Search NASA⌕ Search

SEARCH · Search NASA

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics

Practical methods for simulating photo-induced dynamics of molecules often rely on a trajectory-description of the nuclear dynamics, incorporating approximate electron-nuclear feedback. Through approximating electron–nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. Here, a set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.

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

Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory

We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hybrid eigensolvers for nuclear configuration interaction calculations

We examine and compare several iterative methods for solving large-scale eigenvalue problems arising from nuclear structure calculations. In particular, we discuss the possibility of using block Lanczos method, a Chebyshev filtering based subspace iterations and the residual minimization method accelerated by direct inversion of iterative subspace (RMM-DIIS) and describe how these algorithms compare with the standard Lanczos algorithm and the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm. Although the RMM-DIIS method does not exhibit rapid convergence when the initial approximations to the desired eigenvectors are not sufficiently accurate, it can be effectively combined with either the block Lanczos or the LOBPCG method to yield a hybrid eigensolver that has several desirable properties. We will describe a few practical issues that need to be addressed to make the hybrid solver efficient and robust.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Diabatic Hamiltonian matrix elements made simple

With a view to applying the generator coordinate method to large configuration spaces, we propose a simple approximate formula to compute diabatic many-body matrix elements without having to evaluate two-body interaction matrix elements. Here, the method is illustrated with two analytically solvable Hamiltonians based on the harmonic oscillator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Comparative Convergence Analysis of Nonlinear AMLI-Cycle Multigrid

The purpose of our paper is to provide a comprehensive convergence analysis of the nonlinear algebraic multilevel iteration (AMLI)-cycle multigrid (MG) method for symmetric positive definite problems. We show that the nonlinear AMLI-cycle MG method is uniformly convergent, based on classical assumptions for approximation and smoothing properties. Furthermore, under only the assumption that the smoother is convergent, we show that the nonlinear AMLI-cycle method is always better (or not worse) than the respective V-cycle MG method. Finally, numerical experiments are presented to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Evaluation of Subgrid-scale Hydrometeor Transport Schemes using a High-resolution Cloud-resolving Model

Potential ways of parameterizing vertical turbulent fluxes of hydrometeors are examined using a high-resolution cloud-resolving model. The cloud-resolving model uses the Morrison microphysics scheme, which contains prognostic variables for rain, graupel, ice, and snow. A benchmark simulation with a horizontal grid spacing of 250 m of a deep convection case carried out to evaluate three different ways of parameterizing the turbulent vertical fluxes of hydrometeors: an eddy-diffusion approximation, a quadrant-based decomposition, and a scaling method that accounts for within-quadrant (subplume) correlations. Results show that the down-gradient nature of the eddy-diffusion approximation tends to transport mass away from concentrated regions, whereas the benchmark simulation indicates that the vertical transport tends to transport mass from below the level of maximum to aloft. Unlike the eddy-diffusion approach, the quadri-modal decomposition is able to capture the signs of the flux gradient but underestimates the magnitudes. The scaling approach is shown to perform the best by accounting for within-quadrant correlations, and improves the results for all hydrometeors except for snow. A sensitivity study is performed to examine how vertical transport may affect the microphysics of the hydrometeors. The vertical transport of each hydrometeor type is artificially suppressed in each test. Results from the sensitivity tests show that cloud-droplet-related processes are most sensitive to suppressed rain or graupel transport. In particular, suppressing rain or graupel transport has a strong impact on the production of snow and ice aloft. Lastly, a viable subgrid-scale hydrometeor transport scheme in an assumed probability density function parameterization is discussed.

cloud-resolving model↗
DOE OSTI2024

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗
DOE OSTI2014

Direct X-B mode conversion for high-β national spherical torus experiment in nonlinear regime

Electron Bernstein wave (EBW) can be effective for heating and driving currents in spherical tokamak plasmas. Power can be coupled to EBW via mode conversion of the extraordinary (X) mode wave. The most common and successful approach to study the conditions for optimized mode conversion to EBW was evaluated analytically and numerically using a cold plasma model and an approximate kinetic model. The major drawback in using radio frequency waves was the lack of continuous wave sources at very high frequencies (above the electron plasma frequency), which has been addressed. A future milestone is to approach high power regime, where the nonlinear effects become significant, exceeding the limits of validity for present linear theory. Therefore, one appropriate tool would be particle in cell (PIC) simulation. The PIC method retains most of the nonlinear physics without approximations. In this work, we study the direct X-B mode conversion process stages using PIC method for incident wave frequency f{sub 0} = 15 GHz, and maximum amplitude E{sub 0} = 10{sup 5 }V/m in the national spherical torus experiment (NSTX). The modelling shows a considerable reduction in X-B mode conversion efficiency, C{sub modelling} = 0.43, due to the presence of nonlinearities. Comparison of system properties to the linear state reveals predominant nonlinear effects; EBW wavelength and group velocity in comparison with linear regime exhibit an increment around ∼36% and 17%, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Electric quadrupole transition probabilities for atomic lithium

