Search results

This study aims to propose an aerodynamic force decomposition which, for the first time, allows for thrust/drag bookkeeping in two-dimensional viscous and unsteady flows. Lamb vector-based far-field methods are used at the scope, and the paper starts with extending recent steady compressible formulas to the unsteady regime.

Exact vortical force formulas are derived considering inertial or non-inertial frames, viscous or inviscid flows, fixed or moving bodies. Numerical applications to a NACA0012 airfoil oscillating in pure plunging motion are illustrated, considering subsonic and transonic flow regimes. The total force accuracy and sensitivity to the control volume size is first analysed, then the axial force is decomposed and results are compared to the inviscid force (thrust) and to the steady force (drag).

Two total axial force decompositions in thrust and drag contributions are proposed, providing satisfactory results. An additional force decomposition is also formulated, which is independent of the arbitrary pole appearing in vortical formulas. Numerical inaccuracies encountered in inertial reference frames are eliminated, and the extended formulation also allows obtaining an accurate force prediction in presence of shock waves.

No thrust/drag bookkeeping methodology was actually available for oscillating airfoils in viscous and compressible flows.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

The purpose of the paper is to analyse the performances of closures and compressibility corrections classically used in turbulence models when applied to highly-compressible turbulent boundary layers (TBLs) over flat plates.

A direct numerical simulation (DNS) database of TBLs, covering a wide range of thermodynamic conditions, is presented and exploited to perform a priori analyses of classical and recent closures for turbulent models. The results are systematically compared to the “exact” terms computed from DNS.

The few compressibility corrections available in the literature are not found to capture DNS data much better than the uncorrected original models, especially at the highest Mach numbers. Turbulent mass and heat fluxes are shown not to follow the classical gradient diffusion model, which was shown instead to provide acceptable results for modelling the vibrational turbulent heat flux.

The main originality of the present paper resides in the DNS database on which the a priori tests are conducted. The database contains some high-enthalpy simulations at large Mach numbers, allowing to test the performances of the turbulence models in the presence of both chemical dissociation and vibrational relaxation processes.

This study aims to evaluate blast loads on and the response of submerged structures.

An arbitrary Lagrangian–Eulerian method is developed to model fluid–structure interaction (FSI) problems of close-in underwater explosions (UNDEX). The “fluid” part provides the loads for the structure considers air, water and high explosive materials. The spatial discretization for the fluid domain is performed with a second-order vertex-based finite volume scheme with a tangent of hyperbola interface capturing technique. The temporal discretization is based on explicit Runge–Kutta methods. The structure is described by a large-deformation Lagrangian formulation and discretized via finite elements. First, one-dimensional test cases are given to show that the numerical method is free of mesh movement effects. Thereafter, three-dimensional FSI problems of close-in UNDEX are studied. Finally, the computation of UNDEX near a ship compartment is performed.

The difference in the flow mechanisms between rigid targets and deforming targets is quantified and evaluated.

Cavitation is modeled only approximately and may require further refinement/modeling.

The results demonstrate that the proposed numerical method is accurate, robust and versatile for practical use.

Better design of naval infrastructure [such as bridges, ports, etc.].

To the best of the authors’ knowledge, this study has been conducted for the first time.

As to different angles of attack and nonlinear problems caused by high temperatures in coexisting hypersonic aircraft, people mainly rely on fluid software for research but lack analysis of flow mechanisms. Owing to computational difficulties, few people use numerical algorithms to combine them for discussion. Hence, this study aims to make a deep inquiry into the laminar flow and heat transfer of compressible Newtonian fluid in hypersonic aircraft with small attack angles.

In this paper, on the basis of mass, momentum and energy conservation laws, the governing equations of the hypersonic boundary layer are established. Viscosity, specific heat capacity and thermal conductivity are considered nonlinear functions concerning temperature. In virtue of the MacCormack finite difference method, the stationary numerical solutions are solved directly, and the validity of the algorithm is verified.

