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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Stochastic Least-Squares Petrov--Galerkin Method for Parameterized Linear Systems

We consider the numerical solution of parameterized linear systems where the system matrix, the solution, and the right-hand side are parameterized by a set of uncertain input parameters. We explore spectral methods in which the solutions are approximated in a chosen finite-dimensional subspace. It has been shown that the stochastic Galerkin projection technique fails to minimize any measure of the solution error [A. Mugler and H.-J. Starkloff, ESAIM Math. Model. Numer. Anal., 47 (2013), pp. 1237--1263]. As a remedy for this, we propose a novel stochatic least-squares Petrov--Galerkin (LSPG) method. The proposed method is optimal in the sense that it produces the solution that minimizes a weighted $\ell^2$-norm of the residual over all solutions in a given finite-dimensional subspace. Moreover, the method can be adapted to minimize the solution error in different weighted $\ell^2$-norms by simply applying a weighting function within the least-squares formulation. In addition, a goal-oriented seminorm induced by an output quantity of interest can be minimized by defining a weighting function as a linear functional of the solution. We establish optimality and error bounds for the proposed method, and extensive numerical experiments show that the weighted LSPG method outperforms other spectral methods in minimizing corresponding target weighted norms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are impractical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The empirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which satisfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Backscattering and absorption coefficients for electrons: Solutions of invariant embedding transport equations using a method of convergence

The backscattering coefficient is a magnitude whose measurement is fundamental for the characterization of materials with techniques that make use of particle beams and particularly when performing microanalysis. In this work, we report the results of an analytic method to calculate the backscattering and absorption coefficients of electrons in similar conditions to those of electron probe microanalysis. Starting on a five level states ladder model in 3D, we deduced a set of integro-differential coupled equations of the coefficients with a method know as invariant embedding. By means of a procedure proposed by authors, called method of convergence, two types of approximate solutions for the set of equations, namely complete and simple solutions, can be obtained. Although the simple solutions were initially proposed as auxiliary forms to solve higher rank equations, they turned out to be also useful for the estimation of the aforementioned coefficients. In previous reports, we have presented results obtained with the complete solutions. In this paper, we present results obtained with the simple solutions of the coefficients, which exhibit a good degree of fit with the experimental data. Both the model and the calculation method presented here can be generalized to other techniques that make use of different sorts of particle beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Space-time least-squares Petrov-Galerkin projection in nonlinear model reduction.

Our work proposes a space-time least-squares Petrov-Galerkin (ST-LSPG) projection method for model reduction of nonlinear dynamical systems. In contrast to typical nonlinear model-reduction methods that first apply Petrov-Galerkin projection in the spatial dimension and subsequently apply time integration to numerically resolve the resulting low-dimensional dynamical system, the proposed method applies projection in space and time simultaneously. To accomplish this, the method first introduces a low-dimensional space-time trial subspace, which can be obtained by computing tensor decompositions of state-snapshot data. The method then computes discrete-optimal approximations in this space-time trial subspace by minimizing the residual arising after time discretization over all space and time in a weighted ℓ 2 -norm. This norm can be de ned to enable complexity reduction (i.e., hyper-reduction) in time, which leads to space-time collocation and space-time GNAT variants of the ST-LSPG method. Advantages of the approach relative to typical spatial-projection-based nonlinear model reduction methods such as Galerkin projection and least-squares Petrov-Galerkin projection include: (1) a reduction of both the spatial and temporal dimensions of the dynamical system, (2) the removal of spurious temporal modes (e.g., unstable growth) from the state space, and (3) error bounds that exhibit slower growth in time. Numerical examples performed on model problems in fluid dynamics demonstrate the ability of the method to generate orders-of-magnitude computational savings relative to spatial-projection-based reduced-order models without sacrificing accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An adaptive moments-based interface reconstruction using intersection of the cell with one half-plane, two half-planes and a circle

