Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interpolation”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Comparative collisionless alpha particle confinement in stellarator reactors with the XGC gyrokinetic code

Whole volume simulations of stellarators are necessary to address a number of important physics and engineering issues, including turbulent transport prediction and optimising divertor operation. In this work, advances in development of a whole volume stellarator capability for the global gyrokinetic particle-in-cell code XGC are detailed. A 3D interpolation of the equilibrium magnetic eld to the last close flux surface, calculated using the VMEC ideal magnetohydrodynamic equilibrum code, has been implemented. Orbit tracing with the XGC code is verified for the Wendelstein 7-X stellarator by comparison with the BEAMS3D and EUTERPE codes. The XGC stellarator version has then been used to consider collisionless alpha particle con nement in representative quasi-isodynamic and quasi-axisymmetric designs. It has been shown that, at least without further optimisation for energetic particle confinement, the Wendelstein 7-X-like quasi-isodynamic reactor design exhibits better alpha particle confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Dynamically weighted multireference perturbation theory: Combining the advantages of multi-state and state-averaged methods

In this paper, we introduce two new approaches to compute near-degenerate electronic states based on the driven similarity renormalization group (DSRG) framework. The first approach is a unitary multi-state formalism based on the DSRG (MS-DSRG), whereby an effective Hamiltonian is built from a set of state-specific solutions. The second approach employs a dynamic weighting parameter to smoothly interpolate between the multi-state and the state-averaged DSRG schemes. The resulting dynamically weighted DSRG (DW-DSRG) theory incorporates the most desirable features of both multi-state approaches (ability to accurately treat many states) and state-averaged methods (correct description of avoided crossings and conical intersections). We formulate second-order perturbation theories (PT2) based on the MS- and DW-DSRG and study the potential energy curves of LiF, the conical intersection of the two lowest singlet states of NH 3 , and several low-lying excited states of benzene, naphthalene, and anthracene. The DW-DSRG-PT2 predicts the correct avoided crossing of LiF and avoids artifacts produced by the corresponding state-specific and multi-state theories. Excitation energies of the acenes computed with the DW-DSRG-PT2 are found to be more accurate than the corresponding state-averaged values, showing a small dependence on the number of states computed.

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

MPRAD: A Monte Carlo and ray-tracing code for the proton radiography in high-energy-density plasma experiments

Proton radiography is used in various high-energy-density (HED) plasma experiments. In this paper, we describe a Monte Carlo and raytracing simulation tool called multimegaelectronvolt proton radiography (MPRAD) that can be used for modeling the deflection of proton beams in arbitrary three dimensional electromagnetic fields as well as the diffusion of the proton beams by Coulomb scattering and stopping power. Specifically, the Coulomb scattering and stopping power models in cold matter and fully ionized plasma are combined using interpolation. We discuss the application of MPRAD in a few setups relevant to HED plasma experiments where the plasma density can play a role in diffusing the proton beams and affecting the prediction and interpretation of the proton images. It is shown how the diffusion due to plasma density can affect the resolution and dynamical range of the proton radiography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Insulators at fractional fillings in twisted bilayer graphene partially aligned to hexagonal boron nitride

At partial fillings of its flat electronic bands, magic-angle twisted bilayer graphene (MATBG) hosts a rich variety of competing correlated phases that show sample-to-sample variations. Divergent phase diagrams in MATBG are often attributed to the sublattice polarization energy scale, tuned by the degree of alignment of the hexagonal boron nitride (hBN) substrates typically used in van der Waals devices. Unaligned MATBG exhibits unconventional superconductor and correlated insulator phases, while nearly perfectly aligned MATBG/hBN exhibits zero-field Chern insulating phases and lacks superconductivity. Here we use scanning tunneling microscopy and spectroscopy (STM/STS) to observe gapped phases at partial fillings of the flat bands of MATBG in a new intermediate regime of sublattice polarization, observed when MATBG is only partially aligned (θ Gr-hBN ≈ 1.65°) to the underlying hBN substrate. Under this condition, MATBG hosts not only phenomena that naturally interpolate between the two sublattice potential limits, but also unexpected gapped phases absent in either of these limits. At charge neutrality, we observe an insulating phase with a small energy gap (Δ < 5 meV) likely related to weak sublattice symmetry breaking from the hBN substrate. In addition, we observe new gapped phases near fractional fillings ν = ±1/3 and ν = ±1/6, which have not been previously observed in MATBG. Importantly, energy-resolved STS unambiguously identifies these fractional filling states to be of single-particle origin, possibly a result of the super-superlattice formed by two moiré superlattices. Finally, our observations emphasize the power of STS in distinguishing single-particle gapped phases from many-body gapped phases in situations that could be easily confused in electrical transport measurements, and demonstrate the use of substrate engineering for modifying the electronic structure of a moiré flat-band material.