The results demonstrate that at Mach number 5, compared to the 0° attack angle, the maximum temperature near-wall at the 3° attack angle increases by about 25%. An enjoyable phenomenon is discovered, where the position corresponding to the maximum wall shear force shifts back as the attack angle and Mach number increase. The relationship between the near-wall maximum temperature versus attack angle and Mach number is fitted through numerical calculation results.

Empirical formulas can be used to estimate heat transfer characteristics at small attack angles, which will guide the design of aircraft thermal protection systems.

The mixing of fuel and air plays a pivotal role in enhancing combustion in supersonic regime. Proper mixing stabilizes the flame and prevents blow-off. Blow-off is due to the shorter residence time of fuel and air in the combustor, as the flow is in supersonic regime. The flame is initiated in the local subsonic region created using a flameholder within the supersonic combustor. This study aims to design an effective flameholder which increases the residence time of fuel in the combustor allowing proper combustion preventing blow-off and other instabilities.

The geometry of the strut-based flameholder is altered in the present study to induce a streamwise motion of the fluid downstream of the strut. The streamwise motion of the fluid is initiated by the ramps and grooves of the strut geometry. The numerical simulations were carried out using ANSYS Fluent and are validated against the available experimental and numerical results of cold flow with hydrogen injection using plain strut as the flameholder. In the present study, numerical investigations are performed to analyse the effect on hydrogen injection in strut-based flameholders with ramps and converging grooves using Reynolds-averaged Navier–Stokes equation coupled with Menter’s shear stress transport k-ω turbulence model. The analysis is done to determine the effect of geometrical parameters and flow parameter on the flow structures near the base of the strut where thorough mixing takes place. The geometrical parameters under consideration include the ramp length, groove convergence angle, depth of the groove, groove compression angle and the Mach number. Two different strut configurations, namely, symmetric and asymmetric struts were also studied.

Higher turbulence and complex flow structures are visible in asymmetric strut configuration which develops better mixing of hydrogen and air compared to symmetric strut configuration. The variation in the geometric parameters develop changes in the fluid motion downstream of the strut. The fluid passing through the converging grooves gets decelerated thereby reducing the Mach number by 20% near the base of the strut compared to the straight grooved strut. The shorter ramps are found to be more effective, as the pressure variation in lateral direction is carried along the strut walls downstream of the strut increasing the streamwise motion of the fluid. The decrease in the depth of the groove increases the recirculation zone downstream of the strut. Moreover, the increase in the groove compression angle also increases the turbulence near the base of the strut where the fuel is injected. Variation in the injection port location increases the mixing performance of the combustor by 25%. The turbulence of the fuel jet stream is considerably changed by the increase in the injection velocity. However, the change in the flow field properties within the flow domain is marginal. The increase in fuel mass flow rate brings about considerable change in the flow field inducing stronger shock structures.

The present study identifies the optimum geometry of the strut-based flameholder with ramps and converging grooves. The reaction flow modelling may be performed on the strut geometry incorporating the design features obtained in the present study.

The purpose of this study is to investigate the combustion of the n-Heptane droplets in the supersonic combustor with a cavity-based fuel injection configuration. The focus is on the impacts of the droplet size on combustion efficiency.

The finite volume solver is developed to simulate the two-phase reacting turbulent compressible flow using a single step reaction mechanism as finite rate chemistry. Three different fuel injection settings are studied for the considered physical geometry and flow conditions: the gas fuel injection, small droplet liquid fuel injection and big droplet fuel. The fuel is injected as a slot wall jet from the bottom of the cavity.

The results show that using the small droplet size, the complete fuel consumption and combustion efficiency can be achieved but using the big droplet sizes, most fuel exit the combustor in the liquid phase and gasified unburned fuel. It is also demonstrated that the cavity's temperature distribution of the liquid fuel case is different from the gas fuel, and two flame branches are observed there due to the droplet evaporation and combustion in the cavity.

To the best of the authors’ knowledge, this study is performed for the first time on the combustion of the n-Heptane fuel droplets in scramjet configuration, which is promising propulsion system for the future economic flights.

