Particle Flux Research Articles

Optical measurements of precipitating relativistic electron microbursts during geomagnetic disturbance and pulsating aurora

Abstract To reliably predict the distribution of heat and particle fluxes at the target plates of tokamaks, a comprehensive understanding of turbulence throughout the entire Scrape-Off-Layer (SOL) is imperative. This study examines divertor turbulence systematically across a broad parameter range on the TCV tokamak, including variations in magnetic field direction, plasma current I p ∈ [ 140 , 320 ] kA, edge safety factor q 95 ∈ [ 2.6 , 4.7 ] and Greenwald fraction f G ∈ [ 0.18 , 0.6 ] . The TCV X-point Gas Puff Imaging (GPI) system is used to measure 2D filament properties in the inner and outer divertor region. The fluctuation levels in the divertor are found to strongly increase with density (to 80% over most of the SOL) while remaining insensitive to I p . The previously identified divertor-localized filaments (DLF), located on the bad curvature side of the outer divertor leg, are found to be a common feature on TCV, while no filaments are observed in the PFR. DLFs are present over most of the parameter space and in both field directions. However, they are absent, or appear only closer to the target, for sufficiently large Λ div ≳ 10 or q 95 ≳ 3.7 . Across both I p and f G scans, some clear trends with Λ div are found for divertor filament sizes and velocities, and with target fall-off lengths of density and heat flux profiles at the outer target. This study provides important experimental insights to turbulent transport in the divertor also for comparison with self-consistent, turbulence simulations and extrapolation to future reactor conditions.

The physical mechanisms usually applied to explain the relativistic electron enhancement have been delved into to elucidate non-adiabatic electron acceleration resulting in the ultra-relativistic electron population observed in the outer radiation belt. We considered multisatellite observations of the solar wind parameters, magnetospheric waves, and particle flux to report an unusual local acceleration of ultra-relativistic electrons under a prolonged high-speed solar wind stream (HSS). A corotating interaction region reaches the Earth’s bowshock on August 3, 2016, causing a minor geomagnetic storm. Following this, the magnetosphere was driven for 72 h by a long-term HSS propagating at 600 km/s. During this period, the magnetosphere sustained both ultra-low frequency (ULF) and very-low frequency (VLF) waves in the outer radiation belt region. Besides the waves, the relativistic and ultra-relativistic electron fluxes were enhanced with different time lags regarding the magnetic storm main phase. The efficiency of wave-particle interaction in enhancing ultrarelativistic electrons is evaluated by the diffusion coefficient rates, considering both ULF and VLF waves together with phase space density analyses. Results show that local acceleration by whistler mode chorus waves can occur in a time scale of 2–4 h, whereas ULF waves take around 10’s of hours and magnetosonic waves take a time scale of days. This result is confirmed by the phase space density analysis. Accordingly, it shows that peaks of local acceleration of 1 MeV electrons are consistent with the observation of the highest chorus wave amplitude at the same L-shell and MLT. Thus, we argue that whistler mode chorus waves interacting with relativistic electrons are the main physical mechanisms leading to ultra-relativistic electron enhancement, while ULF and fast magnetosonic waves are found as secondary physical processes. Lastly, our analysis contributes to understanding how whistler and ULF waves can contribute to ultra-relativistic electrons showing up in the inner magnetosphere under the HSS driver.

Abstract Researches on plasma-facing materials/components (PFMs/PFCs) have become a focus in magnetic confinement fusion studies, particularly for advanced tokamak operation scenarios. Similarly, spacecraft surface materials must maintain stable performance under relatively high temperatures and other harsh plasma conditions, making studies of their thermal and ablation resistance critical. Recently, a low-cost, low-energy-storage for superconducting magnets, and compact linear device, HIT-PSI, has been designed and constructed at Harbin Institute of Technology (HIT) to investigate the interaction between stable high heat flux plasma and PFMs/PFCs in scrape-off-layer (SOL) and divertor regions, as well as spacecraft surface materials. The parameters of the argon plasma beam of HIT-PSI are diagnosed using a water-cooled planar Langmuir probe and emission spectroscopy. As magnetic field rises to 2 T, the argon plasma beam generated by a cascaded arc source achieves high density exceeding 1.2×1021 m−3 at a distance of 25 cm from the source with electron temperature surpassing 4 eV, where the particle flux reaches 1024 m−2s−1, and the heat flux loaded on the graphite target measured by infrared camera reaches 4 MW/m2. Combined with probe and emission spectroscopy data, the transport characteristics of the argon plasma beam are analyzed.

