Particle Flux Research Articles

Implications of gravel content on aerodynamic parameters, sand flux, erosion and accumulation during deflation processes over Gobi

The deflationary process of the Gobi Desert contributes to the release of atmospheric dust and may lead to sand hazards. This process is influenced by the presence of gravel on the desert surface. While extensive research has been conducted on the impact of gravel coverage on sand transport, there is lack of studies examining the wind-sand characteristics across varying gravel content in mixed beds of sand and gravel. In a series of wind tunnel experiments, aerodynamic parameters, sand flux, and both erosion and accumulation amount were measured across beds with differing gravel contents. The findings revealed that the maximum roughness and friction velocity were observed at a 35 % gravel content. As gravel content increased, both the sand flux density (q) and the total sand flux within 30 cm (qt30) declined, whereas the saltation height (h50) increased. Both the deflation rate and the sediment entrainment rate showed a decrease as the gravel content increased from 15 % to 55 %. Although particle flux density exhibited significant fluctuations, no clear correlation with gravel content was found. The deflation process is accompanied by distinct stripes of sand accumulation at 15 % and 25 % gravel contents, owing to the saturated sand flow characterized by a high sand entrainment rate and particle flux density. In contrast, a gravel content of 55 % exhibited stability with only minor deflation, corresponding to a low sand entrainment rate and particle flux density, thus approaching an equilibrium state. This observation indicates that increasing gravel content effectively slows the deflation rate.

Spectrometers of Neutrons and Fast Atoms of Tokamak Thermonuclear Plasma Based on CVD Synthesized Diamond Single-Crystal Films

High radiation resistance, chemical inertness, the ability to operate at elevated temperatures, high mobility and efficiency of charge-carrier collection are important properties of diamond for designing detectors and spectrometers of ionizing radiation. Currently, diagnostics of neutrons and neutral particle fluxes based on diamond detectors for the ITER thermonuclear reactor are justified and developed. This work presents the results of a Raman spectroscopy and photoluminescence spectroscopy study of the electronic quality of synthesized epitaxial diamond films obtained by vapor deposition in a hydrogen and methane mixture in the ARDIS reactor on boron-doped single-crystal diamond substrates. To confirm their electronic quality, detectors have been made from films selected by spectrometric methods and the charge collection efficiency and energy resolution have been measured when irradiated with alpha particles from a 241Am source and 14.7 MeV fast neutrons from the ING-07T2 neutron generator.

Determination of threshold values of parameters of electronic irradiation of glass leading to electrostatic discharges

Experimental data are presented on the minimum values of energies and flux densities of electrons, the impact of which on the cover glasses of solar batteries and reflecting elements of thermoradiators of artificial Earth satellites leads to electrostatic discharges. It has been established that the addition of protons to the composition of the particle flux acting on the studied samples can suppress the development of discharges. For a qualitative interpretation of the results obtained, a mathematical model is proposed.

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

Dependence of divertor turbulence on plasma density and current in TCV

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.

Ultra-relativistic electron flux enhancement under persistent high speed solar wind stream

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.

Experimental investigation on high heat flux plasma parameters of HIT-PSI device in argon discharges

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.

Novel Bragg peak characterization method using proton flux measurements on plastic scintillators

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

Laser-accelerated proton beams (LAP) have inspired innovative applications that can leverage proton bunch properties distinct from those of conventionally accelerated proton beams. A new multifunctional LAP beamline has been proposed, utilizing the GEANT4 toolkit, and the solenoids have been simulated using high-precision components.We present the design of the initial segment of the transport beamline and energy selection system for the effective transport of a radiation pressure acceleration source with a high energy spread. This study examines the effects of solenoids on proton beam quality, specifically focusing on profiles, FWHM, and transmission efficiency. The results obtained using our developed GEANT4 toolkit validate the analytical formulas and findings from the DYNAMION code regarding emittance and proton profiles. Additionally, the presented GEANT4 LAP beamline is capable of calculating the flux of secondary particles, such as neutrons and photons. It was observed that careful attention to the precise structure of the solenoids is critical.

Effect of High Heat Flux of Helium and Hydrogen Plasma Jet on the Material Properties of Piezoelectric PZT-Ceramics

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.

State primary standard for units of absorbed dose and absorbed dose rate of photon, electron, proton radiation and in carbon ion beams, quantity, fuence, fux density and energy of particles in proton beams and heavy charged particles GET 38-2024

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.

Sediment transport on rippled beds

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.

Dependence of divertor asymmetries on the toroidal magnetic field

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.

Numerical assessment of ICRF-specific plasma-wall interaction in the new ITER baseline using the SSWICH-SW code

