Numerical simulation of particulate flows using a hybrid of finite difference and boundary integral methods

被引:2
作者
Bhattacharya, Amitabh [1 ]
Kesarkar, Tejas [1 ]
机构
[1] Indian Inst Technol, Dept Mech Engn, Bombay 400076, Maharashtra, India
关键词
DISCRETIZED BOLTZMANN-EQUATION; STOKES EQUATIONS; COMPLEX GEOMETRIES; DYNAMICS; SPHERES; SUSPENSIONS; PHYSALIS; SOLVER; ARRAY;
D O I
10.1103/PhysRevE.94.043309
中图分类号
O35 [流体力学]; O53 [等离子体物理学];
学科分类号
070204 ; 080103 ; 080704 ;
摘要
A combination of finite difference (FD) and boundary integral (BI) methods is used to formulate an efficient solver for simulating unsteady Stokes flow around particles. The two-dimensional (2D) unsteady Stokes equation is being solved on a Cartesian grid using a second order FD method, while the 2D steady Stokes equation is being solved near the particle using BI method. The two methods are coupled within the viscous boundary layer, a few FD grid cells away from the particle, where solutions from both FD and BI methods are valid. We demonstrate that this hybrid method can be used to accurately solve for the flow around particles with irregular shapes, even though radius of curvature of the particle surface is not resolved by the FD grid. For dilute particle concentrations, we construct a virtual envelope around each particle and solve the BI problem for the flow field located between the envelope and the particle. The BI solver provides velocity boundary condition to the FD solver at "boundary" nodes located on the FD grid, adjacent to the particles, while the FD solver provides the velocity boundary condition to the BI solver at points located on the envelope. The coupling between FD method and BI method is implicit at every time step. This method allows us to formulate an O(N) scheme for dilute suspensions, where N is the number of particles. For semidilute suspensions, where particles may cluster, an envelope formation method has been formulated and implemented, which enables solving the BI problem for each individual particle cluster, allowing efficient simulation of hydrodynamic interaction between particles even when they are in close proximity. The method has been validated against analytical results for flow around a periodic array of cylinders and for Jeffrey orbit of a moving ellipse in shear flow. Simulation of multiple force-free irregular shaped particles in the presence of shear in a 2D slit flow has been conducted to demonstrate the robustness of the FD-BI scheme. This method can also be used to solve the full Navier-Stokes equations around particles. In this case, it is shown that, below a threshold particle Reynolds number, which in turn depends on the required resolution of the particle surface, the FD-BI scheme will be more computationally efficient than a pure FD scheme solved on a multiblock grid.
引用
收藏
页数:18
相关论文
共 34 条
[1]   A fast solver for the Stokes equations with distributed forces in complex geometries [J].
Biros, G ;
Ying, LX ;
Zorin, D .
JOURNAL OF COMPUTATIONAL PHYSICS, 2004, 193 (01) :317-348
[2]   A sphere in a uniformly rotating or shearing flow [J].
Bluemink, J. J. ;
Lohse, D. ;
Prosperetti, A. ;
Van Wijngaarden, L. .
JOURNAL OF FLUID MECHANICS, 2008, 600 :201-233
[3]   STOKESIAN DYNAMICS [J].
BRADY, JF ;
BOSSIS, G .
ANNUAL REVIEW OF FLUID MECHANICS, 1988, 20 :111-157
[4]  
Crowe CT, 2011, Multiphase flows with droplets and particies
[5]   A lattice-Boltzmann method with hierarchically refined meshes [J].
Eitel-Amor, G. ;
Meinke, M. ;
Schroeder, W. .
COMPUTERS & FLUIDS, 2013, 75 :127-139
[6]  
Falgout R.D., 2006, Lect. Notes Comput. Sci. Eng., V51, P267, DOI DOI 10.1007/3-540-31619-1_8
[7]   Pursuing scalability for hypre's conceptual interfaces [J].
Falgout, RD ;
Jones, JE ;
Yang, UM .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2005, 31 (03) :326-350
[8]  
Haines E., 1994, GRAPHICS GEMS, P24
[10]  
Hernandez-Ortiz J. P., 2004, THESIS