Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing the CP structure of the top quark Yukawa at the future muon collider

We study the top-Higgs coupling with a CP violating phase ξ at a future multi-TeV muon collider. We focus on processes that are directly sensitive to the top quark Yukawa coupling: $t$$\overline{t}$$h$, $tbhμν$, and $t$$\overline{t}$$h$$v$$\overline{v}$ with $h$ → $b$$\overline{b}$ and semileptonic top decays. At different energies, different processes dominate the cross section, providing complementary information. At and above an energy of $\mathcal{O}$(10) TeV, vector boson fusion processes dominate. As we show, in the Standard Model there is destructive interference in the vector boson fusion processes $t$$\overline{t}$$h$$v$$\overline{v}$ and $tbhμν$ between the top quark Yukawa and Higgs-gauge boson couplings. A CP-violating phase changes this interference, and the cross section measurement is very sensitive to the size of the CP-violating angle. Although we find that the cross sections are measured to $\mathcal{O}$(50%) statistical uncertainty at 1σ, a 10 and 30 TeV muon collider can bound the CP-violating angle |ξ| ≲ 9.0° and |ξ| ≲ 5.4°, respectively. However, cross section measurements are insensitive to the sign of the CP-violating angle. To determine that the coupling is truly CP violating, observables sensitive to CP-violation must be measured. We find in the $t$$\overline{t}$$h$ process the azimuthal angle between the $t$ + $\overline{t}$ plane and the initial state muon+Higgs plane shows good discrimination for ξ = ±0.1π.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1996

A Sundial-Atlas Precursor to the TIMED Mission: A Quick-Response Global Investigation into Coupled Lower Thermospheric, Ionospheric, and Mesospheric Physics

The SUNDIAL-ATLAS effort was a global-scale investigation which responded to the science priorities of the ITM Panel, the Integrated SPD Strategy Implementation Plan as a whole, and the need for potential cost-saving design criteria for the TIMED mission. The investigation focused on coupling processes in the ionospheric-thermospheric system, taking advantage of the timelines of the ATLAS-1 mission (March 1992), and the global-scale ground-based measurement and modeling activities of the SUNDIAL program. The collaborative SUNDIAL-ATLAS activity was the first opportunity for global measurements of the chemistry, kinetics, and electrodynamics which couple the E-, Fl-, and F2-regions into a single interactive system. As such, the program represented an important first step in studying global issues; and accordingly, was an important proof of concept experiment relevant to the strategic mission plans for the ITM community and the upcoming intermediate class satellite program called TIMED. To meet its projected goals, TIMED must perform a number of critical measurements and execute a number of correlations that were to be tried and tested for the first time in the SUNDIAL-ATLAS investigation. This was designed to include global correlations of thermospheric and ionospheric composition during quiet and disturbed conditions and the co-registration of global-scale ground-based measurements with along-track satellite diagnostics. The SUNDIAL component of the current investigation addressed this need by acquiring, reducing, and analyzing a multi-sensor database that complemented and extended that which was generated in the ATLAS mission (Atmospheric Laboratory for Applications and Science). The SUNDIAL data defined the state and condition of the global-scale ionosphere in the altitude range from 100 km to the F2-peak. These data specified the peak heights and densities of the E-, Fl-, and F2-regions, along with the global distributions of intermediate, descending, and sequential layers which play a critical role in the dynamo region of the lower ionospheric-thermospheric domain. The data were collected by the SUNDIAL network of more than 50 ground-based stations utilizing ionosondes, radars, photometers, Fabry-Perot interferometers, and total electron content measurements. The data were acquired during a three-week period centered on the eight-day ATLAS-1 mission, which provided image and photometric sensing of the altitude distributions of the major and minor ions and neutrals in the ITM system. This report focuses on the scientific contributions of the SUNDIAL component of the overall investigation. Specific findings are described in seven papers (attached) published in the Journal of Geophysical Research.

