Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field 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.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Magnetosonic ULF Waves With Anomalous Plasma–Magnetic Field Correlations: Standing Waves and Inhomogeneous Plasmas

Ultra-low frequency (ULF) wave observations across the heliosphere often rely on the sign of correlations between plasma (density/pressure) and magnetic field perturbations to distinguish between fast and slow magnetosonic modes. However, the assumptions behind this magnetohydrodynamic result are not always valid, particularly within the magnetosphere which is inhomogeneous and supports standing waves along the geomagnetic field. Through theory and a global simulation, we find both effects can result in anomalous plasma–magnetic field correlations. The interference pattern in standing waves can lead both body and surface magnetosonic waves to have different cross-phases than their constituent propagating waves. Furthermore, if the scale of gradients in the background are shorter than the wavelength or the waves are near-incompressible, then advection by the wave of inhomogeneities can overcome the wave's inherent sense of compression. These effects need to be allowed for and taken into account when applying the typical diagnostic to observations.

MHD Waves↗
NASA NTRS2009

Recent Developments in Ultra High Temperature Ceramics at NASA Ames

NASA Ames is pursuing a variety of approaches to modify and control the microstructure of UHTCs with the goal of improving fracture toughness, oxidation resistance and controlling thermal conductivity. The overall goal is to produce materials that can perform reliably as sharp leading edges or nose tips in hypersonic reentry vehicles. Processing approaches include the use of preceramic polymers as the SiC source (as opposed to powder techniques), the addition of third phases to control grain growth and oxidation, and the use of processing techniques to produce high purity materials. Both hot pressing and field assisted sintering have been used to make UHTCs. Characterization of the mechanical and thermal properties of these materials is ongoing, as is arcjet testing to evaluate performance under simulated reentry conditions. The preceramic polymer approach has generated a microstructure in which elongated SiC grains grow in the form of an in-situ composite. This microstructure has the advantage of improving fracture toughness while potentially improving oxidation resistance by reducing the amount and interconnectivity of SiC in the material. Addition of third phases, such as Ir, results in a very fine-grained microstructure, even in hot-pressed samples. The results of processing and compositional changes on microstructure and properties are reported, along with selected arcjet results.

Johnson, Sylvia M.↗
NASA NTRS2023

A Surface Radiation Balance Data Set from Siple Dome in West Antarctica for Atmospheric and Climate Model Evaluation

A field campaign at Siple Dome in West Antarctica during the austral summer 2019-2020 offers an opportunity to evaluate climate model performance, particularly cloud microphysical simulation. Over Antarctic ice sheets and ice shelves, clouds are a major regulator of the surface energy balance, and in the warm season their presence occasionally induces surface melt that can gradually weaken an ice shelf structure. This dataset from Siple Dome, obtained using transportable and solar-poweredequipment, includes surface energy balance measurements, meteorology and cloud remote sensing. To demonstrate how these data can be used to evaluate model performance, comparisons are made with meteorological reanalysis known to give generally good performance over Antarctica (ERA5). Surface albedo measurements show expected variability with observed cloud amount, and can be used to evaluate a model's snowpack parameterization. One case study discussed involves a squall with northerly winds, during which ERA5 fails to produce cloud cover throughout one of the days. A second case study illustrates how shortwave spectroradiometer measurements that encompass the 1.6-micron atmospheric window reveal cloud phase transitions associated with cloud lifecycle. Here, continuously precipitating mixed-phase clouds become mainly liquid water clouds from local morning through the afternoon, not reproduced by ERA5. We challenge researchers to run their various regional or global models in a manner that has the large-scale meteorology follow the conditions of this field campaign, compare cloud and radiation simulations with this Siple Dome dataset, and potentially investigate why cloud microphysical simulations or other model components might produce discrepancies with these observations.

cloud remote sensing↗
NASA NTRS1987

Nonuniform sampling techniques for antenna applications

