Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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

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

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

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

Time dependent chlorination of CeO 2 , La 2 O 3 and Nd 2 O 3 by ZrCl 4 dissolved in eutectic LiCl–KCl

Partial chlorination of CeO 2 , La 2 O 3 , and Nd 2 O 3 was achieved using ZrCl 4 dissolved in molten LiCl-KCl at 450°C over a period of 4-5 hr. ZrCl 4 was pre-volatilized into an argon stream to maintain a constant vapor pressure for equilibration with the salt. This approach was effective at maintaining the ZrCl 4 concentration at the solubility limit of 1.3 mol% in the molten salt throughout the duration of two of the chlorination experiments and within about 60% of the solubility limit in the third experiment. Based on the analysis of salt samples, fractional conversion after four hours of reaction for CeO 2 , and five hours for La 2 O 3 , and Nd 2 O 3 were calculated to be 0.97, 0.28, and 0.65 respectively. Assuming a shrinking core diffusion-limited model, the diffusivity of ZrCl 4 through the hypothesized ZrO 2 layer was calculated for the chlorination of CeO 2 , La 2 O 3 , and Nd 2 O 3 were 1.76 x 10 -10 cm 2 /sec, 7.85 x 10 -12 cm 2 /sec, and 7.42 x 10 -12 cm 2 /sec, respectively. Given the extremely low values calculated for diffusivity, it was concluded that the process is chemical reaction rate limited and that a complete reaction of each rare earth oxide could be achieved. Furthermore, a chemically controlled shrinking core model was shown to also fit the experimental data closely.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Energy Spectra of Abundant Cosmic-ray Nuclei in Sources, According to the ATIC Experiment

One of the main results of the ATIC (Advanced Thin Ionization Calorimeter) experiment is a collection of energy spectra of abundant cosmic-ray nuclei: protons, He, C, O, Ne, Mg, Si, and Fe measured in terms of energy per particle in the energy range from 50 GeV to tens of teraelectronvolts. In this paper, the ATIC energy spectra of abundant primary nuclei are back-propagated to the spectra in sources in terms of magnetic rigidity using a leaky-box approximation of three different GALPROP-based diffusion models of propagation that fit the latest B/C data of the AMS-02 experiment. It is shown that the results of a comparison of the slopes of the spectra in sources are weakly model dependent; therefore the differences of spectral indices are reliable data. A regular growth of the steepness of spectra in sources in the range of magnetic rigidity of 50–1350 GV is found for a charge range from helium to iron. This conclusion is statistically reliable with significance better than 3.2 standard deviations. The results are discussed and compared to the data of other modern experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TH-AB-209-03: Overcoming Resolution Limitations of Diffuse Optical Signals in X-Ray Induced Luminescence (XIL) Imaging Via Selective Plane Illumination and 2D Deconvolution

Purpose: To evaluate the resolution and sensitivity of XIL imaging using a surface radiance simulation based on optical diffusion and maximum likelihood expectation maximization (MLEM) image reconstruction. XIL imaging seeks to determine the distribution of luminescent nanophosphors, which could be used as nanodosimeters or radiosensitizers. Methods: The XIL simulation generated a homogeneous slab with optical properties similar to tissue. X-ray activated nanophosphors were placed at 1.0 cm depth in the tissue in concentrations of 10{sup −4} g/mL in two volumes of 10 mm{sup 3} with varying separations between each other. An analytical optical diffusion model determined the surface radiance from the photon distributions generated at depth in the tissue by the nanophosphors. The simulation then determined the detected luminescent signal collected with a f/1.0 aperture lens and back-illuminated EMCCD camera. The surface radiance was deconvolved using a MLEM algorithm to estimate the nanophosphors distribution and the resolution. To account for both Poisson and Gaussian noise, a shifted Poisson imaging model was used in the deconvolution. The deconvolved distributions were fitted to a Gaussian after radial averaging to measure the full width at half maximum (FWHM) and the peak to peak distance between distributions was measured to determine the resolving power. Results: Simulated surface radiances for doses from 1mGy to 100 cGy were computed. Each image was deconvolved using 1000 iterations. At 1mGy, deconvolution reduced the FWHM of the nanophosphors distribution by 65% and had a resolving power is 3.84 mm. Decreasing the dose from 100 cGy to 1 mGy increased the FWHM by 22% but allowed for a dose reduction of a factor of 1000. Conclusion: Deconvolving the detected surface radiance allows for dose reduction while maintaining the resolution of the nanophosphors. It proves to be a useful technique in overcoming the resolution limitations of diffuse optical imaging in tissue. C. S. acknowledges support from the NIH National Institute of General Medical Sciences (Award number R25GM109439, Project Title: University of Chicago Initiative for Maximizing Student Development, IMSD). B. Q. and P. L. acknowledge support from NIH grant R01EB017293.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

