Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Batch”

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 829 records · Page 46

Intensified reactive extraction for the acid-catalyzed conversion of fructose to 5-hydroxymethyl furfural

Here we conduct an experimental study of fructose dehydration to 5-hydroxymethyl furfural (HMF) in a biphasic microreactor as a function of residence time, temperature, and sugar loading using methyl butyl iso-ketone (MIBK) and 2-pentanol as extracting solvents. We demonstrate stable operation with maximum HMF yields of 93% and 87% in the two solvents, respectively, at 200 °C and a residence time of 2 s for a 5 wt% fructose aqueous feed. We report the highest optimal HMF space–time yield of 60 kg/L-hr at 200 °C (10 – 2,500-fold higher than published reports). Unexpectedly, an optimum organic-to-water ratio exists that depends on the solvent. Notably, we observe experimentally an increased fructose rate and HMF yield well above the extraction thermodynamic limit and hypothesize that the solvent plays a dual role, that of an extractant to protect HMF from degradation and a modifier of the fundamental chemistry. We expose mass transfer limitations of microreactors at longer residence times and higher temperatures and provide reactivity maps for their design. We show that batch reactors are unfit for high throughput and distributed manufacturing where small, farm-based systems are necessary.

42 ENGINEERING↗

Rewiring yeast metabolism for producing 2,3-butanediol and two downstream applications: Techno-economic analysis and life cycle assessment of methyl ethyl ketone (MEK) and agricultural biostimulant production

Rising concerns for sustainability and global climate change have driven the development of sustainable pro- duction pathways for biofuels and chemicals from lignocellulosic biomass via integrated biological and chemical processes. We constructed an engineered Saccharomyces cerevisiae capable of producing 2,3-butanediol (2,3-BDO) from glucose without accumulating ethanol and glycerol, which hinder downstream processing of 2,3-BDO, through extensive metabolic reprogramming. Specifically, we introduced heterologous 2,3-BDO biosynthetic enzymes and deleted the major isozymes of ethanol and glycerol biosynthetic enzymes. In addition, we introduced an NAD + regenerating Pyruvate-Malate (PM) cycle and enhanced the NAD + regenerating capability of the PM cycle to resolve the redox imbalance from the deletion of ethanol and glycerol production pathways. The resulting engineered yeast produced 109.9 g/L of 2,3-BDO with a productivity of 1.0 g/L/h and a yield of 0.36 g/ g glucose in a fed-batch fermentation. We also conducted techno-economic analysis (TEA) and life cycle assessment (LCA) of the production of methyl ethyl ketone (MEK) through catalytic dehydration of 2,3-BDO. A TEA based on the experimental results indicated that the minimum product selling price (MPSP) was estimated to be $1.90/kg. Regarding cradle-to-grave LCA, 100-year global warming potential (GWP100) and fossil energy consumption (FEC) were found to be 0.37 kg CO 2 eq/kg and 3.1 MJ/kg, respectively. These results demonstrated the feasibility of cost-competitive and sustainable bio-based MEK production via yeast fermentation. In addition, we explored the possibility of using the fermentation broth containing 2,3-BDO as a biostimulant inducing drought tolerance in plants. As a result, the yeast 2,3-BDO fermentation broth can induce drought tolerance in Arabidopsis thaliana without a complicated purification process.

2,3-butanediol (2,3-BDO)↗

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulated moving bed-inspired method for continuous adsorptive denitrogenation of model fuel

Efficient, continuous routes for removing nitrogen-containing compounds from hydrothermal liquefaction-derived synthetic aviation fuel are needed to enable direct blending with conventional jet fuels. Here, we report a simulated moving bed-inspired process for adsorptive denitrogenation of a model fuel. Unlike conventional simulated moving bed systems, which are designed for sharp separations between similar solutes, this approach was run deliberately outside the classical separation region so that both pyridine and indole were removed together from the hydrocarbon stream. Alcohol solvents were used to regenerate the silica adsorbent, maintaining performance over extended operation and avoiding the downtime and energy demand associated with calcination. Under these conditions, the system demonstrates removal of more than 98% of nitrogen while cutting solvent use by 28% compared to batch operation. Classical modeling tools predicted column concentration profiles even in this nontraditional regime, suggesting a straightforward path to scaling. Together, these results motivate solvent-efficient, continuous denitrogenation strategies that could be integrated with biorefinery processes.

Adsorption↗

Mobilization of trace metals from caprock and formation rocks at the Illinois Basin – Decatur Project demonstration site under geological carbon dioxide sequestration conditions

