Error analysis of hybrid simulation for spacecraft re-entry studies Final report
Error analysis of hybrid simulation for six degree of freedom equations of motion for spacecraft reentry
Jackson, A. S.↗
SEARCH · Search NASA
Results for “hybrid simulation”
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.
Spontaneous Hot Flow Anomalies at Quasi-Parallel Shocks: 2. Hybrid Simulations
Motivated by recent THEMIS observations, this paper uses 2.5-D electromagnetic hybrid simulations to investigate the formation of Spontaneous Hot Flow Anomalies (SHFA) upstream of quasi-parallel bow shocks during steady solar wind conditions and in the absence of discontinuities. The results show the formation of a large number of structures along and upstream of the quasi-parallel bow shock. Their outer edges exhibit density and magnetic field enhancements, while their cores exhibit drops in density, magnetic field, solar wind velocity and enhancements in ion temperature. Using virtual spacecraft in the simulation, we show that the signatures of these structures in the time series data are very similar to those of SHFAs seen in THEMIS data and conclude that they correspond to SHFAs. Examination of the simulation data shows that SHFAs form as the result of foreshock cavitons interacting with the bow shock. Foreshock cavitons in turn form due to the nonlinear evolution of ULF waves generated by the interaction of the solar wind with the backstreaming ions. Because foreshock cavitons are an inherent part of the shock dissipation process, the formation of SHFAs is also an inherent part of the dissipation process leading to a highly non-uniform plasma in the quasi-parallel magnetosheath including large scale density and magnetic field cavities.
Quasi-Parallel↗
Hybrid simulation techniques applied to the earth's bow shock
The application of a hybrid simulation model, in which the ions are treated as discrete particles and the electrons as a massless charge-neutralizing fluid, to the study of the earth's bow shock is discussed. The essentials of the numerical methods are described in detail; movement of the ions, solution of the electromagnetic fields and electron fluid equations, and imposition of appropriate boundary and initial conditions. Examples of results of calculations for perpendicular shocks are presented which demonstrate the need for a kinetic treatment of the ions to reproduce the correct ion dynamics and the corresponding shock structure. Results for oblique shocks are also presented to show how the magnetic field and ion motion differ from the perpendicular case.
Winske, D.↗
Two-dimensional hybrid simulation of a curved bow shock
Results are presented from two-dimensional hybrid simulations of curved collisionless supercritical shocks, retaining both quasi-perpendicular and quasi-parallel sections of the shock in order to study the character and origin of the foreshock ion population. The simulations demonstrate that the foreshock ion population is dominated by ions impinging upon the quasi-parallel side of the shock, while nonlocal transport from the quasi-perpendicular side of the shock into the foreshock region is minimal. Further, it is shown that the ions gain energy by drifting significantly in the direction of the convection electric field through multiple shock encounters.
Thomas, V. A.↗
A hybrid simulation of contact discontinuity
Contact discontinuities in a collisionless plasma are studied by hybrid simulations, in which ions are treated as particles and electrons are considered as a fluid. It is demonstrated that contact discontinuity with a stable density ramp can exist in cases with a finite electron temperature. An electron pressure gradient is present across the contact discontinuity, leading to the presence of a parallel electric field and hence field-aligned potential increase (Delta Phi (sub parallel)) in the transition region. By reflecting ions at the discontinuity, this parallel electric potential peak reduces the interpenetration between hot and cold ions and maintains a stable density ramp across the contact discontinuity. The ratio of the field-aligned electric potential energy to ion thermal energy, e(Delta) Phi(sub parallel)/kT(sub i), is found to be an increasing function of T(sub e)/T(sub i), where T(sub e) and T(sub i) are respectively the electron and ion temperature.
Wu, B. H.↗
Hybrid Simulations of Interstellar Pick-up Ion Acceleration at the Solar Wind Termination Shock
Hybrid simulations have been used to study self consistently the acceleration of interstellar pick-up ions at the solar wind termination shock. Results are presented form one-dimensional simulations of high Mach number oblique shocks with a 10% population of Interstella pick-up hydrogen.
anomalous cosmic rays↗
Using the Apollo hybrid simulator to define astronaut-spacecraft interfaces for extremely low-orbit lunar landmark tracking
Apollo hybrid simulator for man-machine interface in low orbit lunar landmark tracking
Sabionski, G. R.↗
The Interaction of the Solar Wind with Solar Probe Plus - 3D Hybrid Simulation: The Study for the Distance 4.5Rs - Report 1
Our report devotes a 3D numerical hybrid model of the interaction of the solar wind with the Solar Probe spacecraft. The SPP model includes 3 main parts, namely, a non-conducting heat shield, a support system, and cylindrical section or spacecraft bus that contains the particle analysis devices and antenna. One observes an excitation of the low frequency Alfven and whistler type wave directed by the magnetic field with an amplitude of about (0.06-0.6) V/m. The compression waves and the jumps in an electric field with an amplitude of about (0.15-0.7) V/m were also observed. The wave amplitudes are comparable to or greater than previously estimated max wave amplitudes that SPP is expected to measure. The results of our hybrid simulation will be useful for understanding the plasma environment near the SPP spacecraft at the distance 4.5 Rs. Future simulation will take into account the charging of the spacecraft, the charge separation effects, an outgassing from heat shield, a photoionization and an electron impact ionization effects near the spacecraft.
Lipatov, Alexander S.↗
The Interaction of the Solar Wind with Solar Probe Plus - 3D Hybrid Simulation: The Study for the Distance 4.5Rs - Report 1
Our report devotes a 3D numerical hybrid model of the interaction of the solar wind with the Solar Probe spacecraft. The Solar Probe Plus (SPP) model includes 3 main parts, namely, a non-conducting heat shield, a support system, and cylindrical section or spacecraft bus that contains the particle analysis devices and antenna. One observes an excitation of the low frequency Alfven and whistler type wave directed by the magnetic field with an amplitude of about (0.06-0.6) V/m. The compression waves and the jumps in an electric field with an amplitude of about (0.15-0.7) V/m were also observed. The wave amplitudes are comparable to or greater than previously estimated max wave amplitudes that SPP is expected to measure. The results of our hybrid simulation will be useful for understanding the plasma environment near the SPP spacecraft at the distance 4.5 Rs. Future simulation will take into account the charging of the spacecraft, the charge separation effects, an outgassing from heat shield, a photoionization and an electron impact ionization effects near the spacecraft.
Lipatov, Alexander S.↗
Hybrid Simulation of the Shock Wave Trailing the Moon
A standing shock wave behind the Moon was predicted by Michel (1967) but never observed nor simulated. We use 1D hybrid code in order to simulate the collapse of the plasma-free cavity behind the Moon and for the first time to model the formation of this shock. Starting immediately downstream of the obstacle we consider the evolution of plasma expansion into the cavity in the frame of reference moving along with the solar wind. Well-known effects as electric charging of the cavity affecting the plasma flow and counterstreaming ion beams in the wake are reproduced. Near the apex of the inner Mach cone where the plasma flows from the opposite sides of the obstacle meet, a shock wave arises. We expect the shock to be produced at periods of high electron temperature solar wind streams (T(sub i) much less than T(sub e) approximately 100 eV). The shock is produced by the interaction of oppositely directed proton beams in the plane containing solar wind velocity and interplanetary magnetic field vectors. In the direction across the magnetic field and the solar wind velocity, the shock results from the interaction of the plasma flow with the region of the enhanced magnetic field inside the cavity that plays the role of the magnetic barrier. The appearance of the standing shock wave is expected at the distance of approximately 7R(sub M) downstream of the Moon.
Moon↗
Soft X-Ray Imaging of Earth's Dayside Magnetosheath and Cusps Using Hybrid Simulations
Interactions between solar wind ions and neutral hydrogen atoms in Earth's exosphere can lead to the emission of soft X-rays. Upcoming missions such as SMILE and LEXI aim to use soft X-ray imaging to study the global structure of the magnetosphere. Although the magnetosheath and dayside magnetopause can often be driven by kinetic physics, it has typically been omitted from fluid simulations used to predict X-ray emissions. We study the possible results of soft X-ray imaging using hybrid simulations under quasi-radial interplanetary magnetic fields, where ion-ion instabilities drive ultra-low frequency foreshock waves, leading to turbulence in the magnetosheath, affecting the dynamics of the cusp and magnetopause. We simulate soft X-ray emission to determine what may be seen by missions such as LEXI, and evaluate the possibility of identifying kinetic structures. While kinetic structures are visible in high-cadence imaging, current instruments may not have the time resolution to discern kinetic signals.
J Ng↗
A Study of Slow-mode Shocks in the Near-Earth Magnetotail with MMS Observations and Hybrid Simulations
The structure of the magnetic reconnection boundary, particularly the presence of slow-mode shocks in the near-Earth magnetotail was studied by using magnetospheric multiscale (MMS) observations and 2.5D hybrid simulations. A total of 51 crossings of MMS from 2017 to 2021 were analyzed. We found that the detection percentage of slow-mode shocks in the near-Earth magnetotail is 41%–55%. Previous studies have only reported one slow-mode shock event in the near-Earth magnetotail and a slow-mode shock detection percentage of 10% or lower in the mid-to-distant magnetotail. It was observed that if the high-energy beam region data is removed from the slow-mode shock downstream observations then the detection of slow-mode shocks reduces, implying that the kinetic effects play an important role in the detection of slow-mode shocks. For the crossings where the interface was not identified as a slow-mode shock, it was found that the turbulence in those crossings can change the mass flux values and disrupt the detection of slow-mode shock. However, the macroscopic slow-mode shock-like structure stably exists around the magnetic reconnection interface, as most of the conditions for slow-mode shocks were satisfied. This result suggests that slow-mode shocks are a general feature of magnetic reconnection geometry. We find that the lack of detection of slow-mode shocks in previous observations and simulations can be explained by taking into account the kinetic structure of slow-mode shocks and the presence of turbulence.
79 ASTRONOMY AND ASTROPHYSICS↗
Hybrid simulations of the saturated Farley-Buneman instability in the ionosphere
Numerical simulations of the Farley-Buneman instability in 2-1/2 dimensions using particle ions and fluid electrons show the growth, saturation and nonlinear behavior of two-stream waves. This hybrid technique models the saturated state of the instability for a much longer period of time than the pure particle codes that preceded it. While focusing principally on modeling the topside E region equatorial electrojet, many of these results apply to the auroral two-stream instability as well. The following features are seen in all our hybrid simulations: (1) wave growth at an angle offset from the electron drift direction where the angle depends on the strength of the driving electric field, (2) nonlinear coupling to waves traveling perpendicular to the propagation direction of the principal two-stream waves, (3) a saturated wave phase velocity at or above the sound speed but well below the velocity predicted by linear theory and (4) phase velocities which remain almost constant as a simulated radar sweeps from a horizontal direction to nearly vertical. The nonlinear electron motion dominates the behavior of these waves. Further, these simulations indicate that ion kinetic effects are not essential for the saturation of the instability and that electron temperature effects have a minor impact on the final saturated state.
Oppenheim, Meers↗
First Lunar Wake Passage of ARTEMIS: Discrimination of Wake Effects and Solar Wind Fluctuations by 3D Hybrid Simulations
The spacecraft P1 of the new ARTEMIS (Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon's Interaction with the Sun) mission passed the lunar wake for the first time on February 13, 2010. We present magnetic field and plasma data of this event and results of 3D hybrid simulations. As the solar wind magnetic field was highly dynamic during the passage, a simulation with stationary solar wind input cannot distinguish whether distortions were caused by these solar wind variations or by the lunar wake; therefore, a dynamic real-time simulation of the flyby has been performed. The input values of this simulation are taken from NASA OMNI data and adapted to the P1 data, resulting in a good agreement between simulation and measurements. Combined with the stationary simulation showing non-transient lunar wake structures, a separation of solar wind and wake effects is achieved. An anisotropy in the magnitude of the plasma bulk flow velocity caused by a non-vanishing magnetic field component parallel to the solar wind flow and perturbations created by counterstreaming ions in the lunar wake are observed in data and simulations. The simulations help to interpret the data granting us the opportunity to examine the entire lunar plasma environment and, thus, extending the possibilities of measurements alone: A comparison of a simulation cross section to theoretical predictions of MHD wave propagation shows that all three basic MHD modes are present in the lunar wake and that their expansion governs the lunar wake refilling process.
Wiehle, S.↗
On the cross-field diffusion of ions in one- and two-dimensional hybrid simulations of collisionless shocks
It can be demonstrated analytically that under certain geometries used in numerical simulations of collisionless shocks in which there is at least one ignorable spatial coordinate, the transport of particles across the magnetic field is essentially zero. This notion is tested using one- and two-dimensional hybrid simulations (kinetic ions/fluid electrons). We find, as the theorem predicts, the particles treated kinetically are tied to the same field line on which they start.
Giacalone, Joe↗
Large-scale Linear Magnetic Holes with Magnetic Mirror Properties in Hybrid Simulations of Solar Wind Turbulence
Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the long-standing effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this Letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anisotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modeling of MHs.
79 ASTRONOMY AND ASTROPHYSICS↗
Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study
Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.
Convery, P. D.↗
Hybrid simulations of a tangential discontinuity-driven foreshock bubble formation in comparison with a hot flow anomaly formation
Hot flow anomalies (HFAs) and foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere-ionosphere system. Yet, their early formation mechanisms are still not fully understood. To investigate the formation of tangential discontinuity (TD)-driven FBs and HFAs, we use 2-D local hybrid simulations where a reflected or an injected warm foreshock ion beam can interact with a TD whose half-thickness is comparable to the ion inertial scale. We show that the foreshock ions perform a partial gyration within, or across, the TD. Bulk motion differences between partially-gyrating foreshock ions and fluid-electrons lead to the generation of currents. As the trigger, these foreshock-driven currents change the magnetic field topology around the TD and force the frozen-in solar wind plasma to redistribute along with the field lines, shaping the foreshock transient. This confirms a recently proposed kinetic formation model. The extent of the magnetic field direction change across the TD within the foreshock ion gyromotion determines the current profile and thus the type of foreshock transient that forms. For a thin TD, the foreshock ions generate a current that is much stronger on the upstream side than the downstream side, forming an FB with one upstream compressional boundary. For the same foreshock ion gyroradius and magnetic shear, a thick TD yields comparable foreshock-driven currents on the upstream and downstream sides, forming an HFA with two compressional boundaries. Our study suggests that the TD thickness is one of the factors that determine the formation of FBs and HFAs
Andrew Vu↗