Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Emergence of a nocturnal low-level jet from a broad baroclinic zone

An analytical model is presented for the generation of a Blackadar-like nocturnal low-level jet in a broad baroclinic zone. The flow is forced from below (flat ground) by a surface buoyancy gradient and from above (free atmosphere) by a constant pressure gradient force. Diurnally-varying mixing coefficients are specified to increase abruptly at sunrise and decrease abruptly at sunset. With attention restricted to a surface buoyancy that varies linearly with a horizontal coordinate, the Boussinesq-approximated equations of motion, thermal energy, and mass conservation reduce to a system of one-dimensional equations that can be solved analytically. Sensitivity tests with southerly jets suggest that (i) stronger jets are associated with larger decreases of the eddy viscosity at sunset (as in Blackadar theory), (ii) the nighttime surface buoyancy gradient has little impact on jet strength, and (iii) for pure baroclinic forcing (no free-atmosphere geostrophic wind), the nighttime eddy diffusivity has little impact on jet strength, but the daytime eddy diffusivity is very important and has a larger impact than the daytime eddy viscosity. The model was applied to a jet that developed in fair weather conditions over the Great Plains from southern Texas to northern South Dakota on 1 May 2020. The ECMWF Reanalysis v5 (ERA5) for the afternoon prior to jet formation showed that a broad north-south-oriented baroclinic zone covered much of the region. The peak model-predicted winds were in good agreement with ERA5 winds and lidar data from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) central facility in north-central Oklahoma.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Comparison of Adverse Meteorology Inputs for Emergency Management and Safety Basis Dispersion Modeling Applications at SRS

The MELCOR Accident Consequence Calculation System (MACCS), used for safety basis accident analyses at SRS, implements adjustment factors to consider the increased vertical turbulence associated with local surface roughness and the increased lateral turbulence associated with time-based plume meander. The two models used for SRS emergency planning radiological consequence assessment, i.e., Hotspot and Puff/Plume, currently do not have an option to input or apply diffusion adjustments for meander or surface roughness. Adverse meteorology conditions at SRS have been determined as E stability class and 1.3 meters per second wind speed for ground level releases. This is the input used for emergency planning models. Contrastingly safety basis analyses use a different input, representing adverse meteorology values as F stability class and wind speeds of 1.3 meters per second (for a ground level release). This results in an apparent discrepancy between the two approaches. In this report we document the differing modeling assumptions, discuss the logic behind each modeling approach used to represent local roughness and plume meander conditions, and numerically compare dispersion coefficients and concentration estimates obtained from each set of model inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improving a Multilevel Turbulence Closure Model for a Shallow Lake in Comparison With Other 1-D Models

Lakes differ from lands in water availability, heat capacity, albedo, and roughness, which affect local surface-atmospheric interactions. This study modified a multilevel upper ocean model (UOM) for lake applications and evaluated its performance in Lake Taihu (China) with comprehensive measurements against three popular one-dimensional (1-D) lake models. These models were based on different concepts, including the self-similarity (FLake), the wind-driven eddy diffusion (LISSS), the k-ε turbulence closure (SIMSTRAT), and a simplified turbulence closure (UOM). The surface flux scheme in these models was unified to exclude the discrepancies in representing air-lake exchanges. All models in their default formulations presented obvious cold water temperature biases and largely underestimated the lake surface temperature (LST) diurnal range. For each model, these deficiencies were significantly reduced by incorporating new physics schemes or calibrated tunable parameters based on systematic sensitivity tests. The primary modifications for UOM included (1) a new scheme of decreased surface roughness lengths to better characterize the shallow lake, (2) a solar radiation penetration scheme with increased light extinction coefficient and surface absorption fraction to account for the high water turbidity, and (3) turbulent Prandtl number increased by a factor of 20 to reduce the turbulent vertical mixing. All other models were improved in these three aspects (roughness, extinction, and mixing) within their original formulations. Given these improvements, UOM showed superior performance to other models in capturing LST diurnal cycle and daily to seasonal variations, as well as summer-autumn vertical stratification changes. The new UOM is well suited for application in shallow lakes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Conservation laws and effective hadronization models

Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗
DOE OSTI2022

