Search NASA⌕ Search

SEARCH · Search NASA

Results for “London equation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Quantitative Assessment of Fatigue Damage Accumulation in Wavy Slip Metals from Acoustic Harmonic Generation

A comprehensive, analytical treatment is presented of the microelastic-plastic nonlinearities resulting from the interaction of a stress perturbation with dislocation substructures (veins and persistent slip bands) and cracks that evolve during high-cycle fatigue of wavy slip metals. The nonlinear interaction is quantified by a material (acoustic) nonlinearity parameter beta extracted from acoustic harmonic generation measurements. The contribution to beta from the substructures is obtained from the analysis of Cantrell [Cantrell, J. H., 2004, Proc. R. Soc. London A, 460, 757]. The contribution to beta from cracks is obtained by applying the Paris law for crack propagation to the Nazarov-Sutin crack nonlinearity equation [Nazarov, V. E., and Sutin, A. M., 1997, J. Acoust. Soc. Am. 102, 3349]. The nonlinearity parameter resulting from the two contributions is predicted to increase monotonically by hundreds of percent during fatigue from the virgin state to fracture. The increase in beta during the first 80-90 percent of fatigue life is dominated by the evolution of dislocation substructures, while the last 10-20 percent is dominated by crack growth. The model is applied to the fatigue of aluminium alloy 2024-T4 in stress-controlled loading at 276MPa for which experimental data are reported. The agreement between theory and experiment is excellent.

Cantrell, John H.↗
DOE OSTI2020

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Partition function corrections

The vibration-rotation interaction coupling in diatomic molecules is analyzed and first-order expressions for the energy levels are deduced. Corrections to the partition function to account for this coupling are derived. The effects of long-range intermolecular forces due to induced dipole-induced dipole interactions, dipole-induced dipole interactions, charge-dipole interactions, and charge-charge interactions are considered. London's quantized oscillator model of the induced dipole-induced dipole interaction is developed. The partition function including effects of such dispersion forces is derived and is shown to lead to the virial equation of state. For a model combining the hard sphere and long-range attractive type potentials, this is found to lead to the van der Waals equation of state. The equation of state near the critical point is presented in reduced form (thermodynamic variables expressed in units of critical point values). Finally, approximate corrections to the electronic partition functions are derived to account for perturbations of excited electronic states induced by neighboring neutral particles, heavy ions, and fast electron collisions.

Source record↗
DOE OSTI2025

Two-Scale Neural Networks for Partial Differential Equations with Small Parameters

We propose a two-scale neural network method for solving partial differential equations (PDEs) with small parameters using physics-informed neural networks (PINNs). We directly incorporate the small parameters into the architecture of neural networks. The proposed method enables solving PDEs with small parameters in a simple fashion, without adding Fourier features or other computationally taxing searches of truncation parameters. Various numerical examples demonstrate reasonable accuracy in capturing features of large derivatives in the solutions caused by small parameters.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Numerical simulations of the seasonal/latitudinal variations of atomic oxygen and nitric oxide in the lower thermosphere and mesosphere

A 2-Dimensional zonally-averaged thermospheric model and the global University College London (UCL) thermospheric model have been used to investigate the seasonal, solar activity and geomagnetic variation of atomic oxygen and nitric oxide. The 2-dimensional model includes detailed oxygen and nitrogen chemistry, with appropriate completion of the energy equation, by adding the thermal infrared cooling by O and NO. This solution includes solar and auroral production of odd nitrogen compounds and metastable species. This model has been used for three investigations; firstly, to study the interactions between atmospheric dynamics and minor species transport and density; secondly, to examine the seasonal variations of atomic oxygen and nitric oxide within the upper mesosphere and thermosphere and their response to solar and geomagnetic activity variations; thirdly, to study the factor of 7 to 8 peak nitric oxide density increase as solar F sub 10.7 cm flux increases from 70 to 240 reported from the Solar Mesospheric Explorer. Auroral production of NO is shown to be the dominant source at high latitudes, generating peak NO densities a factor of 10 greater than typical number densities at low latitudes. At low latitudes, the predicted variation of the peak NO density, near 110 km, with the solar F sub 10.7 cm flux is rather smaller than is observed. This is most likely due to an overestimate of the soft X-ray flux at low solar activity, for times of extremely low support number, as occurred in June 1986. As observed on pressure levels, the variation of O density is small. The global circulation during solstice and periods of elevated geomagnetic activity causes depletion of O in regions of upwelling, and enhancements in regions of downwelling.

