Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applications development system design”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Methods to Enhance Wellbore Cement Integrity with Microbially-Induced Calcite Precipitation (MICP) (Final Scientific/ Technical Report)

The goal of this project was to develop improved methods for sealing compromised wellbore cement in leaking oil and gas wells, thereby reducing the risk of unwanted upward fluid migration. Novel methods for improving wellbore integrity, such as microbially induced calcite precipitation (MICP), can reduce leakage potential, improve the safety of fossil fuel extraction, improve the public perception of hydraulic fracturing, and promote environmentally-prudent unconventional oil and gas development. Microbes, with the urease enzyme, can catalyze the chemical reaction of urea hydrolysis to induce the precipitation of calcium carbonate which can be used as a cementitious material to seal leakage pathways. In this project, methods to promote robust bio-composite cementitious materials were designed and tested in the laboratory. Scale-up of those methods were tested in meso-scale reactor systems and in field applications. In this report, in Section One, we describe laboratory efforts to develop injection strategies to promote precipitation in wellbore analogs and determine the strength of the bio-composite cements as compared to fine cement. In Section Two, we describe the efforts to scale up the work and study the use of materials that can be used in field application, for example exploring the use of calcium chloride ice melt or urea fertilizer as source chemicals. In Section Three, the three field trials (methods and results) performed as part of the project are described and summarized. At the end of the report is a comprehensive summary and conclusion section which highlights the key findings of the project. The work performed during this project significantly advanced the technology readiness level (TRL) of the MICP wellbore sealing strategy.

03 NATURAL GAS↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2020

Holographic wavefront sensors

A brief historical review of the first designs of holographic wavefront sensors (WFS’s) and the concepts lying in their basis is presented. The main directions in the development of these sensors are highlighted and considered. One of these directions implies a two-stage transition from the use of several separate measuring channels with holograms filtering only one Zernike mode: first to the use of one channel with a multiplexed hologram filtering several first Zernike modes, and then to filtering the entire set of Zernike modes with the aid of one combined multiplexed hologram. Another line of research in this field (related to the first one) is the optimisation of the filter hologram structure in order to reduce cross-modulation interferences, increase multiplexing level, etc. One more line of research implements principles of dynamic holography by introducing spatial light modulators into the WFS composition. Hence, the advantages of time multiplexing of holograms can be used. The approach developed by G. Andersen’s team, aimed at adapting the holographic WFS design for operation as an element of adaptive optical system with a zonal corrector, as well as an approach implying development of hybrid holographic sensors, are considered separately. The results of the authors’ studies of holographic WFS’s with application of the methods of Fourier holography (i.e., holography of focused beams, in particular, using diffuse scatterers in a hologram recording channel) are also reported. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Thermo-Mechanical Modeling and Evaluation of the Cracking Response of Additively Manufactured Monolithic SiC Lattice Structures Subjected to Laser Heating

Increasing operating temperatures of solar receivers is paramount to the efficiency of concentrated solar thermal (CST) and solar power (CSP) systems. Owing to its high temperature stability combined with excellent thermal and optical properties, SiC has been the material of choice for application in high-temperature solar receivers. The state-of-the-art SiC volumetric concentrating solar air receivers such as honeycomb design have been demonstrated in field tests to achieve exit air temperatures approaching 800oC. However, successful application of CST systems for decarbonization requires significant increase in temperature capability of SiC receiver technology. Supported by an award from the Solar Technology Office (SETO), US Department of Energy (DOE), GE Research in collaboration with Heliogen Holdings Inc, is engaged in the development of ultra-High Operating Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR) enabled by additive manufacturing. The program objective is to design and to demonstrate a techno-economically viable SiC air receiver technology to achieve exit air temperature of 1100oC for CST applications. This report summarizes the learnings of a computational study on additively manufactured SiC lattice structures subjected to laser heating. The results from this study aim to help guide the design of a SiC receiver through a better understanding of lattice structure geometric parameters and their implications on the cracking response of the structure. In the study, a heat flux was applied to the surface of a 2”-diameter cylindrical lattice structure to simulate a 4kW CO2 laser. The resulting temperature distribution was applied to a structural model to approximate the stress distribution within the lattice structure and a Weibull analysis was performed to gain insight into the probability of failure and to evaluate the cracking response of the structure. With this approach, the influence of lattice density on temperature, stress, and probability of failure was explored for two different lattice beam spacings. A discussion on modeling assumptions, a comparison with experimental results, and an evaluation of the lattice cracking response is provided.

