Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximation methods”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Optimization of the Jastrow factor using the random-phase approximation and a similarity-transformed Hamiltonian: Application to band-structure calculation for some semiconductors and insulators

Based on the random-phase approximation and the transcorrelated (TC) method, we optimize the Jastrow factor together with one-electron orbitals in the Slater determinant in the correlated wave function with a new scheme for periodic systems. The TC method is one of the promising wave function theories for first-principles electronic structure calculation, where the many-body wave function is approximated as a product of a Slater determinant and a Jastrow factor, and the Hamiltonian is similarity-transformed by the Jastrow factor. Using this similarity-transformed Hamiltonian, we can optimize the one-electron orbitals without evaluating 3N-dimensional integrations for the N-electron system. In contrast, optimization of the Jastrow factor within the framework of the TC method is computationally much more expensive and has not been performed for solid-state calculations before. In this study, we also benefit from the similarity-transformation in optimizing the Jastrow factor. Our optimization scheme is tested in applications to some solids from narrow-gap semiconductors to wide-gap insulators, and it is verified that the band gap of a wide-gap insulator and the lattice constants of some solids are improved by this optimization with reasonable computational cost.

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

Numeric kinetic energy operators for molecules in polyspherical coordinates

Generalized curvilinear coordinates, as, e.g., polyspherical coordinates, are in general better adapted to the resolution of the nuclear Schroedinger equation than rectilinear ones like the normal mode coordinates. However, analytical expressions of the kinetic energy operators (KEOs) for molecular systems in polyspherical coordinates may be prohibitively complicated for large systems. In this paper we propose a method to generate a KEO numerically and bring it to a form practicable for dynamical calculations. To examine the new method we calculated vibrational spectra and eigenenergies for nitrous acid (HONO) and compare it with results obtained with an exact analytical KEO derived previously [F. Richter, P. Rosmus, F. Gatti, and H.-D. Meyer, J. Chem. Phys. 120, 6072 (2004)]. In a second example we calculated {pi}{yields}{pi}* photoabsorption spectrum and eigenenergies of ethene (C{sub 2}H{sub 4}) and compared it with previous work [M. R. Brill, F. Gatti, D. Lauvergnat, and H.-D. Meyer, Chem. Phys. 338, 186 (2007)]. In this ethene study the dimensionality was reduced from 12 to 6 by freezing six internal coordinates. Results for both molecules show that the proposed method for obtaining an approximate KEO is reliable for dynamical calculations. The error in eigenenergies was found to be below 1 cm{sup -1} for most states calculated.

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

A new adaptive modeling strategy based on optimal control for atomic-to-continuum coupling simulations

A new adaptive method based on an optimal control approach is proposed for adaptive modeling of an atomic-to-continuum coupling method constructed from the Arlequin framework. The coupling method provides for an approximation of the solution to a fully atomic model in which errors may arise due to the misplacement of the overlap region. The objective is thus to determine, a posteriori, the best position of this overlap region. More precisely, the problem is to find the optimal size of the atomic region that one needs to consider in the coupled formulation in order to accurately estimate predefined quantities of interest. In this new adaptive process, the position of the overlapping domain between the two models is conveniently parameterized and iteratively determined by searching for the optimal parameters. The performance of the method is demonstrated on one-dimensional and two-dimensional test problems. In particular, it is observed that the approach yields better convergence with respect to quantities of interest when compared to a classical adaptive approach based on a posteriori estimates of the modeling error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Path Integral Monte Carlo Simulations of Warm Dense Matter and Plasmas

New path integral Monte Carlo simulation (PIMC) techniques will be developed and applied to derive the equation of state (EOS) for the regime of warm dense matter and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 10,000 K, this method is not applicable at much higher temperature where electronic excitations become important because the number of partially occupied electronic orbitals reaches intractably large numbers and, more importantly, the use of zero-temperature exchange-correlation functionals introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it has commonly been assumed such methods can only be applied to elements without core electrons like hydrogen and helium, we recently showed how to extend PIMC to heavier elements by performing the first PIMC simulations of carbon and water plasmas [Driver, Militzer, Phys. Rev. Lett. 108 (2012) 115502]. Here we propose to continue this important development to extend the reach of PIMC simulations to yet heavier elements and also lower temperatures. The goal is to provide a robust first-principles simulation method that can accurately and efficiently study materials with excited electrons at solid-state densities in order to access parts of the phase diagram such the regime of warm dense matter and plasmas where so far only more approximate, semi-analytical methods could be applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Calculations of multiquark functions in effective models of strong interaction

