Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsonic 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 289 records · Page 16

Agglomeration Multigrid for an Unstructured-Grid Flow Solver

An agglomeration multigrid scheme has been implemented into the sequential version of the NASA code USM3Dns, tetrahedral cell-centered finite volume Euler/Navier-Stokes flow solver. Efficiency and robustness of the multigrid-enhanced flow solver have been assessed for three configurations assuming an inviscid flow and one configuration assuming a viscous fully turbulent flow. The inviscid studies include a transonic flow over the ONERA M6 wing and a generic business jet with flow-through nacelles and a low subsonic flow over a high-lift trapezoidal wing. The viscous case includes a fully turbulent flow over the RAE 2822 rectangular wing. The multigrid solutions converged with 12%-33% of the Central Processing Unit (CPU) time required by the solutions obtained without multigrid. For all of the inviscid cases, multigrid in conjunction with an explicit time-stepping scheme performed the best with regard to the run time memory and CPU time requirements. However, for the viscous case multigrid had to be used with an implicit backward Euler time-stepping scheme that increased the run time memory requirement by 22% as compared to the run made without multigrid.

Frink, Neal↗

Analysis of boundary conditions for SSME subsonic internal viscous flow analysis

A study was completed of mathematically proper boundary conditions for unique numerical solution of internal, viscous, subsonic flows in the space shuttle main engine. The study has concentrated on well posed considerations, with emphasis on computational efficiency and numerically stable boundary condition statements. The method of implementing the established boundary conditions is applicable to a wide variety of finite difference and finite element codes, as demonstrated.

Baker, A. J.↗

The reversibility theorem for thin airfoils in subsonic and supersonic flow

A method introduced by Munk is extended to prove that the light-curve slope of thin wings in either subsonic flow or supersonic flow is the same when the direction of flight of the wing is reversed. It is also shown that the wing reversal does not change the thickness drag, damping-in-roll parameter or the damping-in-pitch parameter.

Brown, Clinton E↗

Efficient solutions to the Euler equations for supersonic flow with embedded subsonic regions

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two dimensions is described. Convergence of the basic algorithm to the steady state is quadratic for fully supersonic flows and is linear for other flows. This is in contrast to the block alternating direction implicit methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented herein is easily coupled with methods to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, and yields a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, Robert W.↗

Effects of inviscid parallel shear flows on steady and unsteady aerodynamics and flutter

A general theory of planar disturbances in inviscid parallel shear flows, analogous to thin-wing theory in potential flows, has been developed. Integral relations between surface pressure and deformation are obtained which are similar to, and can be solved by the same numerical methods as, those of potential flow. Computed results are shown which illustrate the effects of a model turbulent boundary layer on various lifting and nonlifting surfaces, including an elastic panel in low supersonic flow and an airfoil control surface in subsonic flow.

Williams, M. H.↗

Development of a three-dimensional turbulent duct flow analysis

A method for computing three-dimensional turbulent subsonic flow in curved ducts is described. An approximate set of governing equations is given for viscous flows which have a primary flow direction. The derivation is coordinate invariant, and the resulting equations are expressed in terms of tensors. General tube-like coordinates were developed for a general class of geometries applicable to many internal flow problems. The coordinates are then particularized to pipes having superelliptic cross sections whose shape can vary continuously between a circle and a near rectangle. The analysis is applied to a series of relevant aerodynamic problems including transition from nearly square to round pipes and flow through a pipe with an S-shaped bend.

Eiseman, P. R.↗

An inverse method for subcritical flows

An exact method is presented for two-dimensional subsonic flow within the limitations of the tangent gas approximation. While Woods (1961) studied these equations and proposed iterative methods for the solution of both the iterative and inverse problems, the inverse method here presented is noniterative and exact. It is shown that the direct Euler solution over the designed airfoil is very close to the input speed distribution, and that the constraints necessitated by upstream condition and closure requirements are easily incorporated.

Daripa, P. K.↗

The similarity rules for second-order subsonic and supersonic flow

