Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chemical looping air separation with Sr 0.8 Ca0.2Fe 0.9 Co 0.1 O 3-δ perovskite sorbent: Packed bed modeling, verification, and optimization

Chemical looping air separation (CLAS) represents a promising approach for efficient O 2 production from the air. This present study aims at optimizing the absorber/desorber operations and the separation process with extensive experimental validation. Specifically, a one-dimensional packed bed model was developed to investigate the CLAS operation with a Sr 0.8 Ca 0.2 Fe 0.9 Co 0.1 O 3-δ perovskite sorbent. The redox thermodynamics of perovskite sorbent was measured by TGA and then incorporated into a linear driving force model to describe the O 2 absorption and desorption rates. Both 4-step and 5-step air separation cycle configurations, with various cyclic structures, were performed in a subpilot-scale packed bed. The model predicted O2 purity and productivity were consistent with experimental results, supporting its accuracy and applicability. Parametric analysis and multi-objective optimization were further carried out to assess the performance of CLAS. Both O 2 purity and recovery increased monotonically with the cycle time, airflow rate, steam flow rate, and absorption pressure. Meanwhile, optimal O 2 productivity and power consumption can only be achieved by specific combinations of these parameters. The optimized results showed that CLAS can be highly competitive when compared to conventional pressure swing adsorption (PSA) or cryogenic distillation. The 5-step cycle configuration achieved a minimum power consumption of 118 kW·h for producing 1 ton O 2 with ≥ 95% purity. The maximum O 2 productivity reached 0.0932 g O2 /(g sorbent ·h) with 390 kW·h/ton O 2 of energy consumption (95% pure). The optimization results also indicate that CLAS can potentially be more efficient than cryogenic distillation even when the required O 2 purity is above 99%.

42 ENGINEERING↗
DOE OSTI2025

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗
DOE OSTI2021

Hinderance of C-S-H sheet piling during first drying using a shrinkage reducing agent: A SAXS study

The role of a shrinkage reducing agent (SRA) in the hardened cement paste during a first desorption process was investigated by measuring the small-angle X-ray scattering (SAXS). Fractal disc analysis and Guinier approximation were adopted to interpret the SAXS profiles. We experimentally confirmed that the presence of SRA molecules mitigated the piling of calcium-silicate hydrate (C-S-H) sheets, which could be monitored by the increase in the thickness of the disc, representing the agglomeration of C-S-H sheets. In addition, the basal spacing of C-S-H did not change during the drying. By contrast, the diameter of the discs monotonically decreased and showed a similar trend to that of samples lacking SRA. Based on these findings, we concluded that SRA molecules were on the surface of the C-S-H sheets and mitigated the piling, facilitating water removal from the C-S-H agglomeration in the RH range of 80–40%.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

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

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2020

Rheological and microstructural characterization of aqueous suspensions of carbon black and reduced graphene oxide

Carbon black is often used as a conductivity enhancer in battery electrodes. Reduced graphene oxide is another conducting material, with a high aspect ratio sheet-like morphology, and offers the possibility of generating a conducting network at loadings that are small compared to carbon black. Suspensions of conducting carbon, active materials and binders are typically coated on current collectors during the formation of electrodes. The rheology of these suspensions is an important indicator of the connectivity of particles in the suspension and impacts the electrical conductivity as well as the thickness of the coated layer. In this work, the rheology and microstructures of aqueous suspensions of para-aminobenzoic acid terminated carbon black (CB) with and without reduced graphene oxide (RGO) were investigated. The CB loading varied between 0.05 wt.% -10 wt.%, while the samples containing RGO had an RGO loading of 0.05 wt.%. The pH of the suspension was varied from 7.5 to 3. The carboxyl groups on the CB surface are deprotonated at pH 7.5, making the particles hydrophilic. Protonation of the carboxylate groups at pH 3 makes the particles partially hydrophobic. Suspensions of CB showed Newtonian behavior at all loadings for pH 7.5, and at loadings below 1.5 wt.% for pH 3. They are shear thinning at loadings between 1.5–10 wt.%. Interestingly, the low shear viscosity is non-monotonic with CB loading, rising up to a CB loading of 7.5 wt.%, and then dropping at 10 wt.%. Cryogenic scanning electron microscopy (SEM) showed CB particles forming aggregates at pH 3 for CB loadings greater than 1.5 wt.%. At pH 3, addition of 0.05 wt.% RGO resulted in a viscosity increase for suspensions containing 0.5 wt.% CB, as the RGO formed connections with CB particles that enhanced network formation. Adding RGO had little effect on suspensions that had 1.5 wt.% CB. The combination of SEM imaging and suspension rheology provides insight into the particle connectivity for the carbon-based materials studied here. These insights will be useful in determining coating conditions for slurries used in rechargeable battery electrodes.

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