This paper aims to study passive control techniques for transonic flow over a backward-facing step (BFS) using square-lobed trailing edges. The study investigates the efficacy of upward and downward lobe patterns, different lobe widths and deflection angles on flow separation, aiming for a deeper understanding of the flow physics behind the passive flow control system.

Large Eddy Simulation and Reynolds-averaged Navier–Stokes were used to evaluate the results of the study. The research explores the impact of upward and downward patterns of lobes on flow separation through the effects of different lobe widths and deflection angles. Numerical methods are used to analyse the behaviour of transonic flow over BFS and compared it to existing experimental results.

The square-lobed trailing edges significantly enhance the reduction of mean reattachment length by up to 80%. At Ma = 0.8, the up-downward configuration demonstrates increased effectiveness in reducing the root mean square of pressure fluctuations at a proximity of 5-step height in the wake region, with a reduction of 50%, while the flat-downward configuration proves to be more efficient in reducing the root mean square of pressure fluctuations at a proximity of 1-step height in the near wake region, achieving a reduction of 71%. Furthermore, the study shows that the up-downward configuration triggers early spanwise velocity fluctuations, whereas the standalone flat-downward configuration displays less intense crosswise velocity fluctuations within the wake region.

The findings demonstrate the effectiveness of square-lobed trailing edges as passive control techniques, showing significant implications for improving efficiency, performance and safety of the design in aerospace and industrial systems.

This paper demonstrates that the square-lobed trailing edges are effective in reducing the mean reattachment length and pressure fluctuations in transonic conditions. The study evaluates the efficacy of different configurations, deflection angles and lobe widths on flow and provides insights into the flow physics of passive flow control systems.

The purpose of this paper is to investigate the impact of near-wall treatment approaches, which are crucial parameters in predicting the flow characteristics of open channels, and the influence of different vegetation covers in different layers.

Ansys Fluent, a computational fluid dynamics software, was used to calculate the flow and turbulence characteristics using a three-dimensional, turbulent (k-e realizable), incompressible and steady-flow assumption, along with various near-wall treatment approaches (standard, scalable, non-equilibrium and enhanced) in the vegetated channel. The numerical study was validated concerning an experimental study conducted in the existing literature.

The numerical model successfully predicted experimental results with relative error rates below 10%. It was determined that nonequilibrium wall functions exhibited the highest predictive success in experiment Run 1, standard wall functions in experiment Run 2 and enhanced wall treatments in experiment Run 3. This study has found that plant growth significantly alters open channel flow. In the contact zones, the velocities and the eddy viscosity are low, while in the free zones they are high. On the other hand, the turbulence kinetic energy and turbulence eddy dissipation are maximum at the solid–liquid interface, while they are minimum at free zones.

This is the first study, to the best of the author’s knowledge, concerning the performance of different near-wall treatment approaches on the prediction of vegetation-covered open channel flow characteristics. And this study provides valuable insights to improve the hydraulic performance of open-channel systems.

The purpose of this paper is to investigate the two-phase flow problems involving gas–liquid mixture.

The governed equations consist of a range of conservation laws modeling a classification of two-phase flow phenomena subjected to a velocity nonequilibrium for the gas–liquid mixture. Effects of the relative velocity are accounted for in the present model by a kinetic constitutive relation coupled to a collection of specific equations governing mass and volume fractions for the gas phase. Unlike many two-phase models, the considered system is fully hyperbolic and fully conservative. The suggested relaxation approach switches a nonlinear hyperbolic system into a semilinear model that includes a source relaxation term and characteristic linear properties. Notably, this model can be solved numerically without the use of Riemann solvers or linear iterations. For accurate time integration, a high-resolution spatial reconstruction and a Runge–Kutta scheme with decreasing total variation are used to discretize the relaxation system.

The method is used in addressing various nonequilibrium two-phase flow problems, accompanied by a comparative study of different reconstructions. The numerical results demonstrate the suggested relaxation method’s high-resolution capabilities, affirming its proficiency in delivering accurate simulations for flow regimes characterized by strong shocks.

While relaxation methods exhibit notable performance and competitive features, as far as we are aware, there has been no endeavor to address nonequilibrium two-phase flow problems using these methods.