Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Determination of sin{sup 2}θ{sub W} using ν (ν{sup ¯})-nucleus scattering

We have studied nonisoscalarity and medium effects in the extraction of weak mixing angle using Paschos and Wolfenstein relation in the iron nucleus. Paschos and Wolfenstein(PW) relation is valid for an isoscalar target. We have modified the PW relation for nonisoscalar target as well as incorporated the medium effects like Pauli blocking, Fermi motion, nuclear binding energy and pion rho cloud contributions. In our calculations we have used the relativistic nuclear spectral function which includes nucleon correlations. Finally local density approximation is applied to translate the numerical results to the finite nuclei. We have studied the dependence of sin{sup 2}θ{sub W} on Bjorken variables x and y, four momentum transfer square (Q{sup 2}), energy of the neutrino and antineutrino, and effect of excess neutrons over protons in the nuclear target.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

The rates and time-delay distribution of multiply imaged supernovae behind lensing clusters

Time delays of gravitationally lensed sources can be used to constrain the mass model of a deflector and determine cosmological parameters. We here present an analysis of the time-delay distribution of multiply imaged sources behind 17 strong lensing galaxy clusters with well-calibrated mass models. We find that for time delays less than 1000 days, at z = 3.0, their logarithmic probability distribution functions are well represented by P(log Δt) = 5.3 × 10{sup −4}Δt{sup β-tilde}/M{sub 250}{sup 2β-tilde}, with β-tilde = 0.77, where M{sub 250} is the projected cluster mass inside 250 kpc (in 10{sup 14}M{sub ☉}), and β-tilde is the power-law slope of the distribution. The resultant probability distribution function enables us to estimate the time-delay distribution in a lensing cluster of known mass. For a cluster with M{sub 250} = 2 × 10{sup 14}M{sub ☉}, the fraction of time delays less than 1000 days is approximately 3%. Taking Abell 1689 as an example, its dark halo and brightest galaxies, with central velocity dispersions σ≥500kms{sup −1}, mainly produce large time delays, while galaxy-scale mass clumps are responsible for generating smaller time delays. We estimate the probability of observing multiple images of a supernova in the known images of Abell 1689. A two-component model of estimating the supernova rate is applied in this work. For a magnitude threshold of m{sub AB} = 26.5, the yearly rate of Type Ia (core-collapse) supernovae with time delays less than 1000 days is 0.004±0.002 (0.029±0.001). If the magnitude threshold is lowered to m{sub AB} ∼ 27.0, the rate of core-collapse supernovae suitable for time delay observation is 0.044±0.015 per year.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2015

Variable Neural Adaptive Robust Control: A Switched System Approach

Variable neural adaptive robust control strategies are proposed for the output tracking control of a class of multi-input multi-output uncertain systems. The controllers incorporate a variable-structure radial basis function (RBF) network as the self-organizing approximator for unknown system dynamics. The variable-structure RBF network solves the problem of structure determination associated with fixed-structure RBF networks. It can determine the network structure on-line dynamically by adding or removing radial basis functions according to the tracking performance. The structure variation is taken into account in the stability analysis of the closed-loop system using a switched system approach with the aid of the piecewise quadratic Lyapunov function. The performance of the proposed variable neural adaptive robust controllers is illustrated with simulations.

adaptive↗
DOE OSTI2011

Unified BCS-like model of pairing and alpha-correlations

We construct a BCS-like model that combines nucleonic pairing correlations and possible quartic correlations of alpha-type in a single variational wave function and derive corresponding gap equations. In the approximation of large logarithms typical for the BCS approach, we show that the system reveals two possible types of a condensate which cannot coexist. If the alpha-type condensate prevails at N = Z, the growth of the neutron excess will naturally lead to the first order phase transition to the nucleon condensates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗
DOE OSTI2014

Bound state solution of Dirac equation for Hulthen plus trigonometric Rosen Morse non-central potential using Romanovski polynomial