Objective. Bragg peak measurements play a key role in the beam quality assurance in proton therapy. Used as base data for the treatment planning softwares, the accuracy of the data is crucial when defining the range of the protons in the patient.Approach. In this paper a protocol to reconstruct a Pristine Bragg Peak exploring the direct correlation between the particle flux and the dose deposited by particles is presented. Proton flux measurements at the HollandPTC and FLUKA Monte Carlo simulations are used for this purpose. This new protocol is applicable to plastic scintillator detectors developed for Quality Assurance applications. In order to obtain the Bragg curve using a plastic fiber detector, a PMMA phantom with a decoupled and moveable stepper was designed. The step phantom allows to change the depth of material in front of the fiber detector during irradiations. The Pristine Bragg Peak reconstruction protocol uses the measured flux of particles at each position and multiplies it by the average dose obtained from the Monte Carlo simulation at each position.Main results. The results show that with this protocol it is possible to reconstruct the Bragg Peak with an accuracy of about 470µm, which is in accordance with the tolerances set by the AAPM.Significance. It has the advantage to be able to overcome the quenching problem of scintillators in the high ionization density region of the Bragg peak.

GEANT4 simulation for controlling and focusing of laser-accelerated proton beam for particle therapy using pulsed power solenoids

A set of experimental and measurement techniques to study the influence of a plasma jet on the main material parameters of piezoelectric ceramics has been presented. A series of plasma experiments has been carried out using a pulsed plasma jet system. It allows of a metered-dose exposure to plasma of different composition and fluence with a constant particle flux density of 1021/m2, energy flux density of 0.1 MJ/m2 and average particle energy of 100–200 eV in a pulse duration of 15 μs. The study of the effects that a repeated exposure to an extreme heat flux of helium and hydrogen plasmas has on the near-surface layer structure and basic material parameters of mass-produced piezoelectric ceramic samples has been presented. The main result of the research is an experimental confirmation of the surface micro-structuring starting after just a few cycles of plasma exposure while only a slight decrease of the main material parameters as well as the preservation of polarization has been observed for two types of different compositions of PZT-ceramics. A further increase in the number of exposure pulses leads to practically no change of main material parameters of both ceramics, even showing a tendency for recovery instead.

Atmospheric aerosol particles are essential for forming clouds and precipitation, thereby influencing Earth’s energy budget, water cycle and climate on regional and global scales. However, the origin of aerosol particles over the Amazon rainforest during the wet season is poorly understood. Earlier studies showed new particle formation in the outflow of deep convective clouds and suggested a downward flux of aerosol particles during precipitation events. Here we use comprehensive aerosol, trace gas and meteorological data from the Amazon Tall Tower Observatory to show that rainfall regularly induces bursts of nanoparticles in the nucleation size range. This can be attributed to rain-related scavenging of larger particles and a corresponding reduction of the condensation sink, along with an ozone injection into the forest canopy, which could increase the oxidation of biogenic volatile organic compounds, especially terpenes, and enhance new particle formation. During and after rainfall, the nucleation particle concentrations directly above the canopy are greater than those higher up. This gradient persists throughout the wet season for the nucleation size range, indicating continuous particle formation within the canopy, a net upward flux of newly formed particles and a paradigm shift in understanding aerosol–cloud–precipitation interactions in the Amazon. Particle bursts provide a plausible explanation for the formation of cloud condensation nuclei, leading to the local formation of green-ocean clouds and precipitation. Our findings suggest that an interplay of a rain-related reduction in the condensation sink, primary emissions of gases, mainly terpenes, and particles from the forest canopy, and convective cloud processing determines the population of cloud condensation nuclei in pristine rainforest air.

