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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Block Preconditioners for Stable Mixed Nodal and Edge finite element Representations of Incompressible Resistive MHD

The scalable iterative solution of strongly coupled three-dimensional incompressible resistive magnetohydrodynamics (MHD) equations is quite challenging because disparate time scales arise from the electromagnetics, the hydrodynamics, as well as the coupling between these systems. This study considers a mixed finite element discretization of a dual saddle point formulation of the incompressible resistive MHD equations using a stable nodal (Q2/Q1) discretization for the hydrodynamics and a stable edge-node discretization of a reduced form of the Maxwell equations. This paper introduces new approximate block factorization preconditioners for this system which reduce the system to approximate Schur complement systems that can be solved using algebraic multilevel methods. These preconditioners include a new augmentation-based approximation for the magnetic induction saddle point system as well as efficient approximations of the Schur complements that arise from the complex coupling between the Navier--Stokes equations and the Maxwell equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Model-based predictive control of multi-zone commercial building with a lumped building modelling approach

Here this study investigates the applicability of a lumped building modeling approach to model-based predictive control (MPC) to alleviate the complex modeling process of the grey-box multi-zone building model. Based on experimental data, two building models were estimated in this study. The detailed model as a reference case and a lumped model were estimated with decentralized and conventional approaches, respectively. Then, simulations were performed with two boundary conditions, including the comfort bound and electricity cost structure. The performances of the MPC with the detailed and lumped models were analyzed compared to the feedback control. More savings was achieved with a larger comfort bound and more aggressive electricity cost structure. The savings potential of the proposed lumped model approach was not as high as that of the detailed model. However, the proposed method yields good control performance, whose savings was approximately 8.6% over that of feedback control. These results suggest that the proposed method can be used to facilitate MPC implementation in multi-zone building applications.

42 ENGINEERING↗
DOE OSTI2021

Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model

An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.

Westheimer, Bryce M. (ORCID:0000000267645714)↗
DOE OSTI2018

PV Inverter Fault Response Including Momentary Cessation, Frequency-Watt, and Virtual Inertia

This paper presents two photovoltaic (PV) inverter control methods and an analysis of the two under a significant three-phase transmission line fault contingencies in the Hawaiian island of Oahu power system. Simulations are presented for the two control methods on the island power system with a high penetration of PV generation, approximately 40% of the total. The two control methods discussed are similar: one is a proportional frequency-watt controller, and the second is a proportional-derivative controller, a frequency-watt controller with virtual inertia. Both methods can be beneficial for fast frequency support. The inverter model also emulates inverter fault response including 'momentary cessation' and recovery during low voltage events. The study presented in this paper utilizes a validated power system model of the Hawaiian island of Oahu, modified to include PV resources with the two custom developed control models. Dynamic simulations with PSS/E are presented during a significant transmission line three-phase fault contingency. Simulations are presented with and without PV reserve margin. In addition, a parameter sensitivity analysis is presented for the control methods. Results indicate both methods can significantly improve system response during fault events. Findings indicate that transmission faults can produce severe frequency events, and that fast recovery from momentary cessation is crucial to mitigate severity.

14 SOLAR ENERGY↗
DOE OSTI2016

A new insight into neutrino energy loss by electron capture of iron group nuclei in magnetar surfaces

