Behavior Of Particles Research Articles

The role of near-wall downdraught and asymmetric temperature distribution in dispersion of respiratory aerosols in radiant floor heating systems

In radiant floor heating (RFH) systems, the downdraught due to the asymmetric temperature distribution causes descending airflows towards the floor surface, where the warm air adjacent to the RFH system tends to be driven upwards by the buoyancy force. Hence, this conflict creates a zone that prevents the warm air to ascend affecting the streamlines of the indoor air. The present study aims to investigate the integrated effects of this phenomenon at different levels of the RFH temperature and thermal transmittance on the behaviour of indoor airborne particles (PMs2.5). In this context, a Eulerian-Lagrangian computational fluid dynamic (CFD) code is established, validated against experimental data, to address the dispersion and deposition patterns of PMs2.5. To generate different levels of the downdraught and non-uniform temperature distribution, five scenarios are considered regarding different thermal transmittance (U-value) levels, assessed for four RFH temperatures. Firstly, by introducing the near-wall and zonal spaces, the dispersion of particles in each scenario is evaluated. Then, the role of RFH system temperature in conjunction with each scenario is investigated. Finally, simple correlations are proposed allowing for fast evaluation of the decay rate coefficient of PMs2.5. According to the obtained results, the asymmetric temperature distribution causes a striking disparity in the zonal concentration of suspended PMs2.5, i.e., up to 32 % difference in the number of particles between quarters. It is shown that a faster decay rate of PMs2.5 is associated with a larger value of the characteristic temperature difference and the Rayleigh number (Ra). For a given RFH temperature, thermal performance improvement of the envelope reduces the number of deposited particles on the ceiling surface, whereas it gives a boost to the number of particles adhering to the floor. A sensitivity analysis on the results revealed that a 1 °C increment in the RFH temperature leads to an 8.4 % reduction on average in the number of suspended PMs2.5 in the breathing zone, regardless of the level of thermal transmittance from surrounding walls.

An Innovative Model of English Teaching Based on Particle Swarm Algorithm Construction

English teaching in construction plays a vital role in facilitating effective communication and collaboration within the industry. By equipping construction professionals with proficiency in English, they can navigate global markets, engage with multinational clients, and access a wealth of resources and knowledge available in the English language. English instruction tailored to the construction sector emphasizes industry-specific terminology, communication skills for project management, and comprehension of technical documents and specifications. Additionally, English language proficiency enables construction workers to comply with safety protocols, understand complex engineering plans, and participate in international conferences and training programs. This paper introduces an innovative model of English teaching tailored for construction professionals, utilizing a Particle Swarm Algorithm Construction with Weighted Particle Swarm Classification (W-PSO). Acknowledging the significance of effective communication in the construction industry, especially in multinational contexts, this research endeavors to optimize English language instruction through advanced computational techniques. The proposed model combines traditional teaching methodologies with the W-PSO algorithm, a metaheuristic optimization technique inspired by the behavior of particles in a swarm. W-PSO dynamically adjusts the weightings assigned to various instructional components, such as vocabulary acquisition, grammar instruction, and communication skills development, to optimize learning outcomes for construction professionals. Simulation results a cohort of 100 construction professionals participating in the English teaching program, the average improvement in English language proficiency is measured at 25% after completing the W-PSO-enhanced curriculum. Furthermore, specific language skills exhibit notable enhancements, with vocabulary acquisition increasing by 30%, grammar comprehension by 20%, and communication skills by 35%. The W-PSO algorithm dynamically adjusts the weightings assigned to various instructional components, optimizing the learning process. For instance, vocabulary acquisition receives a weighting of 0.4, grammar comprehension 0.3, and communication skills 0.3, resulting in a balanced and comprehensive approach to English language instruction.

The fractional soliton solutions of dynamical system arising in plasma physics: The comparative analysis

