Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density 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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Numerical solution of the compressible Navier-Stokes equations using density gradients as additional dependent variables

Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.

Kwon, J. H.↗
NASA NTRS2004

A Parametric Study of Jet Interactions with Rarefied Flow

Three-dimensional computational techniques, in particular the uncoupled CFD-DSMC of the present study, are available to be applied to problems such as jet interactions with variable density regions ranging from a continuum jet to a rarefied free stream. When the value of the jet to free stream momentum flux ratio approximately greater than 2000 for a sharp leading edge flat plate forward separation vortices induced by the jet interaction are present near the surface. Also as the free stream number density n (infinity) decreases, the extent and magnitude of normalized pressure increases and moves upstream of the nozzle exit. Thus for the flat plate model the effect of decreasing n (infinity) is to change the sign of the moment caused by the jet interaction on the flat plate surface.

Glass, C. E.↗
NASA NTRS1993

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗
NASA NTRS1985

Applications of variational principles in computing rotational flows

Ecer and Akay (1983) have developed a variational formulation of rotational flow for Euler equations. The present paper provides a summary of these developments. The considered variational formulation provides a transformation of a type considered by Clebsch (1859). In this transformation, a new set of variables replaces the more commonly used primitive variables u(i), rho and p. Here, u(i) denotes the velocity components, while rho is the density, and p the pressure. The employed transformation produces a natural uncoupling of the equations when written in a quasi-linear form. After obtaining the governing equations in terms of the 'Clebsch variables', a solution scheme developed for calculating steady flows is discussed. Attention is given to numerical solutions of Euler equations based on the derived variational principles, and a study of inviscid, separated flows is conducted.

Ecer, A.↗
NASA NTRS1983

Experiments in dilution jet mixing

Experimental results are presented on the mixing of a single row of jets with an isothermal mainstream in a straight duct, with flow and geometric variations typical of combustion chambers in gas turbine engines included. It is found that at a constant momentum ratio, variations in the density ratio have only a second-order effect on the profiles. A first-order approximation to the mixing of jets with a variable temperature mainstream can, it is found, be obtained by superimposing the jets-in-an-isothermal-crossflow and mainstream profiles. Another finding is that the flow area convergence, especially injection-wall convergence, significantly improves the mixing. For opposed rows of jets with the orifice cone centerlines in-line, the optimum ratio of orifice spacing to duct height is determined to be 1/2 of the optimum value for single injection at the same momentum ratio. For opposed rows of jets with the orifice centerlines staggered, the optimum ratio of orifice spacing to duct height is found to be twice the optimum value for single side injection at the same momentum ratio.

Holdeman, J. D.↗
NASA NTRS2018

Investigating the Response of Loop Plasma to Nanoflare Heating Using RADYN Simulations

We present the results of 1D hydrodynamic simulations of coronal loops that are subject to nanoflares, caused by either in situ thermal heating or nonthermal electron (NTE) beams. The synthesized intensity and Doppler shifts can be directly compared with Interface Region Imaging Spectrograph (IRIS) and Atmospheric Imaging Assembly (AIA) observations of rapid variability in the transition region (TR) of coronal loops, associated with transient coronal heating. We find that NTEs with high enough low-energy cutoff (EC) deposit energy in the lower TR and chromosphere, causing blueshifts (up to approximately 20 kilometers per second) in the IRIS Si IV lines, which thermal conduction cannot reproduce. The EC threshold value for the blueshifts depends on the total energy of the events (approximately 5 kiloelectronvolts for 1024 ergs, up to 15 kiloelectronvolts for 1025 ergs). The observed footpoint emission intensity and flows, combined with the simulations, can provide constraints on both the energy of the heating event and EC. The response of the loop plasma to nanoflares depends crucially on the electron density: significant Si IV intensity enhancements and flows are observed only for initially low-density loops (less than 10 (sup 9) per cubic centimeter). This provides a possible explanation of the relative scarcity of observations of significant moss variability. While the TR response to single heating episodes can be clearly observed, the predicted coronal emission (AIA 94 Angstroms) for single strands is below current detectability and can only be observed when several strands are heated closely in time. Finally, we show that the analysis of the IRIS Mg II chromospheric lines can help further constrain the properties of the heating mechanisms.

Sun: Corona↗
DOE OSTI2024

An improved stochastic weighted particle method for boundary driven flows

Here, the stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Zinc dendrite removal in a nickel-zinc battery with flow-through electrodes

