Maxwellian Distribution Function Research Articles

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Analytical irreversible thermodynamics in conduction-radiative heat transfer of neutral gases through a lens of thermal radiation fields

An analytical investigation has been conducted on the steady-state conduction and a diluted gas exposed to a thermal radiation field experiences radiative heat transfer. Analytical solutions have been derived for the transport Boltzmann BGK partial differential equations system. The investigation delves into the irreversible thermodynamic behavior of the system through the utilization of the Liu-Lees model, which employs two-stream Maxwellian distribution functions. The moment approach is utilized to predict the behavior of macroscopic gas parameters like temperature, concentration, and Fourier heat flux. This process entails substituting these parameters into the two-stream Maxwellian distribution functions to analyze the system's non-equilibrium thermodynamic characteristics, which encompass gas particles and the heated plate. Through this analysis, kinetic coefficients, entropy, entropy flux, entropy production, and thermodynamic forces are determined. The investigation assesses the validation of Onsager's reciprocity relation, the Boltzmann H-theorem, the Le Chatelier principle, and the second law of thermodynamics within the system. Furthermore, the Gibbs formula is employed to estimate the ratios between various contributions to internal energy changes, which are determined by the total derivatives of extensive parameters. The non-equilibrium velocity distribution function was calculated for the first time, and its behavior was compared with that of the equilibrium one using 3D graphics. The behavior of the computed variables is predicted through graphical analysis and discussion of the results.

Benchmarking calculations of level energies, transition parameters and electron impact excitation cross sections of S-like Xe38+

Atomic structure parameters for levels corresponding to 2s22p63s23p4, 2s22p63s3p5, 2s22p63s23p23d2, 2s22p63s3p43d and 2s22p63s23p33d configurations of S-like Xe38+ are calculated using the fully relativistic multiconfiguration Dirac-Hatree-Fock (MCDHF) method followed by the subsequent relativistic configuration interaction (RCI) calculations. The parameters evaluated include energies of the lowest 75 levels and E1, M1, E2, and M2 transition parameters among these levels. The effect of the choice of virtual orbitals for generating the wavefunctions is discussed. The accuracy of our line strengths is established through the rigorous calculations of their associated uncertainties using three different methods. Another set of calculations using the many-body perturbation theory (MBPT) to validate the present MCDHF-RCI energies is carried out. Further, the electron impact excitation cross-sections using the relativistic distorted wave (RDW) theory are determined for all transitions from the ground and metastable states in the incident electron energy range from the excitation threshold to 10 keV. For plasma physics applications, the fitting parameters for these cross sections employing two distinct equations tailored for low and high-energy regimes are reported. Moreover, the rate coefficients are determined in the electron temperatures range of 15 eV to 100 eV by taking into account the Maxwellian electron energy distribution function and our computed cross-sections. In addition to addressing the atomic data gap in highly charged Xe ions, the present results are of good accuracy and can serve as benchmark tests for other theoretical evaluations. Moreover, they can also assist in line identification of the complex spectra of highly charged sulphur-like xenon ions.

DC-driven subatmospheric glow discharges in the infrared-stimulated

This paper presents a conceptual framework for experimental research combined with numerical analysis on direct current (DC) glow discharges in microscale planar gas discharge-semiconductor systems (GDSS). In the experimental section, several structural and elemental analyses, including SEM, EDAX, AFM, and near-infrared absorption spectra measurements were carried out for compound semiconductor zinc selenide (ZnSe) cathode sample. Argon (Ar) was charged into the plasma reactor cell of GDSS at pressures of 100 Torr subatmospheric and 760 Torr atmospheric, respectively, by a vacuum pump- gas filling station. Glow discharge light emissions from plasma, excited under three different intensity levels (dark, weak, strong) of infrared beam illumination on ZnSe cathode electrode, were measured by using a phomultiplier tube that is sensitive to UV–Visible wavelengths. In the numerical analysis section, simulation studies were carried out on the two-dimensional gas discharge-semiconductor microplasma system (GDSµPS) cell models using the finite-element method (FEM) solver COMSOL Multiphysics DC plasma program. Calculations and predictions were based on mixture-averaged diffusion drift theory and Maxwellian electron energy distribution function. GDSµPS cell was modeled in a square chamber with planar anode/cathode electrode pair coupled at a 50 μm discharge gap. Single side of ZnSe cathode was finely micro-digitated to increase the effective surface area for enhanced electron emission to the gas discharge cell. The electrical equivalent circuit (EEC) of the proposed model was driven by 1.0 kV DC voltage source. Binary Ar/H2 gas medium in a mixture of 3:2 molar ratio was introduced to the gas discharge chamber at constant 200 Torr subatmospheric pressure. Simulations were run for normal glow discharges to exhibit the electrical fast transient glow discharge behaviours from electron field emission state to self-sustained normal glow discharge state by numerically solving the electron density (ED), electron current density (ECD) and electric potential distribution (EPD) parameters.It is figured out that binary Ar/H2 gas discharge model can undertake a major role in shaping and controlling the spatiotemporal response to transient electro-optical behavior of microplasma-based artificial electromagnetic materials configured for high-efficiency infrared-to-visible wavelength conversion applications.