In light of fractional theory, this paper presents several new effective solitonic formulations for the Langmuir and ion sound wave equations. Prior to this study, no previous research has presented the comparision and obtained the generalized fractional soliton solutions of this kind with power law kernel and Mittag-Leffler kernel. The ion sound and Langmuir wave equations are essential in plasma physics, offering insights into the collective behavior of charged particles in plasmas and enabling diagnostics and control of these complex, ionized gas systems. The two distinct fractional order differential operators are substituted for the traditional order derivative to reshape the examined model. The Atangana-Baleanu non-singular and non-local operator and conformable fractional operator are the fractional-order operators that are used to create the fractional complex system equations for Langmuir waves and ion sound. A constructive approach new auxiliary equation method utilizes to obtain the exact analytical soliton solutions for ion sound and Langmuir wave equation. A wide range of soliton solutions is obtained, including mixed complex solitary shock solutions, singular solutions, mixed shock singular solutions, mixed trigonometric solutions, mixed singular solutions, exact solutions, mixed periodic solutions, and mixed hyperbolic solutions, dark soliton, bright soliton, trigonometric solutions, periodic results, and hyperbolic results. The solitons solution of the ion sound and Langmuir wave equations lies in their ability to maintain wave stability, their role in modeling wave propagation and nonlinear effects, their potential use as diagnostic tools, and their relevance in wave-particle interactions in plasma physics. The solitons provide a valuable framework for understanding the behavior of waves in plasmas and offer insights into the complex dynamics of these charged particle systems. A graphical comparison analysis of a few solutions is also shown here, taking into account appropriate parametric values through the use of the software package. Moreover, the results of this study have important implications for Hamilton's equations and generalized momentum, where solitons are employed in long-range interactions.

INFLUENCE OF SPRAY DISTANCE ON THE POROSITY OF Ni-BASED AMORPHOUS COATINGS: NUMERICAL SIMULATION AND EXPERIMENT

The influence of injection distance on the porosity of Ni-based amorphous coatings (AMCs) prepared with a high-velocity air fuel (HVAF) process is discussed based on a numerical analysis and experimental methods. A computational fluid dynamics model was established to demonstrate the gas flow field and behavior of particles in flight at different spraying distances during HVAF spraying. When analyzing the changes in the particle velocity and temperature, the spraying distance is less than 30 µm. The velocity and temperature changes of small particles have a significant impact, and the optimal spray distance (350 mm) for obtaining a low porosity coating is predicted. The calculation was validated experimentally by producing a Ni-based AMC with a low porosity (1.87 %) that was manufactured using the predicted HVAF optimal spraying distance.

Trajectory dynamics of particles in accelerated toroidal pipe: A computational study using CFD-DPM simulations

In this study, we investigate the behavior of solid particles of varying sizes within an accelerated toroidal filled with fluid using computational fluid dynamics- discrete phase model (CFD-DPM) simulations. A comprehensive grid independence study and validation were conducted to ensure the accuracy of our simulations. We focus on particles with 10-, 50-, and 100-microns diameters. Our simulations indicate that the primary flow (relative tangential velocity) of both the fluid and the particles follow the input acceleration. We also examined the secondary flow within the toroidal. Smaller particles (10 microns) closely follow the fluid flow with minimal lag, while larger particles (50 and 100 microns) exhibit a noticeable lag due to the inertia effect. Furthermore, the displacement from the initial position is also dependent on the particle size, with smaller particles traveling further distances. However, despite the ongoing acceleration of the fluid and the toroidal, the largest particles (100 microns) eventually settle and stop due to their high inertia (mass), and their velocity begins to decrease during the acceleration of the toroidal. These findings highlight the significant role of particle size in determining particle behavior within a toroidal, paving the way for future research and Particle Imaging Velocimetry Gyroscope (PIVG) application.

Influence of the internal structure of straight microchannels on inertial transport behavior of particles

The rapid advancement of Micro-Electro-Mechanical Systems (MEMS) technology has established microfluidics as a pivotal field. This technology marks the onset of a new era in various applications, including drug testing, cell culture, and disease monitoring, underscoring its extensive practicality and potential for future exploration. This research delves into the intricate behavior of particle inertial migration within microchannels, particularly focusing on the impact of different channel structures and Reynolds numbers (Re). Our studies reveal that particles in microchannels with one-sided sharp-cornered microstructures migrate towards the sharp corner at a relative position of 0.4 under low flow rates, and towards the straight wall side at a relative position of 0.8 under high flow rates. The migration pattern of equilibrium positions varies with different arrangements of sharp-corner structures, achieving stability at the channel's center only when the sharp corners are symmetrically arranged on both sides. Our investigation into the shape of microstructures indicates that sharp-cornered structures generate a more stable secondary flow compared to rectangular and semi-circular structures, preventing particle aggregation at the outlet. To address the challenges associated with handling variable cross-section geometries and solid-wall boundaries in dissipative particle dynamics methods effectively, we have developed a dissipative particle dynamics model specifically for analyzing such microchannels. Building upon our previous research, this model introduces a conservative force coefficient for particles within the microstructured region and an interaction zone that only involves repulsive forces, aligning well with experimental outcomes. Through the study of microstructures' geometric shapes, this paper offers guidance for designing microchannels for particle enrichment. Furthermore, the dissipative particle dynamics model established for the particle flow and solid structure interaction within microstructured channels provides insights into the mesoscale dynamics of liquid-solid two-phase flow and particle motion. In conclusion, this paper aims to enhance particle motion sample preparation techniques, thereby broadening the scope of microfluidic applications in the biomedical field.