Based on the relativistic mean-field effective interactions theory, and the Lai dong model, we discuss the influences of superstrong magnetic fields (SMFs) on electron Fermi energy, nuclear blinding energy, and single-particle level structure in magnetar surfaces. Using the Shell-Model Monte Carlo method and the Random Phase Approximation theory, we analyze the neutrino energy loss rates (NELRs) by electron capture for iron group nuclei in SMFs. First, when B {sub 12} < 100, we find that the SMFs have a slight influence on the NELRs for most nuclides at relativistic low temperatures (e.g., T {sub 9} = 0.233); nevertheless, the NELRs increase by more than four orders of magnitude at relativistic high temperatures (e.g., T {sub 9} = 15.53). When B {sub 12} > 100, the NELRs decrease by more than three orders of magnitude (e.g., at T {sub 9} = 15.53 for {sup 52–61}Fe, {sup 55–60}Co, and {sup 56–63}Ni). Second, for a certain value of magnetic field and temperature, the NELRs increase by more than four orders of magnitude when ρ{sub 7}⩽10{sup 3}, but as the density increases (i.e., when ρ{sub 7}>10{sup 3}), there is almost no influence on the density of NELRs. For the density around ρ{sub 7}=10{sup 2}, there is an abrupt increase in NELRs when B {sub 12} ≥ 10{sup 3.5}. Such jumps are an indication that the underlying shell structure has changed due to single-particle behavior by SMFs. Finally, we compare our NELRs with those of Fuller et al. (FFN) and Nabi and Klapdor-Kleingrothaus (NKK). For the case without SMFs, one finds that our rates for certain nuclei are close to about five orders of magnitude lower than FFN and NKK at relativistic low temperatures (e.g., T {sub 9} = 1). However, at a relativistic high temperature (e.g., T {sub 9} = 3), our results are in good agreement with NKK, but about one order of magnitude lower than FFN. For the case with SMFs, our NELRs for some iron group nuclei can be about five orders of magnitude higher than those of FFN and NKK. (Note that B {sub 12}, T {sub 9}, and ρ {sub 7} are in units of 10{sup 12} G, 10{sup 9} K, and 10{sup 7} g cm{sup −3}, respectively.)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Hyperfine interactions for small systems including transition-metal elements using self-interaction corrected density-functional theory

The interactions between the electronic magnetic moment and the nuclear spin moment, i.e., magnetic hyperfine (HF) interactions, play an important role in understanding electronic properties of magnetic systems and in realizing platforms for quantum information science applications. We investigate the HF interactions for atomic systems and small molecules, including Ti or Mn, by using Fermi–Löwdin orbital (FLO) based self-interaction corrected (SIC) density-functional theory. We calculate the Fermi contact (FC) and spin-dipole terms for the systems within the local density approximation (LDA) in the FLO-SIC method and compare them with the corresponding values without SIC within the LDA and generalized-gradient approximation (GGA), as well as experimental data. For the moderately heavy atomic systems (atomic number Z ≤ 25), we find that the mean absolute error of the FLO-SIC FC term is about 27 MHz (percentage error is 6.4%), while that of the LDA and GGA results is almost double that. Therefore, in this case, the FLO-SIC results are in better agreement with the experimental data. For the non-transition-metal molecules, the FLO-SIC FC term has the mean absolute error of 68 MHz, which is comparable to both the LDA and GGA results without SIC. For the seven transition-metal-based molecules, the FLO-SIC mean absolute error is 59 MHz, whereas the corresponding LDA and GGA errors are 101 and 82 MHz, respectively. Therefore, for the transition-metal-based molecules, the FLO-SIC FC term agrees better with experiment than the LDA and GGA results. We observe that the FC term from the FLO-SIC calculation is not necessarily larger than that from the LDA or GGA for all the considered systems due to the core spin polarization, in contrast to the expectation that SIC would increase the spin density near atomic nuclei, leading to larger FC terms.

Chemistry↗
DOE OSTI2017

Optical potential from first principles

Here, we develop a method to construct a microscopic optical potential from chiral interactions for nucleon-nucleus scattering. The optical potential is constructed by combining the Green’s function approach with the coupled-cluster method. To deal with the poles of the Green’s function along the real energy axis we employ a Berggren basis in the complex energy plane combined with the Lanczos method. Using this approach, we perform a proof-of-principle calculation of the optical potential for the elastic neutron scattering on 16 O. For the computation of the ground-state of 16 O, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A ±1 nuclei we use particle-attached/removed equation-of-motion method truncated at two-particle-one-hole and one-particle-two-hole excitations, respectively. We verify the convergence of the optical potential and scattering phase shifts with respect to the model-space size and the number of discretized complex continuum states. We also investigate the absorptive component of the optical potential (which reflects the opening of inelastic channels) by computing its imaginary volume integral and find an almost negligible absorptive component at low-energies. To shed light on this result, we computed excited states of 16 O using equation-of-motion coupled-cluster method with singles-and- doubles excitations and we found no low-lying excited states below 10 MeV. Furthermore, most excited states have a dominant two-particle-two-hole component, making higher-order particle-hole excitations necessary to achieve a precise description of these core-excited states. We conclude that the reduced absorption at low-energies can be attributed to the lack of correlations coming from the low-order cluster truncation in the employed coupled-cluster method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Successful In Situ Remediation of Cr(VI)-Impacted Deeper Vadose Zone Soils for Groundwater Protection - 18603

