SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

DOE OSTI2022

Effect of heavy ion irradiation dose rate and temperature on α' precipitation in high purity Fe-18%Cr alloy

Cr-rich alpha prime precipitates (CrRP) induce hardening and embrittlement of FeCr alloys, but the kinetics of CrRP formation due to particle irradiation are not well understood. In this study, Fe18wt.%Cr alloy in solid solution state and pre-aged to produce relatively coarse CrRP was irradiated with 8 MeV Fe ions. The irradiation conditions involved two midrange doses of 0.37 and 3.7 displacements per atom (dpa), a wide range of dose rates (10 -5 –10 -3 dpa/s) and temperatures (300–450 °C). The distributions of CrRP after irradiation were studied with atom probe tomography (APT). The critical irradiation conditions to suppress CrRP formation were identified as 300 °C and 10-3 dpa/s; CrRP formation occurred readily at lower dose rates or higher temperatures. From 0.37 to 3.7 dpa, CrRP were observed to slightly grow at 350 °C and strongly coarsen at 450 °C. Specimens with pre-existing CrRP evolved into a similar precipitate distribution as detected after ion irradiation on solution annealed specimens at 300–350 °C to 0.37 dpa, indicating that the precipitate microstructure approaches a quasi-equilibrium for doses <1 dpa. Limited shrinking of pre-existing CrRP was observed after irradiation at 450 °C to 0.37 dpa, indicating a higher recovery rate at this temperature. Additionally, the evolution of CrRP is quantitatively explained by employing corrections to the historic Nelson-Hudson-Mazey precipitate stability model, and a radiation modified precipitation mechanism is proposed to account for the competition between radiation enhanced diffusion and ballistic dissolution which results in the modifications on both size and solute concentration of CrRP.

36 MATERIALS SCIENCE↗

DOE OSTI2019

9Cr steel visualization and predictive modeling

The goals of this work were to develop tools to visualize and characterize materials information, to explore new alloy compositions and/or processing, to better predict tensile strength and other mechanical properties. The 9–12 wt% Cr (9Cr) steel test data were compiled from several sources. About 3000 records with 30 predictor variables represent 82 unique steel alloy compositions, variations of thermo-mechanical processing steps and temperatures (homogenization, normalization and tempering cycles), test conditions and outcomes. Detailed data processing steps such as visualization and exploratory analysis, including univariate and bivariate analysis, different clustering techniques used to segment the data, statistical post-hoc analysis to verify significance of the findings, feature engineering to identify predictors of importance, and predictive modeling with cross-validation were performed. The outcome of analysis at each step was reviewed in the context of the domain knowledge of this class of steel, to see if there were underlying physical mechanisms that explain statistical relationships. Data analytics techniques and their parameters were fine-tuned to facilitate interpretation of the results as aligned with the insights from domain experts. Here, the ensemble predictive modeling using Random Forest regressor and post-hoc means comparison corroborated domain knowledge on the role of Co which is known to increase strength of steel alloys through solid-solution strengthening and to affect diffusion of alloying elements and precipitates. Alloys with Co content of 0.7–8 wt% had significantly higher mean strength than the ones without (while having Cr locked within 10–11 wt% range in either group and keeping other compositional element ratios fixed). End products of the computational techniques exploration are presented as the tools that can be used in iterative workflow of materials development and testing.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Use of the Hugoniot Elastic Limit in Laser Shockwave Experiments to Relate Velocity Measurements

---The Idaho National Laboratory (INL) has been developing a Laser Shock system to characterize interface strength in plate fuel for research reactors around the world [1]. This work is in support of the Unites states High Performance Research Reactor program whose goal is to reduce the risk of proliferation of weapons grade fuel while maintaining reactor performance. The INL has been working with National Research Council Canada (NRC) on this project for the last five years. ---One of the concerns of the project is the difficulty of calibrating and standardizing the laser shock technique. The generation of the stress wave by the laser impinging on the plate surface is complex, with amplitude and shape dependence on numerous factors relating to the temporal and spatial distribution of power in the laser spot, surface treatment, and material properties of the plate. Thus the resulting back-face surface velocity, which is used to measure the time-varying stress field within the plate, is difficult to reproduce across different measurement systems. A technique that can generate or indicate a known shock intensity independent of the generating conditions would be valuable and allow us to directly compare results from different systems. ---A Laser Shock round robin test between the system at NRC and the INL was conducted. The goal of the round robin was to determine if there was any bias between the two systems. During the analysis of the data, it was noted that the Hugoniot Elastic Limit (HEL) was consistent within both system data sets independent of the laser power used. An analytical study under development supports that the HEL in metals could be a robust and simple benchmark to compare stresses generated by different laser shock systems. The HEL is a material property that should be insensitive to differences in the shock laser and surface preparation which enhances the reliability of the resulting HEL velocity measurement. 1. J. M. Lacy, J. A. Smith and B. H. Rabin, "Developing a Laser Shockwave Model for Characterizing Diffusion Bonded Interfaces," The 41st Annual Review of Progress in Quantitative Nondestructive Evaluation," Boise, ID, U.S.A., July 2014.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Local measurement of bulk thermal diffusivity using photothermal radiometry

