Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local minimum”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Best Practices for Shielding Analyses of Activated Metals and Spent Resins from Reactor Operation

This report investigated best practices for performing shielding evaluations of Type B waste packages, as defined in 10 CFR Part 71 regulations on packaging and transportation of radioactive material, or packages for which the contents are not defined before loading and may include a broad range of nuclides, geometries, and non-fuel materials. The following non-fuel waste streams were analyzed: Activated metals from decommissioned commercial power reactors, including Type 304 stainless steel, reactor vessel steel, and Inconel, Control blades from boiling water reactors, Neutron-activated corrosion products on surfaces of activated metals, and Spent resins from power plant operations. Measured elemental compositions, including major constituents and impurities, for steel and Inconel samples from commercial power reactors were used in activation calculations to determine radionuclide inventories in activated metals. For a simplified cask model, 60 Co contribution to the total external package dose rate at 30 days after shutdown varied from approximately 60% to 95%, depending on the activated metal, initial cobalt impurity concentration in the metal, and the thickness of the overpack gamma shield. Its maximum contribution to the total external package dose rate of approximately 100% was reached within the time interval of 2 to 5 years after shutdown and was maintained for up to 45 to 60 years after shutdown, depending on material, initial cobalt impurity concentration, and shield thickness. Thereafter, the 60 Co contribution to external package dose rate decreased with increasing decay time. Cobalt-60 is primarily produced by neutron reactions with the cobalt impurity in steel and Inconel. Other important radionuclides in activated metals contributing to package external dose rate are radionuclides with relatively short decay times, including 51 Cr, 59 Fe, 58 Co, and 54 Mn. These radionuclides may be represented as an equivalent 60 Co activity/source because 60 Co gamma ray emissions are bounding in terms of source strength and energy to other important radionuclides identified in the analyzed activated metals. Approaches for modeling the neutron-activated corrosion products that may be attached to activated reactor components were analyzed in this report. It was demonstrated that a surface source is more conservative than a uniform volumetric source for the treatment of neutron-activated corrosion products with respect to external package dose rates. An analysis of the maximum radionuclide loadings reported on spent resins identified the radionuclides 137 Cs, 60 Co, 134 Cs, 65 Zn, and 58 Co as the primary contributors to external package dose rate. For a resin cooled for 3.08 years, the external package dose rate was entirely produced by the reported 137 Cs and 60 Co inventory. The neutron sources from actinides found on spent resins or activated metals produced negligible dose rates and may be ignored in dose rate analyses. Effects of idealized waste material, source geometry, and spatial material/source distributions on external package dose rates were determined based on dose rate results for a simplified cask model under normal conditions of transport. Type 304 stainless steel, zirconium, and aluminum with adjusted mass densities based on a maximum content weight were analyzed for material modeling effects on external package dose rate. These materials produced identical external package dose rates within the statistical uncertainties of the dose rate estimates. Among four different source geometry configurations with homogeneous material of different mass densities, uniform volumetric source distribution, and the same source strength, the geometry configuration with lowest mass density (i.e., minimum self-shielding effect among the four cases) was most conservative. Spatial source distributions that better represent localized peak 60 Co activity values were more conservative than a uniform volumetric source distribution, assuming the same weight and total source strength per package. The increase in external dose rate caused by localized source peaks can be as much as the ratio between source peak activity density to average activity density, depending on the location of the activated metal with peak activity density. Therefore, the shielding analysis may be simplified if localized peak activities and the average activity per package can be measured/determined and documented at the time of cask loading. For simplicity, external package dose rates may be determined based on average source activity and a uniform volumetric source distribution. The dose rate results from that calculation model multiplied by the ratio of peak activity density to average activity density will produce maximum dose rate values for conservative estimates.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Solitary zonal structures in subcritical drift waves: a minimum model

Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same 'equation of state', which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Two-Dimensional Molecular Charge Density Waves in Single-Layer-Thick Islands of a Dirac Fermion System