Anomalous electrohydrodynamic cross-stream particle migration

Particle cross-stream migration in electrohydrodynamic microfluidic systems exhibits intriguing behaviors, which makes it interesting when viewed from a fundamental perspective and promising for nanoparticle focusing and separation applications. So far, particle behavior in such systems has been explained with the slip-induced lift force model (Saffman model), which predicts particle central or side focusing based on the direction of electric field and fluid flow. However, in our previous work, we observed particle migration patterns that did not adhere to the prediction of the Saffman model. In this work, we further studied this novel particle lateral migration behavior, which we termed the “anti-Saffman” behavior. We experimentally investigated how changing the conductivity of the suspending medium influences particle behavior and quantitatively measured the net lateral force experienced by the particles. Then, we compared this net force with the prediction of the relevant lift force models in the literature. We concluded that the anti-Saffman behavior is positively correlated with medium conductivity and shear rate (∝γ̇2). Furthermore, the comparison with the existing force models revealed that none of them can predict the experimentally observed particle lift. The net lift predicted by hydrodynamic lift models indicated that the underlying mechanism behind our experiments also potentially has a hydrodynamic origin. We believe this phenomenon offers the possibility of manipulating and separating nanoparticles suspended in standard aqueous electrolyte solutions, which makes it applicable to various biological samples.

Gas–solid flow mechanism of spherocylindrical particles with various aspect ratios in spouted bed

Spherocylindrical particles encompass elongated spherical particles, such as biomass particles utilized in thermal reactions and capsule pills in the pharmaceutical industry, along with various other particles found in industrial production processes. Understanding the flow behavior of spherocylindrical particles within a spouted bed holds significant practical importance. This study delves into investigating the fluidization characteristics of spherocylindrical particles with varying aspect ratios within a spouted bed. Key parameters studied include particle void fraction, velocity, orientation, and mixing. These were analyzed using the computational fluid dynamics–discrete element method, and the simulation model's accuracy was verified through high-speed photography experiments. The findings revealed distinct void fraction distributions at different locations within the bed. Particles in the higher-bed region exhibited larger void fractions with flatter peaks compared to those in the lower-bed region. Moreover, the correlation between particle aspect ratio and void fraction varied at different heights within the bed. Notably, there was no significant correlation observed between particle aspect ratio and particle velocity. However, the larger the particle aspect ratio/gas velocity, the better the particle mixing quality. Additionally, gas velocity and particle aspect ratio had little effect on particle orientation, but there was a tendency for particle orientation angles to hover around 45°–60° near the wall region.

Study of scalar particles through the Klein–Gordon equation under rainbow gravity effects in Bonnor–Melvin-Lambda space-time

In our investigation, we explore the quantum dynamics of charge-free scalar particles through the Klein–Gordon equation within the framework of rainbow gravity, considering the Bonnor–Melvin-Lambda (BML) space-time background. The BML solution is characterized by the magnetic field strength along the axis of the symmetry direction which is related to the cosmological constant Λ and the topological parameter α of the geometry. The behavior of charge-free scalar particles described by the Klein–Gordon equation is investigated, utilizing two sets of rainbow functions: (i) f(χ)=(eβχ−1)βχ , h(χ) = 1 and (ii) f(χ) = 1, h(χ)=1+βχ2 . Here 0<χ=∣E∣Ep≤1 with E representing the particle’s energy, E p is the Planck’s energy, and β is the rainbow parameter. We obtain the approximate analytical solutions for the scalar particles and conduct a thorough analysis of the obtained results. Afterwards, we study the quantum dynamics of quantum oscillator fields within this BML space-time, employing the Klein–Gordon oscillator. Here also, we choose the same sets of rainbow functions and obtain the approximate eigenvalue solution for the oscillator fields. Notably, we demonstrate that the relativistic approximate energy profiles of charge-free scalar particles and oscillator fields get influenced by the topology of the geometry and the cosmological constant. Furthermore, we show that the energy profiles of scalar particles receive modifications from the rainbow parameter and the quantum oscillator fields by both the rainbow parameter and the frequency of oscillation.

