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,171 records · Page 65

A note on the reliability of goal-oriented error estimates for Galerkin finite element methods with nonlinear functionals

Here, we consider estimating the discretization error in a nonlinear functional J (u) in the setting of an abstract variational problem: find u ϵ $\mathscr{V}$ such that B (u, φ) = L (φ) ∀φ ϵ $\mathscr{V}$, as approximated by a Galerkin finite element method. Here, $\mathscr{V}$ is a Hilbert space, B (. , .) is a bilinear form, and L (∙) is a linear functional. We consider well-known error estimates η of the form J (u) - J (u h ) ≈ η = L (z) - B (u h , z), where u h denotes a finite element approximation to u, and z denotes the solution to an auxiliary adjoint variational problem. We show that there exist nonlinear functionals for which error estimates of this form are not reliable, even in the presence of an exact adjoint solution z. An estimate η is said to be reliable if there exists a constant C ϵ $\mathbb{R}$ >0 independent of u h such that |J (u) - J (u h )| ≤ C|η|. We present several example pairs of bilinear forms and nonlinear functionals where reliability of η is not achieved.

A posteriori↗

Optimizing Storage Operation for a Probabilistic Locational Marginal Pricing Forecast: Preprint

This paper explores the applicability of using dynamic programing (DP) and approximate dynamic programming (ADP) based methods for optimal dispatch of utility scale energy storage systems (ESS). In this study, the effectiveness of these approaches have been tested using the IEEE 13 node test feeder with distributed photovoltaics (PVs) and a utility scale storage system. In this work, a co-simulation based approach has been used to setup the experiment to be able to implement detailed ESS and network models. The results obtained from DP/ADP runs have been compared with three other control strategies both myopic and intelligent. Simulations results show that DP/ADP algorithms are a good candidate for optimal EES dispatch in terms of both solution quality and execution time.

25 ENERGY STORAGE↗

Dirac cones in transition metal doped boron nitride

The transition metal (TM) doped zinc blende boron nitride (c-BN) is studied by using the first principle calculation. TM atoms fill in the interstitials in c-BN and form two-dimensional honeycomb lattice. The generalized gradient approximation and projector augmented wave method are used. The calculated density of states and band structures show that d electrons of TM atoms form impurity bands in the gap of c-BN. When the TM-BN system is in ferromagnetic or non-magnetic state, Dirac cones emerge at the K point in Brillouin zone. When TM is Ti and Co, the Dirac cones are spin polarized and very close to the Fermi level, which makes them promising candidates of Dirac half-metal [H. Ishizuka and Y. Motome, Phys. Rev. Lett. 109, 237207 (2012)]. While TM is Ni and Cu, the system is non-magnetic and Dirac cones located above the Fermi level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Linear Distribution System State Estimation Using Synchrophasor Data and Pseudo-Measurement

State estimation is often a challenging task in distribution systems due to deploying a limited number of measurement devices. Moreover, the integration of precise distributionlevel phasor measurement unites, a.k.a micro-PMUs, along with inaccurate pseudo-measurements in state estimation introduces another challenge. These issues might decrease the efficiency of traditional standard weighted least squares (WLS) for distribution system state estimation. This paper proposes a novel linear distribution system state estimation with taken different types of measurements, including micro-PMU measurements and line current measurements, into account. To involve pseudomeasurements into the linear state estimation, a linearization method based on the Taylor’s approximation is adopted to reformulate pseudo-measurement functions in a linear form. The results obtained via numerical simulations show that the proposed linear state estimation method performs similarly to the standard nonlinear WLS estimator. In addition, the sensitivity analyses results show that our method has a better performance compared to WLS once a limited number of highly precise micro-PMUs are accompanied with inaccurate pseudo-measurements in state estimation.

State estimation Distribution system Micro-PMU Lin↗

Effects of laser radiation field on energies of hydrogen atom in plasmas