Rees, D.↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Atomic Data and Spectral Line Intensities for NI XVII

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for Ni XVII. We include in the calculations the 23 lowest configurations, corresponding to 159 fine-structure levels: 3l3l', 3l4l0'' , and 3s5l0''' , with l,l' = s,p,d, l'' = s,p,d, f, and l''' = s,p,d. Collision strengths are calculated at five incident energies for all transitions at varying energies above the threshold of each transition. One additional energy, very close to the threshold of each transition, has also been included. Calculations have been carried out using the Flexible Atomic Code in the distorted wave approximation. Additional calculations have been performed with the University College London suite of codes for comparison. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates of the present work, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8) - 10(exp 14) / cubic cm and at an electron temperature of logT(sub e)e(K) = 6.5, corresponding to the maximum abundance of Ni XVII. Spectral line intensities are calculated, and their diagnostic relevance is discussed. This dataset will be made available in the next version of the CHIANTI database

Bhatia, A. K.↗
NASA NTRS1983

User's manual for a TEACH computer program for the analysis of turbulent, swirling reacting flow in a research combustor

Described is a computer program for the analysis of the subsonic, swirling, reacting turbulent flow in an axisymmetric, bluff-body research combustor. The program features an improved finite-difference procedure designed to reduce the effects of numerical diffusion and a new algorithm for predicting the pressure distribution within the combustor. A research version of the computer program described in the report was supplied to United Technologies Research Center by Professor A. D. Gosman and his students, R. Benodeker and R. I. Issa, of Imperial College, London. The Imperial College staff also supplied much of the program documentation. Presented are a description of the mathematical model for flow within an axisymmetric bluff-body combustor, the development of the finite-difference procedure used to represent the system of equations, an outline of the algorithm for determining the static pressure distribution within the combustor, a description of the computer program including its input format, and the results for representative test cases.

Chiappetta, L. M.↗
DOE OSTI2022

Two-dimensional parametric design of short pulse laser driven iron opacity experiments

High intensity short pulse lasers are commonly used to create high temperature and high density plasmas. These plasmas are used to study fundamental properties such as the equation of state and opacity. The targets often use small-diameter, thin buried layers of the materials of interest, surrounded by a tamper of low atomic number material such as plastic. Computational modeling is used to design and interpret experiments with short pulse lasers. Most of the modeling to date has assumed one-dimensional plane-parallel geometry. In this paper, the effects of radial gradients in the irradiation of thin planar targets are studied with one- and two-dimensional radiation/hydrodynamic simulations. It is found that the main effect of radial gradients is the averaging of plasma conditions and x-ray emission over the pattern of irradiation. Additionally, differences between one- and two-dimensional simulations arise because the plasma conditions and x-ray emission are, in general, nonlinear functions of the temperature of the plasma and thereby also nonlinear functions of the irradiation intensity. The differences increase with the ratio of the buried layer radius to the laser spot radius. The root mean square difference in the inferred iron L-shell opacity is less than 30% when the ratio is less than 1.0 and about 70% when the ratio is 1.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

An EPR methodology for measuring the London penetration depth for the ceramic superconductors

The use is discussed of electron paramagnetic resonance (EPR) as a quick and easily accessible method for measuring the London penetration depth, lambda for the high T(sub c) superconductors. The method utilizes the broadening of the EPR signal, due to the emergence of the magnetic flux lattice, of a free radical adsorbed on the surface of the sample. The second moment, of the EPR signal below T(sub c) is fitted to the Brandt equation for a simple triangular lattice. The precision of this method compares quite favorably with those of the more standard methods such as micro sup(+)SR, Neutron scattering, and magnetic susceptibility.

Rakvin, B.↗
NASA NTRS1991

Microwave (EPR) measurements of the penetration depth measurements of high-Tc superconductors