Gas Kinetic Scheme Coupled with High-Speed Modifications for Hypersonic Transition Flow Simulations.

The issue of hypersonic boundary layer transition prediction is a critical aerodynamic concern that must be addressed during the aerodynamic design process of high-speed vehicles. In this context, we propose an advanced mesoscopic method that couples the gas kinetic scheme (GKS) with the Langtry-Menter transition model, including its three high-speed modification methods, tailored for accurate predictions of high-speed transition flows. The new method incorporates the turbulent kinetic energy term into the Maxwellian velocity distribution function, and it couples the effects of high-speed modifications on turbulent kinetic energy within the computational framework of the GKS solver. This integration elevates both the transition model and its high-speed enhancements to the mesoscopic level, enhancing the method's predictive capability. The GKS-coupled mesoscopic method is validated through a series of test cases, including supersonic flat plate simulation, multiple hypersonic cone cases, the Hypersonic International Flight Research Experimentation (HIFiRE)-1 flight test, and the HIFiRE-5 case. The computational results obtained from these cases exhibit favorable agreement with experimental data. In comparison with the conventional Godunov method, the new approach encompasses a broader range of physical mechanisms, yielding computational results that closely align with the true physical phenomena and marking a notable elevation in computational fidelity and accuracy. This innovative method potentially satisfies the compelling demand for developing a precise and rapid method for predicting hypersonic boundary layer transition, which can be readily used in engineering applications.

High‐order gas kinetic flux solver for viscous compressible flow simulations

AbstractAlthough the gas kinetic schemes (GKS) have emerged as one of the powerful tools for simulating compressible flows, they exhibit several shortcomings. Since the local solution of continuous Boltzmann equation with the Maxwellian distribution function is used to calculate the numerical fluxes at the cell interface, the flux expression in GKS is usually more complicated. In this paper, a high‐order simplified gas kinetic flux solver (GKFS) is presented for simulating two‐dimensional compressible flows. Circular function‐based GKFS (C‐GKFS), which simplifies the Maxwellian distribution function into the circular function, combined with an improved weighted essentially non‐oscillatory (WENO‐Z) scheme is applied to capture more details of the flow fields with fewer grids. As a result, a simple high‐order accurate C‐GKFS is obtained, which improves the computing efficiency and reduce its complexity to facilitate the practical application of engineering. A series of benchmark‐test problems are simulated and good agreement can be obtained compared with the references, which demonstrate that the high‐order C‐GKFS can achieve the desired accuracy.

Toroidal Alfvén eigenmodes observed in low power JET deuterium–tritium plasmas

The Joint European Torus recently carried out an experimental campaign using a plasma consisting of both deuterium (D) and tritium (T). We observed a high-frequency mode using a reflectometer and an interferometer in a D-T plasma heated with low power neutral beam injection, . This mode was observed at a frequency and was located at major radii . The observed mode was identified as a toroidal Alfvén eigenmode (TAE) using the linear MHD code, MISHKA. Beam ions and fusion-born alpha particles were modelled using the full orbit particle tracking code LOCUST, which produces smooth distribution functions suitable for stability calculations without analytical fits or the use of moments. We calculated the stability of the 21 candidate modes using the HALO code. These calculations revealed that beam ions can drive TAEs with toroidal mode numbers with linear growth rates , while TAEs with n < 8 are damped by the beam ion population. Alpha particles drive modes with significantly smaller linear growth rates, due to the low alpha power generated almost exclusively by beam-thermal fusion reactions. Non-ideal effects were calculated using complex resistivity in the CASTOR code, leading to an assessment of radiative, collisional, and continuum damping for all 21 candidate modes. Ion Landau damping was modelled using Maxwellian distribution functions for bulk D and T ions in HALO. Radiative damping, the dominant bulk damping mechanism, suppresses modes with high toroidal mode numbers. Comparing the drive from energetic particles with damping from thermal particles, we find all but one of the candidate modes are damped. The single net-driven n = 9 TAE with a net growth rate matches experimental observations with a lab frequency and location . The TAE was driven by co-passing particles through the resonance. Both co- and counter-passing alpha particles drive the TAE through the resonance. Additional sideband resonances contribute significant drive for both beam and alpha particles.