On how structures convey non-diffusive turbulence spreading

Abstract We report on comprehensive experimental studies of turbulence spreading in edge plasmas. These studies demonstrate the relation of turbulence spreading and entrainment to intermittent convective density fluctuation events or bursts (i.e. blobs and holes). The non-diffusive character of turbulence spreading is thus elucidated. The turbulence spreading velocity (or mean jet velocity) manifests a linear correlation with the skewness of density fluctuations, and increases with the auto-correlation time of density fluctuations. Turbulence spreading by positive density fluctuations is outward, while spreading by negative density fluctuations is inward. The degree of symmetry breaking between outward propagating blobs and inward propagating holes increases with the amplitude of density fluctuations. Thus, blob-hole asymmetry emerges as crucial to turbulence spreading. These results highlight the important role of intermittent convective events in conveying the spreading of turbulence, and constitute a fundamental challenge to existing diffusive models of spreading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Capacitive radio frequency discharges with a single ring-shaped narrow trench of various depths to enhance the plasma density and lateral uniformity

Spatial structures of the electron density and temperature in ring-shaped hollow cathode capacitive rf plasma with a single narrow trench of 2 mm width have been investigated at various trench depths of D = 5, 8, 10, 12, and 15 mm. It is found that the plasma density is increased in the presence of the trench and that the radial profile of the plasma density has a peak around the narrow hollow trench near the cathode. The density becomes uniform further away from the cathode at all trench depths, whereas the electron temperature distribution remains almost uniform. The measured radial profiles of the plasma density are in good agreement with a theoretical diffusion model for all the trench depths, which explains the local density increase by a local enhancement of the electron heating. Under the conditions investigated, the trench of 10 mm depth is found to result in the highest plasma density at various axial and radial positions. The results show that the radial uniformity of the plasma density at various axial positions can be improved by using structured electrodes of distinct depths rather than planar electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evolving radon diffusion through earthen barriers at uranium waste disposal sites

Field measurements of Rn-222 fluxes from the tops and bottoms of compacted clay radon barriers were used to calculate effective Rn diffusion coefficients (D Rn ) at four uranium waste disposal sites in the western United States to assess cover performance after more than 20 years of service. Values of D Rn ranged from 7.4 × 10 -7 to 6.0 × 10 -9 m 2 /s, averaging 1.42 × 10 -7 . Water saturation (S W ) from soil cores indicated that there was relatively little control of D Rn by S W , especially at higher moisture levels, in contrast to estimates from most steady-state diffusion models. Further, this is attributed to preferential pathways intrinsic to construction of the barriers or to natural process that have developed over time including desiccation cracks, root channels, and insect burrows in the engineered earthen barriers. A modification to some models in which fast and slow pathway D Rn values are partitioned appears to give a good representation of the data; 4% of the fast pathway was needed to fit the data regression. For locations with high S w and highest D Rn (and fluxes) at each site, the proportion of fast pathway ranged from 1.7% to 34%, but for many locations with lower fluxes, little if any fast pathway was needed.

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

