Particle Transport Research Articles

Numerical simulation of the movement of the water surface with macroscopic particles during a partial dam break for various complex reliefs

In this paper, it was considered a numerical simulation of the water surface movement during the partial collapse of the destruction of a dam with complex terrain. The numerical simulations took into account debris after a partial collapse of the dam, which imitates debris and moves downstream with the water flow. Carrying out these calculations brings the results closer to the real scenario. The mathematical model was based on the Navier–Stokes equations and uses the turbulent large eddy method (LES) model describing the flow of an incompressible viscous fluid. To describe the movement of a two-phase fluid, the volume of fluid (VOF) methods for the phase were used, and the Discrete Phase Model (DPM) and Macroscopic Particle Model (MPM) were used to describe the movement of particles. For the numerical solution of this system of equations, the Pressure Implicit Split Operator (PISO) numerical algorithm was chosen. The results show that the model is accurate and capable of handling very complex interactions such as particle transport or hydrodynamic actions on structures if the appropriate scales are reproduced. Also in this work, a three-dimensional (3D) model of the flow of a dam break on uneven terrain was considered and combined problems were performed that are closer to real conditions. For combined problems, flood zones and flood times were determined, knowledge of which will help evacuate people from dangerous areas. The accuracy of the 3D model and the selected numerical algorithm were verified using two natural measurements.

Rotational Taylor dispersion in linear flows

The coupling between advection and diffusion in position space can often lead to enhanced mass transport compared with diffusion without flow. An important framework used to characterize the long-time diffusive transport in position space is the generalized Taylor dispersion theory. In contrast, the dynamics and transport in orientation space remains less developed. In this work we develop a rotational Taylor dispersion theory that characterizes the long-time orientational transport of a spheroidal particle in linear flows that is constrained to rotate in the velocity-gradient plane. Similar to Taylor dispersion in position space, the orientational distribution of axisymmetric particles in linear flows at long times satisfies an effective advection–diffusion equation in orientation space. Using this framework, we then calculate the long-time average angular velocity and dispersion coefficient for both simple shear and extensional flows. Analytic expressions for the transport coefficients are derived in several asymptotic limits including nearly spherical particles, weak flow and strong flow. Our analysis shows that at long times the effective rotational dispersion is enhanced in simple shear and suppressed in extensional flow. The asymptotic solutions agree with full numerical solutions of the derived macrotransport equations and results from Brownian dynamics simulations. Our results show that the interplay between flow-induced rotations and Brownian diffusion can fundamentally change the long-time transport dynamics.

Integrating ytopt and libEnsemble to autotune OpenMC

Ytopt is a Python machine-learning-based autotuning software package developed within the ECP PROTEAS-TUNE project. The ytopt software adopts an asynchronous search framework that consists of sampling a small number of input parameter configurations and progressively fitting a surrogate model over the input-output space until exhausting the user-defined maximum number of evaluations or the wall-clock time. libEnsemble is a Python toolkit for coordinating workflows of asynchronous and dynamic ensembles of calculations across massively parallel resources developed within the ECP PETSc/TAO project. libEnsemble helps users take advantage of massively parallel resources to solve design, decision, and inference problems and expands the class of problems that can benefit from increased parallelism. In this paper we present our methodology and framework to integrate ytopt and libEnsemble to take advantage of massively parallel resources to accelerate the autotuning process. Specifically, we focus on using the proposed framework to autotune the ECP ExaSMR application OpenMC, an open source Monte Carlo particle transport code. OpenMC has seven tunable parameters some of which have large ranges such as the number of particles in-flight, which is in the range of 100,000 to 8 million, with its default setting of 1 million. Setting the proper combination of these parameter values to achieve the best performance is extremely time-consuming. Therefore, we apply the proposed framework to autotune the MPI/OpenMP offload version of OpenMC based on a user-defined metric such as the figure of merit (FoM) (particles/s) or energy efficiency energy-delay product (EDP) on Crusher at Oak Ridge Leadership Computing Facility. The experimental results show that we achieve the improvement up to 29.49% in FoM and up to 30.44% in EDP.

