Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
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↗
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↗
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↗
DOE OSTI2023

The role of oak species in long-term soil P loss in a humid river bottomland

There is a gap in our understanding of if and how bottomland forest type will affect long-term nutrient cycling and loss. This study aims to determine how different forests affect soil hydrologic variability and whole-soil P loss in a humid-subtropical setting. Here, we used replicate-sampling and measured soil physical, chemical, and mineralogical properties at 12 sites in two forest ecosystems, post oak (Quercus stellata) and cherry bark oak (Quercus pagoda) in Clarks River National Wildlife Refuge in Western Kentucky. We hypothesize that wetting–drying events in redox soils of forested bottomlands can cause positive feedback in whole-soil P loss. Trees with greater P demand (e.g., post oak) take up more water creating more frequent and pronounced episodes in swelling and shrinking of expandable clays. Repeated swelling and shrinking of clays occlude the surface of Fe-Mn oxides from further adsorption of P in acidic soil. This can lead to greater loss or plant uptake of available P. Our results show (i) a significant difference in mean whole-soil P loss between the oak species with more loss in soils underlying the post oak forest, and (ii) a difference in the total P found in the sap- and heartwood of the two oak species. Soil analysis reveals that the clay mineralogy of the post and cherry bark oak sites are similar, and thus, may only play a minor role in governing the whole-soil P loss difference. However, the leaf data analysis suggests that the post oak site could be P and nitrogen-limited, while the cherry bark is only nitrogen limited. Our study shows that differences in the oak forest ecosystem may affect the long-term balance in water and nutrient uptake and may alter the redistribution of nutrients in the canopy and the underlying soils.

58 GEOSCIENCES↗
DOE OSTI2020

Water and Ammonia Desorption from 5A Zeolite

Demonstrate adsorption and desorption characteristics of water and ammonia on 5A zeolite during regeneration conditions using a Thermogravimetric analyzer (TGA) on samples previously evaluated by residual gas analyzer (RGA)/mass balance. Molecular sieve (zeolite) adsorbents have the following characteristics: Crystalline alumino-silicates containing pores or 'cages' that adsorb water or other molecules. Polar compound are adsorbed with high loading even at very low concentrations of the fluid. A 5A molecular sieve was chosen for its ability to adsorb both water and ammonia (NH{sub 3}) from a carrier gas stream. Preliminary RGA/mass balance testing indicated: 5A zeolite captured both NH{sub 3} and water as anticipated based on molecular size; Reloading showed limited success after zeolite heated beyond manufacturer recommended operating conditions: Steam dealumination theorized cause; Tests were reprioritized to focus on NH{sub 3} removal and end of life potential holdup of NH{sub 3} and water on 5A zeolite material. TGA measures a sample's weight and volatile composition as it is heated or cooled in a furnace. TGA heats sample to release volatile materials, (change in mass is measured). MS detects/measures low levels of impurities. Perform TGA desorption of zeolite samples taken at the following conditions tested with RGA (RGA results not shown here): Initial Bakeout conditions: hydrogen carrier gas, slowly raise temp to approximately 300 deg. C until desired moisture content <<100 ppm is achieved. Water loading (re-loading): argon carrier gas through bubbler until breach is observed. Ammonia loading (re-loading): hydrogen carrier gas through system. Regenerations/End of life conditions: hydrogen carrier gas, raise temp to approximately 300 deg. C, then slowly to 450 deg. C (end of life). After end of life testing less than 4% of the total mass was residual ammonia or water at bakeout conditions. Peak water and ammonia removal occurred around 275 deg. C and 150 deg. C respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The fate of carbon utilized by the subterranean termite Reticulitermes flavipes

Abstract Subterranean termites are ecosystem engineers that modulate the flow of carbon from dead wood to the atmosphere and soil, yet their contributions to the latter pool are largely unaccounted for in carbon cycling models. The fate of C from wood utilized by Reticulitermes flavipes (Kollar) was determined using a reductionist design in a closed system with δ 13 C labeled wood as a stable isotope tracer. The percentage of wood‐based carbon in termite respiratory gases, tissues, and organic deposits (frass and construction materials) was measured for five colonies to budget wood‐C mass distributed into metabolic and behavioral pathways during a 160‐h incubation period. We found that termites emitted 42% of the C from wood as gas (largely as carbon dioxide), returned 40% to the environment as organic deposits (frass and construction materials), and retained 18% in their tissues (whole alimentary tracts and de‐gutted bodies). Our findings affirm that termites are a source of greenhouse gases but are also ecosystem engineers that return approximately half the C from dead wood as organic deposits into their surrounding environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The utility of Bayesian data reconciliation for separations

Data reconciliation methods for separation processes typically rely on classical statistical approaches to generate estimates of true mass flow rates from measurements. Knowledge regarding the uncertainty of these estimates has value in decision making, but is often not acquired. Bayesian approaches intrinsically quantify uncertainty; however, literature for Bayesian data reconciliation of separation processes is scarce. This publication outlines two Bayesian data reconciliation models and provides details for how the models were implemented for the BayesMassBal (V 1.0.0) software package written in R. To demonstrate the advantages of this approach for data reconciliation, the models were first applied to simulated data and then compared to a classical model through a Monte Carlo experiment. In this example, the Bayesian models were found to provide more accurate estimates of the simulated data, while also providing quantitative information on the estimate uncertainty. To demonstrate the use of the technique in a practical problem, the models were also applied to real data collected from a pilot-scale rare earth solvent extraction process. Here, this publication provides a small window into how Bayesian methods can be used for data reconciliation, but findings suggest Bayesian data reconciliation models for separation processes have distinct advantages over classical alternatives.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