In this paper we present our results of the investigation of multiquark equations in the Nambu-Jona-Lasinio model with chiral symmetry of SU(2) group in the mean-field expansion. To formulate the mean-field expansion we have used an iteration scheme of solution of the Schwinger-Dyson equations with the fermion bilocal source. We have considered the equations for Green functions of the Nambu-Jona-Lasinio model up to third step for this iteration scheme. To calculate the high-order corrections to the mean-field approximation, we propose the method of the Legendre transformation with respect to the bilocal source, which allows effectively to take into account the symmetry constraints related with the chiral Ward identity. We discuss also the problem of calculating the multiquark functions in the mean-field expansion for Nambu-Jona-Lasinio-type models with other types of the multifermion sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Application of a generalized matrix averaging method for the calculation of the effective properties of thin multiferroic layers

It is proposed to use a generalized matrix averaging (GMA) method for calculating the parameters of an effective medium with physical properties equivalent to those of a set of thin multiferroic layers. This approach obviates the need to solve a complex system of magnetoelectroelasticity equations. The required effective characteristics of a system of multiferroic layers are obtained using only operations with matrices, which significantly simplifies calculations and allows multilayer systems to be described. The proposed approach is applicable to thin-layer systems, in which the total thickness is much less than the system length, radius of curvature, and wavelengths of waves that can propagate in the system (long-wave approximation). Using the GMA method, it is also possible to obtain the effective characteristics of a periodic structure with each period comprising a number of thin multiferroic layers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Heavy-ion collisions and fission dynamics with the time-dependent Hartree–Fock theory and its extensions

Microscopic methods and tools to describe nuclear dynamics have considerably been improved in the past few years. They are based on the time-dependent Hartree–Fock (TDHF) theory and its extensions to include pairing correlations and quantum fluctuations. The TDHF theory is the lowest level of approximation of a range of methods to solve the quantum many-body problem, showing its universality to describe many-fermion dynamics at the mean-field level. The range of applications of TDHF to describe realistic systems allowing for detailed comparisons with experiment has considerably increased. For instance, TDHF is now commonly used to investigate fusion, multi-nucleon transfer and quasi-fission reactions. Thanks to the inclusion of pairing correlations, it has also recently led to breakthroughs in our description of the saddle to scission evolution, and, in particular, the non-adiabatic effects near scission. Beyond mean-field approaches such as the time-dependent random-phase approximation (TDRPA) and stochastic mean-field methods have reached the point where they can be used for realistic applications. We review recent progresses in both techniques and applications to heavy-ion collision and fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Plasmon-phonon mechanism of nonradiative transition 2S{sub 1/2} {yields} 2P{sub 1/2} in grazing scattering of a hydrogen atom at the surface of a solid

The probability of nonradiative transition 2S{sub 1/2} {yields} 2P{sub 1/2} in the case of grazing scattering of a hydrogen atom with an energy on the order of 1 keV at an Au and LiF surface is calculated in the semiclassical theory of nonadiabatic transitions. The transitions are induced by the interaction of the atom with the fields of surface plasmons and optical phonons, respectively. A formalism is developed which makes it possible to exceed the approximation based on the method of electrical images of charges. The value of the probability of scattering at LiF is high enough for experimental observation of the effect. The probability of this transition is determined for the motion of the atom parallel to the metal surface.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

e/a classification of Hume–Rothery Rhombic Triacontahedron-type approximants based on all-electron density functional theory calculations

There are three key electronic parameters in elucidating the physics behind the Hume–Rothery electron concentration rule: the square of the Fermi diameter (2kF)2, the square of the critical reciprocal lattice vector and the electron concentration parameter or the number of itinerant electrons per atom e/a. We have reliably determined these three parameters for 10 Rhombic Triacontahedron-type 2/1–2/1–2/1 (N = 680) and 1/1–1/1–1/1 (N = 160–162) approximants by making full use of the full-potential linearized augmented plane wave-Fourier band calculations based on all-electron density-functional theory. We revealed that the 2/1–2/1–2/1 approximants Al13Mg27Zn45 and Na27Au27Ga31 belong to two different sub-groups classified in terms of equal to 126 and 109 and could explain why they take different e/a values of 2.13 and 1.76, respectively. Among eight 1/1–1/1–1/1 approximants Al3Mg4Zn3, Al9Mg8Ag3, Al21Li13Cu6, Ga21Li13Cu6, Na26Au24Ga30, Na26Au37Ge18, Na26Au37Sn18 and Na26Cd40Pb6, the first two, the second two and the last four compounds were classified into three sub-groups with = 50, 46 and 42; and were claimed to obey the e/a = 2.30, 2.10–2.15 and 1.70–1.80 rules, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2018