We present a new adaptive moment-of-fluid (A-MOF) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to construct a shape that approximates the respective material fragment. The new method requires information about the material moments only for the cell under consideration. The adaptive method chooses between shapes obtained by the intersection of the cell with one half-plane, two half-planes, or a circle. The A-MOF method allows to exactly reproduce several convex shapes: corners, filaments, and their concave cell-complements; as well as pieces of the circles and its cell-compliments. Interface reconstruction is formulated as a local (for each cell), non-linear, equality constrained optimization problem, which does not require additional communication and allows for an efficient parallel implementation. In conclusion, we present an extensive set of test problems, both for interface reconstruction on a single cell, and for reconstruction of a variety of shapes on the entire mesh.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Theoretical and experimental study on electron interactions with chlorobenzene: Shape resonances and differential cross sections

In this work, we report theoretical and experimental cross sections for elastic scattering of electrons by chlorobenzene (ClB). The theoretical integral and differential cross sections (DCSs) were obtained with the Schwinger multichannel method implemented with pseudopotentials (SMCPP) and the independent atom method with screening corrected additivity rule (IAM-SCAR). The calculations with the SMCPP method were done in the static-exchange (SE) approximation, for energies above 12 eV, and in the static-exchange plus polarization approximation, for energies up to 12 eV. The calculations with the IAM-SCAR method covered energies up to 500 eV. The experimental differential cross sections were obtained in the high resolution electron energy loss spectrometer VG-SEELS 400, in Lisbon, for electron energies from 8.0 eV to 50 eV and angular range from 7{sup ∘} to 110{sup ∘}. From the present theoretical integral cross section (ICS) we discuss the low-energy shape-resonances present in chlorobenzene and compare our computed resonance spectra with available electron transmission spectroscopy data present in the literature. Since there is no other work in the literature reporting differential cross sections for this molecule, we compare our theoretical and experimental DCSs with experimental data available for the parent molecule benzene.

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

Mutual information-assisted adaptive variational quantum eigensolver

Adaptive construction of ansatz circuits offers a promising route towards applicable variational quantum eigensolvers on near-term quantum hardware. Those algorithms aim to build up optimal circuits for a certain problem and ansatz circuits are adaptively constructed by selecting and adding entanglers from a predefined pool. In this work, we propose a way to construct entangler pools with reduced size by leveraging classical algorithms. Our method uses mutual information between the qubits in classically approximated ground state to rank and screen the entanglers. The density matrix renormalization group method is employed for classical precomputation in this work. We corroborate our method numerically on small molecules. Our numerical experiments show that a reduced entangler pool with a small portion of the original entangler pool can achieve same numerical accuracy. Here, we believe that our method paves a new way for adaptive construction of ansatz circuits for variational quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

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

Mechanical properties of Zircaloy cladding tubes and contributions to M.E.T.A. mechanical property database

To support a multi-laboratory Methodology, Evaluation, Testing, and Analysis (M.E.T.A.) cladding properties database, Oak Ridge National Laboratory’s (ORNL’s) cladding mechanical test geometries were manufactured from several nuclear-relevant cladding alloys and subsequently tested. These geometries were developed as mechanical test specimens to evaluate the properties of tube materials that may be used for irradiation testing at ORNL’s High Flux Isotope Reactor. They may also be used as test articles to be harvested—via in-cell machining—from commercially irradiated fuel rods and later tested. This report explores the differences among axial, hoop, and SSJ tensile geometries with partially recrystallized Zircaloy-2 to test ORNL correlation-based methods on a plate material that approximates, to the greatest extent possible, the characteristics of nuclear industry tubing. Furthermore, several tests were conducted with ORNL’s Zircaloy-4 tube inventory to (1) develop material properties as a standard for future tests, (2) determine the effect of the US Department of Energy’s Advanced Fuels Campaign coating processes on tube mechanical properties, and (3) evaluate the effect of specimen machining methods on the mechanical properties of tube geometries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

A weighted ℓ{sub 1}-minimization approach for sparse polynomial chaos expansions