Charge density waves have been intensely studied in inorganic materials such as transition metal dichalcogenides however its counterpart in organic materials has yet to be explored in detail. Here we report the finding of a robust two dimensional charge density waves in molecular layers formed by α-(BEDT-TTF) 2 -I 3 on a Ag(111) surface. Low temperature scanning tunnelling microscopy images of a multi-layer thick α-(BEDT-TTF) 2 -I 3 on Ag(111) substrate reveal coexistence of 5a0 x 5a0 and √31 a0 x √31 a0 R9° charge density wave patterns commensurate with the underlying molecular lattice at 80 K. Both charge density wave patterns remain in nano-size molecular islands with just a single constituent molecular-layer thickness at 80 K and 5 K. Local tunneling spectroscopy measurements reveal the variation of the gap from 244 meV to 288 meV between the maximum and minimum charge density wave locations. Density functional theory calculations further confirm a vertical positioning of BEDT-TTF molecules in the molecular layer. While the observed charge density wave patterns are stable for the defect sites, they can be reversely switched for one molecular lattice site by means of inelastic tunnelling electron energy transfer with the electron energies exceeding 400 meV using a scanning tunneling microscope manipulation scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The High-resolution Urban Meteorology for Impacts Dataset (HUMID) daily for the Conterminous United States

Many current gridded surface meteorological datasets are inadequate for quantifying near surface spatiotemporal variability because they do not fully represent the impacts of land surface heterogeneity. Of note, explicit representation of the spatial structure and magnitude of local urban warming are usually lacking. Here we enhance the representation of spatial meteorological variability over urban areas in the conterminous United States (CONUS) by employing the High-Resolution Land Data Assimilation System (HRLDAS), which accounts for the fine-scale impacts of spatiotemporally varying land surfaces on weather. We also synthesize in situ meteorological data including local mesonets to create a 1 km grid spacing model-observation fusion product spanning 1981-2018 over the CONUS. Daily maximum, minimum, and mean values for a variety of temperature estimates, humidity, and surface energy budget terms, among others, are included. This High-resolution Urban Meteorology for Impacts Dataset (HUMID) will be useful for studies examining spatial variability of near surface meteorology and the impacts of urban heat islands across many disciplines including epidemiology, ecology, and climatology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The HalfDome multi-survey cosmological simulations: N-body simulations

Upcoming cosmological surveys have the potential to reach groundbreaking discoveries on multiple fronts, including the neutrino mass, dark energy, and inflation. Most of the key science goals require the joint analysis of datasets from multiple surveys to break parameter degeneracies and calibrate systematics. To realize such analyses, a large set of mock simulations that realistically model correlated observables is required. In this paper we present the N-body component of the HalfDome cosmological simulations, designed for the joint analysis of Stage-IV cosmological surveys, such as Rubin LSST, Euclid, SPHEREx, Roman, DESI, PFS, Simons Observatory, CMB-S4, and LiteBIRD. Our 300TB initial data release includes full-sky lightcones and halo catalogs between z = 0–4 for 11 fixed cosmology realizations, as well as an additional run with local primordial non-Gaussianity (f NL = 20). The simulations evolve 61443 particles in a 3.75 h -1 Gpc box, reaching a minimum halo mass of ∼6 × 1012 h -1 M ⊙ and maximum scale of k ∼ h Mpc-1. Our data is publicly available: instructions to access the data and plans for future data releases can be found at https://halfdomesims.github.io.

Bayer, Adrian E↗
DOE OSTI2024

Circumventing superexponential runtimes for hard instances of quantum adiabatic optimization

Classical optimization problems can be solved by adiabatically preparing the ground state of a quantum Hamiltonian that encodes the problem. The performance of this approach is determined by the smallest gap encountered during the evolution. Here, we consider the maximum independent set problem, which can be efficiently encoded in the Hamiltonian describing a Rydberg atom array. We present a general construction of instances of the problem for which the minimum gap decays superexponentially with system size, implying a superexponentially large time to solution via adiabatic evolution. The small gap arises from locally independent choices which cause the system to initially evolve and localize into a configuration far from the solution in terms of Hamming distance. We investigate remedies to this problem. Specifically, we show that quantum quenches in these models can exhibit signatures of quantum many-body scars, which in turn, can circumvent the superexponential gaps. By quenching from a suboptimal configuration, states with a larger ground-state overlap can be prepared, illustrating the utility of quantum quenches as an algorithmic tool. Published by the American Physical Society 2024

