Abstract
In this article, a full-scale computational model for fluid–particle interaction under a magnetic field is developed. In this model, the fluid field is solved by the lattice Boltzmann method, and the hydrodynamic force acting on the particle is computed by the immersed boundary method . The numerical solutions of the magnetic field in the fluid–solid domain are achieved by the immersed interface method with a finite difference scheme, in which the normal and tangential jump conditions of the magnetic field intensity are applied to modify the standard finite difference scheme. The magnetic stress tensor along the fluid–particle interface can be calculated accurately. Unlike the widely used point–dipole model, the magnetic force acting on the particle is determined by the stress integration method. Numerical simulation of several numerical tests are carried out to validate the proposed model. The numerical results demonstrate the validity of the present model. Moreover, the magnetoviscous effect is studied by simulating the motion of elliptical particles under the uniform magnetic field in shear flow.
Introduction
Fluid–particle two-phase flows under a magnetic field can be usually found in nature and engineering applications, ranging from mineral screening (Ku et al., 2015), microfluidic control systems (Cao et al., 2014), magnetorheological fluid (Climent et al., 2004), and other chemical and biological applications. For the modeling of such flows, the two most important issues are to calculate the interaction forces between the particles and fluid/magnetic fields accurately and efficiently.
Fluid–particle interaction is one of the research hotspots in the computational fluid community. The existing models for fluid–particle interaction can fall into two categories: the Eulerian–Eulerian models and the Eulerian–Lagrangian models (Chiesa et al., 2005; Patel et al., 2017). The Eulerian–Eulerian models can capture the collective behavior of particles. However, it needs complex constitutive equations, and it is difficult to predict the details of the flow. Alternatively, the Eulerian–Lagrangian models easily obtain the detailed flowing behavior around the particles. The Eulerian–Lagrangian models can also be classified into two types: the point source model and the full-scale model (Luo et al., 2007; Hu et al., 2018). For the point source model, the empirical drag force formula is utilized to compute the fluid-particle interaction force. However, the point model i lacks enough accuracy for cases with dense or large-sized particles. For full-scale simulation, the fluid is governed by the Navier–Stokes equations, and the particle motion is solved by the Newtonian laws. The interaction force between the fluid and particle is determined by the no-slip condition. From this point of view, the full-scale model is most suitable for the study of flow mechanisms. Within the framework of full-scale models, the fixed mesh methods, which do not require time-consuming mesh generation, have received considerable attention in recent years. Unlike the body-fitted mesh methods, the governing equations of the fluid field are discretized on a fixed mesh, and the boundary is tracked by a set of Lagrangian points or captured by an implicit function. As a result, the computational efficiency is greatly improved. Kang et al. applied the distributed Lagrange multiplier/fictitious domain method to solve flows with suspended paramagnetic particles, in which the no-slip boundary condition on the particle boundary is implemented by Lagrange multipliers (Kang et al., 2008). Kang and Suh proposed the one-stage smoothed profile method for simulation of flows with suspended paramagnetic particles, in which the sigmoid function was used to construct body force by ensuring the rigidity of particles (Kang and Suh, 2011a). Kim and Park presented a level-set method for the analysis of magnetic particle dynamics on a fixed mesh (Kim and Park, 2010). Ke et al. introduced an IBM to simulate the behavior of magnetic particles in a fluid with an external magnetic field (Ke et al., 2017). These works show the effectiveness of fixed mesh methods in handling the fluid–particle interaction under a magnetic field.
In terms of calculation of magnetic force, the calculation models can also be divided into the point–dipole model and the force integration model. For example, Sand et al. developed a point–dipole model to simulate the magnetic particle suspension flow in which the single particle attraction force towards the magnetic pole was computed by a simple Kelvin force formula, and the interaction force between two nearby particles was computed by a magnetic dipole model (Sand et al., 2016). Like the point source model for fluid–particle interaction, the point–dipole model for particle–magnetic field interaction also suffers from inaccurate force estimation in nondilute flow or particle shape-dependent conditions. To this end, some work based on the force integration model has been done. In these models, the governing equation of the magnetic field in the multimedia zone was solved firstly. Then, the magnetic force acting on the particle can be obtained by integrating the force density. The magnetic force density can be calculated by the Helmholtz force density (Kang et al., 2008) or the virtual air gap scheme (Kang and Suh, 2011b). It should be pointed out that the fluid–particle interface is smeared out over several mesh cells in these above-the-force integration models. However, the magnetic stress force at the interface is treated as the continuous smoothed form, which will cause the loss of numerical accuracy.
In this paper, we develop a fully resolved simulation method for the fluid–particle interaction under a magnetic field. The LBM, which is a simple and efficient flow field solver, is adopted. The fluid–particle interaction is handled by the momentum-exchange-based IBM. More importantly, we give a Maxwell stress integration method to calculate the magnetic force based on the IIM with a finite difference scheme. According to the normal and tangential jump conditions of magnetic scalar potential along the medium interface, the discretized difference scheme is modified to ensure second-order accuracy. Then, the Maxwell stress tensor at the interface can be calculated accurately. The magnetic force acting on the particle can be obtained by integrating the magnetic stress force along the interface. Several numerical examples are simulated to validate the present IIM–IBM–LBM coupling model. The numerical results indicate that the present calculated stress force values agree well with the numerical results obtained by the body-fitted mesh. The circular particle sedimentation and motion of elliptical particles in shear flow under the magnetic field are also studied.
Mathematical Model and Numerical Method
In this study, LBM is used to simulate fluid flow, IBM is used to simulate the interaction between fluid and particles, and IIM is used to calculate magnetic field. In this section, we introduce the numerical implementation procedures of LBM, IBM, and IIM.
Lattice Boltzmann Model for Incompressible Fluid Flow
Using a standard uniform Cartesian grid with lattice space , the single relaxation time LB evolution equation with external force term is as follows (Chen and Doolen, 1998):where is the density distribution function. The dimensionless relaxation time readswhere is the coefficient of kinetic viscosity, is the equilibrium distribution function, and is the external force term in discrete velocity space. Then, and could be separately written aswhere is the lattice sound speed and is the body force. In the D2Q9 model, the discrete velocity and weight coefficients separately reads
The density and macrovelocity are written as
Immersed Boundary Method for the Fluid–Particle Interaction
The momentum exchange-based IBM proposed by Niu et al. (2006) is used to calculate the interaction force between fluid and solid particles. In Figure 1, the flow domain is covered with a uniform Cartesian grid, and the fluid–solid interface is divided into a series of Lagrange points with the arc length . The relation between the body force on the Eulerian nodes and the force density on the Lagrange points is expressed aswhere is the two-dimensional discrete Dirac function. In this study, the expression of proposed by Peskin (Peskin, 2002) is adopted:where
FIGURE 1
Configuration of the grids used in the IBM. Euler points (circles) represent the flow field and Lagrange points (crosses) represent the fluid–solid interface.
The calculation of the force density is the key problem in IBM, and it is usually determined by the no-slip condition on the boundary of the rigid particle. For the present momentum exchange-based IBM, the no-slip condition is implemented by the bounce-back rule. Firstly, the density distribution function at Lagrange points is obtained by an interpolation method, i.e.,
Then, the new distribution function could be solved according to the bounce-back rule:where is the velocity at Lagrange point and is the opposite direction of . Next, the force density of the flow field acting on the boundary point is calculated as
Furthermore, we define
The hydrodynamic force and moment acting on the solid particle can be obtained by removing the effect of internal masswhere is the centroid of the solid particle, is the particle density, is the mass of particle, is the moment of inertia, is the translational velocity, and is the rotational velocity, respectively.
Immersed Interface Method for Magnetic Field Calculation
The magnetic field in the multimedia zone can be solved by the Maxwell equations with the interface conditions. As shown in Figure 2, we consider two types of methods to impose the magnetic field: one is the given external uniform magnetic field, and the other is the permanent magnetic field. The Maxwell equations without current are expressed as follows:where is magnetic induction and is magnetic field intensity. To solve the two Maxwell’s equations, the constitutive relationships between magnetic flux density and magnetic field (Eq. 15) are neededwhere is permeability of vacuum, is relative permeability, and is the remanence of permanent magnet. At a current-free interface between two media, the following two continuity conditions are satisfied:where and are the normal and tangential directions, respectively. At the interface between permanent magnet and other media, the normal component of is still continuous across the interface. However, the tangential component of has a jump discontinuity.
FIGURE 2
Interface problems of permanent magnetic field (A) and uniform magnetic field (B).
Due to the irrotationality condition of field, the scalar magnetic potential is introduced:
The equations for solving the magnetic potential can be obtained:
The interface conditions for could be obtained from Formula 19 as follows:
In order to solve the magnetic field and compute magnetic force accurately, the IIM is used. The basic idea of IIM is to adopt the interface jump conditions to modify the finite difference scheme near the interface. As a result, second-order solutions can be achieved in the whole domain. A finite difference scheme for Poisson equation (Eq. 19) can be written asfor use at the point . To sum over k involves points (not more than six points selected near the interface) neighboring point , so the value of each is set to . The coefficients and indices depend on . The correction term is nonzero only if the grid point is irregular. For more details, please refer to LeVeque et al. (1994).
The following formula is used to solve the magnetic field force.
The Maxwell stress tensor iswhere is the identity operator. Moreover, the magnetic field force is solved by the following formula (Blūms et al., 1997):where is the normal vector. So, the magnetic field intensity at the interface needs to be determined.
Here, we introduce a six-point interpolation method to calculate from the correct side of the interface. In Figure 3, for the point P on the interface, we find the nearest point in the Cartesian grid. So, the stencil for six points is , and they transform them into interface coordinates as follows.
FIGURE 3
Selection of six points in the modified six-point method. Select the Euler point 1 closest to the Lagrange point (P), and the sixth point is the Euler point closest to the Lagrange point (P) except points 1 to 5.
Applying the Taylor expansion method for each point, we have
Ignoring the high order term , we can obtain a six-variable linear equation for.
. Once are obtained, we can get at point P on the interface.
Particle Dynamics
In this study, the dynamics of particle adopts Newton’s equation of motion, and the equations controlling particle translation and rotation are as follows:where is the rotation angle, is the gravity force, and and are the hydrodynamics force and magnetic force acting on the particle. To avoid the numerical instability due to the low-density ratio, the following time-stepping scheme for Eq. 29 is utilized:where , , and (Hu et al., 2015). We could easily obtain the center of mass , the rotation angle , the translational velocity , and the rotational velocity of particle at time .
Then, the velocity at the boundary point can be written as
Numerical Results and Discussion
In this section, the present IIM–IBM–LBM model is used to simulate several problems with fluid–particle–magnetic interaction.
Numerical Method Validation
Particle Dynamics Verification of IB–LBM
In order to examine the accuracy of the present IB–LBM, the free settlement of a single particle is simulated, which has been used by some scholars as a benchmark problem (Glowinski et al., 2001; Hu et al., 2015). In Figure 4A, we use a fluid domain with a width of 2 cm and a height of 6 cm, where the fluid density is . A circular particle with a diameter of 0.25 cm falls freely from the position of 1 cm and 4 cm under the effect of gravity. The particle density is set as , and the corresponding kinematic viscosity is , respectively.
FIGURE 4
Model of free settlement of a single particle (A). Model of permanent magnet attracting magnetic conducting cylinder (B).
A grid is used for the numerical simulation. Figures 5A,B show the variations in vertical position and velocity of the particle center with time when and , respectively. Obviously, the particle quickly reaches the terminal velocity, and the velocity becomes small when the particle is close to the bottom wall at about . The results are in good agreement with those of Hu et al. (Hu et al., 2015) and Glowinski et al. (Glowinski et al., 2001). Similarly, we simulate the case with and . As shown in Figures 5C,D, the desired results are obtained.
FIGURE 5
The position of the center of the particle for , (A);, (C), and the velocity of that for , (B); , (D).
Verification of Magnetic Field Calculation of IIM
To test the accuracy of IIM in terms of calculation of magnetic field, the problem of a permanent magnet attracting a magnetic conducting cylinder is studied. The results obtained are compared with those calculated by COMSOL in which the body-fitted mesh is used. As shown in Figures 4A,B, a circular permanent magnet with remanence and a circular magnetic conducting object with are centered at (2 cm, 1.5 cm) and (2 cm, 3 cm), respectively. Both circular solids are 0.8 cm in diameter. The computational domain is a square with a side length of 4 cm.
This numerical simulation is carried out in a grid. Figure 6 displays the distribution of magnetic potential , and it can be seen clearly that the jump of the magnetic field is in the normal direction across the interface. Figure 7 shows the magnetic field distribution. It could be observed that the closer the permanent magnet, the stronger the magnetic field intensity. Also, the contours of across magnetic conducting cylinder interface have a jump. We can see that the present results agree well with those of COMSOL (Figure 7). The horizontal and vertical components of magnetic stress along the interface are plotted in Figure 8, and the symmetrical distribution of them can be found. The numerical results of the total magnetic force acting on the magnetic conducting cylinder are shown in Table 1. As expected, compared with the results obtained by the diffusion interface method, IIM based on the non-body–fitted grid technology can obtain the calculation accuracy similar to that of COMSOL based on the body-fitted grid technology.
FIGURE 6
The magnetic potential of the permanent magnet attracting the particle model is obtained by the IIM method.
FIGURE 7
The permanent magnet attracts the magnetic field of the magnetic conducting cylinder model, which is calculated by COMSOL (A) and IIM (B).
FIGURE 8
Magnetic stress along the interface of the circular conducting magnet.
TABLE 1
| The position of permanent magnet center in the y direction (cm) | IIM | Comsol | Diffuse interface method | |||
|---|---|---|---|---|---|---|
| 1.5 | 0.0084 | −3.0857 | −0.0068 | −3.1410 | −0.9100 | −5.0882 |
| 1.25 | 0.0027 | −1.0788 | −0.0022 | −1.0982 | −0.2429 | −1.8545 |
| 1.0 | 0.00095 | −0.4015 | −0.00081 | −0.4081 | −0.0731 | −0.7075 |
| 0.75 | 0.00034 | −0.1488 | −0.00030 | −0.1508 | −0.0231 | −0.2662 |
The magnetic force on the magnetic conducting cylinder.
Particle Sedimentation Under Permanent Magnetic Field
The diagram of particle sedimentation under a permanent magnetic field is shown in Figure 9. The computational domain for the magnetic field is 8 cm wide and 8 cm high. The computational domain for the flow field is 2 cm wide and 6 cm high. The circular particle with and diameter of 0.25 cm is located (4 cm, 6 cm), and the circular permanent magnet with diameter of 0.8 cm is located (4 cm, 1 cm). The remanence of the permanent magnet is . In this study, the relative permeability of the flow field is set to 1.0.
FIGURE 9
Model of particle sedimentation under the action of permanent magnet.
A grid is used for the present numerical simulation. The flow field simulation is carried out on a grid. Figure 10 shows the variations in vertical position and velocity of the particle center with time when and , respectively. Due to the coexistence of magnetic force and gravity force, the particle settling velocity becomes larger compared with that without the effect of a permanent magnet. The closer the particle is to the permanent magnet, the greater the magnetic force and the faster the particle velocity. When , the acceleration effect of magnetic force on particles is dominant. Figure 11 displays the variety of magnetic stress at the particle interface at . Similarly, we also study the case with and and reached similar conclusions (as shown in Figure 12).
FIGURE 10
Position (A) and velocity (B) of the center of the particle for and .
FIGURE 11
The magnetic stress in x-direction (A) and in y-direction (B) on the particle interface at .
FIGURE 12
Position (A) and velocity (B) of the center of the particle for and .
Shear Viscosity of Suspension Containing Elliptical Particles Under the Magnetic Field
It should be pointed out that the point–dipole model is suitable for spherical-like particles. We also consider an ellipsoidal particle immersed in the two-dimensional shear flows, as shown in Figure 13. In order to verify the reliability of the numerical method, we calculate an example of an elliptic particle rotating in a simple shear flow when Reynolds number and the ratio of major axis to minor axis of ellipse . The curve of angular velocity changing with angular rotation is obtained in Figure 14 and compared with the exact solution calculated by Jeffery (Jeffery, 1922). The simulation results of this paper agree well with Jeffery’s exact solution.
FIGURE 13
An ellipsoidal particle in the two-dimensional shear flows under uniform magnetic field.
FIGURE 14
Comparison of Jeffery solution (Jeffery, 1922) and the present simulation result.
The relation between angular velocity and angle of the elliptic particle in shear flow obtained by Jeffery (Jeffery, 1922) is as follows:where the fluid shear rate and are the length of long and short half axle.
Then, we study the shear viscosity of suspension containing elliptical particles under the magnetic field. The computational domain is 2 cm long and 1 cm wide, in which the velocities of the upper and lower planes are and , respectively. The elliptical particle with length–width ratio and relative permeability is located at the center of the computational domain. Under an external uniform magnetic field , the force acting on the paramagnetic ellipsoid particle is zero, but the torque is not zero.
The shear stress at the fluid node is calculated as
To study the simple rheological properties of a suspension containing elliptical particles under the external magnetic field, the effective viscosity of the suspension is calculated:where is the distance between two plates and is the speed of the top plate. is the average shear force acting on the moving plane wall (Huang et al., 2012).
A grid of is used for the calculation domain, and the left and right boundaries are subjected to the periodic boundary conditions. We set the initial deflection angle of the ellipse particle to . When the flow field becomes steady state, the effective shear viscosity of the suspension containing the elliptical particle under the magnetic field can be obtained. Through the present direct numerical simulation, the effective shear viscosities under different magnetic field strength , Reynolds number , and volume fraction of elliptical particles in suspension are compared, as shown in Figure 15. Obviously, the effective kinematic viscosity increases with an increase in Reynolds number and decreases with an increase in volume fraction . Moreover, the effective kinematic viscosity increases with the increase of magnetic field strength .
FIGURE 15
Variation of effective kinematic viscosity of suspension with Reynolds number in different magnetic field strength when is 5% (A). Variation of effective kinematic viscosity of suspension with volume fraction in different magnetic field strength when is 20 (B).
Conclusion
The fluid–particle–magnetic interactions are modeled using the IIM–IBM–LBM coupling method. The fluid flow simulations are handled by the simple and efficient LBM. The particle motion and the hydrodynamics interaction between the particle and the flow field are computed by the momentum exchange-based IBM. Especially, we use the IIM to solve the magnetic field and calculate the magnetic force with the aid of the interface jump conditions. Unlike the point-source model or point dipole model, the hydrodynamics and magnetic forces acting on the particle are calculated using an integration method, in which the flow details and magnetic distribution around the particle are considered. Two numerical examples are simulated to verify the numerical accuracy of the present full-scale model. Moreover, particle sedimentation under a permanent magnetic field and shear viscosity of suspension containing elliptical particles under the magnetic field are also studied by the present model. The obtained results indicate that the present model has the potential to treat the complex fluid–particle–magnetic interactions.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was supported by the National Natural Science Foundation of China (grant nos. 12172039, 12102228, and 11802159) and Fundamental Research Funds for the Central Universities (grant no. 2020RC201).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
BlūmsE.T͡SebersA. O.CebersA. O.MaiorovM. M. (1997). Magnetic Fluids. Berlin, Germany: Walter de Gruyter.
2
CaoQ.HanX.LiL. (2014). Configurations and Control of Magnetic Fields for Manipulating Magnetic Particles in Microfluidic Applications: Magnet Systems and Manipulation Mechanisms. Lab. Chip14 (15), 2762–2777. 10.1039/c4lc00367e
3
ChenS.DoolenG. D. (1998). Lattice Boltzmann Method for Fluid Flows. Annu. Rev. Fluid Mech.30 (1), 329–364. 10.1146/annurev.fluid.30.1.329
4
ChiesaM.MathiesenV.MelheimJ. A.HalvorsenB. (2005). Numerical Simulation of Particulate Flow by the Eulerian-Lagrangian and the Eulerian-Eulerian Approach with Application to a Fluidized Bed. Comput. Chem. Eng.29 (2), 291–304. 10.1016/j.compchemeng.2004.09.002
5
ClimentE.MaxeyM. R.KarniadakisG. E. (2004). Dynamics of Self-Assembled Chaining in Magnetorheological Fluids. Langmuir20 (2), 507–513. 10.1021/la035540z
6
GlowinskiR.PanT. W.HeslaT. I.JosephD. D.PériauxJ. (2001). A Fictitious Domain Approach to the Direct Numerical Simulation of Incompressible Viscous Flow Past Moving Rigid Bodies: Application to Particulate Flow. J. Comput. Phys.169 (2), 363–426. 10.1006/jcph.2000.6542
7
HuY.LiD.ShuS.NiuX. (2015). Modified Momentum Exchange Method for Fluid-Particle Interactions in the Lattice Boltzmann Method. Phys. Rev. E Stat. Nonlin Soft Matter Phys.91 (3), 033301. 10.1103/PhysRevE.91.033301
8
HuY.LiD.NiuX.ShuS. (2018). Fully Resolved Simulation of Particulate Flows with Heat Transfer by Smoothed Profile-Lattice Boltzmann Method. Int. J. Heat Mass Transf.126, 1164–1167. 10.1016/j.ijheatmasstransfer.2018.05.137
9
HuangH.YangX.KrafczykM.LuX.-Y. (2012). Rotation of Spheroidal Particles in Couette Flows. J. Fluid Mech.692, 369–394. 10.1017/jfm.2011.519
10
JefferyG. B. (1922). The Motion of Ellipsoidal Particles Immersed in a Viscous Fluid. Proc. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character102 (715), 161–179. 10.1098/rspa.1922.0078
11
KangS.SuhY. K. (2011). An Immersed-Boundary Finite-Volume Method for Direct Simulation of Flows with Suspended Paramagnetic Particles. Int. J. Numer. Meth. Fluids67 (1), 58–73. 10.1002/fld.2336
12
KangS.SuhY. K. (2011). Direct Simulation of Flows with Suspended Paramagnetic Particles Using One-Stage Smoothed Profile Method. J. Fluids Struct.27 (2), 266–282. 10.1016/j.jfluidstructs.2010.11.002
13
KangT. G.HulsenM. A.den ToonderJ. M. J.AndersonP. D.MeijerH. E. H. (2008). A Direct Simulation Method for Flows with Suspended Paramagnetic Particles. J. Comput. Phys.227 (9), 4441–4458. 10.1016/j.jcp.2008.01.005
14
KeC.-H.ShuS.ZhangH.YuanH.-Z. (2017). LBM-IBM-DEM Modelling of Magnetic Particles in a Fluid. Powder Technol.314, 264–280. 10.1016/j.powtec.2016.08.008
15
KimY. S.ParkI. H. (2010). FE Analysis of Magnetic Particle Dynamics on Fixed Mesh with Level Set Function. IEEE Trans. Magn.46 (8), 3225–3228. 10.1109/tmag.2010.2045747
16
KuJ.ChenH.HeK.YanQ. (2015). Simulation and Observation of Magnetic Mineral Particles Aggregating into Chains in a Uniform Magnetic Field. Miner. Eng.79, 10–16. 10.1016/j.mineng.2015.05.002
17
LeVequeR. J.LiZ. (1994). The Immersed Interface Method for Elliptic Equations with Discontinuous Coefficients and Singular Sources. SIAM J. Numer. Anal.31 (4), 1019–1044. 10.1137/0731054
18
LuoK.WangZ.FanJ.CenK. (2007). Full-Scale Solutions to Particle-Laden Flows: Multidirect Forcing and Immersed Boundary Method. Phys. Rev. E Stat. Nonlin Soft Matter Phys.76 (6), 066709. 10.1103/PhysRevE.76.066709
19
NiuX. D.ShuC.ChewY. T.PengY. (2006). A Momentum Exchange-Based Immersed Boundary-Lattice Boltzmann Method for Simulating Incompressible Viscous Flows. Phys. Lett. A354 (3), 173–182. 10.1016/j.physleta.2006.01.060
20
PatelR. G.DesjardinsO.KongB.CapecelatroJ.FoxR. O. (2017). Verification of Eulerian-Eulerian and Eulerian-Lagrangian Simulations for Turbulent Fluid-Particle Flows. AIChE J.63 (12), 5396–5412. 10.1002/aic.15949
21
PeskinC. S. (2002). The Immersed Boundary Method. Acta Numer.11, 479–517. 10.1017/s0962492902000077
22
SandA.StenerJ. F.ToivakkaM. O.CarlsonJ. E.PålssonB. I. (2016). A Stokesian Dynamics Approach for Simulation of Magnetic Particle Suspensions. Miner. Eng.90, 70–76. 10.1016/j.mineng.2015.10.015
Summary
Keywords
full-scale simulation, fluid-particle interactions, magnetic field, immersed interface method, immersed boundary method, lattice Boltzmann method
Citation
Peng W, Hu Y, Li D and He Q (2022) Full-Scale Simulation of the Fluid–Particle Interaction Under Magnetic Field Based on IIM–IBM–LBM Coupling Method. Front. Mater. 9:932854. doi: 10.3389/fmats.2022.932854
Received
30 April 2022
Accepted
03 June 2022
Published
05 August 2022
Volume
9 - 2022
Edited by
Xuan Shouhu, University of Science and Technology of China, China
Reviewed by
Huaxia Deng, University of Science and Technology of China, China
Xufeng Dong, Dalian University of Technology, China
Updates
Copyright
© 2022 Peng, Hu, Li and He.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yang Hu, yanghu@bjtu.edu.cn
This article was submitted to Smart Materials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.