Here, an experimental methodology using photothermal radiometry is developed for the accurate measurement of bulk thermal diffusivity of nuclear fuels and materials irradiated to high doses. Under these conditions, nuclear fuels, such as uranium oxide, and moderator materials, such as graphite, become friable, which requires characterization techniques that can accommodate irregularly shaped fragments. Photothermal radiometry, a good candidate for this application, involves locally heating a sample by using a laser and measuring the temperature field by monitoring blackbody radiation. The interaction volume for this study, less than a millimeter, is carefully chosen to sample a statistically significant number of large-scale structural features, such as pores and gas filled bubbles, and is small enough that the sample fragments can be treated as a thermal half-space. The thermal diffusivity standards considered in this study cover a range of thermal diffusivities representative of both fresh and spent nuclear fuels. We also consider a sample having a porous microstructure representative of large-scale structures found in materials irradiated to high doses. Our measurement methodology circumvents complex thermal wave models that address optical diffraction, nonlinear transfer function associated with blackbody radiation, and finite sample size effects. Consequently, the large measurement uncertainty associated with modeling these effects can be avoided. While the emphasis here is on nuclear fuels and materials, this measurement approach is well suited to measure thermal transport in a variety of technologically important materials associated with advanced synthesis techniques. Examples range from small, exotic single crystals grown using hydrothermal growth techniques to additively manufactured components having complex geometries.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Fluid age-based analysis of a lifted turbulent DME jet flame DNS [Residence Time-Based Analysis of a Lifted Turbulent DME Jet Flame]

The link between the distribution of fluid residence time and the distribution of reactive scalars is analysed using Direct Numerical Simulation data. Information about the reactive scalar distribution is needed in order to model the reaction terms that appear in Large Eddy and Reynolds-Averaged simulations of turbulent reacting flows. The lifted flame is simulated taking account of multi-step chemistry for dimethyl-ether fuel and differential diffusion. Due to autoignition and flame propagation, the reaction progress increases with residence time. The variation of fluid residence time is evaluated by solving an Eulerian transport equation for the fluid age. The fluid age is a passive scalar with a spatially-uniform source term, meaning that its moments and dissipation rates in turbulent flows can be modelled using closures already established for conserved scalars such as mixture fraction. In combination with the mixture fraction, the fluid age serves as a useful mapping variable to distinguish younger less-reacted fluid near the inlet from older more-reacted fluid downstream. The local fluctuations of mixture fraction and fluid age have strong negative correlation and, building upon established presumed-pdf models for mixture fraction, this feature can be used to construct an accurate presumed-pdf model for the joint mixture fraction/fluid age pdf. It is demonstrated that the double-conditional first-order moment closure combined with the proposed presumed model for the joint pdf of mixture fraction and fluid age gives accurate predictions for unconditional reaction rates – both for pre-ignition radical species produced by low-temperature processes upstream of the flame base, and for major species that are produced at the flame front.

42 ENGINEERING↗

DOE OSTI2019

Theoretical Determination of Size Effects in Zeolite-Catalyzed Alcohol Dehydration

