Search results

The purpose of this paper is to present the status of the on-going development of the new computerized environment for aircraft synthesis and integrated optimization methods (CEASIOM) and to compare results of different aerodynamic tools. The concurrent design of aircraft is an extremely interdisciplinary activity incorporating simultaneous consideration of complex, tightly coupled systems, functions and requirements. The design task is to achieve an optimal integration of all components into an efficient, robust and reliable aircraft with high performance that can be manufactured with low technical and financial risks, and has an affordable life-cycle cost.

CEASIOM (www.ceasiom.com) is a framework that integrates discipline-specific tools like computer-aided design, mesh generation, computational fluid dynamics (CFD), stability and control analysis and structural analysis, all for the purpose of aircraft conceptual design.

A new CEASIOM version is under development within EU Project AGILE (www.agile-project.eu), by adopting the CPACS XML data-format for representation of all design data pertaining to the aircraft under development.

Results obtained from different methods have been compared and analyzed. Some differences have been observed; however, they are mainly due to the different physical modelizations that are used by each of these methods.

This paper summarizes the current status of the development of the new CEASIOM software, in particular for the following modules: CPACS file visualizer and editor CPACSupdater (Matlab) Automatic unstructured (Euler) & hybrid (RANS) mesh generation by sumo Multi-fidelity CFD solvers: Digital Datcom (Empirical), Tornado (VLM), Edge-Euler & SU2-Euler, Edge-RANS & SU2-RANS Data fusion tool: aerodynamic coefficients fusion from variable fidelity CFD tools above to compile complete aero-table for flight analysis and simulation.

The purpose of this study is to investigate the effects of film holes’ arrangements and jet Reynolds number on flow structure and heat transfer characteristics of jet impingement conjugated with film cooling in a semicylinder double wall channel.

Numerical simulations are used in this research. Streamlines on different sections, skin-friction lines, velocity, wall shear stress and turbulent kinetic energy contours near the concave target wall and vortices in the double channel are presented. Local Nusselt number contours and surface averaged Nusselt numbers are also obtained. Topology analysis is applied to further understand the fluid flow and is used in analyzing the heat transfer characteristics.

It is found that the arrangement of side films positioned far from the center jets helps to enhance the flow disturbance and heat transfer behind the film holes. The heat transfer uniformity for the case of 55° films arrangement angle is most improved and the thermal performance is the highest in this study.

The film holes’ arrangements effects on fluid flow and heat transfer in an impingement cooled concave channel are conducted. The flow structures in the channel and flow characteristics near target by topology pictures are first obtained for the confined film cooled impingement cases. The heat transfer distributions are analyzed with the flow characteristics. The highest heat transfer uniformity and thermal performance situation is obtained in present work.

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.

The purpose of the paper is to predict the aerodynamic performance of a complete scale model H-Darrieus vertical axis wind turbine (VAWT) with end plates at different operating conditions. This paper aims at understanding the flow physics around a model VAWT for three different tip speed ratios corresponding to three different flow regimes.

This study achieves a first three-dimensional hybrid lattice Boltzmann method/very large eddy simulation (LBM-VLES) model for a complete scaled model VAWT with end plates and mast using the solver PowerFLOW. The power curve predicted from the numerical simulations is compared with the experimental data collected at Erlangen University. This study highlights the complexity of the turbulent flow features that are seen at three different operational regimes of the turbine using instantaneous flow structures, mean velocity, pressure iso-contours, blade loading and skin friction plots.

The power curve predicted using the LBM-VLES approach and setup provides a good overall match with the experimental power curve, with the peak and drop after the operational point being captured. Variable turbulent flow structures are seen over the azimuthal revolution that depends on the tip speed ratio (TSR). Significant dynamic stall structures are seen in the upwind phase and at the end of the downwind phase of rotation in the deep stall regime. Strong blade wake interactions and turbulent flow structures are seen inside the rotor at higher TSRs.

The computational cost and time for such high-fidelity simulations using the LBM-VLES remains expensive. Each simulation requires around a week using supercomputing facilities. Further studies need to be performed to improve analytical VAWT models using inputs/calibration from high fidelity simulation databases. As a future work, the impact of turbulent and nonuniform inflow conditions that are more representative of a typical urban environment also needs to be investigated.

The LBM methodology is shown to be a reliable approach for VAWT power prediction. Dynamic stall and blade wake interactions reduce the aerodynamic performance of a VAWT. An ideal operation close to the peak of the power curve should be favored based on the local wind resource, as this point exhibits a smoother variation of forces improving operational performance. The 3D flow features also exhibit a significant wake asymmetry that could impact the optimal layout of VAWT clusters to increase their power density. The present work also highlights the importance of 3D simulations of the complete model including the support structures such as end plates and mast.