One concern for geologic CO 2 sequestration is the potential leakage of CO 2 or CO 2 -saturated brine containing trace metals into overlying aquifers, which poses the risk of adversely affecting underground sources of drinking water. In this work, rock and brine samples collected in the Illinois Basin – Decatur Project (IBDP), a large-scale geological CO 2 sequestration demonstration project in Decatur, Illinois, USA, as well as synthetic brine samples were used to understand the mobilization of trace metals as a result of CO 2 –rock–brine interactions under IBDP-specific conditions (50 °C and 20.6 MPa). Rock sample characterization indicated that, at the IBDP site, trace metal concentrations were greatest in cap rock samples from the Eau Claire Formation, compared with those from formations above and below this Formation. The natural brine samples collected from the Mt. Simon Sandstone reservoir were highly reducing and saline, with trace metal concentrations up to 680 times greater than the U.S. Environmental Protection Agency-prescribed drinking water standards. Batch leaching experiments indicated that both trace metal mobilization from rock and immobilization from the brine occurred when high-pressure CO 2 was introduced into the rock-brine system. The amounts of metals mobilized from the rock generally accounted for <5% of the total metals in the rock, but for some metals, including Ni, Pb, and Tl, up to 63% of the metals in the rock were mobilized to the brine. Leaching of trace metals into synthetic brines was different from that into the natural brine. The results of this study provided current information on the trace element sources and will help in risk simulations and experimental design to further evaluate potential impact of leakage on groundwater quality, which is important for future GCS projects to consider for monitoring and containment assurance purposes.

58 GEOSCIENCES↗

Application of spectroscopic evidence to diffuse layer model (DLM) parameter estimation for cation adsorption onto ferrihydrite in single- and bi-solute systems

The diffuse layer model (DLM) is one of the most frequently used surface complexation models (SCM) for predicting cation adsorption onto ferrihydrite. However, the DLM with its existing database sometimes fails to predict cation adsorption, especially in multi-solute systems. Notably, the failure is often attributed to inaccurate selection of the surface complexes and/or underestimation of the site density for DLM formulation. DLM simulations of batch adsorption data conducted with cadmium, copper, lead, and zinc on ferrihydrite in CO2 free systems demonstrates the ability of the DLM to fit single-solute data reasonably well but highlights its limitations with respect to predicting bi-solute data. Modifications to the DLM were made in two steps. First, spectroscopic evidence was used to guide cation complexation reactions at the ferrihydrite surface. Second, surface characteristics (reactive site density, specific surface area, and acidity constants) of ferrihydrite were updated according to recent literature. The results from this study suggest that when spectroscopic verification is employed and surface site density is increased, a simple two-site DLM can predict cation adsorption on ferrihydrite in single-solute systems for pH data above 3.5 and in most of the bi-solute systems tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sorption and oxidation of Co(II) at the surface of birnessite: Impacts of aqueous Mn(II)

The sorption of aqueous Mn(II) (1 mM) and Co(II) (50 and 200 μM) onto hexagonal birnessite (0.1 g L –1 ) was studied under anoxic conditions at pH 6.5 and 7.5 in binary and ternary experiments using batch kinetic experiments and XRD, ATR-FTIR, and Co K-edge EXAFS analyses. In the binary systems, sorption of Co(II) was accompanied by partial oxidation to Co(III) yielding a mixture of corner-sharing Co(II) and edge-sharing Co(III) complexes at both pH values, while Mn(II)-birnessite interaction resulted in coordination of Mn(II/III) at layer vacancy sites at pH 6.5, and led to reductive transformation of birnessite into feitknechtite at pH 7.5. In the ternary systems, strong mutual interferences between Co(II) and Mn(II) co-sorbates reduced the rate and extent of sorption relative to the binary experiments. The introduction of Mn(II) aq into the high-Co systems (200 μM) halted the slow sorption of Co(II) that was attributed to the incorporation of Co(III) into layer vacancies, while Co(II) aq prevented Mn(II)-driven transformation of birnessite into feitknechtite at pH 7.5. In the low-Co system (50 μM), reductive transformation of birnessite by Mn(II) aq at pH 7.5 produced Co(II)-substituted feitknechtite, a conversion that was accompanied by reduction of sorbed Co(III) to Co(II) which was partially released to solution. The strong effects of co-sorption are attributed to the similarity in sorption mechanisms of Co(II) and Mn(II), which both sorb as a mixture of di- and trivalent species. Finally, the results of this work demonstrate that aqueous Mn(II) may significantly affect the reactivity of phyllomanganate sorbents towards dissolved Co(II), and therefore impact the speciation and solubility of this trace metal in anoxic and suboxic geochemical environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗

Chemical impacts of subsurface CO 2 and brine on shallow groundwater quality