Adaptive wave-particle decomposition in UGKWP method for high-speed flow simulations

With wave-particle decomposition, a unified gas-kinetic wave-particle (UGKWP) method has been developed for multiscale flow simulations. With the variation of the cell Knudsen number, the UGKWP method captures the transport process in all flow regimes without the kinetic solver’s constraint on the numerical mesh size and time step being determined by the kinetic particle mean free path and particle collision time. In the current UGKWP method, the cell Knudsen number, which is defined as the ratio of particle collision time to numerical time step, is used to distribute the components in the wave-particle decomposition. The adaptation of particles in the UGKWP method is mainly for the capturing of the non-equilibrium transport. In this aspect, the cell Knudsen number alone is not enough to identify the non-equilibrium state. For example, in the equilibrium flow regime with a Maxwellian distribution function, even at a large cell Knudsen number, the flow evolution can be still modelled by the Navier-Stokes solver. More specifically, in the near space environment both the hypersonic flow around a space vehicle and the plume flow from a satellite nozzle will encounter a far field rarefied equilibrium flow in a large computational domain. In the background dilute equilibrium region, the large particle collision time and a uniform small numerical time step can result in a large local cell Knudsen number and make the UGKWP method track a huge number of particles for the far field background flow in the original approach. But, in this region the analytical wave representation can be legitimately used in the UGKWP method to capture the nearly equilibrium flow evolution. Therefore, to further improve the efficiency of the UGKWP method for multiscale flow simulations, an adaptive UGKWP (AUGKWP) method is developed with the introduction of an additional local flow variable gradient-dependent Knudsen number. As a result, the wave-particle decomposition in the UGKWP method is determined by both the cell and gradient Knudsen numbers, and the use of particles in the UGKWP method is solely to capture the non-equilibrium flow transport. The current AUGKWP method becomes much more efficient than the previous one with the cell Knudsen number only in the determination of wave-particle composition. Many numerical tests, including Sod shock tube, normal shock structure, hypersonic flow around cylinder, flow around reentry capsule, and an unsteady nozzle plume flow, have been conducted to validate the accuracy and efficiency of the AUGKWP method. Compared with the original UGKWP method, the AUGKWP method achieves the same accuracy, but has advantages in memory reduction and computational efficiency in the simulation for flows with the co-existing of multiple regimes.

Numerical Calculations of Charging Threshold at GEO Altitudes With Two Temperature Non‐Extensive Electrons

AbstractSurface charging at geosynchronous altitude is one of the major concerns for satellites and spacecrafts. Spacecraft anomalies are often associated with extreme surface charging events, especially during substorms in which the GEO plasma is better modeled as two temperatures non‐Maxwellian plasma. In such case, we employ two temperature q‐non‐extensive distribution function to determine the onset of spacecraft surface charging which becomes complex since many parameters control the surface charging. We developed a current balance equation which better explains the charging threshold in comparison to a Maxwellian distribution function. The effect of non‐extensive parameters, temperature and density ratio on the current balance equation has been explained. The modified current balance equation predicts the critical and anti‐critical temperatures for various space‐grade materials both analytically and numerically. A significant change is observed in the quantities characterizing the charging current, average yield and density ratio in the presence of non‐extensive two temperature electrons. The mechanism underlying different charging behaviors at or near the threshold is also indicated at various plasma parametric domains. Furthermore, the general conditions of potential jump are also obtained theoretically which predicts the sudden or smooth potential transition.

Synthetic radio images of structured GRB and kilonova afterglows

ABSTRACT In this paper, we compute and analyse synthetic radio images of gamma-ray bursts (GRBs) and kilonova afterglows. For modelling the former, we consider GRB170817A-inspired set of parameters, while for the latter, we employ ejecta profiles from numerical-relativity simulations. We find that the kilonova afterglow sky map has a doughnut-like structure at early times that becomes more ring-like at late times. This is caused by the fact that the synchrotron emission from electrons following Maxwellian distribution function dominates the early beamed emission, while emissions from electrons following power-law distribution is important at late times. For an on-axis observer, the image flux centroid moves on the image plane, initially away from the observer. The image sizes, we find, are the largest for equal mass merger simulations with the soft equation of state. The presence of a kilonova afterglow affects the properties inferred from the source sky map, even if the GRB afterglow dominates the total flux density. The main effect is the reduction of the mean apparent velocity of the source, and an increase in the source size. However, this effect becomes negligible around the light curve peak of the GRB afterglow. For a far off-axis observer, neglecting the presence of the kilonova afterglow may lead to systematic errors in the inference of GRB properties from the sky map observations. Notably, at the observing angle inferred for GRB170817A, the presence of kilonova afterglow would affect the sky map properties, only at very late times $t\gtrsim 1500\,$ d.

