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

Theoretical study of the relativistic molecular rotational g-tensor

An original formulation of the relativistic molecular rotational g-tensor valid for heavy atom containing compounds is presented. In such formulation, the relevant terms of a molecular Hamiltonian for non-relativistic nuclei and relativistic electrons in the laboratory system are considered. Terms linear and bilinear in the nuclear rotation angular momentum and an external uniform magnetic field are considered within first and second order (relativistic) perturbation theory to obtain the rotational g-tensor. Relativistic effects are further analyzed by carrying out the linear response within the elimination of the small component expansion. Quantitative results for model systems HX (X=F, Cl, Br, I), XF (X=Cl, Br, I), and YH{sup +} (Y=Ne, Ar, Kr, Xe, Rn) are obtained both at the RPA and density functional theory levels of approximation. Relativistic effects are shown to be small for this molecular property. The relation between the rotational g-tensor and susceptibility tensor which is valid in the non-relativistic theory does not hold within the relativistic framework, and differences between both molecular parameters are analyzed for the model systems under study. It is found that the non-relativistic relation remains valid within 2% even for the heavy HI, IF, and XeH{sup +} systems. Only for the sixth-row Rn atom a significant deviation of this relation is found.

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

Generation, dynamics, and correlations of the fission fragments' angular momenta

The generation of angular momentum in fissioning nuclei is not well understood. The predictions of different models disagree, particularly concerning the correlation between the fragments' angular momenta. In this article, a time-dependent collective Hamiltonian model is proposed to treat the generation of the angular momentum in the fission fragments due to the quantum uncertainty principle as well as the dynamics of the collective wave function during and after scission. The model is constructed in the framework of the frozen Hartree-Fock approximation using a Skyrme energy functional to extract deformations of the fission fragments as well as the interactions in a derived collective Hamiltonian. The fission reactions studied are $^{240}$Pu $\rightarrow$ $^{132}$Sn+$^{108}$Ru and $^{240}$Pu $\rightarrow$ $^{144}$Ba+$^{96}$Sr. The model can account for a large part of the angular momentum found in experimental data. The orientation of the angular momentum of each fragment is found to be mainly in the plane perpendicular to the fission axis, in agreement with the experiment. The magnitudes of the angular momenta in the two fragments are nearly uncorrelated, in agreement with the recent experimental data of Wilson et al., Nature (London) 590, 566 (2021). Some of the conclusions of the traditional collective vibration model are supported by the present model but some are not. Surprisingly, it is found that the angular momenta of the fragments are slightly correlated positively as in a wriggling mode. As a result, it is also found that the presence of an octupole deformation in a fragment can significantly increase the generated angular momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Results for aliovalent doping of CeBr{sub 3} with Ca{sup 2+}

Despite the outstanding scintillation performance characteristics of cerium tribromide (CeBr{sub 3}) and cerium-activated lanthanum tribromide, their commercial availability and application are limited due to the difficulties of growing large, crack-free single crystals from these fragile materials. This investigation employed aliovalent doping to increase crystal strength while maintaining the optical properties of the crystal. One divalent dopant (Ca{sup 2+}) was used as a dopant to strengthen CeBr{sub 3} without negatively impacting scintillation performance. Ingots containing nominal concentrations of 1.9% of the Ca{sup 2+} dopant were grown, i.e., 1.9% of the CeBr{sub 3} molecules were replaced by CaBr{sub 2} molecules, to match our target replacement of 1 out of 54 cerium atoms be replaced by a calcium atom. Precisely the mixture was composed of 2.26 g of CaBr{sub 2} added to 222.14 g of CeBr{sub 3}. Preliminary scintillation measurements are presented for this aliovalently doped scintillator. Ca{sup 2+}-doped CeBr{sub 3} exhibited little or no change in the peak fluorescence emission for 371 nm optical excitation for CeBr{sub 3}. The structural, electronic, and optical properties of CeBr{sub 3} crystals were studied using the density functional theory within the generalized gradient approximation. Calculated lattice parameters are in agreement with the experimental data. The energy band structures and density of states were obtained. The optical properties of CeBr{sub 3}, including the dielectric function, were calculated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Semiclassical analysis of the Wigner 12j symbol with one small angular momentum