Statistical Survey of Interchange Events in the Jovian Magnetosphere Using Juno Observations

AbstractInterchange instability is known to drive fast radial transport of electrons and ions in Jupiter's inner and middle magnetosphere. In this study, we conduct a statistical survey to evaluate the properties of energetic particles and plasma waves during interchange events using Juno data from 2016 to 2023. We present representative examples of interchange events followed by a statistical analysis of the spatial distribution, duration and spatial extent. Our survey indicates that interchange instability is predominant at M‐shells from 6 to 26, peaking near 17 with an average duration of minutes and a corresponding M‐shell width of <∼0.05. During interchange events, the associated plasma waves, such as whistler‐mode, Z‐mode, and electron cyclotron harmonic waves exhibit a distinct preferential location. These findings provide valuable insights into particle transport and the source region of plasma waves in the Jovian magnetosphere, as well as in other magnetized planets within and beyond our solar system.

Simulation of special-shaped graded particulate hydraulic transport in deep-sea mining scenarios

Hydrodynamic transport is crucial in deep-sea mining engineering for conveying materials. This study investigates the impact of various particle shapes on slurry transport efficiency and flow characteristics at higher concentrations. A numerical method combining computational fluid dynamics (CFD) and the discrete element method (DEM) within the Euler-Lagrange framework is employed to simulate the flow and hydrodynamic properties of graded heterogeneous coarse particles in a deep-sea hydraulic lift pipe. The reliability of the simulations is verified through comparisons, focusing on feed concentration and particle gradation effects on dynamic characteristics and flow patterns. The study quantitatively analyzes concentration distributions under various flow regimes for different particle shapes and conveying concentrations in a vertical pipeline system. Flow characteristics under clogging conditions are described, considering flow transition and concentration distribution. Additionally, conditions for stable transportation are discussed concerning particle forces and transport reliability.

One-dimensional model for vertical hydraulic transport of high-concentration mineral particles

A novel model is proposed for analyzing high-concentration granular flow systems comprising equally sized spherical particles within vertical, long straight pipelines. This model is specifically tailored for simulating the vertical hydraulic transport of ore particles in marine mining projects. The proposed model treats the granular system akin to a pseudo-fluid and operates through three mechanisms. First, fluid characteristics of the granular system are derived from particle–particle collisions. Second, the resistance exerted by the pipe wall on the granular system is calculated based on the momentum loss of particles during particle–wall collisions. Third, the interaction between individual particles and the surrounding fluid is transformed into an interaction between the carrier fluid and the pseudo-fluid. Additionally, the present work develops a dedicated numerical format and iterative method for solving the one-dimensional two-fluid governing equations. The one-dimensional (1D) model notably enhances computational efficiency and facilitates accurate tracking of high-concentration particles over extended distances within straight pipelines. Notably, the proposed 1D model demonstrates a high degree of predictive accuracy when compared against experimental data as well as results from computational fluid dynamics and discrete element method simulations.

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of E→×B→ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the −E→×B→ direction, where E→ and B→ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

Particulate transport in porous media at pore-scale. Part 1: Unresolved-resolved four-way coupling CFD-DEM

Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) is a powerful approach to simulate particulate flow in porous media at the pore-scale, and hence decipher the complex interplay between particle transport and retention. Two separate CFD-DEM approaches are commonly used in the literature: the unresolved (particle smaller than the grid cell size) and the resolved (particle bigger than the grid cell size) approach. In this paper, we propose a novel CFD-DEM coupling approach that combines both unresolved and resolved coupling. Our new modeling technique allows for the simulation of particulate flows in complex pore morphology characteristic of porous materials. It relies on an efficient searching strategy to find grid cells covered by the particles and on an appropriate calculation of the fluid-solid momentum exchange term. The robustness and efficiency of the computational model are demonstrated using cases for which reference solutions – analytical or experimental – exist. The new unresolved-resolved four-way coupling CFD-DEM is used to investigate pore-clogging and permeability reduction due to the sieving and bridging of particles.

