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 613 records · Page 34

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

We have developed efficient numerical algorithms for the solution of 3D steady-state Poisson-Nernst-Planck equations (PNP) with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are discretized by finite difference scheme and solved iteratively by Gummel method with relaxation. The Nernst-Planck equations are transformed into Laplace equations through the Slotboom transformation. Algebraic multigrid method is then applied to efficiently solve the Poisson equation and the transformed Nernst-Planck equations. A novel strategy for calculating excess chemical potentials through fast Fourier transforms is proposed which reduces computational complexity from O(N2) to O(NlogN) where N is the number of grid points. Integrals involving Dirac delta function are evaluated directly by coordinate transformation which yields more accurate result compared to applying numerical quadrature to an approximated delta function. Numerical results for ion and electron transport in solid electrolyte for Li ion batteries are shown to be in good agreement with the experimental data and the results from previous studies.

Meng, Da↗

Actuator and aerodynamic modeling for high-angle-of-attack aeroservoelasticity

Accurate prediction of airframe/actuation coupling is required by the imposing demands of modern flight control systems. In particular, for agility enhancement at high angle of attack and low dynamic pressure, structural integration characteristics such as hinge moments, effective actuator stiffness, and airframe/control surface damping can have a significant effect on stability predictions. Actuator responses are customarily represented with low-order transfer functions matched to actuator test data, and control surface stiffness is often modeled as a linear spring. The inclusion of the physical properties of actuation and its installation on the airframe is therefore addressed in this paper using detailed actuator models which consider the physical, electrical, and mechanical elements of actuation. The aeroservoelastic analysis procedure is described in which the actuators are modeled as detailed high-order transfer functions and as approximate low-order transfer functions. The impacts of unsteady aerodynamic modeling on aeroservoelastic stability are also investigated in this paper by varying the order of approximation, or number of aerodynamic lag states, in the analysis. Test data from a thrust-vectoring configuration of an F/A-18 aircraft are compared to predictions to determine the effects on accuracy as a function of modeling complexity.

Brenner, Martin J.↗

Scaling laws and stability of nano-sized defect clusters in niobium via atomistic simulations and statistical analysis

The predictions of scaling laws for the structure and properties of defect clusters are generally limited to small defect clusters in their ground-state configurations. We investigated the size and geometrical configuration dependence of nano-sized defect clusters in niobium (Nb) using molecular dynamics. We studied the structure and stability of large clusters of size up to fifty defects for vacancies and one hundred defects for interstitials, as well as the role of helium and metastable configurations on the stability of these clusters. We compared three different interatomic potentials in order to determine the relative stability of these clusters as a function of their size and geometrical configurations. Additionally, we conducted a statistical analysis to predict the formation and binding energies of interstitial clusters as a function of both their size and configuration. We find that the size dependence of vacancy and interstitial clusters can be approximated by functional forms that account for bulk and surface effects as well as some considerations of elastic interactions. We also find that helium and metastable configurations can make vacancy and interstitial clusters thermally stable depending on the configuration. Our parameterized functional forms for the formation and binding energies are valid for a very broad range of defect sizes and configurations making it possible to be used directly in a coarse-grained modeling strategy such as Monte Carlo, cluster dynamics or dislocation dynamics which look at defect accumulation and evolution in microstructures.

36 MATERIALS SCIENCE↗

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Soft mean spherical approximation for dusty plasma liquids: Level of accuracy and analytic expressions

The soft mean spherical approximation is employed for the study of the thermodynamics of dusty plasma liquids, the latter treated as Yukawa one-component plasmas. Within this integral theory method, the only input necessary for the calculation of the reduced excess energy stems from the solution of a single non-linear algebraic equation. Consequently, thermodynamic quantities can be routinely computed without the need to determine the pair correlation function or the structure factor. The level of accuracy of the approach is quantified after an extensive comparison with numerical simulation results. The approach is solved over a million times with input spanning the whole parameter space and reliable analytic expressions are obtained for the basic thermodynamic quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fragment-based {sup 13}C nuclear magnetic resonance chemical shift predictions in molecular crystals: An alternative to planewave methods