The development and deployment of inexpensive energy storage technologies is critical to realizing a clean energy grid. Batteries are being used in this role, but there remains a need for research on systems that are designed specifically for stationary energy storage, with a focus on lowering the overall cost rather than prioritizing the system energy density, specific energy, and power output. Here, we report the development of ultra-thick (1 cm thick) electrodes with engineered flow channels and explore the variables determining how thick these electrodes can feasibly be. Our proof of concept cell, utilizing the alkaline Ni-Zn chemistry, shows stable cycling over the initial 60 cycles but still suffers from the common Zn dendrite growth at the anode. To extend the life of these systems we report our novel methodology to completely remove Zn dendrites by exploiting the flow-through nature of our electrode architecture.

Collins-Wildman, Daniel L↗
NASA NTRS1991

Styles of ejecta emplacement under atmospheric conditions

Laboratory experiments provide essential first-order constraints on processes affecting ballistic ejecta and styles of ejecta emplacement under different atmospheric environments at planetary scales. The NASA-Ames Vertical Gun allows impacting different fine-grained particulate targets under varying atmospheric pressure and density, thereby helping to isolate controlling variables. Further analysis now permits characterizing distinct modes of emplacement that reflect the degree of ejecta entrainment within a turbidity flow created by ejecta curtain movement through the atmosphere.

Schultz, P. H.↗
NASA NTRS2012

Ensemble Averaged Probability Density Function (APDF) for Compressible Turbulent Reacting Flows

In this paper, we present a concept of the averaged probability density function (APDF) for studying compressible turbulent reacting flows. The APDF is defined as an ensemble average of the fine grained probability density function (FG-PDF) with a mass density weighting. It can be used to exactly deduce the mass density weighted, ensemble averaged turbulent mean variables. The transport equation for APDF can be derived in two ways. One is the traditional way that starts from the transport equation of FG-PDF, in which the compressible Navier- Stokes equations are embedded. The resulting transport equation of APDF is then in a traditional form that contains conditional means of all terms from the right hand side of the Navier-Stokes equations except for the chemical reaction term. These conditional means are new unknown quantities that need to be modeled. Another way of deriving the transport equation of APDF is to start directly from the ensemble averaged Navier-Stokes equations. The resulting transport equation of APDF derived from this approach appears in a closed form without any need for additional modeling. The methodology of ensemble averaging presented in this paper can be extended to other averaging procedures: for example, the Reynolds time averaging for statistically steady flow and the Reynolds spatial averaging for statistically homogeneous flow. It can also be extended to a time or spatial filtering procedure to construct the filtered density function (FDF) for the large eddy simulation (LES) of compressible turbulent reacting flows.

Shih, Tsan-Hsing↗
NASA NTRS2016

Uncertainty and Sensitivity Analysis of Afterbody Radiative Heating Predictions for Earth Entry

The objective of this work was to perform sensitivity analysis and uncertainty quantification for afterbody radiative heating predictions of Stardust capsule during Earth entry at peak afterbody radiation conditions. The radiation environment in the afterbody region poses significant challenges for accurate uncertainty quantification and sensitivity analysis due to the complexity of the flow physics, computational cost, and large number of un-certain variables. In this study, first a sparse collocation non-intrusive polynomial chaos approach along with global non-linear sensitivity analysis was used to identify the most significant uncertain variables and reduce the dimensions of the stochastic problem. Then, a total order stochastic expansion was constructed over only the important parameters for an efficient and accurate estimate of the uncertainty in radiation. Based on previous work, 388 uncertain parameters were considered in the radiation model, which came from the thermodynamics, flow field chemistry, and radiation modeling. The sensitivity analysis showed that only four of these variables contributed significantly to afterbody radiation uncertainty, accounting for almost 95% of the uncertainty. These included the electronic- impact excitation rate for N between level 2 and level 5 and rates of three chemical reactions in uencing N, N(+), O, and O(+) number densities in the flow field.

West, Thomas K., IV↗
NASA NTRS1983

Flow in coronal loops with a mass source

This research studies the flow of plasma inside a coronal loop in which an injection of plasma through the lateral surface is permitted. The flow is assumed steady and polytropic. The problem covers two cases: (a) upflow at one footpoint, downflow at the other; (b) downflow at both footpoints. The first case can be shown to be quite similar to that of a mass-conserving flow with variable cross section; the second, instead, is characterized by solutions with a different type of topology; its main new feature is the obvious fact that all the solutions pass through a single point going from negative to positive velocities. In this second case the density ratio between footpoints and top can be much smaller than in a mass conserving flow. This can explain some properties of observed loops.

Noci, G.↗
NASA NTRS1999

Plasma Flow During RF Discharges in VASIMR