The similarity rules for linearized compressible flow theory (Gothert's rule and its supersonic counterpart) are extended to second order. It is shown that any second-order subsonic flow can be related to "nearly incompressible" flow past the same body, which can be calculated by the Janzen-Rayleigh method.

Van Dyke, Milton D↗

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

Solution of transonic flows by an integro-differential equation method

Solutions of steady transonic flow past a two-dimensional airfoil are obtained from a singular integro-differential equation which involves a tangential derivative of the perturbation velocity potential. Subcritical flows are solved by taking central differences everywhere. For supercritical flows with shocks, central differences are taken in subsonic flow regions and backward differences in supersonic flow regions. The method is applied to a nonlifting parabolic-arc airfoil and to a lifting NACA 0012 airfoil. Results compare favorably with those of finite-difference schemes.

Ogana, W.↗

Numerical simulation of separated and vortical flows on bodies at large angles of attack

Navier-Stokes computations were carried out for high-angle-of-attack flow about bodies of revolution in supersonic and subsonic flow. Numerical and physical factors affecting flow predictions are discussed. Accurate prediction of the complex three-dimensional separated flow is found to depend on using grids sufficiently fine to resolve the details of both the viscous boundary layer and the off-surface separated flow structures, and, for turbulent flow cases, use of an eddy-viscosity turbulence model which accounts for the leeward vortical flow structures. An example of high-incidence flow computed about a more complex geometry, the F-18 fighter fuselage forebody and leading edge extension, is also presented.

Schiff, Lewis B.↗

Viscous/Inviscid Interaction On An Oscillating Airfoil

Viscous/inviscid-interaction procedure developed for computation of two-dimensional flow of air about airfoil, angle of attack of which oscillates at high amplitude. Purpose of procedure to study flow effects in dynamic stall. Bulk of flow treated as inviscid, while effects of viscosity treated only on surface of airfoil and in boundary layer. Inviscid part of flow simulated via time-marching, unsteady-subsonic-flow panel method incorporating appropriate distributions of unsteady surface and wake singularities.

Carr, L. W.↗

An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids

An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes. Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3 element and 4 element airfoil configurations. For the calculations, two one equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3 element configuration. For the 4 element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.

Anerson, W. Kyle↗

Application of a Full Reynolds Stress Model to High Lift Flows

A recently developed second-moment Reynolds stress model was applied to two challenging high-lift flows: (1) transonic flow over the ONERA M6 wing, and (2) subsonic flow over the DLR-F11 wing-body configuration from the second AIAA High Lift Prediction Workshop. In this study, the Reynolds stress model results were contrasted with those obtained from one- and two{equation turbulence models, and were found to be competitive in terms of the prediction of shock location and separation. For an ONERA M6 case, results from multiple codes, grids, and models were compared, with the Reynolds stress model tending to yield a slightly smaller shock-induced separation bubble near the wing tip than the simpler models, but all models were fairly close to the limited experimental surface pressure data. For a series of high-lift DLR{F11 cases, the range of results was more limited, but there was indication that the Reynolds stress model yielded less-separated results than the one-equation model near maximum lift. These less-separated results were similar to results from the one-equation model with a quadratic constitutive relation. Additional computations need to be performed before a more definitive assessment of the Reynolds stress model can be made.

Lee-Rausch, E. M.↗

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

Boundary-value problem of configurations with compressible free vortex flow

A self-consistent version of the compressible boundary-value problem for configurations with leading-edge vortex separation is formulated, based on the assumption that the compressible flow field is controlled by the linearized potential equation. The stream surface boundary condition and the zero pressure jump condition of the compressible free vortex flows are analyzed; application of the Goethert rule permits the compressible nonlinear boundary-value problem for the subsonic flow domain to be transformed into an equivalent nonlinear incompressible problem. The compressibility corrections developed are used in numerical calculations of subsonic leading-edge vortex flows about planar wing geometries. The sample calculations, employing an inviscid flow model in which the wing and vortex sheets are represented by piecewise continuous quadratic doublet sheet distributions, are applicable to high subsonic Mach numbers.

Brune, G. W.↗

Computed And Measured Flow In A Turnaround Duct

Report presents theoretical and experimental study of subsonic flow of air in 180 degree-turnaround duct of rectangular cross section. Purpose of study to conduct partial simulation of conditions similar to those in axisymmetric duct of Space Shuttle main engine. Part of continuing effort to develop capability to simulate complicated turbulent internal flows with large curvatures like those in rocket engines and turbomachinery.

Monson, Daryl J.↗