The bound state solutions of Dirac equation for Hulthen and trigonometric Rosen Morse non-central potential are obtained using finite Romanovski polynomials. The approximate relativistic energy spectrum and the radial wave functions which are given in terms of Romanovski polynomials are obtained from solution of radial Dirac equation. The angular wave functions and the orbital quantum number are found from angular Dirac equation solution. In non-relativistic limit, the relativistic energy spectrum reduces into non-relativistic energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Morse–Smale Regression

This paper introduces a novel partition-based regression approach that incorporates topological information. Partition-based regression typically introduces a quality-of-fit-driven decomposition of the domain. The emphasis in this work is on a topologically meaningful segmentation. Thus, the proposed regression approach is based on a segmentation induced by a discrete approximation of the Morse–Smale complex. This yields a segmentation with partitions corresponding to regions of the function with a single minimum and maximum that are often well approximated by a linear model. This approach yields regression models that are amenable to interpretation and have good predictive capacity. Typically, regression estimates are quantified by their geometrical accuracy. For the proposed regression, an important aspect is the quality of the segmentation itself. Thus, this article introduces a new criterion that measures the topological accuracy of the estimate. The topological accuracy provides a complementary measure to the classical geometrical error measures and is very sensitive to overfitting. The Morse–Smale regression is compared to state-of-the-art approaches in terms of geometry and topology and yields comparable or improved fits in many cases. Finally, a detailed study on climate-simulation data demonstrates the application of the Morse–Smale regression. Supplementary Materials are available online and contain an implementation of the proposed approach in the R package msr, an analysis and simulations on the stability of the Morse–Smale complex approximation, and additional tables for the climate-simulation study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Stochastic quasi-steady state approximations for asymptotic solutions of the chemical master equation

In this paper, we propose two methods to carry out the quasi-steady state approximation in stochastic models of enzyme catalytic regulation, based on WKB asymptotics of the chemical master equation or of the corresponding partial differential equation for the generating function. The first of the methods we propose involves the development of multiscale generalisation of a WKB approximation of the solution of the master equation, where the separation of time scales is made explicit which allows us to apply the quasi-steady state approximation in a straightforward manner. To the lowest order, the multi-scale WKB method provides a quasi-steady state, Gaussian approximation of the probability distribution. The second method is based on the Hamilton-Jacobi representation of the stochastic process where, as predicted by large deviation theory, the solution of the partial differential equation for the corresponding characteristic function is given in terms of an effective action functional. The optimal transition paths between two states are then given by those paths that maximise the effective action. Such paths are the solutions of the Hamilton equations for the Hamiltonian associated to the effective action functional. The quasi-steady state approximation is applied to the Hamilton equations thus providing an approximation to the optimal transition paths and the transition time between two states. Using this approximation we predict that, unlike the mean-field quasi-steady approximation result, the rate of enzyme catalysis depends explicitly on the initial number of enzyme molecules. The accuracy and validity of our approximated results as well as that of our predictions regarding the behaviour of the stochastic enzyme catalytic models are verified by direct simulation of the stochastic model using Gillespie stochastic simulation algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Rotating quantum turbulence in the unitary Fermi gas

Quantized vortices carry the angular momentum in rotating superfluids, and are key to the phenomenon of quantum turbulence. Advances in ultracold-atom technology enable quantum turbulence to be studied in regimes with both experimental and theoretical control, unlike the original contexts of superfluid helium experiments. While much work has been performed with bosonic systems, detailed studies of fermionic quantum turbulence are nascent, despite wide applicability to other contexts such as rotating neutron stars. In this paper, we present a large-scale study of quantum turbulence in rotating fermionic superfluids using an accurate time-dependent density functional theory called the superfluid local density approximation. We identify two different modes of turbulent decay in the dynamical equilibration of a rotating fermionic superfluid, and contrast these results with a computationally simpler description provided by the Gross-Pitaevskii equation, which we find can qualitatively reproduce these decay mechanisms if dissipation is explicitly included. These results demonstrate that dissipation mechanisms intrinsic to fermionic superfluids play a key role in differentiating fermionic from bosonic turbulence, which manifests by enhanced damping of Kelvin waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Phonon and Thermodynamic Properties in YNi{sub 2}B{sub 2}C and LuNi{sub 2}B{sub 2}C from First-Principles Study