The problem of ensuring the accuracy and traceability of the measurement results of the absorbed dose in carbon ion beams, as well as the measurement results of the amount, fluence, flux density and energy of particles in proton and heavy charged particles beams is considered. Until now, in practice, these values have been measured only by indirect methods. The lack of approved measuring instruments for the quantities under consideration and the metrological traceability of measurement results of these quantities to standards did not allow achieving consistency of measurement methods used in practice and confirming the reliability of the results obtained. To solve this problem, three measuring complexes have been developed and created, which are included in the State Primary Standard of units of absorbed dose and absorbed dose rate of photon, electron, proton radiation and in carbon ion beams, quantity, fluence, flux density and energy of particles in proton and heavy charged particles beams GET 38-2024. The measuring complex for reproducing the unit of absorbed dose in carbon ion beams consists of an adiabatic calorimeter, a thermostating system, a data collection and processing system and a vacuum pumping station. To reproduce the unit of energy of protons and heavy charged particles, a complex has been implemented, which includes a total absorption calorimeter, a data collection and processing system, a vacuum pumping station and a particle count determination system based on the use of a Faraday cup. To reproduce the units of fluence and particle flux density in proton and heavy charged particles beams, a measuring complex has been created containing a Faraday cup, a set of collimators and a low current meter. The schemes, the principles and the results of studies of the metrological characteristics of the developed measuring complexes are described. The results are relevant for the field of radiation therapy and radiation resistance tests of the electronic component base used in the space industry.

We conduct an Euler-Lagrange, direct numerical simulation of a turbulent channel flow at a shear Reynolds number of Reτ=180 over an erodible particle bed. The particle bed consists of approximately 1.3 × 106 monodisperse particles, resulting in a bed thickness of around 12–13 particles. The particle density and size are chosen to achieve a ratio of 4 for the Shields stress to the critical Shields stress necessary for incipient motion such that particle transport occurs primarily as bedload. The simulation is run long enough for ripples to form. We track the temporal evolution of the particle flux and excess Shields stress for the entire bed as well as for the four regions of a ripple, namely, the crest, trough, lee side, and stoss side. We find that the particle flux and excess Shields stress closely match the Wong and Parker correlation when the particle bed is featureless at early time but diverge from the correlation when ripples form. This deviation primarily arises from particle transport in the trough and lee side regions. Conversely, particle transport in the crest and stoss side regions remains largely consistent with the Wong and Parker correlation. A root mean square-based correction for the bed is proposed to be used in conjunction with the Wong and Parker correlation. Additionally, ripples attain a self-similar profile in the shape and near-bed shear stress when they are sufficiently distant from their upstream neighbor. Any departure from self-similarity occurs when the upstream neighbor gets within close proximity.

A study of the effect of toroidal magnetic field (Bt) on the divertor asymmetries is carried out with a plasma transport code under BOUT++ framework with the magnetic equilibria from EAST and C-Mod. For the simulation cases with drifts, the density is larger at the inner divertor target, while the temperature and heat flux are larger at the outer divertor target. The ratio of the total particle flux at the outer target to that of the inner target increases with increasing Bt, while the ratio of the total heat flux at the outer target to that of the inner target decreases with increasing Bt. The in–out divertor asymmetries of both total particle and heat fluxes get weaker with increasing Bt. Further analysis shows that the edge radial transport induced by drifts is much weaker for the simulation case with higher Bt, indicating that drift-driven divertor asymmetries may be less important for future tokamaks with high Bt.

Abstract Increasing wildfire frequency and intensity indicate the need for improved wildfire Quantitative Risk Assessment (QRA) methodologies for industrial components. Firebrand exposure, the leading cause of infrastructure wildfire ignition, has not been considered in wildfire QRA frameworks. This paper proposes a methodology to simulate and quantify firebrand contact exposure to an industrial storage tank using the Fire Dynamics Simulator (FDS). A low aspect-ratio cylindrical tank obstacle is exposed to constant Lagrangian particles flux (2.5 particles m−2s−1) under varying wind speeds (4ms−1, 7ms−1 and 10ms−1), with varying particle densities (50 kgm−3, 100 kgm−3 and 300 kgm−3). Constant-mass spherical particles simulate firebrands. The highest firebrand landing densities occur at the tank and floor intersection, in two regions located perpendicular to wind direction, and at the tank’s windward recirculation zone edge. The highest landing densities are reached with highest wind speed tested, and with 100 kgm−3 particle density, reaching 7 times the inlet flux in [particles m −2 ]. The highest firebrand amount that reaches the storage tank roof, especially vulnerable to ignition, occurs at 10 ms−1 wind with the lowest particle density tested (50 kgm−3). Results allow quantification of firebrand contact exposure around an industrial critical component; contributing to future comprehensive wildfire NaTech QRA methodologies.