Szuszczewicz, E. P.↗
DOE OSTI2019

Coupling constant corrections in a holographic model of heavy ion collisions with nonzero baryon number density

Sufficiently energetic collisions of heavy ions result in the formation of a droplet of a strongly coupled liquid state of QCD matter known as quark-gluon plasma. By using gauge-gravity duality (holography), a model of a rapidly hydrodynamizing and thermal- izing process like this can be constructed by colliding sheets of energy density moving at the speed of light and tracking the subsequent evolution. In this work, we consider the dual gravitational description of such collisions in the most general bulk theory with a four-derivative gravitational action containing a dynamical metric and a gauge field in five dimensions. Introducing the bulk gauge field enables the analysis of collisions of sheets which carry nonzero “baryon” number density in addition to energy density. Introducing the four-derivative terms enables consideration of such collisions in a gauge theory with finite gauge coupling, working perturbatively in the inverse coupling. While the dynamics of energy and momentum in the presence of perturbative inverse-coupling corrections has been analyzed previously, here we are able to determine the effect of such finite coupling corrections on the dynamics of the density of a conserved global charge, which we take as a model for the dynamics of nonzero baryon number density. In accordance with expectations, as the coupling is reduced we observe that after the collisions less baryon density ends up stopped at mid-rapidity and more of it ends up moving near the lightcone.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Investigation of a model based reconstruction technique for neutron tomography at the PULSTAR reactor

Herein, a model based neutron computed tomography (CT) reconstruction technique is developed and implemented to improve the spatial resolution of the reconstructed images. The spatial resolution of a scintillator and a charge-coupled device (CCD) based digital neutron imaging system is quantified and used as the system response function in the iterative reconstruction process. The response function is included in the forward projection model to provide better estimation of the actual projection, thus improving the spatial resolution of the reconstructed images. Experiments were performed at the neutron imaging facility of the PULSTAR research reactor to corroborate the performance of the developed algorithm with a 50 μm and 250 μa m 6 LiF:ZnS thick scintillator. The reconstruction results of the model based technique have shown an improvement of nearly 40% in achievable spatial resolution as compared to the standard filtered back projection (FBP) technique.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Predicting synthetic mRNA stability using massively parallel kinetic measurements, biophysical modeling, and machine learning

Abstract mRNA degradation is a central process that affects all gene expression levels, though it remains challenging to predict the stability of a mRNA from its sequence, due to the many coupled interactions that control degradation rate. Here, we carried out massively parallel kinetic decay measurements on over 50,000 bacterial mRNAs, using a learn-by-design approach to develop and validate a predictive sequence-to-function model of mRNA stability. mRNAs were designed to systematically vary translation rates, secondary structures, sequence compositions, G-quadruplexes, i-motifs, and RppH activity, resulting in mRNA half-lives from about 20 seconds to 20 minutes. We combined biophysical models and machine learning to develop steady-state and kinetic decay models of mRNA stability with high accuracy and generalizability, utilizing transcription rate models to identify mRNA isoforms and translation rate models to calculate ribosome protection. Overall, the developed model quantifies the key interactions that collectively control mRNA stability in bacterial operons and predicts how changing mRNA sequence alters mRNA stability, which is important when studying and engineering bacterial genetic systems.

Cetnar, Daniel P.↗
NASA NTRS2004

Magnetohydrodynamic Convection in the Outer Core and its Geodynamic Consequences

