Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass balance”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

The H2o Content of the Ureilite Parent Body

The fate of highly volatile elements (H, C, F, Cl and S) during planetary accretion and differentiation is debated. Recent analyses of water in non-carbonaceous chondrites (RC, OC, EC) and achondrites (angrites, eucrites) have been used to argue that inner solar system parent bodies accreted and retained their highly volatile element budgets from their primary feedstock without substantial loss during accretion, metamorphism and differentiation. An alternative model posits that differentiated inner solar system parent bodies (e.g., the angrite parent body, 4 Vesta, Earth) derived the majority of their water from a carbonaceous chondrite-like source, delivered during the final stages of accretion. In order to add new constraints to this debate, we have measured water in nominally anhydrous minerals, melt inclusions, and interstitial glass in ureilites, the largest group of primitive achondrites in the terrestrial meteorite collection. Primitive achondrites did not experience global melting and homogenization. Therefore, these meteorites capture part of the transition from chondritic to achondritic parent bodies, allowing us to constrain the fate of water during the earliest stages of differentiation. Our nano-scale secondary ion mass spectrometry (nanoSIMS) analyses allow us to assess the viability of ureilite-like material as a potential source of terrestrial water. Analyses of pigeonite in main group ureilites yield a range of 2.0 – 6.0 µg/g H2O, and analyses of high-Ca pyroxene and glass (glassy melt inclusions and interstitial glass) in the Almahata Sitta ureilitic trachyandesite yield ranges of 13 – 19 µg/g H2O and 44 – 216 µg/g H2O, respectively. Mass balance, incremental melting, and batch melting calculations yield a preferred ureilite parent body H2O content of 2 – 20 µg/g, similar to previous estimates of water in the eucrite parent body (4 Vesta), but lower than estimates of Earth’s water budget. With these data, we demonstrate that 1) the ureilite parent body is H2O-depleted relative to the Earth; 2) ureilite-like material is unlikely to be a primary source of H2O to the Earth; 3) C and H are not necessarily coupled elements during planetary accretion and thermal processing; and 4) accretion, heating, partial melting, and degassing of rocky planetesimals likely results in significant depletion of H2O.

Primitive achondrites↗
NASA NTRS2024

Global Optimization via Quadratic Disjunctive Programming for Water Networks Design with Energy Recovery

Generalized disjunctive programming (GDP) models with bilinear and concave constraints, often seen in water network design, are challenging optimization problems. This work proposes quadratic and piecewise linear approximations for nonlinear terms to reformulate GDP models into quadratic GDP (QGDP) models that suitable solvers may solve more efficiently. We illustrate the benefits of the quadratic reformulation with a water treatment network design problem in which nonconvexities arise from bilinear terms in the mixers’ mass balances and concave investment cost functions of treatment units. Given the similarities with water network design problems, we suggest quadratic approximation for the GDP model for the optimal design of a large-scale reverse electrodialysis (RED) process. This power technology can recover energy from salinity differences between by-product streams of the water sector, such as desalination brine mixed with regenerated wastewater effluents. The solver Gurobi excels in handling QGDP problems, but weighing the problem’s precision and tractability balance is crucial. The piecewise linear approximation yields more accurate, yet larger QGDP models that may require longer optimization times in large-scale process synthesis problems.

Water Networks↗
NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗
DOE OSTI2021

Advanced Zonal Infiltration Measurement Method for Multifamily Buildings: A Novel Test Procedure to Determine Air Leakage Through External and Internal Surfaces

Measuring air infiltration through the external building envelope as well as between adjacent zones within a building is important since air leakage is related to increased energy use and poor indoor air quality (IAQ). State of the art methods can determine both leakages, however it can be both time consuming and labor intensive as it requires running both guarded and compartmentalization tests. This paper presents our work on developing a simplified testing procedure to determine power law flow equation leakage parameters of both external and internal boundaries of tested zones. The proposed Zonal Multipoint Pressure Test (ZMPT) method modulates airflow in the zone using a variable airflow measurement panel (mounted in either a window or the door), and solving the mass balance equations for the leakage parameters. ZMPT results were compared to conventional test results and found to correlate to within +/-7%.

