Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE SPACE”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Large space telescope, phase A. Volume 1: Executive summary

The Phase A study of the Large Space Telescope (LST) is reported. The study defines an LST concept based on the broad mission guidelines provided by the Office of Space Science (OSS), the scientific requirements developed by OSS with the scientific community, and an understanding of long range NASA planning current at the time the study was performed. The LST is an unmanned astronomical observatory facility, consisting of an optical telescope assembly (OTA), scientific instrument package (SIP), and a support systems module (SSM). The report consists of five volumes. The report describes the constraints and trade off analyses that were performed to arrive at a reference design for each system and for the overall LST configuration. A low cost design approach was followed in the Phase A study. This resulted in the use of standard spacecraft hardware, the provision for maintenance at the black box level, growth potential in systems designs, and the sharing of shuttle maintenance flights with other payloads.

Source record↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

Interaction of ring current ions with the magnetopause

An analysis of phase space density observations of 30-130 keV ring current protons made during a magnetopause crossing of ISEE 1 indicates that the phase space densities of 30-65 keV protons within the magnetosphere are higher for particles streaming parallel to the field than antiparallel to the field. In the magnetosheath, phase space densities parallel to the field are higher than the densities of both outward- and inward-flowing particles parallel to the field within the magnetosphere. The observations are discussed in the terms of reflected and transmitted ring current protons' acceleration in the magneteopause current layer.

Scholer, M.↗
NASA NTRS1988

Plasma electrons as tracers of distant magnetotail structure - ISEE-3

This paper compares the electron spectra and phase space densities measured concurrently by ISEE-3 at 200 R(E), with those measured by DMSP at low altitudes. The field-aligned lobe electron phase space densities above 200 eV at ISEE were found to agree well with the DMSP-measured polar rain phase space densities near the polar cap; the spectral slopes above 200 eV were also similar. Below 100-200 eV, a thermal electron population was measured by ISEE in the distant tail, which arose from local entry of plasma through the distant magnetopause, which is not present at DMSP altitudes. These data show that the suprathermal tail lobe electrons are essentially a test particle population which can move freely along field lines to form polar rain; in contrast, the thermal electrons are bound to the tailward-flowing lobe ion population far down the tail and, thus, cannot reach the polar cap regions.

Baker, D. N.↗
DOE OSTI2025

Neural refinement of sample weights

Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be negative, and several post hoc methods have been developed to eliminate or mitigate the negative values. All of these methods share the common strategy of approximating the average weight as a function of phase space. We introduce an alternative approach, which, instead of reweighting to the average, refines the initial weights with a scaling transformation, utilizing a phase space-dependent factor. Since this new refinement method does not need to model the full weight distribution, it can be more accurate. High-dimensional and unbinned phase space is processed using neural networks for the refinement method. In addition to the refinement method, we introduce a new resampling protocol, which can be used in conjunction with any weight transformation to not only preserve the average weight but also the statistical uncertainties of the initial distribution. Using both realistic and synthetic examples, we show that the new neural refinement method is able to match or exceed the accuracy of similar weight transformations and that the new resampling protocol is simpler in implementation than previous methods while exhibiting equivalent statistical properties.

Artificial neural networks↗
NASA NTRS2002

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗
NASA NTRS2003

The Statistical Mechanics of Ideal MHD Turbulence

Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Shebalin, John V.↗
NASA NTRS1995

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗
NASA NTRS1993

Electron distributions in the inner Jovian magnetosphere: Voyager 1 observations

