SEARCH · Engineering Papers

Results for “numeric profiling”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2026

TRANSP-TGLF core predictive modeling of the JET DT baseline scenario

In recent years, an intense modeling activity has been focused on preparing and analyzing the second JET Deuterium–Tritium (D–T) experimental campaign DTE2. Among the numerous scientific outcomes of this campaign was the unique opportunity to test and validate the state of the art of modeling tools with fusion-relevant DT plasmas using the full metallic ITER-like wall in different scenarios. This work reports on the core predictive modeling of plasma density, electron and ion temperatures (n e , T e , T i ) performed using TRANSP (Pankin et al 2025 Comput. Phys. Commun. 312 109611) coupled with Trapped Gyro-Landau Fluid (TGLF)-SAT2 (Kinsey et al 2008 Phys. Plasmas 15 055908), (Staebler et al 2021 Nucl. Fusion 61 116007) for the JET D–T baseline scenario (I p = 3.5 MA, q 95 = 3, β N < 2, with pellet pacing) (Garzotti et al 2025 Plasma Phys. Control. Fusion 67 075011). The sensitivity to different input parameters, as the $\vec{E}$ x $\vec{B}$ shear parameterization and the values of the kinetic quantities at the boundary of the prediction domain (ρ = 0.85) has been assessed, identifying the confidence interval of the prediction results. In particular, the dependence of the electron density profile on the particle source parameters has been studied, identifying the ionization source as the main cause for the density gradient under-prediction obtained by TRANSP-TGLF, reported in (Hyun-Tae et al 2023 Nucl. Fusion 63 112004).

JET↗

DOE OSTI2020

An Experimental Feasibility Study on Applying Neutron Tomography to Encapsulated Spent Nuclear Fuel - 20050

Visual inspection makes easier to ensure the integrity and safety of spent nuclear fuel (SNF) than any classical techniques. Various classical techniques have been applied but there are no reliable methods to qualitatively and quantitatively verify spent fuel in dry storage. Thus, the present authors have developed the prototype safeguards apparatus for dry storage employing the array of He-4 gas scintillation detectors (S670E, Arktis Radiation Detectors Ltd., Switzerland), newly designed to simultaneously measure thermal and fast neutrons without any moderators. The S670E detector has a cylindrical shape with a diameter of 52 mm and active length of 600 mm (total length: 875 mm). The detector is filled by He-4 gas with an approximate pressure of 180 bar for fast neutron detection, and its inner wall is coated by Li-6 for thermal neutron detection. The scintillation lights generated via Li-6 nuclear reaction and elastic scattering are collected by 24 SiPMs linearly paired at the center of the detector. The detector delivers a TTL (Transistor-Transistor Logic) output for pulse readout and UART (Universal Asynchronous Receiver Transmitter) for device control. In order to assess feasibility of the apparatus, an experimental system has been designed, built, and optimized via computational studies. Cf-252 neutron sources and linearly arrayed detectors, working as a single detector, were occupied for this study due to the difficulties in working with SNF. The laboratory scale cask (diameter: 0.67 m, height: 1.5 m), minimized by a factor of 10 compared to the actual thickness of a commercial TN-32 cask, was also manufactured. The detector array was designed to rotate the lab-scale cask and obtain 36 image profiles at every 10 degrees. All profiles were aligned in single frame image called a sinogram, and the cross-sectional image was then fabricated by the inverse radon transform algorithm. These experiments have been repeated with different configurations and numbers of sources. Some gamma-ray sources were also measured with neutron sources in order to distinguish between neutron and gamma-ray pulses. Basically, a He-4 detector is designed to run on Linux OS so it is difficult to directly apply to Windows-based equipment widely used in S. Korea. Therefore, a new data acquisition board working on Windows OS was designed and built. The board mainly consists of FPGA (Field Programmable Gate Array) and SoC (System on Chip) for TTL pulse readout, sorting measured data, and transferring data to a user interface. In conclusion, the tomographic system at lab-scale has shown considerable potential to detect a partial or gross defect of encapsulated assemblies in dry storage. Next steps of this study will be to 1) repeatedly carry out experiments to demonstrate scientific reliability and validity, and 2) numerically integrate signals with weight factors to enhance the image quality since the suggested system based on passive interrogation method requires longer measurement time. Finally, the system will apply to a commercial dry storage phased out soon in S. Korea. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2022