In the upgrading of biomass pyrolysis vapors to hydrocarbons, dehydration accomplishes a primary objective of removing oxygen, and acidic zeolites represent promising catalysts for the dehydration reaction. Here, we utilized density functional theory calculations to estimate adsorption energetics and intrinsic kinetics of alcohol dehydration over H-ZSM-5, H-BEA, and H-AEL zeolites. The ONIOM (our Own N-layered Integrated molecular Orbital and molecular Mechanics) calculations of adsorption energies were observed to be inconsistent when benchmarked against QM (Quantum Mechanical)/Hartree–Fock and periodic boundary condition calculations. However, reaction coordinate calculations of adsorbed species and transition states were consistent across all levels considered. Comparison of ethanol, isopropanol (IPA), and tert-amyl alcohol (TAA) over these three zeolites allowed for a detailed examination of how confinement impacts on reaction mechanisms and kinetics. The TAA, seen to proceed via a carbocationic mechanism, was found to have the lowest activation barrier, followed by IPA and then ethanol, both of which dehydrate via a concerted mechanism. Barriers in H-BEA were consistently found to be lower than in H-ZSM-5 and H-AEL, attributed to late transition states and either elevated strain or inaccurately estimating long-range electrostatic interactions in H-AEL, respectively. Molecular dynamics simulations revealed that the diffusivity of these three alcohols in H-ZSM-5 were significantly overestimated by Knudsen diffusion, which will complicate experimental efforts to develop a kinetic model for catalytic fast pyrolysis.

09 BIOMASS FUELS↗

DOE OSTI2020

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

DOE OSTI2019

Structure and proton conduction in sulfonated poly(ether ether ketone) semi-permeable membranes: a multi-scale computational approach.

The design of polymeric membranes for proton or ionic exchange highly depends on the fundamental understanding of the physical and molecular mechanisms that control the formation of the conduction channels. There is an inherent relation between the dynamical structure of the polymeric membrane and the electrostatic forces that drive membrane segregation and proton transport. Here, we used a multi-scale computational approach to analyze the morphology of sulfonated poly(ether ether ketone) membranes at the mesoscale. A self-consistent description of the electrostatic phenomenon was adopted, where discrete polymer chains and a continuum proton field were embedded in a continuum fluid. Brownian dynamics was used for the evolution of the suspended polymer molecules, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the proton concentration field. We varied the polymer concentration, the degree of sulfonation and the level of confinement to find relationships between membrane structure and proton conduction. Our results indicate that the reduced mobility of polymer chains, at concentrations above overlap, and a moderate degree of sulfonation -i.e., 30% - are essential elements for membrane segregation and proton domain connectivity. These conditions also ensure that the membrane structure is not affected by size or by potential gradients. Importantly, our analysis shows that membrane conductivity and current are linearly dependent on polymer concentration and quadratically dependent on the degree of sulfonation. We found that the optimal polymeric membrane design requires a polymer concentration above overlap and a degree of sulfonation around 50%. These conditions promote a dynamical membrane morphology with a constant density of proton channels. Our results and measurements agree with previous experimental works, thereby validating our model and observations.

Molina, Jarol↗

DOE OSTI2013

Modeling Solute Thermokinetics in LiCl-KCl Molten Salt for Nuclear Waste Separation (NEUP 10-418)

Recovery of actinides is an integral part of a closed nuclear fuel cycle. Pyrometallurgical nuclear fuel recycling processes have been developed in the past for recovering actinides from spent metallic and nitride fuels. The process is essentially to dissolve the spent fuel in a molten salt and then extract just the actinides for reuse in a reactor. Extraction is typically done through electrorefining, which involves electrochemical reduction of the dissolved actinides and plating onto a cathode. Knowledge of a number of basic thermokinetic properties of salts and salt-fuel mixtures is necessary for optimizing present and developing new approaches for pyrometallurgical waste processing. The properties of salt-fuel mixtures are presently being studied, but there are so many solutes and varying concentrations that direct experimental investigation is prohibitively time consuming and expensive (particularly for radioactive elements like Pu). Therefore, there is a need to reduce the number of required experiments through modeling of salt and salt-fuel mixture properties. This project will develop first-principles-based molecular modeling and simulation approaches to predict fundamental thermokinetic properties of dissolved actinides and fission products in molten salts. The focus of the proposed work is on property changes with higher concentrations (up to 5 mol%) of dissolved fuel components, where there is still very limited experimental data. The properties predicted with the modeling will be density, which is used to assess the amount of dissolved material in the salt; diffusion coefficients, which can control rates of material transport during separation; and solute activity, which determines total solubility and reduction potentials used during electrorefining. The work will focus on La, Sr, and U, which are chosen to include the important distinct categories of lanthanides, alkali earths, and actinides, respectively. Studies will be performed using LiCl-KCl salt at the eutectic composition (58 mol% LiCl, 42 mol% KCl), which is used for treating spent EBR-II fuel. The same process being used for EBRII fuel is currently being studied for widespread international implementation. The methods will focus on first-principles and first- principles derived interatomic potential based simulations, primarily using molecular dynamics. Results will be validated against existing literature and parallel ongoing experimental efforts. The simulation results will be of value for interpreting experimental results, validating analytical models, and for optimizing waste separation by potentially developing new salt configurations and operating conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