Measurement and simulation of top- and bottom-illuminated solar-blind AlGaN metal-semiconductor-metal photodetectors with high external quantum efficiencies

A comprehensive study on top- and bottom-illuminated Al{sub 0.5}Ga{sub 0.5}N/AlN metal-semiconductor-metal (MSM) photodetectors having different AlGaN absorber layer thickness is presented. The measured external quantum efficiency (EQE) shows pronounced threshold and saturation behavior as a function of applied bias voltage up to 50 V reaching about 50% for 0.1 μm and 67% for 0.5 μm thick absorber layers under bottom illumination. All experimental findings are in very good accordance with two-dimensional drift-diffusion modeling results. By taking into account macroscopic polarization effects in the hexagonal metal-polar +c-plane AlGaN/AlN heterostructures, new insights into the general device functionality of AlGaN-based MSM photodetectors are obtained. The observed threshold/saturation behavior is caused by a bias-dependent extraction of photoexcited holes from the Al{sub 0.5}Ga{sub 0.5}N/AlN interface. While present under bottom illumination for any AlGaN layer thickness, under top illumination this mechanism influences the EQE-bias characteristics only for thin layers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development of Corrosion Resistant Coatings for Structural Materials for Liquid Fueled Molten Salts Reactors Applications (Final Report)

Current structural alloys, code certified for the MSRs temperature ranges, contain high levels of chromium, making them highly susceptible to molten halide salt corrosion. One potential solution to circumvent the need for code certification of novel alloys, which is expensive and time consuming, is to design claddings that will protect the underlying code certified materials from corrosion damage during operation. In this work, we examined the corrosion of Ni and Cu electroplated, Ni and CuNi weld overlay, Mo-laser clad, and carburized claddings on SS316H for use in molten salt reactor environments. After characterization of the cladded materials (Task 1), static corrosion tests were used to assess corrosion resistance of the claddings in typical molten fluoride salts (Task 2). The corrosion tests were performed in molten FLiNaK at 700°C up to 1000 hours. Pre- and post-corrosion Scanning Electron Microscopy (SEM), energy-dispersive-spectroscopy (EDS), Transmission Electron Microscopy (TEM), X-ray Diffraction (XRD) and glow-discharge-optical-emission-spectroscopy (GDOES) were performed on cladding cross-sections and surfaces to evaluate degradation. The Cu and Ni electroplated samples, as well as the carburized samples, showed excellent corrosion resistance relative to the bare SS316H. To assess high temperature cladding stability, ageing experiments were performed at temperatures up to 900 C in inert atmosphere for the electroplated and carburized samples and a diffusion model was developed to predict long term cladding behavior (Task 3). It was found that the Cu cladding was basically insensitive to the high temperature ageing, except for small secondary phases forming at the interface. On the other hand, the Ni electroplated cladding experienced significant interdiffusion. Nevertheless, the gain in corrosion resistance for a 100μm Ni electroplated cladding is phenomenal, with more than 50% reduction in chromium dissolution from the substrate material for the first 15 years of salt exposure at 700°C. To assess radiation resistance and phase stability, high-temperature 4MeV Ni heavy ion irradiation was performed across the cladding/substrate interface up to 50 displacement-per-atom (DPA) at 500°C and 700°C to assess the phase stability and irradiation behaviors of the electroplated systems (Task 4). The Ni and Cu electroplated systems did not experience void swelling at the contrary to the SS316H substrate due to their nanocrystalline nature. Some level of recrystallization was observed in the cladding, as well as radiation induced segregation and enhanced diffusion. The interface acted as a potent sink for point defects with the presence of a void denuded zone. The mechanical properties of the claddings were assessed using thermal shock, micro-indentation, nano-indentation, and four-point bend testing experiments (Task 5). These experiments were performed on the claddings in as-received, corroded, irradiated, and thermally aged states to determine the effects of typical molten salt reactor environments on cladding integrity. No significant mass loss was observed after repeated thermal shocks. While the electroplated samples softened after high-temperature ageing, irradiation hardening compensate this effect, such that there is little different with the as-received materials. The results of these experiments suggest that the electroplated (copper/nickel) show the greatest promise for application in molten salt reactor development. These claddings prevented any chromium dissolution from occurring during the static corrosion experiments. Additionally, they demonstrated favorable properties for high temperature diffusion, phase stability, interfacial mechanical properties, and irradiation resistance. Weld-overlay cladding are also of great interest and could reach properties similar to the electroplated samples upon optimization (especially with multiple weld passes). The carburized SS316H showed excellent behavior as well, but more work needs to be done to assess their long-term stability. Finally, the Mo-laser clad system was not pursued further due to the manufacturing challenges. It is also worth noting that the Ni cladding systems should behave relatively well in terms of weldability since the Ni-weld overlay microstructure and associated corrosion rate are sound. While the NiCu weld-overlay has even lower corrosion rates than the Ni-weld overlay, more studies on welding of Cu-electroplated systems are necessary since Cu clusters are known to embrittle steels.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Nucleation and structural growth of cluster crystals