Finite flame thickness effects on Kolmogorov-Petrovsky-Piskunov turbulent burning velocities

KPP (Kolmogorov-Petrovsky-Piskunov) solutions of the reaction-diffusion equation have application in various physical phenomena occurring in biology, ecology, and reacting flows. In particular, these solutions are commonly used in turbulent combustion to scale turbulent burning velocities. Subject to certain conditions on reaction rate profile through the flame brush and turbulent diffusivity, this theory relates the turbulent burning velocity to the derivative of the reaction rate ($\tilde{ω}$) at the leading edge of the flame brush ($d\tilde{ω}$/$d\tilde{c}$ | $\tilde{c} =$0 ). Such waves are often referred to as “pulled fronts.” However, turbulent flames never actually satisfy the KPP conditions for a pulled front, as the turbulent flame brush, parametrized here by the thickness δ t , consists of an ensemble of laminar flamelets of thickness δ, where ε=δ/δ t $\ll$ 1 is very small, but nonzero, and $d\tilde{ω}$ /$d\tilde{c}$ tends to zero at the brush leading edge for high activation energy, combustion-type kinetics. Here, this paper analyzes these effects on KPP wave solutions, parametrized by ε=δ/δ t and Zeldovich number Ze focusing on whether turbulent flames retain their pulled front character and what the correction to the KPP wave speed is. Variational solutions of the reaction-diffusion equation show that the solution can be expanded in powers of 1/|ln⁡ε|. Both numerical and asymptotic results are presented, showing that the wave still exhibits pulled front solutions but with significant corrections to the KPP result. The leading order correction is of the form |ln⁡ε| –2 and independent of Ze. Higher order corrections are function of both ε and Ze. However, the dominant factor influencing the wave speed correction is due to the finite ε, with Ze exhibiting a weaker effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series: Supporting Data

This dataset consists of soil temperature time series that were used to estimate soil thermal diffusivity and its uncertainty trough the probabilistic modelling approach developed and presented in the article named "Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series" and published in Earth Surface Dynamics. There are two compressed (.zip) files that contains synthetic (Synthetic_soiltemp_Teller.zip) and field (Field_soiltemp_Teller.zip) data. There is one MATLAB file that requires MATLAB to execute but any text editor can open it. The Synthetic_soiltemp_Teller.zip file includes 5 comma-delimited data files (.csv) each of which contains soil temperature time series generated through forward modeling (i.e., heat-conduction process in a heterogeneous medium using an explicit finite difference method) to mimic various types of temperature gradients, trend and fluctuations. These synthetic soil temperatures were used to investigate the impact of different environmental conditions on the uncertainty of thermal diffusivity estimates. The Field_soiltemp_Teller.zip file contains 28 comma-delimited data files (.csv) out of which (a) 27 files includes soil temperature time series recorded from 27 temperature probes located in a site along Teller Road about 40 km northwest of Nome, Alaska (64.72°N, 165.94°W), (b) one includes the name and coordinates of the 27 probes. These field soil temperatures were used to infer soil thermal diffusivity at numerous locations and depths in a discontinuous permafrost environment, and to evaluate the links between the estimated soil thermal diffusivity values and soil physical properties. The comma-delimited data files (.csv) of the synthetic and field soil temperature time series includes date and time (UTC) in the first column and soil temperature from 5 cm below the ground surface to 1.05 m depth (with 5 or 10 cm spacing between sensors) in the other columns. The measurements were acquired every 15 minutes. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Interaction between cloud–radiation, atmospheric dynamics and thermodynamics based on observational data from GoAmazon 2014/15 and a cloud-resolving model