The Earth's fluid outer core is in vigorous convection through much of the Earth's history. In addition to generating and maintaining Earth s time-varying magnetic field (geodynamo), the core convection also generates mass redistribution in the core and a dynamical pressure field on the core-mantle boundary (CMB). All these shall result in various core-mantle interactions, and contribute to surface geodynamic observables. For example, electromagnetic core-mantle coupling arises from finite electrically conducting lower mantle; gravitational interaction occurs between the cores and the heterogeneous mantle; mechanical coupling may also occur when the CMB topography is aspherical. Besides changing the mantle rotation via the coupling torques, the mass-redistribution in the core shall produce a spatial-temporal gravity anomaly. Numerical modeling of the core dynamical processes contributes in several geophysical disciplines. It helps explain the physical causes of surface geodynamic observables via space geodetic techniques and other means, e.g. Earth's rotation variation on decadal time scales, and secular time-variable gravity. Conversely, identification of the sources of the observables can provide additional insights on the dynamics of the fluid core, leading to better constraints on the physics in the numerical modeling. In the past few years, our core dynamics modeling efforts, with respect to our MoSST model, have made significant progress in understanding individual geophysical consequences. However, integrated studies are desirable, not only because of more mature numerical core dynamics models, but also because of inter-correlation among the geophysical phenomena, e.g. mass redistribution in the outer core produces not only time-variable gravity, but also gravitational core-mantle coupling and thus the Earth's rotation variation. They are expected to further facilitate multidisciplinary studies of core dynamics and interactions of the core with other components of the Earth.

Kuang, Weijia↗
DOE OSTI2011

Polarization of nuclear spins by a cold nanoscale resonator

A cold nanoscale resonator coupled to a system of nuclear spins can induce spin relaxation. In the low-temperature limit where spin-lattice interactions are ''frozen out,'' spontaneous emission by nuclear spins into a resonant mechanical mode can become the dominant mechanism for cooling the spins to thermal equilibrium with their environment. We provide a theoretical framework for the study of resonator-induced cooling of nuclear spins in this low-temperature regime. Relaxation equations are derived from first principles, in the limit where energy donated by the spins to the resonator is quickly dissipated into the cold bath that damps it. A physical interpretation of the processes contributing to spin polarization is given. For a system of spins that have identical couplings to the resonator, the interaction Hamiltonian conserves spin angular momentum, and the resonator cannot relax the spins to thermal equilibrium unless this symmetry is broken by the spin Hamiltonian. The mechanism by which such a spin system becomes ''trapped'' away from thermal equilibrium can be visualized using a semiclassical model, which shows how an indirect spin-spin interaction arises from the coupling of multiple spins to one resonator. The internal spin Hamiltonian can affect the polarization process in two ways: (1) By modifying the structure of the spin-spin correlations in the energy eigenstates, and (2) by splitting the degeneracy within a manifold of energy eigenstates, so that zero-frequency off-diagonal terms in the density matrix are converted to oscillating coherences. Shifting the frequencies of these coherences sufficiently far from zero suppresses the development of resonator-induced correlations within the manifold during polarization from a totally disordered state. Modification of the spin-spin correlations by means of either mechanism affects the strength of the fluctuating spin dipole that drives the resonator. In the case where product states can be chosen as energy eigenstates, spontaneous emission from eigenstate populations into the resonant mode can be interpreted as independent emission by individual spins, and the spins relax exponentially to thermal equilibrium if the development of resonator-induced correlations is suppressed. When the spin Hamiltonian includes a significant contribution from the homonuclear dipolar coupling, the energy eigenstates entail a correlation specific to the coupling network. Simulations of dipole-dipole coupled systems of up to five spins suggest that these systems contain weakly emitting eigenstates that can trap a fraction of the population for time periods >>100/R{sub 0}, where R{sub 0} is the rate constant for resonator-enhanced spontaneous emission by a single spin 1/2. Much of the polarization, however, relaxes with rates comparable to R{sub 0}. A distribution of characteristic high-field chemical shifts tends to increase the relaxation rates of weakly emitting states, enabling transitions to states that can quickly relax to thermal equilibrium. The theoretical framework presented in this paper is illustrated with discussions of spin polarization in the contexts of force-detected nuclear-magnetic-resonance spectroscopy and magnetic-resonance force microscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Data-Driven Nonlinear Subspace Modeling for Prediction and Control of Molten Iron Quality Indices in Blast Furnace Ironmaking