This work proposes a method for sparse polynomial chaos (PC) approximation of high-dimensional stochastic functions based on non-adapted random sampling. We modify the standard ℓ{sub 1}-minimization algorithm, originally proposed in the context of compressive sampling, using a priori information about the decay of the PC coefficients, when available, and refer to the resulting algorithm as weightedℓ{sub 1}-minimization. We provide conditions under which we may guarantee recovery using this weighted scheme. Numerical tests are used to compare the weighted and non-weighted methods for the recovery of solutions to two differential equations with high-dimensional random inputs: a boundary value problem with a random elliptic operator and a 2-D thermally driven cavity flow with random boundary condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

From Bessel beam to complex-source-point cylindrical wave-function

This investigation shows that a scalar Bessel beam can be transformed into the non-paraxial complex-source-point cylindrical wave (CSPCW). High-order CSPCW solutions, termed here high-order quasi-Gaussian cylindrical beams, which exactly satisfy the Helmholtz equation, are derived analytically. Moreover, partial-derivatives of the high-order CSPCW solutions satisfy the Helmholtz equation. In addition, the CSPCW solutions satisfy the nonrelativistic Schrödinger equation within standard quantum mechanics, thus, the results can be used in the description of elementary particle/matter motion and related applications in quantum scattering theory. Furthermore, the analysis is extended to the case of vector beams in which the components of the electromagnetic (EM) field are obtained based on different polarizations of the magnetic and electric vector potentials, which exactly satisfy Maxwell’s vectorial equations and Lorenz’ gauge condition. An attractive feature of the high-order solutions is the rigorous description of strongly focused (or strongly divergent) cylindrical wave-fields without any approximations, nor the need for numerical methods. Possible applications are in beam-forming design using high-aperture or collimated cylindrical laser/electron quasi-Gaussian beams in imaging microscopy, particle manipulation, optical tweezers, and the study of the scattering, and radiation forces on objects. - Highlights: • Bessel beam is transformed into the non-paraxial cylindrical complex-source-point. • Exact high-order tightly focused solutions are derived without any approximations. • The exact solutions also satisfy the nonrelativistic Schrödinger equation. • Electromagnetic beams are obtained as solutions of Maxwell’s vectorial equations. • Applications are in laser/electron beam imaging, tweezers, and radiation force.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Saddle point approaches in strong field physics and generation of attosecond pulses

Attoscience is the emerging field that accesses the fastest electronic processes occurring at the atomic and molecular length scales with attosecond (1 as = 10 –18 s) time resolution having wide ranging physical, chemical, material science and biological applications. The quintessential and one of the most fundamental processes in this domain is the generation of phase locked XUV attosecond pulses. The theoretical approach to understand the process incorporates a fully quantum or semi classical or relativistic description of coherent charge dynamics in intense ultrashort electromagnetic fields driving a quantum system (an atom, a molecule, solid band gap materials or surface plasmas). Modelling of such physical and dynamical systems in science and also in many other branches often leads to equations represented in terms of complex multi-dimensional integrals. These integrals can often be solved using the stationary phase approximation, which leads to a series of equations identifying the points in the multi-dimensional space, having most significant contributions in their evaluation. These points are usually indicated as saddle points. The description of the dynamics of quantum mechanical or relativistic systems that results from such an approach enables near to classical physics intuitive perceptions of the processes under investigation. Thus, the saddle point methods are very powerful and valuable general theoretical tools to obtain asymptotic expressions of such solutions and help also to gain physical insights on the underlying phenomena. Such techniques developed in the past have been adapted to study the emission of as pulses by different physical systems and have been widely employed in calculating and estimating the response of matter to intense electromagnetic pulses on ultrafast time scales. Here we provide an extensive disposition of the saddle point approaches unifying their ubiquitous applications within the domain of attoscience valid for simple atomic to more complex condensed matter systems undergoing ultrafast dynamics and present current trends and advancements in the field. In this review we would delineate the methodology, present a synthesis of seminal works and describe the state of the art applications. Finally we also address ultrashort time dynamics of novel materials that have gained much attention recently, namely lower dimensional material systems and micro-plasma systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Low-Rank QTT-based Finite Element Method for Elasticity Problems