We study the nucleation of crystalline cluster phases in the generalized exponential model with exponent n = 4. Due to the finite value of this pair potential for zero separation, at high densities the system forms cluster crystals with multiply occupied lattice sites. Here, we investigate the microscopic mechanisms that lead to the formation of cluster crystals from a supercooled liquid in the low-temperature region of the phase diagram. Using molecular dynamics and umbrella sampling, we calculate the free energy as a function of the size of the largest crystalline nucleus in the system, and compare our results with predictions from classical nucleation theory. Employing bond-order parameters based on a Voronoi tessellation to distinguish different crystal structures, we analyze the average composition of crystalline nuclei. We find that even for conditions where a multiply occupied fcc crystal is the thermodynamically stable phase, the nucleation into bcc cluster crystals is strongly preferred. Furthermore, we study the particle mobility in the supercooled liquid and in the cluster crystal. In the cluster crystal, the motion of individual particles is captured by a simple reaction-diffusion model introduced previously to model the kinetics of hydrogen bonds.

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

Report on the AAPM grand challenge on deep generative modeling for learning medical image statistics

Abstract Background The findings of the 2023 AAPM Grand Challenge on Deep Generative Modeling for Learning Medical Image Statistics are reported in this Special Report. Purpose The goal of this challenge was to promote the development of deep generative models for medical imaging and to emphasize the need for their domain‐relevant assessments via the analysis of relevant image statistics. Methods As part of this Grand Challenge, a common training dataset and an evaluation procedure was developed for benchmarking deep generative models for medical image synthesis. To create the training dataset, an established 3D virtual breast phantom was adapted. The resulting dataset comprised about 108 000 images of size 512 512. For the evaluation of submissions to the Challenge, an ensemble of 10 000 DGM‐generated images from each submission was employed. The evaluation procedure consisted of two stages. In the first stage, a preliminary check for memorization and image quality (via the Fréchet Inception Distance [FID]) was performed. Submissions that passed the first stage were then evaluated for the reproducibility of image statistics corresponding to several feature families including texture, morphology, image moments, fractal statistics, and skeleton statistics. A summary measure in this feature space was employed to rank the submissions. Additional analyses of submissions was performed to assess DGM performance specific to individual feature families, the four classes in the training data, and also to identify various artifacts. Results Fifty‐eight submissions from 12 unique users were received for this Challenge. Out of these 12 submissions, 9 submissions passed the first stage of evaluation and were eligible for ranking. The top‐ranked submission employed a conditional latent diffusion model, whereas the joint runners‐up employed a generative adversarial network, followed by another network for image superresolution. In general, we observed that the overall ranking of the top 9 submissions according to our evaluation method (i) did not match the FID‐based ranking, and (ii) differed with respect to individual feature families. Another important finding from our additional analyses was that different DGMs demonstrated similar kinds of artifacts. Conclusions This Grand Challenge highlighted the need for domain‐specific evaluation to further DGM design as well as deployment. It also demonstrated that the specification of a DGM may differ depending on its intended use.