Jayarathne, Tharanga↗
DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗
DOE OSTI2020

Gamma Spectrometry Examination of the AGR-3/4 Irradiation

The results from gamma spectrometry examination of the different components from the combined third and fourth U.S. Advanced Gas Reactor (AGR) TRISO-coated particle fuel irradiation tests (AGR-3/4) have been analyzed. This experiment was designed to provide information about in-pile fission product migration. In each of the 12 capsules, a single stack of four compacts with designed-to-fail (DTF) particles surrounded by inner and outer graphite and/or graphitic matrix rings and a graphite sink ring were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Gamma spectrometry has been used to evaluate the gamma-emitting fission product inventory of compacts from the irradiation and evaluate the burnup of these compacts based on the activity of the radioactive cesium isotopes (i.e., Cs-134 and Cs-137) in the compacts. Burnup from gamma spectrometry compares well with predicted burnup from simulations. The inner and outer rings were also examined by gamma spectrometry to evaluate the total fission product inventory of the rings, and gamma emission computed tomography (GECT) was used to investigate the spatial distribution of gamma-emitting fission products within the rings. This report focuses on results obtained via the Precision Gamma Scanner (PGS) for the compacts and inner and outer rings. These non-destructive gamma measurements are currently being compared to recent destructive measurements of spatial fission product distributions in the rings and to predictions made via mathematical modeling. A separate report discusses the entire AGR-3/4 fission product mass balance for all irradiation capsule components including the inner and outer rings, graphite sink rings, capsule spacers, capsule foils, through tubes, and felts. This report is similar to the conference proceedings: Harp, J.M., Demkowicz, P.A., and Stempien, J.D., “Initial gamma spectrometry examination of the AGR-3/4 irradiation,” International Topical Meeting High Temperature Reactor Technology (HTR 2016), Las Vegas, NV, USA, November 2016, Paper HTR2016-18593, but it has been updated to reflect all the data that has been collected on AGR-3/4 using PGS. The Appendices contain extensive data from the collected gamma spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Impact of Tree Age and Anatomical Fraction on Forest Residue Pyrolysis and Hydrotreating

Utilization of cost-advantaged biomass and waste resources to produce clean, domestic biofuels will be a key factor in decarbonizing transportation in the United States. This study investigates how the distribution of anatomical fractions (stem wood, bark, needles, branches) from pine forest residues of different ages impacts pyrolysis and hydrotreating processes, including product composition, carbon efficiency, economics, and overall sustainability. Commercially sourced 13- and 23-year-old loblolly pine residues were chosen to represent available resources; whole tree thinnings, and tops and branches from merchantable timber, respectively. Detailed characterization of the residue samples, product distributions, bio-oil, char, and gas characterization data are presented for 13 experiments, including the individual fractions, blends, and air-classified samples to reduce ash content. Pyrolysis mass balances averaged 97 +/- 3%. As expected, bio-oil yields were lower for the bark, needles, and blends with high proportions of these (53-56% dry basis), and highest for the debarked stem wood (67-71% dry basis). Lower yields were attributed to the high ash content of the needles, but also to the high lignin content in the bark (>50%), which is hypothesized to have led to increased char formation during pyrolysis. Higher-oil yield was generally correlated with total volatiles, hydrogen, oxygen, glucan, and mannan; somewhat correlated with xylan and acetyl content; and negatively correlated with total carbon, fixed carbon, nitrogen, sulfur, and extractives. Air classification of the residues was effective in reducing ash content (1.5% to 0.64%), extractives (~7% to 2%), and sulfur (0.04% to 0.02%), resulting in a small increase in yield, a significant increase in the GC-detectable proportion of oil (27% to 34%) and a net economic benefit. These experimental results were used to conduct technoeconomic and life cycle analyses, which are also summarized here.

biofuels↗
NASA NTRS1990

The influence of Oort clouds on the mass and chemical balance of the interstellar medium

