SEARCH · Engineering Papers

Results for “dynamical balances”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. Here, we find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Reactive Capture and Conversion of Carbon Dioxide to Methanol with ZnZrO 2 and Alkali-Promoted Mg 3 AlO x Mixed Oxide Catalytic Sorbents

Reactive capture and conversion (RCC) explores the use of a single-unit process to capture CO 2 and produce a product, in this case, methanol (MeOH). In this study, different configurations of a catalytic sorbent (CS) composed of ZnZrO 2 catalyst and Mg 3 AlO x sorbent with and without alkali modification are evaluated for CO 2 adsorption, steady-state catalysis with cofed CO 2 and H 2 , and transient RCC performance. A catalyst composed of a physical mixture of Mg 3 AlO x with ZnZrO 2 resulted in a slight increase in CO 2 uptake, with a low impact on the catalytic activity and RCC of the materials compared to ZnZrO 2 alone. In contrast, Na impregnation significantly increased the level of CO 2 uptake from 0.28 mmol/g (ZnZrO 2 alone) to 0.6 and 1.1 mmol/g for the CS with Na on the catalyst or Mg 3 AlO x , respectively. However, Na impregnation reduced the CO 2 conversion rate and MeOH selectivity during steady-state cofeed experiments at 300 °C and 6 bar. In contrast to steady-state catalysis conditions, RCC, which is a cyclic capture and conversion process, creates dynamic CO 2 and H 2 surface coverages, favoring CH 4 in the early stages of the conversion step and then CO and MeOH as the catalyst CO 2 coverage reduces. The highest MeOH productivity during RCC was achieved with CS that balanced the CO 2 uptake with only moderate catalyst rate reductions caused by Na addition. The optimal material, ZnZrO 2 +10%Na/Mg 3 AlO x , achieved a CO 2 uptake of 0.8 mmol/g and a MeOH productivity of 0.5 mmol/g with 100% selectivity at 260 °C and 6 bar during RCC. This marks the highest RCC MeOH productivity reported to date, although the process needs further optimization and even with optimization, may remain impractical. The results further demonstrate that optimization of catalytic sorbents under steady-state flow conditions does not easily correlate to transient capture and conversion cycles for methanol synthesis from CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Analysis of Transformer Magnetizing Reactance Using Describing Function in Direct Power Control of Back-to-Back Modular Multilevel Converter with Advanced Grid Support

This paper presents a Port Controlled Hamiltonian based direct power control architecture for a back-to-back modular multilevel converter system connecting two ac sources at different frequencies. The system features advanced grid support functionalities based on IEEE 1547-2018, implemented on the inverter side. The rectifier side controller ensures reference following for active and/or reactive powers and maintains the commanded dc bus voltage. The inverter side controller ensures power command following for active and reactive powers.The proposed control architecture is designed to suppress second harmonic oscillations in powers during unbalanced grid voltage sags by dynamically adjusting the currents on each ac side. This also ensures effective elimination of any second harmonic oscillations in the equivalent dc bus voltage. Validation is performed on an OPAL-RT real-time platform with case studies on unbalanced and balanced sags, demonstrating the controller's effectiveness during real time implementation. A reduced-scale laboratory prototype further verifies these case studies, with experimental results for balanced sags due to grid simulator limitations. The results confirm the robustness and efficiency of the proposed control strategy in ensuring stable and reliable operation under various grid conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Development of Electro-chemical Battery Model for Plug-and-Play Eco-system Library

Energy storage components are fundamental to the concept of an Integrated Energy System (IES). They serve to store surplus energy during low-demand periods for later release when other IES components (i.e., Secondary Energy Source, Balance of Plant, etc.) would otherwise have to operate flexibly. This provision for storage avoids high-amplitude power ramps in these components thereby limiting thermal and mechanical stresses to their internals and providing for extended service life. This report describes a dynamic model that has been developed for an electrochemical battery. The lithium-ion (Li-ion) cell was selected as representative technology. The battery model was developed in the Dymola simulation environment and meets the requirements of the ecosystem plug-and-play library. The model accurately describes the electric dynamic response of a Li-ion battery for an imposed charging/discharging power profile. The corresponding physical limitations related to over-power scenarios, and the impact of the residual state of charge are accounted for in the model. A literature review of the major degradation processes affecting Li-ion batteries was performed. Given the purposes of the CTD-IES project, the progressive fade of the installed capacity, the reduction of the round-trip efficiency, and the limits on the number of charging/discharging cycles are aspects that need to be taken into account in techno-economic analyses. The modeling of these degradation phenomena becomes crucial when predictions over long time horizons (capacity expansion) are made. For each one of these phenomena, a brief description is given, and some figures to be implemented in the HERON optimization algorithm are presented.