Observational meteorological data from the field experiment GoAmazon 2014/15 and data from numerical simulations with the cloud-resolving model (CRM) called the System for Atmospheric Modeling (SAM) are used to study the interaction between the cloudiness–radiation as well as the atmospheric dynamics and thermodynamics variables for a site located in the central Amazon region (–3.2°S, –60.6°W) during the wet and dry periods. The main aims are to (a) analyze the temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux as well as (b) to determine the relationship between the integrated cloud fraction, radiative fluxes and large-scale variable anomalies as a function of the previous day's average. The temporal series of the integrated cloud fraction, precipitation rate and downward shortwave flux from SAM simulations showed physical consistency with the observations from GoAmazon 2014/15. Shallow and deep convection clouds show to have a meaningful impact on radiation fluxes in the Amazon region during wet and dry periods. Anomalies of large-scale variables (relative to the previous day's average) are physically associated with cloud formation, evolution and dissipation. SAM consistently simulated these results, where the cloud fraction vertical profile shows a pattern very close to the observed data (cloud type). Additionally, the integrated cloud fraction and large-scale variable anomalies, as a function of the previous day's average, have a good correlation. These results suggest that the memory of the large-scale dynamics from the previous day can be used to estimate the cloud fraction as well as the water content, which is a variable of the cloud itself. In general, the SAM satisfactorily simulated the interaction between cloud–radiation as well as dynamic and thermodynamic variables of the atmosphere during the periods of this study, being able to obtain atmospheric variables that are impossible to obtain in an observational way.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Benzo[a]pyrene (BaP) metabolites predominant in human plasma following escalating oral micro-dosing with [ 14 C]-BaP

Benzo[a]pyrene (BaP) is formed by incomplete combustion of organic materials (petroleum, coal, tobacco, etc.). BaP is designated by the International Agency for Research on Cancer as a group 1 known human carcinogen; a classification supported by numerous studies in preclinical models and epidemiology studies of exposed populations. Risk assessment relies on toxicokinetic and cancer studies in rodents at doses 5–6 orders of magnitude greater than average human uptake. Using a dose–response design at environmentally relevant concentrations, this study follows uptake, metabolism, and elimination of [ 14 C]-BaP in human plasma by employing UPLC - accelerator mass spectrometry (UPLC-AMS). Volunteers were administered 25, 50, 100, and 250 ng (2.7–27 nCi) of [ 14 C]-BaP (with interceding minimum 3-week washout periods) with quantification of parent [ 14 C]-BaP and metabolites in plasma measured over 48 h. [ 14 C]-BaP median T max was 30 min with C max and area under the curve (AUC) approximating dose-dependency. Marked inter-individual variability in plasma pharmacokinetics following a 250 ng dose was seen with 7 volunteers as measured by the C max (8.99 ± 7.08 ng × mL -1 ) and AUC0-48hr (68.6 ± 64.0 fg × hr -1 × mL -1 ). Approximately 3–6% of the [ 14 C] recovered (AUC 0-48 hr ) was parent compound, demonstrating extensive metabolism following oral dosing. Metabolite profiles showed that, even at the earliest time-point (30 min), a substantial percentage of [ 14 C] in plasma was polar BaP metabolites. The best fit modeling approach identified non-compartmental apparent volume of distribution of BaP as significantly increasing as a function of dose (p = 0.004). Bay region tetrols and dihydrodiols predominated, suggesting not only was there extensive first pass metabolism but also potentially bioactivation. AMS enables the study of environmental carcinogens in humans with de minimus risk, allowing for important testing and validation of physiologically based pharmacokinetic models derived from animal data, risk assessment, and the interpretation of data from high-risk occupationally exposed populations.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

The NanoSIMS-HR: The Next Generation of High Spatial Resolution Dynamic SIMS

