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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First-principles study of magnetism and electric field effects in 2D systems

This review article provides a bird's-eye view of what first-principles based methods can contribute to next-generation device design and simulation. After a brief overview of methods and capabilities in the area, the authors focus on published work by their group since 2015 and current work on CrI3. The authors introduce both single- and dual-gate models in the framework of density functional theory and the constrained random phase approximation in estimating the Hubbard U for 2D systems vs their 3D counterparts. A wide range of systems, including graphene-based heterogeneous systems, transition metal dichalcogenides, and topological insulators, and a rich array of physical phenomena, including the macroscopic origin of polarization, field effects on magnetic order, interface state resonance induced peak in transmission coefficients, spin filtration, etc., are covered. For CrI3, the authors present their new results on bilayer systems such as the interplay between stacking and magnetic order, pressure dependence, and electric field induced magnetic phase transitions. The authors find that a bare bilayer CrI3, graphene|bilayer CrI3|graphene, hexagonal boron nitride (h-BN)|bilayer CrI3|h-BN, and h-BN|bilayer CrI3|graphene all have a different response at high field, while at small field, the difference is small except for graphene|bilayer CrI3|graphene. The authors conclude with discussion of some ongoing work and work planned in the near future, with the inclusion of further method development and applications.

Cheng, Hai-Ping (ORCID:0000000159901725)↗
DOE OSTI2022

Hyperspectral imaging of exciton confinement within a moiré unit cell with a subnanometer electron probe

Electronic and optical excitations in two-dimensional systems are distinctly sensitive to the presence of a moiré superlattice. Here we used cryogenic transmission electron microscopy and spectroscopy to simultaneously image the structural reconstruction and associated localization of the lowest-energy intralayer exciton in a rotationally aligned WS 2 -WSe 2 moiré superlattice. In conjunction with optical spectroscopy and ab initio calculations, we determined that the exciton center-of-mass wave function is confined to a radius of approximately 2 nanometers around the highest-energy stacking site in the moiré unit cell. Our results provide direct evidence that atomic reconstructions lead to the strongly confining moiré potentials and that engineering strain at the nanoscale will enable new types of excitonic lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Retroreflection of Optical Radiation from Ellipsoidal-Porous Surfaces

The retroreflection of optical radiation from surfaces with roughness in the form of ellipsoidal pores and two-level roughness in the form of convex ellipsoids with elliptical-porous surfaces is studied using the ROKS-RG Monte-Carlo method software in the geometrical-optics approximation. The dependence of the retroreflection phase function on the parameters of bodies that simulates the rough surface is analyzed. It is shown that the angular dependence of the retroreflection from porous rough surfaces with arbitrarily oriented pores and the aspect ratio higher than unity is a pronounced peak with concave sides at all incidence angles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Ionization-Excitation of helium-like ions by high-energy electrons

The differential and total cross sections are calculated for the ionization of helium-like ions that is accompanied by the excitation of residual ions to ns states owing to an electron impact. Nonrelativistic perturbation theory in electron–electron interaction with Coulomb functions used for a zero-order approximation underlies these calculations. The expressions obtained in this way have a universal character. They are applicable at moderate values of the target charge number Z and high energies of incident electrons. A comparison with total cross sections calculated for the helium atom within various theoretical approaches is performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Kinetics of the Variation in the Magnetic Impurity Ion Concentration in Pb{sub 1–x–y}Sn{sub x}V{sub y}Te Alloys upon Doping

The field and temperature dependences of the magnetization (magnetic fields B ≤ 7.5 T, temperatures T = 2.0–75 K) of samples from a Pb{sub 1–x–y}Sn{sub x}V{sub y}Te (x = 0.08, y = 0.01) single-crystal ingot synthesized by the Bridgman–Stockbarger method. It is established that the sample magnetization contains two main contributions, notably, the paramagnetism of vanadium ions and diamagnetism of the crystal lattice. The field and temperature dependences of the magnetization are approximated by the sum of modified Brillouin functions corresponding to the paramagnetic contributions of vanadium in two different charge states and the diamagnetic contribution linear in terms of field. The concentrations of vanadium ions in two different magnetic states and the character of their variation along the ingot are determined within the scope of the alloy’s electronic-structure rearrangement because of doping. The results are compared with the data of X-ray fluorescence microanalysis and the results of studying the galvanomagnetic properties of the samples.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Electronic Structure of OR-AgInS{sub 2} and TiO{sub 2} Crystals and Its Differential Characteristics

