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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Generic Hilbert space fragmentation in Kogut-Susskind lattice gauge theories

At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in d+1D (d spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Lastly, our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.

Eigenstate thermalization↗

DOE OSTI2022

The black carbon cycle and its role in the Earth system

Black carbon (BC) is produced by incomplete combustion of biomass by wildfires and burning of fossil fuels. BC is environmentally persistent over centuries to millennia, sequestering carbon in marine and terrestrial environments. However, its production, storage and dynamics, and therefore its role in the broader carbon cycling during global change, are poorly understood. Here, in this Review, we discuss BC cycling across the land-to-ocean continuum. Wildfires are the main source of BC, producing 128 ± 84 teragrams per year. Negative climate–BC feedbacks could arise as wildfire increases with anthropogenic warming, producing more BC, which in turn will sequester carbon, but the magnitude of these effects are unknown. Most BC is stored in terrestrial systems with some transported to the ocean via rivers and the atmosphere. However, the oceanic BC budget is not balanced, with known BC removal fluxes exceeding BC inputs. We demonstrate these observed inconsistencies using a simple ocean box model, which highlights key areas of future research. Measurements of BC mineralization and export rates along the land-to-ocean continuum and quantification of previously unexplored sources of oceanic BC are needed to close the global BC budget.

aquatic↗

DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗

DOE OSTI2024

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ± 0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. Furthermore, the An +m /An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Chemical calculations↗

DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗

DOE OSTI2019

Stochastic self-tuning hybrid algorithm for reaction-diffusion systems

Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Furthermore, its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Thermal-fluctuation effects on small-scale statistics in turbulent gas flow

Kolmogorov's theory of turbulence assumes that the small-scale turbulent structures in the energy cascade are universal and are determined by the energy dissipation rate and the kinematic viscosity alone. However, thermal fluctuations, absent from the continuum description, terminate the energy cascade near the Kolmogorov length scale. Here, we propose a simple superposition model to account for the effects of thermal fluctuations on small-scale turbulence statistics. Furthermore, for compressible Taylor–Green vortex flow, we demonstrate that the superposition model in conjunction with data from direct numerical simulation of the Navier–Stokes equations yields spectra and structure functions that agree with the corresponding quantities computed from the direct simulation Monte Carlo method of molecular gas dynamics, verifying the importance of thermal fluctuations in the dissipation range.

42 ENGINEERING↗

DOE OSTI2020

Black hole spin in X-ray binaries: giving uncertainties an $\mathcal{f}$

The two established techniques for measuring black hole spin in X-ray binaries often yield conflicting results, which must be resolved before either method may be deemed robust. In practice, black hole spin measurements based on fitting the accretion disc continuum effectively do not marginalize over the colour-correction factor f col . This factor parametrizes spectral hardening of the disc continuum by the disc atmosphere, whose true properties are poorly constrained. We incorporate reasonable systematic uncertainties in f col into the eight (non-maximal) black hole spin measurements vetted by the disc continuum fitting community. In most cases, an f col uncertainty of ±0.2–0.3 dominates the black hole spin error budget. We go on to demonstrate that plausible departures in f col values from those adopted by the disc continuum fitting practitioners can bring the discrepant black hole spins into agreement with those from iron line modelling. Systematic uncertainties in f col , such as the effects of strong magnetization, should be better understood before dismissing their potentially dominant impact on the black hole spin error budget.

4U 1543–47↗

DOE OSTI2021

Appearance versus disappearance of broad absorption line troughs in quasars

ABSTRACT We present a new set of 84 broad absorption line (BAL) quasars (1.7 < zem < 4.4) exhibiting an appearance of C iv BAL troughs over 0.3–4.8 rest-frame years by comparing the Sloan Digital Sky Survey Data Release (SDSSDR)-7, SDSSDR-12, and SDSSDR-14 quasar catalogues. We contrast the nature of BAL variability in this appearing BAL quasar sample with a disappearing BAL quasar sample studied in the literature by comparing the quasar’s intrinsic, BAL trough, and continuum parameters between the two samples. We find that appearing BAL quasars have relatively higher redshift and smaller probed time-scales as compared to the disappearing BAL quasars. To mitigate the effect of any redshift bias, we created control samples of appearing and disappearing BAL quasars that have similar redshift distribution. We find that the appearing BAL quasars are relatively brighter and have shallower and wider BAL troughs compared to the disappearing BAL sample. The distribution of quasar continuum variability parameters between the two samples is clearly separated, with the appearance of the BAL troughs being accompanied by the dimming of the continuum and vice versa. Spectral index variations in the two samples also point to the anticorrelation between the BAL trough and continuum variations consistent with the ‘bluer when brighter’ trend in quasars. We show that the intrinsic dust model is less likely to be a favourable scenario in explaining BAL appearance/disappearance. Our analysis suggests that the extreme variations of BAL troughs like BAL appearance/disappearance are mainly driven by changes in the ionization conditions of the absorbing gas.

Mishra, Sapna↗

DOE OSTI2024

Entanglement entropy of (2 +1)-dimensional SU(2) lattice gauge theory on plaquette chains

We study the entanglement entropy of Hamiltonian SU(2) lattice gauge theory in 2 +1 dimensions on linear plaquette chains and show that the entanglement entropies of both ground and excited states follow Page curves. The transition of the subsystem size dependence of the entanglement entropy from the area law for the ground state to the volume law for highly excited states is found to be described by a universal crossover function. Quantum many-body scars in the middle of the spectrum, which are present in the electric flux truncated Hilbert space, where the gauge theory can be mapped onto an Ising model, disappear when higher electric field representations are included in the Hilbert space basis. This suggests the continuum (2+1)-dimensional SU(2) gauge theory does not have such scarred states.

Astronomy & Astrophysics↗

DOE OSTI2022

Semileptonic form factors for $$B\rightarrow D^*\ell \nu $$ at nonzero recoil from $$2+1$$-flavor lattice QCD: Fermilab Lattice and MILC Collaborations

Abstract We present the first unquenched lattice-QCD calculation of the form factors for the decay $$B\rightarrow D^*\ell \nu $$ B → D ∗ ℓ ν at nonzero recoil. Our analysis includes 15 MILC ensembles with $$N_f=2+1$$ N f = 2 + 1 flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $$a\approx 0.15$$ a ≈ 0.15 fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence b and c quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $$|V_{cb}|$$ | V cb | . We obtain $$\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}$$ V cb = ( 38.40 ± 0 . 68 th ± 0 . 34 exp ± 0 . 18 EM ) × 10 - 3 . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $$\chi ^2\text {/dof} = 126/84$$ χ 2 /dof = 126 / 84 , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $$R(D^*) = 0.265 \pm 0.013$$ R ( D ∗ ) = 0.265 ± 0.013 , which confirms the current tension between theory and experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Coarse-grained explicit-solvent molecular dynamics simulations of semidilute unentangled polyelectrolyte solutions