Flow topology and mixing in alveolar edema: Unsteady flow in interconnected cavities with moving walls

The study of alveolar fluid mechanics is critical for comprehending respiratory function and lung diseases, particularly in cases of alveolar lesions that result in significant structural and fluid dynamic changes. This study investigates the flow topology and chaotic mixing within both normal and edematous alveoli, where the alveoli in the edematous model are interconnected by pores. To numerically simulate alveolar flow, a mathematical model is developed to ascertain the key parameters of Reynolds number (Re) and alveolar expansion ratio. Subsequently, the flow fields are analyzed to determine wall shear stress (WSS) and to identify WSS critical points and critical points of velocity vector, with a thorough presentation of the various flow topologies corresponding to these critical points. Moreover, a dynamic mode decomposition-based method is introduced to track particle trajectories, and the exploration of chaotic mixing is conducted through tracer advection, Poincare map, and the calculation of finite-time Lyapunov exponents. Results indicate that the edematous model exhibits a higher Re and higher WSS due to the fluid properties. Within the alveoli, high WSS is usually localized at the pores. The pores increase critical points and alter flow topologies, significantly changing chaotic mixing. Additionally, Re and alveolar locations also affect mixing patterns. These findings are crucial for understanding alveolar physiology and designing inhaled drugs for lung diseases, considering the role of chaos in particle transport in the lung acini.

Turbulence characteristics of non-spherical isolated particles with different flatness on a rough bed surface

The transport of non-spherical particles on bed is significantly influenced by flow dynamics and particle shape. To study the impact of changes in the flatness of isolated bed particles on the surrounding flow characteristics, we prefabricated three isovolumetric particles with varying levels of flatness and analyzed the flow field around them using two-dimensional particle image velocimetry measurement methods. The results show that the presence of the particles significantly alters the turbulence parameters of the surrounding flow. Specifically, the influence extends up to 1D upstream of the particles, 0.5D above the particle tops, and throughout most of the downstream region, where D represents the equivalent particle diameter. The particles cause a reduction in Reynolds stress and turbulent kinetic energy in the near-bed upstream region, with a more pronounced effect observed for particles with lower flatness. Increased particle flatness is associated with a smaller difference in the longitudinal time-averaged velocity between upstream and downstream regions. In the wake region, particle flatness shows a negative correlation with Reynolds stress, longitudinal turbulence intensity, and turbulent kinetic energy, with this effect being more pronounced closer to the particles. Quadrant analysis indicates that ejection (Q2) and sweep (Q4) quadrant events remain dominant in the presence of particles, while particle flatness positively correlates with the frequency of outward (Q1) and inward (Q3) events in the wake region and negatively correlates with the frequency of Q2 and Q4 events. This study enhances our understanding of particle-fluid interactions, particularly the response relationship between non-spherical particles and turbulent mechanisms.

The impact of Microbially Induced Calcite Precipitation (MICP) on sand internal erosion resistance: A microfluidic study

Fine particle internal erosion involves the detachment, transport, and deposition of fine particles within soil, significantly impacting agriculture, engineering, and environmental protection. Microbially induced calcite precipitation (MICP) has proven to be an effective method for controlling internal erosion. To optimize MICP protocols for effective erosion control, understanding the microscopic mechanisms by which MICP reduces fine particle erosion is essential but remains unclear. In this study, microfluidic chip experiments were conducted to observe the characteristics of calcium carbonate crystals and fine particles before and after MICP reinforcement and erosion. Through quantitative analysis of the calcium carbonate produced by MICP and eroded fine particles, the effects of bacterial density, concentration of cementation solution, and flow rate on the efficiency of MICP in resisting soil internal erosion were investigated. In addition, three primary mechanisms by which MICP reduces fine particle erosion were identified. Firstly, in situ sand stabilization occurs when calcium carbonate generated by MICP bonds and encapsulates fine particles, forming larger particles that remain stable under erosive flow conditions. Secondly, regional sand stabilization is achieved as MICP-produced calcium carbonate crystals narrow or block the flow channels, preventing extensive migration of fine particles. Lastly, adjacent particle stabilization is facilitated by calcium carbonate crystals, which remain stable during water flow erosion and alter the erosion flow lines, creating zones adjacent to the crystals where fine particle movement is minimized. These findings provide a deeper understanding of the role of MICP in erosion control and can inform the development of more effective erosion mitigation strategies.