The high lateral resolution and sensitivity of the NanoSIMS 50 and 50L series of dynamic SIMS instruments have enabled numerous scientific advances over the past 25 years. Here, in this study, we report on the NanoSIMS-HR, the first major upgrade to the series, and analytical tests in a suite of sample types, including an aluminum sample containing silicon crystals, microalgae, and plant roots colonized with a symbiotic fungus. Significant improvements have been made in the Cs + ion source, high voltage (HV) control, stage reproducibility, and other aspects of the instrument that affect performance. The modified design of the NanoSIMS-HR thermal-ionization Cs + source enables a 5 pA primary ion beam to be focused into a 100 nm spot, a ~2.5-fold increase compared to Cs + sources on previous instruments (~2 pA at 100 nm). The brightness of the new Cs + source enables an ultimate lateral resolution as high as 30 nm and improved detection limits for a given analysis area. Sample stage movement accuracy is higher than 500 nm, enabling many-fold higher throughput automated analyses. With the new HV control, the primary ion beam impact energy can be reduced from 16 to 2 keV, which enables higher depth resolution during depth profiling (a 2-fold improvement), albeit with a 5-fold decrease in lateral resolution. In the NanoSIMS-HR, the secondary ion column and detection system are identical to those used in the previous series, and the isotopic analysis performance is as precise as in previous NanoSIMS instruments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Collapsar Black Holes Are Likely Born Slowly Spinning

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power γ-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with a large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. and GRB observations favor stellar density profiles that yield an accretion rate of $\dot m\approx 10^{-2} M_{\odot}s^{-1}$, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin of a ≃ 0.2, or a ≃ 0.5 for wobbling jets, with the main uncertainty originating in the unknown γ-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anticorrelation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al., we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to a ≲ 0.2, consistent with gravitational-wave observations. We verify our results by performing the first 3D general-relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Plasma Facing Components with Capillary Porous System and Liquid Metal Coolant Flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [Khodak and Maingi (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile, as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a Fusion Nuclear Science Facility [Kessel et al. (2019)]: 10T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

High Heat Flux↗

DOE OSTI2021

Simulation of lithium transport using the BOUT++ framework

A numerical model that calculates the collisional interactions between the lithium atoms from a lithium pellet and the background plasmas has been upgraded. The ion density ($N_t$), electron temperature ($T_e$), ion temperature ($T_i$) and parallel ion velocity ($V_{∥, i}$) are used to characterize the background plasmas. The lithium atom density ($N^{a}_{Li}$) and parallel velocity ($V_{∥,a}$) of lithium atoms evolve with time. For each lithium ion, the density ($N_{Li^{n+}}$), temperature ($T_{Li^{n+}}$) and parallel velocity ($V_{∥, Li^{n+}}$) are self-consistently calculated. A C-mod lower single null equilibrium is used to generate the grid for the BOUT++ simulation. The lithium atoms can be fully ionized to $Li^{3+}$ in ~2 μs. The rapid radial and poloidal expansion of the lithium ions are found in the simulation. After the collision interaction process, the electron temperature rapidly decreases at the pellet location; then, it rapidly poloidally expands, and the temperature at the pellet location starts to recover. The electron pressure increases at the pellet location despite the decrease in electron temperature because of the extra electrons from the lithium ionization. The ion pressure profile decreases in the pellet location due to the decrease in ion temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗

DOE OSTI2023

Earth-system-model evaluation of cloud and precipitation occurrence for supercooled and warm clouds over the Southern Ocean's Macquarie Island

Abstract. Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Portable, heterogeneous ensemble workflows at scale using libEnsemble

libEnsemble is a Python-based toolkit for running dynamic ensembles, developed as part of the DOE Exascale Computing Project. The toolkit utilizes a unique generator–simulator–allocator paradigm, where generators produce input for simulators, simulators evaluate those inputs, and allocators decide whether and when a simulator or generator should be called. The generator steers the ensemble based on simulation results. Generators may, for example, apply methods for numerical optimization, machine learning, or statistical calibration. libEnsemble communicates between a manager and workers. Flexibility is provided through multiple manager–worker communication substrates each of which has different benefits. These include Python’s multiprocessing, mpi4py, and TCP. Multisite ensembles are supported using Balsam or Globus Compute. We overview the unique characteristics of libEnsemble as well as current and potential interoperability with other packages in the workflow ecosystem. We highlight libEnsemble’s dynamic resource features: libEnsemble can detect system resources, such as available nodes, cores, and GPUs, and assign these in a portable way. These features allow users to specify the number of processors and GPUs required for each simulation; and resources will be automatically assigned on a wide range of systems, including Frontier, Aurora, and Perlmutter. Such ensembles can include multiple simulation types, some using GPUs and others using only CPUs, sharing nodes for maximum efficiency. We also describe the benefits of libEnsemble’s generator–simulator coupling, which easily exposes to the user the ability to cancel, and portably kill, running simulations based on models that are updated with intermediate simulation output. We demonstrate libEnsemble’s capabilities, scalability, and scientific impact via a Gaussian process surrogate training problem for the longitudinal density profile at the exit of a plasma accelerator stage. In conclusion, the study uses gpCAM for the surrogate model and employs either Wake-T or WarpX simulations, highlighting efficient use of resources that can easily extend to exascale.

Dynamic ensembles↗

DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