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 541 records · Page 30

Computing Flows Using Chimera and Unstructured Grids

DRAGONFLOW is a computer program that solves the Navier-Stokes equations of flows in complexly shaped three-dimensional regions discretized by use of a direct replacement of arbitrary grid overlapping by nonstructured (DRAGON) grid. A DRAGON grid (see figure) is a combination of a chimera grid (a composite of structured subgrids) and a collection of unstructured subgrids. DRAGONFLOW incorporates modified versions of two prior Navier-Stokes-equation-solving programs: OVERFLOW, which is designed to solve on chimera grids; and USM3D, which is used to solve on unstructured grids. A master module controls the invocation of individual modules in the libraries. At each time step of a simulated flow, DRAGONFLOW is invoked on the chimera portion of the DRAGON grid in alternation with USM3D, which is invoked on the unstructured subgrids of the DRAGON grid. The USM3D and OVERFLOW modules then immediately exchange their solutions and other data. As a result, USM3D and OVERFLOW are coupled seamlessly.

Liou, Meng-Sing↗

Comparison of solutions of the Navier-Stokes and interacting boundary layer equations for separated turbulent flow

Comparisons are presented between numerical solutions of the Navier-Stokes and interacting boundary-layer equations for the subsonic separated turbulent flow over an axisymmetric boattail configuration for which experimental data was available. Comparisons were made using a standard algebraic turbulence model and its modification using a streamwise relaxation model. Grid size studies were performed to provide an approximate assessment of the sensitivity of each method to mesh size. In general, good agreement was observed between the solutions of the Navier-Stokes and interacting boundary-layer equations; however, the interaction approach responded more strongly to the modified turbulence model thereby giving better agreement with experimental data.

Vatsa, V. N.↗

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

Navier-Stokes and viscous shock-layer solutions for radiating hypersonic flows

Results are presented from the Navier-Stokes and viscous shock-layer (VSL) calculations with nonequilibrium and equilibrium chemistry, respectively. These calculations contain coupling to the Aerotherm radiation code RAD. A simplified form of the electron energy equation is used to obtain an electron temperature in the Navier-Stokes calculations. The radiation in the flowfield is calculated using this temperature. The Navier-Stokes code is used at high altitude only, whereas the VSL code is employed for the entire entry period to make estimates of the radiative and convective heating to the Fire II vehicle. Results from the Navier-Stokes code have also been compared with the predictions of Lee and Kawamura, who used gray-gas radiation model and thin-layer Navier-Stokes equations. Quite good agreement is obtained between the measured and computed values of radiative and convective heating from the VSL code in th medium-to-low altitude flight regime of the Fire II vehicle. At high altitudes, the Navier-Stokes calculations considerably overpredict the Fire II flight data for radiative intensity. This is attributed to the deficiencies in the Aerotherm radiation model when used for low-density flight conditions. This model contains the thermal equilibrium assumption and precludes accounting for the collision-limiting phenomenon at high altitudes. Present Navier-Stokes calculations highlight the effect of these assumptions on radiative heating calculations for such conditions.

Gupta, Roop N.↗

Navier-Stokes analysis of two- and three-dimensional flow field in solid rocket motors with segment joints

A multidimensional implicit Navier-Stokes analysis which uses numerical solution of ensemble-averaged Navier-Stokes equations in a nonorthogonal bodyfitted cylindrical-polar coordinate system has been applied to simulation of the internal flow field in solid-propellant rocket motor chambers with segment joints. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. Computations performed to simulate the axisymmetric flow field in the vicinity of the aft field joint in the Space Shuttle SRB using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot. Calculations were also performed for two cases involving asymmetric three-dimensional flow in the vicinity of the aft field joint in the SRB using 721,525 grid points to estimate circumferential velocities and pressure gradients at the joint.

Sabnis, J. S.↗

Internal computational fluid mechanics on supercomputers for aerospace propulsion systems

The accurate calculation of three-dimensional internal flowfields for application towards aerospace propulsion systems requires computational resources available only on supercomputers. A survey is presented of three-dimensional calculations of hypersonic, transonic, and subsonic internal flowfields conducted at the Lewis Research Center. A steady state Parabolized Navier-Stokes (PNS) solution of flow in a Mach 5.0, mixed compression inlet, a Navier-Stokes solution of flow in the vicinity of a terminal shock, and a PNS solution of flow in a diffusing S-bend with vortex generators are presented and discussed. All of these calculations were performed on either the NAS Cray-2 or the Lewis Research Center Cray XMP.