Physics↗
DOE OSTI2023

Eulerian finite element simulations of the drop weight impact test with a dislocation Density-based continuum model

During a drop weight impact test, the kinetic energy of the falling weight is transferred into the sample resting on the anvil. The plastic deformation in the sample is an important mechanism for the dissipation of the input kinetic energy. We use Eulerian finite element analysis to simulate the deformation and temperature evolution in a 1,3,5-trinitro-1,3,5-triazine (RDX) sample consisting of multiple crystals. In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, i.e., advection. In an effort to reduce the advection errors, we use a rate form of a dislocation density-based continuum model by Luscher et al. Here, the simulations predict localization of plastic deformation, and plastic dissipation as a significant source of heat generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗
DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

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

Phase behavior of continuous-space systems: A supervised machine learning approach

The phase behavior of complex fluids is a challenging problem for molecular simulations. Supervised machine learning (ML) methods have shown potential for identifying the phase boundaries of lattice models. In this work, we extend these ML methods to continuous-space systems. We propose a convolutional neural network model that utilizes grid-interpolated coordinates of molecules as input data of ML and optimizes the search for phase transitions with different filter sizes. We test the method for the phase diagram of two off-lattice models, namely, the Widom–Rowlinson model and a symmetric freely jointed polymer blend, for which results are available from standard molecular simulations techniques. The ML results show good agreement with results of previous simulation studies with the added advantage that there is no critical slowing down. We find that understanding intermediate structures near a phase transition and including them in the training set is important to obtain the phase boundary near the critical point. The method is quite general and easy to implement and could find wide application to study the phase behavior of complex fluids.

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

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

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

A generalized forward fit for neutron detectors with energy-dependent response functions

To date, most analysis of neutron time-of-flight data from inertial confinement fusion experiments has focused on the relatively small range of energies corresponding to the primary neutrons from DD and DT fusion, and have therefore employed instrument response functions (IRF’s) corresponding to monoenergetic 2.45-MeV or 14.03-MeV neutrons. For analysis of time-of-flight signals corresponding to broader ranges of neutron energies, accurate treatment of the data requires the use of an energy-dependent IRF. Here, this work describes interpolation of the IRF for neutrons of arbitrary energy, construction of an energy-dependent IRF, and application of this IRF in a forward fit via matrix multiplication. As an example of the application of this method, an analysis of synthetic data relevant to TT fusion experiments at the Omega Laser Facility is discussed. This example is used to illustrate the differences between a forward fit that uses an energy-dependent IRF and a forward fit that uses a monoenergetic IRF. Use of the energy-dependent IRF is shown to result in accurate inference of the fit parameters of interest.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collisional excitation of interstellar PN by H 2 : New interaction potential and scattering calculations

Here, rotational excitation of interstellar PN molecules induced by collisions with H 2 is investigated. We present the first ab initio four-dimensional potential energy surface (PES) for the PN–H 2 van der Waals system. The PES was obtained using an explicitly correlated coupled cluster approach with single, double, and perturbative triple excitations [CCSD(T)-F12b]. The method of interpolating moving least squares was used to construct an analytical PES from these data. The equilibrium structure of the complex was found to be linear, with H 2 aligned at the N end of the PN molecule, at an intermolecular separation of 4.2 Å. The corresponding well-depth is 224.3 cm -1 . The dissociation energies were found to be 40.19 cm -1 and 75.05 cm -1 for complexes of PN with ortho-H 2 and para-H 2 , respectively. Integral cross sections for rotational excitation in PN–H 2 collisions were calculated using the new PES and were found to be strongly dependent on the rotational level of the H 2 molecule. These new collisional data will be crucial to improve the estimation of PN abundance in the interstellar medium from observational spectra.

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

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ice, glass, and solid phases in artificial spin systems with quenched disorder