Schiffer, Benjamin F. (ORCID:0000000189512157)↗
DOE OSTI2021

Onset of selective laser flash sintering of AlN

Abstract Flash sintering uses a combination of heating and electric fields to rapidly densify ceramics. Previously, it has been shown that a scanning laser can be used to initiate flash sintering in localized regions on an yttria‐stabilized zirconia (YSZ) sample in a process known as selective laser flash sintering (SLFS). In this work, we show using a combination of measurements of electric current flowing through the sample and observations of necks formed between powder particles that aluminum nitride (AlN) can also undergo SLFS. Scan conditions required to initiate SLFS are characterized over a range of laser powers and laser scan speeds in a dry nitrogen environment. It is shown that initiation of SLFS in AlN is governed by both the local input energy density per scan and heat dissipation and a numerical model is developed to predict temperatures during SLFS. Assuming the minimum temperature along the conductive path determines the onset of SLFS, the minimum temperature and time required is 450–670 K in 2–0.25 s for the pressed AlN pellets used in this study for laser scan speeds of 33–300 m/s, laser powers of 10–30 W, and an applied electric field of 3000 V/cm.

Gouws, Andre↗
DOE OSTI2021

Techno-Economic Assessment for Opportunities to Integrate Algae Farming with Wastewater Treatment

This report investigates the economic and scalability potential for algal biomass production based on opportunities to couple algae farming with wastewater treatment (WWT). This activity was motivated in part by a recent influx of technology developers who are currently in various stages of planning, construction, and operation of algae facilities for treating municipal wastewater, all making claims of tremendous economic potential for algal WWT relative to other more traditional WWT methods when viewing treated wastewater as the primary product output and algal biomass as a secondary coproduct. In addition to solving for minimum biomass selling price (MBSP) for a farm configured to maximize biomass production while minimizing water footprint (when sourcing makeup water from local groundwater resources at added cost), the algal WWT cases evaluated here also solved for MBSP while assigning a treatment credit as a means to valorize the treated water as a coproduct in the farm model. This enables a consistent comparison against the National Renewable Energy Laboratory’s standard algae farm models.

09 BIOMASS FUELS↗
DOE OSTI2021

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2025

Cultivating clarity: understanding the impact of land cost assumptions on biofuel viability

The transition to sustainable energy has increased interest in biofuel production to reduce greenhouse gas emissions, decrease reliance on imported oil, and ensure energy resilience. Here, this study examines the often-overlooked impact of land cost assumptions on the economic viability of biofuel production. Using a discounted cash flow techno-economic framework, we evaluated three land cost scenarios—no land costs, land rental costs, and land purchase costs—across six bioenergy feedstocks (corn, soybeans, switchgrass, miscanthus, poplar, and microalgae) and three biofuel products (corn ethanol, soybean biodiesel, and sustainable aviation fuel) at the county level for the contiguous United States. The analysis reveals substantial variation in minimum fuel selling prices due to these scenarios. High-yield crops like algae showed low sensitivity to land costs, while low-yield crops such as soybeans were highly sensitive. Geographical differences were significant, with minimum fuel selling price increases most pronounced in high-value land regions like the Corn Belt. Case studies further illustrate the influence of local productivity and land costs on economic outcomes across the United States. These findings emphasize the importance of maintaining consistent land cost assumptions in biofuel economic assessments. By quantifying the interplay between land value, crop productivity, and economic feasibility, this study provides essential insights for policymakers and stakeholders to advance sustainable energy solutions.

09 BIOMASS FUELS↗
DOE OSTI2023

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗
DOE OSTI2024

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Steady-State Thermal-Hydraulic Analysis of the LEU-Fueled Dalat Nuclear Research Reactor