Multispacecraft Observations of Protons and Helium Nuclei in Some Solar Energetic Particle Events toward the Maximum of Cycle 25

The intricate behavior of particle acceleration and transport mechanisms complicates the overall efforts in formulating a comprehensive understanding of solar energetic particle (SEP) events; these efforts include observations of low-energy particles (from tens of keV to hundreds of MeV) by space-borne instruments and measurements by the ground-based neutron monitors of the secondary particles generated in the Earth atmosphere by SEPs in the GeV range. Numerous space-borne missions provided good data on the nature/characteristics of these solar particles in past solar cycles, but more recently—concurrently with the rise toward the maximum of solar cycle 25—the High-Energy Particle Detector (HEPD-01) proved to be well suited for the study of solar physics and space weather. Its nominal 30–300 MeV energy range for protons can enlarge the detection capabilities of solar particles at low Earth orbit, closer to the injection limit of many SEP events. In this work, we characterize three SEP events within the first six months of 2022 through spectral and velocity dispersion analysis, assessing the response of HEPD-01 to >M1 events.

Analytical calculation of the kinetic q factor and resonant response of toroidally confined plasmas

Symmetry-breaking perturbations in axisymmetric toroidal plasma configurations have a drastic impact on particle, energy, and momentum transport in fusion devices, thereby affecting their confinement properties. The perturbative modes strongly affect particles with specific kinetic characteristics through resonant mode–particle interactions. In this work, we present an analytical calculation of the kinetic q factor, enabling the identification of particles with kinetic properties that meet the resonant conditions. This allows us to predict the locations and structures of the corresponding resonant island chains, as well as the existence of transport barriers in the particle phase space. The analytical results, derived for the case of a large aspect ratio configuration, are systematically compared to numerical simulations, and their domain of validity is thoroughly investigated and explained. Our findings demonstrate that calculating the kinetic q factor and its dependence on both particle and magnetic field characteristics provides a valuable tool for understanding and predicting the resonant plasma response to non-axisymmetric perturbations. Moreover, this approach can be semi-analytically applied to generic realistic experimental equilibria, offering a low-computational-cost method for scenario investigations under various multi-scale perturbative modes.

Compression Acceleration of Protons and Heavier Ions at the Heliospheric Current Sheet

Recent observations by the Parker Solar Probe (PSP) suggest that protons and heavier ions are accelerated to high energies by magnetic reconnection at the heliospheric current sheet (HCS). By solving the energetic particle transport equation in large-scale MHD simulations, we study the compression acceleration of protons and heavier ions in the reconnecting HCS. We find that the acceleration of multispecies ions results in nonthermal power-law distributions with a spectral index consistent with the PSP observations. Our study shows that the high-energy cutoff of protons can reach Emax∼0.1 –1 MeV depending on the particle diffusion coefficients. We also study how the high-energy cutoff of different ion species scales with the charge-to-mass ratio Emax∝(Q/M)α . When determining the diffusion coefficients from the quasi-linear theory with a Kolmogorov magnetic power spectrum, we find that α ∼ 0.4, which is somewhat smaller than α ∼ 0.7 observed by PSP.

Distribution of transitional flow deposits in sedimentary environments of mixed sand-mud turbidite system