Shallow soil at a former industrial wood treating facility was impacted with chromium, arsenic, and copper above cleanup goals and shallow soil has been excavated for protection of direct contact exposure pathways. Deeper vadose zone soil was predominantly impacted with chromium and was undergoing remediation for long-term groundwater protection. The 920-meter (m) long Cr(VI) groundwater plume has undergone full-scale remediation with groundwater recirculation and in situ bioremediation. Following approximately 4 years of operation, the Cr(VI) groundwater plume, with former source area groundwater concentrations over 10,000 micrograms per liter (μg/L), has been reduced to below the laboratory method reporting limit (0.50 μg/L). Approximately 23,000 cubic meters (m{sup 3}) of deeper (4.5 m - 11 m below ground surface, bgs) vadose zone soil required Cr(VI) remediation. The remediation strategy relied on reduction of Cr(VI) and geochemical fixation as trivalent chromium (Cr(III)) and iron hydroxide precipitates. Chromium reduction was facilitated by a carbon source (corn syrup and methanol) to establish and maintain biogeochemical reducing conditions and ferrous iron to support direct Cr(VI) chemical reduction. Combined, the injected soil moisture, carbon source and ferrous iron supported fast and sustained ferrous/ferric iron cycling and Cr(VI) reduction. Remedy performance monitoring included time-series water quality data from a network of vacuum lysimeters and confirmation soil borings. Reagent was applied either through liquid injection wells and direct-push injection points, or by deep soil mixing methodology in areas with difficult geology (i.e., fine-grained soils and discontinuous iron-silica hardpan layer). Following both phases of deep soil remediation, soils met long-term protection of groundwater remediation action objectives. In situ treatment of soils provides a cost-effective remedial solution compared to excavation, specifically when impacts are deep and long-term protection of groundwater is a primary remediation goal. The delivery methods employed in this case study can be applied to sites impacted with a broad range of contaminants amendable to geochemical fixation (e.g., chromium, uranium, strontium, nitrate, etc.). (authors)

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

Multigroup Neutron Transport Using a Collision-Based Hybrid Method

A collision-based hybrid algorithm for the discrete ordinates approximation of the neutron transport equation is extended to the isotropic multigroup setting. The algorithm uses discrete energy and angle grids at two different resolutions and approximates the fission and scattering sources on the coarser grids. The coupling of a collided transport equation, discretized on the coarse grid, with an uncollided transport equation, discretized on the fine grid, yields an algorithm that, in most cases, is more efficient than the traditional multigroup approach. In conclusion, the improvement over existing techniques is demonstrated for time-dependent problems with different materials, geometries, and energy groups.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Regularity of Solutions to Quasilinear Parabolic Systems with Time-Nonsmooth Principal Matrix and the Neumann Boundary Condition

We consider a quasilinear parabolic system of equations with nondiagonal principal matrix in a model parabolic cylinder with the Neumann condition on the plane part Γ of the lateral surface of the cylinder. We prove the partial regularity (the Hölder continuity) of the weak solution in a neighborhood of Γ by the method of A(t)-caloric approximations adapted to the problem with the Neumann boundary condition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Unraveling Site Selective Magnetic Properties of Cobalt Sites in Critical Elements Lean RE(TM)5 Magnet Materials

