Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spin-textured Chern bands in AB-stacked transition metal dichalcogenide bilayers

Significance The long-wavelength moiré superlattices in twisted two-dimensional (2D) structures have emerged as a highly tunable platform for strongly correlated electron physics. In this work, we reveal the displacement field–induced topological moiré bands in AB-stacked transition-metal dichalcogenide (TMD) heterobilayer M o T e 2 / W S e 2 . Valley-contrasting Chern bands with nontrivial spin texture are formed from interlayer hybridization between M o T e 2 and W S e 2 bands of nominally opposite spins. The general mechanism for band inversion we uncovered in AB-stacked TMD bilayers is quite robust and can be applied to other material combinations as well. Our study establishes a recipe for creating topological bands in AB-stacked TMD bilayers in general, which provides a highly tunable platform for realizing quantum-spin Hall and interaction-induced quantum anomalous Hall effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗
DOE OSTI2023

Helium bubble size effects on the surface morphological response of plasma-facing tungsten

Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗
DOE OSTI2021

Microscopic investigation of the 8 Li(n,γ)Li 9 reaction

Background: The 8 Li (n, γ) 9 Li reaction plays an important role in several astrophysics scenarios. It cannot be measured directly and indirect experiments have so far provided only cross section limits. Theoretical predictions differ by an order of magnitude. Purpose: In this work we study the properties of 9 Li bound states and low-lying resonances and calculate the 8 Li(n, γ) 9 Li cross section within the no-core shell model with continuum (NCSMC) with chiral nucleon-nucleon and three-nucleon interactions as input. Methods: Here, the NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states well suited to calculate low-energy nuclear scattering and reactions. For the capture cross section calculation, we adjust calculated thresholds to experimental values. Results: Our calculations reproduce the experimentally known bound states as well as the lowest 5/2 – resonance of 9 Li. We predict a 3/2 – spin-parity assignment for the resonance observed at 5.38 MeV. In addition to the very narrow 7/2 – resonance corresponding presumably to the experimental 6.43 MeV state, we find several other broad low-lying resonances. Conclusions: Our calculated 8 Li(n, γ) 9 Li cross section is within the limits derived from the 1998 National Superconducting Cyclotron Laboratory Coulomb-dissociation experiment. However, it is higher than cross sections obtained in recent phenomenological studies. It is dominated by a direct E1 capture to the ground state with a resonant contribution at ≈ 0.2 MeV due to E2/M1 radiation enhanced by the 5/2 – resonance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Ab initio calculation of the β decay from 11 Be to a 10 Be + $\textit{p}$ resonance

In this work, the exotic β -delayed proton emission is calculated in 11 Be from first principles using chiral two- and three-nucleon forces. To investigate the unexpectedly large branching ratio measured in Ayyad et al. [Phys. Rev. Lett. 123, 082501 (2019)], we calculate the proposed (1/2 + , 1/2) proton resonance in 11 B using the no-core shell model with continuum. This calculation helps to address whether this enhancement is caused by unknown dark decay modes or an unobserved proton resonance. We report a branching ratio of $b_p$ = (1.3 ± 0.5) × 10 –6 , suggesting that its unexpectedly large value is caused by an unobserved proton resonance in 11 B.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Search for an s -wave resonance in Li 7 just above the proton-decay threshold

Near-threshold resonances play an important role in nucleosynthesis and applied nuclear science. The study of nuclei removed from stability has greatly extended the list of resonances close to decay thresholds. The no core shell model with continuum (NCSMC) recently predicted an s-wave resonance just above the proton-decay threshold of 7 Li at an excitation energy of 10.2 MeV. This potential case of a near-threshold resonance is dependent on the quantum mechanics of the p + 6 He fragments extended into the continuum. Further, the 6 He(d,n) 7 Li* reaction was employed in an attempt to populate this resonance and search for its proton decay via the invariant-mass method. No evidence of this resonance was found. However, the data collected in this search led to better constraints on the energy and width of the isobaric analog state (IAS) with (J π ,T) = (3/2 - , 3/2) and revealed a new weak resonance just above the IAS in energy, predicted by the NCSMC as a (J π , T ) = (3/2 – , 1/2) state and potentially a part of the antianalog structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Coalescence dynamics of prismatic dislocation loops due to vacancy supersaturation