The electronic structures of orthorhombic silver–indium sulfide and rutile titanium oxide were calculated within the framework of density functional theory using the basis of numerical pseudoatomic orbitals in the LDA and GGA approximations. The features of their electron energy spectra and the character of chemical bonds were considered. The comparative analysis of obtained bandgap widths and effective masses determined for heavy electrons and holes in the Brillouin zone center on the basis of analytical derivatives of one-electron energies was performed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Complete lowest order radiative corrections in semi-inclusive scattering of polarized particles

For this work, the lowest order radiative corrections to the cross section and asymmetries measured in experiments on semi-inclusive deep inelastic scattering of polarized particles were calculated. Both exact and leading log expressions were presented and discussed for the total correction that include the contributions from the processes of i) real photon emission with semi-inclusive processes, ii) loop diagrams, and iii) real photon emission with exclusive processes. Radiative corrections to the Sivers and Collins asymmetries in $\pi^+$ electroproduction were studied numerically within the kinematical conditions of modern experimental environments at Jefferson Laboratory (JLab). The Wandzura-Wilczek approximation for the semi-inclusive structure functions and MAID2007 parameterization for the six amplitudes of exclusive processes were used in numeric analyses. The results show that i) radiative effects can generate a correction comparable to the size of Sivers and Collins asymmetries at the Born level, ii) there is good agreement between the exact and leading-order corrections, iii) external functions (that is, other than the Sivers and Collins functions in the respective asymmetries) can generate a contribution to the radiative correction up to 20%, and iv) there exists a strong dependence of the radiative correction on the models for semi-inclusive and exclusive structure functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Electronic stopping and proton dynamics in InP, GaP, and In 0.5 Ga 0.5 P from first principles

The phosphide-based III-V semiconductors InP, GaP, and In 0.5 Ga 0.5 P are promising materials for solar panels in outer space and radioisotope batteries, for which lifetime is a major issue. In order to understand high radiation tolerance of these materials and improve it further, it is necessary to describe the early stages of radiation damage on fast time and short length scales. In particular, the influence of atomic ordering, as observed e.g. in In 0.5 Ga 0.5 P, on electronic stopping is unknown. We use real-time time-dependent density functional theory and the adiabatic local density approximation to simulate electronic stopping of protons in InP, GaP, and the CuAu-I ordered phase of In 0.5 Ga 0.5 P across a large kinetic energy range. These results are compared to SRIM and we investigate the dependence on the channel of the projectile through the target. We show that stopping can be enhanced or reduced in In 0.5 Ga 0.5 P and explain this using the electron-density distribution. By comparing Ehrenfest and Born-Oppenheimer molecular dynamics, we illustrate the 21 intricate dynamics of a proton on a channeling trajectory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

The fully nonconforming virtual element method for biharmonic problems

Here in this paper, we address the numerical approximation of linear fourth-order elliptic problems on polygonal meshes. In particular, we present a novel nonconforming virtual element discretization of arbitrary order of accuracy for biharmonic problems. The approximation space is made of possibly discontinuous functions, thus giving rise to the fully nonconforming virtual element method. We derive optimal error estimates in a suitable (broken) energy norm and present numerical results to assess the validity of the theoretical estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Efficient Band-to-Trap Tunneling Model Including Heterojunction Band Offset

