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

Hybrid immersed interface-immersed boundary methods for AC dielectrophoresis

Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface–immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of the hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Cohesive energy and structural parameters of binary oxides of groups IIA and IIIB from diffusion quantum Monte Carlo

We have applied the diffusion quantum Monte Carlo (DMC) method to calculate the cohesive energy and the structural parameters of the binary oxides CaO, SrO, BaO, Sc{sub 2}O{sub 3}, Y{sub 2}O{sub 3}, and La{sub 2}O{sub 3}. The aim of our calculations is to systematically quantify the accuracy of the DMC method to study this type of metal oxides. The DMC results were compared with local, semi-local, and hybrid Density Functional Theory (DFT) approximations as well as with experimental measurements. The DMC method yields cohesive energies for these oxides with a mean absolute deviation from experimental measurements of 0.18(2) eV, while with local, semi-local, and hybrid DFT approximations, the deviation is 3.06, 0.94, and 1.23 eV, respectively. For lattice constants, the mean absolute deviations in DMC, local, semi-local, and hybrid DFT approximations are 0.017(1), 0.07, 0.05, and 0.04 Å, respectively. DMC is a highly accurate method, outperforming the DFT approximations in describing the cohesive energies and structural parameters of these binary oxides.

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

Realistic magnetic thermodynamics by local quantization of a semiclassical Heisenberg model

Abstract Classical Monte Carlo simulation of the Heisenberg model poorly describes many thermodynamic phenomena due to its neglect of the quantum nature of spins. Alternatively, we discuss how to semiclassically approach the quantum problem and demonstrate a simple method for introducing a locally approximate form of spin quantization. While the procedure underestimates magnetic short-range order, our results suggest a simple correction for recovering realistic spin–spin correlations above the critical temperature. Moreover, ensemble fluctuations are found to provide reasonably accurate thermodynamics, largely reproducing quantum mechanically calculated heat capacities and experimental magnetometry for ferromagnetic Fe and antiferromagnetic RbMnF 3 . Extensions of the method are proposed to address remaining inaccuracies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Online Adaptive Algorithm for Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method

Here in this research, we propose an online basis enrichment strategy within the framework of a recently developed constraint energy minimizing generalized multiscale discontinuous Galerkin method. Combining the technique of oversampling, one makes use of the information of the current residuals to adaptively construct basis functions in the online stage to reduce the error of multiscale approximation. A complete analysis of the method is presented, which shows the proposed online enrichment leads to a fast convergence from multiscale approximation to the fine-scale solution. The error reduction can be made sufficiently large by suitably selecting oversampling regions and the number of oversampling layers. Further, the convergence rate of the enrichment algorithm depends on a factor of exponential decay regarding the number of oversampling layers and a user-defined parameter. Numerical results are provided to demonstrate the effectiveness and efficiency of the proposed online adaptive algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Analysis of thermo-elastic problems using the improved element-free Galerkin method

The improved element-free Galerkin method (IEFG) is presented to deal with thermo-elastic problems. This mesh-free method is a combination between the element-free Galerkin method and the improved moving least-square approximation. It has not the Kronecker delta property, and the penalty method is used to impose the essential boundary conditions. In this paper, linear and stationary thermo-elasticity is treated. To solve the thermo-elastic problem, this latter is decoupled into two separate parts: first, the heat transfer problem is analyzed to reach the temperature field, which is used as input in the mechanical problem to calculate the displacement field and then the stress fields. Numerical examples with different boundary conditions are illustrated. The performance and the accuracy of the IEFG method are approved when obtained results are compared to finite-element results and analytical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Structural, electronic and elastic properties of RERu{sub 2} (RE=Pr and Nd) Laves phase intermetallic compounds

We have performed the first-principles calculations to study the structural, electronic and elastic properties of RERu{sub 2} (RE = Pr and Nd) Laves phase intermetallic compounds using full-potential linearized augmented plane wave (FP-LAPW) method based on density functional theory (DFT) within the generalized gradient approximation (GGA) for exchange and correlation potential. The optimized lattices constant are in reasonable agreement with available experimental data. The electronic properties are analyzed in terms of band structures, total and partial density of states, which confirm their metallic character. The calculated elastic constants infer that these compounds are mechanically stable in C15 (MgCu{sub 2} type) structure and found to be ductile in nature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