Recovery processes in metals with a high defect concentration are at the heart of a structural alloy's usability in engineering applications. However, such processes take place via nonconservative mechanisms mediated by point defect diffusion and are governed by their interactions with the underlying microstructure. Indeed, one of the main challenges in studying recovery processes is their intrinsic multiscale complexity, which often negates the direct use of experimental techniques. Here, in this Research Letter, we study the mechanism of dislocation loop coalescence in metals by simulating the interaction dynamics between vacancies and prismatic loops using a stochastic framework that seamlessly bridges point defect kinetics with dislocation mechanics. Our simulations reveal two fundamental discoveries: first, that climb processes are governed by local fluctuations that are not captured by continuum elastodiffusion models, and, second, that prismatic loops create internal traps for vacancies, leading to the formation of rings around the perimeter of the loops that affect how coalescence takes place. Our results show an excellent agreement with in situ experimental observations of loop coalescence in Fe under irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic control of molecular transport MXene transistor membranes

Controlled spatial confinement and surface properties of lamellar 2D nanomaterial membranes could enhance many precision separation processes. Traditionally, researchers view channel dimensions, surface properties, and permeation rates of these membranes as intrinsic properties that cannot be modulated in operando. We report that gate voltage applied to the conducting laminar MXene membrane can modulate the permeation rate of ions and neutral solutes, as well as its effective size rejection. In operando wide-angle x-ray scattering measurements reveal that these changes are not driven by electrically induced variations in the d spacing of the MXene layers. Instead, experimental data and continuum electrokinetic modeling reveal that ion transport through the MXene channels is primarily affected by Donnan equilibrium at the membrane-solution interface. We also report a strong increase in the permeation rates through the membrane under a low-frequency ac voltage gating regime that we attribute to diffusioosmotic flow oscillations induced in the membrane. Overall, MXene “transistor” membranes provide a previously unidentified approach to dynamic control of molecular separations.

