Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A semi-analytic model of gas-puff liner-on-target magneto-inertial fusion

A semi-analytic model is presented for the gas-puff Staged Z-pinch, a magneto-inertial fusion concept in which an annular gas-puff liner implodes onto a deuterium or deuterium-tritium target. The one-dimensional model is a modification of the semi-analytic model for MagLIF (SAMM) [R. D. McBride and S. A. Slutz, Phys. Plasmas 22, 052708 (2015)], that addresses the different set of physics inherent to a Staged Z-pinch implosion: azimuthal magnetic field transport, shock heating of the fuel, separate ion and electron energy equations, and a simplified radiation model that approximates the liner transition from optically thin to optically thick. Following the explanation of the model, three sample problems are presented: first, a Staged Z-pinch implosion on the Zebra driver (1 MA, 100 ns) is modeled and compared with the HYDRA simulation results; second, the MagLIF point design is modeled and compared to the original simulation results [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] and results from SAMM; and third, we conduct a simple parameter scan and scaling study for a Staged Z-pinch implosion on the LTD-III driver (0.8 MA, 160 ns). Here, some agreement with HYDRA and SAMM is obtained, and deuterium-deuterium (DD) neutron yield scaling with current is consistent with other existing models and HYDRA simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An analytic expression approximating the Debye heat capacity function

It is important to have analytic expressions for critical functions in the equations of state of materials. The Debye model has been very successful in approximating the thermal energy properties of a variety of solids, but is nonanalytic. Existing approximations suffer from various shortcomings, the most common being lack of applicability over some temperature range. A new analytic and integrable functional form that closely approximates the Debye model for the heat capacity is presented. This form, based on the mean of two Einstein heat capacity functions with a low temperature correction, exhibits deviations from the Debye model smaller than typical experimental scatter in heat capacity data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. I. Interaction via ordinary and anomalous cyclotron resonances

Conditions for net fast ion drive are derived for beam-driven, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes, such as those routinely observed in spherical tokamaks such as NSTX(-U) and MAST. Both co- and counter-propagating CAEs and GAEs are investigated, driven by the ordinary and anomalous Doppler-shifted cyclotron resonance with fast ions. Whereas prior results were restricted to vanishingly narrow distributions in velocity space, broad parameter regimes are identified in this work which enable an analytic treatment for realistic fast ion distributions generated by neutral beam injection. The simple, approximate conditions derived in these regimes for beam distributions of realistic width compare well to the numerical evaluation of the full analytic expressions for fast ion drive. Additionally, previous results in the very narrow beam case are corrected and generalized to retain all terms in ω/ω ci and |k ∥ /k ⊥ |, which are often assumed to be small parameters but can significantly modify the conditions of drive and damping when they are non-negligible. Lastly, favorable agreement is demonstrated between the approximate stability criterion, simulation results, and a large database of NSTX observations of cntr-GAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparison of modular analytical wake models to the Lillgrund wind plant

Efficient and accurate wake models are required for wind plant performance modeling and the suite of engineering processes that support wind plant layout, control, and monitoring. Although many analytical and engineering wake models with low computational costs have been proposed, their ability to represent the power production of large wind plants in a wide range of atmospheric conditions is not completely understood. The following validation study reviews the underlying theory for analytical wake models, outlines quality control procedures for observational data, and compares model results with observational data from the Lillgrund Wind Plant. Lillgrund makes a valuable case study for wake modeling because of its regular arrangement and the relatively close spacing of constituent wind turbines, which lead to regular and significant wake interactions within the wind plant and the development of deep array flow conditions. Formulations for the velocity deficit, wake-added turbulence, and wake superposition methods are considered in a modular sense, yielding many possible configurations to represent wind turbine wakes, of which seven are examined in detail. Velocity deficit models that account for flow conditions in the near wake are better able to reproduce power production for wind turbines in the transitional region of the wind plant, where wind turbines experience as many as five wakes from upstream turbines. Finally, in the deep array, where power production reaches asymptotic values and wake statistics become quasi-periodic, wake superposition schemes become the largest driver in error reduction; using the linear or maximum wake superposition methods can reduce the relative root mean square error by as much as 40% in the deep array.