The contribution of stellar encounters and interstellar erosion to comet cloud mass injection to the ISM is calculated. It is shown that evaporative mass loss from passing stars and SNe results in an average Galactic mass injection rate of up to 10 to the -5th solar mass/yr if such clouds are frequent around solar-type stars. Cometary erosion by interstellar grains produces an injection rate of 10 to the -5th to 10 to the -4th solar mass/yr. An injection rate of 2 x 10 to the -5th solar mass/yr is calculated. Each of these rates could be increased by a factor of about 15 if the comet clouds contain a significant amount of smaller debris. It is concluded that the total mass injection rate of material to the ISM by comet clouds is small compared to other ISM mass injection sources. Comet cloud mass loss to the ISM could be responsible for a sizeable fraction of the metal and dust abundances of the ISM if Oort clouds are common.

Stern, S. Alan↗
NASA NTRS2015

Maintaining Atmospheric Mass and Water Balance Within Reanalysis

This report describes the modifications implemented into the Goddard Earth Observing System Version-5 (GEOS-5) Atmospheric Data Assimilation System (ADAS) to maintain global conservation of dry atmospheric mass as well as to preserve the model balance of globally integrated precipitation and surface evaporation during reanalysis. Section 1 begins with a review of these global quantities from four current reanalysis efforts. Section 2 introduces the modifications necessary to preserve these constraints within the atmospheric general circulation model (AGCM), the Gridpoint Statistical Interpolation (GSI) analysis procedure, and the Incremental Analysis Update (IAU) algorithm. Section 3 presents experiments quantifying the impact of the new procedure. Section 4 shows preliminary results from its use within the GMAO MERRA-2 Reanalysis project. Section 5 concludes with a summary.

IAU↗
DOE OSTI2022

A mass conservative, well balanced, tangency preserving and energy decaying method for the shallow water equations on a sphere

Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Spacebased Observations of Oceanic Influence on the Annual Variation of South American Water Balance

The mass change of South America (SA) continent measured by the Gravity Recovery and Climate Experiment (GRACE) imposes a constraint on the uncertainties in estimating the annual variation of rainfall measured by Tropical Rain Measuring Mission (TRMM) and ocean moisture influx derived from QuikSCAT data. The approximate balance of the mass change rate with the moisture influx less climatological river discharge, in agreement with the conservation principle, bolsters not only the credibility of the spacebased measurements, but supports the characterization of ocean's influence on the annual variation of continental water balance. The annual variation of rainfall is found to be in phase with the mass change rate in the Amazon and the La Plata basins, and the moisture advection across relevant segments of the Pacific and Atlantic coasts agrees with the annual cycle of rainfall in the two basins and the Andes mountains.

rainfall annual cycle↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2022

Effects of thermophoresis on high-pressure binary-species boundary layers with uniform and non-uniform compositions

Direct numerical simulations of high-pressure binary-species temporal boundary layers are performed to investigate the flow physics for three situations: (1) uniform and equal composition, (2) uniform but unequal compositions and (3) non-uniform composition. Both colder- and hotter-wall situations compared with the free stream are simulated. The working fluid is a nitrogen/methane mixture. The analysis is performed at a case-specific self-similar state. Even when the initial composition is uniform, the methane mean mass fraction decreases near the colder wall, whereas it increases near the hotter wall and the mass fraction fluctuates in the entire boundary layer. Analysis of the species-mass diffusion balance and flow structures reveal that both mass-fraction variation and fluctuations are induced by the Soret effect. When the initial composition is non-uniform and the wall is colder, the methane mean mass fraction monotonically increases from the wall akin to its initial profile. However, when the wall is hotter the mean mass fraction decreases near the wall in contrast to its initial profile, a fact traced through the species-mass diffusion balance to the Soret effect being large and enriching methane near the wall. In contrast, the direction of the Soret flux is opposite for the colder wall, thus keeping the methane concentration small. Although the initial magnitude of the difference between the wall and free-stream temperature is the same in all cases, the situation is not symmetric between colder- and hotter-wall cases; the flow structure exhibits much smaller scales when the wall is hotter than when the wall is colder.

Mechanics↗
NASA NTRS1986

Mass, heat and freshwater fluxes in the South Indian Ocean

