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 127 records · Page 7

A Parallel Incompressible Navier-Stokes Solver With a Parallel Multigrid Elliptic Kernal

The paper describes numerical algorithms and parallel implementations of a time-dependent, incompressible Navier-Stokes flow solver and a multigrid elliptic solver, which is also used as a computation kernal in the flow solver. The implemented solvers are numerically stable and computationally efficient, and they scale well to a large number of processors for problems with moderate granularity.

solver Navier-Stokes solver multigrid elliptic sol↗

Flow Simulation of N3-X Hybrid Wing-Body Configuration

System studies show that a N3-X hybrid wing-body aircraft with a turboelectric distributed propulsion system using a mail-slot inlet/nozzle nacelle can meet the environmental and performance goals for N+3 generation transports (three generations beyond the current air transport technology level) set by NASA s Subsonic Fixed Wing Project. In this study, a Navier-Stokes flow simulation of N3-X on hybrid unstructured meshes was conducted, including the mail-slot propulsor. The geometry of the mail-slot propulsor was generated by a CAD (Computer-Aided Design)-free shape parameterization. A body force approach was used for a more realistic and efficient simulation of the turning and loss effects of the fan blades and the inlet-fan interactions. Flow simulation results of the N3-X demonstrates the validity of the present approach.

Kim, Hyoungjin↗

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗

Hydrodynamically induced fluid transfer and non-convective double-diffusion in microgravity sliding solvent diffusion cells

Microgravity can provide a diffusion-dominated environment for double-diffusion and diffusion-reaction experiments otherwise disrupted by buoyant convection or sedimentation. In sliding solvent diffusion cells, a diffusion interface between two liquid columns is achieved by aligning two offset sliding wells. Fluid in contact with the sliding lid of the cavities is subjected to an applied shear stress. The momentum change by the start/stop action of the well creates an additional hydrodynamical force. In microgravity, these viscous and inertial forces are sufficiently large to deform the diffusion interface and induce hydrodynamic transfer between the wells. A series of KC-135 parabolic flight experiments were conducted to characterize these effects and establish baseline data for microgravity diffusion experiments. Flow visualizations show the diffusion interface to be deformed in a sinusoidal fashion following well alignment. After the wells were separated again in a second sliding movement, the total induced liquid transfer was determined and normalized by the well aspect ratio. The normalized transfer decreased linearly with Reynolds number from 3.3 to 4.0% (w/v) for Re = 0.4 (Stokes flow) to a minimum of 1.0% for Re = 23 to 30. Reynolds numbers that provide minimum induced transfers are characterized by an interface that is highly deformed and unsuitable for diffusion measurements. Flat diffusion interfaces acceptable for diffusion measurements are obtained with Reynolds numbers on the order of 7 to 10. Microgravity experiments aboard a sounding rocket flight verified counterdiffusion of different solutes to be diffusion dominated. Ground control experiments showed enhanced mixing by double-diffusive convection. Careful selection of experimental parameters improves initial conditions and minimizes induced transfer rates.

Pollmann, Konrad W.↗

TPS Sizing for Access-to-Space Vehicles

A study was carried out to identify, develop, and benchmark simulation techniques needed for optimum thermal protection system (TPS) material selection and sizing for reusable launch vehicles. Fully viscous, chemically reacting, Navier-Stokes flow solutions over the Langley wing-body single stage to orbit (SSTO) configuration were generated and coupled with an in-depth conduction code. Results from the study provide detailed TPS heat shield materials selection and thickness sizing for the wing-body SSTO. These results are the first ever achieved through the use of a complete, trajectory based hypersonic, Navier-Stokes solution database. TPS designs were obtained for both laminar and turbulent entry trajectories using the Access-to-Space baseline materials such as tailorable advanced blanket insulation. The TPS design effects (materials selection and thickness) of coupling material characteristics to the aerothermal environment are illustrated. Finally, a sample validation case using the shuttle flight database is included.

Henline, William↗

Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows

Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.

Pandya, Mohagna J.↗

An analysis of the flow field in the region of the ASRM field joints