Fully resolved simulation of spherical and non-spherical particles in a turbulent channel flow

This paper investigates the behavior of finite-size particles in a turbulent channel flow using a custom direct numerical simulation solver within the FOAM-Extend framework. The solver integrates the cut-cell immersed boundary module with a Lagrangian particle-tracking subroutine capable of simulating the motion of both spherical and non-spherical particles. The study investigates the complex interactions between particles and turbulent structures, offering insights into how particle shape and orientation affect their behavior within the flow field. Additionally, it examines the collision dynamics of two spherical particles in a turbulent channel. The simulations reveal that particle shape significantly influences particle trajectories, rotation, and their interactions with turbulent structures. The Q-criterion visualization showed the creation of hairpin and vortex ring structures shed by the particles. Finally, the close proximity and collision of particles was shown to significantly modify the flow pattern and particle dynamics.

Eulerian- lagrangian dense discrete phase model (DDPM) of stenotic LAD coronary arteries in comparison with single phase modeling

In computational fluid dynamic studies related to blood flow, investigating the behavior of blood particles is crucial, especially red blood cells as they constitute a significant proportion of blood particles. Additionally, studying red blood cell movements is necessary, especially in stenotic artery geometries. A new multiphase scheme was utilized to demonstrate the effect of red blood cells on hemodynamics in complex coronary arteries and investigate the consequence of their motion. To investigate the effect of red blood cell movement on flow, the dense discrete phase model (DDPM) was used. This simulation was performed in 3D coronary arteries with different degrees of stenosis, utilizing blood pressure as inlet and outlet boundary conditions while assuming the arterial wall to be rigid. The model prediction shows good agreement with experimental data. Velocity values were comparable in both single-phase and two-phase flow simulations, but the shear stress in two-phase modeling had higher values. In the two-phase DDPM modeling, the recirculation areas indicated a higher probability of atherosclerosis plaque re-formation in the pre-stenosis area compared to the stenosis and post-stenosis areas. The DDPM model was found to be more effective in obtaining shear stress values in the artery. Additionally, this model provides good results compared to the single-phase model in investigating the movement of particles along the artery as well as recirculation areas that lead to the deposition of particles.

The transformation behaviour of ZnFe2O4 spinel particles in the CaO-FexO-SiO2-ZnO slag system

Treating zinc hydrometallurgical residues and other Zn-bearing solid waste in direct lead smelting processes has been recognized as the preferred solution for the sustainable development of the lead and zinc industry. However, new technical problems occurred when changing the starting material of lead-making processes, the pyrometallurgical process for treating zinc hydrometallurgical residues and zinc bearing solid waste usually had the CaO-FexO-SiO2-ZnO slag. Along with the increasing of Zn content in the raw material, high ZnO slag formed, and the ZnFe2O4 spinel particles intended to precipitate from the molten slag. The understanding of the particle behaviour of ZnFe2O4 spinel was the core part of monitoring the slag properties. This study investigated the transformation behaviour of zinc ferrite spinel particles and the changes in glassy slag phase through the high-temperature phase equilibrium experimental technique, and the temperature was 1100 °C-1300 °C and the oxygen potential was 10−8 atm – 10−4 atm. It is found that the particle size of ZnFe2O4 spinel in the molten slag increased along with the increasing of the oxygen potential. The Zn/Fe mass ratio in spinel particles decreased when increasing of temperature and it increased when increasing of oxygen potential. The Fe/SiO2 mass ratio and the Zn content in glassy slag phase gradually increased along with increasing temperature and they decreased along with increasing oxygen potential. The growing and decomposition of ZnFe2O4 spinel was essentially a process of material exchange with the glassy slag phase.

Pressureless sintering behavior of ZrC ceramic particles encapsulated by nano-sized carbon layer

