Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximations and Expansions”

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 19 records

Asymptotic expansions and approximations for the Caputo derivative

In this paper, we use the asymptotic expansions of the binomial coefficients and the weights of the L1 approximation to obtain approximations of order 2−α and second-order approximations of the Caputo derivative by modifying the weights of the shifted Grünwald–Letnikov difference approximation and the L1 approximation of the Caputo derivative. A modification of the shifted Grünwald–Letnikov approximation is obtained which allows second-order numerical solutions of fractional differential equations with arbitrary values of the solutions and their first derivatives at the initial point.

97 MATHEMATICS AND COMPUTING↗

Heat of Mixing During Fast Charge/Discharge of a Li-Ion Cell: A Study on NMC523 Cathode

Predicting temperature rise accurately during fast charge/discharge of a Li-ion cell is essential to avoid thermal runaway and extend battery life. While modeling the temperature rise, it is necessary to model the heat generation correctly. Heat of mixing, one of the four main sources of heat generation in a Li-ion battery, has often been considered insignificant and therefore excluded from modeling. When included, it is modeled using the expression from a Taylor expansion approximation. In this work, we have shown the conditions when including the heat of mixing becomes important and quantified the error associated with using the Taylor expansion, especially under high charge/discharge conditions. Consequently, we carry out the calculation of the rate of enthalpy change and the heat generation rate from the most fundamental equation for the rate of the total enthalpy change without simplifying assumptions or approximations. The heat generation rate calculated doing so naturally includes irreversible, reversible and mixing heat. We then exclusively separate out the heat of mixing by subtracting the rate of enthalpy change by reaction from the rate of the total enthalpy change. Results show that the contribution of heat of mixing in the total heat generated increases with the charge/discharge rate and is as large as 23% for a 6 C discharge. This result suggests that while modeling heat generation for fast charge/discharge, it is necessary to include the heat of mixing and avoid calculating it using the Taylor expansion approximation.

25 ENERGY STORAGE↗

Improved Venus ionopause altitude calculation and comparison with measurement

Improved approximations are incorporated into the inviscid fluid method originated by Spreiter et al (1970) in calculations of the altitude of the Venus ionopause. The altitude calculations are then compared with median altitudes measured by the Pioneer Venus retarding potential analyzer. The calculated ionopause shape is found to closely approximate the measured meridian shape in the solar zenith angle (SZA) range of zero-135 deg. The use of improved ionospheric pressure field ionosheath pressure and Prandtl-Meyer expansion approximations lowers the terminator ionopause altitude to approximately half that obtained with the usual Spreiter approximations. It is also determined that the calculated dawn ionopause is about 300 km higher than the dusk ionopause, and that both ionopauses are close to their respective measured meridian values. It is concluded that median ionopause altitude within the SZA angle range may be calculated without inclusion of a viscous interaction in the theory.

Knudsen, W. C.↗

Effects of various assumptions on the calculated liquid fraction in isentropic saturated equilibrium expansions

The saturated equilibrium expansion approximation for two phase flow often involves ideal-gas and latent-heat assumptions to simplify the solution procedure. This approach is well documented by Wegener and Mack and works best at low pressures where deviations from ideal-gas behavior are small. A thermodynamic expression for liquid mass fraction that is decoupled from the equations of fluid mechanics is used to compare the effects of the various assumptions on nitrogen-gas saturated equilibrium expansion flow starting at 8.81 atm, 2.99 atm, and 0.45 atm, which are conditions representative of transonic cryogenic wind tunnels. For the highest pressure case, the entire set of ideal-gas and latent-heat assumptions are shown to be in error by 62 percent for the values of heat capacity and latent heat. An approximation of the exact, real-gas expression is also developed using a constant, two phase isentropic expansion coefficient which results in an error of only 2 percent for the high pressure case.

Bursik, J. W.↗

Scaling Symmetries in Elastic-Plastic Dynamic Cavity Expansion Equations Using the Isovector Method

Cavity-expansion approximations are widely-used in the study of penetration mechanics and indentation phenomena. We apply the isovector method to a well-established model in the literature for elastic-plastic cavity-expansion to systematically demonstrate the existence of Lie symmetries corresponding to scale-invariant solutions. Here we use the symmetries obtained from the equations of motion to determine compatible auxiliary conditions describing the cavity wall trajectory and the elastic-plastic material interface. The admissible conditions are then compared with specific similarity solutions in the literature.

97 MATHEMATICS AND COMPUTING↗

Efficiency of Uncertainty Propagation Methods for Estimating Output Moments

