Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗

A new PNS code for chemical nonequilibrium flows

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic laminar flow of a multicomponent, chemically reacting mixture of thermally perfect gases over two-dimensional and axisymmetric bodies. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, space-marching finite-difference method. The conditions for well-posedness of the space-marching method have been derived from an eigenvalue analysis of the governing equations. The code has been used to compute hypersonic laminar flow of chemically reacting air over wedges and cones. The results of these computations are in good agreement with the results of reacting boundary-layer calculations.

Prabhu, D. K.↗

Use of Navier-Stokes methods for the calculation of high-speed nozzle flow fields

Flows through three reference nozzles have been calculated to determine the capabilities and limitations of the widely used Navier-Stokes solver, PARC. The nozzles examined have similar dominant flow characteristics as those considered for supersonic transport programs. Flows from an inverted velocity profile (IVP) nozzle, an under expanded nozzle, and an ejector nozzle were examined. PARC calculations were obtained with its standard algebraic turbulence model, Thomas, and the two-equation turbulence model, Chien k-epsilon. The Thomas model was run with the default coefficient of mixing set at both 0.09 and a larger value of 0.13 to improve the mixing prediction. Calculations using the default value substantially underpredicted the mixing for all three flows. The calculations obtained with the higher mixing coefficient better predicted mixing in the IVP and underexpanded nozzle flows but adversely affected PARC's convergence characteristics for the IVP nozzle case. The ejector nozzle case did not converge with the Thomas model and the higher mixing coefficient. The Chien k-epsilon results were in better agreement with the experimental data overall than were those of the Thomas run with the default mixing coefficient, but the default boundary conditions for k and epsilon underestimated the levels of mixing near the nozzle exits.

Georgiadis, Nicholas J.↗

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

Measured and predicted vortex-induced leeward heating on a biconic at Mach 6 and 10

Detailed longitudinal and circumferential heating distributions were measured on a spherically blunted, 13 deg/7 deg biconic at angles of attack from 0 deg to 27 deg. The measurements were made in the Langley 20-inch Mach 6 tunnel and 31-inch Mach 10 tunnel to provide heating distributions at Mach 6 and 10 in air. The free-stream Reynolds number based on model length varied from 0.4 to 4.8 million at Mach 6 and 0.3 to 1.3 million at Mach 10. The basic trends in leeward heating are found to be described by a single parameter, the viscous interaction parameter. A decrease in this parameter results in a decrease in leeward heating for attached flow, an increase in heating for separated flow, and a forward movement of the separation region. The parabolized Navier-Stokes code for laminar flow qualitatively predicted the longitudinal and circumferential heating distributions.

Miller, C. G.↗

Hypersonic turbulent wall boundary layer computations

The Baldwin-Lomax (1978) algebraic turbulence model was modified for hypersonic flow conditions. Two coefficients in the outer-layer eddy-viscosity model were determined as functions of Mach number and temperature ratio. By matching the solutions from the Baldwin-Lomax model to those from the Cebeci-Smith (1974) model for a flat plate at hypersonic speed, the new values of the coefficients were obtained. The results show that the values of C(cp) and C(kleb) are functions of both Mach number and wall temperature ratio. The C(cp) and C(kleb) variations with Mach number and wall temperature were used for the calculations of both a 4-deg wedge flow at Mach 18 and an axisymmetric Mach 20 nozzle flow. The Navier-Stokes equations with thin-layer approximation were solved for the above hypersonic flow conditions and the results were compared with existing experimental data. The agreement between the numerical solutions and the existing experimental data were good. The modified Baldwin-Lomax model thus is useful in the computations of hypersonic flows.

Kim, S. C.↗

Transonic Flows About A Fighter Airplane

Report describes computational study of transonic flows about airplane having shape similar to F-16A fighter. Navier-Stokes calculations simulate complicated flow fields and interactions.

Reznick, Steve↗

Aeroloads and secondary flows in a transonic mixed flow turbine stage