A Numerical Study of Forbush Decreases with a 3D Cosmic-Ray Modulation Model Based on an SDE Approach

Based on the reduced diffusion mechanism for producing Forbush decreases (Fds) in the heliosphere, we constructed a three-dimensional (3D) diffusion barrier, and by incorporating it into a stochastic differential equation (SDE) based time-dependent, cosmic-ray transport model, a 3D numerical model for simulating Fds is built and applied to a period of relatively quiet solar activity. This SDE model generally corroborates previous Fd simulations concerning the effects of the solar magnetic polarity, the tilt angle of the heliospheric current sheet (HCS), and cosmic-ray particle energy. Because the modulation processes in this 3D model are multi-directional, the barrier’s geometrical features affect the intensity profiles of Fds differently. We find that both the latitudinal and longitudinal extent of the barrier have relatively fewer effects on these profiles than its radial extent and the level of decreased diffusion inside the disturbance. We find, with the 3D approach, that the HCS rotational motion causes the relative location from the observation point to the HCS to vary, so that a periodic pattern appears in the cosmic-ray intensity at the observing location. Correspondingly, the magnitude and recovery time of an Fd change, and the recovering intensity profile contains oscillation as well. Investigating the Fd magnitude variation with heliocentric radial distance, we find that the magnitude decreases overall and, additionally, that the Fd magnitude exhibits an oscillating pattern as the radial distance increases, which coincides well with the wavy profile of the HCS under quiet solar modulation conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

DOE OSTI2015

Algorithms for propagating uncertainty across heterogeneous domains

We address an important research area in stochastic multi-scale modeling, namely the propagation of uncertainty across heterogeneous domains characterized by partially correlated processes with vastly different correlation lengths. This class of problems arise very often when computing stochastic PDEs and particle models with stochastic/stochastic domain interaction but also with stochastic/deterministic coupling. The domains may be fully embedded, adjacent or partially overlapping. The fundamental open question we address is the construction of proper transmission boundary conditions that preserve global statistical properties of the solution across different subdomains. Often, the codes that model different parts of the domains are black-box and hence a domain decomposition technique is required. No rigorous theory or even effective empirical algorithms have yet been developed for this purpose, although interfaces defined in terms of functionals of random fields (e.g., multi-point cumulants) can overcome the computationally prohibitive problem of preserving sample-path continuity across domains. The key idea of the different methods we propose relies on combining local reduced-order representations of random fields with multi-level domain decomposition. Specifically, we propose two new algorithms: The first one enforces the continuity of the conditional mean and variance of the solution across adjacent subdomains by using Schwarz iterations. The second algorithm is based on PDE-constrained multi-objective optimization, and it allows us to set more general interface conditions. The effectiveness of these new algorithms is demonstrated in numerical examples involving elliptic problems with random diffusion coefficients, stochastically advected scalar fields, and nonlinear advection-reaction problems with random reaction rates.

Cho, Heyrim↗

DOE OSTI2014

Temperature dependent transport characteristics of graphene/n-Si diodes

Realizing an optimal Schottky interface of graphene on Si is challenging, as the electrical transport strongly depends on the graphene quality and the fabrication processes. Such interfaces are of increasing research interest for integration in diverse electronic devices as they are thermally and chemically stable in all environments, unlike standard metal/semiconductor interfaces. We fabricate such interfaces with n-type Si at ambient conditions and find their electrical characteristics to be highly rectifying, with minimal reverse leakage current (<10{sup −10} A) and rectification of more than 10{sup 6}. We extract Schottky barrier height of 0.69 eV for the exfoliated graphene and 0.83 eV for the CVD graphene devices at room temperature. The temperature dependent electrical characteristics suggest the influence of inhomogeneities at the graphene/n-Si interface. A quantitative analysis of the inhomogeneity in Schottky barrier heights is presented using the potential fluctuation model proposed by Werner and Güttler.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

ELM–Wet: Inclusion of a Wet–Landunit With Sub–Grid Representation of Eco–Hydrological Patches and Hydrological Forcing Improves Methane Emission Estimations in the E3SM Land Model (ELM)