Reduced models for ETG transport in the tokamak pedestal

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2021

Local wave-number model for inhomogeneous two-fluid mixing

We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Generation of localized reactor point kinetics parameters using coupled neutronic and thermal fluid models for pebble-bed reactor transient analysis

The systems analysis of anticipated operating occurrences and design basis accidents for pebble-bed reactor systems requires knowledge of neutron point kinetics equations (PKE) parameters. Typically, the generation of PKE parameters is performed in a global manner using standalone neutronics calculations, without the inclusion of thermal fluid distributions. We utilize Griffin and Pronghorn for generating global and local PKE parameters which includes the use of thermal fluid distributions to account for localized effects. This work establishes a methodology for calculating PKE parameters for a pebble bed reactor with a coupled neutronics/thermal fluids analysis. PKE parameters generated on a global and local basis are compared against a diffusion solve for a typical load-following transient. Locally-generated neutron kinetic parameters are able to reduce the maximum error in the transient power level from 5% to below 1.5%; along with this, bulk temperature errors were reduced from 12 K to 4 K.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Brazing characteristics, microstructure, and wettability of laser powder bed fusion additive manufactured GRCop-84 compared to CuCrZr and OFC, and brazing to titanium-zirconium-molybdenum alloy limiters

Laser Powder Bed Fusion (L-PBF) of Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, produces a fully dense, high conductivity alloy with a yield strength of 500 MPa and ultimate tensile strength (UTS) of 740 MPa with 20% elongation; superior to other competing copper alloys. Braze wetting characteristics of GRCop-84 with Ag-Cu-X, and Au-Cu brazes were similar to CuCrZr, but less than oxygen free copper. No difference in wetting was observed between infill and surface contour areas in L-PBF GRCop-84. Wet sanding to 240 grit (R a =0.24 µm) was considered the optimal surface condition. Silver diffusing through GRCop-84 depleted Cr 2 Nb precipitates from the copper grain and deposited agglomerations of coarsened precipitates within silver-rich regions of intergranular diffusion once a density threshold was reached. Microstructure modification was minimized with 50Au-50Cu braze implying that silver caused precipitate coarsening and agglomeration, and not high temperature exposure. Coarsened precipitates were observed on the surface within braze pools implying a contribution to braze wetting. Palcusil-25, Ticusil, CuSil-ABA, and 50Au-50Cu brazes were suitable for brazing to unplated Titanium-Zirconium-Molybdenum (TZM), while sulfamate nickel plating to allows wetting with CuSil or other non-active brazes. Additionally, vacuum brazing techniques were developed to join a 1 mm thick layer of TZM to the front of additive manufactured GRCop-84 waveguides considering the brazing characteristics of both GRCop-84, TZM, and internal stress from the difference in coefficient in thermal expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Synthesis, structure, linear and nonlinear optical properties of noncentrosymmetric quaternary diamond-like semiconductors, Cu 2 ZnGeSe 4 (CZGSe) and the novel Cu 4 ZnGe 2 Se 7

