Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Weak entanglement approximation for nuclear structure

The interacting shell model, a configuration-interaction method, is a venerable approach for low-lying nuclear structure calculations, but it is hampered by the exponential growth of its basis dimension as one increases the single-particle space and/or the number of active particles. Recent, quantum-information-inspired work has demonstrated that the proton and neutron sectors of a nuclear wave function are weakly entangled. Furthermore, the entanglement is smaller for nuclides away from N = Z, such as heavy, neutron-rich nuclides. Here, in this study, we implement a weak entanglement approximation to bipartite configuration-interaction wave functions, approximating low-lying levels by coupling a relatively small number of many-proton and many-neutron states. This truncation scheme, which we present in the context of past approaches, reduces the basis dimension by many orders of magnitude while preserving essential features of nuclear spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Statistical model based iterative reconstruction (MBIR) in clinical CT systems. Part II. Experimental assessment of spatial resolution performance

Purpose: Statistical model based iterative reconstruction (MBIR) methods have been introduced to clinical CT systems and are being used in some clinical diagnostic applications. The purpose of this paper is to experimentally assess the unique spatial resolution characteristics of this nonlinear reconstruction method and identify its potential impact on the detectabilities and the associated radiation dose levels for specific imaging tasks. Methods: The thoracic section of a pediatric phantom was repeatedly scanned 50 or 100 times using a 64-slice clinical CT scanner at four different dose levels [CTDI{sub vol} =4, 8, 12, 16 (mGy)]. Both filtered backprojection (FBP) and MBIR (Veo{sup ®}, GE Healthcare, Waukesha, WI) were used for image reconstruction and results were compared with one another. Eight test objects in the phantom with contrast levels ranging from 13 to 1710 HU were used to assess spatial resolution. The axial spatial resolution was quantified with the point spread function (PSF), while the z resolution was quantified with the slice sensitivity profile. Both were measured locally on the test objects and in the image domain. The dependence of spatial resolution on contrast and dose levels was studied. The study also features a systematic investigation of the potential trade-off between spatial resolution and locally defined noise and their joint impact on the overall image quality, which was quantified by the image domain-based channelized Hotelling observer (CHO) detectability index d′. Results: (1) The axial spatial resolution of MBIR depends on both radiation dose level and image contrast level, whereas it is supposedly independent of these two factors in FBP. The axial spatial resolution of MBIR always improved with an increasing radiation dose level and/or contrast level. (2) The axial spatial resolution of MBIR became equivalent to that of FBP at some transitional contrast level, above which MBIR demonstrated superior spatial resolution than FBP (and vice versa); the value of this transitional contrast highly depended on the dose level. (3) The PSFs of MBIR could be approximated as Gaussian functions with reasonably good accuracy. (4) Thez resolution of MBIR showed similar contrast and dose dependence. (5) Noise standard deviation assessed on the edges of objects demonstrated a trade-off with spatial resolution in MBIR. (5) When both spatial resolution and image noise were considered using the CHO analysis, MBIR led to significant improvement in the overall CT image quality for both high and low contrast detection tasks at both standard and low dose levels. Conclusions: Due to the intrinsic nonlinearity of the MBIR method, many well-known CT spatial resolution and noise properties have been modified. In particular, dose dependence and contrast dependence have been introduced to the spatial resolution of CT images by MBIR. The method has also introduced some novel noise-resolution trade-off not seen in traditional CT images. While the benefits of MBIR regarding the overall image quality, as demonstrated in this work, are significant, the optimal use of this method in clinical practice demands a thorough understanding of its unique physical characteristics.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1987

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions.

Bayliss, A.↗
NASA NTRS1989

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions.

Bayliss, A.↗
DOE OSTI2016

Formulation of the Rice-Walsh equation of state based on Shock Hugoniot data for porous metals

