SEARCH · Search NASA

Results for “generalized spherical functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗

NASA NTRS2006

The Variance of Solar Wind Magnetic Fluctuations: Solutions and Further Puzzles

We study the dependence of the variance directions of the magnetic field in the solar wind as a function of scale, radial distance, and Alfvenicity. The study resolves the question of why different studies have arrived at widely differing values for the maximum to minimum power (approximately equal to 3:1 up to approximately equal to 20:1). This is due to the decreasing anisotropy with increasing time interval chosen for the variance, and is a direct result of the "spherical polarization" of the waves which follows from the near constancy of |B|. The reason for the magnitude preserving evolution is still unresolved. Moreover, while the long-known tendency for the minimum variance to lie along the mean field also follows from this view (as shown by Barnes many years ago), there is no theory for why the minimum variance follows the field direction as the Parker angle changes. We show that this turning is quite generally true in Alfvenic regions over a wide range of heliocentric distances. The fact that nonAlfvenic regions, while still showing strong power anisotropies, tend to have a much broader range of angles between the minimum variance and the mean field makes it unlikely that the cause of the variance turning is to be found in a turbulence mechanism. There are no obvious alternative mechanisms, leaving us with another intriguing puzzle.

Roberts, D. A.↗

NASA NTRS1995

Upper Thermosphere Winds and Temperatures in the Geomagnetic Polar Cap: Solar Cycle, Geomagnetic Activity, and Interplanetary Magnetic Field Dependencies

Ground-based Fabry-Perot interferometers located at Thule, Greenland (76.5 deg. N, 69.0 deg. W, lambda = 86 deg.) and at Sondre Stromfjord, Greenland (67.0 deg. N, 50.9 deg. W, lambda = 74 deg.) have monitored the upper thermospheric (approx. 240-km altitude) neutral wind and temperature over the northern hemisphere geomagnetic polar cap since 1983 and 1985, respectively. The thermospheric observations are obtained by determining the Doppler characteristics of the (OI) 15,867-K (630.0-nm) emission of atomic oxygen. The instruments operate on a routine, automatic, (mostly) untended basis during the winter observing seasons, with data coverage limited only by cloud cover and (occasional) instrument failures. This unique database of geomagnetic polar cap measurements now extends over the complete range of solar activity. We present an analysis of the measurements made between 1985 (near solar minimum) and 1991 (near solar maximum), as part of a long-term study of geomagnetic polar cap thermospheric climatology. The measurements from a total of 902 nights of observations are compared with the predictions of two semiempirical models: the Vector Spherical Harmonic (VSH) model of Killeen et al. (1987) and the Horizontal Wind Model (HWM) of Hedin et al. (1991). The results are also analyzed using calculations of thermospheric momentum forcing terms from the Thermosphere-ionosphere General Circulation Model TGCM) of the National Center for Atmospheric Research (NCAR). The experimental results show that upper thermospheric winds in the geomagnetic polar cap have a fundamental diurnal character, with typical wind speeds of about 200 m/s at solar minimum, rising to up to about 800 m/s at solar maximum, depending on geomagnetic activity level. These winds generally blow in the antisunward direction, but are interrupted by episodes of modified wind velocity and altered direction often associated with changes in the orientation of the Interplanetary Magnetic Field (IMF). The central polar cap (greater than approx. 80 magnetic latitude) antisunward wind speed is found to be a strong function of both solar and geomagnetic activity. The polar cap temperatures show variations in both solar and geomagnetic activity, with temperatures near 800 K for low K(sub p) and F(sub 10.7) and greater than about 2000 K for high K(sub p) and F(sub 10.7). The observed temperatures are significantly greater than those predicted by the mass spectrometer/incoherent scatter model for high activity conditions. Theoretical analysis based on the NCAR TIGCM indicates that the antisunward upper thermospheric winds, driven by upstream ion drag, basically 'coast' across the polar cap. The relatively small changes in wind velocity and direction within the polar cap are induced by a combination of forcing terms of commensurate magnitude, including the nonlinear advection term, the Coriolis term, and the pressure gradient force term. The polar cap thennospheric thermal balance is dominated by horizontal advection, and adiabatic and thermal conduction terms.

