Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity equation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Variation objective analyses for cyclone studies

The objectives were to: (1) develop an objective analysis technique that will maximize the information content of data available from diverse sources, with particular emphasis on the incorporation of observations from satellites with those from more traditional immersion techniques; and (2) to develop a diagnosis of the state of the synoptic scale atmosphere on a much finer scale over a much broader region than is presently possible to permit studies of the interactions and energy transfers between global, synoptic and regional scale atmospheric processes. The variational objective analysis model consists of the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation for a dry hydrostatic atmosphere. Preliminary tests of the model with the SESMAE I data set are underway for 12 GMT 10 April 1979. At this stage of purpose of the analysis is not the diagnosis of atmospheric structures but rather the validation of the model. Model runs for rawinsonde data and with the precision modulus weights set to force most of the adjustment of the wind field to the mass field have produced 90 to 95 percent reductions in the imbalance of the initial data after only 4-cycles through the Euler-Lagrange equations. Sensitivity tests for linear stability of the 11 Euler-Lagrange equations that make up the VASP Model 1 indicate that there will be a lower limit to the scales of motion that can be resolved by this method. Linear stability criteria are violated where there is large horizontal wind shear near the upper tropospheric jet.

Achtemeier, G. L.↗
NASA NTRS1988

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.

Achtemeier, Gary L.↗
NASA NTRS1988

A multivariate variational objective analysis-assimilation method. Part 2: Case study results with and without satellite data

The variational multivariate assimilation method described in a companion paper by Achtemeier and Ochs is applied to conventional and conventional plus satellite data. Ground-based and space-based meteorological data are weighted according to the respective measurement errors and blended into a data set that is a solution of numerical forms of the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation for a dry atmosphere. The analyses serve first, to evaluate the accuracy of the model, and second to contrast the analyses with and without satellite data. Evaluation criteria measure the extent to which: (1) the assimilated fields satisfy the dynamical constraints, (2) the assimilated fields depart from the observations, and (3) the assimilated fields are judged to be realistic through pattern analysis. The last criterion requires that the signs, magnitudes, and patterns of the hypersensitive vertical velocity and local tendencies of the horizontal velocity components be physically consistent with respect to the larger scale weather systems.

Achtemeier, Gary L.↗
NASA NTRS1991

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. Previous analyses of this concept were based on inviscid assumptions. The present results indicate that viscous effects are of primary importance; in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure, and turbulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. In the unstable case, a detonation wave is formed, which propagates upstream and unstarts the ram accelerator. In the stable case, a solution that converges to steady-state is obtained, in which the combustion wave remains stationary with respect to the ram accelerator projectile. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. In addition to these studies, two numerical techniques were tested. These two techniques are vector extrapolation to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices that arise in chemically reacting flows.

Yungster, Shaye↗
NASA NTRS1991

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. The present results indicate that viscous effects are of primary importance in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure and tubulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. Two numerical techniques were tested: vector extrapolation, to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices which arise in chemically reacting flows.

Yungster, Shaye↗
NASA NTRS1995

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗
NASA NTRS2008

Numerical Simulations of a Co-Axial Supersonic-Combusting Free-Jet Experiment

CFD calculations using the Reynolds-averaged Navier-Stokes equations coupled with species continuity equations have been made for a supersonic coaxial-jet CFD-validation experiment to determine the sensitivity of the external flowfield to the main-nozzle exit profile. Four different nozzle exit profiles were used in the study: a uniform profile, one computed using only the nozzle geometry, one computed using the nozzle geometry and part of the upstream facility combustor, and one using the nozzle and the full facility combustor. Two cases were examined using the four profiles: a non-reacting case without coflow and a reacting case with hydrogen coflow. Results show that the nozzle exit profile has a significant effect on the external flowfield. The uniform profile produced the longest jet while the profile created with the full combustor produced the shortest jet. The nozzle-only and part-combustor profiles fell between the other two profiles. The reacting flow was found to be more sensitive to the nozzle exit profile since it affects the downstream mixing and combustion. These calculations indicate the importance of properly setting the nozzle-exit profile for this type of calculation.

Gaffney, Richard L.↗
DOE OSTI2024

A simplified all-frequency stable formulation with an implicit Coulomb gauge

A potential-based finite element formulation has been developed in the past to circumvent the low-frequency breakdown issues commonly encountered in electromagnetic (EM) simulations of low-frequency and multiscale problems. In this formulation, the magnetic vector and electric scalar potentials have been employed to express the electric field and magnetic flux, leading to a set of two equations that represents all four Maxwell’s equations and the current continuity equation. To enforce the Coulomb gauge, an auxiliary potential has been introduced, which results in a total of three equations to be solved simultaneously. To reduce the number of equations and unknowns needed in a simulation, a simplified formulation is proposed in this paper to enforce the Coulomb gauge implicitly without the need for the auxiliary potential. A numerical example is given to demonstrate the accuracy of the proposed formulation.

Mekonnen, Minyechil↗
DOE OSTI2022

Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++

Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Computer program for aerodynamic and blading design of multistage axial-flow compressors

A code for computing the aerodynamic design of a multistage axial-flow compressor and, if desired, the associated blading geometry input for internal flow analysis codes is presented. Compressible flow, which is assumed to be steady and axisymmetric, is the basis for a two-dimensional solution in the meridional plane with viscous effects modeled by pressure loss coefficients and boundary layer blockage. The radial equation of motion and the continuity equation are solved with the streamline curvature method on calculation stations outside the blade rows. The annulus profile, mass flow, pressure ratio, and rotative speed are input. A number of other input parameters specify and control the blade row aerodynamics and geometry. In particular, blade element centerlines and thicknesses can be specified with fourth degree polynomials for two segments. The output includes a detailed aerodynamic solution and, if desired, blading coordinates that can be used for internal flow analysis codes.

Crouse, J. E.↗
NASA NTRS1981

Basic lubrication equations

Lubricants, usually Newtonian fluids, are assumed to experience laminar flow. The basic equations used to describe the flow are the Navier-Stokes equation of motion. The study of hydrodynamic lubrication is, from a mathematical standpoint, the application of a reduced form of these Navier-Stokes equations in association with the continuity equation. The Reynolds equation can also be derived from first principles, provided of course that the same basic assumptions are adopted in each case. Both methods are used in deriving the Reynolds equation, and the assumptions inherent in reducing the Navier-Stokes equations are specified. Because the Reynolds equation contains viscosity and density terms and these properties depend on temperature and pressure, it is often necessary to couple the Reynolds with energy equation. The lubricant properties and the energy equation are presented. Film thickness, a parameter of the Reynolds equation, is a function of the elastic behavior of the bearing surface. The governing elasticity equation is therefore presented.

Hamrock, B. J.↗
NASA NTRS1992

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Averaged Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗
NASA NTRS1993

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Average Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗
NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗
NASA NTRS1999

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub 50 nm MOSFETs: a Statistical 3D 'Atomistic' Simulation Study

A 3D 'atomistic' simulation study of random dopant induced threshold voltage fluctuations and lowering in sub 50 nm MOSFETs is presented. The attention is focused mainly on devices with 30 nm effective channel length which represent the expected level of scaling at the end of the Silicon Roadmap. An efficient algorithm, based on a single 3D ap solution of the Poisson equation and a simplified current continuity equation, is used in the simulations. Large samples of microscopically different devices (typically 200) arc used in order to obtain statistically reliable results. The influence of different aspects of the conventional MOSFET design on the threshold voltage fluctuations and lowering are investigated. Results for fluctuation resistant device architectures based on low-doped epitaxial channel MOSFETs are also presented.

Asenov, Asen↗
NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