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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Fractons with twisted boundary conditions and their symmetries

In this work, we study several exotic systems, including the X-cube model, on a flat three-torus with a twist in the $xy$-plane. The ground state degeneracy turns out to be a sensitive function of various geometrical parameters. Starting from a lattice, depending on how we take the continuum limit, we find different values of the ground state degeneracy. Yet, there is a natural continuum limit with a well-defined (though infinite) value of that degeneracy. We also uncover a surprising global symmetry in $2+1$ and $3+1$ dimensional systems. It originates from the underlying subsystem symmetry, but the way it is realized depends on the twist. In particular, in a preferred coordinate frame, the modular parameter of the twisted two-torus $\tau = \tau_1 + i \tau_2$ has rational $\tau_1 = k / m$. Then, in systems based on $U(1)\times U(1)$ subsystem symmetries, such as momentum and winding symmetries or electric and magnetic symmetries, the new symmetry is a projectively realized $\mathbb{Z}_m\times \mathbb{Z}_m$, which leads to an $m$-fold ground state degeneracy. In systems based on $\mathbb{Z}_N$ symmetries, like the X-cube model, each of these two $\mathbb{Z}_m$ factors is replaced by $\mathbb{Z}_{\gcd(N,m)}$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A continuum damage approach to spallation and the role of microinertia

Spall failure is of interest due to its prevalence in high strain rate problems in which the spallation is driven by the interaction of release waves. In this article, a porosity-based damage model that includes microinertial effects is used to examine spall failure. The model is successfully calibrated to plate impact-driven experiments and then used to evaluate experimental conditions producing more extreme strain rate conditions, such as those in laser-driven experiments. The incorporation of microinertia allows us to better understand the increase in apparent macroscopic spall strength seen at high strain rates. Correspondingly, we conclude that the incorporation of microinertial effects improves the model’s predictive capabilities. Microinertial effects result in more severe local tensile stresses that affect the damage evolution, and microinertia can play a significant role in the approach to the material’s ideal strength at extreme loading rates. Additionally, a preliminary parametric study is also carried out to investigate the role of microstructural aspects such as nucleation volume fraction and initial pore radius. One counter-intuitive result from the microinertial effects is that, for a given nucleation site volume fraction, having larger initial pore nucleation sites can lead to an improved macroscopic spall strength.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The gallium solar neutrino capture cross section revisited

Solar neutrino flux constraints from the legacy GALLEX/GNO and SAGE experiments continue to influence contemporary global analyses of neutrino properties. The constraints depend on the neutrino absorption cross sections for various solar sources. Following recent work updating the 51 Cr and 37 Ar neutrino source cross sections, we reevaluate the 71 Ga solar neutrino cross sections, focusing on contributions from transitions to 71 Ge excited states, but also revising the ground-state transition to take into account new 71 Ge electron-capture lifetime measurements and various theory corrections. The excited-state contributions have been traditionally taken from forward-angle (𝑝, 𝑛) cross sections. Here we correct this procedure for the ≈ 10%–20% tensor operator contribution that alters the relationship between Gamow-Teller and (𝑝, 𝑛) transition strengths. Using state-of-the-art nuclear shell-model calculations to evaluate this correction, we find that it lowers the 8 B and hep neutrino cross sections. However, the addition of other corrections, including contributions from near-threshold continuum states that radiatively decay, leads to an overall increase in the 8 B and hep cross sections of ≈ 10% relative to the values recommended by Bahcall. Uncertainties are propagated using Monte Carlo simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Predictive scale-bridging simulations through active learning

Throughout computational science, there is a growing need to utilize the continual improvements in raw computational horsepower to achieve greater physical fidelity through scale-bridging over brute-force increases in the number of mesh elements. For instance, quantitative predictions of transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, are impossible without accounting for molecular-level interactions. Similarly, inertial confinement fusion simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions. With these two disparate applications in mind, we develop a novel capability which uses an active learning approach to optimize the use of local fine-scale simulations for informing coarse-scale hydrodynamics. Our approach addresses three challenges: forecasting continuum coarse-scale trajectory to speculatively execute new fine-scale molecular dynamics calculations, dynamically updating coarse-scale from fine-scale calculations, and quantifying uncertainty in neural network models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2023

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2021

Extreme Molecular Complexity Resulting in a Continuum of Carbonaceous Species in Biomass Burning Tar Balls from Wildfire Smoke