The effect of CD4 injection on WD molecule sputtering at the divertor target in EAST

Abstract. The central Arabian Sea, a unique tropical basin, is profoundly impacted by monsoon wind reversal affecting its surface circulation and biogeochemistry. Phytoplankton blooms associated with high biological productivity and particle flux occur in the northern part of the central Arabian Sea due to summer-monsoon-induced open-ocean upwelling and winter convection. The core oxygen minimum zone (OMZ) at intermediate water depths is another important feature of the northern central Arabian Sea and fades southward. In this study, we attempt to interlink how these factors collectively impact phytodetrital export to the sediment. Short sediment core-top (1 cm) samples representing the recent particle flux signatures were analysed from five locations (21 to 11° N; 64° E) in the central Arabian Sea. Previously, we used core-top (0–0.5 cm) samples and observed a trend between diatom frustule abundance and diversity with bulk sedimentary parameters indicating a spatial variability in phytodetrital export to the sediment. To verify this observation further, lipid biomarkers of key phytoplankton groups and a sea surface temperature (SST) proxy have been analysed in addition to diatom frustules. The C37 alkenone-based SST proxy indicated cooler SST (27.6 ± 0.25 °C) in the north (21–15° N) mostly due to upwelling (summer) and convective mixing (winter). Warmer SSTs (+0.4 °C) are measured in the south, which usually remains nutrient-poor. This trend was consistent with satellite-derived average SST values (2017–2020). Lipid biomarker analysis suggests that dinoflagellates were likely to be the highest contributor, as indicated by dinosterol and its degradative product dinostanol, followed by brassicasterol and C37 alkenone, likely representing diatoms and coccolithophores, respectively. The north, which largely experiences periodic phytoplankton blooms and is influenced by the thick OMZ, revealed the highest contents of organic matter, diatom frustules (diversity and abundance), dominated by large, thickly silicified cells (e.g. Coscinodiscus and Rhizosolenia) and phytoplankton lipid biomarkers, as well as lower contents of zooplankton biomarkers (cholesterol and cholestanol). In contrast, relatively smaller chain-forming centric (e.g. Thalassiosira) and pennate (e.g. Pseudo-nitzschia, Nitzschia, Thalassionema) diatom frustules along with lower phytoplankton lipid biomarker contents were found in the south, where zooplankton biomarkers and silicious radiolarians were more abundant. The possible impacts of the OMZ on particle flux related to the phytoplankton community, including zooplankton grazing and other factors, have been discussed.

Abstract Large‐scale patterns of the steady‐state magnetosheath plasma flow and their dependence on the interplanetary magnetic field (IMF) have been reconstructed for the first time on the basis of large multi‐year multi‐mission pool of spacecraft observations, concurrent interplanetary data, and an empirical high‐resolution model. The flow model architecture builds upon a recently developed magnetosheath magnetic field representation by flexible expansions of its toroidal and poloidal components in a coordinate system, naturally conformed with the magnetopause and bow shock shapes. The model includes two physics‐based flow symmetry modes: the first one treats the magnetosphere as an axisymmetric unmagnetized obstacle, whereas the second mode takes into account the geodipole tilt, an important factor in the reconnection effects. The spacecraft data pool includes 1‐min average data by Themis (2007–2024), Cluster (2001–2022), and MMS‐1 (2015–2024) missions, as well as OMNI interplanetary data. The model drivers include the solar wind particle flux, IMF components, and the geodipole tilt angle. The model calculations faithfully reproduce the average plasma flow geometry and substantial effects have been found of the IMF orientation and magnitude, a principal factor that defines electromagnetic forces inside the magnetosheath. A strong dependence of the magnetosheath flow patterns on the Earth's dipole tilt indicates an important contribution of reconnection effects at the magnetopause to the solar wind particle transport around the dayside magnetosphere.