Killeen, T. L.↗

NASA NTRS2010

Detecting Thin Cirrus in Multiangle Imaging Spectroradiometer Aerosol Retrievals

Thin cirrus clouds (optical depth (OD) < 03) are often undetected by standard cloud masking in satellite aerosol retrieval algorithms. However, the Mu]tiangle Imaging Spectroradiometer (MISR) aerosol retrieval has the potential to discriminate between the scattering phase functions of cirrus and aerosols, thus separating these components. Theoretical tests show that MISR is sensitive to cirrus OD within Max{0.05 1 20%l, similar to MISR's sensitivity to aerosol OD, and MISR can distinguish between small and large crystals, even at low latitudes, where the range of scattering angles observed by MISR is smallest. Including just two cirrus components in the aerosol retrieval algorithm would capture typical MISR sensitivity to the natural range of cinus properties; in situations where cirrus is present but the retrieval comparison space lacks these components, the retrieval tends to underestimate OD. Generally, MISR can also distinguish between cirrus and common aerosol types when the proper cirrus and aerosol optical models are included in the retrieval comparison space and total column OD is >-0.2. However, in some cases, especially at low latitudes, cirrus can be mistaken for some combinations of dust and large nonabsorbing spherical aerosols, raising a caution about retrievals in dusty marine regions when cirrus is present. Comparisons of MISR with lidar and Aerosol Robotic Network show good agreement in a majority of the cases, but situations where cirrus clouds have optical depths >0.15 and are horizontally inhomogeneous on spatial scales shorter than 50 km pose difficulties for cirrus retrieval using the MISR standard aerosol algorithm..

Pierce, Jeffrey R.↗

DOE OSTI2025

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Simultaneous Measurements of Noncommuting Observables: Positive Transformations and Instrumental Lie Groups

We formulate a general program for describing and analyzing continuous, differential weak, simultaneous measurements of noncommuting observables, which focuses on describing the measuring instrument autonomously, without states. The Kraus operators of such measuring processes are time-ordered products of fundamental differential positive transformations, which generate nonunitary transformation groups that we call instrumental Lie groups. The temporal evolution of the instrument is equivalent to the diffusion of a Kraus-operator distribution function, defined relative to the invariant measure of the instrumental Lie group. This diffusion can be analyzed using Wiener path integration, stochastic differential equations, or a Fokker-Planck-Kolmogorov equation. This way of considering instrument evolution we call the Instrument Manifold Program. We relate the Instrument Manifold Program to state-based stochastic master equations. We then explain how the Instrument Manifold Program can be used to describe instrument evolution in terms of a universal cover that we call the universal instrumental Lie group, which is independent not just of states, but also of Hilbert space. The universal instrument is generically infinite dimensional, in which case the instrument’s evolution is chaotic. Special simultaneous measurements have a finite-dimensional universal instrument, in which case the instrument is considered principal, and it can be analyzed within the differential geometry of the universal instrumental Lie group. Principal instruments belong at the foundation of quantum mechanics. We consider the three most fundamental examples: measurement of a single observable, position and momentum, and the three components of angular momentum. As these measurements are performed continuously, they converge to strong simultaneous measurements. For a single observable, this results in the standard decay of coherence between inequivalent irreducible representations. For the latter two cases, it leads to a collapse within each irreducible representation onto the classical or spherical phase space, with the phase space located at the boundary of these instrumental Lie groups.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2016

Optics of Water Microdroplets with Soot Inclusions: Exact Versus Approximate Results

We use the recently generalized version of the multi-sphere superposition T-matrix method (STMM) to compute the scattering and absorption properties of microscopic water droplets contaminated by black carbon. The soot material is assumed to be randomly distributed throughout the droplet interior in the form of numerous small spherical inclusions. Our numerically-exact STMM results are compared with approximate ones obtained using the Maxwell-Garnett effective-medium approximation (MGA) and the Monte Carlo ray-tracing approximation (MCRTA). We show that the popular MGA can be used to calculate the droplet optical cross sections, single-scattering albedo, and asymmetry parameter provided that the soot inclusions are quasi-uniformly distributed throughout the droplet interior, but can fail in computations of the elements of the scattering matrix depending on the volume fraction of soot inclusions. The integral radiative characteristics computed with the MCRTA can deviate more significantly from their exact STMM counterparts, while accurate MCRTA computations of the phase function require droplet size parameters substantially exceeding 60.

cloud droplets↗

DOE OSTI2022

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

A numerical modeling study of the interaction between the tides and the circulation forced by high-latitude plasma convection

A spectral, time-varying thermospheric general circulation model has been used to study the nonlinear interaction at high latitudes between the tides propagating into the thermosphere from below and the circulation induced by magnetospheric forcing and in situ solar heating. The model is discrete in the vertical with 27 layers spaced by half a scale height. In the horizontal, the fields are expanded in a series of spherical harmonics using a triangular truncation at wave number 31, equivalent to a homogeneous global resolution with a minimum wavelength of 1270 km. A hypothetical uniform grid point model would require a horizontal spacing of 417 km to describe the same minimum wavelength. In the high-latitude F region the tides affect the dusk vortex of the neutral flow very little, but the dawn vortex is either suppressed or amplified dependent upon the universal time and tidal phase. In the E region neutral flow, both the dusk and dawn vortices are shifted in local time by the tides, again as a function of universal time and tidal phase. At dusk a nonlinear amplification of the sunward winds occurs for certain combinations of parameters, and at dawn the winds may be completely suppressed. Below 120 km altitude the magnetospheric forcing creates a single cyclonic vortex which is also sensitive to the high-latitude tidal structure.

Mikkelsen, I. S.↗

DOE OSTI2021

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1997

Topography of the Moon from the Clementine Lidar

Range measurements from the lidar instrument carried aboard the Clementine spacecraft have been used to produce an accurate global topographic model of the Moon. This paper discusses the function of the lidar; the acquisition, processing, and filtering of observations to produce a global topographic model; and the determination of parameters that define the fundamental shape of the Moon. Our topographic model: a 72nd degree and order spherical harmonic expansion of lunar radii, is designated Goddard Lunar Topography Model 2 (GLTM 2). This topographic field has an absolute vertical accuracy of approximately 100 m and a spatial resolution of 2.5 deg. The field shows that the Moon can be described as a sphere with maximum positive and negative deviations of approx. 8 km, both occurring on the farside, in the areas of the Korolev and South Pole-Aitken (S.P.-Aitken) basins. The amplitude spectrum of the topography shows more power at longer wavelengths as compared to previous models, owing to more complete sampling of the surface, particularly the farside. A comparison of elevations derived from the Clementine lidar to control point elevations from the Apollo laser altimeters indicates that measured relative topographic heights generally agree to within approx. 200 in over the maria. While the major axis of the lunar gravity field is aligned in the Earth-Moon direction, the major axis of topography is displaced from this line by approximately 10 deg to the cast and intersects the farside 24 deg north of the equator. The magnitude of impact basin topography is greater than the lunar flattening (approx. 2 km) and equatorial ellipticity (approx. 800 m), which imposes a significant challenge to interpreting the lunar figure. The floors of mare basins are shown to lie close to an equipotential surface, while the floors of unflooded large basins, except for S.P.-Aitken, lie above this equipotential. The radii of basin floors are thus consistent with a hydrostatic mechanism for the absence of significant farside maria except for S.P.-Aitken, whose depth and lack of mare require significant internal compositional and/or thermal heterogeneity. A macroscale surface roughness map shows that roughness at length scales of 10(exp 1) - 10(exp 2) km correlates with elevation and surface age.

Smith, David E.↗

DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗

DOE OSTI2025

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗

NASA NTRS1998

The Use of Langmuir Probes in Non-Maxwellian Space Plasmas

Disturbance of the Maxwellian plasma may occur in the vicinity of a spacecraft due to photoemission, interactions between the spacecraft and thermospheric gases, or electron emissions from other devices on the spacecraft. Significant non-maxwellian plasma distributions may also occur in nature as a mixture of ionospheric and magnetospheric plasmas or secondaries produced by photoionization in the thermosphere or auroral precipitation. The general formulas for current collection (volt-ampere curves) by planar, cylindrical, and spherical Langmuir probes in isotropic and anisotropic non-maxwellian plasmas are examined. Examples are given of how one may identify and remove the non-maxwellian components in the Langmuir probe current to permit the ionospheric parameters to be determined. Theoretical volt-ampere curves presented for typical examples of non-maxwellian distributions include: two-temperature plasmas and a thermal plasma with an energetic electron beam. If the non-ionospheric electrons are Maxwellian at a temperature distinct from that of the ionosphere electrons, the volt-ampere curves can be fitted directly to obtain the temperatures and densities of both electron components without resorting to differenting the current. For an arbitrary isotropic distribution, the current for retarded particles is shown to be identical for the three geometries. For anisotropic distributions, the three probe geometries are not equally suited for measuring the ionospheric electron temperature and density or for determining the distribution function in the presence of non-maxwellian back-round electrons.

Hoegy, Walter R.↗

NASA NTRS2003

Kinetics of Nucleation and Crystal Growth in Glass Forming Melts in Microgravity

This flight definition project has the specific objective of investigating the kinetics of nucleation and crystal growth in high temperature inorganic oxide, glass forming melts in microgravity. It is related to one1 of our previous NASA projects that was concerned with glass formation for high temperature containerless melts in microgravity. The previous work culminated in two experiments which were conducted aboard the space shuttle in 1983 and 1985 and which consisted of melting (at 1500 C) and cooling levitated 6 to 8 mm diameter spherical samples in a Single Axis Acoustic Levitator (SAAL) furnace. Compared to other types of materials, there have been relatively few experiments, 6 to 8, conducted on inorganic glasses in space. These experiments have been concerned with mass transport (alkali diffusion), containerless melting, critical cooling rate for glass formation, chemical homogeneity, fiber pulling, and crystallization of glass forming melts. One of the most important and consistent findings in all of these experiments has been that the glasses prepared in microgravity are more resistant to crystallization (better glass former) and more chemically homogeneous than equivalent glasses made on earth (1g). The chemical composition of the melt appears relatively unimportant since the same general results have been reported for oxide, fluoride and chalcogenide melts. These results for space-processed glasses have important implications, since glasses with a higher resistance to crystallization or higher chemical homogeneity than those attainable on earth can significantly advance applications in areas such as fiber optics communications, high power laser glasses, and other photonic devices where glasses are the key functional materials. The classical theories for nucleation and crystal growth for a glass or melt do not contain any parameter that is directly dependent upon the g-value, so it is not readily apparent why glasses prepared in microgravity should be more resistant to crystallization than equivalent glasses prepared on earth. Similarly, the gravity-driven convection in a fluid melt is believed to be the primary force field that is responsible for melt homogenization on earth. Thus, it is not obvious why a glass prepared in space, where gravity-driven convection is ideally absent, would be more chemically homogeneous than a glass identically prepared on earth. The primary objective of the present research is to obtain experimental data for the nucleation rate and crystal growth rate for a well characterized silicate melt (lithium disilicate) processed entirely in space (low gravity) and compare these rates with the nucleation and crystal growth rates for a similar glass prepared identically on earth (1g).

Day, Delbert E.↗

NASA NTRS2001

Kinetics of Nucleation and Crystal Growth in Glass Forming Melts in Microgravity

Day, Delbert E.↗