Predicting facies distribution in turbidite systems is essential for resource exploration and identifying geohazards from an economic standpoint. Models that describe facies distribution depend heavily on understanding the mechanisms of particle transport and deposition. These processes are closely tied to the volume, concentration, and composition of sediment gravity flows, which display a range of behaviours between turbulent and laminar flow extremes. Recently, there has been a rise in studies on transitional flow deposits, although they remain much less understood than fully turbulent or laminar flows.For the first time, the distribution of Structured sandstone–mudstone associated with transitional flow deposits has been quantitatively presented for various sedimentary environments within the turbidite system. The distribution of Structured sandstone–mudstone was analysed for six areas of the Ropianka Fm (Skole Nappe, Polish Outer Carpathians) across twelve sedimentary environments, including channels, channel-levees, channel-mouths, and sub-environments of the depositional lobe. An increased amount of Structured sandstone–mudstone was observed in proximal settings away from the transport axis and in the distal parts of the turbidite system. It was found that flow transformation can occur in both proximal and distal zones of the turbidite system. Structured sandstone–mudstone in proximal zones is more often deposited from diluted mud-laden flows of small volume, where fine-grained cohesive material likely underwent vertical segregation. In contrast, Structured sandstone–mudstone formed basinward tend to be initially formed by larger flows. In proximal part of depositional lobe setting, the flow transformation due to longitudinal or longitudinal and vertical segregation of fine-grained cohesive material occurs slowly. At this point the velocity of sediment gravity flow is too high and the concentration of cohesive particles is too low for common development of cohesive bonds. Flow transformation accelerates in lobe fringe and lobe distal fringe/interlobe, due to flow deceleration, changes in sand-to-mud ratio, and the time required for development of cohesive bonds and the transition to a transitional flow regime, leading to increased deposition of Structured sandstone–mudstone distally.

Retention of microplastics by biofilms and their ingestion by protists in rivers.

Microplastics (MPs) are released into the environment through human activities and are transported by rivers from land to sea. Biofilms, which are ubiquitous in aquatic ecosystems such as rivers, may play an essential role in the fate of MPs and their ingestion by biofilm protists. To assess this, biofilms were naturally grown on clay tiles in the River Rhine, Germany, and analysed in a combined field and laboratory study. Compared to the ambient river water, biofilms grown for 6, 12, and 18 months in the River Rhine contained up to 10 times more MPs. Between 70% and 78% of all MPs were smaller than 50 μm. In laboratory experiments, clay tiles covered with 1-month-old naturally grown biofilm retained 6-12 times more MPs than clay tiles without biofilm coverage. Furthermore, the ingestion of MPs of 6 and 10 μm by the ciliate Stentor coeruleus was confirmed, and a positive correlation between ingestion rates and ambient MP concentrations was found. The results are relevant for particle transport models in riverine systems, risk assessment of MPs regarding their distribution and fate in the aquatic environment, and the effects of MPs on micro- and macroorganisms.

Full-f gyrokinetic simulations of hydrogen isotope mixing in tokamak plasmas

Hydrogen isotope mixing phenomena in tokamak plasmas are analyzed using global full-f gyrokinetic simulations. Model plasma parameters are chosen based on hydrogen isotope pellet experiments on JET, in which fast hydrogen isotope mixing in the timescale of the energy confinement time τE occurred after injecting deuterium (D) pellets into hydrogen (H) plasmas. Two numerical experiments are conducted using plasma profiles before and after the pellet injection. In both cases, turbulent fluctuations in the plasma core are characterized by ion temperature gradient driven turbulence. In the former case, the density profile of bulk H ions is kept in a quasi-steady state, and the particle confinement time of bulk H ions τH is an order of magnitude longer than τE. In the latter case, the density profiles of bulk H ions and pellet D ions show transient relaxation in the timescale of τE, and both τH and the timescale of particle pinch τD become comparable to τE, indicating the fast hydrogen isotope mixing. In the toroidal angular momentum balance, it is found that the enhanced ion particle transport leading to the fast hydrogen isotope mixing is driven by the toroidal field stress.

How is particulate organic carbon transported through the river-fed submarine Congo Canyon to the deep sea?