Leakage from geologic CO 2 sequestration (GCS) sites to overlying shallow drinking water aquifers is a tangible risk. A primary purpose of this study is to assess the potential impacts of CO 2 leakage into a fresh-water aquifer with associated CO 2 -water-sediment interactions. Here, the study site is the Ogallala aquifer overlying an active demonstration-scale GCS site in north Texas, USA. Using the results of combined batch experiments and reactive transport simulations, we discuss the effects of salinity on potential trace metal release and the potential for groundwater quality recovery after leakage ceases.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic insights into Cr(VI) removal by a combination of zero-valent iron and pyrite

Recent studies have revealed that a combination of zero-valent iron (ZVI) and pyrite (FeS 2 ) can effectively remove (Cr(VI)) from water, but the reasons behind this synergistic effect are still unclear. Our batch experiments showed that dissolved oxygen (DO) is a critical factor in the improved removal of Cr(VI) by ZVI and pyrite. When 0.08 g/L pyrite was combined with 0.5 g/L ZVI in the presence of DO, total Cr was reduced from 10 mg/L to 0.02 mg/L within 6 h. Conversely, in the absence of DO, total Cr was only reduced to 5.6 mg/L. DO oxidation of pyrite produced protons that promote ZVI corrosion, and mixing pyrite with water creates dissolved sulfide, which also contributes to the improved removal of Cr(VI). Here, electron microscopy images and X-ray absorption near edge structure analyses revealed that the presence of dissolved sulfide led to the formation of ferrous sulfide precipitates on the ZVI surface, preventing the formation of a passivating layer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Targeted electrochemical reduction of carcinogenic N -nitrosamines from emission control systems within CO 2 capture plants

N-Nitrosamines are one of the environmentally significant byproducts from aqueous amine-based post-combustion carbon capture systems (CCS) due to their potential risk to human health. Safely mitigating nitrosamines before they are emitted from these CO 2 capture systems is therefore a key concern before widescale deployment of CCS can be used to address worldwide decarbonization goals. Electrochemical decomposition is one viable route to neutralize these harmful compounds. The circulating emission control waterwash system, commonly installed at the end of the flue gas treatment trains to minimize amine solvent emissions, plays an important role to capture N-nitrosamines and control their emission into the environment. The waterwash solution is the last point where these compounds can be properly neutralized before becoming an environmental hazard. In this study, the decomposition mechanisms of N-nitrosamines in a simulated CCS waterwash with residual alkanolamines was investigated using several laboratory-scale electrolyzers utilizing carbon xerogel (CX) electrodes. Hcell experiments revealed that N-nitrosamines were decomposed through a reduction reaction and are converted into their corresponding secondary amines thereby neutralizing their environmental impact. Batch-cell experiments statistically examined the kinetic models of N-nitrosamine removal by a combined adsorption and decomposition processes. The cathodic reduction of the N-nitrosamines statistically obeyed the first-order reaction model. Finally, a prototype flow-through reactor using an authentic waterwash was used to successfully target and decompose N-nitrosamines to below the detectable level without degrading the amine solvent compounds allowing them to be return to the CCS and lower the system operating costs. Furthermore, the developed electrolyzer was able to efficiently remove greater than 98% of N-nitrosamines from the waterwash solution without producing any additional environmentally harmful compounds and offers an effective and safe route to mitigate these compounds from CO 2 capture systems.

20 FOSSIL-FUELED POWER PLANTS↗

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗

Extraction chromatography of 225 Ac and lanthanides on N,N-dioctyldiglycolamic acid /1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide solvent impregnated resin

The alpha-emitter 225 Ac (t 1/2 = 9.92 d) is currently under development for targeted alpha-particle therapy of cancer, and accelerator production of 225 Ac via proton irradiation of thorium targets requires robust separations of 225 Ac from chemically similar fission product lanthanides. Additionally, the lanthanide elements represent critical components in modern technologies, and radiolanthanides such as 140 Nd (t 1/2 = 3.37 d) also have potential application in the field of nuclear medicine. The ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Bmim][NTf 2 ]), combined with the diglycolamide extractant, N,N-dioctyldiglycolamic acid (DODGAA), was adsorbed on macroporous resin support to produce a solvent impregnated resin (SIR) that was investigated for separations of 225 Ac and lanthanides. The equilibrium distribution coefficients (Kd) of the rare earth elements (Sc(III), Y(III), Ln(III)), 225 Ac(III), Th(IV), and U(VI) on the prepared DODGAA/[Bmim][NTf 2 ]-SIR were determined from batch adsorption experiments in HCl and HNO 3 media. The DODGAA/[Bmim][NTf2]-SIR exhibited preferential uptake of the heavier lanthanide elements while allowing for the separation of the lighter lanthanides. Column separations utilizing the DODGAA/[Bmim][NTf 2 ]-SIR were effective at separating the lighter lanthanides from each other, and separating 225 Ac from a mixture of lanthanides, 213 Bi, and 225 Ra without the need for additional complexing agents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a SnO 2 -based 44 Ti/ 44 Sc generator for medical applications