In this work, two quaternary diamond-like semiconductors (DLSs) in the Cu-Zn-Ge-Se system, the known Cu 2 ZnGeSe 4 (also referred to as CZGSe) and the new Cu 4 ZnGe 2 Se 7 , are compared in terms of their crystal structures, electronic structures and physicochemical properties. Both compounds were prepared by high-temperature, solid-state synthesis at 800 degrees C. Single crystal X-ray diffraction was used to determine the structure of Cu 4 ZnGe 2 Se 7 . The structures of both Cu 2 ZnGeSe 4 and Cu 4 ZnGe 2 Se 7 can be considered as derivatives of cubic diamond. Cu 4 ZnGe 2 Se 7 , with a reduced symmetry due to the cation ordering pattern, adopts the Cu 4 NiSi 2 S 7 structure type, with space group C2. The corner-sharing tetrahedra in Cu 4 ZnGe 2 Se 7 are slightly distorted because the charge for some of the S 2 - anions is not compensated by the first-nearest-neighbor cations. Rietveld refinements using synchrotron X-ray powder diffraction data were used to assess the phase purity of the samples and confirm the bulk structural behavior. Diffuse reflectance UV/Vis/NIR spectroscopy shows that Cu 2 ZnGeSe 4 and Cu 4 ZnGe 2 Se 7 have direct optical bandgaps of 1.38 and 0.91 eV, respectively. Electronic structure calculations implementing density functional theory confirm the direct bandgap for Cu 4 ZnGe 2 Se 7 , with a calculated value of 0.62 eV. Both compounds are air stable, thermally stable up to relatively high temperatures, undergo phase transitions and have wide windows of optical clarity. The Kurtz-Perry powder technique was used to determine the second harmonic generation (SHG) responses using a commercial AgGaSe 2 standard. While Cu 4 ZnGe 2 Se 7 displays a weak SHG response, Cu 2 ZnGeSe 4 exhibits a response that is greater than several benchmark materials, with a large SHG coefficient, χ (2) , of 43 ± 6 pm/V at λ = 2900 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solving Inverse Stochastic Problems from Discrete Particle Observations Using the Fokker--Planck Equation and Physics-Informed Neural Networks

The Fokker--Planck (FP) equation governing the evolution of the probability density function (PDF) is applicable to many disciplines, but it requires specification of the coefficients for each case, which can be functions of space-time and not just constants and hence require the development of a data-driven modeling approach. When the data available is directly on the PDF, there exist methods for inverse problems that can be employed to infer the coefficients and thus determine the FP equation and subsequently obtain its solution. Herein, we address a more realistic scenario, where only sparse data are given on the particles' positions at a few time instants, which are not sufficient to accurately construct directly the PDF even at those times from existing methods, e.g., kernel estimation algorithms. To this end, we develop a general framework based on physics-informed neural networks (PINNs) that introduces a new loss function using the Kullback--Leibler divergence to connect the stochastic samples with the FP equation to simultaneously learn the equation and infer the multidimensional PDF at all times. In particular, we consider two types of inverse problems, type I, where the FP equation is known but the initial PDF is unknown, and type II, in which, in addition to the unknown initial PDF, the drift and diffusion terms are also unknown. In both cases, we investigate problems with either Brownian or Lévy noise or a combination of both. Here, we demonstrate the new PINN framework in detail in the one-dimensional (1D) case, but we also provide results for up to five dimensions demonstrating that we can infer both the FP equation and dynamics simultaneously at all times with high accuracy using only very few discrete observations of the particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Ce- and Ni-Codoped Double PrBaMn 2 O 5 Perovskite as a Ceramic SOFC Anode

This study explores the efficacy of cerium introduction, both on its own and in combination with nickel, into PrBaMn 2 O 5+δ (PBM) structures to enhance solid oxide fuel cell (SOFC) anodes. We synthesized Pr 1–x Ce x BaMn 2 O 5+δ compositions for x values of 0.05 (PrCe5) and 0.1 (PrCe10), as well as a nickel-doped variant (PrCe5Ni), assessing their performance under H 2 – 3% H 2 O reducing conditions pertinent to SOFC anode operations. Our findings reveal that the PrCe5 composition exhibits a thermal expansion coefficient (TEC) that not only improves upon that of the Ce-free counterpart but also aligns closely with the TEC standards of prevalent SOFC electrolytes. A notable advancement was achieved with the application of a 3.5 μm gadoliniadoped ceria (GDC) buffer layer through physical vapor deposition, effectively mitigating chemical interactions between PBM-based anodes and yttria-stabilized zirconia (YSZ) electrolytes, a concern highlighted by in situ neutron diffraction analyses. Electrochemical impedance spectroscopy, conducted over 220 h in a H 2 –3% H 2 O atmosphere at 750 °C, demonstrated that the optimal 3.5 μm thickness of the GDC buffer layer significantly minimizes the area-specific resistance (ASR) degradation rate to 0.002 Ω cm 2 /h, markedly outperforming both thinner (1 μm) and thicker (8 μm) GDC layers, which showed higher degradation rates of 0.15–0.2 Ω cm 2 /h due to the diffusion of Ba ions or delamination. Moreover, cerium doping fosters superior microstructural stability and obviates barium diffusion, thereby suggesting an enhanced durability of the doped anodes over their lifespan. Integrating nickel into the PrCe5 structure halved the ASR to 0.5 Ω cm 2 at 750 °C, situating it well within the ideal performance range for SOFC anodes. Furthermore, the enhancement brought about by simultaneously doping PBM with cerium and nickel, which fundamentally relies on the critical contributions of defect chemistry and crystal structure, highlights the significance of these fields in creating sophisticated materials for energy related applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Inference of Experimental Radial Impurity Transport on Alcator C-Mod: Bayesian Parameter Estimation and Model Selection.

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. We find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Unique impacts of strong and westward-extended western Pacific subtropical high on ozone pollution over eastern China