GDB-9-Ex_EOM-CCSD: Dataset containing Equation of Motion Coupled Cluster (EOM-CCSD) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 80,593 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Equation of Motion Coupled Cluster (EOM-CCSD) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The similarity-transformed EOM-CCSD method that used domain-based local pair natural orbitals (DLPNO) approximation which constitutes the STEOM-DLPNO-CCSD method was used. This method is based on the STEOM approach and was found to make accurate predictions of transition energies for organic molecules. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2015

Distribution System Secondary Circuit Parameter Estimation for Model Calibration

To analyze and coordinate the operation of distribution systems with rapidly increasing amounts of PV, more accurate distribution system models are required, especially for the distribution system secondary (low-voltage) circuits down to the point of common coupling for distributed PV. There is a growing need for automated procedures to calibrate the distribution system secondary circuit models that are typically either not modeled at all or are modeled with a lower level of detail than the better modeled medium-voltage systems. This report presents an accurate, flexible, and computationally efficient method to use measurement data to estimate secondary circuit series impedance parameters in existing utility feeder models. The parameter estimation method assumes well-modeled primary circuit models, known secondary circuit topologies, and AMI active power, and reactive power measurements at all the loads in the secondary circuit. The method also requires AMI voltage measurement at most of the loads in the secondary circuit but can handle loads that do not have voltage measurements. No existing secondary circuit model information is needed, except for topology. The method is based on the well-known linearized voltage drop approximation and linear regression. The performance of the method is demonstrated on a three-phase test circuit with ten different secondary circuit topologies and on the Georgia Tech campus distribution system with AMI data. The developed method can be utilized to improve existing utility feeder models for more accurate analysis and operation with ubiquitous distributed PV interconnected on the low-voltage circuits.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Resolution of Tolerance Variations

The connector/adapter assembly (CAA) is one component of a higher-level assembly. A simplified representation of the CAA is illustrated in Figure 1. The location of each CAA is required to meet a specific positional requirement. Several CAAs have been measured and found to be out of tolerance. The connector and the adapter are welded to each other such that the cylindrical axes of the two components are coaxial. The positional tolerance condition requires that all points on the outer surface of the CAA must be within a defined maximum diameter, as illustrated in Figure 1. Two different methods have been developed to measure the CAA positions from coordinate measuring machine (CMM) data. One method (RevD) extrapolates the data from two circles approximately midway from the weld region to each end. The second method (RevF) measures points that are much closer to the actual ends.

42 ENGINEERING↗
DOE OSTI2019

A kinematic comparison of meshfree and mesh-based Lagrangian approximations using manufactured extreme deformation fields

Meshfree methods for solid mechanics have been in development since the early 1990's. Initial motivations included alleviation of the burden of mesh creation and the desire to overcome the limitations of traditional mesh-based discretizations for extreme deformation applications. Here, the accuracy and robustness of both mesh-free and meshbased Lagrangian discretizations are compared using manufactured extreme-deformation fields. For the meshfree discretizations, both moving least squares and maximum entropy are considered. Quantitative error and convergence results are presented for the best approximation in the H 1 norm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method

The Molecular Beam Epitaxial model is derived from the variation of a free energy, that consists of either a fourth order Ginzburg–Landau double well potential or a nonlinear logarithmic potential in terms of the gradient of a height function. One challenge in solving the MBE model numerically is how to develop proper temporal discretization for the nonlinear terms in order to preserve energy stability at the time-discrete level. In this paper, we resolve this issue by developing a first and second order time-stepping scheme based on the “Invariant Energy Quadratization” (IEQ) method. The novelty is that all nonlinear terms are treated semi-explicitly, and the resulted semi-discrete equations form a linear system at each time step. Moreover, the linear operator is symmetric positive definite and thus can be solved efficiently. We then prove that all proposed schemes are unconditionally energy stable. The semi-discrete schemes are further discretized in space using finite difference methods and implemented on GPUs for high-performance computing. Various 2D and 3D numerical examples are presented to demonstrate stability and accuracy of the proposed schemes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Analysis of the Conformally Flat Approximation for Binary Neutron Star Initial Conditions