Six hydrographic sections were used to examine the circulation and property fluxes in the South Indian Ocean from 10 to 32 deg S. The calculations were made by applying an inverse method to the data. In the interior of the South Indian Ocean, the geostrophic flow is generally northward. At 18 deg S, the northward interior mass flux is balanced by the southward Ekman mass flux at the surface, whereas at 32 deg S the northward interior mass flux is balanced by the southward mass flux of the Agulhas Current. There is a weak, southward mass flux of 6 x 10 to the 9th kg/s in the Mozambique Channel. The rate of water exchange between the Pacific Ocean and the Indian Ocean is dependent on the choice of the initial reference level used in the inverse calculation. The choice of 1500 m, the depth of the deep oxygen minimum, has led to a flux of water from the Pacific Ocean to the Indian Ocean at a rate of 6.6 x 10 to the 9th kg/s. Heat flux calculations indicate that the Indian Ocean is exporting heat to the rest of the world's oceans at a rate of -0.69 x 10 to the 15th W at 18 deg S and -0.25 x 10 to the 15th W at 32 deg S (negative values being southward).

Fu, Lee-Lueng↗
NASA NTRS2009

Three-Dimensional Structure and Energy Balance of a Coronal Mass Ejection

UVCS observed Doppler-shifted material of a partial halo coronal mass ejection (CME) on 2001 December 13. The observed ratio of [O VJ/O V] is a reliable density diagnostic important for assessing the state of the plasma. Earlier UVCS observations of CMEs found evidence that the ejected plasma is heated long after the eruption. This paper investigated the heating rates, which represent a significant fraction of the CME energy budget. The parameterized heating and radiative and adiabatic cooling have been used to evaluate the temperature evolution of the CME material with a time-dependent ionization state model. Continuous heating is required to match the UVCS observations. To match the O VI bright knots, a higher heating rate is required such that the heating energy is greater than the kinetic energy.

Lee, J.-Y.↗
NASA NTRS2020

Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves

Ocean-driven basal melting of Antarctica’s floating ice shelves accounts for about half of their mass loss in steady state, where gains in ice-shelf mass are balanced by losses. Ice-shelf thickness changes driven by varying basal melt rates modulate mass loss from the grounded ice sheet and its contribution to sea level, and the changing meltwater fluxes influence climate processes in the Southern Ocean. Existing continent-wide melt-rate datasets have no temporal variability, introducing uncertainties in sea level and climate projections. Here, we combine surface height data from satellite radar altimeters with satellite-derived ice velocities and a new model of firn-layer evolution to generate a high-resolution map of time-averaged (2010–2018) basal melt rates and time series (1994–2018) of meltwater fluxes for most ice shelves. Total basal meltwater flux in 1994 (1,090 ± 150 Gt/yr) was similar to the steady-state value (1,100 ± 60 Gt/yr), but increased to 1,570 ± 140 Gt/yr in 2009, followed by a decline to 1,160 ± 150 Gt/yr in 2018. For the four largest ‘cold-water’ ice shelves, we partition meltwater fluxes into deep and shallow sources to reveal distinct signatures of temporal variability, providing insights into climate forcing of basal melting and the impact of this melting on the Southern Ocean.

Susheel Adusumilli↗
DOE OSTI2022

Influence of the Soret effect on binary-species boundary layers at high pressure

Direct numerical simulations of binary-species temporal boundary layers at high pressure are performed. The main objective is to investigate the influence of the Soret effect on flow physics of binary-species boundary layers where the fluid has a uniform composition. The working fluid is a mixture of 25 % methane and 75 % nitrogen in mass fraction. Although the fluid composition is uniform at the initial condition, the mass fraction of methane increases near the wall when the wall temperature is hotter than the free stream temperature, whereas it decreases when the wall temperature is colder. The non-uniform mass fraction indicates that the uphill diffusion occurs near the wall. Investigation of fluctuations of the mass fraction reveals that the mass fraction fluctuates in the whole boundary layer, indicating that the uphill diffusion occurs even far from the wall. Examination of the species-mass diffusion balance for mean flow fields clarified that the Soret effect flux becomes large near the wall, and the large flux causes the non-uniform profile of the mass fraction near the wall.

Binary Diffusion Coefficient↗