Blast furnace (BF) in ironmaking is a nonlinear dynamic process with complicated physical-chemical reactions, where multi-phase and multi-field coupling and large time delay occur during its operation. In BF operation, the molten iron temperature (MIT) as well as Si, P and S contents of molten iron are the most essential molten iron quality (MIQ) indices, whose measurement, modeling and control have always been important issues in metallurgic engineering and automation field. This paper develops a novel data-driven nonlinear state space modeling for the prediction and control of multivariate MIQ indices by integrating hybrid modeling and control techniques. First, to improve modeling efficiency, a data-driven hybrid method combining canonical correlation analysis and correlation analysis is proposed to identify the most influential controllable variables as the modeling inputs from multitudinous factors would affect the MIQ indices. Then, a Hammerstein model for the prediction of MIQ indices is established using the LS-SVM based nonlinear subspace identification method. Such a model is further simplified by using piecewise cubic Hermite interpolating polynomial method to fit the complex nonlinear kernel function. Compared to the original Hammerstein model, this simplified model can not only significantly reduce the computational complexity, but also has almost the same reliability and accuracy for a stable prediction of MIQ indices. Last, in order to verify the practicability of the developed model, it is applied in designing a genetic algorithm based nonlinear predictive controller for multivariate MIQ indices by directly taking the established model as a predictor. Industrial experiments show the advantages and effectiveness of the proposed approach.

Blast furnace (BF)↗
DOE OSTI2025

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗
DOE OSTI2025

Photoelectrochemical Hydride Generation with Oxide-Coated Silicon

Photoelectrochemical generation of a potent organic hydride donor at silicon is demonstrated. Two different oxide-coated p-type silicon photoelectrodes reduced 1,2,3,5,6-pentamethyl-1H-benzo[d]imidazol-3-ium hexafluorophosphate, [PMBI][PF 6 ], to its corresponding imidazole hydride, PMBIH, in the presence of a proton source. Under 1 sun illumination, in acetonitrile with 1,5-diazabicyclo[4.3.0]non-5-ene (DBN) buffer, the p-Si photoelectrodes convert PMBI + to PMBIH with good Faradaic efficiencies (FEs): 78% FE at −2.3 V vs Fc +/0 for Si|TiO 2 and 83% FE at −2.6 V vs Fc +/0 for Si|SiO 2 (where Si|SiO 2 represents silicon coated with an oxide layer). Generally, the Si|TiO 2 catalyzed the reaction at milder potentials than Si|SiO 2 , but the Si|SiO 2 had better selectivity for PMBIH generation over H 2 evolution than Si|TiO 2 . In light of prior studies of these photoelectrodes, the differences in selectivity and onset potential suggest a difference in mechanism, likely an outer-sphere electron transfer (ET) mechanism at the SiO 2 surface and potentially a proton-coupled ET process at the TiO 2 surface. To help understand reaction efficiency and identify areas of improvement, a thermochemical model for understanding net hydride transfer from the semiconductor to an acceptor in solution was developed. We find that the reactions in the present system are quite downhill. The high overpotentials (even while maintaining selectivity over H 2 evolution) emphasize the need for improved catalysts. Furthermore, this approach to evaluate the thermodynamics of net hydride transfer should be broadly valuable for electrochemical and photoelectrochemical processes.

14 SOLAR ENERGY↗
DOE OSTI2024

Causally‐Informed Deep Learning to Improve Climate Models and Projections