BAM: bias assignment method to generate mock catalogues

Combining the statistics of dark matter tracers from a high-resolution cosmological N-body simulation and the dark matter density field calculated from down-sampled initial conditions using efficient structure formation solvers, we extract the halo-bias relation on a mesh of a 3$h^{-1}$ Mpc cell side resolution as a function of properties of the dark matter density field (e.g. local density, cosmic web type), automatically including stochastic, deterministic, local and non-local components. We use this information to sample the halo density field, accounting for ignored dependencies through an iterative process. By construction, our approach reaches ~1% accuracy in the majority of the k-range up to the Nyquist frequency without systematic deviations for power spectra (about $k ~ 1h$ Mpc -1 ) using either particle mesh or Lagrangian perturbation theory based solvers. When using phase-space mapping to compensate the low resolution of the approximate gravity solvers, our method reproduces the bispectra of the reference within 10% precision studying configurations tracing the quasi-non-linear regime. BAM has the potential to become a standard technique to produce mock halo and galaxy catalogues for future galaxy surveys and cosmological studies being highly accurate, efficient and parameter free.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Toward extracting $$\gamma $$ from $$B\rightarrow DK$$ without binning

Abstract $$B^\pm \rightarrow DK^\pm $$ B ± → D K ± transitions are known to provide theoretically clean information about the CKM angle $$\gamma $$ γ , with the most precise available methods exploiting the cascade decay of the neutral D into CP self-conjugate states. Such analyses currently require binning in the D decay Dalitz plot, while a recently proposed method replaces this binning with the truncation of a Fourier series expansion. In this paper, we present a proof of principle of a novel alternative to these two methods, in which no approximations at the level of the data representation are required. In particular, our new strategy makes no assumptions about the amplitude and strong phase variation over the Dalitz plot. This comes at the cost of a degree of ambiguity in the choice of test statistic quantifying the compatibility of the data with a given value of $$\gamma $$ γ , with improved choices of test statistic yielding higher sensitivity. While our current proof-of-principle implementation does not demonstrate optimal sensitivity to $$\gamma $$ γ , its conceptually novel approach opens the door to new strategies for $$\gamma $$ γ extraction. More studies are required to see if these can be competitive with the existing methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Quantitative Analysis of HAADF–STEM Images of MoVTeTaO M1 Phase Catalyst for Propane Ammoxidation to Acrylonitrile

Abstract High‐angle annular dark‐field scanning transmission electron microscopy (HAADF–STEM) image simulations were performed to determine accurate metal distributions in the MoVTeTaO M1 phase catalyst for propane ammoxidation. QSTEM simulation, based on a multislice method with the frozen phonon approximation, was chosen because of the excellent agreement between experimental and simulated HAADF–STEM images. This simulation method successfully provided accurate metal distributions in the MoVTeTaO M1 phase and indicated significantly lower V occupancies of the linking sites 1, 2, 3, 4, and 7, and higher Ta occupancy in site 9 than the previously reported metal occupancies. This novel methodology for accurate HAADF–STEM image analysis is a valuable tool for accurate analysis of local chemical compositions in a wide range of ordered catalytic materials.

Woo, Jungwon↗
DOE OSTI2016

Coupled wave-packets for non-adiabatic molecular dynamics: a generalization of Gaussian wave-packet dynamics to multiple potential energy surfaces

Accurate simulation of the non-adiabatic dynamics of molecules in excited electronic states is key to understanding molecular photo-physical processes. Here we present a novel method, based on a semiclassical approximation, that is as efficient as the commonly used mean field Ehrenfest or ad hoc surface hopping methods and properly accounts for interference and decoherence effects. This novel method is an extension of Heller's thawed Gaussian wave-packet dynamics that includes coupling between potential energy surfaces. By studying several standard test problems we demonstrate that the accuracy of the method can be systematically improved while maintaining high efficiency. The method is suitable for investigating the role of quantum coherence in the non-adiabatic dynamics of many-atom molecules.

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