In this study, for the first time, the Schrödinger equation with more general exponential cosine screened Coulomb (MGECSC) potential is solved numerically in the presence of laser radiation field within the Ehlotzky approximation using the asymptotic iteration method. The MGECSC potential includes four different potential forms in consideration of different sets of the parameters in the potential. By applying laser field, the total interaction potential of hydrogen atom embedded in plasmas converts to double well-type potential. The plasma screening effects under the influence of laser field as well as confinement effects of laser field on hydrogen atom in Debye and quantum plasmas are investigated by solving the Schrödinger equation with the laser-dressed MGECSC potential. It is resulted that since applying a monochromatic laser field on hydrogen atom embedded in a Debye and quantum plasma causes to shift in the profile of the total interaction potential, the confinement effects of laser field on hydrogen atom in plasmas modeled by the MGECSC potential change localizations of energy states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Smoothness of Functions in Spaces of the Finite Element Method

We find necessary and sufficient conditions for the generalized smoothness of the coordinate functions obtained from the approximate relations. We show that for the coordinate functions the smoothness on their supports is equivalent to that on the boundaries of supports. We obtain conditions for the continuity of Courant type finite element approximations and conditions for the uniqueness of a linear space of such approximations.

97 MATHEMATICS AND COMPUTING↗

Is Biological Subtype Prognostic of Locoregional Recurrence Risk in Women With pT1-2N0 Breast Cancer Treated With Mastectomy?

Purpose: To examine locoregional and distant recurrence (LRR and DR) in women with pT1-2N0 breast cancer according to approximated subtype and clinicopathologic characteristics. Methods and Materials: Two independent datasets were pooled and analyzed. The study participants were 1994 patients with pT1-2N0M0 breast cancer, treated with mastectomy without radiation therapy. The patients were classified into 1 of 5 subtypes: luminal A (ER+ or PR+/HER 2−/grade 1-2, n=1202); luminal B (ER+ or PR+/HER 2−/grade 3, n=294); luminal HER 2 (ER+ or PR+/HER 2+, n=221); HER 2 (ER−/PR−/HER 2+, n=105) and triple-negative breast cancer (TNBC) (ER−/PR−/HER 2−, n=172). Results: The median follow-up time was 4.3 years. The 5-year Kaplan-Meier (KM) LRR were 1.8% in luminal A, 3.1% in luminal B, 1.7% in luminal HER 2, 1.9% in HER 2, and 1.9% in TNBC cohorts (P=.81). The 5-year KM DR was highest among women with TNBC: 1.8% in luminal A, 5.0% in luminal B, 2.4% in luminal HER 2, 1.1% in HER 2, and 9.6% in TNBC cohorts (P<.001). Among 172 women with TNBC, the 5-year KM LRR were 1.3% with clear margins versus 12.5% with close or positive margins (P=.04). On multivariable analysis, factors that conferred higher LRR risk were tumors >2 cm, lobular histology, and close/positive surgical margins. Conclusions: The 5-year risk of LRR in our pT1-2N0 cohort treated with mastectomy was generally low, with no significant differences observed between approximated subtypes. Among the subtypes, TNBC conferred the highest risk of DR and an elevated risk of LRR in the presence of positive or close margins. Our data suggest that although subtype alone cannot be used as the sole criterion to offer postmastectomy radiation therapy, it may reasonably be considered in conjunction with other clinicopathologic factors including tumor size, histology, and margin status. Larger cohorts and longer follow-up times are needed to define which women with node-negative disease have high postmastectomy LRR risks in contemporary practice.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

COMPARISON OF IMPLICIT SCHEMES TO SOLVE EQUATIONS OF RADIATION HYDRODYNAMICS WITH A FLUX-LIMITED DIFFUSION APPROXIMATION: NEWTON–RAPHSON, OPERATOR SPLITTING, AND LINEARIZATION