We derive an asymptotic formula for the Wigner 12j symbol, in the limit of one small and 11 large angular momenta. There are two kinds of asymptotic formulas for the 12j symbol with one small angular momentum. We present the first kind of formula in this paper. Our derivation relies on the techniques developed in the semiclassical analysis of the Wigner 9j symbol [L. Yu and R. G. Littlejohn, Phys. Rev. A 83, 052114 (2011)], where we used a gauge-invariant form of the multicomponent WKB wave functions to derive asymptotic formulas for the 9j symbol with small and large angular momenta. When applying the same technique to the 12j symbol in this paper, we find that the spinor is diagonalized in the direction of an intermediate angular momentum. In addition, we find that the geometry of the derived asymptotic formula for the 12j symbol is expressed in terms of the vector diagram for a 9j symbol. This illustrates a general geometric connection between asymptotic limits of the various 3nj symbols. This work contributes an asymptotic formula for the 12j symbol to the quantum theory of angular momentum, and serves as a basis for finding asymptotic formulas for the Wigner 15j symbol with two small angular momenta.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Hydrogen adsorption in metal-organic frameworks: The role of nuclear quantum effects

The role of nuclear quantum effects on the adsorption of molecular hydrogen in metal-organic frameworks (MOFs) has been investigated on grounds of Grand-Canonical Quantized Liquid Density-Functional Theory (GC-QLDFT) calculations. For this purpose, we have carefully validated classical H{sub 2}-host interaction potentials that are obtained by fitting Born-Oppenheimer ab initio reference data. The hydrogen adsorption has first been assessed classically using Liquid Density-Functional Theory and the Grand-Canonical Monte Carlo methods. The results have been compared against the semi-classical treatment of quantum effects by applying the Feynman-Hibbs correction to the Born-Oppenheimer-derived potentials, and by explicit treatment within the GC-QLDFT. The results are compared with experimental data and indicate pronounced quantum and possibly many-particle effects. After validation calculations have been carried out for IRMOF-1 (MOF-5), GC-QLDFT is applied to study the adsorption of H{sub 2} in a series of MOFs, including IRMOF-4, -6, -8, -9, -10, -12, -14, -16, -18, and MOF-177. Finally, we discuss the evolution of the H{sub 2} quantum fluid with increasing pressure and lowering temperature.

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

Emergence of pygmy dipole resonances: Magic numbers and neutron skins

The pygmy dipole resonances (PDR) for even-even nuclei in 8{<=}Z{<=}40 are studied performing a systematic calculation of the random-phase approximation with the Skyrme functional of SkM*. The calculation is fully self-consistent and does not assume any symmetry in the nuclear shape of the ground state. In every isotopic chain, the PDR emerges by showing a peak of the E1 strength at energies less than 10 MeV. The E1 strength of the PDR strongly depends on the position of the Fermi level and shows a clear correlation with the occupation of the orbits with the orbital angular momenta less than 3({h_bar}/2{pi}) (l{<=}2). We also found a strong correlation between the isotopic dependence of the neutron skin thickness and the pygmy dipole strength. The fraction of the energy-weighted strength exhausted by the PDR and the neutron skin thickness show a linear correlation with the universal rate of about 0.2 fm{sup -1}.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Robust Cross-correlation-based Measurement of Clump Sizes in Galaxies

Stars form in molecular complexes that are visible as giant clouds (∼10{sup 5−6} M{sub ⊙}) in nearby galaxies and as giant clumps (∼10{sup 8−9} M{sub ⊙}) in galaxies at redshifts z≈1–3. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g., Hα) can be recovered statistically using the two-point correlation function (2PCF) if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at ∼20% accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical measurements of Hα clump size agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT): Part I. The clear-sky model