Learning nonlinear level sets for dimensionality reduction in function approximation

We developed a Nonlinear Level-set Learning (NLL) method for dimensionality reduction in high-dimensional function approximation with small data. This work is motivated by a variety of design tasks in real-world engineering applications, where practitioners would replace their computationally intensive physical models (e.g., high-resolution fluid simulators) with fast-to-evaluate predictive machine learning models, so as to accelerate the engineering design processes. There are two major challenges in constructing such predictive models: (a) high-dimensional inputs (e.g., many independent design parameters) and (b) small training data, generated by running extremely time-consuming simulations. Thus, reducing the input dimension is critical to alleviate the over-fitting issue caused by data insufficiency. Existing methods, including sliced inverse regression and active subspace approaches, reduce the input dimension by learning a linear coordinate transformation; our main contribution is to extend the transformation approach to a nonlinear regime. Specifically, we exploit reversible networks (RevNets) to learn nonlinear level sets of a high-dimensional function and parameterize its level sets in low-dimensional spaces. A new loss function was designed to utilize samples of the target functions' gradient to encourage the transformed function to be sensitive to only a few transformed coordinates. The NLL approach is demonstrated by applying it to three 2D functions and two 20D functions for showing the improved approximation accuracy with the use of nonlinear transformation, as well as to an 8D composite material design problem for optimizing the buckling-resistance performance of composite shells of rocket inter-stages.

Zhang, Guannan↗
DOE OSTI2011

Finite-range effects in dilute Fermi gases at unitarity

We develop a theoretical method going beyond the contact-interaction approximation frequently used in mean-field theories of many-fermion systems, based on the low-energy T matrix of the pair potential to rigorously define the effective radius of the interaction. One of the main consequences of our approach is the possibility to investigate finite-density effects, which are outside the range of validity of approximations based on {delta}-like potentials. We apply our method to the calculation of density-dependent properties of an ultracold gas of {sup 6}Li atoms at unitarity, whose two-body interaction potential is calculated using ab initio quantum chemistry methods. We find that density effects will be significant in ultracold gases with densities 1 order of magnitude higher than those attained in current experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Three-dimensional thermo-mechanical simulations of heterogeneous solid propellants

Here in this work, we present a numerical framework that describes thermo-mechanical deformations in a burning heterogeneous solid propellant. These deformations are quasi-static at time scales associated with combustion, and the resulting thermo-mechanical formulation is discretized on a Cartesian grid using a hypoelastic law. We use a weak form of Chorin-type projection method to deal with large difference in shear modulus of the constituent materials. Extending our previous two-dimensional work, grid convergence studies for a three-dimensional propellant configuration are presented for the stress, velocity, and reference map components. Finally, simulations are carried out for a random propellant pack that is coupled to a gas phase, and we present results for the pack undergoing combustion, with and without deformations.

simulations↗
DOE OSTI2020

The effect of atmospheric transmissivity on model and observational estimates of the sea ice albedo feedback

The sea-ice-albedo feedback (SIAF) is the product of the ice sensitivity (IS) – how much the surface albedo in sea-ice regions changes as the planet warms– and the radiative sensitivity (RS) – how much the top of atmosphere radiation changes as the surface albedo changes. We demonstrate that the RS calculated from radiative kernels in climate models is reproduced from calculations using the “approximate partial radiative perturbation” method that uses the climatological radiative fluxes at the top of atmosphere and the assumption that the atmosphere is isotropic to shortwave radiation. This method facilitates the comparison of RS from satellite-based estimates of climatological radiative fluxes with RS estimates across a full suite of coupled climate models and, thus, allows model evaluation of a quantity important in characterizing the climate impact of sea ice concentration changes. The satellite based RS is within the model range of RS that differs by a factor of two across climate models in both the Arctic and Southern Ocean. Observed trends in Arctic sea ice are used to estimate IS which, in conjunction with the satellite-based RS yields an SIAF of 0.16 ± 0.04 W m-2 K-1. This Arctic SIAF estimate suggests a modest amplification of future global surface temperature change by approximately 14% relative to a climate system with no SIAF. We calculate the global albedo feedback in climate models using model specific RS and IS and find a model mean feedback parameter of 0.37 W m-2 K-1 which is 40% larger than the IPCC AR5 estimate based on using RS calculated from radiative kernel calculations in a single climate model.

Arctic, Antarctic, climate change, Climate & Earth↗
DOE OSTI2021

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