Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Active motion of multiphase oil droplets: emergent dynamics of squirmers with evolving internal structure

Synthetic soft matter systems, when driven beyond equilibrium by active processes, offer the potential to achieve dynamical states and functions of a complexity found in living matter. Emulsions offer the basis of a simple yet versatile system for identification of the physicochemical principles underlying active soft matter, but how multiple internal phases within emulsion droplets (e.g., Janus morphologies) organize to impact emergent dynamics is not understood. Here, we create multiphase oil droplets with ultralow interfacial tensions but distinct viscosities, and drive them into motion in aqueous micellar solutions. Preferential solubilization of select components of the oil both drives the droplet motion and yields a progression of internal phase morphological states with distinct symmetries. In this work, we find the active droplets to exhibit five dynamical states during morphogenesis. By quantifying microscopic flow fields, we show that it is possible to map the diverse droplet behaviors to squirmer models of spherical microswimmers in Stokes flow, thus showing that multiphase droplets offer the basis of a versatile platform with which to study and engineer the hydrodynamics of microswimmers.

36 MATERIALS SCIENCE↗
NASA NTRS1995

An unstructured-grid software system for solving complex aerodynamic problems

A coordinated effort has been underway over the past four years to elevate unstructured-grid methodology to a mature level. The goal of this endeavor is to provide a validated capability to non-expert users for performing rapid aerodynamic analysis and design of complex configurations. The Euler component of the system is well developed, and is impacting a broad spectrum of engineering needs with capabilities such as rapid grid generation and inviscid flow analysis, inverse design, interactive boundary layers, and propulsion effects. Progress is also being made in the more tenuous Navier-Stokes component of the system. A robust grid generator is under development for constructing quality thin-layer tetrahedral grids, along with a companion Navier-Stokes flow solver. This paper presents an overview of this effort, along with a perspective on the present and future status of the methodology.

Frink, Neal T.↗
NASA NTRS1992

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗
NASA NTRS1997

Flow of a Gas Turbine Engine Low-Pressure Subsystem Simulated

The NASA Lewis Research Center is managing a task to numerically simulate overnight, on a parallel computing testbed, the aerodynamic flow in the complete low-pressure subsystem (LPS) of a gas turbine engine. The model solves the three-dimensional Navier- Stokes flow equations through all the components within the LPS, as well as the external flow around the engine nacelle. The LPS modeling task is being performed by Allison Engine Company under the Small Engine Technology contract. The large computer simulation was evaluated on networked computer systems using 8, 16, and 32 processors, with the parallel computing efficiency reaching 75 percent when 16 processors were used.

Veres, Joseph P.↗
NASA NTRS2006

High Fidelity Simulations of Unsteady Flow through Turbopumps and Flowliners

High fidelity computations were carried out to analyze the orbiter LH2 feedline flowliner. Computations were performed on the Columbia platform which is a 10,240-processor supercluster consisting of 20 Altix nodes with 512 processor each. Various computational models were used to characterize the unsteady flow features in the turbopump, including the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer, the orbiter manifold and a test article used to represent the manifold. Unsteady flow originating from the orbiter LPFTP inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the gimbal flowliners just upstream of the LPFTP. The flow fields for the orbiter manifold and representative test article are computed and analyzed for similarities and differences. The incompressible Navier-Stokes flow solver INS3D, based on the artificial compressibility method, was used to compute the flow of liquid hydrogen in each test article.

Kiris, Cetin C.↗
NASA NTRS2005

High Fidelity Simulations for Unsteady Flow Through the Orbiter LH2 Feedline Flowliner

High fidelity computations were carried out to analyze the orbiter M2 feedline flowliner. Various computational models were used to characterize the unsteady flow features in the turbopump, including the orbiter Low-Pressure-Fuel-Turbopump (LPFTP) inducer, the orbiter manifold and a test article used to represent the manifold. Unsteady flow originating from the orbiter LPFTP inducer is one of the major contributors to the high frequency cyclic loading that results in high cycle fatigue damage to the gimbal flowliners just upstream of the LPFTP. The flow fields for the orbiter manifold and representative test article are computed and analyzed for similarities and differences. An incompressible Navier-Stokes flow solver INS3D, based on the artificial compressibility method, was used to compute the flow of liquid hydrogen in each test article.

Kiris, Cetin C.↗
NASA NTRS1999

Implementation and Validation of the Chien k-epsilon Turbulence Model in the Wind Navier-Stokes Code