This paper presents results of steady-state thermal-hydraulic analysis for the designed working core of the Dalat Nuclear Research Reactor (DNRR) using the PLTEMP/ANL code. The core was designed to be loaded with 92 low-enriched uranium (LEU) VVR-M2 fuel bundles (FBs) and 12 beryllium rods surrounding a neutron trap at the core center, for replacement of the previous core with 104 high-enriched uranium (HEU) VVR-M2 FBs. Before using this code for thermohydraulic analysis of the designed LEU working core, it was validated by comparing calculation results with experimental data collected from the HEU working core of the DNRR. The discrepancy between calculated results and measured data was at the maximum about 0.8°C and 1.5°C of fuel cladding and outlet coolant temperatures, respectively. In the design calculation, thermohydraulic safety was confirmed through evaluation of the fuel cladding and coolant temperatures, as well as of other safety parameters such as Departure from Nucleate Boiling Ratio (DNBR) and Onset of Nucleate Boiling Ratio (ONBR). The calculation results showed that, in normal operation conditions at full nominal thermal power of 500 kW without uncertainty parameters, the maximum fuel cladding temperature of the hottest FB was about 90.4°C, which is lower than its limit value of 103°C, the minimum DNBR was 32.0, which is much higher than the recommended value of 1.5, and the minimum ONBR was 1.43, which is higher than the recommended value of 1.4 for VVR-M2 LEU fuel type. When the global and local hot channel factors were taken into account, the maximum temperature of fuel cladding at the hottest FB was about 98.4 °C, for global only, and 114.3°C, for global together with local hot channel factors. The calculation results confirm the safety operation of the designed LEU core loaded with 92 fresh VVR-M2 FBs.

Nguyen, Kien-Cuong↗
DOE OSTI2024

Laterally Modulating Carrier Concentration by Ion Irradiation in CdO Thin Films for Mid‐IR Plasmonics

This report demonstrates tunable carrier densities in CdO thin films through local ion irradiation, providing lateral control of mid‐IR optical properties. Ion‐solid interactions produce donor‐like defects that boost electron concentrations from the practical minimum of 2.5 × 10 19 cm −3 to a maximum of 2.5 × 10 20 cm −3 by metered ion exposure. This range is achieved using He, N, Ar, or Au ions at 1–2.8 MeV; when normalized by displacements per atom, all ion species produce comparable results. Since CdO is well‐described by the Drude model, irradiation‐tuned carrier densities directly alter the infrared dielectric function, and in turn, mid‐infrared optical properties. Further, it is demonstrated that by combining irradiation with traditional lithography, CdO films expose to ions in the presence of 3‐µm thick, patterned photoresist exhibit lateral carrier density profiles with ≈400‐nm resolution. Scanning near‐field optical microscopy reveals sharp optical interfaces with almost no companion contrast in surface morphology, microstructure, or crystallinity. Finally, CdO lateral homostructures supporting surface plasmon polaritons (SPPs) are demonstrated whose dispersion relation can be tuned through periodic patterning in a monolithic platform by simple nanofabrication. Numerical simulations show these polaritons result from strong coupling between excitations at CdO plasma frequencies and SPPs supported by the platinum substrate.

cadmium oxide↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2025

Constraints on light QCD and CP-violating axions from the death line of rotation-powered pulsars

For axions that couple to nucleons, the presence of dense nuclear matter can displace the axion from its vacuum minimum, sourcing large field gradients around neutron stars (and, more generally, compact objects). These gradients, which we refer to as axion hair, couple to the local background magnetic field, inducing a large voltage drop near the surface of the star; here, we demonstrate that the presence of axion hair decouples local near-field particle acceleration in the open magnetic field line bundle from the rotational frequency of the pulsar itself. This is significant as the non-observation of old slowly-rotating pulsars is attributed to the fact the rotationally-induced electric fields are not strong enough to sustain $e^\pm$ pair production. In this work, we review the evidence for the existence for `pulsar death', i.e. the threshold at which $e^\pm$ pair production (and thus, by association, coherent radio emission) ceases, and demonstrate using both semi-analytics and particle-in-cell simulations that the existence of axion hair can dramatically extend pulsar lifetimes. We show that the non-observation of extremely old, slowly rotating, pulsars allows for a new probe of light QCD and CP-violating axions. We also demonstrate how the observation of emission from both poles of pulsars with nearly orthogonal rotational and magnetic axes, as seen e.g. in PSR J1906+0746, can be used to set competitive limits on CP-violating axion-nucleon interactions.

Witte, Samuel J. [Oxford U., Theor. Phys.; DESY; H↗
DOE OSTI2021

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗
DOE OSTI2024

The Collaborative Seismic Earth Model: Generation 2

Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.

58 GEOSCIENCES↗