Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Accretion of Gas Giants Constrained by the Tidal Barrier

After protoplanets have acquired sufficient mass to open partial gaps in their natal protostellar disks, residual gas continues to diffuse onto horseshoe streamlines under the effect of viscous dissipation, and to meander in and out of the planets' Hill sphere. Within the Hill sphere, the horseshoe streamlines intercept gas flow in circumplanetary disks. The host star's tidal perturbation induces a barrier across the converging streamlines' interface. Viscous transfer of angular momentum across this tidal barrier determines the rate of mass diffusion from the horseshoe streamlines onto the circumplanetary disks, and eventually the accretion rate onto the protoplanets. Here, we carry out a series of numerical simulations to test the influence of this tidal barrier on superthermal planets. In weakly viscous disks, the protoplanets' accretion rate steeply decreases with their masses above the thermal limit. As their growth timescale exceeds the gas depletion timescale, their masses reach asymptotic values comparable to that of Jupiter. In relatively thick and strongly viscous disks, protoplanets' asymptotic masses exceed several times that of Jupiter. Two-dimensional numerical simulations show that such massive protoplanets strongly excite the eccentricity of nearby horseshoe streamlines, destabilize orderly flow, substantially enhance the diffusion rate across the tidal barrier, and elevate their growth rate until their natal disk is severely depleted. In contrast, eccentric streamlines remain stable in three-dimensional simulations. Based on the upper falloff in the observed mass distribution of known exoplanets, we suggest that their natal disks had relatively low viscosity (α ~ 10 –3 ) and modest thickness (H/R ~ 0.03–0.05).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Optical spectroscopy and modeling of uranium gas-phase oxidation: Progress and perspectives

We report that studies related to U gas-phase oxidation through plasma- and thermo-chemistry are important for many fields, including environmental monitoring, forensic analysis, debris analysis in a weapon detonation event, and nucleation physics. Recently, significant efforts have been made to understand the chemical pathways involved in the progression from U atoms to diatoms (UO) and polyatomic molecules (UxOy), employing optical spectroscopy tools and computational modeling. In many studies, laser ablation of U or a U-containing flow reactor are used as a highly resource-efficient, repeatable, tunable, and lab-scale testbed for studying gas-phase oxidation in U plasmas. The spectroscopic analysis of high-temperature gas-phase oxidation of U is challenging due to the congested U spectra, resolution limitations of instrumentation, and the numerous chemical reaction pathways possible. This article focuses on the current understanding and challenges related to studying U plasma chemistry, specifically U gas-phase oxidation and molecular formation, via optical spectroscopy of plasmas and associated computational and spectral modeling. The physical and chemical processes involved in the evolution from U atoms to U oxide molecules to nanoparticles and agglomerates (i.e., debris) are discussed in the context of optical spectroscopic studies. The article concludes by highlighting opportunities for future research efforts based on existing knowledge published in the literature.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

The Opus radiative opacity code, a new implementation of the super-transition array method

This report describes the development of our new radiative opacity code Opus, submitted as part of a 2020 Level 2 milestone requirement. This code aims to incorporate the capabilities in existing but separate opacity codes previously developed at LLNL (Opal, Topaz, and Tycho). The first phase in the development, which was to implement the super-transition (STA) method of Bar-Shalom and co-workers is now completed. In this work, we describe the underlying physics, implementation, usage, and regression testing as well as providing numerical examples for the opacity tables of boron, chromium, iron, and nickel. The latter are benchmarked against the LLNL legacy STA code Tycho, with an average absolute relative difference of ~1.5 and 0.5% for Planck and Rosseland opacity means, respectively. The cases where differences are larger are well understood in terms of numerical instabilities near discontinuities of the average degree of ionization inherent to the assumed opacity/equation-of-state model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application

In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multiscale plasticity of geomaterials predicted via constrained optimization‐based granular micromechanics

Abstract A general framework to derive nonlinear elastic and elastoplastic material models from granular micromechanics is proposed, where a constraint‐based variational structure is introduced to classical grain contact‐based homogenization methods of hyperelasticity. Like the classical hyperelastic methods, reference solutions for closed‐form hyperelastic material models are analytically derived from the grain‐scale contact mechanics. However, unlike prior methods, the proposed homogenization framework defines closed‐form hyperelastoplastic material models that extend multiscale variational methods to granular plasticity. The proposed framework is used to develop novel granular micromechanics‐based macroscopic models for a Mises type solid, Drucker–Prager type plasticity, and grain‐contact cohesive‐debonding with a deviatorically and volumetrically coupled nonlinearly elastic response. Macroscopic plastic parameters and yield criteria are explicitly related to their microscale counterparts, for example, the friction coefficient governing intergranular slip. Numerical examples and comparison to measurements from the literature, including triaxial compaction of concrete, are provided to investigate model predictions and demonstrate calibration to experimental data.