Wetlands are the largest emitters of biogenic methane (CH 4 ) and represent the highest source of uncertainty in global CH 4 budgets. Here, we aim to improve the realism of wetland representation in the U.S. Department of Energy's Exascale Earth System Model land surface model, ELM, thereby reducing uncertainty of CH 4 flux predictions. We develop an updated version, ELM-Wet, where we activate a separate landunit for wetlands that handles multiple wetland-specific eco-hydrological patch functional types. We introduce more realistic hydrological forcing through prescribing site-level constraints on surface water elevation, which allows resolving different sustained inundation depth for different patches, and if data exists, prescribing inundation depth. We modified the calculation of aerenchyma transport diffusivity based on observed conductance per leaf area for different vegetation types. We use Bayesian Optimization to parameterize CO 2 and CH 4 fluxes in the developed wet-landunit. Site-level simulations of a coastal non-tidal freshwater wetland in Louisiana were performed with the updated model. Eddy covariance observations of CO 2 and CH 4 fluxes from 2012 to 2013 were used to train the model and data from 2021 were used for validation. Patch-specific chamber flux observations and observations of CH 4 concentration profiles in the soil porewater from 2021 were used for evaluation of the model performance. Our results show that ELM-Wet reduces the model's CH 4 emission root mean squared error by up to 33% and is able to represent inter-daily CO 2 and CH 4 flux variability across the wetland's eco-hydrological patches, including during periods of extreme dry or wet conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

TRANSPORT OF SOLAR ENERGETIC PARTICLES ACCELERATED BY ICME SHOCKS: REPRODUCING THE RESERVOIR PHENOMENON

In this work, gradual solar energetic particle (SEP) events observed by multiple spacecraft are investigated with model simulations. Based on a numerical solution of the Fokker-Planck focused transport equation including perpendicular diffusion of particles, we obtained the fluxes of SEPs accelerated by an interplanetary coronal mass ejection driven shock as it propagates outward through the three-dimensional Parker interplanetary magnetic field. The shock is treated as a moving source of energetic particles with an assumed particle distribution function. We look at the time profiles of particle flux as they are observed simultaneously by multiple spacecraft located at different locations. The dependence of particle fluxes on different levels of perpendicular diffusion is determined. The main purpose of our simulation is to reproduce the reservoir phenomenon, during which it is frequently observed that particle fluxes are nearly the same at very different locations in the inner heliosphere, up to 5 AU, during the decay phase of gradual SEP events. The reservoir phenomenon is reproduced in our simulation under a variety of conditions of perpendicular diffusion of particles estimated from the nonlinear guiding center theory (NLGC). As the perpendicular diffusion coefficient increases, the nonuniformity of particle fluxes becomes smaller, making the reservoir phenomenon more prominent. However, if the shock acceleration strength decreases slower than r {sup -2.5} with the radial distance r, the reservoir phenomenon might disappear, with limited perpendicular diffusion constrained by the NLGC theory. Therefore, observation of the reservoir phenomenon in gradual SEP events can be used to test qualitatively theories of particle diffusion and shock acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Opacity Effects on Pulsations of Main-Sequence A-Type Stars

Here, opacity enhancements for stellar interior conditions have been explored to explain observed pulsation frequencies and to extend the pulsation instability region for B-type main-sequence variable stars. For these stars, the pulsations are driven in the region of the opacity bump of Fe-group elements at ~200,000 K in the stellar envelope. Here we explore effects of opacity enhancements for the somewhat cooler main-sequence A-type stars, in which p-mode pulsations are driven instead in the second helium ionization region at ~50,000 K. We compare models using the new LANL OPLIB vs. LLNL OPAL opacities for the AGSS09 solar mixture. For models of two solar masses and effective temperature 7600 K, opacity enhancements have only a mild effect on pulsations, shifting mode frequencies and/or slightly changing kinetic-energy growth rates. Increased opacity near the bump at 200,000 K can induce convection that may alter composition gradients created by diffusive settling and radiative levitation. Opacity increases around the hydrogen and 1st He ionization region (~13,000 K) can cause additional higher-frequency p modes to be excited, raising the possibility that improved treatment of these layers may result in prediction of new modes that could be tested by observations. New or wider convective zones and higher convective velocities produced by opacity increases could also affect angular momentum transport during evolution. More work needs to be done to quantify the effects of opacity on the boundaries of the pulsation instability regions for A-type stars.

79 ASTRONOMY AND ASTROPHYSICS↗