In this study, we present results from explicit-solvent coarse-grained molecular dynamics (MD) simulations of fully charged, salt-free, and unentangled polyelectrolytes in semidilute solutions. The inclusion of a polar solvent in the model allows for a more physical representation of these solutions at concentrations, where the assumptions of a continuum dielectric medium and screened hydrodynamics break down. The collective dynamic structure factor of polyelectrolytes, S(q, t), showed that at q > q*, where q* = 2π/ξ is the polyelectrolyte peak in the structure factor S(q) and ξ is the correlation length, the relaxation time obtained from fits to stretched exponential was $\tau$ KWW ~ q -3 , which describes unscreened Zimm-like dynamics. This is in contrast to implicit-solvent simulations using a Langevin thermostat where $\tau$ KWW ~ q -2 . At q < q*, a crossover region was observed that eventually transitions to another inflection point $\tau$ KWW ~ q -2 at length scales larger than ξ for both implicit- and explicit-solvent simulations. The simulation results were also compared to scaling predictions for correlation length, ξ ~ c$-½\atop{p}$, specific viscosity, η sp ~ c$½\atop{p}$, and diffusion coefficient, D ~ c$0\atop{p}$, where c p is the polyelectrolyte concentration. The scaling prediction for ξ holds; however, deviations from the predictions for η sp and D were observed for systems at higher c p , which are in qualitative agreements with recent experimental results. This study highlights the importance of explicit-solvent effects in molecular dynamics simulations, particularly in semidilute solutions, for a better understanding of polyelectrolyte solution behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Average-atom model with Siegert states

