Particle Collision Rate Research Articles

Ideal Collision Rate of Symmetric Hexapods in a Simple Shear Flow.

An ideal collision rate (ICR) is defined as the average rate of contact between two particles that translate and rotate with the imposed fluid flow in the absence of interparticle interactions. ICR in a simple shear flow provides an estimate of the collision rate in a flowing dilute particle suspension, and its value is known only for traditional convex shapes such as spheres and cylinders. In this work, we compute the ICR for a family of symmetric hexapods (shown in Figure 1) that are particles composed of three orthogonal cylinders with coinciding centers (forming six arms) with at least two cylinders having the same length. We employed the finite-element method to obtain the rotational dynamics of hexapods and Monte Carlo simulations to calculate their ICR. Our results indicate that the ICR for hexapods is not directly proportional to the volume of the particle, in contrast to what is observed for convex shapes such as spheres and cylinders. For hexapods of the same size, the ICR could be similar for particles with order-of-magnitude differences in their volumes and it could also vary by an order of magnitude for particles with similar volumes. Specifically for hexapods termed branched fibers, which have one longer cylinder, the ICR changes by an order of magnitude even when the shorter cylinders are significantly smaller than the longer cylinder. This change is attributed to the increase in the tumbling frequency of the particle due to hydrodynamic forces acting on the shorter arms, in addition to the increased probability of collision afforded by them. Using asymptotic theory for high-aspect ratio cylinders, we showed that the tumbling period of hexapods was proportional to an algebraic power of the ratio of its arm lengths and has a weaker logarithm relationship with the aspect ratio of its arms. The collision cross-section provided the relative cross-streamline particle separations for the most likely binary collisions in the suspensions, and its value was sensitive to the arm lengths, as well. The collision rate of hexapods also could not be estimated within an order of magnitude from existing geometric models, such as the ICR of a circumscribing convex shape, ICR for one of the individual cylinders of the hexapod, or using the particle volume times the shear rate. Our work indicates that collision rate for non-convex particles in shear flows critically depends on the shape of the particle due to nontrivial changes in the particle's orientational dynamics, and the ICR calculation serves as a more reliable method for estimating their true collision rates in the suspension.

Settling and collision of spheroidal particles with an offset mass centre in a quiescent fluid

Gravitational sedimentation and the collision of particles play key roles in various natural and engineering processes. In practice, particles are often non-spherical in shape with non-uniform mass distribution. In this study we investigate how the mass eccentricity influences the settling and gravitational collision of non-spherical particles in a quiescent fluid. Firstly, we theoretically analyse the effect of mass centre offset on the settling motion of a single spheroid under the low-Reynolds-number assumption. We find that the competition of fluid-inertia torque and gravitational torque determines the terminal settling mode of the spheroid. As the mass centre offset increases from zero to a critical value, the orientation of a settling spheroid undergoes a transition from the broad-side-on to narrow-side-on alignment. With an intermediate mass centre offset, the settling spheroid prefers an oblique orientation with a horizontal drift. Secondly, we investigate the gravitational collision rate of settling spheroids. With the change of particle orientation, the collision kernel exhibits a non-monotonic variation with a maximum when particles settle with an intermediate oblique orientation. Therefore, adjusting the mass centre offset to alter particle orientation can indirectly affect the collision rate of settling spheroidal particles. In summary, our findings reveal the significance of the mass eccentricity on particle dynamics in fluid flows, and suggest a potential approach for manipulating the settling motion and collision rate of non-spherical particles by adjusting their mass centre position.

Optimizing Metal‐Surface Water Disinfection: CFD Study on Microorganism Collision Against a Triply Periodic Minimal Surface