25 ENERGY STORAGE↗

DOE OSTI2025

High-Fidelity Building Emulator for Integrated Comfort and Energy Analysis using EnergyPlus and Radiance

The growing need for smart, energy-efficient, and occupant-centric buildings has created a demand for advanced control systems that can optimize building operations to balance energy savings, demand flexibility, and comfort. However, current building energy simulation tools, such as EnergyPlus, have limitations that hinder the development and evaluation of these complex control systems. To address this challenge, we introduce a high-fidelity building emulator that dynamically couples EnergyPlus with Radiance for enhanced daylight modeling. The introduced workflow allows researchers and practitioners to rapidly develop and evaluate innovative control solutions. An example study looking at a south-facing office zone revealed up to 67% deviation in predicted light levels, which can significantly impact building assessment.

Yu, Tammie↗

DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗

DOE OSTI2022

Model files for estimating snow dynamics and stable water isotopes across the East River, CO.

A coupled hydrologic and snowpack stable water isotope model is used to assesses controls on isotopic inputs across the East River, Colorado, a large, mountainous basin. The hydrologic model uses the semi-empirical, spatially distributed and publicly available U.S. Geological Survey numerical code Precipitation-Modelling Runoff System (PRMS). Water and energy are tracked daily through the atmosphere, canopy and subsurface at a 100-m grid resolution. The isotope mass balance model follows previous work by Ala-aho et al. (2017) to track stable isotopes entering the soil system as snowmelt or rain. Water stores and fluxes needed for the isotope model use hydrologic model output for each timestep and model modeled grid location. This data package provides all files related to the hydrologic model (executable, input and output files), the isotopic model calibration and the historic isotopic model (water years 2015-2020). The isotope model source code and executable are provided. The file "readme.txt" includes information on the file structure, steps to run the model, and description of included folders.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Microreactor Automated Control System - Digital Twin Models and Advanced Control Systems Updates

Automation of control systems is expected to be important in the economic and safe operation of microreactors. Therefore, there is a need to develop and demonstrate automated control for microreactors, along with the development of testbeds for this purpose. This report provides updates on the status of a nonnuclear microreactor automated control system (MACS)—a real-time, hardware-in-the-loop testbed for non-nuclear testing of microreactor control system automation. A real-time hardware-in-the-loop testbed incorporates the realistic dynamics of physical systems into control system development and testing. The collaborative effort between Oak Ridge National Laboratory (ORNL) and Idaho National Laboratory (INL) resulted in the development of a prototypic microreactor plant-level digital twin that includes the reactor and a balance of plant system. Advanced control strategies were incorporated to demonstrate testing of control automation solutions. The gRPC communication protocol, which was implemented in the hardware-in-the-loop testbed by INL, was coupled to a digital twin model developed using the TRANsient Simulation Framework of Reconfigurable Models (TRANSFORM) library in Modelica. This digital twin simulation was tested with the ViBRANT hardware for realistic feedback and visual representation of control action in real time. A modular Python client structure was developed to manage functional mock-up unit-based simulation and real-time gRPC communication. Hardware-in-the-loop testing indicated that the modeled reactor—a natural-convection, molten-salt coolant loop configuration—responds well to control of drum positioning for modulation of reactor core power, as well as system-level control and downstream demand changes. Ongoing research is focused on integrating additional control algorithms that utilize data from newly included sensors within the MACS hardware testbed, as well as demonstrating and assessing the performance of the different automated control algorithms on multiple additional operational scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Synthetic Electricity Market Data Generation and HERON Use Case Setup of Advanced Nuclear Reactors Coupled with Thermal Energy Storage Systems

