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

SDSS-IV MaStar: Stellar parameter determination with continuum-supplemented full-spectrum fitting

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent χ 2 fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total χ 2 for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The χ 2 distribution for a given region of model space was sampled using a Markov chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates (T eff , log g, [Fe/H], and [α/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC’s prescribed minimum-χ 2 result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the χ 2 distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Local conservation of energy in fully implicit PIC algorithms

We consider the issue of strict, fully discrete local energy conservation for a whole class of fully implicit local-charge- and global-energy-conserving particle-in-cell (PIC) algorithms. Earlier studies demonstrated these algorithms feature strict global energy conservation. However, whether a local energy conservation theorem exists (in which the local energy update is governed by a flux balance equation at every mesh cell) for these schemes is unclear. In this study, we show that a local energy conservation theorem indeed exists. We begin our analysis with the 1D electrostatic PIC model without orbit-averaging, and then generalize our conclusions to account for orbit averaging, multiple dimensions, and electromagnetic models (Darwin). In all cases, a temporally, spatially, and particle-discrete local energy conservation theorem is shown to exist, proving that these formulations (as originally proposed in the literature), in addition to being locally charge conserving and globally energy conserving, are strictly locally energy conserving as well. In contrast to earlier proofs of local conservation in the literature, which only considered continuum time, our result is valid for the fully implicit time-discrete version of all models considered, including important features such as orbit averaging. We demonstrate the local-energy-conservation property numerically with a paradigmatic numerical example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗
DOE OSTI2021

Towards modelling ghostly damped Ly α s

ABSTRACT We use simple models of the spatial structure of the quasar broad-line region (BLR) to investigate the properties of so-called ghostly damped Ly α (DLA) systems detected in Sloan Digital Sky Survey (SDSS) data. These absorbers are characterized by the presence of strong metal lines but no H i Ly α trough is seen in the quasar spectrum indicating that, although the region emitting the quasar continuum is covered by an absorbing cloud, the BLR is only partially covered. One of the models has a spherical geometry, another one is the combination of two wind flows, whereas the third model is a Keplerian disc. The models can reproduce the typical shape of the quasar Ly α emission and different ghostly configurations. We show that the DLA H i column density can be recovered precisely independently of the BLR model used. The size of the absorbing cloud and its distance to the centre of the AGN are correlated. However, it may be possible to disentangle the two using an independent estimate of the radius from the determination of the particle density. Comparison of the model outputs with SDSS data shows that the wind and disc models are more versatile than the spherical one and can be more easily adapted to the observations. For all the systems, we derive log N(H i)(cm−2) > 20.5. With higher quality data, it may be possible to distinguish between the models.

Laloux, Brivael↗
DOE OSTI2021

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct ab initio calculation of the 4 He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly diagonalizing the nuclear Hamiltonian in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained by directly diagonalizing the ab initio no-core shell model Hamiltonian is achievable numerically. Specifically, we calculate the 4 He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4 He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Astronomy & Astrophysics↗
DOE OSTI2020

Toward Engineering of Solution Microenvironments for the CO 2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional–Continuum Theory

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this study, we apply a recently developed density functional theory (DFT)–continuum approach based on the effective screening medium method and reference interaction site model (ESM–RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Furthermore, applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO 2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

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

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗
DOE OSTI2024

Absence of quantization in the circular photogalvanic effect in disordered chiral Weyl semimetals

The circularly polarized photogalvanic effect (CPGE) is studied in chiral Weyl semimetals with short-range quenched disorder. Without disorder, the topological properties of chiral Weyl semimetals lead to quantization of the CPGE, which is a second-order optical response. Furthermore, using a combination of diagrammatic perturbation theory in the continuum and exact numerical calculations via the kernel polynomial method on a lattice model, we show that disorder perturbatively destabilizes the quantization of the CPGE.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Average-atom calculations of bound-free and free-free cross sections in dense plasmas

Cross sections for photo-induced bound-free and free-free transitions in plasmas are evaluated in the average-atom approximation and applied to determine opacities of dense plasmas of light elements. Parameters characterizing the plasmas (chemical potential, average ionic charge, free electron density, bound and continuum wave functions, and occupation numbers) are obtained from the average-atom model. Lowest-order calculations of the free-free cross sections, which diverge in the low-frequency limit, are regularized by accounting for the finite electron-ion relaxation time. The resulting analysis provides the basis for average-atom studies of plasma opacities. Such studies are presented for dense lithium, beryllium, boron, and carbon. Applications are given to Rosseland mean opacities of dense hydrogen and deuterium plasmas and to comparisons of free-free to bound-free opacities in shock-compressed plasmas. Average-atom cross section and opacity calculations are extended to plasmas consisting of more than one ionic species, boron nitride, polystyrene, and a composite H, He, C plasma.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

SU(3) Topology of Magnon-Phonon Hybridization in 2D Antiferromagnets

Magnon-phonon hybrid excitations are studied theoretically in a two-dimensional antiferromagnet with an easy axis normal to the plane. In this work, we show that two magnon bands and one phonon band are intertwined by the magnetoelastic coupling through a nontrivial SU(3) topology, which can be intuitively perceived by identifying a skyrmion structure in the momentum space. Our results are insensitive to lattice details and generally applicable to two-dimensional antiferromagnets. We show this by developing a continuum theory as the long-wavelength approximation to the tight-binding model. The theoretical results can be probed by measuring the thermal Hall conductance as a function of the temperature and the magnetic field. We envision that the magnetoelastic coupling in antiferromagnets can be a promising venue in search of various topological excitations, which cannot be found in magnetic or elastic models alone.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum-Path-Resolved Attosecond High-Harmonic Spectroscopy

Strong-field ionization of molecules releases electrons which can be accelerated and driven back to recombine with their parent ion, emitting high-order harmonics. This ionization also initiates attosecond electronic and vibrational dynamics in the ion, evolving during the electron travel in the continuum. Revealing this subcycle dynamics from the emitted radiation usually requires advanced theoretical modeling. Here we show that this can be avoided by resolving the emission from two families of electronic quantum paths in the generation process. The corresponding electrons have the same kinetic energy, and thus the same structural sensitivity, but differ by the travel time between ionization and recombination-the pump-probe delay in this attosecond self-probing scheme. We measure the harmonic amplitude and phase in aligned CO 2 and N 2 molecules and observe a strong influence of laser-induced dynamics on two characteristic spectroscopic features: a shape resonance and multichannel interference. This quantum-path-resolved spectroscopy thus opens wide prospects for the investigation of ultrafast ionic dynamics, such as charge migration.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Functional Photoresists for Energy Applications. Final Report

Monolithic ultralow-density porous bulk materials have recently attracted much interest due to many emerging applications in the areas of catalysis, energy storage and conversion, and thermal insulation. They are also important components of high energy density (HED) and inertial confinement fusion (ICF) targets. However, despite tremendous progress that has been made in the synthesis of porous materials, deterministic and independent control over microscopic architecture, density and composition remain key issues, and their integration in high precision devices requires cost and time-intensive mechanical machining that not only reduces reproducibility by generating debris but also limits the complexity of the 3D shapes that can be realized. In this project, we overcame these limitations by developing a universal templating capability that provides deterministic and independent control over density, composition, architecture, and macroscopic sample shape. This was achieved by developing the technology to 1) 3D print ultrahigh resolution, ultra-high precision polymeric micro-lattice templates, 2) coat these templates with the desired materials, and 3) removing the template (Fig. 1a). Atomic layer deposition (ALD) provides the atomic scale coating thickness accuracy required for precisely controlling density. While this templating approach had been demonstrated in prior work, limitations in suitable photoresists, 3D print technologies, print design, and template removal techniques did not allow the fabrication of millimeter-sized high-precision parts with sub-micron resolution. To enable this technology, we developed 1) two-photon polymerization (TPP) print designs that enable the fabrication of millimeter-sized, mechanically robust polymeric templates with sub-micron resolution and 2) a continuum level TPP printing simulation capability for additional print design guidance; 3) atomistic models to study photoresist polymerization kinetics and network topography, 4) refractive index matched polymeric and preceramic TPP photoresists, and 5) functional TPP photoresists including porous voxel structures and self-immolative polymer photoresist chemistries; and 6) damage free template removal techniques that enable the fabrication of defect-free high-precision low-density foam components. We also developed a templating approach for pure carbon microlattices with a unique tube-in-tube ligament morphology. As a test platform, we pursued the fabrication of foam liners that promise to further increase the neutron yield in indirect drive ICF experiments by improving implosion symmetry control and coupling between the laser and the deuterium-tritium fuel. This application requires fabrication and integration of a ultra-high precision, millimeter-sized, thin-walled (200-400 micrometer thick), low-density (10-30 mg/cc), high atomic number (high Z) cylindrical foam tube into the gold hohlraum of an indirect drive ICF target (Fig. 1b). While our hohlraum liner test case will mainly find application in HED and ICF experiments, the underlying science will also directly apply to previously developed nanoparticle and additive manufacturing technologies and will advance those techniques as well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Gamow shell model with realistic interactions: a theoretical framework for ab initio nuclear structure at drip-lines