Using several improvements in the analysis of the observations of the Low Energy Charged Particle (LECP) experiment on Voyager 1, electron phase space densities in the inner Jovian magnetosphere (5 - 10 R(sub J) were first calculated at constant first and second invariants (represented by mu and K, respectively), based on the LECP measurements. The calculated electron phase space density profiles show that in the inner Jovian magnetosphere there exist evident time and longitude variations, energetic electron injections, and present radial transport and distributed losses. To study the radial and pitch angle diffusions of Jovian electrons, we have calculated the phase space densities in the K-L space. It is found that the electron population in the inner Jovian magnetosphere seems to consist of two components: electrons radially diffusing from a main external source and electrons generated from local sources. The radially diffusing electrons have a relatively time stationary and isotropic distribution, while the locally created electrons mainly concentrate around the equatorial plane and have relatively lower energies, in comparison with the inward diffusing electrons. Consequently, the sources of precipitation losses to the ionosphere must be primarily electrons transported from outer sources, and the major precipitations occur in the inner magnetosphere (L less than 7.5 R(sub J). In the inner Jovian magnetosphere (L = 5 to approximately 10 R(sub J)) it is estimated that for electrons with magnetic moment mu = 300 MeV/G, the diffusion coefficient D is roughly 10(exp -8) to approximately 10(exp -6) R(exp 2)(sub J)/s, and the lifetime against the diffusion losses is of the order of 10(exp 4) to approximately 10(exp 6) s.

Ye, G.↗
NASA NTRS1994

Radial diffusion models of energetic electrons and Jupiter's synchrotron radiation. 2: Time variability

We used a radial diffusion code for energetic electrons in Jupiter's magnetosphere to investigate variations in Jupiter's radio emission due to changes in the electron phase space density at L shells between 6 and 50, and due to changes in the radial diffusion parameters. We suggest that the observed variations in Jupiter's radio emission are likely caused by changes in the electron phase space density at some boundary L(sub 1) is greater than 6, if the primary mode of transport of energetic electrons is radial diffusion driven by fluctuating electric and/or magnetic fields induced by upper atmospheric turbulence. We noticed an excellent empirical correlation, both in phase and relative amplitude, between changes in the solar wind ram pressure and Jupiter's synchrotron radiation if the electron phase space density at the boundary L(sub 1) (L(sub 1) is approximately equal to 20-50) varies linearly with the square root of the solar wind ram pressure, f is approximately (N(sub s)nu(exp 2 sub s))(exp 1/2). The calculations were carried out with a diffusion coefficient D(sub LL) = D(sub n)L(exp n) with n = 3. The diffusion coefficient which best fit the observed variations in Jupiter's synchrotron radiation D(sub 3) = 1.3 +/- 0.2 x 10(exp -9)/s is approximately 0.041/yr, which corresponds to a lagtime of approximately 2 years. We further show that the observed short term (days-weeks) variations in Jupiter's radio emission cannot be explained adequately when radial diffusion is taken into account.

De Pater, Imke↗
NASA NTRS1996

Affordable In-Space Transportation Phase 2: An Advanced Concepts Project

The Affordable In-Space Transportation (AIST) program was established by the NASA Office of Space Access to improve transportation and lower the costs from Low Earth Orbit (LEO) to Geostationary Earth Orbit (GEO) and beyond (to Lunar orbit, Mars orbit, inner solar system missions, and return to LEO). A goal was established to identify and develop radically innovative concepts for new upper stages for Reusable Launch Vehicles (RLV) and Highly Reusable Space Transportation (HRST) systems. New architectures and technologies are being identified which have the potential to meet a cost goal of $1,000 to $2,000 per pound for transportation to GEO and beyond for overall mission cost (including the cost to LEO). A Technical Interchange Meeting (TTM) was held on October 16 and 17, 1996 in Huntsville, Alabama to review previous studies, present advanced concepts and review technologies that could be used to meet the stated goals. The TIN4 was managed by NASA-Marshall Space Flight Center (MSFC) Advanced Concepts Office with Mr. Alan Adams providing TIM coordination. Mr. John C. Mankins of NASA Headquarters provided overall sponsorship. The University of Alabama in Huntsville (UAH) Propulsion Research Center hosted the TIM at the UAH Research Center. Dr. Clark Hawk, Center Director, was the principal investigator. Technical support was provided by Christensen Associates. Approximately 70 attendees were present at the meeting. This Executive Summary provides a record of the key discussions and results of the TIN4 in a summary for-mat. It incorporates the response to the following basic issues of the TDVL which addressed the following questions: 1. What are the cost drivers and how can they be reduced? 2. What are the operational issues and their impact on cost? 3. What is the current technology readiness level (TRL) and what will it take to reach TRL 6? 4. What are the key enabling technologies and sequence for their accomplishment? 5 . What is the proposed implementation time frame? See Appendix A for the TIM Agenda and Appendix C for the AIST Program Terms of Reference.

Source record↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
NASA NTRS2012

Vlasov Simulation of Electrostatic Solitary Structures in Multi-Component Plasmas

Electrostatic solitary structures have been observed in the Earth's magnetosheath by the Cluster spacecraft. Recent theoretical work has suggested that these solitary structures are modeled by electron acoustic solitary waves existing in a four-component plasma system consisting of core electrons, two counter-streaming electron beams, and one species of background ions. In this paper, the excitation of electron acoustic waves and the formation of solitary structures are studied by means of a one-dimensional electrostatic Vlasov simulation. The present result first shows that either electron acoustic solitary waves with negative potential or electron phase-space holes with positive potential are excited in four-component plasma systems. However, these electrostatic solitary structures have longer duration times and higher wave amplitudes than the solitary structures observed in the magnetosheath. The result indicates that a high-speed and small free energy source may be needed as a fifth component. An additional simulation of a five-component plasma consisting of a stable four-component plasma and a weak electron beam shows the generation of small and fast electron phase-space holes by the bump-on-tail instability. The physical properties of the small and fast electron phase-space holes are very similar to those obtained by the previous theoretical analysis. The amplitude and duration time of solitary structures in the simulation are also in agreement with the Cluster observation.

magnetohydrodynamics (MHD)↗
NASA NTRS1992

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗
NASA NTRS1971

Modular space station phase B extension: NASA headquarters quarterly review

Progress in the modular space station development program is reported. The space station project office organization and the manner in which the center interfaces with NR are presented along with a sortie and experiment analysis. The operations, subsystems and configurations are reviewed, and future activities are listed.

Cole, E. G.↗
NASA NTRS1989

Radial diffusion and losses of energetic protons in the 5 to 12 R(S) region of Saturn's magnetosphere

Phase space density profiles for nearly equatorial mirroring energetic protons, derived from Pioneer 11 measurements during the 1979 Saturn encounter, were analyzed using solutions of the time-averaged radial diffusion equation in a dipolar planetary magnetic field. Then, several loss models were considered, ranging from a model that takes into account satellite absorption only to the one that accounts for satellite plus ring-E absorption plus added distributed losses, and the corresponding form of the time-averaged radial diffusion coefficient was determined for each case by a minimum-variance fit to the data-derived phase space densiy profiles. The model comparisons clearly indicated the need for distributed losses, in addition to satellite absorption, in order to explain the smooth decline of the experimental energetic proton phase space density profiles in the 5-12-Saturn radius region of Saturn's magnetosphere.

Hood, L. L.↗
DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