AbstractThis research presents a point‐of‐use (POU) water treatment technology utilizing metal/metallic surfaces based on Triply Periodic Minimal Surface (TPMS) structured filtration infills. Computational Fluid Dynamics modeling using Ansys CFX software is employed to analyze the behavior of Escherichia coli bacteria within a continuous liquid phase moving through the filtration infill, assess particle collision dynamics, and evaluate the efficiency of filtration, considering pressure drop as a fundamental factor in process energy consumption. TPMS infill meshes are coded in a Schwarz P shape using Python, and the Darcy‐Forchheimer equation is employed to determine the permeability and resistance loss coefficient of the infill geometry. The results indicate that the TPMS infill efficiently captures particles while introducing a negligible pressure drop into the system. It is found that a single 40 mm infill configuration is the most efficient, exhibiting a higher collision rate compared to the smaller 20 mm infill configuration and a lower pressure drop compared to the two 20 mm infill configurations in the series. Additionally, this study provides insights into the behavior of continuous fluid flow through TPMS infill in view of scaled‐up implementation, including the presence of recirculation zones that can be exploited to further enhance the collision rate of particles.

Turbulence as a Key Driver of Ice Aggregation and Riming in Arctic Low‐Level Mixed‐Phase Clouds, Revealed by Long‐Term Cloud Radar Observations

AbstractTurbulence in clouds is known to enhance particle collision rates, as widely demonstrated for warm rain formation. A similar impact on ice growth processes is expected but a solid observational basis is missing. A statistical analysis of a 15‐month data set of cloud radar observations allows for the first time to quantify the impact of turbulence on ice aggregation and riming in Arctic low‐level mixed‐phase clouds. Increasing eddy dissipation rate (EDR), from below 10−4 to above 10−3 m2 s−3, yields larger ice aggregates, and higher particle concentration, likely caused by increasing fragmentation. In conditions more favorable to riming, higher EDR is associated with dramatically higher particle fall velocities (by up to 125%), under similar liquid water paths, indicative of markedly higher degrees of riming. Our findings thus reveal the key role of turbulence for cold precipitation formation, and highlight the need for an improved understanding of turbulence‐hydrometeor interactions in cold clouds.

Large Area Picosecond Photodetector for the Upgrade II of the LHCb RICH

The Large Area Picosecond Photodetector (LAPPD) is a new generation microchannel plate (MCP) photodetector that is currently capturing the attention of detector physicists thanks to its excellent timing resolution, high gain and low dark count rate. It consists of a photocathode that responds to a single photon by producing a photoelectron that enters a pair of microchannel plates, where it is multiplied to a million electrons in a fast-rising pulse that is easily detected. The speed of this process provides the time resolution essential for high intensity environments such as the high luminosity Large Hadron Collider (LHC) at the CERN laboratory. The LAPPD has a large sensitive area, 200 × 200 mm 2, which reduces the number of sensors that are needed. The increased particle collision rate in the future LHC-HL phase represents a great challenge for the Ring Imaging Cherenkov (RICH) photodetectors in the LHCb experiment. The LHCb collaboration is working for the Upgrade II that will take place during the Long Shutdown 4 of the LHC (>2032). The LAPPD produced and commercialised by INCOM (US) represents one of the possible alternatives for the future RICH photodetectors, thanks to its good time resolution, low dark noise rate and high gain. The LHCb Edinburgh group is involved in the Upgrade II R&D programme, testing and characterising an LAPPD in the laboratory of the University. The LAPPD technology, the current Edinburgh setup and performance from first tests of the photodetector will be presented.

A One‐Dimensional Volcanic Plume Model for Predicting Ash Aggregation

AbstractDuring explosive volcanic eruptions, volcanic ash is ejected into the atmosphere, impacting aircraft safety and downwind communities. These volcanic clouds tend to be dominated by fine ash (<63 μm in diameter), permitting transport over hundreds to thousands of kilometers. However, field observations show that much of this fine ash aggregates into clusters or pellets with faster settling velocities than individual particles. Models of ash transport and deposition require an understanding of aggregation processes, which depend on factors like moisture content and local particle collision rates. In this study, we develop a Plume Model for Aggregate Prediction, a one‐dimensional (1D) volcanic plume model that predicts the plume rise height, concentration of water phases, and size distribution of resulting ash aggregates from a set of eruption source parameters. The plume model uses a control volume approach to solve mass, momentum, and energy equations along the direction of the plume axis. The aggregation equation is solved using a fixed pivot technique and incorporates a sticking efficiency model developed from analog laboratory experiments of particle aggregation within a novel turbulence tower. When applied to the 2009 eruption of Redoubt Volcano, Alaska, the 1D model predicts that the majority of the plume is over‐saturated with water, leading to a high rate of aggregation. Although the mean grain size of the computed Redoubt aggregates is larger than the measured deposits, with a peak at 1 mm rather than 500 μm, the present results provide a quantitative estimate for the magnitude of aggregation in an eruption.