The use is discussed of electron paramagnetic resonance (EPR) as a quick and easily accessible method for measuring the London penetration depth, lambda for the high T sub c superconductors. The method uses the broadening of the EPR signal, due to the emergence of the magnetic flux lattice, of a free radical adsorbed on the surface of the sample. The second moment, of the EPR signal below T sub c is fitted to the Brandt equation for a simple triangular lattice. The precision of this method compares quite favorably with those of the more standard methods such as micro sup(+)SR, neutron scattering, and magnetic susceptibility.

Dalal, N. S.↗
DOE OSTI2025

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗
DOE OSTI2025

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗
NASA NTRS2005

Momentum-Based Dynamics for Spacecraft with Chained Revolute Appendages

An efficient formulation is presented for a sub-class of multi-body dynamics problems that involve a six degree-of-freedom base body and a chain of N rigid linkages connected in series by single degree-of-freedom revolute joints. This general method is particularly well suited for simulations of spacecraft dynamics and control that include the modeling of an orbiting platform with or without internal degrees of freedom such as reaction wheels, dampers, and/or booms. In the present work, particular emphasis is placed on dynamic simulation of multi-linkage robotic manipulators. The differential equations of motion are explicitly given in terms of linear and angular momentum states, which can be evaluated recursively along a serial chain of linkages for an efficient real-time solution on par with the best of the O(N3) methods.

Queen, Steven↗
DOE OSTI2025

A platform to measure isentropes from proton-heated warm dense matter on short pulse laser facilities

We describe the development of an experimental platform that measures the release isentrope of materials heated isochorically to temperatures of a few electron volts, using short-pulse laser-produced protons to heat the sample and long-pulse laser-produced x rays to perform streaked x-ray radiography. The density profiles derived from the radiography data are integrated to generate pressure–density isentropes, independent of prior knowledge of the equation of state of the sample material. In order to understand the sensitivities of isentrope extraction from radiography data, we analyze synthetic radiographs generated by a radiation hydrodynamics code. Noise reduction and high spatial resolution are critical for isentrope reconstruction, as demonstrated by the analysis of a proof-of-principle shot day on the OMEGA-EP facility. In conclusion, the data demonstrate the feasibility of the platform for characterizing isentropes, and we discuss the necessary improvements to enhance precision in differentiating between equation-of-state models.

Equations of state↗
NASA NTRS1992

Einstein coefficients for rotational lines of the (0,0) band of the NO A2sigma(+)-X2Pi system

A summary of the spectroscopic equations necessary for prediction of the molecular transition energies and the Einstein A and B coefficients for rovibronic lines of the gamma(0,0) band of nitric oxide (NO) is presented. The calculated molecular transition energies are all within 0.57/cm of published experimental values; in addition, over 95 percent of the calculated energies give agreement with measured results within 0.25/cm. Einstein coefficients are calculated from the band A00 value and the known Hoenl-London factors and are tabulated for individual rovibronic transitions in the NO A2sigma(+)-X2Pi(0,0) band.

Reisel, John R.↗
DOE OSTI2026

Compression behavior of lithium fluoride up to 80 GPa and 2300 K

Here, the high-temperature compression behavior of lithium fluoride (LiF) has been determined to ∼80 GPa and 2300 K by means of high-resolution synchrotron-based x-ray diffraction in a laser-heated diamond-anvil cell. The room-temperature Vinet equation of state (EOS) for LiF yields an isothermal room-temperature, ambient-pressure lattice parameter a 0 = 4.028 (±0.003) Å, volume V 0 = 65.35 (±0.03) Å 3 , bulk modulus K 0 = 67.57 (±0.34) GPa, and its pressure derivative K 0 ’ = 4.64 (±0.03). This expanded experimental dataset is in good agreement with the LiF thermal equation of state by Myint et al. (2019) up to ∼2300 K and 25 GPa. At higher pressures and temperatures our data shows reduced thermal expansion compared to the prediction and may be particularly important for experiments at more extreme conditions.

High-pressure behavior↗
NASA NTRS2010

GSFC DORIS Contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd 10 submissions included DORIS data from Envisat, TOPEX!Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT -2, SPOT -3, SPOT -4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 x 10-9 m/s(exp 2) with the a priori macromodel to 0.99 x 10-9 m/s2 with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8-0.9 cm in the radial direction, 2.1-3.4 cm cross-track, and 1.7-2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of +/- 5-1O micro s. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993-1994 to 1.5 cm in 2007-2008.

Le Bail, K.↗