The solar energy industry often uses individual steps to empirically compute plane-of-array (POA) irradiance from horizontal irradiance and decompose it to narrow-wavelength bands. Conventional radiative transfer models designed for meteorological applications requires significant computing efforts in practice; however, they provide a physics-based solution of radiance and therefore are capable of computing spectral POA irradiances in a single step. In this study, we integrate the advantages of the current models and develop an innovative radiative transfer model, the Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT), to efficiently compute irradiances on inclined photovoltaics (PV) panels for 2002 narrow-wavelength bands from 0.28 to 4.0 um. This study is reported in two parts. Part I presents the methodology and performance evaluation of the new model under clear-sky conditions. The Simple Model of the Atmospheric Radiative Transfer of Sunshine (SMARTS), which was designed to compute clear-sky irradiances, is employed to rapidly provide the optical properties of a given clear-sky atmosphere. The clear-sky radiances in the narrow-wavelength bands are computed by considering three paths of photon transmission and solving the radiative transfer equation with the single-scattering approximation. The Bi-directional Transmittance Distribution Function (BTDF) of aerosols is given by their single-scattering phase function with a correction using a two-stream approximation. The validation analysis confirms that FARMS-NIT has improved accuracy compared to TMYSPEC as evaluated by both surface observations and a state-of-the-art radiative transfer model. This model substantially improves computational efficiency compared to other radiative transfer models though it uses slightly more computing time than TMYSPEC. Part II of this study addresses the model in cloud-sky conditions and will be published as a companion paper.

14 SOLAR ENERGY↗
DOE OSTI2013

Orientational order and formation of topological defects in two-dimensional systems

The formation dynamics of topological defects and the orientational order in nonideal systems of particles interacting with a screened Coulomb potential is numerically simulated. Calculations are performed in a broad range of parameters corresponding to experimental conditions in a laboratory dusty plasma. A relation is obtained between the number of defects, the shape of the orientational correlation function, and the coupling parameter of the system. New approximations are proposed for the orientational correlation function in the liquid and hexatic phases of structures analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics

Wavefunction extrapolation greatly reduces the number of self-consistent field (SCF) iterations and thus the overall computational cost of Born-Oppenheimer molecular dynamics (BOMD) that is based on the Kohn–Sham density functional theory. Going against the intuition that the higher order of extrapolation possesses a better accuracy, we demonstrate, from both theoretical and numerical perspectives, that the extrapolation accuracy firstly increases and then decreases with respect to the order, and an optimal extrapolation order in terms of minimal number of SCF iterations always exists. We also prove that the optimal order tends to be larger when using larger MD time steps or more strict SCF convergence criteria. By example BOMD simulations of a solid copper system, we show that the optimal extrapolation order covers a broad range when varying the MD time step or the SCF convergence criterion. Therefore, we suggest the necessity for BOMD simulation packages to open the user interface and to provide more choices on the extrapolation order. Another factor that may influence the extrapolation accuracy is the alignment scheme that eliminates the discontinuity in the wavefunctions with respect to the atomic or cell variables. We prove the equivalence between the two existing schemes, thus the implementation of either of them does not lead to essential difference in the extrapolation accuracy.

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

Towards solution and refinement of organic crystal structures by fitting to the atomic pair distribution function

