Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass conservation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Generalized Fluid System Simulation Program

A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

Alok Kumar Majumdar↗
NASA NTRS2004

Numerical Modeling of Conjugate Heat Transfer in Fluid Network

Fluid network modeling with conjugate heat transfer has many applications in Aerospace engineering. In modeling unsteady flow with heat transfer, it is important to know the variation of wall temperature in time and space to calculate heat transfer between solid to fluid. Since wall temperature is a function of flow, a coupled analysis of temperature of solid and fluid is necessary. In cryogenic applications, modeling of conjugate heat transfer is of great importance to correctly predict boil-off rate in propellant tanks and chill down of transfer lines. In TFAWS 2003, the present author delivered a paper to describe a general-purpose computer program, GFSSP (Generalized Fluid System Simulation Program). GFSSP calculates flow distribution in complex flow circuit for compressible/incompressible, with or without heat transfer or phase change in all real fluids or mixtures. The flow circuit constitutes of fluid nodes and branches. The mass, energy and specie conservation equations are solved at the nodes where as momentum conservation equations are solved at the branches. The proposed paper describes the extension of GFSSP to model conjugate heat transfer. The network also includes solid nodes and conductors in addition to fluid nodes and branches. The energy conservation equations for solid nodes solves to determine the temperatures of the solid nodes simultaneously with all conservation equations governing fluid flow. The numerical scheme accounts for conduction, convection and radiation heat transfer. The paper will also describe the applications of the code to predict chill down of cryogenic transfer line and boil-off rate of cryogenic propellant storage tank.

Majumdar, Alok↗
NASA NTRS2005

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗
NASA NTRS2009

Analysis of Stationary, Photovoltaic-based Surface Power System Designs at the Lunar South Pole

Combinations of solar arrays and either batteries or regenerative fuel cells are analyzed for a surface power system module at the lunar south pole. The systems are required to produce 5 kW of net electrical power in sunlight and 2 kW of net electrical power during lunar night periods for a 10-year period between 2020 and 2030. Systems-level models for energy conservation, performance, degradation, and mass are used to compare to various systems. The sensitivities of important and/or uncertain variables including battery specific energy, fuel cell operating voltage, and DC-DC converter efficiency are compared to better understand the system. Switching unit efficiency, battery specific energy, and fuel cell operating voltage appear to be important system-level variables for this system. With reasonably sized solar arrays, the regenerative fuel cell system has significantly lower mass than the battery system based on the requirements and assumptions made herein. The total operational time is estimated at about 10,000 hours in battery discharge/fuel cell mode and about 4,000 and 8,000 hours for the battery charge and electrolyzer modes, respectively. The estimated number of significant depth-of-discharge cycles for either energy storage system is less than 100 for the 10-year period.

Freeh, Joshua E.↗
NASA NTRS2009

Method and Apparatus for Predicting Unsteady Pressure and Flow Rate Distribution in a Fluid Network

A method and apparatus for analyzing steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics, external body forces such as gravity and centrifugal force and conjugate heat transfer. In some embodiments, a graphical user interface provides for the interactive development of a fluid network simulation having nodes and branches. In some embodiments, mass, energy, and specific conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. In some embodiments, contained herein are data objects for computing thermodynamic and thermophysical properties for fluids. In some embodiments, the systems of equations describing the fluid network are solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods.

Majumdar, Alok K.↗
NASA NTRS2009

Evidence for an Intermediate Mass Black Hole in NGC 5408 X-1

We report the discovery with XMM-Newton of correlated spectral and timing behavior in the ultraluminous X-ray source (ULX) NGC 5408 X-1. An approx. 100 ksec pointing with XMM/Newton obtained in January, 2008 reveals a strong 10 mHz QPO in the > 1 keV flux, as well as flat-topped, band limited noise breaking to a power law. The energy spectrum is again dominated by two components, a 0.16 keV thermal disk and a power-law with an index of approx. 2.5. These new measurements, combined with results from our previous January 2006 pointing in which we first detected QPOs, show for the first time in a ULX a pattern of spectral and temporal correlations strongly analogous to that seen in Galactic black hole sources, but at much higher X-ray luminosity and longer characteristic time-scales. We find that the QPO frequency is proportional to the inferred disk flux, while the QPO and broad-band noise amplitude (root mean squared, rms) are inversely proportional to the disk flux. Assuming that QPO frequency scales inversely with black hole mass at a given power-law spectral index we derive mass estimates using the observed QPO frequency - spectral index relations from five stellar-mass black hole systems with dynamical mass constraints. The results from all sources are consistent with a mass range for NGC 5408 X-1 from 1000 - 9000 Stellar mass. We argue that these are conservative limits, and a more likely range is from 2000 - 5000 Stellar mass. Moreover, the recent relation from Gierlinski et al. that relates black hole mass to the strength of variability at high frequencies (above the break in the power spectrum), and the variability plane results of McHardy et al. and Koerding et al., are also suggestive of such a. high mass for NGC 5408 X-1. Importantly, none of the above estimates appears consistent with a black hole mass less than approx. 1000 Stellar mass for NGC 5408 X-1. We argue that these new findings strongly support the conclusion that NGC 5408 X-1 harbors an intermediate mass black hole.