Abstract Climate models are essential to understand and project climate change, yet long‐standing biases and uncertainties in their projections remain. This is largely associated with the representation of subgrid‐scale processes, particularly clouds and convection. Deep learning can learn these subgrid‐scale processes from computationally expensive storm‐resolving models while retaining many features at a fraction of computational cost. Yet, climate simulations with embedded neural network parameterizations are still challenging and highly depend on the deep learning solution. This is likely associated with spurious non‐physical correlations learned by the neural networks due to the complexity of the physical dynamical system. Here, we show that the combination of causality with deep learning helps removing spurious correlations and optimizing the neural network algorithm. To resolve this, we apply a causal discovery method to unveil causal drivers in the set of input predictors of atmospheric subgrid‐scale processes of a superparameterized climate model in which deep convection is explicitly resolved. The resulting causally‐informed neural networks are coupled to the climate model, hence, replacing the superparameterization and radiation scheme. We show that the climate simulations with causally‐informed neural network parameterizations retain many convection‐related properties and accurately generate the climate of the original high‐resolution climate model, while retaining similar generalization capabilities to unseen climates compared to the non‐causal approach. The combination of causal discovery and deep learning is a new and promising approach that leads to stable and more trustworthy climate simulations and paves the way toward more physically‐based causal deep learning approaches also in other scientific disciplines.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Growth rate as a link between microbial diversity and soil biogeochemistry

The growth rate of a microorganism is a simple yet profound way to quantify its impact on the world. The absolute growth rate of a microbial population reflects rates of resource assimilation, biomass production, and element transformation, some of the many ways that organisms affect Earth’s ecosystems and climate. Microbial fitness in the environment depends on the ability to reproduce quickly when conditions are favorable and adopt a survival physiology when conditions worsen, which cells coordinate by adjusting their relative growth rate. At the population level, relative growth rate is a sensitive metric of fitness, linking survival and reproduction to the ecology and evolution of populations. Techniques combining ‘omics and stable isotope probing enable sensitive measurements of growth rates of microbial assemblages and individual taxa in soil. Microbial ecologists can explore how the growth rates of taxa with known traits and evolutionary histories respond to changes in resource availability, environmental conditions, and interactions with other organisms. We anticipate that quantitative and scalable data on the growth rates of soil microorganisms, coupled with measurements of biogeochemical fluxes, will allow scientists to test and refine ecological theory and advance process-based models of carbon flux, nutrient uptake, and ecosystem productivity. Finally, measurements of in situ microbial growth rates provide insights into the ecology of populations and can be used to quantitatively link microbial diversity to soil biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effective field theory for radiative corrections to charged-current processes. II. Axial-vector coupling

Here we discuss the hadronic structure-dependent radiative corrections to the axial-vector coupling that controls single-nucleon weak charged-current processes—commonly denoted by 𝑔 𝐴 . We match the Standard Model at the GeV scale onto chiral perturbation theory at next-to-leading order in the one-nucleon sector, in the presence of electromagnetic and weak interactions. As a result, we provide a representation for the corrections to 𝑔 𝐴 in terms of infrared finite convolutions of simple kernels with the single-nucleon matrix elements of time-ordered products of two and three quark bilinears (vector, axial-vector, and pseudoscalar). We discuss strategies to determine the required nonperturbative input from data, lattice-QCD (+QED), and possibly hadronic models. This work paves the way for a precise comparison of the values of the ratio 𝑔 𝐴 ⁡/𝑔 𝑉 extracted from experiment and from lattice QCD, which constrain physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Competing ionization and dissociation in the H 2 gerade system

A numerically solvable two-dimensional (2D) model, employed by the authors to study the dissociative recombination of H$^+_2$ in the ungerade symmetry, is extended to describe the collision process in the gerade symmetry of H$_2$. In this symmetry, the ionization and dissociation processes are driven primarily by the direct, curve-crossing mechanism. The model is represented by a set of three coupled electronic channels in 2D, in the space of $\textit{s, p, d}$ partial waves of the colliding electron. In this work, we demonstrate that the Born-Oppenheimer properties of the H$_2$ molecule in the relevant range of internuclear distances can be described by such a model. The molecular rotational degrees of freedom are accounted for by the rotational frame transformation. The numerical solution of the model is discussed, and the resulting rovibrationally inelastic and dissociative recombination cross sections are compared with the available data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Achieving Cyber-Resilience for Power Systems using a Learning, Model-Assisted Blockchain Framework

