Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nanodiffraction imaging of ionically driven phase separation in cobaltite heterostructures

Controlling the anionic stoichiometry and distribution of perovskite cobaltites presents a powerful method for tailoring the magnetic and electronic properties for magneto-ionic applications. Nevertheless, the influence of nanoscale morphology, including phase separation and defects in cobaltite heterostructures, remains largely unexplored. We used x-ray nanodiffraction to directly image the nanoscale topotactic phase transformation in La 1−x Sr x CoO 3 (LSCO) heterostructures capped with a reactive aluminum (Al) layer (i.e., oxygen “getter”) of varying thicknesses (0.5, 3, and 7 nm). Nanodiffraction images reveal the formation of highly elongated BM filaments extending over mesoscopic lengths (>67 μm), indicating a strong tendency toward oxygen vacancy ordering in both P and BM phases in all heterostructures. Local strain and strain gradient analysis shows a strong correlation between the depletion of the P phase and strain gradients. Our measurements suggest that elongation of the BM filaments is the predominant mechanism for topotactic transformation in LSCO/Al heterostructures, preferred over coarsening, branching, and fractal growth. Our study provides a detailed nanoscale analysis of ionically driven phase transitions in LSCO/Al heterostructures, which is crucial for tailoring the functional properties of cobaltites for magneto-ionic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A skew quadrupole for the AGS to minimize the polarization losses of the polarized beams

The partial helices installed in the Alternating Gradient Synchrotron (AGS) 9 [1] eliminate both, the imperfection and the vertical intrinsic spin resonances thus yielding a 70% polarized proton beam with rigidity 79.4 [Tm] with 2x10 11 protons/bunch at the end of the AGS acceleration cycle. The initial beam polarization at the exit of the 200 MeV Linac is measured to be 80%. The 10% loss of the polarization at the end of the AGS acceleration cycle is due to the horizontal spin resonances which are caused by the betatron oscillations of the beam in the presence of the partial helices. To reduce the effect of these horizontal spin resonances on the beam polarization, the “jump Quads method” is applied in the AGS [2] which eliminates, almost all, the horizontal spin resonances and increases the final polarization of the proton beam to the value of 70%. A study [3] shows that these horizontal spin resonances can be totally eliminated by introducing into the AGS ring, skew quadrupoles which linearly couple the beam motions to excite new horizontal spin resonances which minimize the horizontal-spin-resonances caused by the partial-helices.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Harvest Intensity Effects on Carbon Stocks and Biodiversity Are Dependent on Regional Climate in Douglas-Fir Forests of British Columbia

Temperate forests provide crucial ecosystems services as living sinks for atmospheric carbon (C) and repositories of biodiversity. Applying harvesting at intensities that minimize losses offers one means for mitigating global change. However, little is known of overstory retention levels that best conserve ecosystem services in different regional climates, and likewise as climate changes. To quantify the effect of harvest intensity on C stocks and biodiversity, we compared five harvesting intensities (clearcutting, seedtree retention, 30% patch retention, 60% patch retention, and uncut controls) across a climatic aridity gradient that ranged from humid to semi-arid in the Douglas-fir (Pseudotsuga menziesii) forests of British Columbia. We found that increased harvesting intensity reduced total ecosystem, above ground, and live tree C stocks 1 year post-harvest, and the magnitude of these losses were negatively correlated with climatic aridity. In humid forests, total ecosystem C ranged from 50% loss following clearcut harvest, to 30% loss following large patch retention harvest. In arid forests this range was 60 to 8% loss, respectively. Where lower retention harvests are sought, the small patch retention treatment protected both C stocks and biodiversity in the arid forests, whereas the seedtree method performed as well or better in the humid forests. Below ground C stocks declined by an average of 29% after harvesting, with almost all of the loss from the forest floor and none from the mineral soil. Of the secondary pools, standing and coarse deadwood declined in all harvesting treatments regardless of cutting intensity or aridity, while C stocks in fine fuels and stumps increased. The understory plant C pool declined across all harvesting intensities in the humid forests, but increased in arid forests. Shannon’s diversity and richness of tree and bryoid species declined with harvesting intensity, where tree species losses were greatest in the humid forests and bryoid losses greatest in arid forests. Shrub and herb species were unaffected. This study showed that the highest retention level was best at reducing losses in C stocks and biodiversity, and clearcutting the poorest, and while partial retention of canopy trees can reduce losses in these ecosystem services, outcomes will vary with climatic aridity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Process-microstructure relationship of laser processed thermoelectric material Bi2Te3