We report here our discovery of crystallographic and interstitial sites and onsite electron correlation propelled intrinsic and derived permanent magnetic properties of critical elements lean RE(TM) 5 (RE = La, Ce and TM = Fe, Co) magnet materials. A full potential linearized augmented plane wave (FP-LAPW) method within the local density approximation (LDA) is used to investigate and analyze the electronic structure and magnetism of these RE(TM) 5 type structures. To better correlate the experimental results, the effective Coulomb (U) and exchange (J) interactions (Hubbard parameters) are crucial at the transition metal sites. Results show that the main propeller of magnetic anisotropy in these compounds is the cobalt atoms at the 2c sites not the 3g sites. This site-specific property and site preference energetics are used to replace 3g sites with non-critical elements such as iron that exhibits a larger magnetic moment. Based on this strategic replacement, we predict two new compounds that have a larger hardness parameter with a relatively large energy product due to the atomic dilution caused by interstitial addition of nitrogen in the compounds: CeCo 2 Fe 3 N 2 and LaCo 2 Fe 3 N 2 .

Rare-earths↗
DOE OSTI2018

A comprehensive first-principles study of solute elements in dilute Ni alloys: Diffusion coefficients and their implications to tailor creep rate

Diffusion regulates a vast number of materials properties and phenomena such as creep, the focus of the present work. However, a deep understanding of the effect of how each alloying element in a Ni-base superalloy affects properties such as diffusion and creep is far from complete. Here, in this study, we report temperature-dependent dilute solute diffusion coefficients and their implications to tailor the creep rate for 26 transition metal solute elements, X's, in fcc Ni using first-principles calculations. Calculations are performed using the five-frequency model for dilute solute diffusion and the nudged elastic band method within the local density approximation. Thermodynamic properties at finite temperatures for all configurations are calculated using the quasi-harmonic Debye model. In general, the fastest diffusing solute elements in Ni are found at the left side of the periodic table and the slowest diffusing solute elements are found in group VIIB. In particular, the present work indicates that the diffusivity of the dilute solute elements is correlated to the compressibility of each solute element on the respective Ni 31 X supercell, and not as strongly to the ionic radius of the solute elements, as previously suggested. Finally, results from the diffusivity study are combined with the previous results of elastic constants and stacking fault energies, and hence, a relative creep rate ratio for these 26 solute elements is modeled. It is shown that in most cases, slower diffusing solute elements provide the most creep resistance. This is true even at higher temperatures, due mainly to the solute's strong bonding with Ni atoms in the host lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The SCALE/AMPX multigroup cross section processing for fast reactor analysis

The SCALE/AMPX multigroup (MG) cross section processing procedure has been updated to minimize reactivity differences for fast reactor designs and boiling water reactors (BWRs) with very high void fractions to provide excellent agreement with continuous-energy reference calculations. SCALE MG calculations are widely applied to thermal spectrum light-water reactor (LWR) systems as well as fast spectrum metallic systems of interest to National Nuclear Security Administration. With growing interest from industry and regulators in applying SCALE for the design of fast spectrum reactors, both sodium and molten salt, and in the licensing of power up rates for BWRs, it is desirable to review the SCALE/AMPX procedure for unresolved resonance self-shielded data and high energy neutron spectra. The data were improved by generating MG unresolved resonance data based on the analytic probability table method with the narrow resonance approximation as well as the use of a very fine group structure typically used in fast system analysis. This study focused on verifying the probability table and the SCALE/AMPX MG cross section processing procedure by performing reaction rate analysis and benchmark calculations for various fast reactor and high void BWR systems. Here, results indicate that the improved SCALE/AMPX MG cross section processing provides excellent results for the fast reactor analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

An efficient computational method for solving nonlinear stochastic Itô integral equations: Application for stochastic problems in physics

Because of the nonlinearity, closed-form solutions of many important stochastic functional equations are virtually impossible to obtain. Thus, numerical solutions are a viable alternative. In this paper, a new computational method based on the generalized hat basis functions together with their stochastic operational matrix of Itô-integration is proposed for solving nonlinear stochastic Itô integral equations in large intervals. In the proposed method, a new technique for computing nonlinear terms in such problems is presented. The main advantage of the proposed method is that it transforms problems under consideration into nonlinear systems of algebraic equations which can be simply solved. Error analysis of the proposed method is investigated and also the efficiency of this method is shown on some concrete examples. The obtained results reveal that the proposed method is very accurate and efficient. As two useful applications, the proposed method is applied to obtain approximate solutions of the stochastic population growth models and stochastic pendulum problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