The spatially conformally flat approximation (CFA) is a viable method to deduce initial conditions for the subsequent evolution of binary neutron stars employing the full Einstein equations. Here in this paper, we analyze the viability of the CFA for the general relativistic hydrodynamic initial conditions of binary neutron stars. We illustrate the stability of the conformally flat condition on the hydrodynamics by numerically evolving ~100 quasicircular orbits. We illustrate the use of this approximation for orbiting neutron stars in the quasicircular orbit approximation to demonstrate the equation of state dependence of these initial conditions and how they might affect the emergent gravitational wave frequency as the stars approach the innermost stable circular orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Conforming and nonconforming virtual element methods for elliptic problems

Here we present, in a unified framework, new conforming and nonconforming virtual element methods for general second-order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and nonsymmetric parts and conditions for stability and accuracy on their discrete counterparts are established. These conditions are shown to lead to optimal H 1 - and L 2 -error estimates, confirmed by numerical experiments on a set of polygonal meshes. The accuracy of the numerical approximation provided by the two methods is shown to be comparable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Spatial locality of electronic correlations in LiFeAs

In this work, we address the question of the degree of spatial nonlocality of the self-energy in the iron-based superconductors, a subject which is receiving considerable attention. Using LiFeAs as a prototypical example, we extract the self-energy from angular-resolved photoemission spectroscopy data. We use two distinct electronic structure references: density functional theory in the local density approximation and linearized quasiparticle self-consistent GW (LQSGW). We find that with the LQSGW reference, spatially local dynamical correlations provide a consistent description of the experimental data, and account for some surprising aspects of the data such as the substantial out-of-plane dispersion of the electron Fermi surface having dominant xz/yz character. Hence, correlations effects can be separated into static nonlocal contributions well described by LQSGW and dynamical local contributions. Hall effect and resistivity data are shown to be consistent with this description.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗
DOE OSTI2018

First-Principle Calculations of Structural, Elastic, and Electronic Properties of Intermetallic Rare Earth R{sub 2}Ni{sub 2}Pb (R = Ho, Lu, and Sm) Compounds

The structural, elastic, and electronic properties of rare earth intermetallic R{sub 2}Ni{sub 2}Pb (where R = Ho, Lu, and Sm) compounds were investigated with the density functional theory (DFT) calculations. The calculations are performed using the full potential-linearized augmented plane wave (FP-LAPW) method within the framework of local density approximation (LDA). The calculated values of the equilibrium lattice constants were in agreement with the available experimental values. The elastic constants (C{sub ij}) were also calculated to understand the mechanical properties and structural stability of the compounds. Furthermore, the density of states and the charge density distributions of the compounds were calculated to understand the nature of the bonding in the material. The calculated results are in accordance with the available data in the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Quantification of uncertainties in thermoelectric properties of materials from a first-principles prediction method: An approach based on Gaussian process regression

We present the electron-phonon averaged via Gaussian process regression (EPA-GPR) method, in which the electron-phonon coupling matrix is represented as a function of two energies and is in turn modeled as a Gaussian process. The EPA-GPR method can be used as an efficient method to estimate thermoelectric properties of materials for fast-screening applications, comparable to the original electron-phonon averaged (EPA) method and the electron-phonon averaged via moving-least-squares (EPA-MLS) method. The EPA-GPR method does not require specification of any open parameter, unlike the other EPA-related methods, since all the hyperparameters in the model can be unambiguously estimated within the type II maximum likelihood (ML-II) approximation. Thus, the EPA-GPR method is a parameter-free estimation method. Additionally, the concept of Gaussian processes in the EPA-GPR method allows us to quantify the uncertainty in estimated properties of thermoelectric materials. One can randomly realize the electron-phonon coupling coefficients from the identified Gaussian process, and those realized samples can be further analyzed in the solution process of the semiclassical Boltzmann transport equation for charge carriers. The results of the semiclassical Boltzmann transport equation provide the statistical properties of the thermoelectric properties of interest. The means, standard deviations, histograms, and confidence intervals of the Seebeck coefficient, the electrical conductivity, and the power factor can be constructed and analyzed. Furthermore the proposed EPA-GPR method is applied to a p-type half-Heusler compound, i.e., HfCoSb, as a case example, the results of which clearly present the advantages of the method.

30 DIRECT ENERGY CONVERSION↗