A two-dimensional sampling technique, which can employ irregularly spaced samples (amplitude and phase) in order to generate the complete far-field patterns is presented. The technique implements a matrix inversion algorithm, which depends only on the nonuniform sampled data point locations and with no dependence on the actual field values at these points. A powerful simulation algorithm is presented to allow a real-life simulation of many reflector/feed configurations and to determine the usefulness of the nonuniform sampling technique for the copolar and cross-polar patterns. Additionally, an overlapped window concept and a generalized error simulation model are discussed to identify the stability of the technique for recovering the field data among the nonuniform sampled data. Numerical results are tailored for the pattern reconstruction of a 20-m offset reflector antenna operating at L-band. This reflector is planned to be used in a proposed measurement concept of large antenna aboard the Space Shuttle, whereby it would be almost impractical to accurately control the movement of the Shuttle with respect to the RF source in prescribed directions in order to generate uniform sampled points. Also, application of the nonuniform sampling technique to patterns obtained using near-field measured data is demonstrated. Finally, results of an actual far-field measurement are presented for the construction of patterns of a reflector antenna from a set of nonuniformly distributed measured amplitude and phase data.

Rahmat-Samii, Yahya↗
NASA NTRS1991

Far-field pattern of a coherently combined beam from large-scale laser diode arrays

The far-field pattern of a large-scale amplifier array (LSAA) consisting of a large number (2000) of diode laser amplifiers is numerically simulated, and the power collection efficiencies are determined. Random distributions of phase mismatches, misorientations, and element failures in the LSAA system are considered. Phase mismatches and misorientations of the element amplifiers are found to be the most critical parameters of those affecting the power-collection efficiency. Errors of 0.2 wavelength and 25 percent for phase and diffraction angle, respectively, cause a 10 percent reduction in power-collection efficiency. The results are used to evaluate the concept of space-laser power transmission. It is found that an overall transmission efficiency of 80 percent could be realized with a 5-m-diam. receiver at a distance of 10,000 km when an LSAA transmitter 6 m in diam. is aimed with state-of-the-art pointing accuracy.

Kwon, Jin H.↗
NASA NTRSDate not supplied

Towards a more realistic representation of NASA CERES-derived surface radiative fluxes during polar night: a comparison with the MOSAiC field campaign

The Arctic is one of the most sensitive regions of Earth to climate change, and yet remains one of the more difficult regions to both observe and simulate. The high latitude of the Arctic gives it some unusual features compared with much of Earth, such as large and variable sea ice cover, extensive mixed phase cloud cover, and the polar day/night cycle. These features significantly influence the radiative environment of the Arctic as well as the ability to observe it. The remoteness of the Arctic from most of civilization makes large volumes of in situ observations difficult to collect, and so satellite observations are a critical tool for observing the region, such as those from Clouds and the Earth’s Radiant Energy System (CERES). But validating the satellite radiative flux estimates is difficult because of the lack of in situ measurements. A useful remedy is utilizing measurements from field campaigns, such as the year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition during 2019-2020. The extensive high-quality surface radiative flux and meteorological measurements collected from MOSAiC provide a useful check on flux retrievals from CERES instruments. Previous work compared MOSAiC and CERES surface radiative fluxes during the 2020 melt season (April-September). Errors in CERES fluxes were related to errors in estimated atmospheric optical depth and surface albedo. We conduct a similar assessment of CERES surface flux retrievals using MOSAiC data, but for the months of October through March instead. As this time period is mostly during polar night and twilight, the focus of the study is on longwave surface fluxes. To better understand the reasons for errors in CERES fluxes, we use the large set of meteorological measurements also collected by MOSAiC, including surface temperature, thermodynamic vertical profiles, surface turbulent fluxes, and cloud properties. For example, one notable possible source of error in the CERES estimate of surface downwelling longwave flux is the estimate of cloud base height when compared with MOSAiC W-band radar measurements – an overestimate of cloud base height yields an underestimate of the downwelling flux on the order of tens of W m2.

J. Brant Dodson↗
NASA NTRS1993

A numerical study of wave propagation in a confined mixing layer by eigenfunction expansions