We assess the quality of fragment-based ab initio isotropic {sup 13}C chemical shift predictions for a collection of 25 molecular crystals with eight different density functionals. We explore the relative performance of cluster, two-body fragment, combined cluster/fragment, and the planewave gauge-including projector augmented wave (GIPAW) models relative to experiment. When electrostatic embedding is employed to capture many-body polarization effects, the simple and computationally inexpensive two-body fragment model predicts both isotropic {sup 13}C chemical shifts and the chemical shielding tensors as well as both cluster models and the GIPAW approach. Unlike the GIPAW approach, hybrid density functionals can be used readily in a fragment model, and all four hybrid functionals tested here (PBE0, B3LYP, B3PW91, and B97-2) predict chemical shifts in noticeably better agreement with experiment than the four generalized gradient approximation (GGA) functionals considered (PBE, OPBE, BLYP, and BP86). A set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided based on these benchmark calculations. Statistical cross-validation procedures are used to demonstrate the robustness of these fits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamical systems and nonlinear transient rheology of the far-from-equilibrium Bjorken flow

In relativistic kinetic theory, the one-particle distribution function is approximated by an asymptotic perturbative power series in the Knudsen number which is divergent. For the Bjorken flow, we expand the distribution function in terms of its moments and study their nonlinear evolution equations. The resulting coupled dynamical system can be solved for each moment consistently using a multiparameter transseries which makes the constitutive relations inherit the same structure. A new nonperturbative dynamical renormalization scheme is born out of this formalism that goes beyond the linear response theory.We show that there is a Lyapunov function, also known as dynamical potential, which is, in general, a function of the moments and time satisfying Lyapunov stability conditions along renormalization group flows connected to the asymptotic hydrodynamic fixed point. As a result, the transport coefficients get dynamically renormalized at every order in the time-dependent perturbative expansion by receiving nonperturbative corrections present in the transseries. The connection between the integration constants and the UV data is discussed using the language of dynamical systems. Furthermore, we show that the first dissipative correction in the Knudsen number to the distribution function is not only determined by the known effective shear viscous term but also a new high-energy nonhydrodynamic mode. It is demonstrated that the survival of this new mode is intrinsically related to the nonlinear mode-to-mode coupling with the shear viscous term. Finally, we comment on some possible phenomenological applications of the proposed nonhydrodynamic transport theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Atmospheric refraction errors in laser ranging systems

The effects of horizontal refractivity gradients on the accuracy of laser ranging systems were investigated by ray tracing through three dimensional refractivity profiles. The profiles were generated by performing a multiple regression on measurements from seven or eight radiosondes, using a refractivity model which provided for both linear and quadratic variations in the horizontal direction. The range correction due to horizontal gradients was found to be an approximately sinusoidal function of azimuth having a minimum near 0 deg azimuth and a maximum near 180 deg azimuth. The peak to peak variation was approximately 5 centimeters at 10 deg elevation and decreased to less than 1 millimeter at 80 deg elevation.

Gardner, C. S.↗

Image reconstruction in higher dimensions: myocardial perfusion imaging of tracer dynamics with cardiac motion due to deformation and respiration

Myocardial perfusion imaging (MPI) using slow rotating large field of view cameras requires spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration. In vivo, MPI contains additional degrees of freedom involving unavoidable motion of the heart due to quasiperiodic beating and the effects of respiration, which can severely degrade the quality of the images. This work develops a technique for a single photon emission computed tomography (SPECT) that reconstructs the distribution of the radiotracer concentration in the myocardium using a tensor product of different sets of basis functions that approximately describe the spatiotemporal variation of the radiotracer concentration and the motion of the heart. In this study the temporal B-spline basis functions are chosen to reflect the dynamics of the radiotracer, while the intrinsic deformation and the extrinsic motion of the heart are described by a product of a discrete set of Gaussian basis functions. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it deforms due to cardiac beating, and is displaced due to respiratory motion. We find these results are compared with the conventional 4D-spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. The higher dimensional reconstruction method proposed here improves bias, yet the signal-to-noise ratio (SNR) decreases slightly due to redistribution of the counts over the cardiac-respiratory gates. Finally, there is a trade-off between the number of gates and the number of projections per gate to achieve high contrast images.

60 APPLIED LIFE SCIENCES↗

A fast semi-implicit algorithm for problems of mixed type

Certain physical processes are modeled by partial differential equations which are parabolic over part of the domain and elliptic over the remainder. A family of semi-implicit algorithms which are well suited to initial-boundary value problems of this mixed type is discussed. One important feature of these algorithms is the use of an approximate inverse for the solution of the implicit linear system. A strong error analysis results in an estimate of the total error as a function of approximate inverse error e and time step h.

Frederickson, P. O.↗