This study evaluates and optimizes advanced nuclear reactors coupled with thermal energy storage (TES) systems in an Integrated Energy System (IES) architecture to enable advanced nuclear power plants (A NPP) to participate in multi-commodity markets, thus enhancing their economic competitiveness. Nuclear-TES coupling scenarios studied herein are designed attenuate the nuclear heat dynamics and defer energy delivery to a later time, enabling the nuclear reactor to continue operating at or near steady-state design conditions as usual while also enabling flexible generation. Three A-NPPs, namely, an advanced light-water reactor (A LWR), a high temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating the technoeconomic of thermally balanced energy storage coupling design for thermal power extraction. Each of the reactor technologies were evaluated in two different electricity markets. Stochastic optimization approach was adopted which included the evaluation of price signals from the Pennsylvania-New Jersey-Maryland (PJM) market, and Electric Reliability Council of Texas (ERCOT), using an autoregressive moving average (ARMA) model. Risk Analysis Virtual Environment (RAVEN) tool and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON), were used to perform dispatch and capacity optimization, using the price data provided by the ARMA models. The results from the Nuclear-TES use cases will be used to design and characterize dynamic integrated system behavior and feedback.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Influence of Steric and Electronic Properties of P2 Groups on Covalent Inhibitor Binding to SARS-CoV-2 Main Protease

The main protease (MPro) of SARS-CoV-2 is a critical enzyme required for viral replication, making it a prime target for antiviral drug development. Covalent inhibitors, which form a stable interaction with the catalytic C145, have demonstrated strong inhibition of MPro, but the influence of steric and electronic properties of P2 substituents, designed to engage the S2 substrate-binding subsite within the MPro active site, on inhibitor binding affinity remains underexplored. Here, in this study, we design and characterize two hybrid covalent inhibitors, BBH-3 and BBH-4, and present their X-ray crystallographic structures in complex with MPro, providing molecular insights into how their distinct P2 groups, a dichlorobenzyl moiety in BBH-3 and an adamantyl substituent in BBH-4, affect binding conformation and active site adaptability. Comparative structural analyses with previously characterized inhibitors, including BBH-2 and Mcule-5948770040, reveal how the P2 bulkiness and electronic properties influence active site dynamics, particularly through interactions with the S2 and S5 subsites. The P2 group of BBH-3 induces conformational shifts in the S2 helix and the S5 loop, while BBH-4 displaces M49, stabilizing its binding through hydrophobic interactions. Isothermal titration calorimetry further elucidates the impact of P2 modifications on inhibitor affinity, revealing a delicate balance between enthalpic and entropic contributions. The data demonstrate that BBH-3 exhibits less favorable binding, affirming that dichlorobenzyl substitution at the P2 position has a more negative impact on the affinity for MPro than bulky saturated cyclic groups. This underscores the feature that MPro active site malleability may be accompanied by a conformational strain.

SARS-CoV-2↗

DOE OSTI2024

Nutrient Dynamics in a Coupled Terrestrial Biosphere and Land Model (ELM-FATES-CNP)

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Quantitative probing of the fast particle motion during the solidification of battery electrodes

Synchrotron x-ray tomography represents one of the most powerful imaging techniques that can visualize the three-dimensional microstructural complexity with a high spatial resolution. To date, however, the temporal resolution of this technique has been handicapped by the intrinsic limitations in the data acquisition speed and the lack of a suitable data processing method. The present study aims to fundamentally resolve this challenge by implementing a data analysis method that is based on differentiating discrete tomography reconstruction from a continuous scan reconstruction. With this methodological development, we capture and quantify the particle motion during the solidification of battery electrodes and reveal the statistics of the dynamically evolving motion in the drying process, which has been challenging to resolve. We discover that the particle motion exhibits a strong dependence on its geometric location within a drying electrode. Our results also imply that the final electrode quality can be controlled by balancing the solvent evaporation rate and the particle mobility in the region close to the drying surface.

25 ENERGY STORAGE↗

DOE OSTI2020

Comparison of aircraft measurements during GoAmazon2014/5 and ACRIDICON-CHUVA