Biomass burning emits a wide range of carbonaceous particles into the atmosphere and has negative impacts on human health and the Earth’s radiative balance. Nonvolatile spherical organic aerosol particles, commonly known as tar balls, represent one of the most abundant particles in aged biomass burning smoke. However, the detailed molecular-level composition of ambient tar balls is largely unknown but critical to access their environmental impacts. Ambient aerosol samples collected during a wildfire event, which were ~90% tar balls by number fraction, were analyzed using ultrahigh resolution Orbitrap Elite mass spectrometry with four complementary ionization modes. Our results show the molecular composition of tar balls to be complex, composed of over 10,000 molecular formulas. Model estimated saturation mass concentrations and relative humidity dependent glass transition temperatures were consistent with low volatility and solid morphology as expected for tar balls. Room-temperature evaporation kinetics showed that these particles were composed of ~90% non-volatile species. The molecular complexity detected here signifies a continuum of carbonaceous species, ranging from C3 to C45 with continuous ranges of oxygenation and hydrogen saturation for each Cn. Approximately 24% of molecular formulas were estimated the be highly aromatic, which could indicate chemical compounds with negative health effects, and which may contribute to visible light absorption. The carbon continuum observed here has significant implications for the molecular characterization of atmospheric organic matter. The level of complexity detected here should not be ignored in future studies and we demonstrate that multiple analytical methods may be required to suitably interpret this complexity on a molecular level.

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

Improving evapotranspiration computation with electrical resistivity tomography in a maize field

Abstract Hydrogeophysical methods have been increasingly used to study subsurface soil–water dynamics, yet their application beyond the soil compartment or the quantitative link to soil hydraulic properties remains limited. To examine how these methods can inform model‐based evapotranspiration (ET) calculation under varying soil water conditions, we conducted a pilot‐scale field study at an experimental maize plot with manipulated irrigation treatments. Our goal was to develop a workflow for (1) acquiring and inverting field electrical resistivity tomography (ERT) data, (2) correlating ERT to soil hydraulic properties, (3) spatially characterizing soil water stress that feeds into ET modeling (the FAO‐56 model), and (4) evaluating the performance of ERT‐based ET computation. Our results showed that ERT was able to capture decimeter‐scale soil water content (SWC) dynamics from root water uptake and irrigation manipulation and the contrast of soil water stress between deficiently and fully irrigated maize. We also demonstrated the flexibility of using ERT to spatially integrate soil water stress in the soil volume of interest, which could be adjusted based on different crops and plot layouts. The integration of the ERT datasets into ET modeling provided insights into the spatial heterogeneity of the subsurface that has been challenging for point‐based sensing, which can further our understanding of the hydraulic dynamics in the soil‐plant‐atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reviews and syntheses: Ironing out wrinkles in the soil phosphorus cycling paradigm

Soil phosphorus (P) management remains a critical challenge for agriculture worldwide, and yet we are still unable to predict soil P dynamics as confidently as that of carbon (C) or nitrogen (N). This is due to both the complexity of inorganic P (P i ) and organic P (P o ) cycling and the methodological constraints that have limited our ability to trace P dynamics in the soil–plant system. In this review, we describe the challenges for building parsimonious, accurate, and useful biogeochemical models that represent P dynamics and explore the potential of new techniques to usher P biogeochemistry research and modeling forward. We conclude that research efforts should focus on the following: (1) updating the McGill and Cole (1981) model of P o mineralization by clarifying the role and prevalence of biochemical and biological P o mineralization, which we suggest are not mutually exclusive and may co-occur along a continuum of P o substrate stoichiometry; (2) further understanding the dynamics of phytate, a six C compound that can regulate the poorly understood stoichiometry of soil P; (3) exploring the effects of C and P o saturation on P sorption and P o mineralization; and (4) resolving discrepancies between hypotheses about P cycling and the methods used to test these hypotheses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Status of Lattice QCD Calculations of One and Two-Nucleon Matrix Elements

Lattice QCD calculations of basic hadronic quantities are now often performed with systematic control over the physical pion mass, continuum and infinite volume limits. This progress signals a new era in which lattice QCD will be used to compute basic properties of hadrons, nucleons and light nuclei directly from the Standard Model (SM). In particular, we will be able to use lattice QCD, combined with Effective Field Theory, to make quantitative statements about the interaction of SM matter and potential beyond the SM physics ranging from direct dark matter detection to permanent electric dipole moments. I will briefly review the status and challenges of lattice QCD calculations of select one- and two-nucleon matrix elements. These calculations are significantly more challenging than the basic quantities computed with all systematics controlled, but substantial progress is being made.

Walker-Loud, Andre↗
DOE OSTI2021

The O K -2 V spectrum of CO: the influence of the second core-hole

Using synchrotron radiation in the tender X-ray regime, a photoelectron spectrum showing the formation of single site double-core-hole pre-edge states, involving the K shell of the O atom in CO, has been recorded by means of high-resolution electron spectroscopy. The experimentally observed structures have been simulated, interpreted and assigned, employing state-of-the-art ab initio quantum chemical calculations, on the basis of a theoretical model, accounting for their so-called direct or conjugate character. Features appearing above the double ionization threshold have been reproduced by taking into account the strong mixing between multi-excited and continuum states. The shift of the σ* resonance below the double ionization threshold, in combination with the non-negligible contributions of multi-excited configurations in the final states reached, gives rise to a series of avoided crossings between the different potential energy curves.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Hydrogen reionization ends by z = 5.3: Lyman-α optical depth measured by the XQR-30 sample