Improvements to Stochastic Approaches for Simulating Agglomeration in Particle‐Laden Flows

AbstractDirect numerical simulation, facilitated by the spectral element method, has been used to study collisions and agglomeration in particle‐laden fluid flows through a channel at shear Reynolds number. The particulate phase is simulated by deterministic Lagrangian particle tracking alongside a stochastic technique to resolve collisions. Interestingly, the implementation of agglomeration in the deterministic approach caused a major reduction in the particle collision rate. To determine the cause of this effect, analysis of the collision rate and location across the channel for different particle properties was performed. For systems without agglomeration, collisions between particles tend to be located on consistent streamlines, occurring between single pairs of particles. A numerical analysis of this effect confirmed the result, where inter‐particle collisions across the channel mostly took place between consistent particle pairs in close proximity. Repeat collisions are shown to be almost eliminated with the addition of the agglomeration mechanism. The impact and implications of this effect on the accuracy of the stochastic technique is discussed, and modifications are suggested to offer improvements while accounting for the repeat collisions.

Vertically Aligned One-Dimensional Crystal-Structured Sb2Se3 for High-Efficiency Flexible Solar Cells via Regulating Selenization Kinetics

Recently, antimony selenide (Sb2Se3) has exhibited an exciting potential for flexible photoelectric applications due to its unique one-dimensional (1D) chain-type crystal structure, low-cost constituents, and superior optoelectronic properties. The 1D structure endows Sb2Se3 with a strong anisotropy in carrier transport and a lasting mechanical deformation tolerance. The control of the crystalline orientation of the Sb2Se3 film is an essential requirement for its device performance optimization. However, the current state-of-the-art Sb2Se3 devices suffer from unsatisfactory orientation control, especially for the (001) orientation, in which the chains stand vertically. Herein, we achieved an unprecedented control of the (001) orientation for the growth of the Sb2Se3 film on a flexible Mo-coated mica substrate by balancing the collision rate and kinetic energy of Se vapor particles with the surface of Sb film by regulating the selenization kinetics. Based on this (001)-oriented Sb2Se3 film, a high efficiency of 8.42% with a record open-circuit voltage (VOC) of 0.47 V is obtained for flexible Sb2Se3 solar cells. The vertical van der Waals gaps in the (001) orientation provide favorable diffusion paths for Se atoms, which results in a Se-rich state at the bottom of the Sb2Se3 film and promotes the in situ formation of the MoSe2 interlayer between Mo and Sb2Se3. These phenomena contribute to a back-surface field enhanced absorber layer and a quasi-Ohmic back contact, improving the device's VOC and the collection of carriers. This method provides an effective strategy for the orientation control of 1D materials for efficient photoelectric devices.

CFD-DEM study on the interaction between triboelectric charging and fluidization of particles in gas-solid fluidized beds

The triboelectric charging of particles has a substantial effect on the fluidization process and flow dynamics. Conversely, the fluidization state has an impact on the charge transfer of particles. In this study, we added models to the Computational Fluid Dynamic coupled Discrete Element Method (CFD-DEM) that could describe particle triboelectric charging and electrostatic interactions. According to numerical simulations, it was found for the first time that the behavior of bipolar charged particles undergoes a chaining stage and an agglomerating stage as the particle size decreases, where the interactions between triboelectric charging and fluidization are different. The effect of the wall is positively correlated to the ratio of the particle-wall collision rate to the total particle collision rate, and the fluidization conditions of particles affect the charging process of particles via this correlation. This work contributes to further understanding of the interaction between triboelectric charging and fluidization of particles and helps to optimize fluidized bed design and fluidization conditions.

Study on the slagging trends of the pre-combustion chamber in industrial pulverized coal boiler under different excess air coefficients by CFD numerical simulation