A numerical simulation of a transonic mixed flow turbine stage has been carried out using an average passage Navier-Stokes analysis. The mixed flow turbine stage considered here consists of a transonic nozzle vane and a highly loaded rotor. The simulation was run at the design pressure ratio and is assessed by comparing results with those of an established throughflow design system. The three-dimensional aerodynamic loads are studied as well as the development and migration of secondary flows and their contribution to the total pressure loss. The numerical results indicate that strong passage vortices develop in the nozzle vane, mix out quickly, and have little impact on the rotor flow. The rotor is highly loaded near the leading edge. Within the rotor passage, strong spanwise flows and other secondary flows exist along with the tip leakage vortex. The rotor exit loss distribution is similar in character to that found in radial inflow turbines. The secondary flows and non-uniform work extraction also tend to significantly redistribute a non-uniform inlet total temperature profile by the exit of the stage.

Kirtley, K. R.↗

Navier-Stokes Simulation Of Turbine Rotor/Stator Interaction

Patched-grid method used to effect three-dimensional Navier-Stokes numerical simulation of flow interacting with both rotor and stator in turbine. Involves use of two or more grids or sets of grids to solve finite-difference versions of Navier-Stokes or other equations of fluid flow. Ability to perform such calculations contributes to understanding of complicated turbomachinery flows and to consequent improvements in designs of turbomachinery.

Rai, Man Mohan↗

End wall flows in rotors and stators of a single stage compressor

A computer code for solving the parabolized Navier-Stokes equations for internal flows was developed. Oscillations that develop in the calculation procedure are discussed. The measurements made in the hub and annulus wall boundary layers are summarized. The flow in the hub wall boundary layer, starting ahead of the inlet guide vanes to the inlet of the rotor is traced.

Govindan, T. R.↗

Workshop report - A validation study of Navier-Stokes codes for transverse injection into a Mach 2 flow

A CFD (Computational Fluid Dynamics) competition was held at the Third Scramjet Combustor Modeling Workshop to assess the current state-of-the-art in CFD codes for the analysis of scramjet combustors. Solutions from six three-dimensional Navier-Stokes codes were compared for the case of staged injection of air behind a step into a Mach 2 flow. This case was investigated experimentally at the University of Virginia and extensive in-stream data was obtained. Code-to-code comparisons have been made with regard to both accuracy and efficiency. The turbulence models employed in the solutions are believed to be a major source of discrepancy between the six solutions.

Eklund, Dean R.↗

Complete flowfields over low and wide angle AOTV conceptual configurations

Flowfields over two conceptual aeroassisted orbital transfer vehicles, a moderate lift-to-drag ratio biconic, and an axisymmetric zero lift aerobrake are investigated. The solution procedures employ a compressible Navier-Stokes time-asymptotic technique for blunt body flow and a Parabolized Navier-Stokes spatial marching technique for the primarily supersonic shock layer flow over low-angle bodies. Emphasis is placed on the laminar convective heating predictions and comparisons to experimental data where available. Code robustness and solution times are discussed and guidelines are given for the specification of various numerical parameters. The Parabolized Navier-Stokes technique has a high probability of solution success for the range of applications considered herein, and very good comparisons with experimental heat-transfer data have been obtained. The blunt body Navier-Stokes code is more sensitive to parameter specification for the aerobrake calculations. Still, coarse grid solutions can be obtained quickly to get a good first cut at flowfield definition. Fair to good comparisons have been obtained between experimental heating and pressure data and coarse and fine grid solutions.

Gnoffo, P. A.↗

Solutions for hypersonic viscous flow over delta wings

A conical Navier-Stokes solver is developed to study hypersonic flow over delta wings. The differencing scheme used is van Leer flux vector splitting with MUSCL differencing. A flux limiter which involves a second difference in pressure was found to provide sharp, nonoscillatory shocks in windward and leeward flow regions. The governing equations are cast in finite-volume form and are solved using an explicit multistage method. A von Neumann stability analysis of a one-dimensional flux split model is used to predict the stability characteristics of the conical solver. The stability behavior of the scheme is found to strongly depend on the stage coefficients of the multistage scheme and on the order of accuracy of the upwind discretization.