The secure integration and management of distributed energy resources (DER) and power aggregators in the electric grid requires secure communications and a physics-aware Command and Control (C2) strategy. A Blockchain (BC)-based overlay network was developed to provide a security layer for the existing power grid network that mitigates risks in current and legacy network and C2 protocols. By integrating a Model-Assisted Machine Learning (MAML) framework with a Secure Blockchain Overlay Network (SBON) a defense-in-depth strategy was achieved. In our approach, the MAML framework leveraged a smart contract framework to gather network data and learn the dynamics of DER to develop detection strategies for attacks targeting sensors and actuators used by DER. The MAML framework learned dynamical systems models for individual DERs to detect sensor attacks. For DER we utilized a Digital Twin (DT) to accelerate the learning process for a model resistant to stealthy attacks. The project created DT for PV inverters and BESS. The DTs were coupled with a model-assisted, data-driven learning of DER behavior. Specifically, we evaluated architectures for model-based learning with model-free fine-tuning. Additionally, differential privacy techniques were used to obfuscate data, while still allowing the computation of attack detection results based on obfuscated data. The SBON developed leverages a private permissioned blockchain network orchestrated with the Hyperledger Fabric framework. To connect the cyber world, which orchestrates the blockchain fabric, and the physical world where the power network resides, we developed a system implementation to enable the secure interaction of the physical world and the abstracted blockchain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Influence of an Interlayer on Dual Hydraulic Fractures Propagation

Multi-cluster hydraulic fracturing of long-range horizontal wells is an approach for enhancing the productivity of low-permeability shale reservoirs. In this study, RFPA-Petrol (rock failure process analysis on petroleum problems) is applied for modeling hydraulic fracture propagation in multilayered formations. RFPA-Petrol based on coupled hydraulic-mechanical-damage (HMD) modeling was first tested by modeling a laboratory scale experiment on a physical (cement) model with a single completion. The modeling demonstrated the capability of RFPA-Petrol for simulating hydraulic fracture propagation. Then, we used RFPA-Petrol to investigate how the difference in material properties between oil-bearing layers and interlayers and the fracturing fluid properties influence the propagation of dual fractures in multilayered laboratory-scale models. In this case, the models with geological discontinuities in the vertical direction are strongly heterogeneous and RFPA-Petrol simulations successfully modeled the fracture configurations.

58 GEOSCIENCES↗
DOE OSTI2020

Development of aerosol activation in the double-moment Unified Model and evaluation with CLARIFY measurements

Representing the number and mass of cloud and aerosol particles independently in a climate, weather prediction or air quality model is important in order to simulate aerosol direct and indirect effects on radiation balance. Here we introduce the first configuration of the UK Met Office Unified Model in which both cloud and aerosol particles have “double-moment” representations with prognostic number and mass. The GLObal Model of Aerosol Processes (GLOMAP) aerosol microphysics scheme, already used in the Hadley Centre Global Environmental Model version 3 (HadGEM3) climate configuration, is coupled to the Cloud AeroSol Interacting Microphysics (CASIM) cloud microphysics scheme. We demonstrate the performance of the new configuration in high-resolution simulations of a case study defined from the CLARIFY aircraft campaign in 2017 near Ascension Island in the tropical southern Atlantic. We improve the physical basis of the activation scheme by representing the effect of existing cloud droplets on the activation of new aerosol, and we also discuss the effect of unresolved vertical velocities. We show that neglect of these two competing effects in previous studies led to compensating errors but realistic droplet concentrations. While these changes lead only to a modest improvement in model performance, they reinforce our confidence in the ability of the model microphysics code to simulate the aerosol–cloud microphysical interactions it was designed to represent. Capturing these interactions accurately is critical to simulating aerosol effects on climate.

54 ENVIRONMENTAL SCIENCES↗