The dimensionless material parameter R introduced by Wu and Jing into the Rice-Walsh equation of state (EOS) has been deduced from the LASL shock Hugoniot data for porous Al and Cu. It was found that the parameter R/p decays smoothly with shock pressure p and displays small experimental scatter in the high pressure region. This finding led to the conclusion that the parameter has only a weak temperature dependence and is well approximated by a function of pressure alone, and the Grüneisen parameter should be temperature dependent under compression. The thermodynamic formulation of the Rice-Walsh EOS for Al and Cu was realized using the empirically determined function R(p) for each material and their known shock Hugoniot. It was then possible to reproduce porous shock Hugoniot for these metals. For most degrees of porosity, agreement between the porous data and the calculated Hugoniots using the empirical function described was very good. However, slight discrepancies were seen for Hugoniots with very high porosity. Two new thermal variables were introduced after further analysis, which enabled the calculation of the cold compression curve for these metals. The Grüneisen parameters along full-density and porous Hugoniot curve were calculated using a thermodynamic identity connecting R and the Grüneisen parameter. It was shown that the Grüneisen parameter is strongly temperature dependent. The present analysis suggested that the Rice-Walsh type EOS is a preferable choice for the analysis with its simple form, pressure-dependent empirical Wu-Jing parameter, and its compatibility with porous shock data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Reduced density-matrix functional theory: Correlation and spectroscopy

In this work, we explore the performance of approximations to electron correlation in reduced density-matrix functional theory (RDMFT) and of approximations to the observables calculated within this theory. Our analysis focuses on the calculation of total energies, occupation numbers, removal/addition energies, and spectral functions. We use the exactly solvable Hubbard dimer at 1/4 and 1/2 fillings as test systems. This allows us to analyze the underlying physics and to elucidate the origin of the observed trends. For comparison, we also report the results of the GW approximation, where the self-energy functional is approximated, but no further hypothesis is made concerning the approximations of the observables. In particular, we focus on the atomic limit, where the two sites of the dimer are pulled apart and electrons localize on either site with equal probability, unless a small perturbation is present: this is the regime of strong electron correlation. In this limit, using the Hubbard dimer at 1/2 filling with or without a spin-symmetry-broken ground state allows us to explore how degeneracies and spin-symmetry breaking are treated in RDMFT. We find that, within the used approximations, neither in RDMFT nor in GW, the signature of strong correlation is present, when looking at the removal/addition energies and spectral function from the spin-singlet ground state, whereas both give the exact result for the spin-symmetry broken case. Moreover, we show how the spectroscopic properties change from one spin structure to the other.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Closed-Form Approximation of the Total Variation Proximal Operator

Total variation (TV) is a widely used function for regularizing imaging inverse problems that is particularly appropriate for images whose underlying structure is piecewise constant. TV regularized optimization problems are typically solved using proximal methods, but the way in which they are applied is constrained by the absence of a closed-form expression for the proximal operator of the TV function. A closed-form approximation of the TV proximal operator has previously been proposed, but its accuracy was not theoretically explored in detail. Here, we address this gap by making several new theoretical contributions, proving that the approximation leads to a proximal operator of some convex function, it is equivalent to a gradient descent step on a smoothed version of TV, and that its error can be fully characterized and controlled with its scaling parameter. We experimentally validate our theoretical results on image denoising and sparse-view computed tomography (CT) image reconstruction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2015

EVOLUTION OF FAST MAGNETOACOUSTIC PULSES IN RANDOMLY STRUCTURED CORONAL PLASMAS

We investigate the evolution of fast magnetoacoustic pulses in randomly structured plasmas, in the context of large-scale propagating waves in the solar atmosphere. We perform one-dimensional numerical simulations of fast wave pulses propagating perpendicular to a constant magnetic field in a low-β plasma with a random density profile across the field. Both linear and nonlinear regimes are considered. We study how the evolution of the pulse amplitude and width depends on their initial values and the parameters of the random structuring. Acting as a dispersive medium, a randomly structured plasma causes amplitude attenuation and width broadening of the fast wave pulses. After the passage of the main pulse, secondary propagating and standing fast waves appear. Width evolution of both linear and nonlinear pulses can be well approximated by linear functions; however, narrow pulses may have zero or negative broadening. This arises because narrow pulses are prone to splitting, while broad pulses usually deviate less from their initial Gaussian shape and form ripple structures on top of the main pulse. Linear pulses decay at an almost constant rate, while nonlinear pulses decay exponentially. A pulse interacts most efficiently with a random medium with a correlation length of about half of the initial pulse width. This detailed model of fast wave pulses propagating in highly structured media substantiates the interpretation of EIT waves as fast magnetoacoustic waves. Evolution of a fast pulse provides us with a novel method to diagnose the sub-resolution filamentation of the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A high-order time-parallel scheme for solving wave propagation problems via the direct construction of an approximate time-evolution operator