The Variable Specific Impulse Magnetoplasma Rocket (VASIMR) plasma source consists of a helical antenna, driven at frequencies of 4 to 19 MHz with powers up to 1 kW, in a magnetic field up to 3 kG. Helium is the current test gas, and future experiments with hydrogen are planned. Plasma density and temperature profiles were measured by a reciprocating Langmuir probe, and plasma flow profiles were measured with a reciprocating Mach probe. Both probes were located about 0.5 m downstream from the helical antenna. The plasma source operated in capacitive and inductive modes in addition to a helicon mode. During capacitive and inductive modes, densities were low and plasma flow was < 0.5 Cs. When the plasma operated in a helicon mode, the densities measured downstream from the source were higher [10(exp 12) / cubic cm ] and plasma flow along the magnetic field was of the order Mach 1. Details of the measurements will be shown.

Jacobson, V. T.↗
NASA NTRS1991

A near-wall two-equation model for compressible turbulent flows

A near-wall two-equation turbulence model of the K - epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K - omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

Zhang, H. S.↗
NASA NTRS1991

Solar control of the earth's emission of energetic O(+)

Energetic (0.1-16 keV/e) O(+) data obtained in the earth's plasma sheet (between 10 and 23 RE) by an ion mass spectrometer on the ISEE-1 spacecraft are compared statistically with published data on the concurrent solar wind and IMF. The most strongly variable parameter of the plasma sheet O(+) is its density, which is found to be well correlated with certain solar wind parameters, especially with the solar wind flow speed and the IMF component perpendicular to the flow vector. When those two solar wind parameters are combined to form an electric field (-v x B), both the number density and the energy density of the O(+) are found to vary in proportion to the square of that electric field, on average, suggesting that the emission of energetic O(+) ions from the earth may be powered by that same field. Based on this and on the previously published correlation with solar activity, it is argued that the emission of O(+) is controlled by a combination of HF (ionizing) and quasi-static (accelerating) solar electromagnetic fields.

Lennartsson, W.↗
NASA NTRS1992

A near-wall two-equation model for compressible turbulent flows

A near-wall two-equation turbulence model of the k-epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K-omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

Zhang, H. S.↗
NASA NTRS1994

Bubble Formation at a Submerged Orifice in Reduced Gravity

The dynamic regime of gas injection through a circular plate orifice into an ideally wetting liquid is considered, when successively detached bubbles may be regarded as separate identities. In normal gravity and at relatively low gas flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole time of its evolution. If the flow rate exceeds a certain threshold value, another stage of the detachment process takes place in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical stem that lengthens as the bubble rises above the plate. The bubble shape resembles then that of a mushroom and the upper envelope continues to grow until the gas supply through the stem is completely cut off. Such a stage is always present under conditions of sufficiently low gravity, irrespective of the flow rate. Two major reasons make for bubble detachment: the buoyancy force and the force due to the momentum inflow into the bubble with the injected gas. The former force dominates the process at normal gravity whereas the second one plays a key role under negligible gravity conditions. It is precisely this fundamental factor that conditions the drastic influence on bubble growth and detachment that changes in gravity are able to cause. The frequency of bubble formation is proportional to and the volume of detached bubbles is independent of the gas flow rate in sufficiently low gravity, while at normal and moderately reduced gravity conditions the first variable slightly decreases and the second one almost linearly increases as the flow rate grows. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Buyevich, Yu A.↗
NASA NTRS1995

Effects of Gravity on Bubble Formation at a Plate Orifice

A model of the dynamic regime of gas injection through a submerged plate orifice into an ideally wetting liquid is developed in the circumstance when successively detached bubbles may be regarded as independent objects. Two major factors favor bubble detachment: the buoyancy force and a force due to the momentum inflow into the bubble with injected gas. In normal and moderately reduced gravity, the first factor dominates. At relatively low flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole period of its evolution till detachment. If the flow rate exceeds a critical value, another stage of bubble evolution occurs in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical gaseous stem that lengthens as the bubble rises above the plate. The bubble continues to grow until the gas supply through the stem is completely cut off, after which back flow of gas into the stem from the bubble becomes possible. In microgravity, the second factor prevails, and the latter stage is always present irrespective of the flow rate. However, the gas envelope rises and the stem lengthens very slowly. This difference in the underlying physical mechanisms provides for key properties of bubble growth and detachment being drastically different in appreciable and sufficiently reduced gravity. The frequency of bubble formation slightly decreases with and the detachment volume is almost proportional to the gas flow rate in the first case, in accordance with familiar relations. In the second case, the first variable is proportional to the flow rate whereas the second one is independent of it. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Webbon, Bruce W.↗