Search NASA⌕ Search

SEARCH · Search NASA

Results for “INCOMPRESSIBLE FLUID”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Toward Better Modeling of Supercritical Turbulent Mixing

study was done as part of an effort to develop computational models representing turbulent mixing under thermodynamic supercritical (here, high pressure) conditions. The question was whether the large-eddy simulation (LES) approach, developed previously for atmospheric-pressure compressible-perfect-gas and incompressible flows, can be extended to real-gas non-ideal (including supercritical) fluid mixtures. [In LES, the governing equations are approximated such that the flow field is spatially filtered and subgrid-scale (SGS) phenomena are represented by models.] The study included analyses of results from direct numerical simulation (DNS) of several such mixing layers based on the Navier-Stokes, total-energy, and conservation- of-chemical-species governing equations. Comparison of LES and DNS results revealed the need to augment the atmospheric- pressure LES equations with additional SGS momentum and energy terms. These new terms are the direct result of high-density-gradient-magnitude regions found in the DNS and observed experimentally under fully turbulent flow conditions. A model has been derived for the new term in the momentum equation and was found to perform well at small filter size but to deteriorate with increasing filter size. Several alternative models were derived for the new SGS term in the energy equation that would need further investigations to determine if they are too computationally intensive in LES.

Selle, Laurent↗
NASA NTRS1992

Low Mach number two-dimensional hydrodynamic turbulence - Energy budgets and density fluctuations in a polytropic fluid

Theory suggests that three distinct types of turbulence can occur in the low Mach number limit of polytropic flow: nearly incompressible flows dominated by vorticity, nearly pure acoustic turbulence dominated by compression, and flows characterized by near statistical equipartition of vorticity and compressions. Distinctions between these kinds of turbulence are investigated here by direct numerical simulation of two-dimensional compressible hydrodynamic turbulence. Dynamical scalings of density fluctuations, examination of the ratio of transverse to longitudinal velocity fluctuations, and spectral decomposition of the fluctuations are employed to distinguish the nature of these low Mach number solutions. A strong dependence on the initial data is observed, as well as a tendency for enhanced effects of compressibility at later times and at higher wave numbers, as suggested by theories of nearly incompressible flows.

Ghosh, S.↗
NASA NTRS1988

Domain decomposition algorithms and computation fluid dynamics

In the past several years, domain decomposition was a very popular topic, partly motivated by the potential of parallelization. While a large body of theory and algorithms were developed for model elliptic problems, they are only recently starting to be tested on realistic applications. The application of some of these methods to two model problems in computational fluid dynamics are investigated. Some examples are two dimensional convection-diffusion problems and the incompressible driven cavity flow problem. The construction and analysis of efficient preconditioners for the interface operator to be used in the iterative solution of the interface solution is described. For the convection-diffusion problems, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is discussed.

Chan, Tony F.↗
NASA NTRS1996

Thermal analysis of HGFQ using FIDAP(trademark): Solidification front motion

The High Gradient Furnace with Quench (HGFQ) is being designed by NASA/MSFC for flight on the International Space Station. The furnace is being designed specifically for solidification experiments in metal and metallic alloy systems. The HGFQ Product development Team (PDT) has been active since January 1994 and their effort is now in early Phase B. Thermal models have been developed both by NASA and Sverdrup (support contractor) to assist in the HGFQ design effort. Both these models use SINDA as a solution engine, but the NASA model was developed using PATRAN and includes more detail than the Sverdrup model. These models have been used to guide design decisions and have been validated through experimentation on a prototypical 'Breadboard' furnace at MSFC. One facet of the furnace operation of interest to the designers is the sensitivity of the solidification interface location to changes in the furnace setpoint. Specifically of interest is the motion (position and velocity) of the solidification front due to a small perturbation in the furnace temperature. FIDAP(TM) is a commercially available finite element program for analysis of heat transfer and fluid flow processes. Its strength is in solution of the Navier-Stokes equations for incompressible flow, but among its capabilities is the analysis of transient processes involving radiation and solidification. The models presently available from NASA and Sverdrup are steady-state models and are incapable of computing the motion of the solidification front. The objective of this investigation is to use FIDAP(TM) to compute the motion of the solidification interface due to a perturbation in the furnace setpoint.

Woodbury, Keith A.↗
NASA NTRS2004

Numerical Simulations of Vortex Shedding in Hydraulic Turbines

Turbomachines for rocket propulsion applications operate with many different working fluids and flow conditions. Oxidizer boost turbines often operate in liquid oxygen, resulting in an incompressible flow field. Vortex shedding from airfoils in this flow environment can have adverse effects on both turbine performance and durability. In this study the effects of vortex shedding in a low-pressure oxidizer turbine are investigated. Benchmark results are also presented for vortex shedding behind a circular cylinder. The predicted results are compared with available experimental data.

Dorney, Daniel↗
NASA NTRS2002

A Novel Approach for Modeling Chemical Reaction in Generalized Fluid System Simulation Program