Industrial pulverized coal boiler (IPCB) usually employ a pre-combustion chamber (PCC) to ensure stable combustion. However, there are potential slagging problems in the PCC, which threaten the safe operation of the boiler. Air-staged combustion technology is an effective method to reduce NOx emission, and different excess air coefficients affect slagging. CFD simulation is an important method for studying slagging in the furnace. In this paper, Using CFD simulation with slagging model embedded by user-defined function, the slagging trends in the PCC are studied under different excess air coefficients. The particles' collision rate (CR), deposition rate (DR), and sticking rate (SR) on the PCC walls are obtained by CFD simulation. As the excess air coefficient in the PCC (αPCC) increases from 0.3 to 1.1, the CR of particles on all walls of PCC (CRPCC) increases exponentially, indicating the abrasion degree on the PCC walls is growing rapidly. In addition, due to the sizeable annular recirculation zone when αPCC = 0.3, the DR of particles on all walls of PCC (DRPCC) and the SR of particles on all walls of PCC (SRPCC) are much larger than those under other conditions, the slagging trend is the worse. As the αPCC increases from 0.3 to 1.1, the slagging trends in the PCC first decrease and then increase. To improve the safety of boiler operation, it is recommended that the αPCC be adjusted between 0.5 and 0.9.

Kinematic simulations as a subgrid-scale model for particle motion: a priori LES of homogeneous isotropic turbulence

In this work we investigated the impact of filtering and subgrid-scale modeling on particle settling velocity and collision-related statistics in a turbulent flow. To reduce the complexity of this task we first studied the motion of inertial particles in the low-pass filtered homogeneous and isotropic turbulence, which was subsequently enriched with the subgrid-scale velocity components obtained from a frozen high-pass filtered velocity field. Particular emphasis has been put on the radial distribution function and the radial relative velocity of nearly touching particles both in the presence or absence of the gravitational settling. These statistics are the key input parameters to the kinematic collision kernel which is of crucial importance in determining the collision rate of inertial particles in a turbulent flow. Furthermore, kinematic simulations were selected as a means of enhancing the fluid velocity at particle locations. We analyzed a wide range of Stokes numbers, i.e. a measure of particle inertia, and, in contrast to scientific premises found in the literature, we observed no improvement of particle statistics when the low-pass filtered fluid velocity was enriched with both a synthetic or spectrally-filtered small-scale structures. We discuss the shortages of any frozen-velocity-based subgrid-scale model in predicting both single- and two-point particle statistics. We also indicate that in some cases, particularly concerning the collision rate of particles suspended in homogeneous and isotropic turbulence, subgrid-scale contribution in the particle equation of motion can be neglected.

Tethered hard spheres: A bridge between the fluid and solid phases.

The thermodynamics of hard spheres tethered to a Face-Centered Cubic (FCC) lattice is investigated using event-driven molecular-dynamics. The particle-particle and the particle-tether collision rates are related to the phase space geometry and are used to study the FCC and fluid states. In tethered systems, the entropy can be determined by at least two routes: (i) through integration of the tether collision rates with the tether length rT or (ii) through integration of the particle-particle collision rates with the hard-sphere diameter σ (or, equivalently, the density). If the entropy were an entirely analytic function of rT and σ, these two methods for calculating the entropy should lead to the same results; however, a non-analytic region exists as an extension of the solid-fluid phase transition of the untethered hard-sphere system, and integration paths that cross this region will lead to values for the entropy that depend on the particular path chosen. The difference between the calculated entropies appears to be related to the communal entropy, and the location of the non-analytic region appears to be related to conditions where the regions of phase space associated with the FCC configuration become separated from those associated with the disordered fluid. The non-analytic region is finite in extent, vanishing below rT/a ≈ 0.55, where a is the lattice spacing, and there are many continuous paths that connect the fluid and solid phases that can be used to determine the crystal free energy with respect to the fluid.

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

Intricate relations among particle collision, relative motion and clustering in turbulent clouds: computational observation and theory