Electric quadrupole transition probabilities for atomic lithium have been calculated using the weakest bound electron potential model theory (WBEPMT). We have employed numerical non-relativistic Hartree–Fock wavefunctions for expectation values of radii and the necessary energy values have been taken from the compilation at NIST. The results obtained with the present method agree very well with the Coulomb approximation results given by Caves (1975). Moreover, electric quadrupole transition probability values not existing in the literature for some highly excited levels have been obtained using the WBEPMT.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Surrogate Model Guided Optimization of Expensive Black-Box Multi-Objective Problems: A Posteriori Methods

Many engineering applications require the simultaneous optimization of multiple conflicting objective functions. Often, these objective functions are evaluated using highly accurate computer simulations that are computationally too expensive to be evaluated hundreds or thousands of times during optimization. Thus, the goal is to find good approximations of the Pareto front using as few of these expensive simulations as possible. Here, we describe an optimization approach based on surrogate models and diverse sampling strategies to accelerate the search for the Pareto solutions. We use a separate surrogate model for approximating each objective function and then we use the surrogate models to inform where additional expensive simulations should be run. The surrogate models are updated in an active learning framework whenever new information from the expensive simulations becomes available. The sampling strategies aim at balancing local improvements of the approximate Pareto front and global exploration to identify the extrema and fill in large gaps of the approximate Pareto front. We demonstrate on a large set of benchmark problems the effectiveness of the method for finding good approximations of the Pareto front.

MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Describing long-range charge-separation processes with subsystem density-functional theory

Long-range charge-transfer processes in extended systems are difficult to describe with quantum chemical methods. In particular, cost-effective (non-hybrid) approximations within time-dependent density functional theory (DFT) are not applicable unless special precautions are taken. Here, we show that the efficient subsystem DFT can be employed as a constrained DFT variant to describe the energetics of long-range charge-separation processes. A formal analysis of the energy components in subsystem DFT for such excitation energies is presented, which demonstrates that both the distance dependence and the long-range limit are correctly described. In addition, electronic couplings for these processes as needed for rate constants in Marcus theory can be obtained from this method. It is shown that the electronic structure of charge-separated states constructed by a positively charged subsystem interacting with a negatively charged one is difficult to converge — charge leaking from the negative subsystem to the positive one can occur. This problem is related to the delocalization error in DFT and can be overcome with asymptotically correct exchange–correlation (XC) potentials or XC potentials including a sufficiently large amount of exact exchange. We also outline an approximate way to obtain charge-transfer couplings between locally excited and charge-separated states.

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

Uncertainty Quantification of Spatial Correction Factor for Sjoestrand Method due to Cross-section Data