The indirect effect of atmospheric aerosol particles on the Earth's radiation balance remains one of the most uncertain components affecting climate change throughout the industrial period. The large uncertainty is partly due to the incomplete understanding of aerosol–cloud interactions. One objective of the GoAmazon2014/5 and the ACRIDICON (Aerosol, Cloud, Precipitation, and Radiation Interactions and Dynamics of Convective Cloud Systems)-CHUVA (Cloud Processes of the Main Precipitation Systems in Brazil) projects was to understand the influence of emissions from the tropical megacity of Manaus (Brazil) on the surrounding atmospheric environment of the rainforest and to investigate its role in the life cycle of convective clouds. During one of the intensive observation periods (IOPs) in the dry season from 1 September to 10 October 2014, comprehensive measurements of trace gases and aerosol properties were carried out at several ground sites. In a coordinated way, the advanced suites of sophisticated in situ instruments were deployed aboard both the US Department of Energy Gulfstream-1 (G1) aircraft and the German High Altitude and Long-Range Research Aircraft (HALO) during three coordinated flights on 9 and 21 September and 1 October. Here, we report on the comparison of measurements collected by the two aircraft during these three flights. Such comparisons are challenging but essential for assessing the data quality from the individual platforms and quantifying their uncertainty sources. Similar instruments mounted on the G1 and HALO collected vertical profile measurements of aerosol particle number concentrations and size distribution, cloud condensation nuclei concentrations, ozone and carbon monoxide mixing ratios, cloud droplet size distributions, and downward solar irradiance. We find that the above measurements from the two aircraft agreed within the measurement uncertainties. The relative fraction of the aerosol chemical composition measured by instruments on HALO agreed with the corresponding G1 data, although the total mass loadings only have a good agreement at high altitudes. Furthermore, possible causes of the discrepancies between measurements on the G1 and HALO are examined in this paper. Based on these results, criteria for meaningful aircraft measurement comparisons are discussed.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Dual Energy Electron Storage Ring Cooler Design for Relativistic Ion Beams

Collider experiments demand small beam emittances in order to achieve high luminosity. For light particles such as electrons, there exists a natural synchrotron radiation damping resulting in low emittance beams at equilibrium. In the case of heavy particle beams such as proton or ion beams, there is no significant synchrotron radiation damping effect and some cooling mechanism is needed to get to low emittance beams. A dual energy storage ring cooler is a novel concept proposed to cool hadron beams at higher energies. The design consists of two rings: a low energy ring and a high energy ring connected by the energy recovering linac (ERL) that provides the necessary energy difference. The low energy ring has an electron cooler where the cooling interaction takes place between the electron and hadron beams whereas electron beam going through the high energy ring undergoes synchrotron radiation damping. In this document, we present a possible design of a dual energy storage ring-based electron cooler for high energy proton beam cooling suitable for the Electron Ion Collider (EIC) to be built at Brookhaven National laboratory (BNL). A special feature of a dual energy storage ring cooler design is that the electron beam energy in the low energy section must be tuned to match the hadron beam velocity in the cooling section, and the electron beam energy in the high energy section is chosen to provide an adequate synchrotron radiation damping. In addition, the ring optics design is carried out considering chromaticity correction, dynamic aperture, momentum aperture, beam lifetime, radiation damping and the intrabeam scattering effect. Finally, the cooling performance is simulated for 275 GeV proton beams at the top energy of the EIC. In such a collider it is desirable to cool the hadron beams to balance the emittance growth rates due to intra-beam scattering (IBS) to maintain a high luminosity during the collision runs. In such a situation, a dual energy storage ring based cooler could be a practical approach to be used for hadron beam cooling.

Dhital, Bhawin↗

DOE OSTI2023

Balancing Impedance and Controllability in Response Reconstruction with $\mathrm{TS-IMMAT}$

In a prior work the authors proposed a variant on the Impedance Matched Multi-Axis Test (IMMAT) in which a fixture is defined, called the Transmission Simulator (TS), and the desired environment is matched at a set of sensors on the TS. If the motion of the TS is matched then the response of the rest of the component will also match, provided that the attached component has the same dynamics as it did when the environment was measured. Hence, one would like the TS to be flexible so that it reproduces the boundary conditions that the component of interest experiences during flight, but the more flexible the TS, the more shakers might be needed to control its response. This work presents a derivation that gives expressions for these two potential error sources in TS-IMMAT. Then, various case studies are presented, both on simulated and real hardware, to understand the importance of each error term in practical testing. The theory explains the phenomena that were observed when using measurements from a component that flew on a sounding rocket. The environmental response was measured and then various fixtures were attached, each comprising more of the next assembly, or the hardware to which the component was attached in flight. MIMO testing was repeated with each fixture and the results were compared to seek to understand the role of the impedance match in this type of testing. Here the results show that the number of modes that are active in the transmission simulator is also very important, and so the best solution balances these two considerations. An improved method of simulating the MIMO test is then proposed, so simulations can be used to predict what fixture, or transmission simulator, will give the best results in a TS-IMMAT test.

42 ENGINEERING↗

DOE OSTI2020

Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

58 GEOSCIENCES↗