Sensitivity of synchrotron radiation to the superthermal electron population in mildly relativistic plasma

Synchrotron radiation has markedly different behavior in ∼ 10 keV and in ∼ 100 keV plasma. We show that high-energy electrons that occupy the tail of velocity distribution function have disproportionate impact on power loss of ∼ 100 keV plasma. If electrons with energy more than cutoff energy are redistributed while keeping the Maxwellian distribution function below cutoff energy intact, both emission and absorption of synchrotron radiation act to decrease the lost power. These novel radiation transport effects in non-equilibrium plasma suggest large utility in the deconfinement of high-energy electrons to reduce synchrotron radiation in applications where the radiation is deleterious.

Third-order accurate 13-moment equations for non-continuum transport phenomenon

The derivation of analytical equations of non-continuum macroscopic transport phenomena is underpinned by approximate descriptions of the particle distribution function and is required due to the inability of the Navier–Stokes equations to describe flows at high Knudsen number (Kn ∼ 1). In this paper, we present a compact representation of the second-order correction to the Maxwellian distribution function and 13-moment transport equations that contain fewer terms compared to available moment-based representations. The intrinsic inviscid and isentropic assumptions of the second-order accurate distribution function are then relaxed to present a third-order accurate representation of the distribution function, using which corresponding third-order accurate moment transport equations are derived. Validation studies performed for Grad’s second problem and the force-driven plane Poiseuille flow problem at non-zero Knudsen numbers for Maxwell molecules highlight an improvement over results obtained by using the Navier–Stokes equations and Grad’s 13-moment (G13) equations. To establish the ability of the proposed equations to accurately capture the bulk behavior of the fluid, the results of Grad’s second problem have been validated against the analytical solution of the Boltzmann equation. For the planar Poiseuille flow problem, validations against the direct simulation Monte Carlo method data reveal that, in contrast to G13 equations, the proposed equations are capable of accurately capturing the Knudsen boundary layer.

Three-dimensional coupling of electron cyclotron drift instability and ion–ion two stream instability

Electron cyclotron drift instability (ECDI) and ion–ion two stream instability (IITSI) are both kinetic instabilities that can be present in low-temperature, partially magnetized plasmas. The coupling of instabilities in a three-dimensional configuration leads to the existence of more than one unstable roots to the kinetic dispersion relation. In this paper, a generalized method has been developed for numerically evaluating solutions to the three-dimensional dispersion relation for coupled ECDI and IITSI, assuming cold singly and doubly charged ions and a Maxwellian velocity distribution function for the electrons. The present study demonstrates the coupling between ECDI and IITSI that affects the most unstable mode as a function of the wavenumbers in three dimensions and various plasma properties, including the applied electric field, magnetic field, electron temperature, ion velocities, and plasma density. One of the most notable results is that, while the most unstable mode with the largest growth rate is in the direction of the E×B drift in the two-dimensional cases, the most unstable mode for the three-dimensional configuration occurs in the oblique direction between the applied electric field and the E×B drift. This agrees with experimental observations in cross field plasma sources using coherent Thomson scattering.

Wave–particle interaction in contact with a chaotic thermostat

This numerical investigation explores the properties of the wave–particle interaction in a situation in which the particle is simultaneously in contact with a chaotic thermostat [G. J. Morales, Phys. Rev. E 99, 062218 (2019)]. The role of the thermostat is to establish a Maxwellian velocity distribution function through deterministic chaotic orbits. The particle response is quantified by calculating the complex mobility, μk, ω, from the numerically obtained orbits for a wave of constant amplitude, with wave number k and frequency ω. It is found that in the limit of weak coupling to the thermostat, the behavior is that predicted by the plasma dispersion function, which implies collisionless Landau damping. As the coupling to the thermostat is increased (equivalent to increasing collisionality), the behavior follows the generalized collisional plasma dispersion function [Fried et al., Phys. Fluids 9, 292 (1966)]. For strong coupling, the response agrees with the Braginskii mobility. The nonlinear mobility associated with intermittent particle trapping is obtained for the various collisional regimes.

ON PERFORMANCE OF ACCEPTANCE SAMPLING PLANS USING SEQUENTIAL PROBABILITY RATIO TEST BASED ON TRUNCATED LIFE TESTS USING SOME DISTRIBUTION FUNCTION