Repetitively pulsed nanosecond discharge plasma decay in propane–oxygen gas mixture in the presence of a heating electric field

Plasma decay in the afterglow of a repetitively pulsed nanosecond discharge in a stoichiometric propane–oxygen mixture was experimentally investigated when a weak heating DC electric field was applied and in its absence. The discharge was ignited at room gas temperature and a pressure of 1–2 Torr and was characterized by low specific energy inputs (<0.004 eV per molecule in one pulse). Using microwave interferometry, the temporal evolution of the electron density during plasma decay was studied, and the effective recombination coefficients were obtained from data processing. It was shown that the rate of plasma decay behaved in a non-monotonic manner with increasing degree of propane oxidation; at first the decay rate grew, then passed through a maximum, fell and saturated in the limit of a large (~2000) number of pulses. In this limit, the effect of the heating DC electric field on the plasma decay decreased with approaching chemical equilibrium. Numerical simulation of the observed effects was performed for low and high oxidation degrees of propane taking into account changes in the composition of positive ions in the plasma. Good agreement was obtained between measurements and calculations of the electron density during plasma decay in these cases. Here it was shown that the formation of cluster ions in the discharge afterglow plays a fundamental role. The plasma decay was controlled by electron recombination with hydrocarbon cluster ion at low oxidation degree of propane and with water cluster (hydrated) ions at high oxidation degree. A hypothesis was proposed to explain the observed nonmonotonic behavior of the plasma decay rate with an increase in the propane oxidation in the discharge, based on the formation of hydrated hydrocarbon ions C x H y + (H 2 O) k at moderate oxidation degrees.

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

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗
DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗
DOE OSTI2025

Plasma assisted NH 3 /H 2 /air ignition in nanosecond discharges with non-equilibrium energy transfer

Ammonia (NH 3 ), with its high energy density and easiness to store and transport as a hydrogen carrier, has become a promising alternative green fuel. However, its adoption in power generation is hindered by challenges such as low burning velocity, slow low-temperature oxidation, high NO x emissions, and ignition difficulty. Here, this work computationally investigates the effects of non-equilibrium energy transfer by nanosecond discharges on NH 3 ignition and flame propagation in an NH 3 /H 2 /air flow at 700 K and 1 atm. The simulation results demonstrate that NH 3 /air mixtures require a large ignition energy due to their large critical ignition radius. It is shown that adding 30 % hydrogen significantly reduces the critical ignition radius and minimum ignition energy. Two-dimensional modeling further shows a non-monotonic dependence of ignition kernel volume on the applied voltage and reduced electric field. The optimum ignition enhancement occurs at 200 Td where the generation of electronically excited species and radicals including N 2 (B), O( 1 D) and OH becomes most efficient. Higher voltages divert electron energy toward ionization, which makes it less effective for NH 3 ignition. The study also identifies an optimal electrode gap size for a given pulse energy. Smaller gap sizes increase deposited energy density, raising temperature and radical concentrations. However, excessive reduction of the gap distance reduces flame propagation speed due to the flame stretch effect in rich mixtures with the effective Lewis number greater than unity. A nonlinear relationship between pulse repetition frequency and ignition kernel volume is observed in a nanosecond pulsed high frequency discharge (NPHFD). An optimal frequency range of 200 kHz to 2 MHz is found when two pulses are used. In addition, an optimal number of pulses exists for each pulse repetition frequency, with higher frequencies requiring more pulses to maximize the overlap region. These findings provide critical insights on developing controlled plasma discharge techniques for efficient NH 3 ignition in reactive flows within internal combustion engines and gas turbines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Study on ammonia and dimethyl ether oxidation and kinetic interaction up to 100 atm

