Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservation of mass”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Mass Conservation in a Chemical Transport Model and its Effect on CO2 and SF6 Simulations

Chemical transport models (CTMs) must conserve mass to be useful for applications involving assessment of the effect of various pollutants on the troposphere and stratosphere. Furthermore, calculations of the evolution of constituents such as SF6 are used to evaluate overall model transport, and interpretation of such simulations is clouded if mass conservation is not assured. For realistic simulations or predictions, it is crucial that constituents are not produced or lost by transport or other processes in the CTMs. Analysis of CO2 and SF6 experiments using a CTM shows that problems with mass conservation can seriously degrade the simulations. Failure to conserve mass results from inconsistency of the surface pressure tendency and the divergence of horizontal mass flux when the model is forced by assimilated meteorological data. We have developed an effective method to eliminate the inconsistency by modifying the divergent part of the wind field. The changes in the wind fields are quite small but the impact on mass conservation is large. Parameterizations of physical processes such as convection or turbulent transport can also affect mass conservation. The lack of conservation is small but accumulates when integrations are lengthy such as required for SF6. This lack of conservation is found using winds from either a GCM or from an assimilation system. A simple adjustment removes much of the inaccuracy in the convective parameterization. A CO2 simulation using assimilated winds from the most recent version of the Goddard Earth Observing System Data Assimilation System will be used to illustrate the impact of these transport improvements.

Zhu, Z.↗
NASA NTRS2013

Mass Conservation and Positivity Preservation with Ensemble-type Kalman Filter Algorithms

Maintaining conservative physical laws numerically has long been recognized as being important in the development of numerical weather prediction (NWP) models. In the broader context of data assimilation, concerted efforts to maintain conservation laws numerically and to understand the significance of doing so have begun only recently. In order to enforce physically based conservation laws of total mass and positivity in the ensemble Kalman filter, we incorporate constraints to ensure that the filter ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. We show that the analysis steps of ensemble transform Kalman filter (ETKF) algorithm and ensemble Kalman filter algorithm (EnKF) can conserve the mass integral, but do not preserve positivity. Further, if localization is applied or if negative values are simply set to zero, then the total mass is not conserved either. In order to ensure mass conservation, a projection matrix that corrects for localization effects is constructed. In order to maintain both mass conservation and positivity preservation through the analysis step, we construct a data assimilation algorithms based on quadratic programming and ensemble Kalman filtering. Mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate constraints. Some simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. The results show clear improvements in both analyses and forecasts, particularly in the presence of localized features. Behavior of the algorithm is also tested in presence of model error.

EnKR↗
DOE OSTI2025

A Mass Conservation Relaxed (MCR) LSTM Model for Streamflow Simulation Across CONUS

The recent development of the physics-aware Mass-Conserving Long Short-Term Memory network (MC-LSTM) provides an alternative to other data-driven Deep Learning (DL) models in hydrology. Mass-Conserving Long Short-Term Memory incorporates mass conservation directly into the LSTM architecture. Despite the theoretical advancements, studies have reported a surprisingly limited performance of the MC-LSTM in streamflow simulation. We hypothesize that such a limitation is due to the unrealistic mass conservation scheme in MC-LSTM, which overlooks unobserved incoming water fluxes beyond precipitation. As an attempt to verify this hypothesis, we propose a Mass Conservation Relaxed LSTM (MCR-LSTM), which incorporates a bi-directional mass relaxation (MR) component to account for potential incoming water fluxes beyond precipitation. We train and test the proposed MCR-LSTM model across 531 watersheds in the contiguous United States (CONUS) against three baseline models: the Sacramento Soil Moisture Accounting, LSTM, and MC-LSTM. Our results show that MCR-LSTM outperforms MC-LSTM despite its underperformance compared to LSTM. Specifically, MCR-LSTM's advantage over MC-LSTM is mainly seen in the Plains and Western U.S., where the newly incorporated MR component better simulates water loss and suggests the likely existence of additional incoming water fluxes beyond precipitation, respectively. The novelty and contribution of this study are twofold: firstly, it introduces an alternative physics-aware DL tool (i.e., MCR-LSTM) in hydrology with higher accuracy in specific regions compared to MC-LSTM. Secondly, it provides a diagnosis of regions where strict, precipitation-based mass conservation constraints may be unrealistic in streamflow simulation.

deep learning↗
NASA NTRS2011

Improving Simulated Soil Moisture Fields Through Assimilation of AMSR-E Soil Moisture Retrievals with an Ensemble Kalman Filter and a Mass Conservation Constraint

Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.

Li, Bailing↗
NASA NTRS1985

A 3-dimensional mass conserving element for compressible flows

A variety of finite element schemes has been used in the numerical approximation of compressible flows particularly in underwater acoustics. In many instances instabilities have been generated due to the lack of mass conservation. Two- and three-dimensional elements are developed which avoid these problems.

Fix, G.↗
NASA NTRS2009

Mass Conservation in Modeling Moisture Diffusion in Multi-Layer Carbon Composite Structures

Moisture diffusion in multi-layer carbon composite structures is difficult to model using finite difference methods due to the discontinuity in concentrations between adjacent layers of differing materials. Applying a mass conserving approach at these boundaries proved to be effective at accurately predicting moisture uptake for a sample exposed to a fixed temperature and relative humidity. Details of the model developed are presented and compared with actual moisture uptake data gathered over 130 days from a graphite epoxy composite sandwich coupon with a Rohacell foam core.

Nurge, Mark A.↗
NASA NTRS1994

The coupling between flame surface dynamics and species mass conservation in premixed turbulent combustion

Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Trouve, A.↗
NASA NTRS1982

Development of a locally mass flux conservative computer code for calculating 3-D viscous flow in turbomachines

