Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solving 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 271 records · Page 15

CryoTran user's manual, version 1.0

The development of cryogenic fluid management systems for space operation is a major portion of the efforts of the Cryogenic Fluids Technology Office (CFTO) at the NASA Lewis Research Center. Analytical models are a necessary part of experimental programs which are used to verify the results of experiments and are also used as a predictor for parametric studies. The CryoTran computer program is a bridge to obtain analytical results. The object of CryoTran is to coordinate these separate analyses into an integrated framework with a user-friendly interface and a common cryogenic property database. CryoTran is an integrated software system designed to help solve a diverse set of problems involving cryogenic fluid storage and transfer in both ground and low-g environments.

Cowgill, Glenn R.↗

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity have been included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary and the subadjacent fluid motion; it also features Gauss' method for solving the nonlinear inverse problem associated with steady motional induction. The trade-off between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and fluid flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Goerte V.↗

Three-Fluid Magnetohydrodynamic Modeling of the Solar Wind in the Outer Heliosphere

We have developed a three-fluid, fully three-dimensional magnetohydrodynamic model of the solar wind plasma in the outer heliosphere as a co-moving system of solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Our approach takes into account the effects of electron heat conduction and dissipation of Alfvenic turbulence on the spatial evolution of the solar wind plasma and interplanetary magnetic fields. The turbulence transport model is based on the Reynolds decomposition of physical variables into mean and fluctuating components and uses the turbulent phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. We solve the coupled set of the three-fluid equations for the mean-field solar wind and the turbulence equations for the turbulence energy, cross helicity, and correlation length. The equations are written in the rotating frame of reference and include heating by turbulent dissipation, energy transfer from interstellar pickup protons to solar wind protons, and solar wind deceleration due to the interaction with the interstellar hydrogen. The numerical solution is constructed by the time relaxation method in the region from 0.3 to 100 AU. Initial results from the novel model are presented.

Usmanov, Arcadi V.↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Fluid-plasticity of thin cylindrical shells

Dynamic plastic response of a thin cylindrical shell, immersed in a potential fluid initially at rest and subjected to internal pressure pulse of arbitrary shape and duration, is examined. The shell is assumed to respond as a rigid-perfectly plastic material while the fluid is taken as inviscid and incompressible. The fluid back pressure is incorporated into the equation of motion of the shell as an added mass term. Since arbitrary pulses can be reduced to equivalent rectangular pulses, the equation of motion is solved only for a rectangular pulse. The influence of the fluid in reducing the final plastic deformation is demonstrated by a numerical example.

Krajcinovic, D.↗

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗

A Multiblock Approach for Calculating Incompressible Fluid Flows on Unstructured Grids

A multiblock approach is presented for solving two-dimensional incompressible turbulent flows on unstructured grids. The artificial compressibility form of the governing equations is solved by a vertex-centered, finite-volume implicit scheme which uses a backward Euler time discretization. Point Gauss-Seidel relaxations are used to solve the linear system of equations at each time step. This work introduces a multiblock strategy to the solution procedure, which greatly improves the efficiency of the algorithm by significantly reducing the memory requirements while not increasing the CPU time. Results presented in this work shows that the current multiblock algorithm requires 70% less memory than the single block algorithm.

Sheng, Chunhua↗

On the motion of compact objects in relativistic viscous fluids

We present a world-line effective field theory of compact objects moving relativistically through a viscous fluid. The theory is valid when velocity gradients are small compared to the inverse size of the object. Working within the EFT eliminates the need to solve a boundary value problem by turning all interactions between the fluid and the object into a source term in the action. We use the EFT to derive the relativistic equations of motion for a compact object immersed in a viscous fluid in a curved background, when the relative velocity of the object and the fluid is small compared to the speed of light.

astrophysical black holes↗

Software Performs Complex Design Analysis

Designers use computational fluid dynamics (CFD) to gain greater understanding of the fluid flow phenomena involved in components being designed. They also use finite element analysis (FEA) as a tool to help gain greater understanding of the structural response of components to loads, stresses and strains, and the prediction of failure modes. Automated CFD and FEA engineering design has centered on shape optimization, which has been hindered by two major problems: 1) inadequate shape parameterization algorithms, and 2) inadequate algorithms for CFD and FEA grid modification. Working with software engineers at Stennis Space Center, a NASA commercial partner, Optimal Solutions Software LLC, was able to utilize its revolutionary, one-of-a-kind arbitrary shape deformation (ASD) capability-a major advancement in solving these two aforementioned problems-to optimize the shapes of complex pipe components that transport highly sensitive fluids. The ASD technology solves the problem of inadequate shape parameterization algorithms by allowing the CFD designers to freely create their own shape parameters, therefore eliminating the restriction of only being able to use the computer-aided design (CAD) parameters. The problem of inadequate algorithms for CFD grid modification is solved by the fact that the new software performs a smooth volumetric deformation. This eliminates the extremely costly process of having to remesh the grid for every shape change desired. The program can perform a design change in a markedly reduced amount of time, a process that would traditionally involve the designer returning to the CAD model to reshape and then remesh the shapes, something that has been known to take hours, days-even weeks or months-depending upon the size of the model.

Source record↗

Electromagnetic ion beam instabilities - Comparison of oneand two-dimensional simulations

Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principal features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.

Winske, D.↗

Analytical modeling of helium II in forced flow conditions

A numerical model was developed to examine the flow behavior of He II under forced flow conditions. The model basically solves the governing equations in the two-fluid model of He II, and can be used to assess the pump requirement for transferring He II in spaceborne systems. Predictions of liquid temperature and pressure drop agree well with empirical data obtained from small scale tests. At high transfer rates, the model predicts that the forced convection dominates the Gorter-Mellink conduction in heat transfer, and the classical turbulent flow correlation is sufficient for calculating the pressure drop across the transfer line.