The two equation k-epsilon turbulence model of Chien has been implemented in the WIND Navier-Stokes flow solver. Details of the numerical solution algorithm, initialization procedure, and stability enhancements are described. Results obtained with this version of the model are compared with those from the Chien k-epsilon model in the NPARC Navier-Stokes code and from the WIND SST model for three validation cases: the incompressible flow over a smooth flat plate, the incompressible flow over a backward facing step, and the shock-induced flow separation inside a transonic diffuser. The k-epsilon model results indicate that the WIND model functions very similarly to that in NPARC, though the WIND code appears to he slightly more accurate in the treatment of the near-wall region. Comparisons of the k-epsilon model results with those from the SST model were less definitive, as each model exhibited strengths and weaknesses for each particular case.

Yoder, Dennis A.↗
NASA NTRS2017

Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow

To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54 turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJkg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilizing the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly under predicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJkg condition, but with up to a factor of three over prediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJkg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.

Earth entry↗
NASA NTRS2017

Experimental and Numerical Investigation of Air Radiation in Superorbital Expanding Flow

To investigate air radiation in expanding flows and provide experimental data for validating associated computational models, experiments were conducted in the X2 expansion tunnel facility at the Centre for Hypersonics of the University of Queensland. A 54deg turning angle wedge model was employed to generate steady expanding flows with in flow total enthalpies of 50.7, 63.4 and 75.4 MJ/kg. VUV spectra from 118 to 180 nm were acquired across the wedge at three equispaced distances away from the top of the model, as well as through its top surface. High speed filtered images were also obtained by coupling a Shimadzu 1 MHz high speed camera to a bandpass filter to obtain calibrated images of the 777 nm oxygen triplet. Both the across-wedge VUV spectra and filtered images of the 777 nm atomic oxygen were compared with NEQAIR simulations, which were performed using flow field data from two-dimensional CFD simulations with two-temperature 11-species air chemistry utilising the in-house Navier-Stokes flow solver Eilmer3. Data extracted from consecutive frames of the filtered high speed images confirmed up to 8 s of available test time for the flow conditions tested. For the strongly radiating 149 and 174 nm atomic nitrogen lines, large disagreement between experimental data and NEQAIR predictions can be observed from the start of the expansion fan where the electron-ion recombination process commences. The spatial extent, or spans of the radiance profiles of the 149 and 174 nm N lines are significantly underpredicted by NEQAIR, and are very close to those of N, N+ and electron number density profiles, which follow that of flow density. The electron-ion recombination process is proposed as the main reason for these discrepancies. The comparisons between NEQAIR simulations and filtered images of the 777 nm oxygen triplet show good agreement in the post-shock compression region and the start of the expansion fan for the 63.4 MJ/kg condition, but with up to a factor of three overprediction by NEQAIR further downstream, which is attributed to electron-impact excitation. Similar trends are found with the 75.4 MJ/kg condition, with reduced level of agreement in the compression region, which can be due to uncertainties in inflow condition.

expansion tube↗
NASA NTRS1995

The 3D 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 three-dimensional (3D) Navier-Stokes flow solver. A two-equation turbulence model, and a porous bleed model based on unchoked bleed hole discharge coefficients were used. Comparisons were made with experimental data, inviscid theory, and two-dimensional Navier-Stokes analyses. The main objective was to gain insight into the inlet fluid dynamics. Examination of the computational results along with the experimental data suggest that the cowl shock-sidewall boundary layer interaction near the leading edge caused a substantial separation in the wind tunnel inlet model. As a result, the inlet performance may have been compromised by increased spillage and higher bleed mass flow requirements. The internal flow contained substantial waves that were not in the original inviscid design. 3D effects were fairly minor for this inlet at on-design conditions. Navier-Stokes analysis appears to be an useful tool for gaining insight into the inlet fluid dynamics. It provides a higher fidelity simulation of the flowfield than the original inviscid design, by taking into account boundary layers, porous bleed, and their interactions with shock waves.

Mizukami, M.↗
NASA NTRS2005

Heart Pump Design for Cleveland Clinic Foundation