Accelerator-driven System (ADS), which couples a subcritical core with a spallation neutron source, is proposed as a nuclear transmutation technology to reduce the radioactive half-life of high-level waste. To maintain a safe operation in ADS, the core should be kept in subcritical state. Thus, the subcritical state of the core must be quantified through measurement of subcriticality, which is a quantitative margin from critical state. The area ratio method, or the Sjoestrand method, is one of the subcriticality measurement techniques applicable to the ADS. In this method, the absolute value of subcriticality can be obtained by analyzing the temporal variation of neutron density in the core due to the periodic injection of neutron pulses. In practice, the value of subcriticality is measured by the temporal variation of neutron count rate detected by a detector in the core. Time integrals of count rates due to prompt and delayed neutrons A{sub p} and A{sub d}, are measured. Then, the ratio of A{sub p} to A{sub d}, which is called as the area ratio, yields negative reactivity -ρ, i.e. subcriticality, in dollar units: A{sub p}/A{sub d} ≅ -ρ/β{sub eff'} (1) where β{sub eff} is the effective delayed neutron fraction. Equation (1) represents the basic idea of the Sjoestrand method. Note, however, that Eq. (1) is approximately derived using the point kinetics model. In an actual subcritical system, the measured area ratio is not necessarily equal to the exact value of subcriticality in dollar units, due distortion by spatial higher-order modes. Hence, in the actual measurement, this spatial effect is corrected by a spatial correction factor such as the Bell and Glasstone correction factor (the Bell factor), or by extraction of the fundamental mode (so-called the extrapolated area ratio method). In addition, there is an effect due to the neutron detector materials (e.g. {sup 3}He(n,p) and {sup 235}U(n,f)). These effects must be considered in the correction of the measured area ratio. In this study, we focused on the uncertainty of Bell factor. The Bell factor f is defined by the ratio of the area ratio to the subcriticality as follows: f ≡ -[ρ/β{sub eff}]/[A{sub p}/A{sub d}] (2) The Bell factor f can be evaluated by a core analysis. By multiplying the measured area ratio by the Bell factor?, the spatial higher-order mode effect in the area ratio can be corrected. However, the uncertainty of the Bell factor due to uncertainties of cross-section data (covariance data) has not been quantitatively evaluated so far. In this study, we propose a method to quantify the uncertainty of the Bell factor in an actual core analysis. As a feasibility study, the uncertainty quantification is carried out in a simple one-dimensional subcritical core system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Direct Delta-MBPT(2) method for ionization potentials, electron affinities, and excitation energies using fractional occupation numbers

A direct method (D-Delta-MBPT(2)) to calculate second-order ionization potentials (IPs), electron affinities (EAs), and excitation energies is developed. The Delta-MBPT(2) method is defined as the correlated extension of the Delta-HF method. Energy differences are obtained by integrating the energy derivative with respect to occupation numbers over the appropriate parameter range. This is made possible by writing the second-order energy as a function of the occupation numbers. Relaxation effects are fully included at the SCF level. This is in contrast to linear response theory, which makes the D-Delta-MBPT(2) applicable not only to single excited but also higher excited states. We show the relationship of the D-Delta-MBPT(2) method for IPs and EAs to a second-order approximation of the effective Fock-space coupled-cluster Hamiltonian and a second-order electron propagator method. We also discuss the connection between the D-Delta-MBPT(2) method for excitation energies and the CIS-MP2 method. Finally, as a proof of principle, we apply our method to calculate ionization potentials and excitation energies of some small molecules. For IPs, the Delta-MBPT(2) results compare well to the second-order solution of the Dyson equation. For excitation energies, the deviation from EOM-CCSD increases when correlation becomes more important. When using the numerical integration technique, we encounter difficulties that prevented us from reaching the Delta-MBPT(2) values. Most importantly, relaxation beyond the Hartree Fock level is significant and needs to be included in future research.

Beste, Ariana↗
DOE OSTI2021

van der Waals corrected density functionals for cylindrical surfaces: Ammonia and nitrogen dioxide adsorbed on a single-walled carbon nanotube

In this work, we extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, NH3 and NO2 molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70–115 meV for the system of CNT + NH 3 and 300–500 meV for the system of CNT + NO 2 . The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW–dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to –4 to –4.5 at distances D about 20–50 Å. At even larger distances, the vdW power-law exponent approaches –5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

The influence of defects on magnetic properties of fcc-Pu

The influence of vacancies and interstitial atoms on magnetism in Pu is considered in the framework of the density functional theory. The crystal structure relaxation arising due to various types of defects is calculated using the molecular dynamics method with a modified embedded atom model. The local density approximation with explicit inclusion of Coulomb and spin-orbit interactions is applied in matrix invariant form to describe correlation effects in Pu with these types of defects. The calculations show that both vacancies and interstitials give rise to local moments in the f-shell of Pu in good agreement with experimental data for aged Pu. Magnetism appears due to the destruction of a delicate balance between spin-orbit and exchange interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

BHResist [SWR-25-57]

BHResist is a pure python library for computing thermal resistance within single-u, double-u, and coaxial grouted borehole heat exchangers. For single and double u-tube configurations, the methods use the 1st-order closed-form multipole approximations, which typically produces results with less than 1% error when compared to the 10th-order multipole method. Coaxial borehole methods apply a simple 1D resistance network method. This is intended to be a lightweight library that can be easily imported into any other Python tool, with no bulky dependencies.

Mitchell, Matt [National Renewable Energy Laborato↗
DOE OSTI2015

Development of a New 47-Group Library for the CASL Neutronics Simulators

The CASL core simulator MPACT is under development for the neutronics and thermal-hydraulics coupled simulation for the pressurized light water reactors. The key characteristics of the MPACT code include a subgroup method for resonance self-shielding, and a whole core solver with a 1D/2D synthesis method. The ORNL AMPX/SCALE code packages have been significantly improved to support various intermediate resonance self-shielding approximations such as the subgroup and embedded self-shielding methods. New 47-group AMPX and MPACT libraries based on ENDF/B-VII.0 have been generated for the CASL core simulator MPACT of which group structure comes from the HELIOS library. The new 47-group MPACT library includes all nuclear data required for static and transient core simulations. This study discusses a detailed procedure to generate the 47-group AMPX and MPACT libraries and benchmark results for the VERA progression problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

THE BEHAVIOR OF TRANSVERSE WAVES IN NONUNIFORM SOLAR FLUX TUBES. I. COMPARISON OF IDEAL AND RESISTIVE RESULTS

Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius of the tube, i.e., the so-called thin boundary approximation. Here, we develop a general analytic method to compute the dispersion relation and the eigenfunctions of ideal MHD waves in pressureless flux tubes with transversely nonuniform layers of arbitrary thickness. Results for kink waves are produced and compared with fully numerical resistive MHD eigenvalue computations in the limit of small resistivity. We find that the frequency and resonant damping rate are the same in both ideal and resistive cases. The actual results for thick nonuniform layers deviate from the behavior predicted in the thin boundary approximation and strongly depend on the shape of the nonuniform layer. The eigenfunctions in ideal MHD are very different from those in resistive MHD. The ideal eigenfunctions display a global character regardless of the thickness of the nonuniform layer, while the resistive eigenfunctions are localized around the resonance and are indistinguishable from those of ordinary resistive Alfvén modes. Consequently, the spatial distribution of wave energy in the ideal and resistive cases is dramatically different. This poses a fundamental theoretical problem with clear observational consequences.

79 ASTRONOMY AND ASTROPHYSICS↗