17 WIND ENERGY↗
DOE OSTI2020

Analytic investigation of combined effects of anisotropic thermal transport and energetic particles on stability of resistive plasma resistive wall mode

The combined effects of anisotropic thermal transport and trapped energetic particles (EPs) on the stability of the resistive plasma resistive wall mode (RPRWM) are investigated by an energy-principle based analytical model. The results qualitatively confirm that of a recent toroidal modeling study [Bai et al., Phys. Plasmas 27, 072502 (2020)], in which the thermal transport can stabilize the RPPWM depending on the parameters of both the plasma equilibrium and energetic particles. The analytical model predicts a complete stabilization of the RPRWM in highly resistive plasmas, at sufficiently high EPs’ pressure, and finite plasma flow. The stabilizing effect of thermal transport originates from its enhancement of energy dissipations associated with both the resistive layer and the trapped energetic particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic evaluation of energy first derivatives for spin–orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties

The paper reports a formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin–orbit coupling included at the orbital level and an implementation for evaluation of first-order properties. Here, the new development extends the applicability of relativistic coupled cluster methods for heavy-element chemistry and spectroscopy. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the usefulness of the present analytic scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A semi-analytical model for laser induced densification in fused silica at elevated temperatures

A semi-analytical model for deterministic surface smoothing of fused silica below the evaporation temperature threshold is presented. Using this model, a method for calculating deterministic track depths for laser polishing fused silica is derived and validated with experiments. In this work, model predictions were in good agreement with experiment across varying laser preheating temperature and power. With the semi-analytical method presented here, deterministic laser smoothing can be achieved in the densification regime to close the process loop of laser polishing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Analytical model of diffuse reflectance spectrum of skin tissue

We have derived simple analytical expressions that enable highly accurate calculation of diffusely reflected light signals of skin in the spectral range from 450 to 800 nm at a distance from the region of delivery of exciting radiation. The expressions, taking into account the dependence of the detected signals on the refractive index, transport scattering coefficient, absorption coefficient and anisotropy factor of the medium, have been obtained in the approximation of a two-layer medium model (epidermis and dermis) for the same parameters of light scattering but different absorption coefficients of layers. Numerical experiments on the retrieval of the skin biophysical parameters from the diffuse reflectance spectra simulated by the Monte Carlo method show that commercially available fibre-optic spectrophotometers with a fixed distance between the radiation source and detector can reliably determine the concentration of bilirubin, oxy- and deoxyhaemoglobin in the dermis tissues and the tissue structure parameter characterising the size of its effective scatterers. We present the examples of quantitative analysis of the experimental data, confirming the correctness of estimates of biophysical parameters of skin using the obtained analytical expressions. (biophotonics)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Wavefunction stability analysis without analytical electronic Hessians: application to orbital-optimised second-order Møller–Plesset theory and VV10-containing density functionals

Wavefunction stability analysis is commonly applied to converged self-consistent field (SCF) solutions to verify whether the electronic energy is a local minimum with respect to second-order variations in the orbitals. By iterative diagonalisation, the procedure calculates the lowest eigenvalue of the stability matrix or electronic Hessian. However, analytical expressions for the electronic Hessian are unavailable for most advanced post-Hartree-Fock (HF) wave function methods and even some Kohn-Sham (KS) density functionals. To address such cases, here we formulate the Hessian-vector product within the iterative diagonalisation procedure as a finite difference of the electronic gradient with respect to orbital perturbations in the direction of the vector. As a model application, following the lowest eigenvalue of the orbital-optimised second-order Møller-Plesset perturbation theory (OOMP2) Hessian during H 2 dissociation reveals the surprising stability of the spin-restricted solution at all separations, with a second independent unrestricted solution. We show that a single stable solution can be recovered by using the regularised OOMP2 method (δ-OOMP2), which contains a level shift. Internal and external stability analyses are also performed for SCF solutions of a recently developed range-separated hybrid density functional, ωB97X-V, for which the analytical Hessian is not yet available due to the complexity of its long-range non-local VV10 correlation functional.