Radiation is an important process of energy transport, a force, and a basis for synthetic observations, so radiation hydrodynamics (RHD) calculations have occupied an important place in astrophysics. However, although the progress in computational technology is remarkable, their high numerical cost is still a persistent problem. In this work, we compare the following schemes used to solve the nonlinear simultaneous equations of an RHD algorithm with the flux-limited diffusion approximation: the Newton–Raphson (NR) method, operator splitting, and linearization (LIN), from the perspective of the computational cost involved. For operator splitting, in addition to the traditional simple operator splitting (SOS) scheme, we examined the scheme developed by Douglas and Rachford (DROS). We solve three test problems (the thermal relaxation mode, the relaxation and the propagation of linear waves, and radiating shock) using these schemes and then compare their dependence on the time step size. As a result, we find the conditions of the time step size necessary for adopting each scheme. The LIN scheme is superior to other schemes if the ratio of radiation pressure to gas pressure is sufficiently low. On the other hand, DROS can be the most efficient scheme if the ratio is high. Although the NR scheme can be adopted independently of the regime, especially in a problem that involves optically thin regions, the convergence tends to be worse. In all cases, SOS is not practical.

79 ASTRONOMY AND ASTROPHYSICS↗

Method for Describing the Angular Distribution of Optical Radiation Scattered by a Monolayer of Ordered Spherical Particles (Normal Illumination)

We have developed a method for describing the angular distribution of intensity of radiation scattered by a monolayer of homogeneous spatially ordered monodisperse spherical particles normally illuminated by a plane circularly polarized electromagnetic wave. The method is based on the quasicrystalline approximation (QCA) of the theory of multiple scattering of waves (TMSW) using the multipole expansion of fields and the tensor Green function in vectorial spherical wavefunctions. The method is applied for analyzing the characteristics of radiation scattered by a partially ordered monolayer and a monolayer with a nonideal lattice. The results of calculations are compared with the available experimental data on the position of the first-order diffraction peak on the angular and spectral dependences of the intensity of radiation scattered by a closely packed monolayer with a nonideal triangular lattice of SiO{sub 2} particles. Good conformity of the results has been established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗

How to remove the spurious resonances from ring polymer molecular dynamics

Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here, we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the established features of RPMD: thermostatted RPMD is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here, we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Fourth-Order Symmetry Energy of Finite Nuclei

The fourth-order symmetry energy E{sub sym,4}(A) of heavy nuclei is investigated based on the Skyrme energy density functional in combination with a local density approximation. Unlike some previous works, in our method, the interferences from the other energy terms are removed since it is completely isolated from the rest of energy terms. The calculated E{sub sym,4}(A) is much less than that extracted from nuclear masses. The underlying reason for the big difference is discussed. The Brueckner theory also gives a small fourth-order symmetry energy coefficient of nuclear matter, which is also different from recent conclusions with another methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mapping Thermal Conductivity at the Atomic Scale: A Step toward the Thermal Design of Materials

We describe a spatial decomposition of the thermal conductivity, termed site-projected thermal conductivity, which quantifies the thermal conduction activity at each atomic site—a critical parameter for the thermal design of materials. The method is based on the Green–Kubo formula and the harmonic approximation and requires the force-constant and dynamical matrices, as well as a relaxed structural model. Throughout the study, it uses high quality models previously tested and compared to many experiments. It discusses the method and underlying approximations for amorphous silicon, carries the detailed analysis for amorphous silicon, and then examines an amorphous-crystal silicon interface and representative carbon materials. This study identifies the sites and local structures that reduce heat transport, and quantify these (estimate the spatial range) over which these “thermal defects” are effective. It identifies filamentary structures in the amorphous silicon network which impact heat transport and electronic structure (the Urbach edge) and electronic transport.

36 MATERIALS SCIENCE↗

Half-metallicity and stability of the rock salt BaC and SrC (111) surfaces: A density functional study