Abstract Ab initio approaches are among the most advanced models to solve the nuclear many-body problem. In particular, the no-core–shell model and many-body perturbation theory have been recently extended to the Gamow shell model framework, where the harmonic oscillator basis is replaced by a basis bearing bound, resonance and scattering states, i.e. the Berggren basis. As continuum coupling is included at basis level and as configuration mixing takes care of inter-nucleon correlations, halo and resonance nuclei can be properly described with the Gamow shell model. The development of the no-core Gamow shell model and the introduction of the Q ¯ ˆ -box method in the Gamow shell model, as well as their first ab initio applications, will be reviewed in this paper. Peculiarities compared to models using harmonic oscillator bases will be shortly described. The current power and limitations of ab initio Gamow shell model will also be discussed, as well as its potential for future applications.

Physics↗
DOE OSTI2025

Correlation effects in magic-angle twisted bilayer graphene: An auxiliary-field quantum Monte Carlo study

Magic-angle twisted bilayer graphene (MATBG) presents a fascinating platform for investigating the effects of electron interactions in topological flat bands. The Bistritzer-MacDonald (BM) model provides a simplified quantitative description of the flat bands. Introducing long-range Coulomb interactions leads to an interacting BM (IBM) Hamiltonian, a momentum-space continuum description which offers a very natural starting point for many-body studies of MATBG. Accurate and reliable many-body computations in the IBM model are challenging, however, and have been limited mostly to special fillings or smaller lattice sizes. We employ a state-of-the-art auxiliary-field quantum Monte Carlo (AFQMC) method to study the IBM model, which constrains the sign problem to enable accurate treatment of large system sizes. We determine ground-state properties and quantify errors compared to mean-field theory calculations. Our calculations identify correlated metal states and their competition with the insulating Kramers intervalley-coherent state at both half-filling and charge neutrality. Additionally, we investigate one- and three-quarter fillings, and examine the effect of many-body corrections beyond single Slater determinant solutions. We discuss the effect that details of the IBM Hamiltonian have on the results, including different forms of double-counting corrections, and the need to establish and precisely specify many-body Hamiltonians to allow more direct and quantitative comparisons with experiments in MATBG. Published by the American Physical Society 2025

