Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Collisions of localized shocks and quantum circuits

We study collisions between localized shockwaves inside a black hole interior. We give a holographic boundary description of this process in terms of the overlap of two growing perturbations in a shared quantum circuit. The perturbations grow both exponentially as well as ballistically. Due to a competition between different physical effects, the circuit analysis shows dependence on the transverse locations and exhibits four regimes of qualitatively different behaviors. On the gravity side we study properties of the post-collision geometry, using exact calculations in simple setups and estimations in more general circumstances. We show that the circuit analysis offers intuitive and surprisingly accurate predictions about gravity computations involving non-linear features of general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Influence of electronic transport on electrochemical performance of (Cu,Mn) 3 O 4 solid oxide fuel cell cathodes

Alkaline Earth free spinel oxides provide a potential benefit over Sr-doped perovskite-based materials commonly used as electrodes in high-temperature electrochemical energy conversion devices, e.g., solid oxide fuel cells (SOFCs). Sr-segregation is a known issue leading to performance degradation. In this study, Cu x Mn 3-x O 4 (x = 1, 1.2, and 1.5) porous electrodes were examined as SOFC cathodes using electrochemical impedance spectroscopy to investigate the oxygen reduction reaction (ORR) kinetics in relation to the material's intrinsic conductivity, the extrinsic electrode structure, and the cell test design. Similar to the electronic conducting (La,Sr)MnO 3 SOFC cathodes, the ORR kinetics of Cu x Mn 3-x O 4 spinel electrodes was governed by the oxygen adsorption and diffusion at the particle surface as well as the charge transfer at the triple phase boundaries. The overall electrode polarization resistance was highly dependent on contact density with the metallic current collector, active material particle connectivity, electrode thickness, and the intrinsic electronic materials conductivity. Here, we describe the importance of effective electronic charge transport parallel to the electrode surface in maximizing the electrochemically active electrode volume and enhancing electrode performance. We discuss an approach to optimize cell and electrode design with respect to active materials properties. This aspect is critical to ensure reliable evaluation of new materials, since laboratory-scale button-cells typically exhibit a high degree of electrode microstructure (e.g. porosity, thickness) and electrical contact density variation from sample to sample.

(Cu,Mn)3O4↗
DOE OSTI2023

Topological band gap in intercalated epitaxial graphene

The study of functional manipulation of graphene is a critical subject, both for fundamental research and practical applications. In this study, we present that the intercalation of 5d transition metals into epitaxial graphene on SiC is a promising strategy for the realization of topologically nontrivial phases with a finite band gap in graphene. Employing first-principles calculations, grounded in density functional theory, we demonstrate that Re- and Ta-intercalated graphene evolve into two-dimensional topological insulators. These exhibit linear Dirac cones and quadratic bands with topological band gaps, respectively. The emergence of these topological states is attributed to the strong spin–orbit coupling strength of the intercalants. We show that the corresponding topological edge states persist within the finite bulk band gap, aligning with the bulk-boundary correspondence. Furthermore, we explore the spin splitting of the band structure, brought about by the inversion symmetry breaking and the spin–orbit coupling. Finally, our study underscores that the intercalation of graphene is an effective and a feasible approach for manipulating the band gap and the topological nature of graphene. Such intercalated graphene systems hold potential utility for spintronics and low-dimensional quantum device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Lithium Spatial Distribution and Split-Off Electronic Bands at Nanoscale V 2 O 5 /LiPON Interfaces

A combination of depth-resolved cathodoluminescence spectroscopy (DRCLS) and X-ray photoemission depth profiling (XPS) measured the pronounced changes in both the electronic density of states and lithium composition near the nanoscale LixV2O5/LiPON interface. DRCLS studies of electrochemically lithiated bare V2O5 and the sputterdeposited V2O5 plus LiPON overlayer electrochemically lithiated in stages both showed that in the bulk the luminescence intensity of the “split-off” hybridized bonding density of states was anticorrelated with XPSmeasured Li content, decreasing as the Li content increased. However, the LiPON overlayer was found to modify the band structure of the underlying LixV2O5 (LVO) to a depth of at least 30 nm beneath the V2O5 interface. DRCLS spectra near the electrochemically lithiated LiPON/ LVO interface showed a significant intensity of the split-off band, implying a low Li content. However, XPS depth profiling revealed a pronounced negative gradient of Li extending from a maximum Li content at the intimate LiPON boundary to its lowest content of ∼30 nm into the V2O5 in the same region, indicating a strong interaction between band structure and Li electrochemical potential near this heterojunction. These results provide evidence for substantial effects on the local band structure near an electrolyte/cathode interface and insights into the electrochemical interface behavior of solid-state batteries in general.

25 ENERGY STORAGE↗
DOE OSTI2023

Summertime Near-Surface Temperature Biases Over the Central United States in Convection-Permitting Simulations

Convection-Permitting Model (CPM) simulations of the Central United States climate for the summer of 2011 are studied to understand the causes of warm biases in 2-m air temperature (T 2m ) and related underestimates of precipitation including that from mesoscale convective systems (MCSs). Based on 10 CPM simulations and 9 coarser-resolution model simulations, we quantify contributions from evaporative fraction (EF) and radiation to the T 2m bias with both types of models overestimating T 2m largely because they underestimate EF. The performance of CPMs in capturing MCS characteristics (frequency, rainfall, propagation) varies. The pre-summer precipitation bias has large correlation with mean summertime T 2m bias but the relationship between summertime MCS mean rainfall bias and T 2m bias is non-monotonic. Analysis of lifting condensation level deficit and convective available potential energy suggests that models with T 2m warm biases and low EF have too dry and stable boundary layers, inhibiting the formation of clouds, precipitation and MCSs. Among the CPMs with differing model formulations (e.g., transpiration, infiltration, cloud macrophysics and microphysics), evidence suggests that altering the land-surface model is more effective than altering the atmospheric model in reducing T 2m biases. In conclusion, these results demonstrate that land-atmosphere interactions play a very important role in determining the summertime climate of the Central United States.