Abstract. Considering turbulent clouds containing small inertial particles, we investigate the effect of particle collision, in particular collision–coagulation, on particle clustering and particle relative motion. We perform direct numerical simulation (DNS) of coagulating particles in isotropic turbulent flow in the regime of small Stokes number (St=0.001–0.54) and find that, due to collision–coagulation, the radial distribution functions (RDFs) fall off dramatically at scales r∼d (where d is the particle diameter) to small but finite values, while the mean radial component of the particle relative velocity (MRV) increases sharply in magnitude. Based on a previously proposed Fokker–Planck (drift-diffusion) framework, we derive a theoretical account of the relationship among particle collision–coagulation rate, RDF and MRV. The theory includes contributions from turbulent fluctuations absent in earlier mean-field theories. We show numerically that the theory accurately accounts for the DNS results (i.e., given an accurate RDF, the theory could produce an accurate MRV). Separately, we also propose a phenomenological model that could directly predict MRV and find that it is accurate when calibrated using fourth moments of the fluid velocities. We use the model to derive a general solution of RDF. We uncover a paradox: the past empirical success of the differential version of the theory is theoretically unjustified. We see a further shape-preserving reduction of the RDF (and MRV) when the gravitational settling parameter (Sg) is of order O(1). Our results demonstrate strong coupling between RDF and MRV and imply that earlier isolated studies on either RDF or MRV have limited relevance for predicting particle collision rate.

Gasdynamic electron cyclotron ion sources: Basic physics, applications, and diagnostic techniques.

The gasdynamic electron cyclotron resonance (ECR) ion source is a type of the device in which the ionization efficiency is achieved primarily due to a high plasma density. Because of a high particle collision rate, the confinement is determined by a gasdynamic plasma outflow from a magnetic trap. Due to high efficiency of resonant heating, electrons gain energy significantly higher than that in inductively or capacitively coupled plasmas. As a consequence of such a parameter combination, the gasdynamic ECR plasma can be a unique source of low to medium charged ions, providing a high current and an ultimate quality of an ion beam. One of the most demanded directions of its application today is a development of high-current proton injectors for modern accelerators and neutron sources of different intensities. Special plasma parameters allow for the use of diagnostic techniques, traditional for multiply charged ECR plasmas as well as for other types of discharges with a high plasma density. Among the additional techniques, one can mention the methods of numerical simulation and reconstruction of the plasma density and temperature from the parameters of the extracted ion beams. Another point is that the high plasma density makes it possible to measure it from the Stark broadening of hydrogen lines by spectroscopy of plasma emission in the visible range, which is a fairly convenient non-invasive diagnostic method. The present paper discusses the main physical aspects of the gasdynamic ECR plasma, suitable diagnostic techniques, and possibilities and future prospects for its various applications.

Relativistic Magnetized Astrophysical Plasma Outflows in Black-Hole Microquasars

In this work, we deal with collimated outflows of magnetized astrophysical plasma known as astrophysical jets, which have been observed to emerge from a wide variety of astrophysical compact objects. The latter systems can be considered as either hydrodynamic (HD) or magnetohydrodynamic (MHD) in nature, which means that they are governed by non-linear partial differential equations. In some of these systems, the velocity of the jet is very high and they require relativistic MHD (RMHD) treatment. We mainly focus on the appropriate numerical solutions of the MHD (and/or RMHD) equations as well as the transfer equation inside the jet and simulate multi-messenger emissions from specific astrophysical compact objects. We use a steady state axisymmetric model assuming relativistic magnetohydrodynamic descriptions for the jets (astrophysical plasma outflows) and perform numerical simulations for neutrino, gamma-ray and secondary particle emissions. By adopting the existence of such jets in black hole microquasars (and also in AGNs), the spherical symmetry of emissions is no longer valid, i.e., it is broken, and the system needs to be studied accordingly. One of the main goals is to estimate particle collision rates and particle energy distributions inside the jet, from black-hole microquasars. As concrete examples, we choose the Galactic Cygnus X-1 and the extragalactic LMC X-1 systems.

Hybrid quantum classical graph neural networks for particle track reconstruction