Xiao, Zhi-Yu (ORCID:0000000219531579)↗
DOE OSTI2023

Quantum spin ice response to a magnetic field in the dipole-octupole pyrochlore Ce 2 Zr 2 O 7

The pyrochlore magnet Ce 2 Zr 2 O 7 has attracted much attention as a quantum spin ice candidate whose novelty derives in part from the dipolar-octupolar nature of the Ce 3+ pseudospin-1/2 degrees of freedom it possesses. We report heat capacity measurements on single crystal samples of Ce 2 Zr 2 O 7 down to T~0.1K in a magnetic field along the [$1,\bar{1},0$] direction. These measurements show that the broad hump in the zero-field heat capacity moves higher in temperature with increasing field strength and is split into two separate humps by the [$1,\bar{1},0$] magnetic field at ~2T. These separate features are due to the decomposition of the pyrochlore lattice into effectively decoupled chains for fields in this direction: One set of chains (α chains) is polarized by the field while the other (β chains) remains free. This situation is similar to that observed in the classical spin ices Ho 2 Ti 2 O 7 and Dy 2 Ti 2 O 7 , but with the twist that here the strong transverse exchange interactions produce substantial quantum effects. Our theoretical modeling suggests that the β chains are close to a critical state, with nearly-gapless excitations. We also report elastic and inelastic neutron scattering measurements on single crystal Ce 2 Zr 2 O 7 in [$1,\bar{1},0$] and [0,0,1] magnetic fields at temperatures down to T=0.03K. The elastic scattering behaves consistently with the formation of independent chains for a [$1,\bar{1},0$] field, while the [0,0,1] field produces a single field-induced elastic magnetic Bragg peak at (0,2,0) and equivalent wavevectors, indicating a polarized spin ice state for fields above ~3T. For both [$1,\bar{1},0$] and [0,0,1] magnetic fields, our inelastic neutron scattering results show an approximately dispersionless continuum of scattering that increases in both energy and intensity with increasing field strength. By modeling the complete set of experimental data using numerical linked cluster and semiclassical molecular dynamics calculations, we demonstrate the dominantly multipolar nature of the exchange interactions in Ce 2 Zr 2 O 7 and the smallness of the parameter θ, which controls the mixing between dipolar and octupolar degrees of freedom. Finally, these results support previous estimates of the microscopic exchange parameters and place strong constraints on the theoretical description of this prominent spin ice candidate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