The thermodynamic properties of the tetragonal RNi{sub 2}B{sub 2}C (R=Y, Lu) are investigated by means of the density-functional theory (DFT) and density-functional perturbation theory (DFPT) within the quasiharmonic approximation (QHA). The phonon dispersion curves and phonon total and partial density of states are studied, and it is found that the phonon branches can be divided into three obvious regions: a low-frequency region (LFR), an intermediate-frequency region (IFR), and a high-frequency region (HFR). The temperature dependence of the various quantities including isothermal bulk modulus, thermal expansion coefficient, and heat capacity is obtained, and the variation features of these quantities with temperature are discussed in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Frequency dispersion of the first hyperpolarizabilities of reference molecules for nonlinear optics

The frequency dispersion of the hyper-Rayleigh scattering first hyperpolarizabilities (β{sub HRS}) of five reference molecules for nonlinear optics, namely, carbon tetrachloride, chloroform, dichloromethane, acetonitrile, and trichloroacetonitrile, is described using the coupled-cluster singles and doubles quadratic response function (CCSD-QRF) as well as approximate schemes. Comparisons to approximate schemes in which the frequency dispersion is evaluated as either a multiplicative or an additive correction to the static hyperpolarizability yield the following observations: (i) errors of the order of 10% or less are usually encountered when using the multiplicative scheme for photon energies far from the lowest dipole-allowed excitation energies, (ii) spurious cases cannot be excluded as evidenced by carbon tetrachloride where the multiplicative scheme predicts a decrease of β{sub HRS} in contradiction to the increase obtained using the CCSD-QRF method, and (iii) the additive scheme is at best as reliable as the multiplicative approximation. The two-state approximation presents the advantage of correcting the wrong behavior of the additive and multiplicative schemes for carbon tetrachloride, but it is not an improved solution for the other compounds, while the question of selecting the appropriate dominant excited state remains unanswered. Finally, a new β{sub xyz} value of 18.9 a.u. is proposed for carbon tetrachloride in gas phase at λ = 1064 nm, to be compared with the measured 16.9 ± 1.4 a.u. value due to Shelton.

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

Generalized Enhanced Multivariance Product Representation for Data Partitioning: Constancy Level

Enhanced Multivariance Product Representation (EMPR) method is used to represent multivariate functions in terms of less-variate structures. The EMPR method extends the HDMR expansion by inserting some additional support functions to increase the quality of the approximants obtained for dominantly or purely multiplicative analytical structures. This work aims to develop the generalized form of the EMPR method to be used in multivariate data partitioning approaches. For this purpose, the Generalized HDMR philosophy is taken into consideration to construct the details of the Generalized EMPR at constancy level as the introductory steps and encouraging results are obtained in data partitioning problems by using our new method. In addition, to examine this performance, a number of numerical implementations with concluding remarks are given at the end of this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Communication: Testing and using the Lewin-Lieb bounds in density functional theory

Lewin and Lieb have recently proven several new bounds on the exchange-correlation energy that complement the Lieb-Oxford bound. We test these bounds for atoms, for slowly-varying gases, and for Hooke’s atom, finding them usually less strict than the Lieb-Oxford bound. However, we also show that, if a generalized gradient approximation is to guarantee satisfaction of the new bounds for all densities, new restrictions on the exchange-correlation enhancement factor are implied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

LyMAS: Predicting large-scale Lyα forest statistics from the dark matter density field