Bryant, E. C.↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
DOE OSTI2023

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Soil microbial community response to corrinoids is shaped by a natural reservoir of vitamin B12

Abstract Soil microbial communities perform critical ecosystem services through the collective metabolic activities of numerous individual organisms. Most microbes use corrinoids, a structurally diverse family of cofactors related to vitamin B12. Corrinoid structure influences the growth of individual microbes, yet how these growth responses scale to the community level remains unknown. Analysis of metagenome-assembled genomes suggests that corrinoids are supplied to the community by members of the archaeal and bacterial phyla Thermoproteota, Actinobacteria, and Proteobacteria. Corrinoids were found largely adhered to the soil matrix in a grassland soil, at levels exceeding those required by cultured bacteria. Enrichment cultures and soil microcosms seeded with different corrinoids showed distinct shifts in bacterial community composition, supporting the hypothesis that corrinoid structure can shape communities. Environmental context influenced both community- and taxon-specific responses to specific corrinoids. These results implicate corrinoids as key determinants of soil microbiome structure and suggest that environmental micronutrient reservoirs promote community stability.

58 GEOSCIENCES↗
DOE OSTI2026

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Emergency Voltage Regulation in Power Systems via Ripple-Type Control

With increasing penetrations of volatile renewable generation and cyber-physical disruptions, ensuring the safe operation of bulk power systems has become unprecedentedly challenging. Because communication and computational costs restrict centralized system dispatch to being called upon every few minutes, and because purely local schemes are shown to be insufficient, distributed controls have been advocated for handling unanticipated system conditions in real time. The applicability of distributed control schemes, however, is fundamentally limited by their need for widespread communication and model cognizance. In this context, we put forth a hybrid, low-communication, saturation-driven protocol for the coordination of control agents that are distributed over a physical system and are allowed to communicate with peers over a "hotline" communication network. Under this protocol, when agents observe a constraint violation based on local measurements, they respond locally until their control resources saturate, in which case they send a beacon for assistance to peer agents. The scheme ensures that minor violations are efficiently mitigated via fast local controls, whereas severe violations can be handled by collaboration among a relatively small set of agents. We evaluate the performance of this scheme via numerical tests on the IEEE 14-bus test feeder, where agents act upon noisy measurements under diverse scenarios of load variations and severe low-/high-voltage events.

collaboration↗
DOE OSTI2022

Numerical simulations of noisy quantum circuits for computational chemistry

The opportunities afforded by near-term quantum computers to calculate the ground-state properties of small molecules depend on the structure of the computational ansatz as well as the errors induced by device noise. Here we investigate the behavior of these noisy quantum circuits using numerical simulations to estimate the accuracy and fidelity of the prepared quantum states relative to the ground truth obtained by conventional means. We implement several different types of ansatz circuits derived from unitary coupled cluster theory for the purposes of estimating the ground-state energy of sodium hydride using the variational quantum eigensolver algorithm. We show how relative error in the energy and the fidelity scale with the levels of gate-based noise, the internuclear configuration, the ansatz circuit depth, and the parameter optimization methods.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Compare low-gain 1D and 3D models toward a storage ring-based X-ray free-electron-laser oscillator

Fourth-generation diffraction-limited synchrotron light sources (DL-SLS) present a compelling platform for advanced light source applications due to their ultra-low emittance (around tens of picometers), sub-micrometer beam stability, and high repetition rates (greater than 1 MHz). Despite these advantages, fundamental limitations, such as long bunch durations (tens of picoseconds), large relative energy spread (~10⁻³), and low peak current (<300 A), have hindered their effectiveness in free-electron laser (FEL) applications. Using a low-gain three-dimensional (3D) model based on NSLS-II upgrade lattice parameters, we demonstrate that a 6% single-pass gain is achievable with an optimized transverse gradient undulator (TGU). Concurrently, recent advances in X-ray optics—specifically, a four-mirror bow-tie resonator utilizing diamond mirrors and beryllium compound refractive lenses (CRLs)—enable total round-trip losses below 4% and an output coupling efficiency near 1%, rendering a Storage Ring based X-ray FEL oscillator (XFELO) a realistic possibility. In the 3D framework, the TGU effectively mitigates the detrimental effects of large energy spread by spatially dispersing the beam, leading to enhanced gain relative to conventional one-dimensional (1D) models. We establish both analytical and numerical connections between the 1D and 3D descriptions, leveraging GENESIS simulations to map equivalent gain conditions. While 1D resonance is defined in the energy domain, the 3D TGU system shifts the resonance dependence to transverse gradient and dispersion parameters. Once optimized, this framework allows predictive modeling using 1D formulas, which can then be validated through full 3D simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Technical Evaluation of the Department of Energy Research and Development Activities in Underground Research Laboratories - 20287