Towards application of 44 Sc for diagnostic nuclear medicine, a 44 Ti/ 44 Sc generator based on an inorganic resin has been evaluated. Unlike other radionuclide generators used for medical applications, the long-term retention of the parent 44 Ti is vital due to its long half life. In this work, tin dioxide (SnO 2 ), a robust inorganic-based resin, has been synthesized and used as the stationary phase for a 44 Ti/ 44 Sc generator. The sorption behavior of 44 Ti/ 44 Sc was tested on SnO 2 with varying acids, concentrations, and times. Preliminary batch study results showed >88 % 44 Ti retention to the resin at lower acid concentrations (0.05 M HNO 3 and 0.05 M HCl). A pilot generator was evaluated for a year, demonstrating 85.3 ± 2.8 % 44 Sc elution yields and 0.71 ± 0.14 % 44 Ti breakthrough in 5 M HNO 3 . Based on capacity studies, a 7.4 MBq (200 µCi) upscaled generator system was constructed for further evaluation of the SnO 2 resin stability and the efficacy of the eluted 44 Sc for radiolabeling. 44 Sc could be regularly eluted from this generator in 5 M HNO 3 with an overall average radiochemical yield 84.7 ± 9.5 %. Post-elution processing of the 44 Sc with DGA-normal resin removed all 44 Ti present and allowed for high 44 Sc-DOTA labeling yields of 94.2 ± 0.5 %. Overall, SnO 2 has been shown to be a viable material for a 44 Ti/ 44 Sc generator.

07 ISOTOPE AND RADIATION SOURCES↗

Phyllosilicate mineral dissolution upon alkaline treatment under aerobic and anaerobic conditions

The dissolution of phyllosilicate minerals exposed to high-pH environments has been studied to determine the influence of alkaline treatments and variable redox conditions on clay dissolution and the potential formation of new phases via incongruent dissolution phenomena. The present study focused on the dissolution of phyllosilicate minerals (illite, muscovite, and montmorillonite) under anaerobic and aerobic conditions using comparative solutions (sodium hydroxide and ammonium hydroxide) at similar ionic strength in the presence and absence of oxygen. Our batch data show that there is a rapid decrease in aluminum dissolution (< 240 hrs) and slow increase in silica over time (up to 1440 hrs). This trend was particularly evident for montmorillonite which shows the greatest dissolution for ammonium hydroxide, likely due to intercalation of the polyatomic cation ammonium into the mineral’s expandable layers. When comparing alkaline treatments, the strong base, sodium hydroxide, dissolved more of the mica minerals, illite and muscovite, likely due to ion-pairing between the silicate tetrahedra [SiO4]n- and Na+ cations in solution. In addition, redox conditions affect the systems similarly with treatment, although the sodium hydroxide treatment results in greater variability in solution redox conditions. For all investigated phyllosilicates, the calculated aqueous aluminum and silicon ratios over time are significantly different from the minerals’ stoichiometric ratios. As a result, we conclude that incongruent dissolution occurred and suggest formation of secondary precipitates. Understanding the potential for clay mineral alterations from interaction with alkaline solutions has implications for in situ remediation, mining operations, and waste interactions with the subsurface. This research shows that incongruent dissolution is likely and will lead to secondary precipitation which may have long term physical and chemical impacts in the subsurface.

mineral dissolution, phyllosilicate, alkaline, amm↗

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗

Bayesian sequential optimal experimental design for nonlinear models using policy gradient reinforcement learning

We present a mathematical framework and computational methods for optimally designing a finite sequence of experiments. This sequential optimal experimental design (sOED) problem is formulated as a finite-horizon partially observable Markov decision process (POMDP) under a Bayesian setting and with information-theoretic utilities. The formulation is general and may accommodate continuous random variables, non-Gaussian posteriors, and nonlinear forward models. The sOED design policy incorporates elements of feedback and lookahead simultaneously, and we show it to generalize the commonly-used batch and greedy design strategies. We solve for the sOED policy using the policy gradient (PG) method from reinforcement learning, and provide a derivation for the PG expression in the sOED context. Adopting an actor-critic approach, the policy and value functions are parameterized using deep neural networks and improved via PG estimates produced from simulated episodes of designs and observations. The new PG-sOED algorithm is first validated on a linear-Gaussian benchmark, and then compared against other design baselines on a sensor movement problem for contaminant source inversion in a convection-diffusion field. As a result, we provide explanation for the policy behaviors using knowledge of the underlying physical process.

97 MATHEMATICS AND COMPUTING↗