Ruffin, Stephen M.↗

Radiative interactions in chemically reacting supersonic internal flows

The two-dimensional, elliptic Navier-Stokes equations are used to investigate supersonic flows with finite-rate chemistry and radiation for hydrogen-air systems. The chemistry source terms in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The explicit, unsplit MacCormack finite-difference scheme is used to advance the governing equations in time, until convergence is achieved. The specific problem considered is the premixed flow in a channel with a ten-degree compression ramp. Three different chemistry models are used, accounting for increasing number of reactions and participating species. Two chemistry models assure nitrogen as inert, while the third model accounts for nitrogen reactions and NO(x) formation. The tangent slab approximation is used in the radiative flux formulation. A pseudo-gray model is used to represent the absorption-emission characteristics of the participating species. Results obtained for specific conditions indicate that the radiative interactions vary substantially, depending on reactions involving HO2 and NO species and that this can have a significant influence on the flowfield.

Tiwari, S. N.↗

Analysis of viscous transonic flow over airfoil sections

A full Navier-Stokes solver has been used to model transonic flow over three airfoil sections. The method uses a two-dimensional, implicit, conservative finite difference scheme for solving the compressible Navier-Stokes equations. Results are presented as prescribed for the Viscous Transonic Airfoil Workshop to be held at the AIAA 25th Aerospace Sciences Meeting. The NACA 0012, RAE 2822 and Jones airfoils have been investigated for both attached and separated transonic flows. Predictions for pressure distributions, loads, skin friction coefficients, boundary layer displacement thickness and velocity profiles are included and compared with experimental data when possible. Overall, the results are in good agreement with experimental data.

Huff, Dennis L.↗

Analysis of viscous transonic flow over airfoil sections

A full Navier-Stokes solver has been used to model transonic flow over three airfoil sections. The method uses a two-dimensional, implicit, conservative finite difference scheme for solving the compressible Navier-Stokes equations. Results are presented as prescribed for the Viscous Transonic Airfoil Workshop to be held at the AIAA 25th Aerospace Sciences Meeting. The NACA 0012, RAE 2822 and Jones airfoils have been investigated for both attached and separated transonic flows. Predictions for pressure distributions, loads, skin friction coefficients, boundary layer displacement thickness and velocity profiles are included and compared with experimental data when possible. Overall, the results are in good agreement with experimental data.

Huff, Dennis L.↗

Computational studies of an impulsively started viscous flow

Progress in validating incompressible Navier-Stokes codes is described using a predictor/corrector scheme. The flow field under study is the impulsive start of a circular cylinder and the unsteady evolution of the separation bubble. In the current code, a uniform asymptotic expansion is used as an initial condition in order to correctly capture the initial growth of the vortex sheet. Volocity fields at selected instants of time are decomposed into vectorial representations of Navier-Stokes equations which are then used to analyze dominant contributions in the boundary-layer region.

Davis, Sanford S.↗

Vortical flows over delta wings and numerical prediction of vortex breakdown

Navier-Stokes solutions of subsonic vortical flow over a 75 deg sweep delta wing with a sharp leading edge are presented. The sensitivity of the solution to the numerical scheme is examined using both a partially upwind scheme and a scheme with central differencing in all directions. At moderate angles of attack, no vortex breakdown is observed, whereas the higher angle-of-attack cases exhibit breakdown. The effect of numerical grid density is investigated, and solutions that are obtained with various grid densities are compared with experimental data. An embedded grid approach is implemented to enable higher resolution in selected isolated flow regions, such as the leeward-side surface, the leading-edge vortical flow, and the vortex breakdown region.

Ekaterinaris, J. A.↗