Additive manufacturing allows fabrication of custom-shaped thermoelectric materials while minimizing waste, reducing processing steps, and maximizing integration compared to conventional methods. Establishing the process-structure-property relationship of laser additive manufactured thermoelectric materials facilitates enhanced process control and thermoelectric performance. This research focuses on laser processing of bismuth telluride (Bi 2 Te 3 ), a well-established thermoelectric material for low temperature applications. Single melt tracks under various parameters (laser power, scan speed and number of scans) were processed on Bi 2 Te 3 powder compacts. A detailed analysis of the transition in the melting mode, grain growth, balling formation, and elemental composition is provided. Rapid melting and solidification of Bi 2 Te 3 resulted in fine-grained microstructure with preferential grain growth along the direction of the temperature gradient. Experimental results were corroborated with simulations for melt pool dimensions as well as grain morphology transitions resulting from the relationship between temperature gradient and solidification rate. Samples processed at 25 W, 350 mm/s with 5 scans resulted in minimized balling and porosity, along with columnar grains having a high density of dislocations.

Oztan, Cagri↗
DOE OSTI2026

Radial Profiles of the Binary Fraction in Elliptical Galaxies

The radial profile of the binary fraction may vary with environment and is of significant importance for studying the formation mechanisms of binary stars and their dynamical evolution within globular clusters and galaxies. However, existing studies remain limited to the Milky Way and its neighboring galaxies. Leveraging the method proposed by Zhang et al. for estimating the variation in the binary fraction from integrated spectral features, we analyze a sample of 513 elliptical galaxies drawn from the Mapping Nearby Galaxies at Apache Point Observatory survey to measure their radial binary fraction profiles. Our results show that after accounting for the effect induced by radial variations in the stellar population (SP), the median SP-subtracted binary fraction, $r^{\textrm{med}}_{\textrm{b,sub}}$, becomes approximately flat. For nearly all elliptical galaxies in our sample, the variation in the binary fraction relative to the galaxy center at 1R e is less than 5%. No clear correlation is found between the ​​​​​gradient of the binary fraction and the gradients of the SP properties. Moreover, we also compare differences between ultraviolet (UV) upturn and non-UV upturn galaxies. The overall binary fraction profiles and the SP properties of the non-UV upturn galaxies in our sample are comparable to those of the UV upturn galaxies. This similarity may arise from the presence of residual star formation in the non-UV upturn systems.

