Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “GEOMETRIC FACTOR”

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 OSTI2024

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Aims.In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods.We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳ ap 3 〉 in a tomographic set-up. We validated our pipeline withN-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFITfor the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳ ap 3 〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results.We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩ m < 0.02 σ Ω m , and of the structure growth parameter, ΔS 8 < 0.05 σ S 8 . Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ω m = 0.248 −0.055 +0.062 andS 8 = σ 8 √(Ω m /0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ω m andS 8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.

Astronomy & Astrophysics↗
DOE OSTI2023

Thermodynamic characterization of synthetic lead-arsenate apatites with different halogen substitutions

Abstract Thermodynamic parameters have been measured for synthetic analogs of the mimetite-group minerals Pb5(AsO4)3X (X = OH, Cl, Br, I) belonging to the apatite supergroup. Phases precipitated from aqueous solutions under ambient conditions with well characterized structures and compositions were studied. For each phase, dissolution enthalpy was experimentally determined by oxide melt drop solution calorimetry in a molten solvent of sodium molybdate (3Na2O·4MoO3) at 976 K. The enthalpy of formation from the elements ΔHf,elo was calculated using thermochemical cycles and was −3030.6 ± 11.5, −3026.6 ± 15.8, −2967.6 ± 25.0, and −2993.1 ± 12.2 kJ/mol for Pb5.00(AsO4)3.00OH0.86(CO3)0.07, Pb5.00(AsO4)3.00Cl0.80(CO3)0.10, Pb5.00(AsO4)3.00Br0.80(CO3)0.10, and Pb5.00(AsO4)3.00I0.45OH0.35(CO3)0.10, respectively. These ΔHf,elo values exhibit typical trends for apatites: they increased (were less negative) with the increasing molar mass and ionic radius of X and decreased with the electronegativity and ionization energy of X. The compilation and comparison of data for Ca-, Pb-, P-, and As-apatites revealed correlations indicating that thermodynamic enthalpic stability is largely influenced by chemical factors (e.g., differences in electronegativities of the elements, ionization energy, or ionic characteristics of the bonds) and to a lesser extent by physical and geometric parameters in the crystal structure related to the mass and size of the X anion. Using the correlations, it was possible to estimate the value of hitherto unknown ΔHf,elo for Pb5(AsO4)3F, −3144.3 ± 66.5 kJ/mol. The observed relationships apply to the entire apatite supergroup and can be used to predict the values of ΔHf,elo for phases that have not been studied experimentally. The new data on environmentally significant phases will contribute to the modeling of mineral-water interactions, particularly for potential use in the remediation of soils and wastes contaminated with Pb and As and in the immobilization of radioactive waste containing I-129.

Geochemistry & Geophysics↗
DOE OSTI2021

Twisted bilayer graphene. V. Exact analytic many-body excitations in Coulomb Hamiltonians: Charge gap, Goldstone modes, and absence of Cooper pairing

We find exact analytic expressions for the energies and wave functions of the charged and neutral excitations above the exact ground states (at rational filling per unit cell) of projected Coulomb Hamiltonians in twisted bilayer graphene. Our exact expressions are valid for any form of the Coulomb interaction and any form of A A and A B / B A tunneling. The single charge excitation energy is a convolution of the Coulomb potential with a quantum geometric tensor of the TBG bands. The neutral excitations are (high-symmetry group) magnons, and their dispersion is analytically calculated in terms of the form factors of the active bands in TBG. The two-charge excitation energy and wave functions are also obtained, and a sufficient condition on the graphene eigenstates for obtaining a Cooper pair from Coulomb interactions is obtained. For the actual TBG bands at the first magic angle, we can analytically show that the Cooper pair binding energy is zero in all such projected Coulomb models, implying that either phonons and/or nonzero kinetic energy are needed for superconductivity. Since Vafek and Kang [Phys. Rev. Lett. 125, 257602 (2020)] showed that the kinetic energy bounds on the superexchange energy are less 10 -3 in Coulomb units, the phonon mechanism becomes then very likely. If nonetheless the superconductivity is due to kinetic terms which render the bands nonflat, one prediction of our theory is that the highest T c would not occur at the highest DOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.4 is a program that obtains a steady-state flow and temperature solution, in geometrical extent, for a nuclear reactor core, or a single fuel assembly, including the effect of manufacturing and modelling uncertainties by the means of hot channel factors. It is based on an evolutionary sequence of codes originally used for calculating plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over the past 30 years. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability for select design and safety parameters is available to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Beyond the acceptance limit of DRAGON: The case of the 6 Li(α,γ) 10 B reaction

Radiative capture reactions play a pivotal role for our understanding of the origin of the elements in the cosmos. Recoil separators provide an effective way to study these reactions, in inverse kinematics, and take advantage of the use of radioactive ion beams. However, a limiting factor in the study of radiative capture reactions in inverse kinematics is the momentum spread of the product nuclei, which can result in an angular spread larger than the geometric acceptance of the separator. The DRAGON facility at TRIUMF is a versatile recoil separator, designed to study radiative capture reactions relevant to astrophysics in the A~10–30 region. In this work we present the first attempt to study with DRAGON a reaction, 6 Li(α, γ) 10 B, for which the recoil angular spread exceeds DRAGON’s acceptance. Furthermore, our result is in good agreement with the literature value, showing that DRAGON can measure resonance strengths of astrophysically important reactions even when not all the recoils enter the separator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spin dynamics of the antiferromagnetic Heisenberg model on a kagome bilayer

We research the spin dynamics of a classical Heisenberg antiferromagnet with nearest-neighbor interactions on a quasi-two-dimensional kagome bilayer. This geometrically frustrated lattice consists of two kagome layers connected by a triangular-lattice linking layer. By combining Monte Carlo with precessional spin dynamics simulations, we compute the dynamical structure factor of the classical spin liquid in the kagome bilayer and investigate the thermal and dilution effects. While the low-frequency and long-wavelength dynamics of the cooperative paramagnetic phase is dominated by spin diffusion, weak magnon excitations persist at higher energies, giving rise to the half-moon pattern in the dynamical structure factor. In the presence of spin vacancies, the dynamical properties of the diluted system can be understood within the two-population picture. The spin diffusion of the “correlated” spin clusters is mainly driven by the zero-energy weather-vane modes, giving rise to an autocorrelation function that decays exponentially with time. On the other hand, the diffusive dynamics of the quasifree “orphan” spins leads to a distinctive longer-time power-law tail in the autocorrelation function. We discuss the implications of our work for the glassy behaviors observed in the archetypal frustrated magnet SrCr 9p Ga 12-9 pO 19 (SCGO).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Modeling Nanoscale Ohmics in Carbon Supports of Fuel Cell Cathodes

Here, reducing platinum (Pt) loading in polymer electrolyte fuel cells (PEFCs) while meeting performace requirements is critical to their widespread deployment. However, significant polarization losses manifest at higher current densities in cathodes with lower Pt content. The morphology of the carbon supports in PEFC cathodes affects the location of Pt deposition into the micropits or onto the surface of the carbon support, translating into different kinetic and transport resistances. In this work, we derive an agglomerate scale model that differentiates the sink terms for Pt on the surface and in the pits of carbon supports. We develop an approach to assess the impact of nanoscale ohmic resistance to Pt in the micropits arising from weakly ionic solution in the carbon support on PEFC performance. Effectiveness factors relating the actual reaction rate to the maximum reaction rate (had all the Pt been exposed) are derived and embedded into a one-dimensional catalyst layer model. Parameters in the catalyst layer model are tuned based on experimental local oxygen transport resistances. Subsequently, we estimate bounds for the micropore resistances based on geometric and physical arguments. Lastly, polarization curves are simulated to assess the effect of the micropore resistance in fully-humidified and oxygen-rich environments.

25 ENERGY STORAGE↗
DOE OSTI2023

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding

Accurate protein synthesis is determined by the two-subunit ribosome’s capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site’s physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

AI-Assisted Conceptual Development of a Pre-Geometric Cosmological Model - An Exercise in AI-Assisted Conceptual Framework Generation, Paper II: Local Geometry and Metric Structure

This paper develops the geometric sector of the replication-driven cosmogenesis framework introduced in Paper I. Starting from a pre-geometric spectral substrate and a minimal set of replication axioms, we show how coherent self-replicating units generate a spatial adjacency graph whose continuum limit acquires an effective Riemannian structure. The replication dynamics determines a characteristic correlation length that seeds the local metric, while overlap relations among coherent units produce an isotropic neighborhood geometry with an emergent dimensionality $d_{\rm eff}\simeq 3$ across a broad range of replication factors. As replication slows and causal order stabilizes, a limiting signal speed $c_\ast$ appears, providing the basis for the Lorentzian structure of spacetime without assuming a pre-existing light cone. We derive conditions under which the adjacency graph converges to a smooth three-dimensional manifold, describe the transition from Euclidean to Lorentzian propagation, and identify geometric invariants controlled by the replication parameters. This work establishes the geometric and causal layer of the replication cosmogenesis program, bridging the spectral axioms of Paper I to the cosmological dynamics explored in Paper III.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2023

Measurement of the Z boson production cross-section in proton-lead collisions at $\sqrt{s_{\textrm{NN}}}$ = 8.16 TeV

This article presents the first measurement of the differential Z-boson production cross-section in the forward region using proton-lead collisions with the LHCb detector. The dataset was collected at a nucleon-nucleon centre-of-mass energy of $\sqrt{s_{NN}}$ = 8.16 TeV in 2016, corresponding to an integrated luminosity of 30.8 nb –1 . The forward-backward ratio and the nuclear modification factors are measured together with the differential cross-section as functions of the Z boson rapidity in the centre-of-mass frame, the transverse momentum of the Z boson and a geometric variable $\phi$*. The results are in good agreement with the predictions from nuclear parton distribution functions, providing strong constraining power at small Bjorken-x.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Investigating the Relationship Between Bolide Entry Angle and Apparent Direction of Infrasound Signal Arrivals

Infrasound sensing offers critical capabilities for detecting and geolocating bolide events globally. However, the observed back azimuths, directions from which infrasound signals arrive at stations, often differ from the theoretical expectations based on the bolide’s peak brightness location. For objects with shallow entry angles, which traverse longer atmospheric paths, acoustic energy may be emitted from multiple points along the trajectory, leading to substantial variability in back azimuth residuals. This study investigates how the entry angle of energetic bolides affects the back azimuth deviations, independent of extrinsic factors such as atmospheric propagation, station noise, and signal processing methodologies. A theoretical framework, the Bolide Infrasound Back-Azimuth EXplorer Model (BIBEX-M), was developed to compute predicted back azimuths solely from geometric considerations. The model quantifies how these residuals vary as a function of source-to-receiver distance, revealing that bolides entering at shallow angles, e.g., 10°, can produce average residuals of 20°, with deviations reaching up to 46° at distances below 1000 km, and remaining significant even at 5000 km (up to 8°). In contrast, bolides with steeper entry angles, e.g., > 60°, show smaller deviations, typically under 5° at 1000 km and diminishing to less than ~1° beyond 5000 km. These findings attest to the need for careful interpretation when evaluating signal detections and estimating bolide locations. This work is not only pertinent to bolides but also to other high-energy, extended-duration atmospheric phenomena such as space debris and reentry events, where similar geometric considerations can influence infrasound arrival directions.

Acoustics↗
DOE OSTI2024

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol

Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.

Aqueous chemistry↗
DOE OSTI2026

Linear complexity

We present factorization and solution phases for a new linear complexity direct solver designed for concurrent batch operations on fine-grained parallel architectures, for matrices amenable to hierarchical representation. We focus on the strong-admissibility-based $\mathscr{H}^{2}$ format, where strong recursive skeletonization factorization compresses remote interactions. We build upon previous implementations of $\mathscr{H}^{2}$ matrix construction for efficient factorization and solution algorithm design, which are illustrated graphically in stepwise detail. The algorithms are ‘blackbox’ in the sense that the only inputs are the matrix and right-hand side, without analytical or geometrical information about the origin of the system. We demonstrate linear complexity scaling in both time and memory on four representative families of dense matrices up to one million in size. Parallel scaling up to 16 threads is enabled by a multi-level matrix graph coloring and avoidance of dynamic memory allocations thanks to prefix-sum memory management. An experimental backward error analysis is included. We break down the timings of different phases, identify phases that are memory-bandwidth limited, and discuss alternatives for phases that may be sensitive to the trend to employ lower precisions for performance.

Boukaram, Wajih↗
DOE OSTI2021

Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes

Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

Valdivia, Pedro↗
DOE OSTI2020

Precautions of Using Three-Electrode Configuration to Measure Electrode Overpotential in Solid Oxide Electrochemical Cells: Insights from Finite Element Modeling

Accurate determination of electrode overpotentials is essential to assess the performance of the electrode and understand the rate-limiting steps involved. Three-electrode configuration with the use of a reference electrode at a fixed potential is a standard way to measure overpotential of a specific electrode in liquid electrochemical systems. However, application of such three-electrode configuration to solid electrochemical cells for overpotential determination is not straightforward and requires extra caution. Here we report a theoretical Finite Element Analysis on the geometrical requirements for which the reference electrode can be applied to anode- or electrolyte-supported solid oxide button cells. The modeling results suggest that the symmetry of the working and counter electrodes is the key factor determining if a reference electrode is suitable to use. For anode-supported fuel cells with asymmetrical working and counter electrodes, reference electrode cannot be used under all circumstances. To use reference electrode for overpotential measurements, electrolyte-supported cells with symmetrical semicircular-shaped electrodes are preferred. A data processing procedure has also been presented to obtain the electrode overpotential from the measured potential using the three-electrode scheme in solid oxide electrochemical cells.

25 ENERGY STORAGE↗
DOE OSTI2024

Upsampling Monte Carlo Reactor Simulation Tallies in Depleted Sodium-Cooled Fast Reactor Assemblies Using a Convolutional Neural Network

The computational demand of neutron Monte Carlo transport simulations can increase rapidly with the spatial and energy resolution of tallied physical quantities. Convolutional neural networks have been used to increase the resolution of Monte Carlo simulations of light water reactor assemblies while preserving accuracy with negligible additional computational cost. Here, we show that a convolutional neural network can also be used to upsample tally results from Monte Carlo simulations of sodium-cooled fast reactor assemblies, thereby extending the applicability beyond thermal systems. The convolutional neural network model is trained using neutron flux tallies from 300 procedurally generated nuclear reactor assemblies simulated using OpenMC. Validation and test datasets included 16 simulations of procedurally generated assemblies, and a realistic simulation of a European sodium-cooled fast reactor assembly was included in the test dataset. We show the residuals between the high-resolution flux tallies predicted by the neural network and high-resolution Monte Carlo tallies on relative and absolute bases. The network can upsample tallies from simulations of fast reactor assemblies with diverse and heterogeneous materials and geometries by a factor of two in each spatial and energy dimension. The network’s predictions are within the statistical uncertainty of the Monte Carlo tallies in almost all cases. This includes test assemblies for which burnup values and geometric parameters were well outside the ranges of those in assemblies used to train the network.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