Pendse, Aaditya [Lawrence Livermore National Labor↗
DOE OSTI2023

Co-Design of Multijunction Photoelectrochemical Devices for Unassisted CO 2 Reduction to Multicarbon Products

Photoelectrochemical (PEC) CO 2 reduction (PEC CO 2 R) is a prospective approach for utilizing solar energy to synthesize a variety of carbon-containing chemicals and fuels, the most valuable of which are multicarbon (C 2+ ) products, such as ethylene and ethanol. While these products can be produced with high faradaic efficiency using Cu, this occurs over a relatively narrow potential range, which, in turn, imposes constraints on the design of a device for PEC CO 2 R. Herein, we used continuum-scale modeling to simulate the solar-to-C 2+ (STC 2+ ) efficiency of PEC CO 2 R devices fed with CO 2 -saturated, 0.1 M CsHCO 3 . We then explored how cell architecture and the use of single or dual photoelectrode(s) alters the optimal combination of photoelectrode bandgaps for high STC 2+ efficiency. Ultimately, this work provides guidance for the co-design of the device architecture and photoelectrode bandgaps required to achieve high STC 2+ efficiency. The insights gained are then used to identify systems that yield the highest amount of C 2+ products throughout the day and year.

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

Solid Electrolyte Interphase Formation at the Ionic Liquid Electrolyte—Lithium-Metal Interface Using an Ab Initio Molecular Dynamics Approach

An interfacial study is performed using ab initio molecular dynamics (AIMD) simulations to elucidate the electrochemical interfacial stability of an ionic liquid (IL) as component of an electrolyte in contact with a Li-metal anode when an additive, vinylene carbonate (VC), is added to the electrolyte. The IL electrolyte is composed of 1-butyl-1-methylpyrrolidinium bis(fluoro-sulfonyl)-imide (PYR14 + FSI − ) with lithium bis(trifluoromethane-sulfonyl)imide salt (Li + TFSI − ) in a 8:2 concentration ratio. Reactions observed at the interface electrolyte/Li-metal during the AIMD simulations are verified by electronic structure calculations using density functional theory calculations considering a continuum solvent model. We found that the addition of VC to the IL prevents further reduction of the counter-ion TFSI − and retards reduction of FSI − in comparison to the case when the additive VC is absent. The PYR14 + reduction is never observed with or without the additive. The 1st, 2nd, and 3rd electron affinities of the PYR14 + are −0.08, −0.36 and −3.65 eV, indicating little to large repulsive behavior against electrons. The opposite is observed for the two negative ions, FSI − and TFSI − , which yield 1st electron affinities of 1.73 and 1.83 eV, respectively. Reactions obtained under open circuit conditions. Dissociation mechanisms are presented for all IL components reduced at the IL/Li-metal interface.

Electrochemistry↗
DOE OSTI2024

Communication—First-Principles Simulations of LiPF 6 Decomposition in Ethylene Carbonate-Based Electrolytes

We revisit a theoretical result by Okamoto (2013 Journal of The Electrochemical Society , 160 , A404) who calculated the energy barrier for the decomposition of lithium hexafluorophosphate (LiPF 6 ) into LiF + PF 5 when solvated in Ethylene carbonate (EC)-based electrolyte. Using different numerical techniques to discretize the Density Functional Theory (DFT) equations, and different continuum solvation models with the same dielectric constant, our results largely confirm the original calculation. However, simulations with a higher dielectric permittivity value, closer to that of EC, show a lower energy barrier. More importantly, First-Principles simulations with an explicit solvent show a substantially lower energy barrier.

25 ENERGY STORAGE↗
DOE OSTI2025

Role of Pb/PbSO 4 Morphologies in Dynamic Charge Acceptance of Lead Acid Batteries

The dynamic charge acceptance (DCA) of negative pastes in lead acid batteries is studied using a continuum level model that is validated against in situ X-ray diffraction experiments. The model matches the current transient during a ten-minute, 2.4 V charge at 80% state of charge (SOC) using a distribution of PbSO 4 particles. <2% by volume of small, reactive particles is shown to provide ∼50% of the capacity during the first ten seconds. Increases in DCA with increasing carbon content are captured with increasing lead surface area. The increase in DCA when the 2.4 V charge is done after discharging to 80% SOC (vs charging to 80%) is captured by reducing the PbSO 4 particle sizes. It is hypothesized that the lower DCA after charging is caused by the absence of smaller particles that were consumed during the initial charge to 80% SOC, which is confirmed from X-ray measurements. Finally, two-particle case studies show that reducing the PbSO 4 particle size is more effective for increasing DCA than increasing the PbSO 4 dissolution rate alone. For a given lead surface area, the optimal DCA occurs when the particle sizes are large enough to avoid blocking the surface but small enough to provide dissolution/diffusion improvements.

Knehr, Kevin W. [Argonne National Laboratory (ANL)↗
DOE OSTI2026

Optimized Tandem Catalyst Patterning for CO 2 Reduction Flow Reactors

Tandem catalysis involves two or more catalysts arranged in proximity within a single reaction vessel, with the aim of synergistically aligning the catalysts’ reaction pathways to maximize overall system performance. This study presents a proof of concept showing the integration of continuum transport modeling with design optimization in a simplified two-dimensional flow reactor setup for electrochemical CO 2 reduction. Ag catalysts provide the CO 2 ⟶ CO reaction capability, and Cu catalysts provide the CO ⟶ high-value products reaction capability. Given a set of input parameters, the optimization algorithm uses adjoint methods to modify the Ag/Cu surface patterning in order to maximize the current density toward high-value products, such as ethylene. The optimized designs yield significant performance enhancement especially at more negative applied voltages (i.e., stronger surface reactions) and for larger numbers of patterning sections. For an applied voltage of −1.7 V vs. SHE, the 12-section optimized design increases the current density toward ethylene by up to 65% compared to the unoptimized 2-section design. For the optimized cases, observed differences in the production and consumption of CO (the key intermediate species) and minimized zones of low CO reactant surface concentration on Cu sections explain the improved reactor performance.

CO2 reduction↗