The typical high-density zirconium carbide (ZrC) ceramic requires high temperature and high pressure, leading to high manufacturing cost for complex shapes. Herein, nano-sized carbon layer encapsulated ZrC particles were pressureless sintered to achieve near-net preparation of high-density and complex-shape ZrC ceramic. The results indicate that reaction between carbon layer and the ZrCxOy phase in ZrC was improved because the carbon layer shortened the oxygen diffusion pathway. The oxygen content decreased from 3.52 wt% to 0.82 wt% while the thickness of carbon layer decreased from 4.0 nm to 0.75 nm after sintering at 1950 °C for 1 h. Furthermore, the migration of oxygen from ZrCxOy phase to the carbon layer effectively facilitate the formation and growth of the sintering neck. The resulting ZrC ceramic exhibited a high relative density of 92.3 % and excellent mechanical properties. The Vickers hardness, indentation fracture toughness, micro hardness, elastic modulus, and flexural strength were 17.3 ± 0.7 GPa, 4.0 ± 0.5 MPa m1/2, 22.2 ± 2.4 GPa, 429.9 ± 20.7 GPa, and 324 ± 8 MPa, respectively. This work provides a novel perspective to fabricate near-net, high-density and complex-shape ZrC ceramic at lower temperature without pressure.

Design of closed loop for reproduction of gas-dominated hydrate flow.

The study presents a novel setup for measuring the flow regime of hydrate particles in a gas-dominated flow, which is of interest for applications such as natural gas transportation. A closed-flow loop, driven by a novel internal fan, enables continuous observation of hydrate particle behavior in a gas flow. The experimental setup allows the production and insertion of HFC134a gas hydrate particles with diameters of 10-50μm into the gas flow loop via a bypass loop. The performance curve of the internal fan is validated, and its suitability for achieving the required flow speed (5m/s) is demonstrated. Through an observation window using camera systems, the flow regime of glass beads is successfully visualized and analyzed. To validate the experimental data, a coupled computational fluid dynamics-discrete element method model is used to simulate the particle flow density distribution. The study findings demonstrate the effectiveness of the experimental setup in characterizing the flow regime of hydrate particles in a gas-dominated flow.

Analyzing sediment transport in dam-break-driven swash processes: Insights from laboratory experiments using a high-speed camera

Dam-break-driven experimental tests are conducted to explore the intricate dynamics of water and sediment movement during a swash process. Utilizing a combined approach of a high-speed camera and an acoustic doppler velocity profiler, detailed records of flow velocity and sediment movements are captured. This study offers novel insights regarding flow shear stress, layer thickness and velocity of sediment transport in the backwash process, as well as behaviors of individual sand particles. Linear relationships serve to depict alterations in bed shear stress over time. The changing rate of bed shear stress can be influenced by the factors such as the beach slope, sediment grain size, and location on the beach. Increasing layer thickness of sediment transport corresponds with increasing outer layer velocity, exhibiting a linear-exponential relationship. A tentative interpretation of the inflection point (of linear and exponential relationship) designates it as a division between sheet load and a conjunction of sheet load and suspended load. With the onset of a bore and under the influence of turbulence and bubbles, sediment predominantly moves upwards, leading to sediment pickup. However, during both the late uprush phase and the backwash phase, the primary direction of sediment movement reverses downwards, favoring sediment settling. Empirical equations are formulated to estimate the parameters on the bed shear stress and the sediment transport layer. This study enhances understanding of water and sediment movement in swash processes, supplying validation data for numerical models.

The Separation Behavior of TiB2 during Cl2-Free Degassing Treatment of 5083 Aluminum Melt

Utilizing titanium diboride (TiB2) inoculation for grain-refining purposes is a widely established practice in aluminum casthouses and foundries. Since this inoculation is usually implemented jointly with or between routine melt treatment steps ahead of casting, it is important to know whether and how other melt treatment processes affect the fade of TiB2 particles. For the present study, we investigated the influence of degassing process on the separation behavior of TiB2 particles in aluminum melt. Multiple sampling methods were employed and the samples were analyzed via spectrometer analysis. The removal efficiency of TiB2 during the gas-purging process of 5083 aluminum melt was confirmed to be significant over 10 min of treatment time. The rate at which the TiB2 content decays was found to increase with the impeller rotary speed from 400 rounds per minute (rpm) to 700 rpm. The separation rate of TiB2 particles was obtained to be 0.05–0.08 min−1 by fitting the experimental data. Particle mapping results suggest that the TiB2 particles were separated to a dross layer. The obtained experimental results were used to quantitatively evaluate the conventional deterministic flotation model. The deviation between the conventional model and the experimental data was explained through the entrainment–entrapment (EE) model. Suggestions were made for future analytical and experimental works which may validate the EE model.