In 2023, switching the material on the first wall of ITER to tungsten (W) was recommended. In magnetic Fusion devices, waves in the Ion Cyclotron Range of Frequencies (ICRF) interact with the Scrape-Off Layer (SOL) via RF-sheath rectification. This contribution re-assesses this phenomenon close to the ITER ICRF antenna, focusing on the ICRF-specific gross erosion of W from the antenna port sides. Our quantitative estimates rely on predictive multi-2D numerical simulations of the ICRF antenna environment using the SSWICH-SW code. They combine Slow Wave propagation from the antenna mouth to the SOL, the excitation of RF oscillations in the sheath voltages at the antenna port sides and a subsequent DC biasing of the SOL. Maps of the parallel RF electric field at the antenna mouth, from the antenna code TOPICA, excite the system. Our simulations cover more than four decades in the local densities near the antenna. Since both the sputtering and the local heat loads are proportional to the local particle fluxes, the most intense Plasma-Wall Interaction is found for high local density, with or without ICRF waves. In these conditions, larger margins also exist for coupling the ICRF power. We tested several operational trade-offs between these two constraints. The simulated target plasma contains 2% of neon ions. These are efficient at sputtering W, already at low accelerating voltages. Consequently, although the RF-sheath rectification sufficiently amplifies the local sputtering at the antenna port for a detection using visible spectroscopy, the ICRF-induced increment of the gross W production represents at worse 22% of the W source expected from thermal sheaths over the eighteen out-board mid-plane ports. An upper bound, independent of our main assumptions, is proposed for this enhancement factor. This moderate expected global increase questions the ability to detect ICRF-specific W contamination of the plasma core, even at the planned maximal ICRF power.

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

Mechanisms of formation and losses of radiation belts are the most important questions of magnetospheric physics, especially in a subsecond temporal scale. Energetic particles release their energy in the atmosphere producing fluorescent emission in characteristic wavelength bands. This emission is measurable and can be an additional information source on the spatiotemporal structure of particle fluxes and spectra. Here we present the world’s first measurements of UV-microbursts during geomagnetic disturbance and pulsating aurora caused by high-energy electron precipitation. It demonstrates that fundamental questions of magnetospheric physics mentioned above can be addressed by using the optical measurements by highly sensitive photometers with high temporal resolution. Such a photometer was installed at Verkhnetulomsky observatory at Kola Peninsula and measured a series of short (less than 0.5 s) pulses of emission with an angular size of bright spot ∼0.2 rad. Simultaneous measurements of high-energy electron fluxes made by the NOAA-19 satellite and fine temporal structure of geomagnetic pulsations demonstrate a magnetospheric origin of the observed events.

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

Tungsten (W) is one of the preferred plasma-facing materials for future fusion reactors because of the advantages of high melting point, low physical sputtering rate and low fuel retention. The lifetime and service performance of W materials are limited by erosion processes. In recent years, tungsten deuteride molecular (WD) sputtering has been discovered in tokamak. Studying WD sputtering helps to better understand the erosion of tungsten materials in tokamak. EAST adopts the method of injecting impurity gas in the divertor to alleviate the particle flux and heat flux bombarding the tungsten divertor, and CD4 is one of the impurities used. CD4 injection experiments on EAST show that after CD4 is injected from the upper outer (UO) divertor target, despite that the electron temperature (Te) around the strike point on the UO divertor target drops to about 8 eV, strong WD sputtering occurs in the far scrape-off layer (SOL) region. The divertor probe data shows that after CD4 injection, the particle flux from the upstream plasma onto the divertor split into two bands, creating a new particle flux strike point in the far SOL region. The energy of incident particles at the original strike point is significantly reduced, while the energy of incident particles at the far SOL region is enhanced. The edge electron density profile at mid-plane indicates that CD4 injection causes an increase of density in the SOL, which could improve the coupling of Lower-Hybrid Wave (LHW) with the edge plasma as evidenced by the decreasing reflected power of LHW. The increase in absorbed power of LHW could change the particle and energy transport in the edge plasma, thus plays an important role in the enhanced WD sputtering in the far SOL region after CD4 injection.

Sedimentary organic matter signature hints at the phytoplankton-driven biological carbon pump in the central Arabian Sea

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.

Magnetosheath Plasma Flow and Its Response to IMF and Geodipole Tilt as Obtained From the Data‐Based Modeling

AbstractLarge‐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.

Magnetospheric response on impact of solar wind diamagnetic structures borne by eruptive prominence

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.

Mass and momentum balance during particle migration in the pressure-driven flow of frictional non-Brownian suspensions

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.