Strohmayer, Tod E.↗
NASA NTRS2011

Magnetic Reconnection in a Compressible MHD Plasma

Using steady-state resistive MHD, magnetic reconnection is reinvestigated for conditions of high resistivity/low magnetic Reynolds number, when the thickness of the diffusion region is no longer small compared to its length. Implicit expressions for the reconnection rate and other reconnection parameters are derived based on the requirements of mass, momentum, and energy conservation. These expressions are solved via simple iterative procedures. Implications specifically for low Reynolds number/high resistivity are being discussed

Hesse, Michael↗
NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗
NASA NTRS2017

Advection of Potential Temperature in the Atmosphere of Irradiated Exoplanets: A Robust Mechanism to Explain Radius Inflation

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a two-dimensional steady-state atmospheric circulation model, the validity of which is assessed by comparison to three-dimensional calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occurring in the Earth's atmosphere or oceans. In the deep atmosphere, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD 209458b, but also reproduce the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non-uniform atmospheric heating thus provides a robust mechanism to explain the inflated radii of irradiated hot Jupiters.

tmospheric effects; methods: numerical; planets an↗
NASA NTRS2018

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗
NASA NTRS2002

Lidar Measurements of Wind, Moisture, and Boundary Layer Evolution in a Dry Line during 1HOP 2002

Variability in the convective boundary layer moisture, wind and temperature fields and their importance in the forecasting and understanding of storms have been discussed in the literature. These . variations have been reported in relation to frontal zones, stationary boundaries and during horizontal convective rolls. While all three vary substantially in the convective boundary layer, moisture poses a particular challenge. Moisture or water vapor concentration (expressed as a mass mixing ratio, g/kg), is conserved in all meteorological processes except condensation and evaporation. The water vapor mixing ratio often remains distinct across an air-mass boundary even when the temperature difference is indistinct. These properties make it an ideal choice in visualizing and understanding many of the atmosphere's dynamic features. However, it also presents a unique measurement challenge because water vapor content can vary by more than three orders of magnitude in the troposphere. Characterization of the 3D-distribution of water vapor is also difficult as water vapor observations can suffer from large sampling errors and substantial variability both in the vertical and horizontal. This study presents ground-based measurements of wind, boundary layer structure and water vapor mixing ratio measurements observed by three co-located lidars. This presentation will focus on the evolution and variability of moisture and wind in the boundary layer during a dry line event that occurred on 22 May 2002. These data sets and analyses are unique in that they combine simultaneous measurements of wind, moisture and CBL structure to study the detailed thermal variability in and around clear air updrafts during a dryline event. It will quantify the variation caused by, in and around buoyant plumes and across a dryline. The data presented here were collected in the panhandle of Oklahoma as part of the International H2O Project (MOP-2002), a field experiment that took place over the Southern Great Plains (SGP) of the United States from 13 May to 30 June 2002. The chief goal of MOP-2002 is to improve characterization of the four-dimensional (4-D) distribution of water vapor and its application to improving the understanding and prediction of convection

Demoz, Belay↗
DOE OSTI2025

Stemflow Hydrodynamics