In this paper, we present an efficient band-to-trap tunneling model based on the Schenk approach, in which an analytic density-of-states (DOS) model is developed based on the open boundary scattering method. The new model explicitly includes the effect of heterojunction band offset, in addition to the well-known field effect. Its analytic form enables straightforward implementation into TCAD device simulators. It is applicable to all one-dimensional potentials, which can be approximated to a good degree such that the approximated potentials lead to piecewise analytic wave functions with open boundary conditions. The model allows for simulating both the electric-field-enhanced and band-offset-enhanced carrier recombination due to the band-to-trap tunneling near the heterojunction in a heterojunction bipolar transistor (HBT). Simulation results of an InGaP/GaAs/GaAs NPN HBT show that the proposed model predicts significantly increased base currents, due to the hole-to-trap tunneling enhanced by the emitter-base junction band offset. Finally, the results compare favorably with experimental observation.

42 ENGINEERING↗
DOE OSTI2018

BiEntropy for Python v. 1.0

This Python package provides high-performance implementations of the functions and examples presented in "BiEntropy - The Approximate Entropy of a Finite Binary String" by Grenville J. Croll, presented at ANPA 34 in 2013. https://arxiv.org/abs/1305.0954 According to the paper, BiEntropy is "a simple algorithm which computes the approximate entropy of a finite binary string of arbitrary length" using "a weighted average of the Shannon Entropies of the string and all but the last binary derivative of the string."

Helinski, Ryan↗
DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

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

The EMC Effect in A=3 Nuclei

The European Muon Collaboration (EMC) discovered an unexpected and puzzling result in 1983 when comparing the deep inelastic scattering(DIS) nuclear structure functions of Deuterium and Iron. The per-nucleon structure functions were found to be different for the two nuclei, rather than a simple average over the proton and neutron structure functions. In subsequent experiments, this phenomenon was confirmed for additional nuclei, with the magnitude approximately scaling with the density. The exact mechanism leading to this in- medium structure function modification has not been identified. I will discuss recent results from JLab exploring the EMC effect by using an electron beam to probe two mirror nuclei, Helium-3 and Tritium. iv

Bane, Jason↗
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Magnetic Pumping as a Source of Particle Heating and Power-Law Distributions in the Solar Wind

Based on the rate of expansion of the solar wind, the plasma should cool rapidly as a function of distance to the Sun. Observations show this is not the case. In this work, a magnetic pumping model is developed as a possible explanation for the heating and the generation of power-law distribution functions observed in the solar wind plasma. Most previous studies in this area focus on the role that the dissipation of turbulent energy on microscopic kinetic scales plays in the overall heating of the plasma. However, with magnetic pumping, particles are energized by the largest-scale turbulent fluctuations, thus bypassing the energy cascade. In contrast to other models, we include the pressure anisotropy term, providing a channel for the large-scale fluctuations to heat the plasma directly. A complete set of coupled differential equations describing the evolution, and energization, of the distribution function are derived, as well as an approximate closed-form solution. Numerical simulations using the VPIC kinetic code are applied to verify the model's analytical predictions. The results of the model for realistic solar wind scenario are computed, where thermal streaming of particles are important for generating a phase shift between the magnetic perturbations and the pressure anisotropy. In turn, averaged over a pump cycle, the phase shift permits mechanical work to be converted directly to heat in the plasma. Here, the results of this scenario show that magnetic pumping may account for a significant portion of the solar wind energization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multimodal photon and neutron imaging devices with fast event-by-event multisampling readout and optimized particle type discrimination

A particle imaging method for distinguishing between types of incident particles, such as neutrons, photons, and alphas, and improving the position resolution of particle imaging devices with matrix readout. The method includes high frequency multisampling readout electronics that provides the sequences of multiple measurements for each detected event, resulting in recorded detailed waveform information describing the signals. Such detailed information is used to approximate each signal waveform with a parameterized function in which the extracted parameter sets determine the type of the incident particle in an optimized fashion. The detailed event-by-event multisampling information for each signal readout channel in the matrix readout of the radiation imaging devices improves and optimizes the position resolution for variable shapes of the signals. Such devices can be used in mixed radiation fields, creating a new class of multimodal photon and neutron imagers.

Degtiarenko, Pavel V.↗
DOE OSTI2024

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