Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗
DOE OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
DOE OSTI2025

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
DOE OSTI2022

Electron temperature gradient driven transport model for tokamak plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the multi-mode anomalous transport module [Rafiq et al., Phys Plasmas 20, 032506 (2013)] to predict a time-dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, β and gradient of β (β'), collisionality, safety factor, and toroidicity, where β is the ratio of the plasma pressure to the magnetic pressure. Furthermore, the electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Machine Learning Automated Analysis of Enormous Synchrotron X-ray Diffraction Datasets

X-ray diffraction (XRD) data analysis can be a time-consuming and laborious task. Deep neural network (DNN) based models trained with synthetic XRD patterns have been proven to be a highly efficient, accurate, and automated method for analyzing common XRD data collected from solid samples in ambient environments. However, it remains unclear whether synthetic XRD-based models can be effective in solving micro(μ)-XRD mapping data for in situ experiments involving liquid phases, which always have lower quality and significant artifacts. In this study, we collected μ-XRD mapping data from a LaCl 3 -calcite hydrothermal fluid system and trained two categories of models to analyze the experimental XRD patterns. Here, the models trained solely with synthetic XRD patterns showed low accuracy (as low as 64%) when solving experimental μ-XRD mapping data. However, the accuracy of the DNN models significantly improved (90% or above) when we trained them with a data set containing both synthetic and a small number of labeled experimental μ-XRD patterns. This study highlights the importance of labeled experimental patterns in training DNN models to solve μ-XRD mapping data from in situ experiments involving liquid phases.

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

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Enhanced Models for the Simulation of Two-Phase Boiling Flow and Critical Heat Flux with the STAR-CCM+ and NEK-2P CFD Codes: Executive Summary (Final Report)

The goal of the project was to develop, verify and validate an advanced two-phase boiling flow modeling capability that enhances the existing two-phase flow modeling of the STAR-CCM+ Computational Fluid Dynamics (CFD) code. The project focused on the simulation of two phase flow and heat transfer phenomena that occur in Boiling Water Reactor (BWR) and Pressurized Water Reactor (PWR) fuel bundles under various operating conditions. The extended STAR-CCM+ two-phase flow models implemented during the project follow the approach similar to the Extended Boiling Framework [1-6] previously developed at Argonne but include more fundamental physical models of boiling phenomena. The development and validation of the extended two-phase models is supported by the U.S. DOE GAIN program and performed in collaboration with Framatome U.S. Inc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Lagrangian and Eulerian drag models that are consistent between Euler-Lagrange and Euler-Euler (two-fluid) approaches for homogeneous systems

The undisturbed flow of a particle is of fundamental importance since it controls both the undisturbed flow force and the perturbation force (which includes quasisteady, added-mass, and history forces). Here we use the pairwise interaction extended point particle framework to evaluate the undisturbed flow of each particle through superposition of the perturbation flow induced by all its neighbors. This approach allows calculation of various statistics related to undisturbed fluid velocity under conditions of both stationary and nonstationary particles. In a random distribution of stationary particles, while the macroscale undisturbed flow is slowly varying, the microscale undisturbed flow that arises due to the perturbation flow of neighbors varies substantially from one particle to another and this in turn leads to large variation in the hydrodynamic force exerted on the particles. The effect of particle motion is generally to increase the particle-to-particle variation in the undisturbed fluid velocity of the particles. We observe that this increase is greater for the transverse component than for the streamwise component. As a result, with increasing random particle motion, the distribution of undisturbed fluid-velocity fluctuation becomes isotropic. Three different normalized forces are defined: $Φ^L$ is the Lagrangian normalized force on an individual particle suitable for application in a microscale-informed Euler-Lagrange simulation, $Φ^E$ is the Eulerian normalized average force suitable for application in an Euler-Euler simulation, and $Φ^{LE}$ is the Lagrangian normalized force on an individual particle suitable for application in the standard Euler-Lagrange simulation. In this work, we establish precise relations between these different definitions. The drag laws developed based on particle-resolved direct numerical simulation results and experiments are appropriate for application only as the Eulerian normalized average force. We introduce the force consistency relation and use it to obtain an expression for $Φ^L$, which when applied to each particle and averaged over all the particles equals $Φ^E$. The results are first obtained in the limit of stationary particles and then extended to the general case of nonstationary particles.

42 ENGINEERING↗
DOE OSTI2021

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model

Geologic carbon storage is one of the promising strategies to mitigate climate change by reducing the emission of carbon dioxide to the atmosphere. As part of the National Risk Assessment Partnership (NRAP), a systems-level stochastic analysis tool called the open source integrated assessment model, NRAP-Open-IAM, has been developed to estimate and manage the risk of containment loss at a geological carbon sequestration site. NRAP-Open-IAM contains several wellbore leakage model components that estimate the fluid leak rate that may occur through compromised legacy wells due to the increase in pressure resulting from CO 2 injection activities. Coupled to a reservoir component model, these components estimate the leakage of CO 2 and/or brine from a storage reservoir to overlying aquifer layers and the atmosphere through legacy wells. This report presents the theoretical framework and quality testing of the multisegmented wellbore reduced-order model. The model allows for segmenting of the legacy wells passing through the overlying stratigraphy into several intervals to simulate a site’s specific stratigraphic and hydrogeologic properties. For quality assurance, the analytical model is validated against numerical reservoir flow simulations for single and multiple aquifer(s) models. The results indicate that the model accurately predicts the transport of two-phase fluids (brine and injected CO 2 ) through the well over time. A detailed description of the model helps users to understand the model and provides a basis for future improvements.

58 GEOSCIENCES↗