Abstract. The transfer of carbon from land to the near-coastal ocean is increasingly being recognized in global carbon budgets. However, a more direct transfer of terrestrial organic carbon to the deep sea is comparatively overlooked. Among systems that connect coastal to deep-sea environments, the submarine Congo Canyon is of particular interest since the canyon head starts 30 km into the Congo River estuary, which delivers ∼7 % of the dissolved and particulate organic carbon from the world's rivers. However, sediment and particulate organic carbon transport mechanisms that operate in the Congo Canyon and submarine canyons more globally are poorly constrained compared to rivers because monitoring of deep-sea canyons remains challenging. Using a novel array of acoustic instruments, sediment traps, and cores, this study seeks to understand the hydrodynamic processes that control delivery of particulate organic carbon via the submarine Congo Canyon to the deep sea. We show that particulate organic carbon transport in the canyon axis is modulated by two processes. First, we observe periods where the canyon dynamics are dominated by tides, which induce a background oscillatory flow (speeds of up to 0.15 m s−1) through the water column, keeping muds in suspension, with a net upslope transport direction. Second, fast-moving (up to 8 m s−1) turbidity currents occur for 35 % of the time during monitoring periods and transport particulate organic carbon with mud and sand at an estimated transit flux that is more than 3 to 6 times the flux induced by tides. Organic carbon transported and deposited in the submarine canyon has a similar isotopic composition to organic carbon in the Congo River and in the deep-sea fan at 5 km of water depth. Episodic turbidity currents thus promote efficient transfer of river-derived particulate organic carbon in the Congo submarine fan, leading to some of the highest terrestrial carbon preservation rates observed in marine sediments globally.

Soil Erosion Analysis with GIS: A Case Study of Oued Sly Catchment, Northern Algeria

Soil erosion is a process that involves the detachment and transportation of soil particles by water, presenting notable challenges of the depletion of soil fertility, diminished agricultural productivity, and consequences for water reservoirs. The integration of Geographic Information System (GIS) and modeling techniques provides a valuable platform for assessing and addressing soil water erosion. This study specifically examines the risk of soil water erosion in the Oued Sly catchment within the Middle Cheliff watershed in northern Algeria. Employing the Universal Soil Loss Equation (USLE) and thematic maps, the research evaluates the key factors influencing erosion. Results reveal that the catchment generally experiences low to moderate erosion levels, averaging 23.4 to 40 t ha-1 yr-1 across 90.7% of the total area. However, around 9.3% of the catchment exhibits higher erosion rates exceeding 40 t ha-1 yr-1. These findings offer valuable insights for informed watershed management and strategies to mitigate soil water erosion in comparable regions

A Study of the Cereal Benches Based System for Protection of Soil in Algerian Oued Mellegue Tebessa Dry Lands

The study was intended to mitigate the destructive after-effects of erosion, where runoff water erodes the arable layers through the increased transport of soil particles. The study was conducted in Ain Zerga across 40 plots, covering an area of 405 hectares, during 2023 and 2024, with a focus on the dimensional evaluation of cereal beds. The evaluation mainly focused (a) on the durability of the structure, including the measurement of settlement, which indicates distribution and strength of resistance to climatic hazards, and (b) proportional sustainability and the interbank aspects. The height of benches, at 0.481± 0.086 m, and the length of 82.583 ± 24.203 m were influenced by height differences and the intensity of water flows. The retention of runoff water was influenced by the slope toward the natural spillway. Furthermore, the destructive intensity of the runoff, amplified by the sharp increase in altitude, was confirmed by statistically significant differences in the length, height, and number of cracks between the study sites (p-value = 0.000). The number of cracks per breach at the Gastel site (4.25 units) and Kef Mebdoua (3.917 units) was higher compared to Ain Zerga, Berket Frass, Ouled Mebarek, and Zlassde, where the number ranged from 0.750 to 1.583 units. The statistical significance in the number of cracks (p<0.001) signified that variation in cracks between the locations was due to the topographic steep slopes on site. The non-significant difference in width suggests that the observed differences could be due to the phenomenon of compaction following the annual infiltration of water on the ridges of cereal bench, plus the act of trampling constantly exercised by sheep and cattle herds when passing to the paths. This work proposed in this article is an attempt to structure a chain of impact analysis of risks and vulnerability of intensity via the sustainability of cereal bench carried out in the Tebessa region.