We present an efficient and robust numerical algorithm for solving the linear elasticity problem that combines the Quantized Tensor Train format and a domain partitioning strategy. This approach makes it possible to solve the linear elasticity problem on a computational domain that is more general than a square. By integrating Z-ordering and subdomain concatenation, our method substantially decreases memory usage and achieves a notable reduction in rank compared to established Finite Element implementations like the FEniCS platform. This efficiency is maintained while still guaranteeing exponential convergence with respect to the number of degrees of freedom. This performance gain, however, requires a fundamental rethinking of how core finite element operations are implemented. This includes changes to mesh discretization, node and degree of freedom ordering, stiffness matrix and internal nodal force assembly, and the execution of algebraic matrix-vector operations. In this work, we discuss all these aspects in detail and assess the method’s performance in the numerical approximation of three representative test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Simple Models of Lateral Heterostructures

General analytical expressions for densities of states of a lateral heterostructure, formed by the contact of two square semi-infinite lattices with single-band and two-band spectra, were obtained in the tight-binding approximation by the Green’s function method. The semi-elliptical density of states was used for numerical estimates, and the model of two interacting dimers was proposed to estimate the charge transfer. Application of this approach to description of lateral epitaxial and graphene-like heterostructures is discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Analysis of the shifted boundary method for the Stokes problem

The analysis of stability and accuracy of the shifted boundary method is developed for the Stokes flow equations. The key feature of the shifted boundary method, an embedded finite element method, is the shifting of the location where boundary conditions are applied from the true to a surrogate boundary, and an appropriate modification (shifting) of the value of the boundary conditions. An inf–sup condition is proved for the variational formulation associated to the shifted boundary method and we derive, by way of Strang’s second lemma, an optimal error estimate in the natural SBM norm. Here we also derive an L 2 -error estimate for the velocity field, by means of an extension of the Aubin–Nitsche approach to embedded, non-consistent, mixed finite element methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Verification and Large Deformation Analysis Using the Reproducing Kernel Particle Method

The reproducing kernel particle method (RKPM) is a meshless method used to solve general boundary value problems using the principle of virtual work. RKPM corrects the kernel approximation by introducing reproducing conditions which force the method to be complete to arbritrary order polynomials selected by the user. Effort in recent years has led to the implementation of RKPM within the Sierra/SM physics software framework. The purpose of this report is to investigate convergence of RKPM for verification and validation purposes as well as to demonstrate the large deformation capability of RKPM in problems where the finite element method is known to experience difficulty. Results from analyses using RKPM are compared against finite element analysis. A host of issues associated with RKPM are identified and a number of potential improvements are discussed for future work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Comparing machine learning and interpolation methods for loop-level calculations

The need to approximate functions is ubiquitous in science, either due to empirical constraints or high computational cost of accessing the function. In high-energy physics, the precise computation of the scattering cross-section of a process requires the evaluation of computationally intensive integrals. A wide variety of methods in machine learning have been used to tackle this problem, but often the motivation of using one method over another is lacking. Comparing these methods is typically highly dependent on the problem at hand, so we specify to the case where we can evaluate the function a large number of times, after which quick and accurate evaluation can take place. We consider four interpolation and three machine learning techniques and compare their performance on three toy functions, the four-point scalar Passarino-Veltman D_0 D 0 function, and the two-loop self-energy master integral M. We find that in low dimensions (d = 3), traditional interpolation techniques like the Radial Basis Function perform very well, but in higher dimensions (d=5, 6, 9) we find that multi-layer perceptrons (a.k.a neural networks) do not suffer as much from the curse of dimensionality and provide the fastest and most accurate predictions.

97 MATHEMATICS AND COMPUTING↗