Our manuscript demonstrates a technique for efficiently solving the classical wave equation, the shallow water equations, and, more generally, equations of the form ∂u/∂t=Lu∂u/∂t=Lu, where LL is a skew-Hermitian differential operator. The idea is to explicitly construct an approximation to the time-evolution operator exp(τL)exp(τL) for a relatively large time-step ττ. Recently developed techniques for approximating oscillatory scalar functions by rational functions, and accelerated algorithms for computing functions of discretized differential operators are exploited. Principal advantages of the proposed method include: stability even for large time-steps, the possibility to parallelize in time over many characteristic wavelengths and large speed-ups over existing methods in situations where simulation over long times are required. Numerical examples involving the 2D rotating shallow water equations and the 2D wave equation in an inhomogenous medium are presented, and the method is compared to the 4th order Runge–Kutta (RK4) method and to the use of Chebyshev polynomials. The new method achieved high accuracy over long-time intervals, and with speeds that are orders of magnitude faster than both RK4 and the use of Chebyshev polynomials.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Spline Classification Methods

The use of spline functions in the development of classification algorithms is discussed. A method is formulated for producing spline approximations to univariate density functions when each density function is described by a histogram of measurements. The resulting approximations are then incorporated into a Bayesian classification procedure for which the probability of misclassification can be readily computed. Some preliminary numerical results are presented to illustrate the method.

Guseman, L. F., Jr.↗
NASA NTRS2006

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗
DOE OSTI2016

Theoretical prediction of the band offsets at the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions using the self-consistent ab initio DFT/GGA-1/2 method

The band offsets of the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions are calculated using the density functional theory/generalized gradient approximation (DFT/GGA)-1/2 method, which takes into account the self-energy corrections and can give an approximate description to the quasiparticle characteristics of the electronic structure of semiconductors. We present the results of the ionization potential (IP)-based and interfacial offset-based band alignments. In the interfacial offset-based band alignment, to get the natural band offset, we use the surface calculations to estimate the change of reference level due to the interfacial strain. Based on the interface models and GGA-1/2 calculations, we find that the valence band maximum and conduction band minimum of ZnO, respectively, lie 0.64 eV and 0.57 eV above those of anatase TiO{sub 2}, while lie 0.84 eV and 1.09 eV below those of GaN, which agree well with the experimental data. However, a large discrepancy exists between the IP-based band offset and the calculated natural band offset, the mechanism of which is discussed. Our results clarify band alignment of the ZnO/anatase TiO{sub 2} heterojunction and show good agreement with the GW calculations for the GaN/ZnO heterojunction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

A Rigorous Framework for Optimization of Expensive Functions by Surrogates

The goal of the research reported here is to develop rigorous optimization algorithms to apply to some engineering design problems for which design application of traditional optimization approaches is not practical. This paper presents and analyzes a framework for generating a sequence of approximations to the objective function and managing the use of these approximations as surrogates for optimization. The result is to obtain convergence to a minimizer of an expensive objective function subject to simple constraints. The approach is widely applicable because it does not require, or even explicitly approximate, derivatives of the objective. Numerical results are presented for a 31-variable helicopter rotor blade design example and for a standard optimization test example.

Booker, Andrew J.↗
DOE OSTI2014

Data approximation using a blending type spline construction

