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

Testing the Perey effect

Here, the effects of non-local potentials have historically been approximately included by applying a correction factor to the solution of the corresponding equation for the local equivalent interaction. This is usually referred to as the Perey correction factor. In this work we investigate the validity of the Perey correction factor for single-channel bound and scattering states, as well as in transfer (p, d) cross sections. Method: We solve the scattering and bound state equations for non-local interactions of the Perey-Buck type, through an iterative method. Using the distorted wave Born approximation, we construct the T-matrix for (p,d) on 17 O, 41 Ca, 49 Ca, 127 Sn, 133 Sn, and 209 Pb at 20 and 50 MeV. As a result, we found that for bound states, the Perey corrected wave function resulting from the local equation agreed well with that from the non-local equation in the interior region, but discrepancies were found in the surface and peripheral regions. Overall, the Perey correction factor was adequate for scattering states, with the exception of a few partial waves corresponding to the grazing impact parameters. These differences proved to be important for transfer reactions. In conclusion, the Perey correction factor does offer an improvement over taking a direct local equivalent solution. However, if the desired accuracy is to be better than 10%, the exact solution of the non-local equation should be pursued.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2012

Linear source approximation in CASMO5

A Linear Source (LS) approximation has been implemented in the two-dimensional Method of Characteristics (MOC) transport solver in a prototype version of CASMO5. The LS approximation, which relies on the computation of trajectory-based spatial moments over source regions to obtain the linear source expansion coefficients, improves the solution accuracy relative to the 'flat' or constant source approximation. In addition, the LS formulation is capable of treating arbitrarily-shaped source regions and is compatible with standard Coarse-Mesh Finite Difference (CMFD) acceleration. Numerical tests presented in this paper for the C5G7 MOX benchmark show that, for comparable accuracy with respect to the reference solution, the LS approximation can reduce the run time by a factor of four and the memory requirements by a factor often relative to the FS scheme. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2012

Lead Slowing-Down Spectrometry Time Spectral Analysis for Spent Fuel Assay: FY12 Status Report

Executive Summary Developing a method for the accurate, direct, and independent assay of the fissile isotopes in bulk materials (such as used fuel) from next-generation domestic nuclear fuel cycles is a goal of the Office of Nuclear Energy, Fuel Cycle R&D, Material Protection and Control Technology (MPACT) Campaign. To meet this goal, MPACT supports a multi-institutional collaboration, of which PNNL is a part, to study the feasibility of Lead Slowing Down Spectroscopy (LSDS). This technique is an active nondestructive assay method that has the potential to provide independent, direct measurement of Pu and U isotopic masses in used fuel with an uncertainty considerably lower than the approximately 10% typical of today’s confirmatory methods. This document is a progress report for FY2012 PNNL analysis and algorithm development. Progress made by PNNL in FY2012 continues to indicate the promise of LSDS analysis and algorithms applied to used fuel assemblies. PNNL further refined the semi-empirical model developed in FY2011 based on singular value decomposition (SVD) to numerically account for the effects of self-shielding. The average uncertainty in the Pu mass across the NGSI-64 fuel assemblies was shown to be less than 3% using only six calibration assemblies with a 2% uncertainty in the isotopic masses. When calibrated against the six NGSI-64 fuel assemblies, the algorithm was able to determine the total Pu mass within <2% uncertainty for the 27 diversion cases also developed under NGSI. Two purely empirical algorithms were developed that do not require the use of Pu isotopic fission chambers. The semi-empirical and purely empirical algorithms were successfully tested using MCNPX simulations as well applied to experimental data measured by RPI using their LSDS. The algorithms were able to describe the 235U masses of the RPI measurements with an average uncertainty of 2.3%. Analyses were conducted that provided valuable insight with regard to design requirements (e.g. Pb stack size, neutron source location) of an LSDS for the purpose of assaying used fuel assemblies. Sensitivity studies were conducted that provide insight as to how the LSDS instrument can be improved by making it more sensitive to the center of the fuel assemblies. In FY2013, PNNL will continue efforts to develop and refine design requirements of an LSDS for the ultimate purpose of assaying used fuel assemblies. Future efforts will be directed toward more extensive experimental benchmarking of currently implemented time-spectra analysis algorithms.

Kulisek, Jonathan A.↗
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↗