Ng, Y. S.↗

An analytic model of axisymmetric mantle plume due to thermal and chemical diffusion

An analytic model of axisymmetric mantle plumes driven by either thermal diffusion or combined diffusion of both heat and chemical species from a point source is presented. The governing equations are solved numerically in cylindrical coordinates for a Newtonian fluid with constant viscosity. Instead of starting from an assumed plume source, constraints on the source parameters, such as the depth of the source regions and the total heat input from the plume sources, are deduced using the geophysical characteristics of mantle plumes inferred from modelling of hotspot swells. The Hawaiian hotspot and the Bermuda hotspot are used as examples. Narrow mantle plumes are expected for likely mantle viscosities. The temperature anomaly and the size of thermal plumes underneath the lithosphere can be sensitive indicators of plume depth. The Hawaiian plume is likely to originate at a much greater depth than the Bermuda plume. One suggestive result puts the Hawaiian plume source at a depth near the core-mantle boundary and the source of the Bermuda plume in the upper mantle, close to the 700 km discontinuity. The total thermal energy input by the source region to the Hawaiian plume is about 5 x 10(10) watts. The corresponding diameter of the source region is about 100 to 150 km. Chemical diffusion from the same source does not affect the thermal structure of the plume.

Liu, Mian↗

A spectral multi-domain technique applied to high-speed chemically reacting flows

The first applications of a spectral multidomain method for viscous compressible flow is presented. The method imposes a global flux balance condition at the interface so that high-order continuity of the solution is preserved. The global flux balance is imposed in terms of a spectral integral of the discrete equations across adjoining domains. Since the discretized equations interior to each domain solved are uncoupled from each other, and since the interface relation has a block structure, the solution scheme can be adapted to the particular requirements of each subdomain. The spectral multidomain technique presented is well-suited for the multiple scales associated with the chemically reacting and transition flows in hypersonic research. A nonstaggered multidomain discretization is used for the chemically reacting flow calculation, and the first implementation of a staggered multidomain mesh is presented for accurately solving the stability equation for a viscous compressible fluid.

Macaraeg, Michele G.↗

Experimental and analytical studies of flow through a ventral and axial exhaust nozzle system for STOVL aircraft

Flow through a combined ventral and axial exhaust nozzle system was studied experimentally and analytically. The work is part of an ongoing propulsion technology effort at NASA Lewis Research Center for short takeoff, vertical landing (STOVL) aircraft. The experimental investigation was done on the NASA Lewis Powered Lift Facility. The experiment consisted of performance testing over a range of tailpipe pressure ratios from 1 to 3.2 and flow visualization. The analytical investigation consisted of modeling the same configuration and solving for the flow using the PARC3D computational fluid dynamics program. The comparison of experimental and analytical results was very good. The ventral nozzle performance coefficients obtained from both the experimental and analytical studies agreed within 1.2 percent. The net horizontal thrust of the nozzle system contained a significant reverse thrust component created by the flow overturning in the ventral duct. This component resulted in a low net horizontal thrust coefficient. The experimental and analytical studies showed very good agreement in the internal flow patterns.

Esker, Barbara S.↗

Experimental and analytical studies of flow through a ventral and axial exhaust nozzle system for STOVL aircraft

Flow through a combined ventral and axial exhaust nozzle system was studied experimentally and analytically. The work is part of an ongoing propulsion technology effort at NASA Lewis Research Center for short takeoff, vertical landing (STOVL) aircraft. The experimental investigation was done on the NASA Lewis Powered Lift Facility. The experiment consisted of performance testing over a range of tailpipe pressure ratios from 1 to 3.2 and flow visualization. The analytical investigation consisted of modeling the same configuration and solving for the flow using the PARC3D computational fluid dynamics program. The comparison of experimental and analytical results was very good. The ventral nozzle performance coefficients obtained from both the experimental and analytical studies agreed within 1.2 percent. The net horizontal thrust of the nozzle system contained a significant reverse thrust component created by the flow overturning in the ventral duct. This component resulted in a low net horizontal thrust coefficient. The experimental and analytical studies showed very good agreement in the internal flow patterns.

Esker, Barbara S.↗

Computerized reduction of elementary reaction sets for combustion modeling

If the entire set of elementary reactions is to be solved in the modeling of chemistry in computational fluid dynamics, a set of stiff ordinary differential equations must be integrated. Some of the reactions take place at very high rates, requiring short time steps, while others take place more slowly and make little progress in the short time step integration. The goal is to develop a procedure to automatically obtain sets of finite rate equations, consistent with a partial equilibrium assumptions, from an elementary set appropriate to local conditions. The possibility of computerized reaction reduction was demonstrated. However, the ability to use the reduced reaction set depends on the ability of the CFD approach in incorporate partial equilibrium calculations into the computer code. Therefore, the results should be tested on a code with partial equilibrium capability.

Wikstrom, Carl V.↗

High speed propeller acoustics and aerodynamics - A boundary element approach

The Boundary Element Method (BEM) is applied in this paper to the problems of acoustics and aerodynamics of high speed propellers. The underlying theory is described based on the linearized Ffowcs Williams-Hawkings equation. The surface pressure on the blade is assumed unknown in the aerodynamic problem. It is obtained by solving a singular integral equation. The acoustic problem is then solved by moving the field point inside the fluid medium and evaluating some surface and line integrals. Thus the BEM provides a powerful technique in calculation of high speed propeller aerodynamics and acoustics.

Farassat, F.↗

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity were included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary (CMB) and the sub-adjacent fluid motion; it also features Gauss' method for solving the non-linear inverse problem associated with steady motional induction. The tradeoff between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models (DGRF's) is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition, a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Coerte V.↗