In this study, a method towards the solution and refinement of organic crystal structures by fitting to the atomic pair distribution function (PDF) is developed. Approximate lattice parameters and molecular geometry must be given as input. The molecule is generally treated as a rigid body. The positions and orientations of the molecules inside the unit cell are optimized starting from random values. The PDF is obtained from carefully measured X-ray powder diffraction data. The method resembles `real-space' methods for structure solution from powder data, but works with PDF data instead of the diffraction pattern itself. As such it may be used in situations where the organic compounds are not long-range-ordered, are poorly crystalline, or nanocrystalline. The procedure was applied to solve and refine the crystal structures of quinacridone (β phase), naphthalene and allopurinol. In the case of allopurinol it was even possible to successfully solve and refine the structure in P1 with four independent molecules. As an example of a flexible molecule, the crystal structure of paracetamol was refined using restraints for bond lengths, bond angles and selected torsion angles. In all cases, the resulting structures are in excellent agreement with structures from single-crystal data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Conductive Ink with Circular Life Cycle for Printed Electronics

Electronic waste carries energetic costs and an environmental burden rivaling that of plastic waste due to the rarity and toxicity of the heavy-metal components. Recyclable conductive composites are introduced for printed circuits formulated with polycaprolactone (PCL), conductive fillers, and enzyme/protectant nanoclusters. Circuits can be printed with flexibility (breaking strain ≈80%) and conductivity (≈2.1 × 10 4 S m -1 ). These composites are degraded at the end of life by immersion in warm water with programmable latency. Approximately 94% of the functional fillers can be recycled and reused with similar device performance. The printed circuits remain functional and degradable after shelf storage for at least 7 months at room temperature and one month of continuous operation under electrical voltage. Here, the present studies provide composite design toward recyclable and easily disposable printed electronics for applications such as wearable electronics, biosensors, and soft robotics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A simple finite element method for the Stokes equations

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Effects of La Doping and Zn or O Vacancy on the Magnetic Property of ZnO

Research on the magnetic origin and mechanism of La-doped ZnO system has been controversial. To solve this problem, we studied the effects of La doping and Zn or O vacancy on the magnetic property of ZnO using geometry optimization and energy calculation on the basis of first-principle generalized gradient approximation method of density functional theory. The electronic structures and magnetic properties of undoped ZnO, La-doped ZnO system, and systems with coexisting La doping and Zn or O vacancies were calculated. Results showed that the systems of La-doped ZnO and the coexistence of La doping and O vacancy are nonmagnetic. In addition, the system of La replacing Zn and one Zn vacancy exhibits long-range orderly ferromagnetism, and the Curie temperature of the doping system can achieve room temperature. The magnetism source of the systems with La doping and Zn vacancy coexisting in ZnO demonstrates strong hybrid coupling electron exchange effects existing among Zn-4s, Zn-3p, O-2p, and La-5s orbits that are nearest to the Zn vacancy. The next nearest distance between doping and Zn vacancy leads to the lowest formation energy and highest stability under the same doping condition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effect of Fe Doping and Point Defects (V{sub O} and V{sub Sn}) on the Magnetic Properties of SnO{sub 2}

We investigated the effects of Fe doping and point defects (oxygen and tin vacancies, V{sub O} and V{sub Sn}, respectively) on the ferromagnetism of SnO{sub 2} through first-principles calculations by using generalized gradient approximation + U under density functional theory. Calculation results show that Fe-doped SnO{sub 2} with V{sub O} or V{sub Sn} exhibit long-range ordered ferromagnetism and Curie temperature above room temperature. The system has different magnetic properties with different doping modes, which are beneficial to the magnetic adjustment of dilute magnetic semiconductor (DMS). The double-exchange effects among the electrons of O–2p and Fe–3d states are the origin of ferromagnetism in Sn{sub 15}FeO{sub 31}. Similarly, ferromagnetism in Sn{sub 14}FeO{sub 32} is derived from the double-exchange effects among the electrons of O–2p and Fe–3d states. These effects are produced by hole carriers after complexes form by Fe doping and V{sub O} or V{sub Sn}. Simultaneously, Sn{sub 14}FeO{sub 32} shows a half-metallic behavior, thus resulting in 100% conduction hole polarizability, which is highly beneficial to DMSs as a hole-injection source. Moreover, results also show that the decreased distance between Fe and vacancies (V{sub O} or V{sub Sn}) lowers formation energy and thus increases the system stability.

36 MATERIALS SCIENCE↗