Search results

An aerothermodynamic design code for axisymmetric projectiles has been developed using a viscous‐inviscid interaction scheme. Separate solution procedures for the inviscid and the viscous (boundary layer) fluid dynamic equations are coupled by an iterative solution procedure. Non‐equilibrium, equilibrium and perfect gas boundary layer equations are included. The non‐equilibrium gas boundary layer equations assume a binary mixture (two species; atoms and molecules) of chemically reacting perfect gases. Conservation equations for each species include finite reaction rates applicable to high temperature air. The equilibrium gas boundary layer equations assume infinite rate reactions, while the perfect gas equations assume no chemical reactions. Projectile near‐wall and surface flow profiles (velocity, pressure, density, temperature and heat transfer) representing converged solutions to both the inviscid and viscous equations can be obtained in less than two minutes on minicomputers. A technique for computing local reverse flow regions is included. Computations for yawed projectiles are accomplished using a coordinate system transformation technique that is valid for small angle‐of‐attack. Computed surface pressure, heat transfer rates and aerodynamic forces and moments for 1.25 &le Mach No. &le 10.5 are compared to wind tunnel and free flight measurements on flat plate, blunt‐cone, and projectile geometries such as a cone‐cylinder‐flare.

The purpose of this paper is to introduce a unique technique to couple the two-integral boundary layer solutions to a generic inviscid solver in an iterative fashion.

The boundary layer solution is obtained using the two-integral method to solve displacement thickness point by point with a local Newton method, at a fraction of the cost of a conventional mesh-based, full viscous solution. The boundary layer solution is coupled with an existing inviscid solver. Coupling occurs by moving the wall to a streamline at the computed boundary layer thickness and treating it as a slip boundary, then solving the flow again and iterating. The Goldstein singularity present when solving boundary layer equations is overcome by solving an auxiliary velocity equation along with the displacement thickness.

The proposed method obtained favourable results when compared with the analytical solutions for flat and inclined plates. Further, it was applied to modelling the flow around a NACA0012 airfoil and yielded results similar to those of the widely used XFOIL code.

A unique method is proposed for coupling of the boundary layer solution to the inviscid flow. Rather than the traditional transpiration boundary condition, mesh movement is employed to simulate the boundary layer thickness in a more physically meaningful way. Further, a new auxiliary velocity equation is presented to circumvent the Goldstein singularity.

Two numerical methods, based on high order finite volume formulations and upwind schemes, are used to compute the two‐ and three‐dimensional flow field in a transonic nozzle. The influence of numerical diffusivity, boundary treatment and mesh structure is explored for inviscid and turbulent configurations. First order computations provide significantly different inviscid results. However, high order methods lead to similar solutions. An explanation of the error generated through the shockwave is proposed in this case. The two‐dimensional interaction of the shock with the thin turbulent boundary layer developing on the bump wall is also presented. Good agreement between both approaches is obtained considering the rapid thickening of the boundary layer due to the shock. Furthermore, the downstream velocity recovery is almost identical. Only slight discrepancies occur in the main flow, near the outer edge of the boundary layer. These seem to be related to the way the turbulence model deals with the free stream turbulence. Finally, preliminary three‐dimensional unstructured turbulent results are presented and discussed.

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.

The first compressible flow solution based solely on the locally analytical method is developed. This is accomplished by developing the flow model and locally analytical solution for inviscid subsonic compressible flow. The stream function for irrotational, compressible flow without body forces was chosen as the governing differential equation. To demonstrate the modelling and locally analytical solution, this analysis is then applied to predict the flow in convergent nozzles, both planar and axially symmetric, for different back pressures. Results are presented which demonstrate the effectiveness of this technique.

An upwind unstructured grid cell‐centred scheme for the solution of the compressible Euler and Navier‐Stokes equations in two dimensions is presented. The algorithm employs a finite volume formulation. Calculation of the inviscid fluxes is based on the approximate Riemann solver of Roe. Viscous fluxes are obtained from solution gradients computed by a variational recovery procedure. Higher order accuracy is achieved through performing a monotonic linear reconstruction of the solution over each cell. The steady state is obtained by a point implicit time integration of the unsteady equations using local time stepping. For supersonic inviscid flow an alternative space marching algorithm is proposed. This latter approach is applicable to supersonic flow fields containing regions of local subsonic flow. Numerical results are presented to show the performance of the proposed scheme for inviscid and viscous flows.

The purpose of this study is to examine theoretically the axisymmetric flow of a steady free-surface jet emerging from a tube for high inertia flow and moderate surface tension effect.

The method of matched asymptotic expansion is used to explore the rich dynamics near the exit where a stress singularity occurs. A boundary layer approach is also proposed to capture the flow further downstream where the free surface layer has grown significantly.