The Large Hadron Collider (LHC) at the European Organisation for Nuclear Research (CERN) will be upgraded to further increase the instantaneous rate of particle collisions (luminosity) and become the High Luminosity LHC (HL-LHC). This increase in luminosity will significantly increase the number of particles interacting with the detector. The interaction of particles with a detector is referred to as “hit”. The HL-LHC will yield many more detector hits, which will pose a combinatorial challenge by using reconstruction algorithms to determine particle trajectories from those hits. This work explores the possibility of converting a novel graph neural network model, that can optimally take into account the sparse nature of the tracking detector data and their complex geometry, to a hybrid quantum-classical graph neural network that benefits from using variational quantum layers. We show that this hybrid model can perform similar to the classical approach. Also, we explore parametrized quantum circuits (PQC) with different expressibility and entangling capacities, and compare their training performance in order to quantify the expected benefits. These results can be used to build a future road map to further develop circuit-based hybrid quantum-classical graph neural networks.

Tethered-particle model: The calculation of free energies for hard-sphere systems

Two methods for computing the entropy of hard-sphere systems using a spherical tether model are explored, which allow the efficient use of event-driven molecular-dynamics simulations. An intuitive derivation is given, which relates the rate of particle collisions, either between two particles or between a particle and its respective tether, to an associated hypersurface area, which bounds the system's accessible configurational phase space. Integrating the particle-particle collision rates with respect to the sphere diameter (or, equivalently, density) or the particle-tether collision rates with respect to the tether length then directly determines the volume of accessible phase space and, therefore, the system entropy. The approach is general and can be used for any system composed of particles interacting with discrete potentials in fluid, solid, or glassy states. The entropies calculated for the liquid and crystalline hard-sphere states using these methods are found to agree closely with the current best estimates in the literature, demonstrating the accuracy of the approach.

Effect of magnetic field-dependent effective thermal conductivity of melted layer on nanosecond laser ablation of copper and formation of nanoparticles at atmospheric air pressure

For a nanosecond laser ablation of metals, the key physical phenomena involved are thermal evaporation, melt ejection, instability of the molten metal, etc., which depend on the initial temperature evolution in the metal. Understanding the evolution of temperature of the metal needs an effective simulation. In the present paper, we report on the finite element method-based simulation of nanosecond laser ablation of copper in the absence and presence of the magnetic field. Our studies showed that the effective thermal conductivity of the melted layer on the copper surface in the presence of the magnetic field affects the viscosity of the layer, mass ablation rate, instability, and then particle formation. The calculations showed that the condensed nuclei of large critical size are produced in the magnetic field. It is attributed to an increase in the collision rate of plasma particles in the magnetically confined plasma. The simulations are in good agreement with the experimentally measured values.

Structural subgrid scale model based on wavelet filter for large eddy simulation of particle-laden turbulence

Coherent vortex structures in subgrid scale (SGS) motions are important for preferential concentration and collision of particles with small and intermediate Stokes number (Xiong et al., 2019). In this study, a new large eddy simulation (LES) strategy that combines wavelet filtered large eddy simulation with a differential filter SGS model is proposed. First, a wavelet filtered direct numerical simulation (WFDNS) enables good preservation of SGS structures due to the high compressibility and local fidelity of the wavelet filter. Second, a differential filter (DF) model, containing only a parameter related to the nominal filter width, is used to dynamically reconstruct unresolved eddies. The SGS model presented here is verified using direct numerical simulation (DNS) data for particle-laden homogenous isotropic turbulence. Compared to the classical spectral-filtered DNS (FDNS) model, the new SGS model enables to achieve better agreement with DNS results related to the dispersed-phase statistics, such as particle acceleration, particle-seen fluid kinetic energy, particle-seen Lagrangian integral time etc. Furthermore, obtained results also exhibit the advantages of a proposed model over the stochastic Langevin model in the prediction of particle dynamics, especially collision-related statistics such as radial distribution function, radial relative velocities, and particle collision rates. Moreover, the model also demonstrates good performance in application to particle-pair related statistics at different Reynolds numbers and different filter depths. This is attributed to the ability of this model to multispectrally enhance coherent vortex structures in SGS, which is promising to apply for recovering SGS effects in the LES of turbulent particle-laden flows.