Generalized expo-rational B-splines (GERBS) is a blending type spline construction where local functions at each knot are blended together by C{sup k}-smooth basis functions. One way of approximating discrete regular data using GERBS is by partitioning the data set into subsets and fit a local function to each subset. Partitioning and fitting strategies can be devised such that important or interesting data points are interpolated in order to preserve certain features. We present a method for fitting discrete data using a tensor product GERBS construction. The method is based on detection of feature points using differential geometry. Derivatives, which are necessary for feature point detection and used to construct local surface patches, are approximated from the discrete data using finite differences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2014

Low-lying excited states and nonradiative processes of 9-methyl-2-aminopurine

The UV spectrum of the adenine analogue 9-methyl-2-aminopurine (9M-2AP) is investigated with one- and two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm{sup −1} resolution in a supersonic jet. The electronic origin at 32 252 cm{sup −1} exhibits methyl torsional subbands that originate from the 0A{sub 1}{sup ′′} (l = 0) and 1E{sup ″} (l = ±1) torsional levels. These and further torsional bands that appear up to 0{sub 0}{sup 0}+230 cm{sup −1} allow to fit the threefold (V{sub 3}) barriers of the torsional potentials as |V{sub 3}{sup ′′}|=50 cm{sup −1} in the S{sub 0} and |V{sub 3}{sup ′}|=126 cm{sup −1} in the S{sub 1} state. Using the B3LYP density functional and correlated approximate second-order coupled cluster CC2 methods, the methyl orientation is calculated to be symmetric relative to the 2AP plane in both states, with barriers of V{sub 3}{sup ′′}=20 cm{sup −1} and V{sub 3}{sup ′}=115 cm{sup −1}. The 0{sub 0}{sup 0} rotational band contour is 75% in-plane (a/b) polarized, characteristic for a dominantly long-axis {sup 1}ππ{sup *} excitation. The residual 25% c-axis polarization may indicate coupling of the {sup 1}ππ{sup *} to the close-lying {sup 1}nπ{sup *} state, calculated at 4.00 and 4.01 eV with the CC2 method. However, the CC2 calculated {sup 1}nπ oscillator strength is only 6% of that of the {sup 1}ππ{sup *} transition. The {sup 1}ππ{sup *} vibronic spectrum is very complex, showing about 40 bands within the lowest 500 cm{sup −1}. The methyl torsion and the low-frequency out-of-plane ν{sub 1}{sup ′} and ν{sub 2}{sup ′} vibrations are strongly coupled in the {sup 1}ππ{sup *} state. This gives rise to many torsion-vibration combination bands built on out-of-plane fundamentals, which are without precedence in the {sup 1}ππ{sup *} spectrum of 9H-2-aminopurine [S. Lobsiger, R. K. Sinha, M. Trachsel, and S. Leutwyler, J. Chem. Phys. 134, 114307 (2011)]. From the Lorentzian broadening needed to fit the 0{sub 0}{sup 0} contour of 9M-2AP, the {sup 1}ππ{sup *} lifetime is τ ⩾ 120 ps, reflecting a rapid nonradiative transition.

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

A unified and constructive framework for the universality of neural networks

Abstract One of the reasons why many neural networks are capable of replicating complicated tasks or functions is their universal approximation property. Though the past few decades have seen tremendous advances in theories of neural networks, a single constructive and elementary framework for neural network universality remains unavailable. This paper is an effort to provide a unified and constructive framework for the universality of a large class of activation functions including most of the existing ones. At the heart of the framework is the concept of neural network approximate identity (nAI). The main result is as follows: any nAI activation function is universal in the space of continuous functions on compacta. It turns out that most of the existing activation functions are nAI, and thus universal. The framework induces several advantages over the contemporary counterparts. First, it is constructive with elementary means from functional analysis, probability theory, and numerical analysis. Second, it is one of the first unified and constructive attempts that is valid for most of the existing activation functions. Third, it provides new proofs for most activation functions. Fourth, for a given activation and error tolerance, the framework provides precisely the architecture of the corresponding one-hidden neural network with a predetermined number of neurons and the values of weights/biases. Fifth, the framework allows us to abstractly present the first universal approximation with a favorable non-asymptotic rate. Sixth, our framework also provides insights into the developments, and hence providing constructive derivations, of some of the existing approaches.

97 MATHEMATICS AND COMPUTING↗