Influence of the Variations of Orthotropy Parameters on the Stress State of Hollow Cylinders with Concave Corrugated Cross Sections

In the three-dimensional statement, by the method of approximation of functions by discrete Fourier series, we analyze the influence of orthotropy parameters on the stress state of hollow cylinders with cross sections in the form of connected concave semicorrugations subjected to the action of internal pressure under certain conditions on the ends. The results of investigations of the stress state are presented in the form of plots and a table.

36 MATERIALS SCIENCE↗

Accurate Electron Affinities and Orbital Energies of Anions from a Nonempirically Tuned Range-Separated Density Functional Theory Approach

The treatment of atomic anions with Kohn–Sham density functional theory (DFT) has long been controversial because the highest occupied molecular orbital (HOMO) energy, E HOMO , is often calculated to be positive with most approximate density functionals. We assess the accuracy of orbital energies and electron affinities for all three rows of elements in the periodic table (H–Ar) using a variety of theoretical approaches and customized basis sets. Among all of the theoretical methods studied here, we find that a nonempirically tuned range-separated approach (constructed to satisfy DFT-Koopmans’ theorem for the anionic electron system) provides the best accuracy for a variety of basis sets, even for small basis sets where most functionals typically fail. Previous approaches to solve this conundrum of positive E HOMO values have utilized non-self-consistent methods; however, electronic properties, such as electronic couplings/gradients (which require a self-consistent potential and energy), become ill-defined with these approaches. In contrast, the nonempirically tuned range-separated procedure used here yields well-defined electronic couplings/gradients and correct E HOMO values because both the potential and resulting electronic energy are computed self-consistently. Orbital energies and electron affinities are further analyzed in the context of the electronic energy as a function of electronic number (including fractional numbers of electrons) to provide a stringent assessment of self-interaction errors for these complex anion systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy spectra and wave function of trigonometric Rosen-Morse potential as an effective quantum chromodynamics potential in D-dimensions

The Energy Spectra and Wave Function of Schrodinger equation in D-Dimensions for trigonometric Rosen-Morse potential were investigated analytically using Nikiforov-Uvarov method. This potential captures the essential traits of the quark-gluon dynamics of Quantum Chromodynamics. The approximate energy spectra are given in the close form and the corresponding approximate wave function for arbitrary l-state (l ≠ 0) in D-dimensions are formulated in the form of differential polynomials. The wave function of this potential unnormalizable for general case. The wave function of this potential unnormalizable for general case. The existence of extra dimensions (centrifugal factor) and this potential increase the energy spectra of system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Piecewise linear approximation with minimum number of linear segments and minimum error: A fast approach to tighten and warm start the hierarchical mixed integer formulation

In several areas of economics and engineering, it is often necessary to fit discrete data points or approximate nonlinear functions with continuous functions. Piecewise linear (PWL) functions are a convenient way to achieve this. PWL functions can be modeled in mathematical problems using only linear and integer variables. Moreover, there is a computational benefit in using PWL functions that have the least possible number of segments. This work proposes a novel hierarchical mixed integer linear programming (MILP) formulation that identifies a continuous PWL approximation with minimum number of linear segments for a given target maximum error. The proposed MILP formulation also identifies the solution with the least maximum error among the solutions with minimum number of segments. Then, this work proposes a fast iterative algorithm that identifies non necessarily continuous PWL approximations by solving O(S log N) linear programming (LP) problems, where N is the number of data points and S is the minimum number of segments in the non necessarily continuous case. This work demonstrates that tight bounds for the MILP problem can be derived from these approximations. Next, a fast algorithm is introduced to transform a non necessarily continuous PWL approximation into a continuous one. Finally, the tight bounds and the continuous PWL approximations are used to tighten and warm start the MILP problem. The tightened formulation is shown in experimental results to be more efficient, especially for large data sets, with a solution time that is up to two orders of magnitude less than the existing literature.

97 MATHEMATICS AND COMPUTING↗

Fermi orbital derivatives in self-interaction corrected density functional theory: Applications to closed shell atoms

A recent modification of the Perdew-Zunger self-interaction-correction to the density-functional formalism has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Löwdin orthonormalized Fermi-orbitals which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested, here, on atoms. Total energies and ionization energies in closed-shell singlet atoms, where correlation is less important, using the Perdew-Wang 1992 Local Density Approximation (PW92) functional, are in good agreement with experiment and non-relativistic quantum-Monte-Carlo results albeit slightly too low.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