Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic electromagnetic model”

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 55 records · Page 3

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi Ion Magnetosphere - 1

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2003] is presented. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate spatial, temporal, and spectral evolutions of the ring current and electromagnetic ion cyclotron waves. To demonstrate the effects of EMIC wave propagation and refraction on the EMIC wave energy distributions and evolution we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, due to the density gradient at the plasmapause, the net wave refraction is suppressed, and He(+)-mode grows preferably at plasmapause. This result is in a total agreement with the previous ray tracing studies, and very clear observed in presented B-field spectrograms. Second, comparison the global wave distributions with the results from other ring current model [Kozyra et al., 1997] reveals that our model provides more intense and higher plasmapause organized distributions during the May, 1998 storm period. Finally, the found He(+)-mode energy distributions are not Gaussian distributions, and most important that wave energy can occupy not only the region of generation, i. e. the region of small wave normal angles, but the entire wave normal angle region and even only the region near 90 degrees. The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping, and subsequent downward heat transport and excitation of stable auroral red arcs.

Khazanov, G. V.↗

Transport at high ${\beta_p}$ and development of candidate steady state scenarios for ITER

On DIII-D, the high β p scenario has an internal transport barrier (ITB), β N ~β p ~3,q 95 ~10, and very high normalized confinement H 98,y2 ~1.6. Recently, plasmas starting with these conditions have been dynamically driven to q 95 ~6 and β p ~2, where we find the ITB and high performance persist for five energy confinement times. These conditions are projected to meet the ITER steady-state goal of Q = 5. The ITB is maintained at lower β p with a strong reverse shear, consistent with predictions that negative central shear can lower the β p threshold for the ITB. There are two observed confinement states in the high β p scenario: H-mode confinement state with a high edge pedestal, and an enhanced confinement state with a low pedestal and an ITB. It has been observed in a scan of external resonant magnetic perturbation amplitude that when there are no large type-I ELMs, there is no transition to enhanced confinement. This is consistent with the proposed mechanism for ITB formation being a type-I ELM. Quasilinear gyro-Landau fluid predictive modeling of ITER suggests that only a modest reverse shear is required to achieve the ITB formation necessary for Q=5 when electromagnetic physics including the kinetic ballooning mode (KBM) is incorporated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves (Khazanov et al., 2003) is presented In order to adequately take into account wave propagation and refraction in a multi-ion magnetosphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate the spatial, temporal, and spectral evolution of the ring current and of electromagnetic ion cyclotron waves To demonstrate the effects of EMIC wave propagation and refraction on the wave energy distribution and evolution, we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, owing to the density gradient at the plasmapause, the net wave refraction is suppressed, and He+-mode grows preferably at the plasmapause. This result is in total agreement with previous ray tracing studies and is very clearly found in presented B field spectrograms. Second, comparison of global wave distributions with the results from another ring current model (Kozyra et al., 1997) reveals that this new model provides more intense and more highly plasmapause-organized wave distributions during the May 1998 storm period Finally, it is found that He(+)-mode energy distributions are not Gaussian distributions and most important that wave energy can occupy not only the region of generation, i.e., the region of small wave normal angles, but all wave normal angles, including those to near 90 . The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping and subsequent downward heat transport and excitation of stable auroral red arcs.

Khazanov, G. V.↗

Absorber Optimization Study for the Terahertz Intensity Mapper (TIM)

We discuss the design and optimization of the absorber for the long-wavelength arm of the Terahertz Intensity Mapper (TIM), a balloon-borne spectrometer employing Kinetic Inductance Detectors (KIDs). Electromagnetic simulations of our design indicate in-band absorption efficiency over 80% in both linear polarization modes. By developing custom transmission line model and mode-matching calculations, we find the absorption efficiency is affected by the absorber’s reactive part and overall shape. These insights into the operation of this design provide guidance for its optimization for low-resistance absorber. materials.

Nie, Rong↗

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radial drift of plasma blobs in a toroidal magnetic field with fully kinetic and reduced fluid models

Abstract In curved magnetic geometries, field-aligned regions of enhanced plasma pressure and density, termed ‘blobs,’ move as coherent filaments across the magnetic field lines. Coherent blobs account for a significant fraction of transport at the edges of magnetic fusion experiments and arise in naturally-occurring space plasmas. This work examines the dynamics of blobs with a fully kinetic electromagnetic particle-in-cell code and with a drift-reduced fluid code. In low-beta regimes with moderate blob speeds, good agreement is found in the maximum blob velocity between the two simulation schemes and simple analytical estimates. The fully kinetic code demonstrates that blob speeds saturate near the initial sound speed, which is a regime outside the validity of the reduced fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.

Max, C. E.↗

Cluster/Peace Electrons Velocity Distribution Function: Modeling the Strahl in the Solar Wind

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the electron heat-flux and temperature anisotropy. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by fitting the VDF's to a spectral spherical harmonic model from which the moments are derived directly from the model coefficients.

Figueroa-Vinas, Adolfo↗

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic ICWs. Initial Results: Waves and Precipitation Fluxes

Initial results from the new developed model of the interacting ring current ions and ion cyclotron waves are presented. The model described by the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another one gives wave evolution. Such system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. Calculating ion-wave relationships, on a global scale under non steady-state conditions during May 2-5, 1998 storm, we presented the data at three time cuts around initial, main, and late recovery phases of May 4, 1998 storm phase. The structure and dynamics of the ring current proton precipitating flux regions and the wave active ones are discussed in detail.

Khazanov, G. V.↗

Dark matter, dark radiation and gravitational waves from mirror Higgs parity

