View source for CFD-DEM

The '''CFD-DEM model''', or '''Computational Fluid Dynamics / Discrete Element Method model''', is a process used to model or simulate systems combining fluids with solids or particles. In CFD-DEM, the motion of discrete solids or particles phase is obtained by the [[Discrete Element Method]] (DEM) which applies [[Newton's laws of motion]] to every particle, while the flow of continuum fluid is described by the local averaged [[Navier–Stokes equations]] that can be solved using the traditional [[Computational Fluid Dynamics]] (CFD) approach. The interactions between the fluid phase and solids phase is modeled by use of [[Newton's third law]].

The direct incorporation of CFD into DEM to study the gas fluidization process so far has been attempted by Tsuji et al.<ref>{{cite journal | last1=Tsuji | first1=Y. | last2=Kawaguchi | first2=T. | last3=Tanaka | first3=T. | title=Discrete particle simulation of two-dimensional fluidized bed | journal=Powder Technology | publisher=Elsevier BV | volume=77 | issue=1 | year=1993 | issn=0032-5910 | doi=10.1016/0032-5910(93)85010-7 | pages=79–87}}</ref><ref>{{cite journal | last1=Tsuji | first1=Y. | last2=Tanaka | first2=T. | last3=Ishida | first3=T. | title=Lagrangian numerical simulation of plug flow of cohesionless particles in a horizontal pipe | journal=Powder Technology | publisher=Elsevier BV | volume=71 | issue=3 | year=1992 | issn=0032-5910 | doi=10.1016/0032-5910(92)88030-l | pages=239–250}}</ref> and most recently by Hoomans et al.,<ref>{{cite journal | last1=Hoomans | first1=B.P.B. | last2=Kuipers | first2=J.A.M. | last3=Briels | first3=W.J. | last4=van Swaaij | first4=W.P.M. | title=Discrete particle simulation of bubble and slug formation in a two-dimensional gas-fluidised bed: A hard-sphere approach | journal=Chemical Engineering Science | publisher=Elsevier BV | volume=51 | issue=1 | year=1996 | issn=0009-2509 | doi=10.1016/0009-2509(95)00271-5 | pages=99–118| s2cid=17460834 | url=https://research.utwente.nl/en/publications/discrete-particle-simulation-of-bubble-and-slug-formation-in-a-twodimensional-gasfluidised-bed-a-hardsphere-approach(fa42a6ac-0c72-4892-9a60-822f247cfbf4).html }}</ref> Deb et al.<ref>{{cite journal | last1=Deb | first1=Surya | last2=Tafti | first2=Danesh | title=Investigation of flat bottomed spouted bed with multiple jets using DEM–CFD framework | journal=Powder Technology | publisher=Elsevier BV | volume=254 | year=2014 | issn=0032-5910 | doi=10.1016/j.powtec.2014.01.045 | pages=387–402}}</ref> and Peng et al.<ref>{{Cite journal|doi=10.1002/aic.14421|title=Influence of void fraction calculation on fidelity of CFD-DEM simulation of gas-solid bubbling fluidized beds|year=2014|last1=Peng|first1=Z.|last2=Doroodchi|first2=E.|last3=Luo|first3=C.|last4=Moghtaderi|first4=B.|journal=AIChE J|volume=60|issue=6|page=2000}}</ref> A recent overview over fields of application was given by Kieckhefen et al.<ref>{{cite journal | last1=Kieckhefen | first1=P. | last2=Pietsch | first2=S. | last3=Dosta | first3=M.| last4=Heinrich| first4=S. | title=Possibilities and Limits of Computational Fluid Dynamics–Discrete Element Method Simulations in Process Engineering: A Review of Recent Advancements and Future Trends | journal=Annual Review of Chemical and Biomolecular Engineering | volume=11 | year=2020 | pages=397–422 | doi=10.1146/annurev-chembioeng-110519-075414| pmid=32169000 | s2cid=212707563 }}</ref>

==Parallelization==

[[OpenMP]] has been shown to be more efficient in performing coupled CFD-DEM calculations in parallel framework as compared to [[Message Passing Interface|MPI]] by Amritkar et al.<ref>{{Cite journal|doi=10.1016/j.jcp.2013.09.007|title=Efficient parallel CFD-DEM simulations using OpenMP|year=2014|last1=Amritkar|first1=Amit|last2=Deb|first2=Surya|last3=Tafti|first3=Danesh|journal=Journal of Computational Physics|volume=256|page=501|bibcode=2014JCoPh.256..501A}}</ref> Recently, a multi-scale parallel strategy<ref>{{Cite journal |last1=Lu |first1=L. |last2=Xu |first2=J. |last3=Ge |first3=W.|last4=Gao |first4=G.|last5=Jiang |first5=Y. |last6=Zhao |first6=M.|last7=Liu |first7=X. |last8=Li |first8=J. |doi=10.1016/j.ces.2016.08.013 |title=Computer virtual experiment on fluidized beds using a coarse-grained discrete particle method—EMMS-DPM |journal=Chemical Engineering Science |volume=155 |pages=314–337 |year=2016 }}</ref> is developed. Generally, the simulation domain is divided into many sub-domains and each process calculates only one sub-domain using MPI passing boundary information; for each sub-domain, the CPUs are used to solve the fluid phase while the general purpose GPUs are used to solve the movement of particles. However, in this computation method CPUs and GPUs work in serial. That is, the CPUs are idle while the GPUs are calculating the solid particles, and the GPUs are idle when the CPUs are calculating the fluid phase. To further accelerate the computation, the CPU and GPU computing can be overlapped using the [[shared memory]] of a Linux system. Thus, the fluid phase and particles can be calculated at the same time.

==Reducing computation cost using Coarse Grained Particles==

The computation cost of CFD-DEM is huge due to a large number of particles and small time steps to resolve particle-particle collisions. To reduce computation cost, many real particles can be lumped into a Coarse Grained Particle (CGP).<ref>{{Cite journal |last1=Lu |first1=L. |last2=Yoo |first2=K. |last3=Benyahia |first3=S. |doi=10.1021/acs.iecr.6b02688 |title=Coarse-Grained-Particle Method for Simulation of Liquid–Solids Reacting Flows |journal=Industrial & Engineering Chemistry Research |volume=55 |pages=10477–10491 |year=2016 |issue=39 |osti=1477859 }}</ref><ref name="Lu2014">{{Cite journal |last1=Lu |first1=L. |last2=Xu |first2=J. |last3=Ge |first3=W.|last4=Yue |first4=Y. |last7=Liu |first7=X. |last8=Li |first8=J. |doi=10.1016/j.ces.2014.08.004 |title=EMMS-based discrete particle method (EMMS–DPM) for simulation of gas–solid flows |journal=Chemical Engineering Science |volume=120 |pages=67–87 |year=2014 }}</ref> The diameter of the CGP is calculated by the following equation:
: <math> d_{CGP} = d_{realParticle} \cdot W^{1/3},</math>
where <math> W </math> is the number of real particles in CGP. Then, the movement of CGPs can be tracked using DEM.
In simulations using Coarse Grained Particles, the real particles in a CGP are subjected to the same drag force, same temperature and same species mass fractions. The momentum, heat and mass transfers between fluid and particles are firstly calculated using the diameter of real particles and then scaled by <math> W </math> times. The value of <math> W </math> is directly related to computation cost and accuracy.<ref>{{Cite journal |last1=Lu |first1=L. |last2=Konan |first2=A. |last3=Benyahia |first3=S. |doi=10.1016/j.cej.2017.06.002 |title=Influence of grid resolution, parcel size and drag models on bubbling fluidized bed simulation |journal=Chemical Engineering Journal |volume=326 |pages=627–639 |year=2017 |osti=1404697 |doi-access=free }}</ref> When <math> W </math> is equal to unity, the simulation becomes DEM-based achieving results that are of the highest possible accuracy. As this ratio increases, the speed of the simulation increases drastically but its accuracy deteriorates. Apart from an increase in speed, general criteria for selecting a value for this parameter is not yet available. However, for systems with distinct mesoscale structures, like bubbles and clusters, the parcel size should be small enough to resolve the deformation, aggregation, and breakage of bubbles or clusters. The process of lumping particles together reduces the collision frequency, which directly influences the energy dissipation. To account for this error, an effective restitution coefficient was proposed by Lu et al.,<ref name="Lu2014" /> based on kinetic theory of granular flow, by assuming the energy dissipation during collisions for the original system and the coarse grained system are identical.
: <math> e_{CGP} = \sqrt{1+\left(e_{realParticle}^2-1\right)W^{1/3}}</math>

== References ==
{{Reflist}}

{{DEFAULTSORT:Cfd-Dem}}
[[Category:Computational physics]]