The jet is found to always contract near the tube exit. In contrast to existing numerical studies, the author explores the strength of upstream influence and the flow in the wall layer, resulting from jet contraction. This influence becomes particularly evident from the nonlinear pressure dependence on the upstream distance, as well as the pressure undershoot and overshoot at the exit for weak and strong gravity levels, respectively. The approach is validated against existing experimental and numerical data for the jet profile and centerline velocity where good agreement is obtained. Far from the exit, the author shows how the solution in the diffusive region can be matched to the inviscid far solution, providing the desired appropriate initial condition for the inviscid far flow solution. The location, at which the velocity becomes uniform across the jet, depends strongly on the gravity level and exhibits a non-monotonic behavior with respect to gravity and applied pressure gradient. The author finds that under weak gravity, surface tension has little influence on the final jet radius. The work is a crucial supplement to the existing numerical literature.

Given the presence of the stress singularity at the exit, the work constitutes a superior alternative to a computational approach where the singularity is typically and inaccurately smoothed over. In contrast, in the present study, the singularity is entirely circumvented. Moreover, the flow details are better elucidated, and the various scales involved in different regions are better identified.

This paper aims to introduce physics-informed neural networks (PINN) applied to the two-dimensional steady-state laminar Navier–Stokes equations over a flat plate with roughness elements and specified local heating. The method bridges the gap between asymptotics theory and three-dimensional turbulent flow analyses, characterized by high costs in analysis setups and prohibitive computing times. The results indicate the possibility of using surface heating or wavy surface to control the incoming flow field.

The understanding of the flow control mechanism is normally caused by the unsteady interactions between the aircraft structure and the turbulent flows as well as some studies have shown, surface roughness can significantly influence the fluid dynamics by inducing perturbations in the velocity profile.

The description of the boundary-layer flow, based upon a triple-deck structure, shows how a wavy surface and a local surface heating generate an interaction between the inviscid region and the viscous region near the flat plate.

To the best of the authors’ knowledge, the presented approach is especially original in relation to the innovative concept of PINN as a solver of the asymptotic triple-deck method applied to the viscous–inviscid boundary layer interaction.

This study aims to feature the application of the artificial compressibility method (ACM) for the numerical prediction of two-dimensional (2D) axisymmetric swirling flows.

The respective academic numerical solver, named IGal2D, is based on the axisymmetric Reynolds-averaged Navier–Stokes (RANS) equations, arranged in a pseudo-Cartesian form, enhanced by the addition of the circumferential momentum equation. Discretization of spatial derivative terms within the governing equations is performed via unstructured 2D grid layouts, with a node-centered finite-volume scheme. For the evaluation of inviscid fluxes, the upwind Roe’s approximate Riemann solver is applied, coupled with a higher-order accurate spatial reconstruction, whereas an element-based approach is used for the calculation of gradients required for the viscous ones. Time integration is succeeded through a second-order accurate four-stage Runge-Kutta method, adopting additionally a local time-stepping technique. Further acceleration, in terms of computational time, is achieved by using an agglomeration multigrid scheme, incorporating the full approximation scheme in a V-cycle process, within an efficient edge-based data structure.

A detailed validation of the proposed numerical methodology is performed by encountering both inviscid and viscous (laminar and turbulent) swirling flows with axial symmetry. IGal2D is compared against the commercial software ANSYS fluent – by using appropriate metrics and characteristic flow quantities – but also against experimental measurements, confirming the proposed methodology’s potential to predict such flows in terms of accuracy.

This study provides a robust methodology for the accurate prediction of swirling flows by combining the axisymmetric RANS equations with ACM. In addition, a detailed description of the convective flux Jacobian is provided, filling a respective gap in research literature.

The purpose of this paper is to explore the characteristics of hypersonic flow past a blunt body.

The implicit finite volume schemes are derived from axisymmetric Navier–Stokes equations by means of AUSM+ and LU-SGS methods, and programmed in FORTRAN. Based on the verified result that a 2D axisymmetric chemical equilibrium flow has a good agreement with the literature, the characteristics of hypersonic flow past a sphere are simulated by using four different models which involve four factors, namely, viscous, inviscid, equilibrium and calorically perfect gas.

Compared with the calorically perfect gas under hypervelocity condition, the shock wave of the equilibrium gas is more close to the blunt body, gas density and pressure become bigger, but gas temperature is lower due to the effect of real gas. Viscous effects are not obvious in the calculations of the equilibrium gas or the calorically perfect gas. In a word, the model of equilibrium gas is more suitable for hypersonic flow and the calculation of viscous flow has a smaller error.

The computer codes are developed to simulate the characteristics of hypersonic flows, and this study will be helpful for the design of the thermal protection system in hypersonic vehicles.