An exact parity replicates the Standard Model giving a Mirror Standard Model, SM ↔ SM ' . This “Higgs Parity” and the mirror electroweak symmetry are spontaneously broken by the mirror Higgs, ( H ' ) = v ' >> ( H ), yielding the Standard Model Higgs as a Pseudo-Nambu-Goldstone Boson of an approximate SU (4) symmetry, with a quartic coupling λ SM ( v ' ) ~ 10 - 3 . Mirror electromagnetism is unbroken and dark matter is composed of e ' and \( {\overline{e}}^{\prime } \) . Direct detection may be possible via the kinetic mixing portal, and in unified theories this rate is correlated with the proton decay rate. With a high reheat temperature after inflation, the e t dark matter abundance is determined by freeze-out followed by dilution from decays of mirror neutrinos, ν ' → ℓH . Remarkably, this requires v ' ~ (10 8 –10 10 ) GeV, predicting a Higgs mass of 123 ± 3 GeV at 1 σ and a Standard Model neutrino mass of (10 - 2 –10 - 1 ) eV, consistent with observed neutrino masses. The mirror QCD sector exhibits a first order phase transition producing gravitational waves that may be detected by future observations. Mirror glueballs decay to mirror photons giving dark radiation with Δ N eff ~ 0 . 03–0 . 4. With a low reheat temperature after inflation, the e ' dark matter abundance is determined by freeze-in from the SM sector by either the Higgs or kinetic mixing portal.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic Ion Cyclotron Waves, Initial Results: Waves and Precipitating Fluxes

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. These equations for the ion phase space distribution function and for the wave power spectral density were solved on aglobal magnetospheric scale undernonsteady state conditions during the 2-5 May 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the ion cyclotron wave-active zones during extreme geomagnetic disturbances on 4 May 1998 are presented and discussed in detail.

Khazanov, G. V.↗

The Nonlinear Coupling of Electromagnetic Ion Cyclotron and Lower Hybrid Waves in the Ring Current Region: The Magnetic Storm May 1-7 1998

The excitation of lower hybrid waves (LHWs) is a widely discussed mechanism of interaction between plasma species in space, and is one of the unresolved questions of magnetospheric multi-ion plasmas. In this paper we present the morphology, dynamics, and level of LHW activity generated by electromagnetic ion cyclotron (EMIC) waves during the May 2-7, 1998 storm period on the global scale. The LHWs were calculated based on our newly developed self-consistent model that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes the evolution of EMIC waves. It is found that the LHWs are excited by helium ions due to their mass dependent drift in the electric field of EMIC waves. The level of LHW activity is calculated assuming that the induced scattering process is the main saturation mechanism for these waves. The calculated LHWs electric fields are consistent with the observational data.

Khazanov, G. V.↗

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large- and Small-Scale Ring Current Electrodynamic Coupling

In this talk we will address the two primary issues of ring current (RC) electrodynamic coupling: 1. RC self-consistent magnetosphere-ionosphere coupling that includes calculation of the magnetospheric electric field (large scale electrodynamic coupling); and 2. RC self-consistent coupling with electromagnetic ion cyclotron (EMIC) waves (small scale electrodynamic coupling). Our study will be based on two RC models that we have recently developed in our group. The first model by Khazanov et al. [2002] couples the system of two kinetic equations: one equation which describes the RC ion dynamics and another equation which describes the energy density evolution of EMIC waves. The second model by Khazanov et al. [2003] deals with large scale electrodynamic coupling processes and provides a self-consistent simulation of RC ions and the magnetospheric electric field. There is presently no model that addresses both of these issues simultaneously in a self-consistent calculation. However, the need exists for such a model, because these two processes directly influence each other, with the mesoscale coupling changing the drift paths of the thermal and energetic particle populations in the inner magnetosphere, thereby changing the wave interactions, and the microscale coupling altering the pitch angle distributions and ionospheric conductivities (through increased precipitation), thus changing the field-aligned currents and electric potential structure. The initial thrust of the work will be the development of a combined kinetic model of micro- and meso-scale RC electrodynamic coupling processes and to examine their interactions with each other on a global scale.

Khazanov, G. V.↗

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗

On the Ion Precipitation due to Field Line Curvature (FLC) and EMIC Wave Scattering and Their Subsequent Impact on Ionospheric Electrodynamics

Both field line curvature (FLC) and electromagnetic ion cyclotron (EMIC) wave scattering are believed to be associated with energetic ion precipitation, but their relative contributions to ionospheric ion precipitation and subsequent effects are still unclear. Here, in this study, by using a kinetic ring current model, we investigate their impact on the ionosphere from two aspects: the global distribution of ion precipitation and resulting ionospheric conductance. Our results show that the intensity and coverage of ion precipitation due to EMIC waves are larger than that due to FLC scattering, while the latter mostly contributes to ion precipitation in outer regions (L > 4–5). We then estimate the conductance with empirical models using the simulated ion precipitation energy flux. When the EMIC wave associated proton precipitation is included, the conductance is significantly enhanced in the dusk-to-midnight sector and has a wide magnetic latitude (MLAT) range from around 52° to 62°, considerably altering the electric potential in the dusk sector, which further influences particle dynamics in the magnetosphere. On the contrary, the proton precipitation caused by FLC scattering only occurs at higher latitudes above MLAT = 60° and the corresponding conductance is slightly enhanced at midnight, with negligible influence on the convective electric potential. Although electron precipitation-associated conductance is predominant globally, the results show that proton precipitation can also play an important role in ionospheric electrodynamics, especially when EMIC wave-associated precipitating proton energy flux exceeds that of electrons in the dusk sector, a region where many subauroral coupling processes occur.

79 ASTRONOMY AND ASTROPHYSICS↗