The VANS successive approximation numerical method was extended to the computation of three dimensional, viscous, transonic flows in turbomachines. A cross-sectional computer code, which conserves mass flux at each point of the cross-sectional surface of computation was developed. In the VANS numerical method, the cross-sectional computation follows a blade-to-blade calculation. Numerical calculations were made for an axial annular turbine cascade and a transonic, centrifugal impeller with splitter vanes. The subsonic turbine cascade computation was generated in blade-to-blade surface to evaluate the accuracy of the blade-to-blade mode of marching. Calculated blade pressures at the hub, mid, and tip radii of the cascade agreed with corresponding measurements. The transonic impeller computation was conducted to test the newly developed locally mass flux conservative cross-sectional computer code. Both blade-to-blade and cross sectional modes of calculation were implemented for this problem. A triplet point shock structure was computed in the inducer region of the impeller. In addition, time-averaged shroud static pressures generally agreed with measured shroud pressures. It is concluded that the blade-to-blade computation produces a useful engineering flow field in regions of subsonic relative flow; and cross-sectional computation, with a locally mass flux conservative continuity equation, is required to compute the shock waves in regions of supersonic relative flow.

Walitt, L.↗
NASA NTRS1980

Gas bearings

The present work deals with the fundamentals of gas lubrication theory, which forms the foundation of all analytical design tools for gas bearings. Most of the hard lessons learned in the past are outlined with reference to dry contact, debris ingestion, sliding speed, and chemical stability of lubricant. The mathematical theory of gas lubrication is described for scaling rules in thin-film viscous flow, momentum conservation, mass conservation, energy conservation, isothermal gas bearing theory, coupling effects, and global bearing characteristics. Particular attention is given to the governing differential equations for common bearing configurations. Also discussed are representative solutions of self-acting gas bearings, externally pressurized bearings, and time-dependent effects.

Pan, C. H. T.↗
NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗
NASA NTRS2014

Conservation of Mass and Preservation of Positivity with Ensemble-Type Kalman Filter Algorithms

This paper considers the incorporation of constraints to enforce physically based conservation laws in the ensemble Kalman filter. In particular, constraints are used to ensure that the ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. In certain situations filtering algorithms such as the ensemble Kalman filter (EnKF) and ensemble transform Kalman filter (ETKF) yield updated ensembles that conserve mass but are negative, even though the actual states must be nonnegative. In such situations if negative values are set to zero, or a log transform is introduced, the total mass will not be conserved. In this study, mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate non-negativity constraints. Simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. In two examples, an update that includes a non-negativity constraint is able to properly describe the transport of a sharp feature (e.g., a triangle or cone). A number of implementation questions still need to be addressed, particularly the need to develop a computationally efficient quadratic programming update for large ensemble.

EnKF↗
NASA NTRS1978

Prediction of plasma properties in mercury ion thrusters

A simplified theoretical model was developed which obtains to first order the plasma properties in the discharge chamber of a mercury ion thruster from basic thruster design and controllable operating parameters. The basic operation and design of ion thrusters is discussed, and the important processes which influence the plasma properties are described in terms of the design and control parameters. The conservation for mass, charge and energy were applied to the ion production region, which was defined as the region of the discharge chamber having as its outer boundary the surface of revolution of the innermost field line to intersect the anode. Mass conservation and the equations describing the various processes involved with mass addition and removal from the ion production region are satisfied by a Maxwellian electron density spatial distribution in that region.

Longhurst, G. R.↗
DOE OSTI2024

Structure and dynamics of magneto-inertial, differentially rotating laboratory plasmas

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}^{-1}$ . The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$ . As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}^{-1}$ . We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.

Valenzuela-Villaseca, V. (ORCID:0000000167721441)↗
DOE OSTI2024

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Calculation of flow about two-dimensional bodies by means of the velocity-vorticity formulation on a staggered grid

A method for calculating the incompressible viscous flow about two-dimensional bodies, utilizing the velocity-vorticity form of the Navier-Stokes equations using a staggered-grid formulation is presented. The solution is obtained by employing an alternative-direction implicit method for the solution of the block tridiagonal matrix resulting from the finite-difference representation of the governing equations. The boundary vorticity and the conservation of mass are calculated implicitly as a part of the solution. The mass conservation is calculated to machine zero for the duration of the computation. Calculations for the flow about a circular cylinder, a 2-pct thick flat plate at 90-deg incidence, an elliptic cylinder at 45-deg incidence, and a NACA 0012, with and without a deflected flap, at - 90-deg incidence are performed and compared with the results of other numerical investigations.

Stremel, Paul M.↗
DOE OSTI2025

Towards smarter green infrastructure: Fusing bark ecology and stemflow hydrodynamics on tree stems

A wide array of bark surfaces sheath wooded plants in rural and urban areas alike. Much work has examined the function and role of bark in different contexts and different environs, including urban areas, finding that it is rich in life and can play a role in the transfer of water and matter to the ground surface. Accordingly, this paper presents a first step to weld and fuse bark ecology and stemflow hydrodynamics. It is an effort to develop a physically-based understanding of the transport of water and matter (e.g., solutes, particulates, microorganisms) along tree stems using relevant equations to allow a more informed consideration of bark in green infrastructure initiatives. In particular, the hydrodynamical equations are based on the conservation of water mass, conservation of momentum, and conservation of scalar mass. These equations, coupled with contemplation of corticular life, underpin and substantiate bark’s unifying role as a modulator and cultivator. By elucidating the ‘black box’ of the tree stem and utilizing the formulations set forth in this paper, urban foresters and planners can develop green infrastructure to help advance ecosystem services and sustainability development goals (SDG), especially SDG 11 and SDG 15.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