We address the sequence of Sun-to-Earth phenomena, that enables to study the mechanism for geoefficiency of eruptive prominences propagating from the Sun inside coronal mass ejections (CMEs). An eruptive prominence ejected in the solar wind (SW) moves at the SW velocity Earthward like adiamagnetic structure of eruptive prominence (DSEP).The key feature of the latter is a largesharp plasma concentration jump N inside the DSEP at a simultaneous sharp drop in the interplanetary magnetic field (IMF) modulus B. It is the anti-correlation between the N and B profiles in DSEP, due to which its contact with the magnetosphere may lead not only to magnetosphere compression, but also to penetration of DSEP substance into the magnetosphere. The duration of the magnetospheric disturbance (in the form of dayside auroras), global increase in the current systems, charged particle flux enhancement in the radiation belts, and generation of the irregular Pi2-3 oscillations aredetermined by the DSEP size. We present statistical investigations into DSEPs observed in different years of solar activity and builta qualitative modelfor DSEP geoefficiency.

The transient shear-induced particle migration of frictional non-Brownian suspensions is studied using particle-resolved simulations. The numerical method – the fictitious domain method – is well suited to heterogeneous flows thanks to a frame-invariant formulation of the subgrid (lubrication) corrections that does not involve the ambient flow (Orsi et al., J. Comput. Phys., vol. 474, 2023, 111823). The paper aims to give an accurate quantitative picture of the mass and momentum balance during the flow. The various assumptions and local constitutive laws that together form the suspension balance model (SBM) are thoroughly examined. To this purpose, the various quantities of interest are locally averaged in space and time, and their profile across the channel is extensively studied, with specific attention to the time evolution of the different contributions, either hydrodynamic in nature or from contact interactions, to the shear and normal stresses. The latter, together with the velocity gradient in the wall-normal direction and the volume fraction profile, yield the local constitutive laws, which are compared with their counterpart obtained in homogeneous shear flow. A fair agreement is observed except in a layering area at the boundaries and at the very centre of the channel. In addition, the main assumption of the SBM, i.e. the local relation between the hydrodynamic force on the particles and the particle flux, is meticulously investigated. The hydrodynamic force is found to be mainly a drag, except in the lower range of the probed volume fractions, where a non-drag contribution is observed.

Abstract The parameterization of suspended sediments in vegetated flows presents a significant challenge, yet it is crucial across various environmental and geophysical disciplines. This study focuses on modeling suspended sediment concentrations (SSC) in vegetated flows with a canopy density of avH ∈ [0.3, 1.0] by examining turbulent dispersive flux. While conventional studies disregard dispersive momentum flux for avH > 0.1, our findings reveal significant dispersive sediment flux for large particles with a diameter‐to‐Kolmogorov length ratio when dp/η > 0.1. Traditional Rouse alike approaches therefore must be revised to account for this effect. We introduce a hybrid methodology that combines physical modeling with machine learning to parameterize dispersive flux, guided by constraints from diffusive and settling fluxes, characterized using recent covariance and turbulent settling methods, respectively. The model predictions align well with reported SSC data, demonstrating the versatility of the model in parameterizing sediment‐vegetation interactions in turbulent flows.

Abstract The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths $\lamq$, input powers, and particle fuelling locations. Under H-mode scenario conditions with an upstream separatrix density of 1\e20 m$^{-3}$, the most conservative range of $\lamq$ extrapolations ($\approx$ 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power ($P_{SOL}$ = 29 MW) conditions, and 2/3 field with $P_{SOL}$ = 10 MW. Increasing the cross-field SOL diffusivities by 2-10x reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. 

A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fuelling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced.

The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

Three-dimensional (3D) equilibrium calculations, including the plasma rotation shielding effect to resonant magnetic perturbations (RMPs) produced by the island divertor (ID) coils, were carried out using the HINT and MARS-F codes on J-TEXT. Validation of 3D equilibrium calculations with experimental observations demonstrates that the shielding effect will prevent the penetration of the edge m/n = 3/1 mode component when the ID coil current is 4 kA, while change the size of magnetic islands once the current exceeds the penetration threshold. This indicates that equilibrium calculations including the plasma rotation shielding effect to RMPs can lead to better agreements with experimental observations compared to the vacuum approximation method. Additionally, the magnetic topology at the boundary undergoes changes, impacting the interaction between the plasma and the target plate. These results may be important in understanding RMP effects on edge transport and magnetohydrodynamic (MHD) instability control, as well as divertor heat and particle flux distribution control.