Search results

A higher-order implicit shock-capturing scheme is presented for the Euler equations based on time linearization of the implicit flux vector rather than the residual vector.

The flux vector is linearized through a truncated Taylor-series expansion whose leading-order implicit term is an inner product of the flux Jacobian and the vector of differences between the current and previous time step values of conserved variables. The implicit conserved-variable difference vector is evaluated at cell faces by using the reconstructed states at the left and right sides of a cell face and projecting the difference between the left and right states onto the right eigenvectors. Flux linearization also facilitates the construction of implicit schemes with higher-order spatial accuracy (up to third order in the present study). To enhance the diagonal dominance of the coefficient matrix and thereby increase the implicitness of the scheme, wave strengths at cell faces are expressed as the inner product of the inverse of the right eigenvector matrix and the difference in the right and left reconstructed states at a cell face.

The accuracy of the implicit algorithm at Courant–Friedrichs–Lewy (CFL) numbers greater than unity is demonstrated for a number of test cases comprising one-dimensional (1-D) Sod’s shock tube, quasi 1-D steady flow through a converging-diverging nozzle, and two-dimensional (2-D) supersonic flow over a compression corner and an expansion corner.

The algorithm has the advantage that it does not entail spatial derivatives of flux Jacobian so that the implicit flux can be readily evaluated using Roe’s approximate Jacobian. As a result, this approach readily facilitates the construction of implicit schemes with high-order spatial accuracy such as Roe-MUSCL.

A novel finite-volume-based higher-order implicit shock-capturing scheme was developed that uses time linearization of fluxes at cell interfaces.

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.

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.

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.

This study aims to investigate the morphing concepts able to manipulate the dynamics of the downstream unsteadiness in the separated shear layers and, in the wake, be able to modify the upstream shock–boundary layer interaction (SBLI) around an A320 morphing prototype to control these instabilities, with emphasis to the attenuation or even suppression of the transonic buffet. The modification of the aerodynamic performances according to a large parametric study carried out at Reynolds number of 4.5 × 10⁶, Mach number of 0.78 and various angles of attack in the range of (0, 2.4)° according to two morphing concepts (travelling waves and trailing edge vibration) are discussed, and the final benefits in aerodynamic performance increase are evaluated.

This article examines through high fidelity (Hi-Fi) numerical simulation the effects of the trailing edge (TE) actuation and of travelling waves along a specific area of the suction side starting from practically the most downstream position of the shock wave motion according to the buffet and extending up to nearly the TE. The present paper studies through spectral analysis the coherent structures development in the near wake and the comparison of the aerodynamic forces to the non-actuated case. Thus, the physical mechanisms of the morphing leading to the increase of the lift-to-drag ratio and the drag and noise sources reduction are identified.

This study investigates the influence of shear-layer and near-wake vortices on the SBLI around an A320 aerofoil and attenuation of the related instabilities thanks to novel morphing: travelling waves generated along the suction side and trailing-edge vibration. A drag reduction of 14% and a lift-to-drag increase in the order of 8% are obtained. The morphing has shown a lift increase in the range of (1.8, 2.5)% for angle of attack of 1.8° and 2.4°, where a significant lift increase of 7.7% is obtained for the angle of incidence of 0° with a drag reduction of 3.66% yielding an aerodynamic efficiency of 11.8%.

This paper presents results of morphing A320 aerofoil, with a chord of 70cm and subjected to two actuation kinds, original in the state of the art at M = 0.78 and Re = 4.5 million. These Hi-Fi simulations are rather rare; a majority of existing ones concern smaller dimensions. This study showed for the first time a modified buffet mode, displaying periodic high-lift “plateaus” interspersed by shorter lift-decrease intervals. Through trailing-edge vibration, this pattern is modified towards a sinusoidal-like buffet, with a considerable amplitude decrease. Lock-in of buffet frequency to the actuation is obtained, leading to this amplitude reduction and a drastic aerodynamic performance increase.