The flow field in the region of a solid rocket motor field joint is very important since fluid dynamic and mechanical propellant stresses can couple to cause a motor failure at a joint. Presented here is an examination of the flow field in the region of the Advanced Solid Rocket Motor (ASRM) field joints. The analyses were performed as a first step in assessing the design of the ASRM forward and aft field joints in order to assure the proper operation of the motor prior to further development of test firing. The analyses presented here were performed by employing a two-dimensional axisymmetric assumption. Fluent/BFC, a three dimensional full Navier-Stokes flow field code, was used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical algorithm. Wall functions are used to determine the character of the laminar sublayer, and a standard kappa-epsilon turbulence model is used to close the fluid dynamic equations. The analyses performed to this date verify that the ASRM field joint design operates properly. The fluid dynamic stresses at the field joints are small due to the inherent design of the field joints. A problem observed in some other solid rocket motors is that large fluid dynamic stresses are generated at the motor joint on the downstream propellant grain due to forward facing step geometries. The design of the ASRM field joints are such that this is not a problem as shown by the analyses. Also, the analyses of the inhibitor stub left protruding into the port flow from normal propellant burn back show that more information is necessary to complete these analyses. These analyses were performed as parametric analyses in relation to the height of the inhibitor stub left protruding into the motor port. A better estimate of the amount of the inhibitor stub remaining at later burn times must be determined since the height which the inhibitor stub protrudes into the port flow drastically affects the fluid dynamics induced stresses on the propellant grain at the field joints.

Dill, Richard A.↗

Analysis and Inverse Design of the HSR Arrow Wing Configuration with Fuselage, Wing, and Flow Through Nacelles

The design process for developing the natural flow wing design on the HSR arrow wing configuration utilized several design tools and analysis methods. Initial fuselage/wing designs were generated with inviscid analysis and optimization methods in conjunction with the natural flow wing design philosophy. A number of designs were generated, satisfying different system constraints. Of the three natural flow wing designs developed, the NFWAc2 configuration is the design which satisfies the constraints utilized by McDonnell Douglas Aerospace (MDA) in developing a series of optimized configurations; a wind tunnel model of the MDA designed OPT5 configuration was constructed and tested. The present paper is concerned with the viscous analysis and inverse design of the arrow wing configurations, including the effects of the installed diverters/nacelles. Analyses were conducted with OVERFLOW, a Navier-Stokes flow solver for overset grids. Inverse designs were conducted with OVERDISC, which couples OVERFLOW with the CDISC inverse design method. An initial system of overset grids was generated for the OPT5 configuration with installed diverters/nacelles. An automated regridding process was then developed to use the OPT5 component grids to create grids for the natural flow wing designs. The inverse design process was initiated using the NFWAc2 configuration as a starting point, eventually culminating in the NFWAc4 design-for which a wind tunnel model was constructed. Due to the time constraints on the design effort, initial analyses and designs were conducted with a fairly coarse grid; subsequent analyses have been conducted on a refined system of grids. Comparisons of the computational results to experiment are provided at the end of this paper.

Krist, Steven E.↗

An Iron-Rain Model for Core Formation on Asteroid 4 Vesta

Asteroid 4 Vesta is differentiated into a crust, mantle, and core, as demonstrated by studies of the eucrite and diogenite meteorites and by data from NASA's Dawn spacecraft. Most models for the differentiation and thermal evolution of Vesta assume that the metal phase completely melts within 20 degrees of the eutectic temperature, well before the onset of silicate melting. In such a model, core formation initially happens by Darcy flow, but this is an inefficient process for liquid metal and solid silicate. However, the likely chemical composition of Vesta, similar to H chondrites with perhaps some CM or CV chondrite, has 13-16 weight percent S. For such compositions, metal-sulfide melting will not be complete until a temperature of at least 1350 degrees Centigrade. The silicate solidus for Vesta's composition is between 1100 and 1150 degrees Centigrade, and thus metal and silicate melting must have substantially overlapped in time on Vesta. In this chemically and physically more likely view of Vesta's evolution, metal sulfide drops will sink by Stokes flow through the partially molten silicate magma ocean in a process that can be envisioned as "iron rain". Measurements of eucrites show that moderately siderophile elements such as Ni, Mo, and W reached chemical equilibrium between the metal and silicate phases, which is an important test for any Vesta differentiation model. The equilibration time is a function of the initial metal grain size, which we take to be 25-45 microns based on recent measurements of H6 chondrites. For these sizes and reasonable silicate magma viscosities, equilibration occurs after a fall distance of just a few meters through the magma ocean. Although metal drops may grow in size by merger with other drops, which increases their settling velocities and decreases the total core formation time, the short equilibration distance ensures that the moderately siderophile elements will reach chemical equilibrium between metal and silicate before metal drop merger becomes important. In this model, there must be at least 30 percent melting of the silicate phase when metal melting is complete, corresponding to a crust thickness of at least 30 kilometers on Vesta, consistent with Dawn gravity observations. Greater degrees of silicate melting and a correspondingly thicker crust are possible if Vesta accreted sufficiently rapidly.

