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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Establishing the Role of Operating Potential and Mass Transfer in Multicarbon Product Generation for Photoelectrochemical CO 2 Reduction Cells Using a Cu Catalyst

There is increasing interest in the possibility of photoelectrochemical (PEC) reduction of CO 2 to C 2+ products; however, the criteria for maximizing PEC solar-to-C 2+ (STC 2+ ) rates are not well understood. We report here a continuum-scale model of PEC CO 2 reduction (CO 2 R) on Cu in 0.1 M CsHCO 3 and use it to optimize the design and operating conditions for generating C 2+ products. Furthermore, we demonstrate that the potential-dependent product distribution of CO 2 R on Cu requires operating near the potential that maximizes C 2+ generation rates ($V$ id ), unlike PEC water splitting, which desires operation at the maximum photocurrent density. Because of this requirement, the criterion for a high STC 2+ rate includes high-photocurrent semiconductors with photovoltages near $V$ id and low series resistance. The STC 2+ rate in these systems is enhanced by optimal CO 2 transport and exhibits low sensitivity to dirunal solar irradiance variations.

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

Surface Charge in Electrical Double Layer as a Kinetic Descriptor of Electrocatalytic Reactions

The successful commercialization of electrochemical energy-conversion systems hinges on a deeper understanding of electrocatalytic reaction kinetics. Despite extensive research, a key descriptor that characterizes electrolyte effects on reaction kinetics remains elusive. Here, surface charge in electrical double layers (EDLs) is introduced as a descriptor for electrolyte-dependent kinetics. The surface charge is calculated with a continuum EDL model parameterized by density-functional theory. The model is validated by reproducing the anomalously low slope of Pt(111) in Parsons-Zobel plots. Strong correlations are observed between calculated surface charge and experimental kinetic currents for hydrogen evolution, oxygen reduction, and CO 2 -reduction reactions across various pH levels and cationic species. These correlations can be either promotional or inhibitory, depending on solute-intermediate interactions. In acidic media, incorporating adsorbate charge captures specific adsorption effects in oxygen reduction reaction. In conclusion, these findings establish surface charge density as a key descriptor for electrolyte-dependent kinetics, which will guide the design of the electrode/electrolyte interface.

Adsorption↗
DOE OSTI2020

Quadrupole topological photonic crystals

Quadrupole topological phases, exhibiting protected boundary states that are themselves topological insulators of lower dimensions, have recently been of great interest. Extensions of these ideas from current tight binding models to continuum theories for realistic materials require the identification of quantized invariants describing the bulk quadrupole order. Here we identify the analog of quadrupole order in Maxwell’s equations for a gyromagnetic photonic crystal (PhC) through a double-band-inversion process. The quadrupole moment is quantized by the simultaneous presence of crystalline symmetry and broken time-reversal symmetry, which is confirmed using three independent methods: analysis of symmetry eigenvalues, numerical calculations of the nested Wannier bands and the expectation value of the quadrupole operator. Furthermore, we reveal the boundary manifestations of quadrupole phases as quantized edge polarizations and fractional corner charges. The latter are the consequence of a filling anomaly of energy bands as first predicted in electronic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Microstructure dependent shear localization in 316L SS

Shear localization is characterized by the concentration of plastic deformation within narrow bands of material and constitutes the precursor to failure in materials subjected to high strain rates. These bands have been observed experimentally, but classical continuum material models have not been able to fully capture their development. The aim of this research is to investigate the role of microstructural differences in the dynamic shear response of wrought and additively manufactured (AM) 316L SS. Compact Forced Simple Shear specimens for Digital Image Correlation (CFSS-DIC) of wrought and additively manufactured 316L SS were tested on a Split Hopkinson Pressure Bar (SHPB) in conjunction with infrared thermometry measurements. Results indicate that the microstructures of wrought and AM 316L SS exhibit differences that affect the shear response of this materials. In conclusion, CFSS-DIC samples tested on the SHPB did not exhibit recrystallization inside the shear region, fact confirmed by thermometry measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗
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↗