Radiology, Nuclear Medicine & Medical Imaging↗
NASA NTRSDate not supplied

Edge-Based Viscous Method for Mixed-Element Node-Centered Finite-Volume Solvers

A novel, efficient, edge-based viscous (EBV) discretization method has been recently developed, implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver, and applied to viscous-kernel computations that include evaluations of meanflow viscous fluxes, turbulence-model and chemistry-model diffusion terms, and the corresponding Jacobian contributions. Initially, the EBV method had been implemented for tetrahedral grids and demonstrated multifold acceleration of all viscous-kernel computations. This paper presents an extension of the EBV method for mixed-element grids. In addition to the primal edges of a given mixed-element grid, virtual edges are introduced to connect cell nodes that are not connected by a primal edge. The EBV method uses an efficient loop over all (primal and virtual) edges and features a compact discretization stencil based on the nearest neighbors. This study verifies the EBV method and assesses its efficiency on mixed-element grids by comparing the EBV solution accuracy and iterative convergence with those of well-established solutions obtained using a cell-based viscous (CBV) discretization method. The EBV solver’s memory footprint is optimized and often smaller than the memory footprint of the CBV solver. A multifold speedup is demonstrated for all viscous-kernel computations resulting in significant reduction of the time to solutions for several benchmark mixed-element-grid computations, including simulations of a flow around NASA’s juncture-flow model and a hypersonic, chemically reacting flow around a blunt body.

CFD↗
NASA NTRSDate not supplied

Edge-Based Viscous Method for Mixed-Element Node-Centered Finite-Volume Solvers

A novel, efficient, edge-based viscous (EBV) discretization method has been recently developed, implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver, and applied to viscous-kernel computations that include evaluations of meanflow viscous fluxes, turbulence-model and chemistry-model diffusion terms, and the corresponding Jacobian contributions. Initially, the EBV method had been implemented for tetrahedral grids and demonstrated multifold acceleration of all viscous-kernel computations. This paper presents an extension of the EBV method for mixed-element grids. In addition to the primal edges of a given mixed-element grid, virtual edges are introduced to connect cell nodes that are not connected by a primal edge. The EBV method uses an efficient loop over all (primal and virtual) edges and features a compact discretization stencil based on the nearest neighbors. This study verifies the EBV method and assesses its efficiency on mixed-element grids by comparing the EBV solution accuracy and iterative convergence with those of well-established solutions obtained using a cell-based viscous (CBV) discretization method. The EBV solver’s memory footprint is optimized and often smaller than the memory footprint of the CBV solver. A multifold speedup is demonstrated for all viscous-kernel computations resulting in significant reduction of the time to solutions for several benchmark mixed-element-grid computations, including simulations of a flow around NASA’s juncture-flow model and a hypersonic, chemically reacting flow around a blunt body.

Edge-based viscous method↗
DOE OSTI2023

Computational diagnostics for flame acceleration and transition to detonation in a hydrogen/air mixture

A new computational diagnostic method for pressure-induced compressibility is proposed by projecting its local contribution to the chemical explosive mode (CEM) in the chemical explosive mode analysis (CEMA) framework. The new method is validated for the study of detonation development during the deflagration-to-detonation transition (DDT) process. The flame characteristics are identified through the quantification of individual CEM contributions of chemical reaction, diffusion, and pressure-induced compressibility. Numerical simulations are performed to investigate the DDT processes in a stoichiometric hydrogen-air mixture. A Godunov algorithm, fifth-order in space, and third-order in time are used to solve the fully compressible Navier-Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model (CDM) described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The new diagnostic method is first applied to onedimensional (1D) canonical flame configurations followed by two-dimensional (2D) simulations of DDT in an obstructed channel where different detonation initiation scenarios are examined using the new CEMA projection formulation. Detailed examinations of the idealized configuration of detonation initiation through shock focusing mechanism at a flame front are also studied using the new formulation. A comparison of the currently proposed CEMA projection and the original formulation by the authors suggests that including the pressure-induced compressibility is essential for the use of CEMA in DDT process. The results also show that the new formulation of CEMA projection can successively capture the detonation initiation through either a gradient mechanism or a direct initiation mechanism, and therefore can be used as an effective local analytical tool for the computational diagnostics of detonation initiation in a DDT process. It was found that detonation development is characterized by a strong contribution of chemistry role to the CEM which is pivotal to the initiation of detonation. The role of compressibility is found enhanced at the edge of the detonation front where diffusion was found to have minimal effects on detonation development.

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

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗
DOE OSTI2025