Brueckner orbitals↗
DOE OSTI2019

Data analytics approach for melt-pool geometries in metal additive manufacturing

Modern data analytics was employed to understand and predict physics-based melt-pool formation by fabricating Ni alloy single tracks using powder bed fusion. An extensive database of melt-pool geometries was created, including processing parameters and material characteristics as input features. Correlation analysis provided insight for relationships between process parameters and melt-pools, and enabled the development of meaningful machine learning models via the use of highly correlated features. We successfully demonstrated that data analytics facilitates understanding of the inherent physics and reliable prediction of melt-pool geometries. This approach can serve as a basis for the melt-pool control and process optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2015

PROBING THE ROLE OF DYNAMICAL FRICTION IN SHAPING THE BSS RADIAL DISTRIBUTION. I. SEMI-ANALYTICAL MODELS AND PRELIMINARY N-BODY SIMULATIONS

We present semi-analytical models and simplified N-body simulations with 10{sup 4} particles aimed at probing the role of dynamical friction (DF) in determining the radial distribution of blue straggler stars (BSSs) in globular clusters. The semi-analytical models show that DF (which is the only evolutionary mechanism at work) is responsible for the formation of a bimodal distribution with a dip progressively moving toward the external regions of the cluster. However, these models fail to reproduce the formation of the long-lived central peak observed in all dynamically evolved clusters. The results of N-body simulations confirm the formation of a sharp central peak, which remains as a stable feature over time regardless of the initial concentration of the system. In spite of noisy behavior, a bimodal distribution forms in many cases, with the size of the dip increasing as a function of time. In the most advanced stages, the distribution becomes monotonic. These results are in agreement with the observations. Also, the shape of the peak and the location of the minimum (which, in most of cases, is within 10 core radii) turn out to be consistent with observational results. For a more detailed and close comparison with observations, including a proper calibration of the timescales of the dynamical processes driving the evolution of the BSS spatial distribution, more realistic simulations will be necessary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ANALYTIC FORMS OF THE PERPENDICULAR DIFFUSION COEFFICIENT IN NRMHD TURBULENCE

In the past different analytic limits for the perpendicular diffusion coefficient of energetic particles interacting with magnetic turbulence were discussed. These different limits or cases correspond to different transport modes describing how the particles are diffusing across the large-scale magnetic field. In the current paper we describe a new transport regime by considering the model of noisy reduced magnetohydrodynamic turbulence. We derive different analytic forms of the perpendicular diffusion coefficient, and while we do this, we focus on the aforementioned new transport mode. We show that for this turbulence model a small perpendicular diffusion coefficient can be obtained so that the latter diffusion coefficient is more than hundred times smaller than the parallel diffusion coefficient. This result is relevant to explain observations in the solar system where such small perpendicular diffusion coefficients have been reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CALIBRATION OF SEMI-ANALYTIC MODELS OF GALAXY FORMATION USING PARTICLE SWARM OPTIMIZATION

We present a fast and accurate method to select an optimal set of parameters in semi-analytic models of galaxy formation and evolution (SAMs). Our approach compares the results of a model against a set of observables applying a stochastic technique called Particle Swarm Optimization (PSO), a self-learning algorithm for localizing regions of maximum likelihood in multidimensional spaces that outperforms traditional sampling methods in terms of computational cost. We apply the PSO technique to the SAG semi-analytic model combined with merger trees extracted from a standard Lambda Cold Dark Matter N-body simulation. The calibration is performed using a combination of observed galaxy properties as constraints, including the local stellar mass function and the black hole to bulge mass relation. We test the ability of the PSO algorithm to find the best set of free parameters of the model by comparing the results with those obtained using a MCMC exploration. Both methods find the same maximum likelihood region, however, the PSO method requires one order of magnitude fewer evaluations. This new approach allows a fast estimation of the best-fitting parameter set in multidimensional spaces, providing a practical tool to test the consequences of including other astrophysical processes in SAMs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