The Generalized Fluid System Simulation Program (GFSSP) is a computer code developed at NASA Marshall Space Flight Center for analyzing steady state and transient flow rates, pressures, temperatures, and concentrations in a complex flow network. The code, which performs system level simulation, can handle compressible and incompressible flows as well as phase change and mixture thermodynamics. Thermodynamic and thermophysical property programs, GASP, WASP and GASPAK provide the necessary data for fluids such as helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, water, a hydrogen, isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, several refrigerants, nitrogen trifluoride and ammonia. The program which was developed out of need for an easy to use system level simulation tool for complex flow networks, has been used for the following purposes to name a few: Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump Secondary Flow Circuits, Axial Thrust Balance of the Fastrac Engine Turbopump, Pressurized Propellant Feed System for the Propulsion Test Article at Stennis Space Center, X-34 Main Propulsion System, X-33 Reaction Control System and Thermal Protection System, and International Space Station Environmental Control and Life Support System design. There has been an increasing demand for implementing a combustion simulation capability into GFSSP in order to increase its system level simulation capability of a liquid rocket propulsion system starting from the propellant tanks up to the thruster nozzle for spacecraft as well as launch vehicles. The present work was undertaken for addressing this need. The chemical equilibrium equations derived from the second law of thermodynamics and the energy conservation equation derived from the first law of thermodynamics are solved simultaneously by a Newton-Raphson method. The numerical scheme was implemented as a User Subroutine in GFSSP.

Sozen, Mehmet↗
DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. 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. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Conjugated heat transfer in a laminar boundary layer with heat source at the wall

This work presents the solution of the temperature field in a two-dimensional laminar incompressible flow over a conducting solid plate with a line heat source located at the fluid-solid interface perpendicular to the flow direction. A numerical scheme was used to obtain the temperature profiles as a function of the source strength, and of the properties of the fluid and the solid. The heat conduction and forced convection in the fluid and the heat conduction in the solid were solved for the case of moderate temperature rise, where the assumption of constant properties applies. The model enables the improvement of an instrument for the detection of boundary layer separation. It was found that for the actual parameters of the separation detector, a distance of 4 to 24 mm between the sensors gives an indication of 70 percent of the maximum temperature difference.

Brosh, A.↗
NASA NTRS1985

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modelled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗
NASA NTRS1985

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later Shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modeled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗
NASA NTRS1988

Theoretical studies of solar lasers and converters

The previously constructed one dimensional model for the simulated operation of an iodine laser assumed that the perfluoroalkyl iodide gas n-C3F7I was incompressible. The present study removes this simplifying assumption and considers n-C3F7I as a compressible fluid.

Heinbockel, John H.↗
NASA NTRS1987

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗
NASA NTRS1992

High order parallel numerical schemes for solving incompressible flows

The use of parallel computers for numerically solving flow fields has gained much importance in recent years. This paper introduces a new high order numerical scheme for computational fluid dynamics (CFD) specifically designed for parallel computational environments. A distributed MIMD system gives the flexibility of treating different elements of the governing equations with totally different numerical schemes in different regions of the flow field. The parallel decomposition of the governing operator to be solved is the primary parallel split. The primary parallel split was studied using a hypercube like architecture having clusters of shared memory processors at each node. The approach is demonstrated using examples of simple steady state incompressible flows. Future studies should investigate the secondary split because, depending on the numerical scheme that each of the processors applies and the nature of the flow in the specific subdomain, it may be possible for a processor to seek better, or higher order, schemes for its particular subcase.

Lin, Avi↗
NASA NTRS1996

Impulsive Start of a Symmetric Airfoil at High Angle of Attack

The fluid dynamic phenomena following the impulsive start of a NACA 0015 airfoil were studied by using a time accurate solution of the incompressible laminar Navier-Stokes equations. Angle of attack was set at 10 deg to simulate steady-state poststall conditions at a Reynolds number of 1.2 x 10(exp 4). The calculation revealed that large initial lift values can be obtained, immediately following the impulsive start, when a trapped vortex develops above the airfoil. Before the buildup of this trapped vortex and immediately after the airfoil was set into motion, the fluid is attached to the airfoil's surface and flows around the trailing edge, demonstrating the delay in the buildup of the classical Kutta condition. The transient of this effect is quite short and is followed by an attached How event that leads to the trapped vortex that has a longer duration. The just described initial phenomenon eventually transits into a fully developed separated flow pattern identifiable by an alternating, periodic vortex shedding.

Katz, Joseph↗
NASA NTRS1989

A theoretical study of fluid forces on a centrifugal impeller rotating and whirling in a vaned diffuser

The fluid forces on a centrifugal impeller rotating and whirling in a vaned diffuser are analyzed on the assumption that the number of impeller and diffuser vanes is so large that the flows are perfectly guided by the vanes. The flow is taken to be two dimensional, inviscid, and incompressible, but the effects of impeller and diffuser losses are taken into account. It is shown that the interaction with the vaned diffuser may cause destabilizing fluid forces. From these discussions, it is found that the whirling forces are closely related to the steady head-capacity characteristics of the impeller. This physical understanding of the whirling forces can be applied also to the cases with volute casings. At partial capacities, it is shown that the impeller forces change greatly when the flow rate and whirl velocity are near to the impeller or vaned diffuser attributed rotating stall onset capacity, and the stall propagation velocity, respectively. In such cases the impeller forces may become destabilizing for impeller whirl.

Tsujimoto, Yoshinobu↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗
NASA NTRS1985

Turbulent solutions of equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