14 SOLAR ENERGY↗
DOE OSTI2024

Analyzing inference workloads for spatiotemporal modeling

Ensuring power grid resiliency, forecasting climate conditions, and optimization of transportation infrastructure are some of the many application areas where data is collected in both space and time. Spatiotemporal modeling is about modeling those patterns for forecasting future trends and carrying out critical decision-making by leveraging machine learning/deep learning. Once trained offline, field deployment of trained models for near real-time inference could be challenging because performance can vary significantly depending on the environment, available compute resources and tolerance to ambiguity in results. Users deploying spatiotemporal models for solving complex problems can benefit from analytical studies considering a plethora of system adaptations to understand the associated performance-quality trade-offs. To facilitate the co-design of next-generation hardware architectures for field deployment of trained models, it is critical to characterize the workloads of these deep learning (DL) applications during inference and assess their computational patterns at different levels of the execution stack. In this paper, we develop several variants of deep learning applications that use spatiotemporal data from dynamical systems. We study the associated computational patterns for inference workloads at different levels, considering relevant models (Long short-term Memory, Convolutional Neural Network and Spatio-Temporal Graph Convolution Network), DL frameworks (Tensorflow and PyTorch), precision (FP16, FP32, AMP, INT16 and INT8), inference runtime (ONNX and AI Template), post-training quantization (TensorRT) and platforms (Nvidia DGX A100 and Sambanova SN10 RDU). Overall, our findings indicate that although there is potential in mixed-precision models and post-training quantization for spatiotemporal modeling, extracting efficiency from contemporary GPU systems might be challenging. Instead, co-designing custom accelerators by leveraging optimized High Level Synthesis frameworks (such as SODA High-Level Synthesizer for customized FPGA/ASIC targets) can make workload-specific adjustments to enhance the efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Computational Modeling, Simulation, and Potential Applications of Optical Stochastic Cooling