To cool is to accrete: Analytic scalings for nebular accretion of planetary atmospheres

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic expression for how the atmospheric mass grows with time t as a function of the underlying core mass M{sub core} and nebular conditions, including the gas metallicity Z. Planets accrete as much gas as can cool: an atmosphere's doubling time is given by its Kelvin–Helmholtz time. Dusty atmospheres behave differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass ratio (GCR) of a dusty atmosphere scales as GCR ∝ t{sup 0.4}M{sub core}{sup 1.7}Z{sup −0.4}μ{sub rcb}{sup 3.4}, where μ{sub rcb}∝1/(1−Z) (for Z not too close to 1) is the mean molecular weight at the innermost radiative–convective boundary. This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their radiative–convective boundaries, which regulate cooling, are not set by the external environment, but rather by the internal microphysics of dust sublimation, H{sub 2} dissociation, and the formation of H{sup −}. By contrast, dust-free atmospheres have their radiative boundaries at temperatures T{sub rcb} close to nebular temperatures T{sub out}, and grow faster at larger orbital distances where cooler temperatures, and by extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR ∝ t{sup 0.4}T{sub rcb}{sup −1.9}M{sub core}{sup 1.6}Z{sup −0.4}μ{sub rcb}{sup 3.3}, while beyond 1 AU in a gas-rich nebula, GCR ∝ t{sup 0.4}T{sub rcb}{sup −1.5}M{sub core}{sup 1}Z{sup −0.4}μ{sub rcb}{sup 2.2}. We confirm our analytic scalings against detailed numerical models for objects ranging in mass from Mars (0.1M{sub ⊕}) to the most extreme super-Earths (10–20M{sub ⊕}), and explain why heating from planetesimal accretion cannot prevent the latter from undergoing runaway gas accretion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Analytical collisionless damping rate of geodesic acoustic mode

Collisionless damping of geodesic acoustic mode (GAM) is analytically investigated by considering the finite-orbit-width (FOW) resonance effect to the 3rd order in the gyro-kinetic equations. Here, a concise and transparent expression for the damping rate is presented for the first time. Good agreement is found between the analytical damping rate and the previous TEMPEST simulation result (Xu 2008 et al Phys. Rev. Lett. 100 215001) for systematic q scans. Our result also shows that it is of sufficient accuracy and has to take into account the FOW effect to the 3rd order.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Convenient analytical solution for vibrational distribution function of molecules colliding with a wall

We study formation of the Vibrational Distribution Function (VDF) in a molecular gas at low pressure, when vibrational levels are excited by electron impact and deactivated in collisions with walls and show that this problem has a convenient analytical solution that can be used to obtain VDF and its dependence on external parameters. The VDF is determined by excitation of vibrational levels by an external source and deactivation in collisions with the wall. Deactivation in wall collisions is little known process. However, we found that the VDF is weakly dependent on the functional form of the actual form of probability ${\gamma }_{\upsilon ^{\prime} \to \upsilon }$ for a vibrational number $\upsilon ^{\prime} $ to transfer into a lower level $\upsilon $ at the wall. Because for a given excitation source of vibrational states, the problem is linear the solution for VDF involves solving linear matrix equation. The matrix equation can be easily solved if we approximate probability, in the form: ${\gamma }_{\upsilon ^{\prime} \to \upsilon }=(1/\upsilon ^{\prime} )\theta (\upsilon ^{\prime} -\upsilon ).$ In this case, the steady-state solution for VDF($\upsilon $) simply involves a sum of source rates for levels above $\upsilon ,$ with a factor of $1/(\upsilon +1).$ As an example of application, we study the vibrational kinetics in a hydrogen gas and verify the analytical solution by comparing with a full model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