Accurate predictions of power performance for Darrieus VAWTs could help in better siting of wind turbines thus improving return of investment and reducing levelized cost of energy. It could promote the development of onsite electricity generation, especially for industrial sites/urban areas and renew interest for VAWT wind farms.

A first high-fidelity simulation of a complete VAWT with end plates and supporting structures has been performed using the LBM approach and compared with experimental data. The 3D flow physics has been analyzed at different operating regimes of the turbine. These physical insights and prediction capabilities of this approach could be useful for commercial VAWT manufacturers.

Reynolds-averaged Navier–Stokes (RANS) models often perform poorly in shock/turbulence interaction regions, resulting in excessive wall heat load and incorrect representation of the separation length in shockwave/turbulent boundary layer interactions. The authors suggest that this can be traced back to inadequate numerical treatment of the inviscid fluxes. The purpose of this study is an extension to the well-known Harten, Lax, van Leer, Einfeldt (HLLE) Riemann solver to overcome this issue.

It explicitly takes into account the broadening of waves due to the averaging procedure, which adds numerical dissipation and reduces excessive turbulence production across shocks. The scheme is derived based on the HLLE equations, and it is tested against three numerical experiments.

Sod’s shock tube case shows that the scheme succeeds in reducing turbulence amplification across shocks. A shock-free turbulent flat plate boundary layer indicates that smooth flow at moderate turbulence intensity is largely unaffected by the scheme. A shock/turbulent boundary layer interaction case with higher turbulence intensity shows that the added numerical dissipation can, however, impair the wall heat flux distribution.

The proposed scheme is motivated by implicit large eddy simulations that use numerical dissipation as subgrid-scale model. Introducing physical aspects of turbulence into the numerical treatment for RANS simulations is a novel approach.

This paper aims to analyze the propulsive performance of small-amplitude pitching foils at very high frequencies with double objectives: to find out scaling laws for the time-averaged thrust and propulsive efficiency at very high frequencies; and to characterize the Strouhal number above which the effect of turbulence on the mean values cannot be neglected.

The thrust force and propulsive efficiency of a pitching NACA0012 foil at high reduced frequencies (k) and a Reynolds number Re = 16 000 are analyzed using accurate numerical simulations, both assuming laminar flow and using a transition turbulence model. The time-averaged results are validated with available experimental data for k up to about 12 (Strouhal number, St, up to 0.6). This study also compares the present numerical results with the predictions of theoretical models and existing numerical results. For a foil pitching about its quarter chord with amplitude α₀ = 8^o, the reduced frequency is varied here up to k = 30 (St up to 2), much higher than in any previous numerical or experimental work.

For this pitch amplitude, turbulence effects are found negligible for St ≲ 0.8, and affecting less than 10% to the time-averaged thrust coefficient CT¯ for larger St Linear potential theory fails for very large k, even for the small pitch amplitude considered, particularly for the power coefficient, and therefore for the propulsive efficiency. It is found that CT¯ ∼ St² for large St, in agreement with recent models, and the propulsive efficiency decays as 1/k, in disagreement with the linear potential theory.

Pitching foils are increasingly studied as efficient propellers and energy harvesting devices. Their performance at very high reduced frequencies has not been sufficiently analyzed before. The authors provide accurate numerical simulations to discern when turbulence is relevant for the computation of the time-averaged thrust and efficiency and how their scaling with the reduced frequency is affected in relation to the laminar-flow predictions. This is relevant because some small-amplitude theoretical models predict high propulsive efficiency of pitching foils at very high frequencies over certain ranges of the structural parameters, and only very accurate numerical simulations may decide on these predictions.

One of the major challenges in the design of thermal equipment is to minimize the entropy production and enhance the thermal dissipation rate for improving energy efficiency of the devices. In several industrial applications, the structure of thermal device is cylindrical shape. In this regard, this paper aims to explore the impact of isothermal cylindrical solid block on nanofluid (Ag – H₂O) convective flow and entropy generation in a cylindrical annular chamber subjected to different thermal conditions. Furthermore, the present study also addresses the structural impact of cylindrical solid block placed at the center of annular domain.

The alternating direction implicit and successive over relaxation techniques are used in the current investigation to solve the coupled partial differential equations. Furthermore, estimation of average Nusselt number and total entropy generation involves integration and is achieved by Simpson and Trapezoidal’s rules, respectively. Mesh independence checks have been carried out to ensure the accuracy of numerical results.

