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 595 records · Page 33

Investigation of upwind, multigrid, multiblock numerical schemes for three dimensional flows. Volume 1: Runge-Kutta methods for a thin layer Navier-Stokes solver

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.

Cannizzaro, Frank E.↗

Flow visualization and pressure distributions for an all-body hypersonic aircraft

A CFD code-validation effort has been conducted at the NASA-Ames 3.5-ft hypersonic wind tunnel, using a generic, 'all-body' hypersonic aircraft configuration model. The CFD methods to be validated encompass approximate inviscid ones and the upwind parabolized Navier-Stokes solver code. Flow visualizations and pressure distributions are obtained for the cases of zero and 15 deg angles of attack. A complex leeward flow is observed at angle-of-attack with crossflow separation and vortices, and significant changes are noted in the transition from the forebody's conical to the afterbody's nonconical flows.

Lockman, William K.↗

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗

Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver

A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Kim, S.-W.↗

Rotational CARS application to simultaneous and multiple-point temperature and concentration determination in a turbulent flow

Coherent anti-Stokes Raman scattering (CARS) from the pure rotational Raman lines of N2 is employed to measure the instantaneous (approximately 10 ns) rotational temperature of N2 gas at room temperature and below with good spatial resolution (0.2 x 0.2 x 3.0 cu mm). A broad bandwidth dye laser is used to obtain the entire rotational spectrum from a single laser pulse; the CARS signal is then dispersed by a spectrograph and recorded on an optical multichannel analyzer. A best fit temperature is found in several seconds with the aid of a computer for each experimental spectrum by a least squares comparison with calculated spectra. The model used to calculate the theoretical spectra incorporates the temperature and pressure dependence of the pressure-broadened rotational Raman lines, includes the nonresonant background susceptibility, and assumes that the pump laser has a finite linewidth. Temperatures are fit to experimental spectra recorded over the temperature range of 135 to 296 K, and over the pressure range of .13 to 15.3 atm.

Snow, J. B.↗

Rotational cars application to simultaneous and multiple-point temperature and concentration determination in a turbulent flow

Coherent Anti-stokes Raman Scattering (CARS) from the pure rotational Raman lines of N2 is employed to measure the instantaneous rotational temperature of N2 gas at room temperature and below with good spatial resolution. A broad-bandwidth dye laser is used to obtain the entire rotational spectrum from a signal laser pulse; the CARS signal is then dispersed by a spectrograph and recorded on an optical multichannel analyzer. A best-fit temperature is found in several seconds with the aid of a computer for each experimental spectrum by a least squares comparison with calculated spectra. The model used to calculate the theoretical spectra incorporates the temperature and pressure dependence of the pressure-broadened rotational Raman lines, includes the nonresonant background susceptibility, and assumes that the pump laser has a finite linewidth. Temperatures are fit to experimental spectra recorded over the temperature range of 135 to 296K, and over the pressure range of 0.13 to 15.3 atm. In addition to the spatially resolved single point work, we have used multipoint CARS to obtain information from many spatially resolved volume elements along a cylindrical line (0.1 x 0.1 x 2.0 mm). We also obtained qualitative information on the instantaneous species concentration and temperature at 20 spatially resolved volume elements (0.1 x 0.1 x 0.1 mm) along a line.

Snow, J. B.↗

Test code for the assessment and improvement of Reynolds stress models

An existing two-dimensional, compressible flow, Navier-Stokes computer code, containing a full Reynolds stress turbulence model, was adapted for use as a test bed for assessing and improving turbulence models based on turbulence simulation experiments. To date, the results of using the code in comparison with simulated channel flow and over an oscillating flat plate have shown that the turbulence model used in the code needs improvement for these flows. It is also shown that direct simulation of turbulent flows over a range of Reynolds numbers are needed to guide subsequent improvement of turbulence models.

Rubesin, M. W.↗

Navier-Stokes dynamics and aeroelastic computations for vortical flows, buffet and aeroelastic applications

The accomplishments achieved during the period include conference and proceedings publications, journal papers, and abstracts which are either published, accepted for publication or under review. Conference presentations and NASA highlight publications are also included. Two of the conference proceedings publications are attached along with a Ph.D. dissertation abstract and table of contents. In the first publication, computational simulation of three-dimensional flows around a delta wing undergoing rock and roll-divergence motions is presented. In the second publication, the unsteady Euler equations and the Euler equations of rigid body motion, both written in the moving frame of reference, are sequetially solved to simulate the limit-cycle rock motion of slender delta wings. In the dissertation abstract, unsteady flows around rigid or flexible delta wings with and without oscillating leading-edge flaps are considered.

Kandil, Osama A.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗

A Comparison of Three Navier-Stokes Solvers for Exhaust Nozzle Flowfields

A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer.

Georgiadis, Nicholas J.↗

Confined separated-swirling flows in diffusing ducts

The objective is to establish a detailed experimental data base for evaluation of Navier-Stokes codes for confined separated flows in diffusing s-ducts. The computational thrusts include the following: (1) extension and validation of the LeRC parabolized Navier-Stokes solver, PEPSIG, into the separated flow regime using 'flare' type approximations; (2) evaluation and extensions of state-of-the-art turbulence models for confined separated flow with and without swirl; and (3) evaluation and validation of LeRC time marching 3-D Navier-Stokes code, PROTEUS, into confined separate flow regime. Various aspects of the study are presented in viewgraph form.

Anderson, Bernhard H.↗

Thrust chamber performance using Navier-Stokes solution

The viscous, axisymmetric flow in the thrust chamber of the space shuttle main engine (SSME) was computed on the CRAY 205 computer using the general interpolants method (GIM) code. Results show that the Navier-Stokes codes can be used for these flows to study trends and viscous effects as well as determine flow patterns; but further research and development is needed before they can be used as production tools for nozzle performance calculations. The GIM formulation, numerical scheme, and computer code are described. The actual SSME nozzle computation showing grid points, flow contours, and flow parameter plots is discussed. The computer system and run times/costs are detailed.

Chan, J. S.↗