Shape changes of a single hairy particle with mobile ligands at a liquid-liquid interface

We investigate rearrangements of a single hairy particle at a liquid-liquid interface using coarse-grained molecular dynamics simulations. We consider the particles with the same (symmetrical interactions) and different (asymmetrical interactions) affinity to the liquids. We show how ligand mobility affects the behavior of the hairy particle at the liquid-liquid interface. We found that such a hairy particle can take various shapes at the interface. For example, a Janus-like snowman consisting of a segment cluster and a bare part of the core, Saturn-like structures, and the core with a wide “plume” on one side. A configuration of the particle at the interface is characterized by the vertical displacement distance and the orientation of the particle relative to the phase boundary. The selected descriptors are used to characterize the shape of the segment cloud. We found that the shape of a particle and its localization at the interface can be determined by tuning the interactions with the liquids.

Dimerizing hard spherocylinders in porous media

This research focuses on the unique phase behavior of non-spherical patchy colloids in porous environments. Based on the theory of scaled particle (SPT), methods have been refined and applied to analyze the thermodynamic properties of non-spherical patchy particles in a disordered porous medium. Utilizing the associative theory of liquids in conjunction with SPT, we investigated the impact of associative interactions and connections between the functional nodes of particles on the formation of the nematic phase. Calculations of orientational and spatial distributions were conducted, which helped to understand the phase behavior of particles during the transition from isotropic to nematic phase under the spatial constraints imposed by the disordered matrix of the porous medium.

Fabrication and characterization of particle-stacking microporous nickel using laser powder bed fusion

Nickel (Ni) and nickel-based materials with specific pore configurations, particularly interconnected micropores as channels, have significant potential in new energy technologies. Laser 3D printing is a simple, rapid, and effective method for one-step formation; however, challenges remain in the fabrication of micropores (< 50 μm). This study successfully designed a particle-to-particle stacking microporous nickel structure using laser energy to regulate the melting and fusion behaviors of nickel particles. While laser power had a minimal effect on the pore size, it significantly impacted the porosity. The formation of interconnected pores depended on a large hatch distance and high laser scanning speed, while the latter also contributed to the creation of uniform micropores characterized by narrow sintered bands. The particle-stacking microporous nickel exhibited a homogeneous pore distribution, with pore sizes ranging from approximately 10 to 50 μm in both mean pore size and median pore size (D50), along with a high percentage size distribution in 1 – 10 μm. Furthermore, it displayed controllable porosities between 18.44% and 52.94%, along with remarkable resistance to deformation and superior compressive performance, featuring a high yield strength exceeding 50 MPa compared to the conventionally-processed porous nickels.

Modeling and analysis of particle behavior in fluidized bed bioreactors during non-Newtonian sewage treatment

The stability of biofilms attached to support particles plays a crucial role in fluidized bed bioreactors (FBBRs) when treating non-Newtonian sewage. Biofilm instability is governed by particle-particle (P-P) and particle-wall (P-W) collisions. This paper presents a numerical study of particle behavior in a FBBR using a combined approach of computational fluid dynamics (CFD) and discrete element method (DEM) tailored for non-Newtonian flow systems. Following validation, the CFD-DEM model is utilized to quantify the FBBR performance regarding bed expansion, void fraction distribution, and P-P/W interactions. The impacts of three variables are examined, including sewage superficial velocity and two rheological properties, i.e., sewage temperature and solid concentration. The results reveal that increasing sewage superficial velocity and sewage solid concentration or decreasing sewage temperature, leads to a transition in particle flow from symmetrical to asymmetrical in the radial direction, along with a circulating flow at the lower part of the bed. The relationships between the variables and the bed expansion, uniformity, and P-P/W interactions are also established, identifying optimum conditions to enhance bed uniformity and reduce P-P/W interactions. Overall, the findings demonstrate the value of the CFD-DEM model for understanding, designing, and optimizing FBBRs.