Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate optimization”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Retracted Article: Distributed energy management for networked microgrids in a three-phase unbalanced distribution network

Owing to increased penetration of three-phase and single-phase microgrids, distributed energy resources (DERs), and responsive loads, the maintenance of a three-phase balance by distribution networks is a significant challenge. Existing literature on distributed energy management for networked microgrids generally neglects the distribution network or employs a simplified phase balanced distribution network; thus, these evaluations are not applicable. Further, the underlying mutual coupling between the different phases of distribution feeders results in a more challenging situation. Here, to solve this issue, this study sought to propose distributed energy management based on a three-phase unbalanced distribution network. Various three-phase or single-phase microgrids, utility-owned DERs, and responsive loads were coordinated through iteratively adjusted price signals. Based on the price signals received, the microgrid controllers (MCs) and distribution management system (DMS) updated the schedules of the DERs and responsive loads under their jurisdiction separately. The price signals were then updated according to the generation-load mismatch at each node and distributed to the corresponding MCs and DMS for the next iteration. The iteration continued until a sufficiently small generation-load mismatch was achieved at all nodes, that is, a balanced generation and load at all nodes under the agreed price signals. Considering a three-phase unbalanced distribution network, the price signals were determined per phase per node. Overall, the proposed distributed energy management coordinates microgrids, utility-owned DERs, responsive loads with guaranteed network constraints, and preserves the privacy of microgrid customers. This distributed energy management method was further demonstrated through various case studies on a three-phase networked microgrid test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Evaluating the Impact of Managed EV Charging for Reliable Operation of Bulk Power Systems with High Non-Dispatchable Generation

The growth of electric vehicles (EVs) and variable-generation (VG) sources introduces new challenges for power-system operations. This study introduces a modeling framework and evaluates five EV charging strategies under projected 2040 grid conditions in the Evergy service territory with high non-dispatchable generation. Using realistic EV behavior and generation models, their impacts on system peak demand, ramp rate, and reserve capacity are evaluated. Results show that only the peak-avoidance strategy effectively reduces system peak demand, while decentralized strategies-particularly cost based dynamic charging-can exacerbate peaks due to synchronized user behavior. However, ramp-rate minimization strategy significantly reduce the stress on dispatchable generation achieving the lowest maximum absolute ramp rate (MARR) (56.81 MW) and lowest reserve requirement (2.39 GW). In contrast, unmanaged and TOU random strategies increase the stress on dispatchable generation sources with increased MARR and reserve requirements. These findings highlight the importance of coordinated, system-aware managed charging strategies to ensure reliable and affordable grid operation in the presence of EVs and VG sources.

14 - SOLAR ENERGY↗
NASA NTRS2000

Validation of a High-Order Prefactored Compact Scheme on Nonlinear Flows with Complex Geometries

Three benchmark problems are solved using a sixth-order prefactored compact scheme employing an explicit 10th-order filter with optimized fourth-order Runge-Kutta time stepping. The problems solved are the following: (1) propagation of sound waves through a transonic nozzle; (2) shock-sound interaction; and (3) single airfoil gust response. In the first two problems, the spatial accuracy of the scheme is tested on a stretched grid, and the effectiveness of boundary conditions is shown. The solution stability and accuracy near a shock discontinuity is shown as well. Also, 1-D nonlinear characteristic boundary conditions will be evaluated. In the third problem, a nonlinear Euler solver will be used that solves the equations in generalized curvilinear coordinates using the chain rule transformation. This work, continuing earlier work on flat-plate cascades and Joukowski airfoils, will focus mainly on the effect of the grid and boundary conditions on the accuracy of the solution. The grids were generated using a commercially available grid generator, GridPro/az3000.

Hixon, Ray↗
DOE OSTI2021

Power Quality and Stability in a Cluster of Microgrids with Coordinated Power and Energy Management