NH 3 and DME dual fuel oxidation and kinetic coupling are experimentally studied by using a supercritical pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm, over a temperature range of 500–900 K, and at fuel-lean and stoichiometric conditions with NH 3 to DME molar ratios of 4 and 0.62, respectively. An HP-Mech model for high-pressure NH3/DME oxidation is developed based on our previous studies and it shows generally better performance than other models in the literature on high-pressure oxidation experiments. Due to DME’s strong low-temperature reactivity, it dramatically enhances the oxidation of NH3 at low temperature. However, the effect of NH 3 on DME oxidation varies with temperature. At low temperatures, NH 3 inhibits low-temperature reactivity by consuming OH radicals. In addition, the NH 2 and NO x formation from NH 3 further suppresses the low-temperature DME reactivity by reducing alkylperoxyl and O 2 QOOH radicals via RO 2 + NH 2 = RO + H 2 NO, RO 2 + NO = RO + NO 2 and O 2 QOOH + NO = 2CH 2 O + HO 2 + NO 2 . At intermediate temperatures, due to enhanced kinetic coupling of NH 2 /NO x /HO 2 chemistry, DME oxidation is significantly promoted. It is found that there are two major NH 2 /NO x /HO 2 coupling pathways for OH radical production brought by NH 3 : (1) NH 2 + HO 2 = H 2 NO + OH and (2) NH 2 + NO 2 = H 2 NO + NO with NO + HO 2 = OH + NO 2 . Moreover, the resulting H 2 NO will further contribute to OH production via the H 2 NO/NO x /HONO coupling pathways: H 2 NO + NO 2 = HNO + HONO, HNO + NO 2 = HONO + NO, and HONO (+M) = OH + NO (+M). These new NH 2 /NO x /HO 2 and H 2 NO/NO x /HONO pathways play a critical role in promoting DME oxidation at intermediate temperatures. Novelty and Significance Statement: Blending NH 3 with DME can significantly enhance the reactivity of NH 3 and facilitate its practical application in advanced internal combustion engines. NH 3 /DME oxidation is studied in a novel supercritical pressure jet-stirred reactor up to 100 atm, a much higher pressure than previous studies. An interesting non-monotonic effect of NH3 addition on fuel oxidation has been found. Unique kinetic couplings of NH2/HO2/NOx and H2NO/NOx/HONO brough by NH3 and DME blending are identified to play important roles in radical generation and reactivity promotion. The valuable experimental data and key kinetics revealed in the current work can tremendously improve our understanding of NH 3 /DME oxidation via NO x /RO 2 and NH 2 / H 2 NO/HO 2 /NO x kinetic couplings at low to intermediate temperatures.

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

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗
DOE OSTI2020

Investigation of the discontinuous precipitation of U-Nb alloys via thermodynamic analysis and phase-field modeling