Andersen, Bernhard H.↗

A numerical study of viscous vortex rings using a spectral method

Viscous, axisymmetric vortex rings are investigated numerically by solving the incompressible Navier-Stokes equations using a spectral method designed for this type of flow. The results presented are axisymmetric, but the method is developed to be naturally extended to three dimensions. The spectral method relies on divergence-free basis functions. The basis functions are formed in spherical coordinates using Vector Spherical Harmonics in the angular directions, and Jacobi polynomials together with a mapping in the radial direction. Simulations are performed of a single ring over a wide range of Reynolds numbers (Re approximately equal gamma/nu), 0.001 less than or equal to 1000, and of two interacting rings. At large times, regardless of the early history of the vortex ring, it is observed that the flow approaches a Stokes solution that depends only on the total hydrodynamic impulse, which is conserved for all time. At small times, from an infinitely thin ring, the propagation speeds of vortex rings of varying Re are computed and comparisons are made with the asymptotic theory by Saffman. The results are in agreement with the theory; furthermore, the error is found to be smaller than Saffman's own estimate by a factor square root ((nu x t)/R squared) (at least for Re=0). The error also decreases with increasing Re at fixed core-to-ring radius ratio, and appears to be independent of Re as Re approaches infinity). Following a single ring, with Re=500, the vorticity contours indicate shedding of vorticity into the wake and a settling of an initially circular core to a more elliptical shape, similar to Norbury's steady inviscid vortices. Finally, we consider the case of leapfrogging vortex rings with Re=1000. The results show severe straining of the inner vortex core in the first pass and merging of the two cores during the second pass.

Stanaway, S. K.↗

Exploring a Method for Improving Turbulent Separated-Flow Predictions with kappa-omega Models

A particular failing of Reynolds-averaged Navier-Stokes separated turbulent flow computations is addressed within the context of a kappa-omega two-equation turbulence model. The failing is the tendency for turbulence models to under-predict turbulent shear stress in the shear layers of some separation bubbles, yielding late boundary layer reattachment and recovery. Inspired by unpublished work of Volker, Langtry, and Menter, the author undertook an independent investigation in an attempt to improve the ability of the Menter shear stress transport (SST) model to predict flowfield characteristics in and downstream of separation bubbles. The fix is an ad hoc term that is a function of the local ratio of turbulent production to dissipation; it is used to multiply the omega-destruction term, increasing eddy viscosity in separated regions. With this fix, several flowfields are investigated. Results show that, although the "separation fix" can provide dramatic improvement in some cases, it is not consistently good for all flows. Thus, although it may prove helpful in many situations in its current form, this model may benefit from further refinements, including better sensitization to the energetics of turbulence in the separated region.

Rumsey, Christopher L.↗

Navier-Stokes analysis and experimental data comparison of compressible flow in a diffusing S-duct

Full three-dimensional Navier-Stokes computational results are compared with new experimental measurements for the flowfield within a round diffusing S-duct. The present study extends previous computational and experimental results for a similar smaller scale S-duct. Predicted results are compared with the experimental static and total pressure fields, and velocity vectors. Additionally, wall pressures, velocity profiles in wall coordinates, and skin friction values are presented. The CFD results employ algebraic and k-epsilon turbulence models. The CFD computed and experimentally determined separated flowfield is carefully examined.

Harloff, Gary J.↗

Navier-Stokes analysis and experimental data comparison of compressible flow in a diffusing S-duct

Full three-dimensional Navier-Stokes computational results are compared with new experimental measurements for the flowfield within a round diffusing S-duct. The present study extends previous computational and experimental results for a similar smaller scale S-duct. Predicted results are compared with the experimental static and total pressure fields, and velocity vectors. Additionally, wall pressures, velocity profiles in wall coordinates, and skin friction values are presented. The CFD results employ algebraic and k-epsilon turbulence models. The CFD computed and experimentally determined separated flowfield is carefully examined.

Harloff, Gary J.↗

Computation of aircraft component flow fields at transonic Mach numbers using a three-dimensional Navier-Stokes algorithm