Congress created the U.S. Nuclear Waste Technical Review Board (NWTRB) in the 1987 Nuclear Waste Policy Amendments Act (Public Law 100-203) to evaluate the technical and scientific validity of activities undertaken by the Secretary of Energy to implement the Nuclear Waste Policy Act. Since 2012, DOE has collaborated in research conducted in several underground research laboratories (URLs) located in Europe and Asia. According to DOE, these international collaborations have been beneficial to its spent nuclear fuel (SNF) and high-level waste radioactive (HLW) disposal research program, particularly after the termination of the Yucca Mountain repository program when DOE began generic research on alternative host rocks (crystalline, clay, and salt) and repository environments very different from those at Yucca Mountain. In accordance with its mandate, the NWTRB is reviewing the DOE research and development (R and D) activities related to URLs. The NWTRB's review is documented in a report to Congress and the Secretary of Energy that will be released in January 2020, and this paper summarizes the NWTRB review and findings. The NWTRB held a fact-finding meeting with DOE and subsequently held a workshop on international URL collaborations in April of 2019. Based on the presentations and discussions at the workshop and at the fact-finding meeting, as well as information from reports published by DOE and others, there are four principal findings related to DoE's URL-related R and D activities. First, DOE participation in URL-related international research greatly benefits the U.S. geologic disposal R and D program by furthering its understanding of generic and site-specific disposal issues relevant to alternative repository host rocks and environments. DOE-funded R and D activities also are benefiting the URL-related research of other countries, especially in the area of complex analytical and numerical model/software development. Second, the more developed repository programs in other countries have focused on creating and strengthening their safety cases and making them transparent to the public. Repository programs in other countries use URLs to explain the technical bases underlying their safety cases, periodically reassess knowledge gaps and define new activities to strengthen the technical bases, and demonstrate the technology that will allow implementation of the proposed safety concept. Third, countries with more developed geologic disposal programs have found domestic URLs essential to their repository programs. DOE needs domestic URLs to advance geologic disposal efforts over the next decades and further its ability to train the next generation of scientists, engineers, and skilled technical workers. Fourth, DoE's international URL collaborations have advanced its generic disposal R and D program, including development of modeling capabilities recognized internationally as state-of-the-art, but further work on its coupled thermal-hydrological-mechanical-chemical models and URL- and laboratory-based research can strengthen its program. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

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

Influence of Nanoconfinement and Elevated Temperatures on Geocolloid-Facilitated Transport of Energy-related Contaminants

The main objective of this project is to obtain a better understanding of the transport behavior and interactions of geocolloids in the presence of energy- related contaminants under bulk and nanoconfined conditions. In that, we first established solid protocols to synthesize and fabricate numerous types of geocolloids with a. various sizes ranging from 40 – 800 nm and b. different hydrophile-lipophile balance (HLB) ratios in the range of 25:75 to 75:25. Direct force measurements with geocolloids having different degrees of surface coverages (i.e. HLB to realistically mimic adsorption of energy related contaminants) were conducted using the Surface Forces Apparatus (SFA) over a distance regime starting from 8μm all the way down to molecular contact. Repulsive forces were observed on approach starting from > 3μm, followed by an exponential increase of which magnitude appears to be larger than a decay length obtained from Derjaguin–Landau–Verwey–Overbeek (DLVO) theory in pure water. When the geocolloids were confined in salted water, the magnitude of onset of repulsion was varied as a function of salinity in solution, which can significantly alter a purely repulsive screened electrostatic (coulombic) interaction arising from, among geocolloids as well as between geocolloids and geosurfaces. The viscosity and flow characteristics of geocolloidal suspensions at different degrees of confinement were also investigated where we identified highly discontinuous rheological behaviors below a critical nanoconfinement level. We anticipate that the knowledge gained through this study will enable the scientists and researchers to better assess transport and fate behaviors of geocolloidal dispersions that can carry energy-related contaminants under realistically emulated geosystem.

03 NATURAL GAS↗