Equilibrium Fe isotope fractionation between olivine, pyroxene, spinel and MORB glass: Implications for mantle partial melting to generate MORBs

Primitive mid-ocean ridge basalts (MORBs) exhibit Fe isotopic compositions heavier than the upper mantle by +0.074 ± 0.028 ‰ for δ 56 Fe. The processes responsible for this isotopic difference remain unclear. Modeling of Fe isotope fractionation during mantle partial melting requires reliable equilibrium Fe isotope fractionation factors between minerals and melts, for which consistent data are still lacking. Here, in this study, we used Nuclear Resonant Inelastic X-ray Scattering (NRIXS) technique to measure Fe force constants for a MORB glass (ALV 519-4-1) and natural mantle minerals (olivine, orthopyroxene, clinopyroxene, and spinel) to determine the equilibrium Fe isotope fractionation factors between them. The force constants determined in this study, in increasing order, are 167 ± 26 N/m for spinel, 175 ± 17 N/m for olivine, 176 ± 20 N/m for MORB glass, 205 ± 26 N/m for clinopyroxene, and 219 ± 36 N/m for orthopyroxene. We evaluated the previously proposed mechanisms for the heavy Fe isotopic composition of MORBs, including (i) mantle partial melting, (ii) mantle lithological heterogeneity, with pyroxenite in the source, (iii) mantle metasomatism by low-degree melts, and (iv) fractional crystallization of olivine from melts. For (i), we used the pMELTS program to simulate adiabatic decompression melting of mantle peridotites, and calculated Fe isotope fractionation based on Fe 3+ –Fe 2+ equilibrium-controlled fractionation, where Fe 3+ forms stronger bonds and is more incompatible than Fe 2+ . At 10 wt% peridotite melting, corresponding to MORB generation, only +0.03 ‰ Fe isotope fractionation between the melt and the original bulk composition (Δ 56 Fe = δ 56 Fe melt - δ 56 Fe 0 ) was produced, insufficient to account for the observed MORB-upper mantle difference. For (ii), melting of pyroxenites yields smaller Fe isotope fractionation than melting of peridotites, making it unlikely the cause for the MORB-upper mantle isotopic difference. For (iii), both the Fe 3+ /ΣFe ratio and the δ 56 Fe of melts increase with the degree of partial melting, indicating that low-degree melts are not isotopically heavy enough to significantly alter the isotopic composition of lithospheric mantle through metasomatism. For (iv), equilibrium isotope fractionation between olivine and melt is near zero. These results suggest that equilibrium Fe isotope fractionation alone cannot explain the MORB isotopic signature, highlighting the potential role of kinetic isotope fractionation. Using a diffusion model, we calculated kinetic Fe and Mg isotope fractionations associated with (iv) olivine crystallization from a melt, and found that the predicted Fe and Mg isotope fractionations were inconsistent with observations in MORBs. Qualitatively, two processes could have induced kinetic Fe isotope fractionation during MORB generation: (a) Fe-Mg interdiffusion between melt and solid during melt migration and (b) reactive melt-rock interactions during melt focusing. However, a quantitative understanding of their role in modifying the melt isotopic composition remains limited and requires further investigation.

Fe isotopes↗