We describe Lyα Mass Association Scheme (LyMAS), a method of predicting clustering statistics in the Lyα forest on large scales from moderate-resolution simulations of the dark matter (DM) distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the 'Horizon-MareNostrum' simulation, a 50 h {sup –1} Mpc comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F{sub s} |δ {sub s}) of the transmitted flux F{sub s} , smoothed (one-dimensionally, 1D) over the spectral resolution scale, on the DM density contrast δ {sub s}, smoothed (three-dimensionally, 3D) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at z = 2.5, and we find optimal results for a DM smoothing length σ = 0.3 h {sup –1} Mpc (comoving). In its simplest form, LyMAS draws randomly from the hydro-calibrated P(F{sub s} |δ {sub s}) to convert DM skewers into Lyα forest pseudo-spectra, which are then used to compute cross-sightline flux statistics. In extended form, LyMAS exactly reproduces both the 1D power spectrum and one-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum DM field, LyMAS accurately predicts the two-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 h {sup –1} Mpc, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ({sup F}luctuating Gunn-Peterson Approximation{sup )}, often used for large-volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024{sup 3} N-body simulations of a 300 h {sup –1} Mpc and 1.0 h {sup –1} Gpc cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3D Lyα forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, DM, space geometry, and intergalactic medium physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Photoionization and photofragmentation of singly charged positive and negative Sc 3 N @ C 80 endohedral fullerene ions

Photoprocesses of the endohedral fullerene ions $\mathrm{Sc_3N@C_{80}^{+}}$ and $\mathrm{Sc_3N@C_{80}-}$ in the gas phase have been investigated in the photon energy ranges 30-50 eV and 280-420 eV. Single and double ionization as well as single ionization accompanied by the release of a C 2 dimer were observed as a function of the photon energy for the positive parent ion and double detachment was measured for the negative parent ion. The emphasis of the experiments was on the specific effects of the encapsulated trimetallic nitride cluster Sc 3 N on the observed reactions. Clear evidence of photoexcitation near the Sc L edge is obtained with the dominating contributions visible in the one- and two-electron-removal channels. K-vacancy production in the encapsulated central nitrogen atom is seen in the single ionization of $\mathrm{Sc_3N@C_{80}^{+}}$ but is much less pronounced in the photoionization-with-fragmentation channel. Comparison of the cross sections near the carbon K edge with the corresponding channels measured previously in the photoionization of $\mathrm{Lu_3N@C_{80}^+}$ reveal strong similarities. Previously predicted sharp resonance features in the ionization of $\mathrm{Sc_3N@C_{80}^{+}}$ ions below the Sc M edge are not confirmed. Finally, the experiments are accompanied by quantum-chemistry calculations in the Hartree-Fock approximation and by model calculations employing density functional theory (DFT).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Stable half-metallic monolayers of FeCl{sub 2}

The structural, electronic, and magnetic properties of single layers of Iron Dichloride (FeCl{sub 2}) were calculated using first principles calculations. We found that the 1T phase of the single layer FeCl{sub 2} is 0.17 eV/unit cell more favorable than its 1H phase. The structural stability is confirmed by phonon calculations. We found that 1T-FeCl{sub 2} possess three Raman-active (130, 179, and 237 cm{sup −1}) and one infrared-active (279 cm{sup −1}) phonon branches. The electronic band dispersion of the 1T-FeCl{sub 2} is calculated using both gradient approximation of Perdew-Burke-Ernzerhof and DFT-HSE06 functionals. Both functionals reveal that the 1T-FeCl{sub 2} has a half-metallic ground state with a Curie temperature of 17 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Epitaxial Carbyne: Analytical Results

The Green’s function method in the tight-binding approximation is used to obtain analytical expressions for the electronic spectra and densities of states for two carbyne structural modifications (cumulene and polyyne). Metal and semiconductor substrates are considered and charge-transfer estimates are deduced. Criteria for the occurrence of a charge-density wave in free-standing and epitaxial carbynes are proposed taking into account the intra- and interatomic Coulomb repulsion. Analytical expressions for the phonon spectrum of free-standing carbyne are presented. The data obtained are compared with the results of numerical calculations performed by other researchers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