It is well known that the growth rate of instability waves of a two-dimensional free shear layer is reduced greatly at supersonic convective Mach numbers. In previous works, it has been shown that new wave modes exist when the shear layers are bounded by a channel due to the coupling effect between the acoustic wave modes and the motion of the mixing layer. The present work studies the simultaneous propagation of multiple stability waves using numerical simulation. It is shown here that the coexistence of two wave modes in the flow field can lead to an oscillatory growth of disturbance energy with each individual wave mode propagating linearly. This is particularly important when the growth rates of the unstable waves are small. It is also shown here that the propagation of two neutrally stable wave modes can lead to a stationary periodic structure of rms fluctuations. In the numerical simulations presented here the forced wave modes are propagating at same frequency, but with different phase velocities. In order to track the growth of each wave mode as it propagates downstream, a numerical method that can effectively detect and separate the contribution of the individual wave is given. It is demonstrated that by a least square fitting of the disturbance field with eigenfunctions the amplitude of each wave mode can be found. Satisfactory results as compared to linear theory are obtained.

Hu, Fang Q.↗
NASA NTRS1996

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗
NASA NTRS2005

Interstellar PAH in the Laboratory and in Space. What have we Learned from the New Generation of Laboratory and Observational Studies?

Polycyclic Aromatic Hydrocarbons (PAHs) are an important and ubiquitous component of carbon-bearing materials in space. PAHs are the best-known candidates to account for the IR emission bands (UIR bands) and PAH spectral features are now being used as new probes of the ISM. PAHs are also thought to be among the carriers of the diffuse interstellar absorption bands (DIBs). In the model dealing with the interstellar spectral features, PAHs are present as a mixture of radicals, ions and neutral species. PAH ionization states reflect the ionization balance of the medium while PAH size, composition, and structure reflect the energetic and chemical history of the medium. A major challenge for laboratory astrophysics is to reproduce (in a realistic way) the physical conditions that exist in the emission and/or absorption interstellar zones. An extensive laboratory program has been developed at NASA Ames to assess the physical and chemical properties of PAHs in such environments and to describe how they influence the radiation and energy balance in space and the interstellar chemistry. In particular, laboratory experiments provide measurements of the spectral characteristics of interstellar PAH analogs from the ultraviolet and visible range to the infrared range for comparison with astronomical data. This paper will focus on the recent progress made in the laboratory to measure the direct absorption spectra of neutral and ionized PAHs in the gas phase in the near-UV and visible range in astrophysically relevant environments. These measurements provide data on PAHs and nanometer-sized particles that can now be directly compared to astronomical observations. The harsh physical conditions of the IS medium - characterized by a low temperature, an absence of collisions and strong VUV radiation fields - are simulated in the laboratory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature (of the order of 100 K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase. Intrinsic band profiles and band positions of cold gas-phase PAHs can now be measured with high-sensitivity spectroscopy and directly compared to the astronomical data. The electronic bands measured for ionized PAH are found to be intrinsically broad (about 20 cm(sup -1)) while the bands associated with the neutral precursors are narrower (of the order of 2 - 10 cm(sup -1)). The laboratory data are discussed and compared with recent astronomical spectra of large and narrow DIBs and with the spectra of circumstellar environments of selected carbon stars and the implications for the interstellar PAH population are derived. Preliminary results also show that carbon nanoparticles are formed during the short residence time of the precursors in the plasma.

Salama, Farid↗
NASA NTRS1982

Computer simulation of inner magnetospheric dynamics for the magnetic storm of July 29, 1977

The Rice University convection model is applied to the early main phase of the July 29, 1977 magnetic storm through a computer implementation that self-consistently calculates electric fields, currents, and plasma distributions and velocities in the inner magnetosphere/ionosphere system. On the basis of solar wind parameters and AL index as inputs, the model predicts the injection of plasma sheet plasma to form a substantial storm time ring current whose total predicted strength agrees with the observed Dst index. The possibility that the magnetic field may be sufficiently inflated to make 60 deg field lines extend to the outer magnetosphere is examined. In the model, distortion of the plasma sheet inner edge by magnetospheric compression associated with the sudden commencement temporarily disrupts the normal Birkeland current pattern. Normal Birkeland currents and shielding reassert themselves in about an hour.

Wolf, R. A.↗
NASA NTRS1994

Electron acceleration at nearly perpendicular collisionless shocks. 3: Downstream distributions

Spacecraft observations at the Earth's bow shock and at interplanetary shocks have established that the largest fluxes of accelerated suprathermal electrons occur in so-called shock spike events immediately downstream of the shock ramp. Previous theoretical efforts have mainly focused on explaining upstream energetic electron beams. Here we investigate the general motion and acceleration of energetic electrons in a curved, nearly perpendicular shock by numerically integrating the orbits of solar wind halo electrons in shock fields generated by a hybrid simulation (core electron fluid and kinetic ions). Close to the angle Theta(sub Bn) = 90 degs between the upstream magnetic field and shock normal, the calculations result in a (perpendicular) temperature increase proportional to the magnetic field ratio and give the highest phase space densities in the overshoot. For a steep distribution, the temperature change can correspond to an enhancement of the distribution by several orders of magnitude. These results are in agreement with predictions from adiabatic mapping. With smaller angles Theta(sub Bn), the overshoot and downstream densities fall off quickly, because the adiabatic energy gain is less and fewer electrons transmit. The shock curvature also leads to an accumulation of electrons close to 90 degs. Without pitch angle scattering, energization is only significant within a few (approximately 5 to 10 degs) degrees of the point of tangency. However, shock spike events appear to be observed more easily and farther away from 90 degs. Given that over a region of several degrees around 90 degs the theory gives enhancements of up to approximately 4 orders of magnitude, such electrons could in principle account for the typically observed enhancements of 1 to 2 orders of magnitude, if they were distributed over Theta(sub Bn). To test the idea that scattering could efficiently redistribute the energetic electrons, we have conducted test particle simulations in which artificial pitch angle scattering is included. We find that such a process can indeed be very efficient and can explain observations of shock spike events far away from Theta(sub Bn) is approximately 90 degs. The scattering naturally leads to much higher phase space densities at smaller Theta(sub Bn) than what a local one-dimensional mapping would predict and thus can account for an observed discrepancy with adiabatic theory stated in the literature.

Krauss-Varban, D.↗
NASA NTRS2005

PAHs and the Diffuse Interstellar Bands. What have we Learned from the New Generation of Laboratory and Observational Studies?

Polycyclic Aromatic Hydrocarbons (PAHs) are an important and ubiquitous component of carbon-bearing materials in space. PAHs are the best-known candidates to account for the IR emission bands (UIR bands) and PAH spectral features are now being used as new probes of the ISM. PAHs are also thought to be among the carriers of the diffuse interstellar absorption bands (DIBs). In the model dealing with the interstellar spectral features, PAHs are present as a mixture of radicals, ions and neutral species. PAH ionization states reflect the ionization balance of the medium while PAH size, composition, and structure reflect the energetic and chemical history of the medium. A major challenge for laboratory astrophysics is to reproduce (in a realistic way) the physical conditions that exist in the emission and/or absorption interstellar zones, An extensive laboratory program has been developed at NASA Ames to characterize the physical and chemical properties of PAHs in astrophysical environments and to describe how they influence the radiation and energy balance in space and the interstellar chemistry. In particular, laboratory experiments provide measurements of the spectral characteristics of interstellar PAH analogs from the ultraviolet and visible range to the infrared range for comparison with astronomical data. This paper will focus on the recent progress made in the laboratory to measure the direct absorption spectra of neutral and ionized PAHs in the gas phase in the near-W and visible range in astrophysically relevant environments. These measurements provide data on PAHs and nanometer-sized particles that can now be directly compared to astronomical observations. The harsh physical conditions of the IS medium - characterized by a low temperature, an absence of collisions and strong V W radiation fields - are simulated in the laboratory by associating a molecular beam with an ionizing discharge to generate a cold plasma expansion. PAH ions are formed from the neutral precursors in an isolated environment at low temperature (of the order of 100 K). The spectra of neutral and ionized PAHs are measured using the high sensitivity methods of cavity ring down spectroscopy (CRDS). These experiments provide unique information on the spectra of free, cold large carbon molecules and ions in the gas phase. Intrinsic band profiles and band positions of cold gas-phase PAHs can now be measured with high- sensitivity spectroscopy and directly compared to the astronomical data. The electronic bands measured for ionized PAH are found to be intrinsically broad (about 20/cm) while the bands associated with the neutral precursors are narrower (of the order of 2 - 10/cm).

Salama, Farid↗
NASA NTRS2015

Core Analysis Combining MT (TIPPER) and Dielectric Sensors (Sans EC) in Earth and Space

On terrestrial planets and moons of our solar system cores reveal details about a geological structure's formation, content, and history. The strategy for the search for life is focused first on finding water which serves as a universal solvent, and identifying the rocks which such solvent act upon to release the constituent salts, minerals, ferrites, and organic compounds and chemicals necessary for life. Dielectric spectroscopy measures the dielectric properties of a medium as a function of frequency. Reflection measurements in the frequency range from 300 kHz to 300 MHz were carried out using RF and microwave network analyzers interrogating SansEC Sensors placed on clean geological core samples. These were conducted to prove the concept feasibility of a new geology instrument useful in the field and laboratory. The results show that unique complex frequency spectra can be acquired for a variety of rock core samples. Using a combination of dielectric spectroscopy and computer simulation techniques the magnitude and phase information of the frequency spectra can be converted to dielectric spectra. These low-frequency dielectric properties of natural rock are unique, easily determined, and useful in characterizing geology. TIPPER is an Electro-Magnetic Passive-Source Geophysical Method for Detecting and Mapping Geothermal Reservoirs and Mineral Resources. This geophysical method uses distant lightning and solar wind activity as its energy source. The most interesting deflections are caused by the funneling of electrons into more electrically conductive areas like mineralized faults, water or geothermal reservoirs. We propose TIPPER to be used with SansEC for determining terrain/ocean chemistry, ocean depth, geomorphology of fracture structures, and other subsurface topography characteristics below the ice crust of Jovian moons. NASA envisions lander concepts for exploration of these extraterrestrial icy surfaces and the oceans beneath. One such concept would use a nuclear powered heated tip for melting through the ice sheath of Europa and inserting a down hole SansEC with TIPPER interface. NASA's Juno space probe already on the way to Jupiter as part of the Exploration New Frontiers Program and the planned Europa mission will conduct detailed reconnaissance of Jupiter's moon Europa and investigate whether the icy moon could harbor conditions suitable for life. It has already been observed that Jovian moons have auroras that may serve as naturally occurring active energy sources for a TIPPER instrument.

Mound, Michael C.↗
NASA NTRS1989

The formation of topological defects in phase transitions

It was argued, and fought through numerical work that the results of non-dynamical Monte Carlo computer simulations cannot be applied to describe the formation of topological defects when the correlation length at the Ginzburg temperature is significantly smaller than the horizon size. To test the current hypothesis that infinite strings at formation are essentially described by Brownian walks of size the correlation length at the Ginzburg temperature, fields at the Ginzburg temperature were equilibrated. Infinite structure do not exist in equilibrium for reasonable definitions of the Ginzburg temperature, and horizons must be included in a proper treatment. A phase transition, from small-scale to large-scale string or domain wall structure, is found to occur very close to the Ginzburg temperature, in agreement with recent work. The formation process of domain walls and global strings were investigated through the breaking of initially ordered states. To mimic conditions in the early Universe, cooling times are chosen so that horizons exist in the sample volume when topological structure formation occurs. The classical fields are evolved in real-time by the numerical solution of Langevin equations of motion on a three dimensional spatial lattice. The results indicate that it is possible for most of the string energy to be in small loops, rather than in long strings, at formation.

Hodges, Hardy M.↗
NASA NTRS2004

A Critical Review of the Transport and Decay of Wake Vortices in Ground Effect

This slide presentation reviews the transport and decay of wake vortices in ground effect and cites a need for a physics-based parametric model. The encounter of a vortex with a solid body is always a complex event involving turbulence enhancement, unsteadiness, and very large gradients of velocity and pressure. Wake counter in ground effect is the most dangerous of them all. The interaction of diverging, area-varying, and decaying aircraft wake vortices with the ground is very complex because both the vortices and the flow field generated by them are altered to accommodate the presence of the ground (where there is very little room to maneuver) and the background turbulent flow. Previous research regarding vortex models, wake vortex decay mechanisms, time evolution within in ground effect of a wake vortex pair, laminar flow in ground effect, and the interaction of the existing boundary layer with a convected vortex are reviewed. Additionally, numerical simulations, 3-dimensional large-eddy simulations, a probabilistic 2-phase wake vortex decay and transport model and a vortex element method are discussed. The devising of physics-based, parametric models for the prediction of (operational) real-time response, mindful of the highly three-dimensional and unsteady structure of vortices, boundary layers, atmospheric thermodynamics, and weather convective phenomena is required. In creating a model, LES and field data will be the most powerful tools.

Sarpkaya, T.↗
NASA NTRS1983

Cumulus clouds - Numerical models, observations and entrainment

The first computer simulation of the organization phase of a buoyant atmospheric thermal is described. Although crude, it showed the spontaneous development of a rounded tight-gradient 'cap' and internal vortical circulation. The complexities involved in these 'field of motion' models in part motivated the development of entity models, based upon laboratory thermals. These one-dimensional models and their uses with observations are briefly described as well as their limitations. Finally, an application of Schlesinger's three-dimensional model to a GATE cumulus situation clarifies many apparently conflicting observations and postulates, thereby raising further challenging questions to be addressed jointly by the more sophisticated measuring and modeling tools available in the 1980's.

Simpson, J.↗
NASA NTRS2004

Part 1 of a Computational Study of a Drop-Laden Mixing Layer

This first of three reports on a computational study of a drop-laden temporal mixing layer presents the results of direct numerical simulations (DNS) of well-resolved flow fields and the derivation of the large-eddy simulation (LES) equations that would govern the larger scales of a turbulent flow field. The mixing layer consisted of two counterflowing gas streams, one of which was initially laden with evaporating liquid drops. The gas phase was composed of two perfect gas species, the carrier gas and the vapor emanating from the drops, and was computed in an Eulerian reference frame, whereas each drop was tracked individually in a Lagrangian manner. The flow perturbations that were initially imposed on the layer caused mixing and eventual transition to turbulence. The DNS database obtained included transitional states for layers with various liquid mass loadings. For the DNS, the gas-phase equations were the compressible Navier-Stokes equations for conservation of momentum and additional conservation equations for total energy and species mass. These equations included source terms representing the effect of the drops on the mass, momentum, and energy of the gas phase. From the DNS equations, the expression for the irreversible entropy production (dissipation) was derived and used to determine the dissipation due to the source terms. The LES equations were derived by spatially filtering the DNS set and the magnitudes of the terms were computed at transitional states, leading to a hierarchy of terms to guide simplification of the LES equations. It was concluded that effort should be devoted to the accurate modeling of both the subgridscale fluxes and the filtered source terms, which were the dominant unclosed terms appearing in the LES equations.

Okong'o, Nora A.↗
NASA NTRS2023

Parameter Dependencies of Early-Stage Tangential Discontinuity-Driven Foreshock Bubbles in Local Hybrid Simulations

Foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere-ionosphere system. In the kinetic formation model, foreshock ions interact with the discontinuity by performing partial gyrations to generate currents that change the magnetic field topology around the discontinuity. However, how different foreshock ion properties affect the growth of the field variations is not well understood. Therefore, we use 2-D local hybrid simulations to study the effects of different foreshock ion distributions and properties on the growth of tangential discontinuity (TD)-driven FBs. We discover that for a gyrophase-bunched distribution with an initial phase where the guiding center is on the other side of the TD, the foreshock ions gyrate together across the TD, causing more foreshock ions to cross the TD and leading to a faster expansion of the structure than for a Maxwellian distribution. A ring distribution also yields higher expansion speeds because of the higher projected velocity into the new perpendicular direction. For Maxwellian distributions, there are positive and linear correlations of the FB expansion speeds with the initial foreshock ion densities, thermal speeds, parallel speeds, and sine of the TD magnetic shear angles. These parameter dependencies grow in strength as the structures evolve with time. The foreshock ion distributions and properties that lead to stronger currents produce more significant magnetic field variations and higher expansion speeds. Our study helps quantify the formation and expansion of FBs to forecast their space weather effects and contribution to shock acceleration.

Andrew Vu↗