Li, Xiejin [Chinese Academy of Sciences; Universit↗
DOE OSTI2022

Closed-loop pressure retarded osmosis draw solutions and their regeneration processes: A review

Pressure-Retarded Osmosis (PRO) is an osmotic process that has been used to harvest energy from salinity gradients using a semi permeable membrane. A comparison between open-loop PRO (OLPRO) and closed-loop PRO (CLPRO) was made regarding their performance and costs. In CLPRO, where the diluted draw solution is re-concentrated in the regeneration system to be reutilized in the process, has recently received an intensive focus as the most viable configuration for a standalone power plant. The choice of the PRO draw solution in CLPRO is crucial to garner a high osmotic pressure as the key for the feasibility of the process. Here, in this review, the draw solutions are critically evaluated in the literature in terms of energy output as well as the method of regeneration used to recirculate them. A set of practical criteria has been suggested to appraise the adequacy of the solution for CLPRO application. It was concluded that NH 3 – CO 2 theoretically can produce 170 W/m 2 of power density. Inorganic draw solutes such as NaCl can generate high power density up to 87 W/m 2 . Organic draw solutes with their remarkably low reverse salt flux (RSF) have promising potential for future application in PRO. Similarly, the regeneration systems of the diluted draw solutions have also been reviewed and discussed. How the energy consumption of the regeneration process affects the feasibility of CLPRO is explained. For the specific case of osmotic heat engines (OHEs), when the energy of the regeneration process is supplied by heat waste, the range of applicability of the heat waste in CLPRO in terms of efficiency is defined and compared to Organic Rankine Cycle (ORC). The results showed that CLPRO has better efficiency than ORC for temperatures T < 80 °C, which makes it a promising process or low-grade heat energy recovery. In addition, a PRO-RO hybrid system coupled with solar power can reduce the net specific energy consumption (SEC) to 0.39 kWh/m 3 . The conditions that regeneration processes should operate under to make PRO viable are discussed in the last section. Overall, the study indicates the key factors for optimizing the performance of CLPRO process.

42 ENGINEERING↗
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Climate, food and humans predict communities of mammals in the United States

Abstract Aim The assembly of species into communities and ecoregions is the result of interacting factors that affect plant and animal distribution and abundance at biogeographic scales. Here, we empirically derive ecoregions for mammals to test whether human disturbance has become more important than climate and habitat resources in structuring communities. Location Conterminous United States. Time Period 2010–2021. Major Taxa Studied Twenty‐five species of mammals. Methods We analysed data from 25 mammal species recorded by camera traps at 6645 locations across the conterminous United States in a joint modelling framework to estimate relative abundance of each species. We then used a clustering analysis to describe 8 broad and 16 narrow mammal communities. Results Climate was the most important predictor of mammal abundance overall, while human population density and agriculture were less important, with mixed effects across species. Seed production by forests also predicted mammal abundance, especially hard‐mast tree species. The mammal community maps are similar to those of plants, with an east–west split driven by different dominant species of deer and squirrels. Communities vary along gradients of temperature in the east and precipitation in the west. Most fine‐scale mammal community boundaries aligned with established plant ecoregions and were distinguished by the presence of regional specialists or shifts in relative abundance of widespread species. Maps of potential ecosystem services provided by these communities suggest high herbivory in the Rocky Mountains and eastern forests, high invertebrate predation in the subtropical south and greater predation pressure on large vertebrates in the west. Main Conclusions Our results highlight the importance of climate to modern mammals and suggest that climate change will have strong impacts on these communities. Our new empirical approach to recognizing ecoregions has potential to be applied to expanded communities of mammals or other taxa.

Kays, Roland↗
DOE OSTI2022

Pressure Barrier Applicability to Polymer Flood Design

Effective oil displacement from a reservoir requires adequate and properly directed pressure gradients in areas of high oil saturation. If the polymer bank is too large or too viscous during a polymer flood, the pressure drop from the injection well to the polymer front may act as a pressure barrier by usurping most of the downstream driving force for oil displacement. Polymer injection pressures must be limited. The maximum allowable injection pressure is commonly constrained by caprock integrity, injection equipment, and/or regulations, even though fractures can be beneficial to polymer injectivity (and even sweep efficiency in some cases). This paper examines when the pressure-barrier concept limits the size and viscosity of the polymer bank during a polymer flood. Both analytical and numerical methods are used to address this issue. We examine the relevance of the pressure barrier concept for a wide variety of circumstances, including oil viscosities ranging from 10-cp (like at Daqing, China) to 1650-cp (like at Pelican Lake, Alberta), vertical wells (like at Tambaredjo, Suriname) versus horizontal wells (like at Milne Point, Alaska), single versus multiple layered reservoirs, permeability contrast, and with versus with crossflow between layers. We also examine the relation between the pressure-barrier concept and fractures and fracture extension during polymer injection. We demonstrate that in reservoirs with single layers, the pressure-barrier concept only limits the optimum viscosity of the injected polymer if the mobility of the polymer bank is less than the mobility of the displaced oil bank. The same is true for multi-zoned reservoirs with no crossflow between layers. Thus, for these cases, the optimum polymer viscosity is likely to be dictated by the mobility of the oil bank, unless other factors (like fracture extension) intervene. For multi-zoned reservoirs with free crossflow between layers, the situation is different. A compromise must be reached between injected polymer viscosity and the efficiency of oil recovery. The relevance of our findings is applied to operations for several existing polymer floods. This work is particularly relevant to viscous-oil reservoirs (like Pelican Lake and others) where the injected polymer viscosities are substantially lower than the oil viscosity

Wang, Dongmei↗
DOE OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Particle transport in Markov stochastic mixtures with spatial gradients

Markov geometries are a prototype class of stochastic media that are widely used to model complex disordered systems. In a series of recent works, we have shown that spatially homogeneous (but possibly non-isotropic) Poisson tessellations can be conveniently used in order to generate three-dimensional realizations of Markov media and thus characterize the features of particle transport e.g. in fragmented fuel elements following severe nuclear accidents or Rayleigh-Taylor turbulent mixing for inertial confinement fusion. In this paper we show that Poisson tessellations can be extended to take into account the presence of spatial gradients, which commonly occur in real-world applications due to material stratification. For this purpose, we will provide an explicit construction of spatially-non-homogeneous Poisson tessellations, fully characterize their statistical properties, and finally illustrate the behaviour of particle transport in such random media. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Ion separations with membranes

Ion separations are important for resource recovery, water treatment, and energy production and storage. Techniques such as chemical precipitation, selective adsorption, and solvent extraction are effective, but membranes may separate ions continuously with less waste and lower energy costs. Separation of monovalent and multivalent ions with nanofiltration or electrodialysis membranes already enables water softening and edible salt purification. Similar membranes are attractive as separators in vanadium redox flow batteries. Selective partitioning of divalent counter-ions into ion-exchange membranes even allows transport of these ions against their concentration gradients in salt mixtures. However, separations of ions with the same charge is more challenging. Recent research demonstrated highly selective ion “sieving” at small scales. Here, separations using electrical potentials and differences in ion electrophoretic mobilities are promising, but relatively unexplored. Carrier-mediated transport affords high selectivity in liquid membranes, but these systems are not very stable, and selective transport via hopping between anchored carriers has proven elusive. Finally, this paper discusses how concentration polarization decreases selectivities in many membrane processes. Although development of selective, inexpensive ion-separation membranes is a work in progress, successes in water softening and edible salt purification suggests that future selective membranes will serve as complementary methods to traditional purification techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2024

Efficient generation of grids and traversal graphs in compositional spaces towards exploration and path planning

Abstract Diverse disciplines across science and engineering deal with problems related to compositions, which exist in non-Euclidean simplex spaces, rendering many standard tools inaccurate or inefficient. This work explores such spaces conceptually in the context of materials discovery, quantifies their computational feasibility, and implements several essential methods specific to simplex spaces through a new high-performance open-source library . Most significantly, we derive and implement an algorithm for constructing a novel n-dimensional simplex graph data structure, containing all discretized compositions and possible neighbor-to-neighbor transitions. Critically, no distance or neighborhood calculations are performed, instead leveraging pure combinatorics and order in procedurally generated simplex grids, keeping the algorithm $${\mathcal{O}}(N)$$ O ( N ) , with minimal memory, enabling rapid construction of graphs with billions of transitions in seconds. Additionally, we demonstrate how such graph representations can be combined to homogeneously express complex path-planning problems, while facilitating efficient deployment of existing high-performance gradient descent, graph traversal, and other optimization algorithms.

Krajewski, Adam M. (ORCID:0000000222660099)↗
DOE OSTI2022

Enhanced analysis of experimental x-ray spectra through deep learning

X-ray spectroscopic data from high-energy-density laser-produced plasmas has long required thorough, time-consuming analysis to extract meaningful source conditions. There are often confounding factors due to rapidly evolving states and finite spatial gradients (e.g., the existence of multi-temperature, multi-density, multi-ionization states, etc.) that make spectral measurements and analysis difficult. Here, in this paper, we demonstrate how deep learning can be applied to enhance x-ray spectral data analysis in both speed and intricacy. Neural networks (NNs) are trained on ensemble atomic physics simulations so that they can subsequently construct a model capable of extracting plasma parameters directly from experimental spectra. Through deep learning, the models can extract temperature distributions as opposed to single or dual temperature/density fits from standard trial-and-error atomic modeling at a significantly reduced computational cost compared to traditional trial-and-error methods. These NNs are envisioned to be deployed with high repetition rate x-ray spectrometers in order to provide detailed real-time analysis of experimental spectra.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