As a subtropical anticyclonic high-pressure system that typically forms over the northwestern Pacific Ocean in summer, the Western Pacific subtropical high (WPSH) affects meteorological conditions and ozone pollution in China. The relationship between maximum daily 8-h average ozone (MDA8 O 3 ) concentrations and the extremely strong and westward-extended WPSH occurred in 2022 is investigated using observations, reanalysis data and GEOS-Chem model simulations. The intensity of WPSH has a significant positive correlation with MDA8 O3 over southern China during July-August in 2022, with a correlation coefficient of +0.44, but the correlation is negative (–0.40) in northern China. During the strong WPSH days, MDA8 O 3 increased by 16.5µgm -3 (16.4% relative to July-August average) over southern China and decreased by 19.0µgm -3 (14.5%) in northern China compared to the weak WPSH days. The unique dipole pattern in the relationship between ozone levels and the WPSH in 2022 exhibited a contrast to that during 2015–2021. The difference is primarily due to the extremely strong WPSH intensity and its unusual westward expansion in 2022. In this case, an anomalous anticyclone at 500 hPa dominates over southern China, which creates conditions conducive for ozone formation and accumulation. The anticyclone weakened horizontal winds and reduced the dispersion of ozone, alongside a high temperature and low relative humidity, which favored the chemical production of ozone. In contrast, abnormal northerly winds enhanced ozone diffusion in northern China and the low temperature reduced ozone chemical production. Here, this study reveals the mechanism for the significant impact of strong and westward-extended WPSH on ozone concentrations over China, emphasizing the role of the WPSH location in modulating meteorology and ozone levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Inference of experimental radial impurity transport on Alcator C-Mod: Bayesian parameter estimation and model selection

We present a fully Bayesian approach for the inference of radial profiles of impurity transport coefficients and compare its results to neoclassical, gyrofluid and gyrokinetic modeling. Using nested sampling, the Bayesian Impurity Transport InferencE (BITE) framework can handle complex parameter spaces with multiple possible solutions, offering great advantages in interpretative power and reliability with respect to previously demonstrated methods. BITE employs a forward model based on the pySTRAHL package, built on the success of the well-known STRAHL code [Dux, IPP Report, 2004], to simulate impurity transport in magnetically-confined plasmas. In this paper, we focus on calcium (Ca, Z=20) Laser Blow-Off injections into Alcator C-Mod plasmas. Multiple Ca atomic lines are diagnosed via high-resolution X-ray Imaging Crystal Spectroscopy and Vacuum Ultra-Violet measurements. We analyze a sawtoothing I-mode discharge for which neoclassical and turbulent (quasilinear and nonlinear) predictions are also obtained. Here, we find good agreement in diffusion across the entire radial extent, while turbulent convection and density profile peaking are estimated to be larger in experiment than suggested by theory. Efforts and challenges associated with the inference of experimental pedestal impurity transport are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