With the rising demand for intense particles beams, much research is being conducted in the area of particle beam cooling. One of these methods, called Stochastic Cooling (SC) (developed at CERN in the 1970's), delivered a feedback method to improve the quality and lifetime of circulating proton beams by reducing their 6D phase-space and has been widely implemented in a number of hadron machines. However, traditional stochastic cooling schemes are limited by the bandwidth of microwave frequency systems. Optical Stochastic Cooling (OSC) is a promising extension of the stochastic cooling beam cooling technique. OSC instead uses optical wavelengths which allows for improved control and increased cooling but creates its own technical challenges. This lays out work conducted toward the experimental demonstration of OSC at Fermilab's Integrable Optics Test Accelerator (IOTA) storage ring. This includes the design and characterization of parts of the optical delay system, the development and validation of a high-fidelity computational model of the OSC process, and the investigations into possible applications of the OSC mechanism to advanced beam manipulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Development of an Unmanned Mobile Current Turbine Platform: Preprint

The design and development of a prototype unmanned mobile floating platform, equipped with a custom low-flow marine current turbine for autonomously seeking and harnessing tidal/coastal currents is described. The platform is an unmanned surface vehicle in the form of a catamaran with two electric outboard motors and capabilities for autonomous navigation. An undershot water wheel, aided by a custom flow concentrator, has been selected as the basic design for the marine current turbine, which is mounted on the stern of the unmanned surface vehicle platform. The concept of operation is that the platform would navigate to a designated marine current resource, autonomously anchor at the location, align itself in the current and deploy the current turbine using a custom cable-lift deployment mechanism. As the turbine harnesses the local current, an onboard power-take-off device converts the harnessed mechanical energy to electricity which is stored in onboard batteries. Considerations of deployment in tidal and coastal waters required obtaining the necessary environmental permits for conducting in-water testing; developing required mitigation measures in protecting local wildlife and their habitats; and identifying potential in-water test sites and surveying them for their suitability in terms of the current resource, bottom type, water depth, and local boat traffic. The design and development of the turbine and the results of initial in-water testing are discussed. The potential for scaling up the system for extended capacity is presented.

ENGINEERING,TIDAL AND WAVE POWER↗
DOE OSTI2021

Exploration ToolKit (ExTK)

The Exploration Toolkit (ExTK) is a reusable Extended Reality (XR) system developed for incorporating and exploring 3-Dimensional (3-D) computer aided design (CAD) models in XR, with a primary focus on Augmented Reality (AR). The ExTK consists of a Developer Mode and a User Mode. In Developer Mode, ExTK provides developers with the ability to easily import 3-D CAD models and activate desired exploration functionality and layout. Multiple models can be added to a single instantiation of the ExTK using Unity's Scene capability. Exploration functions include scaling, rotating, explode/contract, animations, hiding parts, submodules, and measurement functions. In User Mode, ExTK provides a menu system that allows users to select models and initiate exploration functions. ExTK is architected for reusability and developers can customize the ExTK layout and functions according to application needs. ExTK is designed to be hardware agnostic, although initial development focused on the Microsoft HoloLens as the primary deployment platform. The ExTK is developed using the Unity Game Engine Development Platform. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-1506 O

Klein, BrandonThorin↗
DOE OSTI2024

MODAQ 2.0 (Modular Ocean Data AcQuisition System v2.0) [SWR-24-138]

MODAQ 2.0 is the continuation of the Modular Ocean/Offshore Data AcQuisition system which offers new benefits of being based on generally available hardware and open-source development tools. The MODAQ 2.0 project is the amalgamation of ROS packages, tools, and guidance for developing data acquisition and control applications for marine energy devices using a ROS2 based architecture. The MODAQ 2.0 Reference Design includes several packages that are targeted for the marine energy sector and allow most developers with a basic background in programming to spin up a high-quality DAQ and control system. It is comprised of multiple software repositories which can be found at the MODAQ2 GitHub Organization: https://github.com/NREL-MODAQ2

Nichols, Casey↗
DOE OSTI2021

EXAGRAPH: Graph and combinatorial methods for enabling exascale applications

Combinatorial algorithms in general and graph algorithms in particular play a critical enabling role in numerous scientific applications. However, the irregular memory access nature of these algorithms makes them one of the hardest algorithmic kernels to implement on parallel systems. With tens of billions of hardware threads and deep memory hierarchies, the exascale computing systems in particular pose extreme challenges in scaling graph algorithms. The codesign center on combinatorial algorithms, ExaGraph, was established to design and develop methods and techniques for efficient implementation of key combinatorial (graph) algorithms chosen from a diverse set of exascale applications. Algebraic and combinatorial methods have a complementary role in the advancement of computational science and engineering, including playing an enabling role on each other. In this paper, we survey the algorithmic and software development activities performed under the auspices of ExaGraph from both a combinatorial and an algebraic perspective. In particular, we detail our recent efforts in porting the algorithms to manycore accelerator (GPU) architectures. We also provide a brief survey of the applications that have benefited from the scalable implementations of different combinatorial algorithms to enable scientific discovery at scale. We believe that several applications will benefit from the algorithmic and software tools developed by the ExaGraph team.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fundamental Research Aimed at Diverting Excess Reducing Power in Photosynthesis to Orthogonal Metabolic Pathways

Photosystems are incredible biological machines that use sunlight to drive the conversion of carbon dioxide to sugar. The amount of sunlight available for photosynthesis sometimes exceeds the amount of energy plants can use. This excess energy has to be safely dissipated through non-productive biological processes. The ultimate goal of this project is to understand whether we can utilize that otherwise unused excess energy. In our previous work, we showed that, in principle, it is possible to attach a catalyst to photosystem I and generate H2 using light. Our current strategy is to genetically fuse parts of the photosystem I complex with a recently discovered oxygen-tolerant [FeFe] hydrogenase. Our rationale is that such chimeric proteins may potentially result in the natural incorporation of the photosystem I-hydrogenase link using the inherent genetic machinery of the cell. In this project, we aim to verify that light-driven hydrogen production in this construction is possible. Throughout the project, we designed nanoconstructs that showcase the plausibility of this technology, at least in vitro. We take advantage of these constructs to investigate details of the coupling between photosystem I and a H2-producing enzyme called [FeFe] hydrogenase. This part of the project reveals details of the electron transfer between photosystem I and the attached hydrogenase, providing information that can lead to new strategies for improved biological photocatalysis. We also researched efficient and robust tethering of the [FeFe] hydrogenase to photosystem I in cyanobacteria. This work will highlight successful design strategies to guide the future development of photosynthetic biohybrids. Uncovering the principles governing the utilization of otherwise unusable energy significantly further our understanding of cyanobacterial photosynthesis. The work proposed establishes the feasibility of diverting excess energy under high light conditions to orthogonal enzymatic pathways and set design rules for efficient utilization of such a strategy for scientific and industrial applications in biosensing, renewable energy, and high-value chemicals production. The work addresses the DOE-BES Photosynthetic Systems program goal to develop a multidimensional understanding of photosystems that would provide specific metrics that instruct strategies for improving biological photosynthesis and for guiding the future development of bioreactors and biomimetic energy systems.

Photosynthetic systems, hydrogenase, cyanobacteria↗
DOE OSTI2020

The emergent field of high entropy oxides: Design, prospects, challenges, and opportunities for tailoring material properties

A new class of ceramics, called entropy stabilized oxides, High Entropy Oxides (HEOs), multicomponent oxides, compositionally complex oxides, or polycation oxides, has generated considerable research interest since the first report in 2015. This multicomponent approach has created new opportunities for materials design and discovery. This Perspective will highlight some current research developments and possible applications while also providing an overview of the many successfully synthesized HEO systems to date. The polycation approach to composition development will be discussed along with a few case studies, challenges, and future possibilities afforded by this novel class of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2026

Flashlamp drive system for a high-energy 10–100 kHz pulse-burst Nd:YAG laser

We present the design of the flashlamp drive system for the NG100 laser, a 1 J/pulse, 10–100 kHz, pulse-burst Nd:YAG laser system being developed for application in a Thomson scattering plasma diagnostic. This flashlamp drive system is under active development, with a prototype now being constructed. The flashlamp drive system is modular, with each module capable of driving a series pair of linear flashlamps. Each drive module contains and is controlled by a dedicated Analog Regulator Controller. Thus each module is independently operable and controllable. This modular approach imposes no intrinsic limit to the number of modules that may be applied to drive the flashlamp pairs in a laser system. Each flashlamp drive module has a switch-regulated topology. An 1800 V main capacitor bank provides 25 kJ of energy storage, while a lower voltage output capacitor bank provides filtering and the initial energy delivered to the flashlamps at the start of the drive pulse. The main bank is recharged after each flashlamp drive pulse. As energy is drawn from the output bank by the flashlamps, an IGBT switching regulator feeds current from the main bank through an inductor to replenish the output capacitor bank. The rate of replenishment is feedback-controlled to maintain a regulated supply of power to the flashlamp load, with a setpoint range of 0.07 to 1.65 MW. An Analog Regulator Controller produces two-state variable pulse width feedback switching of the regulator IGBT. The switching frequency is ≤ 20 kHz, dynamically adjusted to limit ripple of the flashlamp power to ±3% statistical standard deviation of mean. For development or troubleshooting, each module is operable independent of the laser digital control system (microcontroller and FPGA).

Plasma diagnostics - interferometry↗
DOE OSTI2025

Visualizing energy transfer between redox-active colloids

Redox-active colloids (RACs) represent a novel class of energy carriers that exchange electrical energy upon contact. Understanding contact-mediated electron transfer dynamics in RACs offers insights into physical contact events in colloidal suspensions and enables quantification of electrical energy transport in nonconjugated polymers. Redox-based electron transport was directly observed in monolayers of micron-sized RACs containing ethyl-viologen side groups via fluorescence microscopy through an unexpected nonlinear electrofluorochromism that is quantitatively coupled to the redox state of the colloid. Via imaging studies, using this electrofluorochromism, the apparent charge transfer diffusion coefficient D CT of the RAC was easily determined. The visualization of energy transport within suspensions of redox-active colloids was also demonstrated. Our work elucidates fundamental mechanisms of energy transport in colloidal systems, informs the development of next-generation redox flow batteries, and may inspire new designs of smart active soft matter including conductive polymers for applications ranging from electrochemical sensors and organic electronics to colloidal robotics.

Qu, Alan Subing [Univ. of Illinois at Urbana-Champ↗
DOE OSTI2026

Design of Controller Hardware-In-the-Loop Model of Microgrid with Modular Building Blocks and Automated Design Script

The scalability of controller hardware-in-the-loop (CHIL) simulation is critical for validating control coordination and energy management in microgrids with distributed energy resources, especially as these modern systems become more complex and decentralized. This paper presents a CHIL modeling methodology that combines modular building blocks with an automated design script to streamline the development of high-fidelity microgrid models. Standardized subsystem templates for resources, converters, and buses are integrated with a Python-based script that compiles structured JSON configuration files into simulation-ready initialization code. The proposed approach reduces development time, improves model consistency, and enhances simulation fidelity. The methodology is validated on a Typhoon HIL604 platform and is broadly applicable to real-time simulation of complex, networked microgrid systems. This framework establishes a foundation for automated, scalable CHIL validation and accelerates the design of next-generation distributed energy systems.

Kim, Namwon [ORNL] (ORCID:0000000200438489)↗
DOE OSTI2022

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