Through a Lewis CommTech Program project with the Cleveland Clinic Foundation, the NASA Lewis Research Center is playing a key role in the design and development of a permanently implantable, artificial heart pump assist device. Known as the Innovative Ventricular Assist System (IVAS), this device will take on the pumping role of the damaged left ventricle of the heart. The key part of the IVAS is a nonpulsatile (continuous flow) artificial heart pump with centrifugal impeller blades, driven by an electric motor. Lewis is part of an industry and academia team, led by the Ohio Aerospace Institute (OAI), that is working with the Cleveland Clinic Foundation to make IVAS a reality. This device has the potential to save tens of thousands of lives each year, since 80 percent of heart attack victims suffer irreversible damage to the left ventricle, the part of the heart that does most of the pumping. Impeller blade design codes and flow-modeling analytical codes will be used in the project. These codes were developed at Lewis for the aerospace industry but will be applicable to the IVAS design project. The analytical codes, which currently simulate the flow through the compressor and pump systems, will be used to simulate the flow within the blood pump in the artificial heart assist device. The Interdisciplinary Technology Office heads up Lewis' efforts in the IVAS project. With the aid of numerical modeling, the blood pump will address many design issues, including some fluid-dynamic design considerations that are unique to the properties of blood. Some of the issues that will be addressed in the design process include hemolysis, deposition, recirculation, pump efficiency, rotor thrust balance, and bearing lubrication. Optimum pumping system performance will be achieved by modeling all the interactions between the pump components. The interactions can be multidisciplinary and, therefore, are influenced not only by the fluid dynamics of adjacent components but also by thermal and structural effects. Lewis-developed flow-modeling codes to be used in the pump simulations will include a one-dimensional code and an incompressible three-dimensional Navier-Stokes flow code. These codes will analyze the prototype pump designed by the Cleveland Clinic Foundation. With an improved understanding of the flow phenomena within the prototype pump, design changes to improve the performance of the pump system can be verified by computer prior to fabrication in order to reduce risks. The use of Lewis flow modeling codes during the design and development process will improve pump system performance and reduce the number of prototypes built in the development phase. The first phase of the IVAS project is to fully develop the prototype in a laboratory environment that uses a water/glycerin mixture as the surrogate fluid to simulate blood. A later phase of the project will include testing in animals for final validation. Lewis will be involved in the IVAS project for 3 to 5 years.

Source record↗
NASA NTRS1992

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and liquid oxygen (LOX) flow distributions through the LOX posts have great effect on the combustion phenomenon inside the main combustion chamber. An advanced computational fluid dynamics (CFD) analysis will help to provide more accurate and efficient characterization of this type of flow field. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time and still maintain good accuracy in describing the flow features, the concept of porosity is employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. A validated non-isotropic porosity model is developed and incorporated into an existing Navier-Stokes flow solver (FDNS). Volume and surface porosity parameters, which are based on the configurations of local LOX post clustering, are introduced into the governing equations.

Cheng, Gary C.↗
NASA NTRS1992

On the drag of model dendrite fragments at low Reynolds number

An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with the primary dendrite arm aspect ratio and a second dimensionless shape parameter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.

Zakhem, R.↗
NASA NTRS1993

On the drag of model dendrite fragments at low Reynolds number

An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with primary dendrite arm aspect ratio and a second dimensionless shape paremeter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.

Zakhem, R.↗
NASA NTRS1993

Heat transfer in rocket combustion chambers

Complexities of liquid rocket engine heat transfer which involve the injector faceplate and film cooled walls are being investigated by computational analysis. A conjugate heat transfer analysis was used to describe localized heating phenomena associated with particular injector configurations and film coolant flows. These components were analyzed, and the analyses verified when appropriate test data were available. The component analyses are being synthesized into an overall flowfield/heat transfer model. A Navier-Stokes flow solver, the FDNS code, was used to make the analyses. Particular attention was given to the representation of the thermodynamic properties of the fluid streams. Unit flow models of specific coaxial injector elements have been developed and are being used to describe the flame structure near the injector faceplate.

Anderson, P.↗
NASA NTRS2013

Lessons Learned in the Selection and Development of Test Cases for the Aeroelastic Prediction Workshop: Rectangular Supercritical Wing

The Aeroelastic Prediction Workshop brought together an international community of computational fluid dynamicists as a step in defining the state of the art in computational aeroelasticity. The Rectangular Supercritical Wing (RSW) was chosen as the first configuration to study due to its geometric simplicity, perceived simple flow field at transonic conditions and availability of an experimental data set containing forced oscillation response data. Six teams performed analyses of the RSW; they used Reynolds-Averaged Navier-Stokes flow solvers exercised assuming that the wing had a rigid structure. Both steady-state and forced oscillation computations were performed by each team. The results of these calculations were compared with each other and with the experimental data. The steady-state results from the computations capture many of the flow features of a classical supercritical airfoil pressure distribution. The most dominant feature of the oscillatory results is the upper surface shock dynamics. Substantial variations were observed among the computational solutions as well as differences relative to the experimental data. Contributing issues to these differences include substantial wind tunnel wall effects and diverse choices in the analysis parameters.

Heeg, Jennifer↗