Distributed Energy Resources (DER) such as photovoltaic (PV) systems and battery energy storage systems (BESS) can be operated collectively as microgrids. Microgrids can be effective in meeting local load requirements as well as in improving the power quality and stability in a modern power distribution system. Multiple microgrids can also be operated in a coordinated manner as a cluster, to improve the resiliency of the power distribution system. In the event of an outage caused by a transmission system failure, the microgrids in a cluster can use their distributed generation capacity and energy storage resources to recover and extend power availability to the critical loads in the system. In this paper, an illustrative cluster of two microgrids based on the IEEE 13-bus model with DERs has been used, to demonstrate the operation of the concept using real-time simulations. Using real-time simulations, capabilities such as switching reconfiguration following faults, identification of optimal DER placement for effective power quality management, and power electronic controller validation have been demonstrated. For the two-microgrid cluster, this paper also presents improvements in load voltage quality and

Chowdhury, Prithwiraj Roy↗
DOE OSTI2024

Optimized Substrate Positioning Enables Switches in the C–H Cleavage Site and Reaction Outcome in the Hydroxylation–Epoxidation Sequence Catalyzed by Hyoscyamine 6β-Hydroxylase

Hyoscyamine 6β-hydroxylase (H6H) is an Fe(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase that catalyzes the last two steps in the biosynthesis of scopolamine, a prolifically administered anti-nausea drug. After its namesake first reaction, H6H couples the newly installed C6-bonded oxygen to C7 to form the epoxide of scopolamine. Oxoiron(IV) (ferryl) intermediates initiate both reactions by cleaving C–H bonds, but it remains unclear how the enzyme switches target site and promotes (C6)O–C7 coupling in preference to C7 hydroxylation in the second step. In one possible epoxidation mechanism, the C6 oxygen would – analogously to mechanisms proposed for the Fe/2OG halogenases and, in the preceding paper, N-acetylnorloline synthase (LolO) – coordinate as alkoxide to the C7–H-cleaving ferryl intermediate to enable alkoxyl coupling to the ensuing C7 radical. Here we provide structural and kinetic evidence that H6H instead exploits the distinct spatial dependencies of competitive C–H-cleavage (C6 vs C7) and C–O-coupling (oxygen rebound vs cyclization) steps to promote the two-step sequence without substrate coordination or repositioning for the epoxidation step. Structural comparisons of ferryl-mimicking vanadyl complexes of wild-type H6H and a variant that preferentially hydroxylates C7 of 6-hydroxyhyoscyamine suggest that only a modest (~ 10°) shift in the Fe–O–H(C7) approach angle is sufficient to change the outcome. Finally, the observation that, in wild-type H6H, 2 H 2 O solvent also increases the C7-hydroxylation:epoxidation ratio by ~ 8-fold implies that the latter outcome requires cleavage of the alcohol O-H bond, which, unlike in the LolO oxacyclization, is not accomplished in advance of C–H cleavage.

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

Multi-Agent Deep Reinforcement Learning for Realistic Distribution System Voltage Control Using PV Inverters

Over the last few decades, the deployment of distributed solar photovoltaic (PV) systems has increased consistently. High PV penetration could cause adverse effects on the grid, such as voltage violations. This paper proposes a new distributed soft actor-critic based multi-agent deep reinforcement learning (SAC-MADRL) control solution to minimize the PV real power curtailment while keeping the grid voltage in an acceptable range. New reward functions have been designed to coordinate different agents during the learning process, yielding improved convergence. Comparison results with other control methods on a real feeder in western Colorado U.S. with 80% penetration of PVs demonstrate that the proposed method has better capability of effectively regulating voltage while minimizing the PV real power curtailment.

distribution system↗
DOE OSTI2025

AI Applications to Physics Experiments at Jefferson Lab

We survey how AI/ML is being deployed across Jefferson Lab's experimental and accelerator programs. In EPSCI, Hydra applies computer vision to automate real-time data-quality monitoring across all four experimental halls, replacing manual inspection of hundreds to thousands of histograms per shift. AIEC (AI Experiment Controls) uses ML to stabilize drift chamber gains and is now part of standard CEBAF production running, while AI Optimized Polarization (AIOP) targets autonomous control of polarized targets and photon beam angular alignment. In CASA, cavity fault classification models identify faulted cavities and trip types from waveform data with ~85% and ~78% agreement to labeled data, respectively, and are deployed in production; a separate effort applies LLMs and hybrid search to make the CEBAF operations logbook AI-ready. QCD-focused work includes transformer- and GAN-based generative models for particle-level event simulation, with distributed GAN training scaling studies on Polaris. Additional efforts span ML-on-FPGA for the EIC and a new Data Science Department coordinating anomaly detection, uncertainty quantification, and HPC-scalable ML lab-wide. Collectively, these projects illustrate AI's growing role in improving efficiency across JLab's nuclear physics mission.

Mei, Xinxin [Thomas Jefferson National Accelerator↗
NASA NTRS2006

4-Vinyl-1,3-Dioxolane-2-One as an Additive for Li-Ion Cells

Electrolyte additive 4-vinyl-1,3-dioxolane-2-one has been found to be promising for rechargeable lithium-ion electrochemical cells. This and other additives, along with advanced electrolytes comprising solutions of LiPF6 in various mixtures of carbonate solvents, have been investigated in a continuing effort to improve the performances of rechargeable lithium-ion electrochemical cells, especially at low temperatures. In contrast to work by other researchers who have investigated the use of this additive to improve the high-temperature resilience of Li-ion cells, the current work involves the incorporation of 4-vinyl-1,3-dioxolane-2-one into quaternary carbonate electrolyte mixtures, previously optimized for low-temperature applications, resulting in improved low-temperature performance. The benefit afforded by 4-vinyl-1,3- dioxolane-2-one can be better understood in the light of relevant information from a number of prior NASA Tech Briefs articles about electrolytes and additives for such cells. To recapitulate: The loss of performance with decreasing temperature is attributable largely to a decrease of ionic conductivity and the increase in viscosity of the electrolyte. What is needed to extend the lower limit of operating temperature is a stable electrolyte solution with relatively small lowtemperature viscosity, a large electric permittivity, adequate coordination behavior, and appropriate ranges of solubilities of liquid and salt constituents. Whether the anode is made of graphitic or non-graphitic carbon, a film on the surface of the anode acts as a solid/electrolyte interface (SEI), the nature of which is critical to low-temperature performance. Desirably, the surface film should exert a chemically protective (passivating) effect on both the anode and the electrolyte, yet should remain conductive to lithium ions to facilitate intercalation and de-intercalation of the ions into and out of the carbon during discharging and charging, respectively. The additives investigated previously include alkyl pyrocarbonates. Those additives help to improve low-temperature performances by giving rise to the formation of SEIs having desired properties. The formation of the SEIs is believed to be facilitated by products (e.g., CO2) of the decomposition of these additives. These decomposition products are believed to react to form Li2CO3-based films on the carbon electrodes. The present additive, 4-vinyl-1,3-dioxolane-2-one, also helps to improve lowtemperature performance by contributing to the formation of SEIs having desired properties, but probably in a different manner: It is believed that, as part of the decomposition process, the compound polymerizes on the surfaces of carbon electrodes.

Smart, Marshall↗
NASA NTRSDate not supplied

Model Predictive Control Study for an Electrified Turbofan Engine

A control sensitivity study is conducted for electrified aircraft propulsion concepts using a publicly available turbofan engine simulation. A Proportional-Integral (PI) controller serves as a baseline for comparison against an advanced Model Predictive Control (MPC) approach. MPC is used in the control scheme to allow for subsystem coordination of engine and electric power loads. The electrified turbofan engine includes electric machines to augment the amount of torque being applied to the spools of the engine during transient phases of operation. Power augmentation allows for energy to be stored or exhausted in a manner that increases the fuel efficiency of the engine, reducing aircraft emissions and noise. The study compares the fan speed, high pressure compressor stall margin, and low pressure compressor stall margin against the baseline PI control system, an MPC system without electric machines, and an MPC system with electric machines. The MPC system with electric machines is optimized using the sample time, control horizon and prediction horizon. The MPC design provided a stable damped response as desired and allowed operability margins to be maintained

Electrified Aircraft Propulsion↗
DOE OSTI2026

A Survey on the Expanding Scope and Interdisciplinary Opportunities for Processing-in-Memory Techniques

Processing-in-Memory (PIM) is emerging as a practical path to overcome the limitations of traditional von Neumann architectures. At its core, PIM systems implement computing primitives such as logic operations and multiply-accumulate acceleration through compute-in-memory, near-memory processing, or hybrid designs. The role of memory cells varies widely across technologies, acting as inputs, outputs, or analog accumulators through bit-lines and sense amplifiers. This diversity creates trade-offs in precision, bandwidth, latency, and programmability, making it difficult to build a unified understanding on the progress of the field. In this survey, we organize recent advances of PIM into three areas. First, we discuss the progress on the architectural optimizations of PIM and its integration with both DRAM and emerging non-volatile memories. Second, we examine how PIM is being used to accelerate key computing domains, including generative AI workloads and high-performance kernels, along with new approaches. Third, we highlight the growing adoption of PIM in computational sciences, where it is being applied to solve interdisciplinary problems such as genome analysis, mRNA quantification, mass spectrometry, quantum circuit simulation, wave modeling, and secure computation. Finally, we synthesize the major challenges that continue to slow PIM adoption, including manufacturing constraints, power delivery, thermal reliability, data consistency, runtime and memory-management coordination, and the difficulty of building portable software abstractions without sacrificing commercial viability. This work provides an updated, structured perspective on PIM’s potential across computing and computational sciences and the barriers that must be solved for it to reach its full impact.

Asifuzzaman, Kazi [Oak Ridge National Laboratory (↗
DOE OSTI2025

Metal–Oxide Interface Sites Created Using Atomic Layer Deposition and Tested for CO Oxidation

The performance of catalysts made out of Pt supported on TiO 2 thin films grown on SBA-15 (a silica mesoporous material) by atomic layer deposition (ALD) was characterized systematically by combining in situ infrared absorption spectroscopy (IR) with other techniques including electron microscopy and adsorption−desorption isothermal measurements. The titania films in the resulting high-surface-area catalysts were evenly distributed throughout the inner surface of the SBA-15 mesopores, and their thickness could be controlled at a submonolayer level, with 3 to 4 TiO 2 ALD cycles needed for the complete coverage of the silica sites. The titania films could be deposited either before or after adding the metal (Pt), which was dispersed in the form of small nanoparticles (NPs) approximately 4−6 nm in diameter, in order to exert some control on the density and nature of the Pt/TiO 2 interface sites. One important lesson deriving from this work is that such an order of deposition leads to significantly different catalysts in spite of the fact that most of their structural properties are similar. If the Pt is deposited on the titania films, the resulting metal NPs are slightly smaller than those grown on silica and display CO adsorption sites with lower surface Pt coordination numbers. On the other hand, when TiO 2 is deposited on the Pt/SBA-15 starting material, some titania grows on the metal and partially blocks its surface while also creating new interface sites where CO binds more weakly and displays lower C−O stretching frequencies. In terms of catalytic performance, the results from in situ IR CO site titration and kinetic measurements combined suggest a mechanism where CO first adsorbs on Pt atop sites and then migrates to Pt/TiO 2 interface sites, where oxidation takes place. Both types of sites appear to be similar in all the catalysts tested, but catalytic performance could be optimized by tuning their surface densities. Maximum catalytic activity was obtained when the TiO2 films were deposited first and with TiO 2 coverages of at least half a monolayer, that is, after at least 2 ALD cycles.

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

Rational solvent molecule tuning for high-performance lithium metal battery electrolytes

Electrolyte engineering improved cycling of Li metal batteries and anode-free cells at low current densities; however, high-rate capability and tuning of ionic conduction in electrolytes are desirable yet less-studied. In this study, we design and synthesize a family of fluorinated-1,2-diethoxyethanes as electrolyte solvents. The position and amount of F atoms functionalized on 1,2-diethoxyethane were found to greatly affect electrolyte performance. Partially fluorinated, locally polar –CHF 2 is identified as the optimal group rather than fully fluorinated –CF 3 in common designs. Paired with 1.2 M lithium bis(fluorosulfonyl)imide, these developed single-salt-single-solvent electrolytes simultaneously enable high conductivity, low and stable overpotential, >99.5% Li||Cu half-cell efficiency (up to 99.9%, ±0.1% fluctuation) and fast activation (Li efficiency >99.3% within two cycles). Combined with high-voltage stability, these electrolytes achieve roughly 270 cycles in 50-μm-thin Li||high-loading-NMC811 full batteries and >140 cycles in fast-cycling Cu||microparticle-LiFePO 4 industrial pouch cells under realistic testing conditions. The correlation of Li + –solvent coordination, solvation environments and battery performance is investigated to understand structure–property relationships.

25 ENERGY STORAGE↗
NASA NTRS1978

A coordinated X-ray, optical, and microwave study of the flare star Proxima Centauri

Results are reported for a three-day coordinated observing program to monitor the flare star Proxima Centauri in the X-ray, optical, and radio spectrum. During this interval 30 optical flares and 12 possible radio bursts were observed. The SAS 3 X-ray satellite made no X-ray detections. An upper limit of 0.08 on the X-ray/optical luminosity ratio is derived for the brightest optical flare. The most sensitive of the radio telescopes failed to detect 6-cm emission during one major and three minor optical flares, and on this basis an upper limit on the flare radio emission (1 hundred-thousandth of the optimal luminosity) is derived.

Haisch, B. M.↗
DOE OSTI2024

Probing Basal and Prismatic Planes of Graphitic Materials for Metal Single Atom and Subnanometer Cluster Stabilization

Abstract Supported metal single atom catalysis is a dynamic research area in catalysis science combining the advantages of homogeneous and heterogeneous catalysis. Understanding the interactions between metal single atoms and the support constitutes a challenge facing the development of such catalysts, since these interactions are essential in optimizing the catalytic performance. For conventional carbon supports, two types of surfaces can contribute to single atom stabilization: the basal planes and the prismatic surface; both of which can be decorated by defects and surface oxygen groups. To date, most studies on carbon‐supported single atom catalysts focused on nitrogen‐doped carbons, which, unlike classic carbon materials, have a fairly well‐defined chemical environment. Herein we report the synthesis, characterization and modeling of rhodium single atom catalysts supported on carbon materials presenting distinct concentrations of surface oxygen groups and basal/prismatic surface area. The influence of these parameters on the speciation of the Rh species, their coordination and ultimately on their catalytic performance in hydrogenation and hydroformylation reactions is analyzed. The results obtained show that catalysis itself is an interesting tool for the fine characterization of these materials, for which the detection of small quantities of metal clusters remains a challenge, even when combining several cutting‐edge analytical methods.

Vidal, Mathieu↗
DOE OSTI2025

Synthesis of UO 2 nanoparticles via coulometric titration: influence of electrolytes and ligands on synthesized particle properties

Here, this study investigates the controlled synthesis of uranium oxide nanoparticles (UO 2 NPs) via coulometric titration under ambient conditions, focusing on the impact of electrolyte anions, salt concentrations, and strongly complexing ligands on particle formation and properties. Using dilute solutions of mineral acids (HNO 3 , HCl, and HClO 4 ), we demonstrate that the choice of electrolyte anion affects particle size, polydispersity, and surface charge. Particle characterization using dynamic light scattering and transmission electron microscopy shows that particles synthesized in perchlorate are largest and show a high degree of polydispersity. In comparison, particles synthesized in chloride are smaller and more uniform in size. Synthesis in nitrate yields a wide variety of particle sizes, with the major fraction (>65 %) having a smaller size than that obtained in the other two electrolytes. The presence of more strongly coordinating ligands such as sulfate and acetate modulate hydrolysis and condensation reactions, with sulfate suppressing nanoparticle formation across a wide concentration range ([SO 4 2- ] > 5 mM) and acetate enabling stable colloidal suspensions at [HOAc] ≤ 0.1 M. Synthesis in concentrated electrolytes (2 M NaNO 3 , NaCl, or NaClO 4 ) accelerates reaction kinetics but introduces challenges, as particles showed increased polydispersity and aggregation, and were more prone to oxidation. Electrolyte effects on actinide oxide nanoparticle formation are discussed and a short comparison to established nanoparticle syntheses is drawn. This work underscores the importance of tailoring synthesis parameters to achieve desired nanoparticle properties, providing valuable insights for optimizing UO 2 NP production for various applications.

Actinide hydrolysis↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2022

Consistent responses of vegetation gas exchange to elevated atmospheric CO 2 emerge from heuristic and optimization models

Elevated atmospheric CO 2 concentration is expected to increase leaf CO 2 assimilation rates, thus promoting plant growth and increasing leaf area. It also decreases stomatal conductance, allowing water savings, which have been hypothesized to drive large-scale greening, in particular in arid and semiarid climates. However, the increase in leaf area could reduce the benefits of elevated CO 2 concentration through soil water depletion. The net effect of elevated CO 2 on leaf- and canopy-level gas exchange remains uncertain. To address this question, we compare the outcomes of a heuristic model based on the Partitioning of Equilibrium Transpiration and Assimilation (PETA) hypothesis and three model variants based on stomatal optimization theory. Predicted relative changes in leaf- and canopy-level gas exchange rates are used as a metric of plant responses to changes in atmospheric CO 2 concentration. Both model approaches predict reductions in leaf-level transpiration rate due to decreased stomatal conductance under elevated CO 2 , but negligible (PETA) or no (optimization) changes in canopy-level transpiration due to the compensatory effect of increased leaf area. Leaf- and canopy-level CO 2 assimilation is predicted to increase, with an amplification of the CO 2 fertilization effect at the canopy level due to the enhanced leaf area. The expected increase in vapour pressure deficit (VPD) under warmer conditions is generally predicted to decrease the sensitivity of gas exchange to atmospheric CO 2 concentration in both models. The consistent predictions by different models that canopy-level transpiration varies little under elevated CO 2 due to combined stomatal conductance reduction and leaf area increase highlight the coordination of physiological and morphological characteristics in vegetation to maximize resource use (here water) under altered climatic conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Embedded training of neural-network subgrid-scale turbulence models

We report that the weights of a deep neural-network model are optimized in conjunction with the governing flow equations to provide a model for subgrid-scale stresses in a temporally developing plane turbulent jet at Reynolds number Re 0 = 6000 . The objective function for training is first based on the instantaneous filtered velocity fields from a corresponding direct numerical simulation, and the training is by a stochastic gradient descent method, which uses the adjoint Navier-Stokes equations to provide the end-to-end sensitivities of the model weights to the velocity fields. In-sample and out-of-sample testing on multiple dual-jet configurations show that its required mesh density in each coordinate direction for prediction of mean flow, Reynolds stresses, and spectra is half that needed by the dynamic Smagorinsky model for comparable accuracy. The same neural-network model trained directly to match filtered subgrid-scale stresses, without the constraint of being embedded within the flow equations during the training, fails to provide a qualitatively correct prediction. The coupled formulation is generalized to train based only on mean-flow and Reynolds stresses, which are more readily available in experiments. The mean-flow training provides a robust model, which is important, though a somewhat less accurate prediction for the same coarse meshes, as might be anticipated due to the reduced information available for training in this case. The anticipated advantage of the formulation is that the inclusion of resolved physics in the training increases its capacity to extrapolate. This is assessed for the case of passive scalar transport, for which it outperforms established models due to improved mixing predictions.

42 ENGINEERING↗