We present a numerical study on a disordered artificial spin-ice system, which interpolates between the long-range ordered square ice and the fully degenerate shakti ice. Starting from the square-ice geometry, disorder is implemented by adding vertical/horizontal magnetic islands to the center of some randomly chosen square plaquettes of the array at different densities. When no island is added, we have ordered square ice. When all square plaquettes have been modified, we obtain shakti ice, which is disordered yet in a topological phase corresponding to the Rys F-model. In between, geometrical frustration due to these additional center spins disrupts the long-range Ising order of square ice, giving rise to a spin-glass regime at low temperatures. The artificial spin system proposed in our work provides an experimental platform to study the interplay between quenched disorder and geometrical frustration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sublimation temperature sensor for temperature locale

The sublimation temperature sensor (or “sublime sensor”) provides a continuum of measurement locations in which certain maximum temperatures are achieved during a heat up/cool down cycle. A predetermined material is encapsulated within a vacuum-sealed, non-volatile long tube (i.e., both ends capped and L$\gg$ D). This assembly is then inserted and centered into a heated zone, such as a furnace, exhaust pipe, or reactor. As the temperature increases, the material will sublimate (i.e., a process of having both the solid and gaseous states of matter simultaneously present) and will begin to fill the void - moving outward in both directions toward the ends of the tube. Once beyond the elevated temperatures, the gas will de-sublimate (i.e., deposition) onto the inner wall of the tube. The desired result of the sensor is the ring of material that develops over a relatively short period of time. This material deposit can be equated with temperature at an exact location. There is no need to interpolate and/or extrapolate for the desired measurement. Accuracy has been recorded for temperature locations on the range of ±2 mm over a 1 m span. Likewise, the precision of the measurement is ±0.2% the overall sensor domain. Furthermore, individual tubes with unique materials and pressures can be bundled together to provide a complete temperature profile of the heated zone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Influence of strain and dislocations on GaSb/GaAs quantum dots: From nested to staggered band alignment

We investigate the influence of strain and dislocations on band alignment in GaSb/GaAs quantum dot systems. Composition profiles from cross-sectional scanning tunneling microscopy images are interpolated onto a finite element mesh in order to calculate the distribution of local elastic strain, which is converted to a spatially varying band alignment using deformation potential theory. Our calculations predict that dislocation-induced strain relaxation and charging lead to significant local variations in band alignment. Furthermore, misfit strain induces a transition from a nested (type I) to a staggered (type II) band alignment. Although dislocation-induced strain relaxation prevents the type I to type II transition, electrostatic charging at dislocations induces the staggered band alignment once again.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Pressure, temperature, and orientation dependent thermal conductivity of pentaerythritol tetranitrate (PETN)

We use reverse non-equilibrium molecular dynamics to determine the thermal conductivity tensor, κ, of tetragonal pentaerythritol tetranitrate (PETN). The most stable form under standard ambient conditions (PETN I, with space group P4¯2 1 c) is considered in the temperature and pressure intervals 200–500 K and 0–8 GPa, which covers approximately the stability range for this polymorph. We compute κ along the [100] and [001] directions, which is enough to construct the full thermal conductivity tensor for a system with tetragonal symmetry. In addition, we also determine κ along the [010] direction and confirm that κ 100 ~ κ 010 , with less than 4% average absolute error between the two quantities. We observe an anisotropic response for κ, with κ 100 > κ 001 across the whole (T,P) interval, and 37% difference at 300 K, 0 GPa. Furthermore, we provide analytical functions to interpolate κ(P,T) within the fitting interval and serve as input for continuum-scale simulations.

36 MATERIALS SCIENCE↗