Kiefer, Walter S.↗

Parametrics on 2D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a two-dimensional (2D) Navier-Stokes flow solver. Parametric studies were performed on turbulence models, computational grids and bleed models. The computer flowfield was substantially different from the original inviscid design, due to interactions of shocks, boundary layers, and bleed. Good agreement with experimental data was obtained in many aspects. Many of the discrepancies were thought to originate primarily from 3D effects. Therefore, a balance should be struck between expending resources on a high fidelity 2D simulation, and the inherent limitations of 2D analysis. The solutions were fairly insensitive to turbulence models, grids and bleed models. Overall, the k-e turbulence model, and the bleed models based on unchoked bleed hole discharge coefficients or uniform velocity are recommended. The 2D Navier-Stokes methods appear to be a useful tool for the design and analysis of supersonic inlets, by providing a higher fidelity simulation of the inlet flowfield than inviscid methods, in a reasonable turnaround time.

Mizukami, M.↗

Solution of three-dimensional afterbody flow using reduced Navier-Stokes equations

The flow over afterbody geometries was investigated using the reduced Navier-Stokes (RNS) approximation. Both pressure velocity flux-split and composites velocity primitive variable formulations were considered. Pressure or pseudopotential relaxation procedures are combined with sparse matrix or coupled strongly implicit algorithms to form a three-dimensional solver for general non-orthogonal coordinates. Three-dimensional subsonic and transonic viscous/inviscid interacting flows were evaluated. Solutions with and without regions of recirculation were obtained.

Almahroos, H. M. H.↗

Progress in incompressible Navier-Stokes computations for propulsion flows and its dual-use applications

Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.

Kiris, Cetin↗

Approximate solutions to the Navier-Stokes initial value problem

A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.

Rosen, G.↗

Broadness, decay, and correlation functions of isotropic homogeneous turbulence

General theoretical relationships consistent with experiments are obtained for isotropic homogeneous incompressible-fluid turbulence governed by the Navier-Stokes flow equation. A key quantity that structures the turbulence is the 'broadness' B of the probability measure over velocity fields. Both the decay law and the longitudinal correlation function for small r appear as simple functions of this quasiconstant parameter B.

Rosen, G.↗

Analytical estimates of turbulent MHD transport coefficients

Turbulent transfer rates from small-scale MHD excitations to large-scale Fourier modes are calculated algebraically, using the method of Biskamp and Welter. Three cases are considered: two-dimensional Navier-Stokes flows, two-dimensional incompressible MHD, and the weakly three-dimensional Strauss equations. In all cases, an initially large spectral gap between the small-scale and large-scale excitations is assumed, and attention focusses on the initial values of the back-transfer rates. The sign of the transfer is determined by the sign of an analytically calculable eddy viscosity and/or anomalous resistivity. We are able to confirm the results of Biskamp and Welter for the case of two-dimensional MHD, but find some differences for the case of the Strauss equations. It is argued that the Strauss equations may not exhibit an inverse cascade phenomenon for the spatially periodic case unless their initial spectra are such that the behavior is essentially that of two-dimensional MHD.

Montgomery, D.↗

Navier-Stokes computations for exotic airfoils

An efficient hyperbolic grid generator with improvements for handling sharp corners and concave surfaces is combined with an efficient and accurate Navier-Stokes flow solver. This combination is applied to some rather complex two-dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 degrees are presented and compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Innovative graphical techniques for both the static and unsteady display of the flow fields are presented and discussed.

Barth, T. J.↗

Navier-Stokes computations for exotic airfoils

A hyperbolic grid generator with improvements for handling sharp corners and concave surfaces combined with a Navier-Stokes flow solver is applied to complex two dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 deg are compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Graphical techniques for the static and unsteady display of the flow fields are discussed.

Barth, T. J.↗