The electronic structure and magnetic properties of relaxed (111) surfaces of the alkaline-earth monocarbides BaC and SrC in the stable rock salt structure, are calculated on the basis of first principle density functional theory within the framework of self-consistent field plane wave pseudo-potential method, using the generalized gradient approximation for the exchange-correlation functional. The results of this study reveal that the C-terminated (111) surfaces retain the bulk half-metallic property in both BaC and SrC. The half-metallicity of the C-terminated BaC surface is found to be more robust compared to the bulk BaC due to the larger half-metallic energy gap. In contrast, the half-metallic energy gap of the C-terminated SrC surface is found to be smaller than that of the bulk. The Ba-terminated surface of BaC and the Sr-terminated surface of SrC, however, lose their bulk half-metallicity due to the formation of surface states in the majority spin band gap. The calculations also show that the atomic magnetic moments at the half-metallic C-terminated surfaces in both BaC and SrC increase considerably with respect to the corresponding bulk values, which is explained in terms of an increase in the number of unpaired 2p electrons of the carbon atom at the surface. We also discuss the stability of the surfaces via the calculated bulk formation energies. The bulk formation energies for both BaC and SrC in the rock salt structure are found to be positive, which indicate that the surfaces are not stable at normal pressure and temperature conditions, and non-equilibrium growth techniques may be required for the realization of BaC and SrC thin films.

74 ATOMIC AND MOLECULAR PHYSICS↗

Half-Metallic and Half-Semiconductor Gaps in Cr-Based Chalcogenides: DFT + U Calculations

Using full-potential linear-augmented plane waves plus local orbitals (FP-LAPW) method, within the generalized gradient approximation (GGA) and GGA + U (Hubbard Coulomb onsite correction) as well as the modified Becke-Johnson approach (mBJ-GGA and mBJ-GGA + U), we have carried out a systematic investigation of the structural, electronic, and magnetic properties of CrS, CrSe, and CrTe compounds in two competing phases: the zinc-blende (ZB) and wurtzite (WZ) structures. We obtained half-metallic (HM) gaps of 1.35 (2.32), 1.65 (2.49), and 1.52 (1.80) for the CrS, CrSe, and CrTe compounds, respectively, in the ZB phase within the GGA + U (mBJ-GGA + U) approach. The HM gaps are 1.78, 2.11, and 1.53 for CrS, CrSe, and CrTe compounds, respectively, in the WZ phase using the GGA + U. Notably, the CrSe and CrTe using the mBJ-GGA + U approach are a half-semiconductor (HSC) in the WZ phase. The HSC gaps are 0.94 and 0.25 eV for CrSe and CrTe compounds, respectively, while CrS has a HM gap of 2.30 eV. We observed a HSC-to-HM transition for the CrSe and CrTe crystal structures on application of strain. The CrS, CrSe, and CrTe compounds have a total magnetic moment of 4.0 µ{sub B} per formula unit with the magnetization predominantly from the Cr ion. The main features observed from the density of states and strain evolution should motivate further experimental exploration on the possible application of the CrS, CrSe, and CrTe compounds for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vibrational frequency scaling factors for correlation consistent basis sets and the methods CC2 and MP2 and their spin-scaled SCS and SOS variants

We present scaling factors for vibrational frequencies calculated within the harmonic approximation and the correlated wave-function methods coupled cluster singles and doubles model (CC2) and Møller-Plesset perturbation theory (MP2) with and without a spin-component scaling (SCS or spin-opposite scaling (SOS)). Frequency scaling factors and the remaining deviations from the reference data are evaluated for several non-augmented basis sets of the cc-pVXZ family of generally contracted correlation-consistent basis sets as well as for the segmented contracted TZVPP basis. We find that the SCS and SOS variants of CC2 and MP2 lead to a slightly better accuracy for the scaled vibrational frequencies. The determined frequency scaling factors can also be used for vibrational frequencies calculated for excited states through response theory with CC2 and the algebraic diagrammatic construction through second order and their spin-component scaled variants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