Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Well-posedness of a model for structural acoustic coupling in a cavity enclosed by a thin cylindrical shell

A fully coupled mathematical model describing the interactions between a vibrating thin cylindrical shell and enclosed acoustic field is presented. Because the model will ultimately be used in control applications involving piezoceramic actuators, the loads and material contributions resulting from piezoceramic patches bonded to the shell are included in the discussion. Theoretical and computational issues lead to the consideration of a weak form of the modeling set of partial differential equations (PDE's) and through the use of a semigroup formulation, well-posedness results for the system model are obtained.

Banks, H. T.↗
NASA NTRS1996

An aeroelastic analysis of helicopter rotor blades incorporating piezoelectric fiber composite twist actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consist of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for numerical integration. The twist actuation responses for three conceptual full-scale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗
NASA NTRS2001

Carbon Dioxide Clouds at High Altitude in the Tropics and in an Early Dense Martian Atmosphere

We use a time dependent, microphysical cloud model to study the formation of carbon dioxide clouds in the Martian atmosphere. Laboratory studies by Glandor et al. show that high critical supersaturations are required for cloud particle nucleation and that surface kinetic growth is not limited. These conditions, which are similar to those for cirrus clouds on Earth, lead to the formation of carbon dioxide ice particles with radii greater than 500 micrometers and concentrations of less than 0.1 cm(exp -3) for typical atmospheric conditions. Within the current Martian atmosphere, CO2 cloud formation is possible at the poles during winter and at high altitudes in the tropics during periods of increased atmospheric dust loading. In both cases, temperature perturbations of several degrees below the CO2 saturation temperature are required to nucleate new cloud particles suggesting that dynamical processes are the most common initiators of carbon dioxide clouds rather than diabatic cooling. The microphysical cloud model, coupled to a two-stream radiative transfer model, is used to reexamine the impact of CO2 clouds on the surface temperature within a dense CO2 atmosphere. The formation of carbon dioxide clouds leads to a warmer surface than what would be expected for clear sky conditions. The amount of warming is sensitive to the presence of dust and water vapor in the atmosphere, both of which act to dampen cloud effects. The radiative warming associated with cloud formation, as well as latent heating, work to dissipate the clouds when present. Thus, clouds never last for periods much longer than several days, limiting their overall effectiveness for warming the surface. The time average cloud optical depth is approximately unity leading to a 5-10 K warming, depending on the surface pressure. However, the surface temperature does not rise about the freezing point of liquid water even for pressures as high as 5 bars, at a solar luminosity of 75% the current value.

Colaprete, Anthony↗
NASA NTRS2001

Coupling of Armchair and Zigzag Tubes to a Free Electron Metal

The effect of nanotube chirality is of prime importance in determining its electronic properties. We address the issue of how chirality affects the coupling of a nanotube to metal contacts. We model coupling of armchair and zigzag nanotubes to metal contacts, in both the side- and end-contacted geometries. In the side-contacted geometry, the coupling of armchair and metallic-zigzag nanotubes to a free electron metal are significantly different. Namely, it is possible to drive a larger current through a metallic-zigzag nanotube. The predicted difference holds good when both (a) the entire circumference and (b) only a finite sector of the nanotube makes contact to the metal electrode. It might be possible to observe the predicted difference between armchair and zigzag nanotubes using gold contacts.

Anantram, M. P.↗
NASA NTRS2003

Solar Effects on Climate and the Maunder Minimum: Minimum Certainty

The current state of our understanding of solar effects on climate is reviewed. As an example of the relevant issues, the climate during the Maunder Minimum is compared with current conditions in GCM simulations that include a full stratosphere and parameterized ozone response to solar spectral irradiance variability and trace gas changes. The GISS Global Climate/Middle Atmosphere Model coupled to a q-flux/mixed layer model is used for the simulations, which begin in 1500 and extend to the present. Experiments were made to investigate the effect of total versus spectrally-varying solar irradiance changes; spectrally-varying solar irradiance changes on the stratospheric ozone/climate response with both pre-industrial and present trace gases; and the impact on climate and stratospheric ozone of the preindustrial trace gases and aerosols by themselves. The results showed that: (1) the Maunder Minimum cooling relative to today was primarily associated with reduced anthropogenic radiative forcing, although the solar reduction added 40% to the overall cooling. There is no obvious distinguishing surface climate pattern between the two forcings. (2)The global and tropical response was greater than 1 C, in a model with a sensitivity of 1.2 C per W m-2. To reproduce recent low-end estimates would require a sensitivity 1/4 as large. (3) The global surface temperature change was similar when using the total and spectral irradiance prescriptions, although the tropical response was somewhat greater with the former, and the stratospheric response greater with the latter. (4) Most experiments produce a relative negative phase of the NAO/AO during the Maunder Minimum, with both solar and anthropogenic forcing equally capable, associated with the tropical cooling and relative poleward EP flux refraction. (5) A full stratosphere appeared to be necessary for the negative AO/NAO phase, as was the case with this model for global warming experiments, unless the cooling was very large, while the ozone response played a minor role and did not influence surface temperature significantly. (6) Stratospheric ozone was most affected by the difference between present day and preindustrial atmospheric composition and chemistry, with increases in the upper and lower stratosphere during the Maunder Minimum. While the estimated UV reduction led to ozone decreases, this was generally less important than the anthropogenic effect except in the upper middle stratosphere, as judged by two different ozone photochemistry schemes. (7) The effect of the reduced solar irradiance on stratospheric ozone and on climate was similar in Maunder Minimum and current atmospheric conditions.

Rind, David↗
NASA NTRS1996

Aeroelastic Analysis of Helicopter Rotor Blades Incorporating Anisotropic Piezoelectric Twist Actuation

A simple aeroelastic analysis of a helicopter rotor blade incorporating embedded piezoelectric fiber composite, interdigitated electrode blade twist actuators is described. The analysis consists of a linear torsion and flapwise bending model coupled with a nonlinear ONERA based unsteady aerodynamics model. A modified Galerkin procedure is performed upon the rotor blade partial differential equations of motion to develop a system of ordinary differential equations suitable for dynamics simulation using numerical integration. The twist actuation responses for three conceptual fullscale blade designs with realistic constraints on blade mass are numerically evaluated using the analysis. Numerical results indicate that useful amplitudes of nonresonant elastic twist, on the order of one to two degrees, are achievable under one-g hovering flight conditions for interdigitated electrode poling configurations. Twist actuation for the interdigitated electrode blades is also compared with the twist actuation of a conventionally poled piezoelectric fiber composite blade. Elastic twist produced using the interdigitated electrode actuators was found to be four to five times larger than that obtained with the conventionally poled actuators.

Wilkie, W. Keats↗
NASA NTRS2009

Radiative Forcing and Temperature Response to Changes in Urban Albedos and Associated CO2 Offsets

The two main forcings that can counteract to some extent the positive forcings from greenhouse gases from pre-industrial times to present-day are the aerosol and related aerosol-cloud forcings, and the radiative response to changes in surface albedo. Here, we quantify the change in radiative forcing and surface temperature that may be obtained by increasing the albedos of roofs and pavements in urban areas in temperate and tropical regions of the globe. Using the catchment land surface model (the land model coupled to the GEOS-5 Atmospheric General Circulation Model), we quantify the response of the total outgoing (outgoing shortwave+longwave) radiation to urban albedo changes. Globally, the total outgoing radiation increased by 0.5 W/square m and temperature decreased by -0.008 K for an average 0.003 increase in albedo. For the U.S. the total outgoing total radiation increased by 2.3 W/square meter, and temperature decreased by approximately 0.03 K for an average 0.01 increase in albedo. These values are for the boreal summer (Tune-July-August). Based on these forcings, the expected emitted CO2 offset for a plausible 0.25 and 0.15 increase in albedos of roofs and pavements, respectively, for all global urban areas, was found to be approximately 57 Gt CO2 . A more meaningful evaluation of the impacts of urban albedo increases on climate and the expected CO2 offsets would require simulations which better characterizes urban surfaces and represents the full annual cycle.

Menon, Surabi↗
NASA NTRS2014

Global Warming and 21st Century Drying

Global warming is expected to increase the frequency and intensity of droughts in the twenty-first century, but the relative contributions from changes in moisture supply (precipitation) versus evaporative demand (potential evapotranspiration; PET) have not been comprehensively assessed. Using output from a suite of general circulation model (GCM) simulations from phase 5 of the Coupled Model Intercomparison Project, projected twentyfirst century drying and wetting trends are investigated using two offline indices of surface moisture balance: the Palmer Drought Severity Index (PDSI) and the Standardized Precipitation Evapotranspiration Index (SPEI). PDSI and SPEI projections using precipitation and Penman- Monteith based PET changes from the GCMs generally agree, showing robust cross-model drying in western North America, Central America, the Mediterranean, southern Africa, and the Amazon and robust wetting occurring in the Northern Hemisphere high latitudes and east Africa (PDSI only). The SPEI is more sensitive to PET changes than the PDSI, especially in arid regions such as the Sahara and Middle East. Regional drying and wetting patterns largely mirror the spatially heterogeneous response of precipitation in the models, although drying in the PDSI and SPEI calculations extends beyond the regions of reduced precipitation. This expansion of drying areas is attributed to globally widespread increases in PET, caused by increases in surface net radiation and the vapor pressure deficit. Increased PET not only intensifies drying in areas where precipitation is already reduced, it also drives areas into drought that would otherwise experience little drying or even wetting from precipitation trends alone. This PET amplification effect is largest in the Northern Hemisphere mid-latitudes, and is especially pronounced in western North America, Europe, and southeast China. Compared to PDSI projections using precipitation changes only, the projections incorporating both precipitation and PET changes increase the percentage of global land area projected to experience at least moderate drying (PDSI standard deviation of < or = -1) by the end of the twenty-first century from 12 to 30 %. PET induced moderate drying is even more severe in the SPEI projections (SPEI standard deviation of > or = -1; 11 to 44 %), although this is likely less meaningful because much of the PET induced drying in the SPEI occurs in the aforementioned arid regions. Integrated accounting of both the supply and demand sides of the surface moisture balance is therefore critical for characterizing the full range of projected drought risks tied to increasing greenhouse gases and associated warming of the climate system.

moisture↗
NASA NTRSDate not supplied

Climate Change Impact on Surface Ozone Based on CMIP6 Earth System Models

It is presented an analysis of the effect of climate change on surface ozone (O3) discussing the related penalties and benefits around the globe from the global modeling perspective based on simulations with five CMIP6 (Coupled Model Intercomparison Project Phase 6) Earth System Models. All models conducted simulation experiments considering future climate (ssp370SST) and present-day climate (ssp370pdSST) under the same future emissions scenario (SSP3-7.0). Over regions remote from pollution sources, there is a robust decline in mean surface ozone concentration varying spatially from -0.2 to -2 ppbv oC-1, with strongest decline over tropical oceanic regions, which is mainly linked to the dominating role of enhanced ozone chemical loss with higher water vapour abudances under a warmer climate. However, ozone increases over regions close to anthropogenic pollution sources or close to enhanced natural Biogenic Volatile Organic Compounds (BVOC) emission sources with a rate ranging regionally from 0.2 to 2 ppbv oC-1, implying a regional surface ozone penalty due to global warming. The individual models show this robustly for south-eastern China and India as well as for regions of Africa but there are inter-model differences in areas within Europe and the United States (US) as well as in South America. The future climate change enhances the efficiency of precursor emissions to generate surface ozone in polluted regions and thus the magnitude of this effect depends on the regional emission changes considered in this study within the SSP3_7.0 scenario. The comparison of the climate change impact effect on surface ozone versus the combined effect of climate and emission changes indicates the dominant role of precursor emission changes in projecting surface ozone concentrations under future climate change scenarios.

climate change↗
NASA NTRS2024

NPCC4: Tail Risk, Climate Drivers of Extreme Heat, and New Methods for Extreme Event Projections

We summarize historic New York City (NYC) climate change trends and provide the latest scientific analyses on projected future changes based on a range of global greenhouse gas emissions scenarios. Building on previous NPCC assessment reports, we describe new methods used to develop the projections of record for sea level rise, temperature, and precipitation for NYC, across multiple emissions pathways and analyze the issue of the “hot models” associated with the 6th phase of the Coupled Model Intercomparison Project (CMIP6) and their potential impact on NYC's climate projections. We describe the state of the science on temperature variability within NYC and explain both the large-scale and regional dynamics that lead to extreme heat events, as well as the local physical drivers that lead to inequitable distributions of exposure to extreme heat. We identify three areas of tail risk and potential for its mischaracterization, including the physical processes of extreme events and the effects of a changing climate. Finally, we review opportunities for future research, with a focus on the hot model problem and the intersection of spatial resolution of projections with gaps in knowledge in the impacts of the climate signal on intraurban heat and heat exposure.

NPCC4↗
NASA NTRS2024

How Can Future Climate Change Affect the Corn Production System in Nebraska, USA?

The increasing temperatures and changing precipitation patterns have affected both crop yields and water requirements challenging agricultural systems in sustaining high crop yields and ensuring water and food security. In the present study, we evaluated the impact of climate change on the corn production under irrigated and rainfed cropping systems in Nebraska considering four climate zones until 2100 as a case study for the US Corn Belt. Additionally, suggests ways to reduce the impact of climate change on corn yield and water resources. Five global circulation model datasets from the sixth phase of the Coupled Model Intercomparison Project (CMIP6) interconnected to Shared Socioeconomic Pathways (SSPs), were selected. Three future periods: near (2020–2039), mid (2040–2069), and far future (2070–2100) under SSP226, SSP370, and SSP585 scenarios were used to predict the climate changes, crop cycle length, crop yield, and crop water requirement. Additionally, we input annual CO2 concentration scenarios for each SSP scenario to model changes from 2015 to 2100 using the CERES-Maize model. Temperatures are forecasted to rise, precipitation to decline, and crop cycle length to decrease, particularly for SSP370 (up to 49% yield decrease) and SSP585 (up to 67% yield decrease) by 2070–2100. Yields are expected to decline across all scenarios and periods for irrigated and rainfed fields due to the reduction in crop cycle length. Crop water demand is projected to increase, with rainfed systems facing heightened impact due to reduced precipitation. Addressing climate change in Nebraska requires actions like developing new corn varieties to longer the crop season, transitioning to efficient irrigation systems, improving irrigation management for water and food security, and reducing the risk of conflicts over its use. These results open possibilities for discussion about the actions that can be taken to cope with climate change and guarantee water and food security in the US Corn Belt.

CMIP6↗
NASA NTRSDate not supplied

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗
NASA NTRS1991

An equilibrium model for the coupled ocean-atmosphere boundary layer in the tropics

An atmospheric convective boundary layer (CBL) model is coupled to an ocean mixed-layer (OML) model in order to study the equilibrium state of the coupled system in the tropics, particularly in the Pacific region. The equilibrium state of the coupled system is solved as a function of sea-surface temperature (SST) for a given surface wind and as a function of surface wind for a given SST. It is noted that in both cases, the depth of the CBL and OML increases and the upwelling below the OML decreases, corresponding to either increasing SST or increasing surface wind. The coupled ocean-atmosphere model is solved iteratively as a function of surface wind for a fixed upwelling and a fixed OML depth, and it is observed that SST falls with increasing wind in both cases. Realistic gradients of mixed-layer depth and upwelling are observed in experiments with surface wind and SST prescribed as a function of longitude.

Sui, C.-H.↗
NASA NTRS2023

Modeling Sea Ice Albedo and Transmittance Measurements With A Fully-Coupled Radiative Transfer Model

A rigorous treatment of the sea ice medium has been incorporated in the advanced Coupled Ocean-Atmosphere Radiative Transfer (COART) model. The inherent optical properties (IOPs) of brine pockets and air bubbles over the 0.25-4.0 μm spectral region are parameterized as a function of the vertical profile of the sea ice physical properties (temperature, salinity and density). We test the model performance using available albedo and transmittance measurements collected during the Impacts of Climate on the Ecosystems and Chemistry of the Arctic Pacific Environment (ICESCAPE) and the Surface Heat Budget of the Arctic Ocean (SHEBA) field campaigns. The observations are adequately simulated when at least three layers are used to represent bare (first-year and multi-year) ice, including a thin top layer characterized by low density and high scattering. Two layers can be sufficient to model isolated cases of multi-year ice, and apply well to ponded ice except for shallow ponds over thick ice. The albedo and transmittance of ponded ice in the visible are mainly determined by the optical properties of the ice underlying the water layer used to model the pond. Sensitivity results indicate that the air volume or ice density has the largest impact on the simulated fluxes. Possible contamination from light-absorbing impurities, such as black carbon or ice algae, is also implemented in the model and is able to effectively reduce the albedo and transmittance in the visible spectrum to further improve the model-observation agreement. The agreement between the observed and modeled spectra validates the parameterization of the sea ice IOPs, and endorses COART as an accurate tool for radiation studies in the cryosphere.

Coupled Ocean-9 Atmosphere Radiative Transfer (COA↗
DOE OSTI2016

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. Lastly, HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

58 GEOSCIENCES↗
DOE OSTI2013

Stimulation at Desert Peak -modeling with the coupled THM code FEHM

Numerical modeling of the 2011 shear stimulation at the Desert Peak well 27-15. This submission contains the FEHM executable code for a 64-bit PC Windows-7 machine, and the input and output files for the results presented in the included paper from ARMA-213 meeting.

15 GEOTHERMAL ENERGY↗
NASA NTRS2016

Quantifying PM2.5-Meteorology Sensitivities in a Global Climate Model

Climate change can influence fine particulate matter concentrations (PM2.5) through changes in air pollution meteorology. Knowledge of the extent to which climate change can exacerbate or alleviate air pollution in the future is needed for robust climate and air pollution policy decision-making. To examine the influence of climate on PM2.5, we use the Geophysical Fluid Dynamics Laboratory Coupled Model version 3 (GFDL CM3), a fully-coupled chemistry-climate model, combined with future emissions and concentrations provided by the four Representative Concentration Pathways (RCPs). For each of the RCPs, we conduct future simulations in which emissions of aerosols and their precursors are held at 2005 levels while other climate forcing agents evolve in time, such that only climate (and thus meteorology) can influence PM2.5 surface concentrations. We find a small increase in global, annual mean PM2.5 of about 0.21 micro-g/cu m3 (5%) for RCP8.5, a scenario with maximum warming. Changes in global mean PM2.5 are at a maximum in the fall and are mainly controlled by sulfate followed by organic aerosol with minimal influence of black carbon. RCP2.6 is the only scenario that projects a decrease in global PM2.5 with future climate changes, albeit only by -0.06 micro-g/cu m (1.5%) by the end of the 21st century. Regional and local changes in PM2.5 are larger, reaching upwards of 2 micro-g/cu m for polluted (eastern China) and dusty (western Africa) locations on an annually averaged basis in RCP8.5. Using multiple linear regression, we find that future PM2.5 concentrations are most sensitive to local temperature, followed by surface wind and precipitation. PM2.5 concentrations are robustly positively associated with temperature, while negatively related with precipitation and wind speed. Present-day (2006-2015) modeled sensitivities of PM2.5 to meteorological variables are evaluated against observations and found to agree reasonably well with observed sensitivities (within 10e50% over the eastern United States for several variables), although the modeled PM2.5 is less sensitive to precipitation than in the observations due to weaker convective scavenging. We conclude that the hypothesized "climate penalty" of future increases in PM2.5 is relatively minor on a global scale compared to the influence of emissions on PM2.5 concentrations.

PM2.5↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