Uncertainty propagation methods are used to estimate the distribution of model outputs resulting from a set of uncertain model outputs. There are a number of uncertainty propagation methods available in literature. This paper compares six non-intrusive uncertainty propagation methods, Latin Hypercube Sampling, Full Factorial Integration, Univariate Dimension Reduction, Halton series, Sobol series, and Polynomial Chaos Expansion, in terms of their efficiency for estimating the first four moments of the output distribution using computational experiments. Here, the results suggest employing FFNI if there are few uncertain inputs, up to three. Uncertainty propagation methods that utilize Halton and Sobol series are found to be robust for estimating output moments as the number of uncertain inputs increases. In general, higher order polynomial chaos expansion approximations (3rd-5th order) obtained accurate estimates of model outputs with fewer model evaluations.

97 MATHEMATICS AND COMPUTING↗

Ionic Asymmetry and Solvent Excluded Volume Effects on Spherical Electric Double Layers: A Density Functional Approach

In this article we present a classical density functional theory for electrical double layers of spherical macroions that extends the capabilities of conventional approaches by accounting for electrostatic ion correlations, size asymmetry and excluded volume effects. The approach is based on a recent approximation introduced by Hansen-Goos and Roth for the hard sphere excess free energy of inhomogeneous fluids (J. Chem. Phys. 124, 154506). It accounts for the proper and efficient description of the effects of ionic asymmetry and solvent excluded volume, especially at high ion concentrations and size asymmetry ratios including those observed in experimental studies. Additionally, we utilize a leading functional Taylor expansion approximation of the ion density profiles. In addition, we use the Mean Spherical Approximation for multi-component charged hard sphere fluids to account for the electrostatic ion correlation effects. These approximations are implemented in our theoretical formulation into a suitable decomposition of the excess free energy which plays a key role in capturing the complex interplay between charge correlations and excluded volume effects. We perform Monte Carlo simulations in various scenarios to validate the proposed approach, obtaining a good compromise between accuracy and computational cost. We use the proposed computational approach to study the effects of ion size, ion size asymmetry and solvent excluded volume on the ion profiles, integrated charge, mean electrostatic potential, and ionic coordination number around spherical macroions in various electrolyte mixtures. Our results show that both solvent hard sphere diameter and density play a dominant role in the distribution of ions around spherical macroions, mainly for experimental water molarity and size values where the counterion distribution is characterized by a tight binding to the macroion, similar to that predicted by the Stern model.

Medasani, Bharat↗

New Alloy for Glass-to-Metal Seals

Coefficient of thermal expansion approximates that of glass more closely. Alloy composed of about 60 percent iron, 40 percent nickel, and traces of six other elements. Developed as replacement for Kovar Fe/Ni/Co alloy in ferrule-and-tube assembly, new alloy has same strength, solderability, and compatibility with fuel as does Kovar. Used in glass-to-metal seals without excessive residual stresses. Potential for other applications in which low thermal expansion important; mechanical measuring devices and precise sliding parts that must function over wide temperature ranges.

Schmuck, A. J.↗

Development of Pore-Free Ti-Si-C MAX/Al-Si Composite Materials Manufactured by Squeeze Casting Infiltration

Open-porous MAX phase skeletons from Ti{sub 3}SiC{sub 2} were manufactured by Microwave-Assisted Self-propagating High-temperature Synthesis (MASHS) and subsequently subjected to squeeze casting infiltration with an Al-Si lightweight casting alloy (EN AC-44200). This alloy was chosen due to its high flowability, corrosion resistance and good machinability. The manufactured composites, together with a reference sample of the original alloy, underwent testing of thermal properties, including thermal conductivity and diffusivity, specific heat and thermal expansion in the temperature range 50-500 °C, which corresponds to the expected working temperatures of the material. The fabricated AlSi/Ti{sub 3}SiC{sub 2} composites have significantly increased thermal stability, with coefficients of thermal expansion (approximately 10-11 × 10{sup −6} °C{sup −1}) half that of the original alloy. As regards mechanical properties, the instrumental Young’s modulus and Vickers hardness of the composite materials are 170.8 and 8.5 GPa, respectively. Moreover, the microstructure and phase composition, structural defects and potential impacts between constituents of the manufactured composites were characterized using SEM, TEM and STEM microscopy and EDS and XRD analysis.

36 MATERIALS SCIENCE↗

Frozen flux violation, electron demagnetization and magnetic reconnection