Sequential sampling procedure constitute a powerful tool to achieve experimental efficiency a sequence of sample are taken from the lot and allow the number of samples to be determined entirely by the result of the sampling process, it forms the basis of numerous sequential technique for different applications which was the special case for multiple acceptance plans, the performance of sequential probability ratio test provides better detection for a wide range of source strength and useful in choosing the best statistical test, which take less time to make decision, was developed on the software reliability measure based on failure intensity for the life time random variable for a truncated life test. In this paper, we considered the performance of acceptance sampling plans using sequential probability ration (SPRT) based on truncated life tests on Maxwell distribution and exponential distribution function in terms of obtaining the minimum number of sample sizes necessary to obtain specified average life time under a given consumer’s and producer’s risk. The SPRT provide an alternative to fix the plans that can help diminish producer and consumer risk of reaching at a wrong decision.

Theoretical analysis of thermo-responsive behavior of microgels loaded with silver nanoparticles

Temperature is one of the extremely significant factors which is responsible for modifications mostly in physical but also in chemical properties of some particular matter. In present study, we carried out theoretical analysis of thermo-responsive behavior of Ag-p (N-isopropylmethacrylamide) hybrid microgels. We employed the classical Maxwellian distribution function with respect to the velocity of spherical silver nanoparticles inside the microgel’s network. The behavior of distribution function relative to the temperature represents that the distribution function still obeys the Maxwell-Boltzmann law even at low temperatures. Further, the diffusion coefficient, probability density, surface tension and surface energy of silver nano-particles are numerically evaluated and then graphically illustrated.

Gyrokinetic modelling of the Alfvén mode activity in ASDEX Upgrade with an isotropic slowing-down fast-particle distribution

In the present paper, the evolution of the Alfvén modes (AMs) is studied in a realistic ASDEX Upgrade equilibrium by analyzing the results of simulations with the global, electromagnetic, gyrokinetic particle-in-cell code ORB5. The energetic particles (EPs) are modelled both via the newly implemented isotropic slowing-down and with Maxwellian distribution functions. The comparison of the numerical results shows that modelling the EPs with the equivalent Maxwellian rather than with the slowing-down, does not significantly affect the frequency of the driven AM, while its growth rate appears to be underestimated with a quantitative difference as large as almost 30%. Additionally the choice of the isotropic slowing-down allows a better description of the nonlinear modification of the dominant AM frequency, while an equivalent Maxwellian underestimates it. A good comparison with the experimental spectrogram is found.

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

AbstractSpacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E &lt; tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E &lt; 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 &lt; E &lt; 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

Bernstein–Greene–Kruskal Ion Modes in Dusty Space Plasmas Application in Saturn’s Magnetosphere

Frequent observations of ion beams moving out from Saturn’s plasma environment hints at the generation of ion Bernstein–Greene–Kruskal (BGK) modes. As the plasma environments of Saturn and its moon Enceladus are characterized by the ubiquitous presence of massive negatively charged dust particles, the existing BGK theory for electron-ion plasma models cannot address this scenario. This manuscript develops a theoretical model for studying ion BGK modes in dusty plasmas. The analysis reveals that the presence of dust in the plasma enhances the stability of BGK modes. As the dust density increases, the effect of other parameters on stability, such as the electron temperature, becomes negligible. The model is developed by assuming that electrons and ions follow a kappa distribution, featuring a long tail trend in the superthermal component, in agreement with observations. Different scenarios with either electrons or ions obeying a Maxwell or kappa distribution function have been considered. A thorough analysis of the trapped ion distribution function considering various combinations indicates that a plasma where electrons are in thermal equilibrium and ions follow kappa distribution is the least favorable system for the generation of BGK modes.

Influence of the electron kinetics on Ar/NF3 inductively coupled plasma

A self-consistent two-dimensional axisymmetric fluid plasma model coupled with Maxwell's equations is used to analyze the plasma kinetics of inductively coupled Ar/NF3 low-pressure plasma. We have developed a plasma chemical reaction mechanism for Ar/NF3 plasma and validated the same against the experimental data found in the literature. We analyze two sets of the electron-neutral reaction rate coefficients obtained for assumed Maxwellian and non-Maxwellian electron energy distribution functions. We find that the electron density in the discharge is sensitive to the choice of the electron energy distribution, while the fluorine atoms density is less sensitive. Since the density of atoms is controlled by a balance between production by dissociation reactions and diffusion to the walls, their density is very sensitive to the wall sticking probability.