Computations have been performed to analyze the simultaneous multiple influences, such as different thermal conditions, size and aspect ratio of the hot obstacle, Rayleigh number and nanoparticle shape on buoyancy-driven nanoliquid movement, heat dissipation, irreversibility distribution, cup-mixing temperature and performance evaluation criteria in an annular chamber. The computational results reveal that the nanoparticle shape and obstacle size produce conducive situation for increasing system’s thermal efficiency. Furthermore, utilization of nonspherical shaped nanoparticles enhances the heat transfer rate with minimum entropy generation in the enclosure. Also, greater performance evaluation criteria has been noticed for larger obstacle for both uniform and nonuniform heating.

The current numerical investigation can be extended to further explore the thermal performance with different positions of solid obstacle, inclination angles, by applying Lorentz force, internal heat generation and so on numerically or experimentally.

A pioneering numerical investigation on the structural influence of hot solid block on the convective nanofluid flow, energy transport and entropy production in an annular space has been analyzed. The results in the present study are novel, related to various modern industrial applications. These results could be used as a firsthand information for the design engineers to obtain highly efficient thermal systems.

The typical approach of modeling the aerodynamics of an aircraft is to develop a complete database through testing or computational fluid dynamics (CFD). The database will be huge if it has a reasonable resolution and requires an unacceptable CFD effort during the conceptional design. Therefore, this paper aims to reduce the computing effort required via establishing a general aerodynamic model that needs minor parameters.

The model structure was a preconfigured polynomial model, and the parameters were estimated with a recursive method to further reduce the calculation effort. To uniformly disperse the sample points through each step, a unique recursive sampling method based on a Voronoi diagram was presented. In addition, a multivariate orthogonal function approach was used.

A case study of a flying wing aircraft demonstrated that generating a model with acceptable precision (0.01 absolute error or 5% relative error) costs only 1/54 of the cost of creating a database. A series of six degrees of freedom flight simulations shows that the model’s prediction was accurate.

This method proposed a new way to simplify the model and recursive sampling. It is a low-cost way of obtaining high-fidelity models during primary design, allowing for more precise flight dynamics analysis.

To reduce the heat load of a gas turbine blade, its surface is covered with an outer layer of ceramics with high thermal resistance. The purpose of this paper is the selection of ceramics with such a low heat conduction coefficient and thickness, so that the permissible metal temperature is not exceeded on the metal-ceramics interface due to the loss ofmechanical properties.

Therefore, for given temperature changes over time on the metal-ceramics interface, temperature changes over time on the inner side of the blade and the assumed initial temperature, the temperature change over time on the outer surface of the ceramics should be determined. The problem presented in this way is a Cauchy type problem. When analyzing the problem, it is taken into account that thermophysical properties of metal and ceramics may depend on temperature. Due to the thin layer of ceramics in relation to the wall thickness, the problem is considered in the area in the flat layer. Thus, a one-dimensional non-stationary heat flow is considered.

The range of stability of the Cauchy problem as a function of time step, thickness of ceramics and thermophysical properties of metal and ceramics are examined. The numerical computations also involved the influence of disturbances in the temperature on metal-ceramics interface on the solution to the inverse problem.

The computational model can be used to analyze the heat flow in gas turbine blades with thermal barrier.

A number of inverse problems of the type considered in the paper are presented in the literature. Inverse problems, especially those Cauchy-type, are ill-conditioned numerically, which means that a small change in the inputs may result in significant errors of the solution. In such a case, regularization of the inverse problem is needed. However, the Cauchy problem presented in the paper does not require regularization.

Joule heating effect is a pervasive phenomenon in electro-osmotic flow because of the applied electric field and fluid electrical resistivity across the microchannels. Its effect in electro-osmotic flow field is an important mechanism to control the flow inside the microchannels and it includes numerous applications.

This research article details the numerical investigation on alterations in the profile of stream wise velocity of simple Couette-electroosmotic flow and pressure driven electro-osmotic Couette flow by the dynamic viscosity variations happened due to the Joule heating effect throughout the dielectric fluid usually observed in various microfluidic devices.

The advantages of the Joule heating effect are not only to control the velocity in microchannels but also to act as an active method to enhance the mixing efficiency. The results of numerical investigations reveal that the thermal field due to Joule heating effect causes considerable variation of dynamic viscosity across the microchannel to initiate a shear flow when EDL (Electrical Double Layer) thickness is increased and is being varied across the channel.

This research work suggest how joule heating can be used as en effective mechanism for flow control in microfluidic devices.