A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.

Shrewsbury, George D.↗

Study of hypersonic inlet flow fields with a three-dimensional Navier-Stokes code

A three-dimensional Navier-Stokes code is used to study the flow through an inlet with a rectangular capture area, but with circular or essentially circular internal cross-sections. The inlet may have a center body extending partially or fully over the length of the inlet. Calculations are made under a variety of geometrical constraints imposed on the inlet cross-sections to study the flow quality in the inlet. It is observed that for this class of inlets, the boundary layer has a tendency to thicken and separate readily due to the adverse pressure gradient caused by the continuous compression in the inlet. Since no experimental results are yet available, this paper has presented only the numerical results from the present analysis.

Kumar, A.↗

Application of CFD Analysis to Design Support and Problem Resolution for ASRM and RSRM

The use of Navier-Stokes CFD codes to predict the internal flow field environment in a solid rocket motor is a very important analysis element during the design phase of a motor development program. These computational flow field solutions uncover a variety of potential problems associated with motor performance as well as suggesting solutions to these problems. CFD codes have also proven to be of great benefit in explaining problems associated with operational motors such as in the case of the pressure spike problem with the STS-54B flight motor. This paper presents results from analyses involving both motor design support and problem resolution. The issues discussed include the fluid dynamic/mechanical stress coupling at field joints relative to significant propellant deformations, the prediction of axial and radial pressure gradients in the motor associated with motor performance and propellant mechanical loading, the prediction of transition of the internal flow in the motor associated with erosive burning, the accumulation of slag at the field joints and in the submerged nozzle region, impingement of flow on the nozzle nose, and pressure gradients in the nozzle region of the motor. The analyses presented in this paper have been performed using a two-dimensional axisymmetric model. Fluent/BFC, a three dimensional Navier-Stokes flow field code, has been used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical pressure-velocity coupling algorithm. Wall functions are used to represent the character of the viscous sub-layer flow, and an adjusted k-epsilon turbulence model especially configured for mass injection internal flows, is used to model the growth of turbulence in the motor port. Conclusions discussed in this paper consider flow field effects on the forward, center, and aft propellant grains except for the head end star grain region of the forward propellant segment. The field joints and the submerged nozzle are discussed as well. Conclusions relative to both the design evaluation of the ASRM and the RSRM scenarios explaining the pressure spikes were based on the flow field solutions presented in this paper.

Dill, Richard A.↗

The computation of flow past an oblique wing using the thin-layer Navier-Stokes equations

Essential aspects are presented for computing flow past an oblique wing with the thin-layer Navier-Stokes equations. A new method is developed for generating a grid system around a realistic wing. This method utilizes a series of conformal transformations. The thin-shear-layer approximation and an algebraic eddy-viscosity turbulence model are used to simplify the Reynolds-averaged Navier-Stokes equations. An implicit, factored numerical scheme and the concept of pencil data structure are utilized. For the first time, some flow fields caused by the oblique wing in a supersonic free stream are discussed, emphasizing the separated vortex flows associated with such a wing.

Mehta, Unmeel↗

Dual Stokes cars measurements in supersonic combusting flow

Coherent antiStokes Raman spectroscopy (CARS) was used to make simultaneous measurements of temperature, nitrogen and oxygen density in a reacting supersonic flow. A supersonic burner (SSB) was designed to provide supersonic flow in which combustion can be studied in the laboratory. Measurements made with the CARS system will be used for validation of computational fluid dynamic (CFD) codes. Preliminary results of the CFD calculations are presented.

Jarrett, Olin, Jr.↗

Computation of Navier-Stokes equations for three-dimensional flow separation

Supersonic flows over a sharp and a flat-faced blunt fin mounted on a flat plate are simulated numerically. Several basic issues involved in the resultant three-dimensional steady flow separation are studied. Using the same number of grid points, different grid spacings are employed to investigate the effects of a grid resolution on the origin of the line of separation. Various shock strengths are used to study the so-called separated and unseparated boundary layer and to establish the existence or absence of secondary separation. The length of separation ahead of the flat-faced blunt fin, bifurcation of a horseshoe vortex, and the accessibility of a closed-type separation are investigated. The usual interpretation of the flow field from previous studies and new interpretations arising from the present simulation are discussed.

Hung, Ching-Mao↗