An automated multi-component gas adsorption system (MC GAS)

The knowledge gap on adsorption of complex mixtures in the literature relative to single component data represents a persistent obstacle to developing accurate process models for adsorption separations. The collection of mixed gas adsorption data is an imminent need for improved understanding of the behavior of adsorbent systems in these diverse adsorption applications. Current approaches to understanding mixture adsorption using predictive theories based on pure component adsorption experiments often fail to capture the behavior of more complex, non-ideal systems. In this work, we present an automated volumetric instrument for the measurement of mixed gas adsorption isotherms. This instrument was validated by comparison to other in-house instruments and data available in the literature, and the binary adsorption measurements were found to be thermodynamically consistent. The automation of this instrument allows for rapid collection of high-quality mixture adsorption data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Previously unaccounted atmospheric mercury deposition in a midlatitude deciduous forest

Significance Direct measurements of atmospheric deposition of gaseous elemental mercury (GEM) over a temperate forest showed a pronounced annual deposition of 25.1 µg ⋅ m −2 , which dominated as a source of mercury. GEM deposition was five times greater than wet deposition and three times greater than litterfall deposition, which has been used as a proxy for GEM deposition until now. Measured GEM deposition is driven by combined plant GEM uptake and underlying forest floor GEM uptake. Global forests may be a much larger global GEM sink than currently assumed, which may explain high mercury levels in soils across forests. Forest mercury mobilizes via watershed runoff and bioaccumulates in aquatic biota, ultimately leading to mercury exposures in wildlife and humans.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Modeling scale-up of particle coating by atomic layer deposition

Atomic layer deposition (ALD) is a promising technique to functionalize particle surfaces for energy applications including energy storage, catalysis, and decarbonization. In this work, we present a set of models of ALD particle coating to explore the transition from lab scale to manufacturing. Our models encompass the main particle coating manufacturing approaches including rotary bed, fluidized bed, and continuously vibrating reactors. These models provide key metrics, such as throughput and precursor utilization, required to evaluate the scalability of ALD manufacturing approaches and their feasibility in the context of energy applications. Our results show that designs that force the precursor to flow through fluidized particles transition faster to a transport-limited regime where throughput is maximized. They also exhibit higher precursor utilization. In the context of continuous processes, our models indicate that it is possible to achieve self-extinguishing processes with almost 100% precursor utilization. A comparison with past experimental results of ALD in fluidized bed reactors shows excellent qualitative and quantitative agreement.

25 ENERGY STORAGE↗
DOE OSTI2023

Crystalline versus Glassy Nature of Iron Phosphate Waste Forms Subjected to Different Slow Cooling Curves

This report provides experimental details and associated results on the elemental distributions (compositions), crystalline/amorphous phase distributions, and microstructures for iron phosphate reference materials (DPF5-336) made with different cooling curves and starting from different temperatures (i.e., 1050°C and 1200°C). The primary goal of this work was to evaluate the properties of the DPF5- 336 reference material after it was melted and cooled at different rates to simulate the cooling profile of different types of canisters. Crystal fractions across the SCC#2, SCC#3, and SCC#5 samples were all similar ranging in crystal content of 35.66–41.14 mass% with amorphous fractions being the balance (64.34–58.86 mass%, respectively). This shows that the heat treatment profile did not seems to greatly affect either the total crystal content or the phase distributions. The quenched sample was almost completely amorphous (99.55 mass%) with very small peaks identified as Li 3 Fe 2 (PO 4 ) 3 whereas the SCC-treated materials showed peaks for what appear to be ten separate phases present in various concentrations. Another goal of these experiments was to see how the phase distribution affected chemical durability of these waste forms, but that information is being collected at Argonne National Laboratory and will be published in a separate report

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Development of Steady-State and Dynamic Mass and Energy Constrained Neural Networks for Distributed Chemical Systems Using Noisy Transient Data

The paper presents the development of algorithms for mass and energy constrained neural network models that can exactly conserve the overall mass and energy of distributed chemical process systems, even though the noisy transient data used for optimal model training violate the same. In contrast to approximately satisfying mass and energy balance constraints of a system by soft penalization of objective function, algorithms have been developed for solving equality-constrained nonlinear optimization problems, thus providing the guarantee of exactly satisfying the system mass and energy conservation laws. For developing dynamic mass-energy constrained network models for distributed systems, hybrid series and parallel dynamic-static neural networks have been leveraged. The developed algorithms for solving both the training and forward problems are validated using both steady-state and dynamic data in the presence of various noise characteristics. The developed data-driven algorithms are flexible to exactly satisfy mass and energy balance constraints for dynamic chemical processes if the system holdup information is available. The proposed network structures and algorithms are applied to the development of data-driven lumped and distributed models of an adiabatic superheater/reheater system, a nonisothermal continuous stirred tank reactor, as well as an electrically heated plug-flow reactor system where one form of energy gets transformed to another. It has been observed that the mass-energy constrained neural networks yield a root mean squared error of <1% with respect to the system truth for the case studies evaluated in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