2-m air temperature↗
DOE OSTI2022

Emergence of Fermi arcs due to magnetic splitting in an antiferromagnet

The Fermi surface plays an important role in controlling the electronic, transport and thermodynamic properties of materials. As the Fermi surface consists of closed contours in the momentum space for well-defined energy bands, disconnected sections known as Fermi arcs can be signatures of unusual electronic states, such as a pseudogap. Another way to obtain Fermi arcs is to break either the time-reversal symmetry or the inversion symmetry of a three-dimensional Dirac semimetal, which results in formation of pairs of Weyl nodes that have opposite chirality, and their projections are connected by Fermi arcs at the bulk boundary. Here, we present experimental evidence that pairs of hole- and electron-like Fermi arcs emerge below the Neel temperature (TN) in the antiferromagnetic state of cubic NdBi due to a new magnetic splitting effect. The observed magnetic splitting is unusual, as it creates bands of opposing curvature, which change with temperature and follow the antiferromagnetic order parameter. This is different from previous theoretically considered and experimentally reported cases of magnetic splitting, such as traditional Zeeman and Rashba, in which the curvature of the bands is preserved. Therefore, our findings demonstrate a type of magnetic band splitting in the presence of a long-range antiferromagnetic order that is not readily explained by existing theoretical ideas.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic study of silicon alloying in the spallation behavior of nanocrystalline aluminum systems

Here, the effect of concentration of the alloying element silicon (Si) on the shock and spall response of aluminum (Al) is studied using classical molecular dynamics (MD) simulations. Silicon is distributed as individual grains and as grain boundaries and the dynamic response under shock is studied for an impact velocity of 1000 m/s. A general trend observed here is that the spallation is almost exclusively intergranular and limited to Al-Si interfaces. The spall strengths are observed to a) decrease with an increase in alloying of Si as grains and b) show little increase for up to 50% grain boundary segregation. An increase in grain size of the whole system showed negligible differences in spall strength, highlighting the significance of concentration of the alloying species rather than its grain size. These results provide an atomistic understanding of light-weight metallic systems under shock compression and pave the way for designing multiphase metal matrix alloys and composites for defense/armor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Can a transport model predict inverse signatures in Lithium metal batteries without modifying kinetics?

In this study, a one-dimensional transport model is developed and analyzed that predicts the inverse signatures obtained in lithium metal batteries. This simple approach’s ability to predict inverse signatures stems from the competing interplay between moving boundary rates and mass transfer limitations to yield inverse signatures without modifying the kinetics. The numerical scheme used for the present model simulations is presented in detail which has been further used to study the effect of design parameters on the prevalence and strength of inverse signatures. It was found that the proposed model and the analysis is more pertinent to thick symmetric cells, commonly used for in-depth fundamental studies.

model, transport, dendrite, lithium metal battery,↗
DOE OSTI2022

An operator-based approach to topological photonics

Abstract Recently, the study of topological structures in photonics has garnered significant interest, as these systems can realize robust, nonreciprocal chiral edge states and cavity-like confined states that have applications in both linear and nonlinear devices. However, current band theoretic approaches to understanding topology in photonic systems yield fundamental limitations on the classes of structures that can be studied. Here, we develop a theoretical framework for assessing a photonic structure’s topology directly from its effective Hamiltonian and position operators, as expressed in real space, and without the need to calculate the system’s Bloch eigenstates or band structure. Using this framework, we show that nontrivial topology, and associated boundary-localized chiral resonances, can manifest in photonic crystals with broken time-reversal symmetry that lack a complete band gap, a result that may have implications for new topological laser designs. Finally, we use our operator-based framework to develop a novel class of invariants for topology stemming from a system’s crystalline symmetries, which allows for the prediction of robust localized states for creating waveguides and cavities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of microstructure and neutron irradiation defects on deuterium retention in SiC

Retention of hydrogen isotopes is a critical concern for operating fusion reactors as retained tritium both activates components and removes scarce fuel from the fuel cycle. Radiation-induced displacement damage in SiC influences the retention of hydrogen isotopes compared to pristine SiC. Deuterium retention in neutron irradiated high purity SiC has been compared to different microstructures of non-irradiated high purity SiC using thermal desorption spectroscopy after gas charging and low energy ion implantation. Experimental results show lower deuterium retention in single crystal SiC than in polycrystal SiC indicating that grain boundaries are key trapping features in unirradiated SiC. Deuterium is released at lower temperatures in neutron irradiated polycrystal SiC compared to pristine polycrystal SiC, suggesting weaker trapping by radiation-induced defects compared to grain boundary trapping sites in the pristine materials. Low energy ion implantation caused a high deuterium release temperature, highlighting the sensitivity of deuterium release behaviour to radiation defect characteristics. First principles calculations have been conducted to identify energetically favourable trapping sites in SiC at the H ABc V Si and H TSi V C complexes, and migration barriers between interstitial sites. This helps interpret experimental results and derive effective diffusivity of hydrogen isotopes in SiC in the presence of vacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2023

MST-8 Colloquium: The Effect of Microstructure on the Spall Strength and Damage of 1045 Steel [Slides]

Three different microstructures of the same steel alloy were investigated to study the influence of microstructure on spall strength and damage. Microstructure influences spall behavior within the explored dynamic regime. Second phase morphology, interphase interface characteristics, and grain boundary character shown to influence spall behavior. Continued understanding of microstructure’s influence on spall and damage needed for targeted material design.

36 MATERIALS SCIENCE↗