We argue that the analogue in collisionless plasma of the collisional diffusion region of magnetic reconnection is properly defined in terms of the demagnetization of the plasma electrons that enable “frozen flux” slippage to occur. This condition differs from the violation of the “frozen-in” condition, which only implies that two fluid effects are involved, rather than the necessary slippage of magnetic flux as viewed in the electron frame. Using 2D Particle In Cell (PIC) simulations, this approach properly finds the saddle point region of the flux function. Our demagnetization conditions are the dimensionless guiding center approximation expansion parameters for electrons which we show are observable and determined locally by the ratio of non-ideal electric to magnetic field strengths. Proxies for frozen flux slippage are developed that (a) are measurable on a single spacecraft, (b) are dimensionless with theoretically justified threshold values of significance, and (c) are shown in 2D simulations to recover distinctions theoretically possible with the (unmeasurable) flux function. A new potentially observable dimensionless frozen flux rate, Λ{sub Φ}, differentiates significant from anecdotal frozen flux slippage. A single spacecraft observable, ϒ, is shown with PIC simulations to be essentially proportional to the unobservable local Maxwell frozen flux rate. This relationship theoretically establishes electron demagnetization in 3D as the general cause of frozen flux slippage. In simple 2D cases with an isolated central diffusion region surrounded by separatrices, these diagnostics uniquely identify the traditional diffusion region (without confusing it with the two fluid “ion-diffusion” region) and clarify the role of the separatrices where frozen flux violations do occur but are not substantial. In the more complicated guide and asymmetric 2D cases, substantial flux slippage regions extend out along, but inside of, the preferred separatrices, demonstrating that Λ{sub Φ} ≠ 0 violations are present over significant distances (in ion inertial units) from the separator identified by the 2D flux function; these violations are, however, generally weaker than seen at known separators in 2D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A New Method for Interpreting Nonstationary Running Correlations and Its Application to the ENSO‐EAWM Relationship

Abstract We here propose a new statistical method to interpret nonstationary running correlations by decomposing them into a stationary part and a first‐order Taylor expansion approximation for the nonstationary part. Then, this method is applied to investigate the nonstationary behavior of the El Niño–Southern Oscillation (ENSO)‐East Asian winter monsoon (EAWM) relationship, which exhibits prominent multidecadal variations. It is demonstrated that the first‐order approximation of the nonstationary part can be expressed to a large extent by the impact of the nonlinear interaction between the Atlantic Multidecadal Oscillation (AMO) and ENSO (AMO*Niño3.4) on the EAWM. Therefore, the nonstationarity in the ENSO‐EAWM relationship comes predominantly from the impact of an AMO modulation on the ENSO‐EAWM teleconnection via this key nonlinear interaction. This general method can be applied to investigate nonstationary relationships that are often observed between various different climate phenomena.

Geng, Xin↗

Towards an accurate description of perovskite ferroelectrics: exchange and correlation effects

Using the van der Waals density functional with C09 exchange (vdW-DF-C09), which has been applied to describing a wide range of dispersion-bound systems, we explore the physical properties of prototypical ABO 3 bulk ferroelectric oxides. Surprisingly, vdW-DF-C09 provides a superior description of experimental values for lattice constants, polarization and bulk moduli, exhibiting similar accuracy to the modified Perdew-Burke-Erzenhoff functional which was designed specifically for bulk solids (PBEsol). The relative performance of vdW-DF-C09 is strongly linked to the form of the exchange enhancement factor which, like PBEsol, tends to behave like the gradient expansion approximation for small reduced gradients. These results suggest the general-purpose nature of the class of vdW-DF functionals, with particular consequences for predicting material functionality across dense and sparse matter regimes.

36 MATERIALS SCIENCE↗

Neutral atomic jet generation by laser ablation of copper targets

This work aimed the obtainment of a neutral atomic jet departing from a plume generated by laser ablation of copper targets. A pair of electrodes together with a transducer pressure sensor was used to study the ablated plume charge composition and also to measure the ion extraction from the plasma plume. The neutral beam was produced with this setup and the relative abundance of neutrals in the plasma was measured, it decreases from 30% to 8% when the laser fluence is varied from 20 J/cm{sup 2} to 32 J/cm{sup 2}. The necessary voltage to completely remove the ions from the plume varied from 10 V to 230 V in the same fluence range. TOF analysis resulted in center of mass velocities between 3.4 and 4.6 km/s, longitudinal temperature in the range from 1 × 10{sup 4} K to 2.4 × 10{sup 4} K and a Mach number of M = 2.36, calculated using purely hydrodynamic expansion approximation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Density functional theory of material design: fundamentals and applications—II

Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.

36 MATERIALS SCIENCE↗

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗

A mathematical analysis of the theory of interplanetary scintillation in the weak scattering approximation

A simplified analytical technique is presented for modeling the interplanetary scintillation of radio sources of finite angular size with a power-law electron-density-fluctuation power spectrum. The simplification results from the representation of the scintillation spectrum in confluent hypergeometric functions. The approximations presented allow fast numerical evaluation of a spectrum for a weakly scattering but extended medium with less than 10% error over the entire spectrum. Parameters describing anisotropic electron irregularities as well as anisotropic source structure are included, and the dependence of the spectrum normalization on the scales of the medium is derived explicitly. The parametric description of the domains of convergence of the approximate expansions also provides a simple conceptualization of the relative contributions of the scattered radiation along the line of sight to the observed spectrum. This is particularly useful for sources of finite angular size. This technique is applied to previously published observations.

Mitchell, D. G.↗