In plasmas, electronic states can be well-localized bound states or itinerant free states, or something in between. In self-consistent treatments of plasma electronic structure such as the average-atom model, all states must be accurately resolved in order to achieve a converged numerical solution. Furthermore, this is a challenging numerical and algorithmic problem in large part due to the continuum of free states which is relatively expensive and difficult to resolve accurately. Siegert states are an appealing alternative. They form a complete eigenbasis with a purely discrete spectrum while still being equivalent to a representation in terms of the usual bound states and free states. However, many of their properties are unintuitive, and it is not obvious that they are suitable for self-consistent plasma electronic structure calculations. Here it is demonstrated that Siegert states can be used to accurately solve an average-atom model and offer advantages over the traditional finite-difference approach, including a concrete physical picture of pressure ionization and continuum resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Sequential spectral line analysis for accurate density and temperature diagnosis of laboratory opacity measurements

The accuracy of iron opacity calculated in stellar interiors has been questioned since the discovery of the “solar problem” and the discrepancies between the measured and modeled iron opacity reported in 2015. Experimental opacity benchmarks require accurate temperature and density measurements, which were inferred by analyzing tracer magnesium spectra in those experiments. Could the observed discrepancy be explained by insufficient accuracy in the inferred temperature, density, and their uncertainties? Previous analyses may have yielded biased results due to three limitations: (1) simultaneous multi-line fitting, (2) approximations in line-shape models, and (3) exclusion of certain spectral lines due to insufficient background characterization. Notably, the first issue is a common concern for many inversion methods, including Bayesian inferences. We present a refined analysis method that overcomes these limitations, applied to three categories of iron opacity experiments (Anchor 1, 2, and 3). In particular, the sequential fitting method yields unbiased results with more realistic uncertainties by accounting for line inconsistencies in the parameter uncertainties. The average electron temperature and density values are 162 ± 6 eV and (7.0 ± 1.9) × 10 21 cm −3 for six Anchor 1 experiments, 189 ± 7 eV and (3.4 ± 0.3) × 10 22 cm −3 for 21 Anchor 2 experiments, and 201 ± 6 eV and (4.8 ± 1.1) × 10 22 cm −3 for nine Anchor 3 experiments. These results show ∼4% temperature and ∼20% density reproducibility over a decade, which also aligns with the inferred parameter uncertainties. In conclusion, the resulting temperature and density uncertainties lead to a quasi-continuum iron opacity variation of ±4%–7% for wavelengths below 9.5 Å, which is insufficient to explain the significant model-data discrepancies reported in 2015.

Absorption spectroscopy↗

DOE OSTI2025

Optical Activity of Chiral Excitons

Recent activity in the area of chiroptical phenomena has been focused on the connection between structural asymmetry, electron spin configuration and light/matter interactions in chiral semiconductors. In these systems, spin‐splitting phenomena emerge due to inversion symmetry breaking and the presence of extended electronic states, yet the connection to chiroptical phenomena is lacking. Here, in this study, we develop an analytical effective mass model of chiral excitons, parameterized by density functional theory. The model accounts for parity mixing of the band edge Bloch functions resulting from polar distortions, resulting in allowed magnetic dipole transitions. Through the study of a prototypical chiral 2D hybrid perovskite semiconductor, we show that circular dichroism of the chiral exciton and its interband continuum emerges from spin‐splitting via cross‐coupling of Rashba‐like and chiral/helical spin‐texture components. To demonstrate the generality of our approach, and as a counterpoint, we apply our model to describe chiroptical properties of three‐dimensionally confined excitons in perovskite nanocrystals that occur without chiral lattice distortions.

14 SOLAR ENERGY↗

DOE OSTI2021

Quantitative comparison of opacities calculated using the R -matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Delahaye, F.↗

DOE OSTI2023

Enhanced negative energy with a massless Dirac field

Motivated by traversable wormhole constructions that require large amounts of negative energy, we explore constraints on the amount of negative energy that can be carried by a free Dirac field in a slab-shaped region between two parallel spatial planes. Specifically, we ask what is the minimum possible uniform energy density that can exist at some time, considering all possible states and all possibilities for the physics outside the slab. The vacuum state where we identify the two sides of the slab with antiperiodic boundary conditions gives one possible state with uniform negative energy, but we argue that states with more negative energy exist above 1+1 dimensions. Technically, we reduce the problem to studying a massive Dirac field on an interval in 1+1 dimensions and numerically search for states with uniform energy density in a lattice regulated model. We succeed in finding states with enhanced negative energy (relative to the antiperiodic vacuum) which also appear to have a sensible continuum limit. Our results for the mass-dependence of the minimum uniform energy density in 1+1 dimensions suggest that for a 3+1 dimensional massless Dirac fermion, it is possible to have states with arbitrarily large uniform negative energy density in an arbitrarily wide slab.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