Stemflow hydrodynamics is the study of water movement along the exterior surface area of plants. Its primary goal is to describe water velocity and water depth along the stem surface area. Its significance in enriching the rhizosphere with water and nutrients is not in dispute. Yet, the hydrodynamics of stemflow have been entirely overlooked. This review seeks to fill this knowledge gap by drawing from thin film theories to seek outcomes at the tree scale. The depth‐averaged conservation equations of water and solute mass are derived at a point. These equations are then supplemented with the conservation of momentum that is required to describe water velocities or relations between water velocities and water depth. Relevant forces pertinent to momentum conservation are covered and include body forces (gravitational effects), surface forces (wall friction), line forces (surface tension), and inertial effects. The inclusion of surface tension opens new vistas into the richness and complexity of stemflow hydrodynamics. Flow instabilities such as fingering, pinching of water columns into droplets, accumulation of water within fissures due to surface tension and their sudden release are prime examples that link observed spatial patterns of stemflow fronts and morphological characteristics of the bark. Aggregating these effects at the tree‐ and storm‐ scales are featured using published experiments. The review discusses outstanding challenges pertaining to stemflow hydrodynamics, the use of dynamic similarity and 3D printing to enable the interplay between field studies and controlled laboratory experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

Lin, Shian-Jiann↗
NASA NTRS2014

Cabin Environment Physics Risk Model

This paper presents a Cabin Environment Physics Risk (CEPR) model that predicts the time for an initial failure of Environmental Control and Life Support System (ECLSS) functionality to propagate into a hazardous environment and trigger a loss-of-crew (LOC) event. This physics-of failure model allows a probabilistic risk assessment of a crewed spacecraft to account for the cabin environment, which can serve as a buffer to protect the crew during an abort from orbit and ultimately enable a safe return. The results of the CEPR model replace the assumption that failure of the crew critical ECLSS functionality causes LOC instantly, and provide a more accurate representation of the spacecraft's risk posture. The instant-LOC assumption is shown to be excessively conservative and, moreover, can impact the relative risk drivers identified for the spacecraft. This, in turn, could lead the design team to allocate mass for equipment to reduce overly conservative risk estimates in a suboptimal configuration, which inherently increases the overall risk to the crew. For example, available mass could be poorly used to add redundant ECLSS components that have a negligible benefit but appear to make the vehicle safer due to poor assumptions about the propagation time of ECLSS failures.

Risk↗
DOE OSTI2025

CUPID, the Cuore upgrade with particle identification

CUPID, the CUORE Upgrade with Particle Identification, is a next-generation experiment to search for neutrinoless double beta decay (0 νββ ) and other rare events using enriched Li 2 100 MoO 4 scintillating bolometers. It will be hosted by the CUORE cryostat located at the Laboratori Nazionali del Gran Sasso in Italy. The main physics goal of CUPID is to search for 0 νββ of 100 Mo with a discovery sensitivity covering the full neutrino mass regime in the inverted ordering scenario, as well as the portion of the normal ordering regime with lightest neutrino mass larger than 10 meV. With a conservative background index of 10 −4 cts/(keV·kg·yr), 240 kg isotope mass, 5 keV FWHM energy resolution at 3 MeV and 10 live-years of data taking, CUPID will have a 90% C.L. half-life exclusion sensitivity of 1.8 · 10 27 yr, corresponding to an effective Majorana neutrino mass ( m ββ ) sensitivity of 9–15 meV, and a 3σ discovery sensitivity of 1 · 10 27 yr, corresponding to an m ββ range of 12–21 meV.

Alfonso, K. [Virginia Polytechnic Institute and St↗
DOE OSTI2025

Modeling the differential rate for signal interactions in coincidence with noise fluctuations or large rate backgrounds

The characteristic energy of a relic dark matter interaction with a detector scales strongly with the putative dark matter mass. Consequently, experimental search sensitivity at the lightest masses will always come from interactions whose size is similar to noise fluctuations and low energy backgrounds in the detector. In this paper, we correctly calculate the net change in measured differential event rate due to low rate signal interactions that overlap in time with noise and backgrounds, accounting for both periods of time when the signal is coincident with noise/backgrounds and for the decreased amount of time in which only noise/backgrounds occur and we show that the introduction of random simulated signal events into the continuous raw data stream (a form of “salting”) provides a correct and practical implementation to estimate this net linear signal sensitivity. Unfortunately this does not apply to the situation with non negligible pileups as we show through explicit examples. Consequently, this exclusion technique should be complemented by less sensitive techniques that can separately exclude large rate, high pileup signal parameter space. An analysis threshold above the Gaussian-like noise component of the measured differential rate spectrum should also likely produce conservative limits. Though previous light mass dark matter searches did not correctly account for decreased background only live time effects in their signal sensitivity, conservative analysis choices likely protected their published limits from being nonconservative.

Li, Xinran↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