ABSTRACT The presence of excess scatter in the Ly-α forest at z ∼ 5.5, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionization. However, low data quality and systematic uncertainties complicate the use of Ly-α transmission as a precision probe of reionization’s end stages. In this paper, we assemble a sample of 67 quasar sightlines at z > 5.5 with high signal-to-noise ratios of >10 per ≤15 km s−1 spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR > 20 per pixel) spectra of 30 quasars at z > 5.7. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-α systems. We present improved measurements of the mean Ly-α transmission over 4.9 < z < 6.1. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-α transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to z ≤ 5.2 (<1σ), and mild tension (∼2.5σ) at z = 5.3. Homogeneous UVB models are ruled out by excess Ly-α transmission scatter at z ≥ 5.4 with high confidence (>3.5σ). Our results indicate that reionization-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least z = 5.3 (t = 1.1 Gyr after the big bang). This is further evidence for a late end to reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Connecting Low- and High-redshift Weak Emission-line Quasars via Hubble Space Telescope Spectroscopy of Lyα Emission

We present ultraviolet spectroscopy covering the Lyα + N v complex of six candidate low-redshift (0.9 < z < 1.5) weak emission-line quasars (WLQs) based on observations with the Hubble Space Telescope. The original systematic searches for these puzzling Type 1 quasars with intrinsically weak broad emission lines revealed an N ≈ 100 WLQ population from optical spectroscopy of high-redshift (z > 3) quasars, defined by a Lyα + N v rest-frame equivalent width (EW) threshold <15.4 Å. Identification of lower-redshift (z < 3) WLQ candidates, however, has relied primarily on optical spectroscopy of weak broad emission lines at longer rest-frame wavelengths. With these new observations expanding existing optical coverage into the ultraviolet, we explore unifying the low- and high-z WLQ populations via EW[Lyα+N v]. Two objects in the sample unify with high-z WLQs, three others appear consistent with the intermediate portion of the population connecting WLQs and normal quasars, and the final object is consistent with typical quasars. The expanded wavelength coverage improves the number of available line diagnostics for our individual targets, allowing a better understanding of the shapes of their ionizing continua. The ratio of EW[Lyα+N v] to EW[Mg ii] in our sample is generally small but varied, favoring a soft ionizing continuum scenario for WLQs, and we find a lack of correlation between EW[Lyα+N v] and the X-ray properties of our targets, consistent with a "slim-disk" shielding gas model. We also find indications that weak absorption may be a more significant contaminant in low-z WLQ populations than previously thought.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Asymptotic consistency of the WSINDy algorithm in the limit of continuum data

In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.

asymptotic consistency↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗
DOE OSTI2020

Mobile orbitons in Ca 2 CuO 3 : Crucial role of Hund's exchange

We investigate the Cu L 3 edge resonant inelastic x-ray scattering (RIXS) spectra of a quasi-1D antiferromagnet Ca 2 CuO 3 . In addition to the magnetic excitations, which are well-described by the two-spinon continuum, we observe two dispersive orbital excitations, the 3d xy and the 3d yz orbitons. We carry out a quantitative comparison of the RIXS spectra, obtained with two distinct incident polarizations, with a theoretical model. We show that any realistic spin-orbital model needs to include a finite, but realistic, Hund's exchange J H ≈ 0.5 eV . Its main effect is an increase in orbiton velocities, so that their theoretically calculated values match those observed experimentally. Even though Hund's exchange also mediates some interaction between spinon and orbiton, the picture of spin-orbit separation remains intact and describes orbiton motion in this compound.

1-dimensional spin chains↗
DOE OSTI2025

Possible Observation of Quadrupole Waves in Spin Nematics

Discovery of new states of matter is a key objective in modern condensed matter physics, which often leads to revolutionary technological advancements such as superconductivity. Quantum spin nematic, a “hidden order” that evades conventional magnetic probes, is one such state. Na 2 ⁢BaNi⁢(PO 4 ) 2 is a potential spin nematic material, suggested by the observation of a two-magnon Bose-Einstein condensation from above the saturation field. However, direct confirmation of the spin nematicity remains elusive. Here, this Letter presents inelastic neutron scattering spectra from the putative spin nematic phases of Na 2⁢ BaNi⁢(PO 4 ) 2 , revealing low-energy quadrupole waves that are absent in the neighboring conventional magnetic phases. A spin-one model quantitatively captures the full details of the spin excitation spectra across all low-temperature phases, providing direct evidence of the spin nematic orders. Additionally, we show evidence of the three-magnon continuum and two-magnon bound states in the 1/3-magnetization plateau, revealing condensation of the two-magnon bound state as the origin of the low-field spin nematic supersolid phase.

Sheng, Jieming [Southern University of Science and↗