U-Nb’s discontinuous precipitation, $\gamma^{bcc}_{matrix}$, $\alpha^{orth}_{cellular}$, $\gamma\prime^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $\gamma\prime$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb’s discontinuous precipitation hypothesized that the energy of $\alpha$ forms an intermediate common tangent with the first potential of the double-well energy of $\gamma$ at the $\gamma\prime$ composition. In this work, we examine different possible mechanisms by which the discontinuous precipitation product in the U-Nb system can be stabilized. We put forward a mechanism by which the bulk free energy of the $\gamma$ can develop a non-monotonic dependency with composition. Additionally we posit that local contributions due to lattice mismatch between the $\alpha$ and $\gamma$ phases may be responsible for the generation of metastable states that may stabilized by thermodynamics as well as by kinetics. Our work suggests that local misfit strain tends to play a crucial role in the growth of the discontinuous precipitation product. Depending on the magnitude of strain developed at the $\alpha$/$\gamma\prime$ interfaces, either an increasing $\gamma\prime$ composition or a random distribution of $\gamma\prime$ compositions around the equiatomic value with respect to increasing temperature could be expected. Moreover, we show how it is possible to stabilize the discontinuous precipitation front through highly anisotropic and fast interface diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Phase-field modeling of radiation-induced composition redistribution: An application to additively manufactured austenitic Fe–Cr–Ni

Multicomponent alloys undergoing irradiation damage develop radiation-induced composition redistribution at point defect sinks such as grain boundaries (GBs) and dislocations. Such redistribution results in undesired changes to their mechanical behavior and corrosion resistance. Additively manufactured alloys proposed for future nuclear applications are expected to demonstrate a distinct response to irradiation owing to their unique microstructure with as-solidified dislocation density and chemical microsegregation. To capture the composition redistribution in such systems, we develop a mesoscale model with coupled evolution of atomic and point defect components in the presence of dislocation density, dislocation heterogeneity, and thermodynamic interactions at the GB. The model is parameterized for an FCC Fe–Cr–Ni alloy as a representative system for austenitic stainless steels, and simulations are performed in 1D and 2D as a function of irradiation temperature, dose, dislocation density, and grain size. Radiation-induced segregation (RIS) characterized by Cr depletion and Ni enrichment is predicted at both the GB and the dislocation cell wall, with RIS being lower in magnitude but wider at the cell wall. Strongly biased absorption of self-interstitials by dislocations is found to suppress Ni enrichment but slightly enhance Cr depletion under certain conditions. Thermodynamic segregation at the GB is predicted to be narrower and opposite in sign to RIS for both Cr and Ni. Importantly, non-monotonic segregation is found to occur when both thermodynamic and RIS mechanisms are considered, providing a novel physical interpretation of experimental observations. The model is expected to serve as a key tool in accelerated qualification of irradiated materials.

additively manufactured microstructure↗
DOE OSTI2025

Impact of amorphous pockets on displacement damage evolution in silicon

Silicon has long been known to exhibit amorphization in response to heavy particle bombardment. For doses below the total amorphization threshold, partial amorphization is observed in the form of scattered amorphous pockets. While extensive research has gone into modeling the formation and evolution of amorphous pockets in response to irradiation, no studies yet investigate their impact on the evolution of other damage such as interstitial supersaturation and clustering. In this study, we survey the impact of amorphous pockets on defect evolution in silicon when treated as static sinks. MD is first used to show that amorphous pockets provide energetically favorable sites for point defects relative to the crystalline bulk, supporting the hypothesis that they act as sinks. A 0-D cluster dynamics model is then constructed, taking an interstitial clustering model from the literature and including amorphous pockets as a sink species. We conduct our survey for temperatures between 30 and 400 °C and sink strengths between 1 to 6 x 10 10 cm −2 . Both implantation- and radiation-induced damage states are investigated using interstitial and vacancy concentrations as initial condition variables. We find that, due to the differing migration rates of the interstitial and the vacancy, amorphous pockets have a non-monotonic impact on the final damage state depending on the effective sink strength of the amorphous pockets, resulting in increased damage formation in regimes of intermediate amorphization. In conclusion, this result emphasizes the important role of amorphous pockets in governing the evolution of damage in partially amorphized crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