A fictitious domain method for dynamic simulation of particle sedimentation in Bingham fluids

Yu, Zhaosheng; Wachs, Anthony

doi:10.1016/j.jnnfm.2007.02.007

Cited by 78 publications

(78 citation statements)

References 34 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In addition, we assume that the fluid is Newtonian and that the flow is isothermal (see our previous work on the flow of particles in a non-Newtonian fluid [14,19] or with heat transfer [13]). However, both smoothly shaped (circular cylinder in 2D and sphere in 3D) and angular (polygon in 2D and polyhedron in 3D) are considered in the following simulations.…”

Section: Resultsmentioning

confidence: 99%

“…Compared to [3], we add the correction step 6, as suggested by Glowinski et al [1] and used by us in [13,14] and we solve sub-problems 2 and 3 together as an advection/diffusion one when the CFL condition is fulfilled with the constant time step employed:…”

Section: Computational Featuresmentioning

confidence: 99%

“…As originally suggested by Glowinski et al in [1] as well as in our previous work [3,14], the matrix system (20) (or equivalently (31)) is solved by an Uzawa/CG algorithm. Details of this iterative procedure in basic serial implementation is given in Appendix A.…”

Section: Mpi Solution Of the Dlm/fd Saddle-pointmentioning

confidence: 99%

“…In this perspective, a fully MPI implementation running on large clusters is a real breakthrough. Compared to our previous work [3,13,14], we present here our new code PeliGRIFF. This new version of our DEM-DLM/FD approach is fully MPI and enables us to extend significantly the number of particles present in our systems.…”

mentioning

confidence: 99%

See 3 more Smart Citations

PeliGRIFF, a parallel DEM-DLM/FD direct numerical simulation tool for 3D particulate flows

Wachs

2010

J Eng Math

Self Cite

View full text Add to dashboard Cite

The problem of particulate flows at moderate to high concentration and finite Reynolds number is addressed by parallel direct numerical simulation. The present contribution is an extension of the work published in Computers & Fluids 38:1608(2009, where systems of moderate size in a 2D geometry were examined. At the numerical level, the suggested method is inspired by the framework established by Glowinski et al. (Int J Multiph Flow 25:755, 1999) in the sense that their Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) formulation and their operator-splitting idea are employed. In contrast, particle collisions are handled by an efficient Discrete Element Method (DEM) granular solver, which allows one to consider both smoothly (sphere) and non-smoothly (angular polyhedron) shaped particles. From a computational viewpoint, a basic though efficient strategy has been developed to implement the DLM/FD method in a domain decomposition/distributed fashion. To achieve this goal, the serial code, GRIFF (GRains In Fluid Flow; see Comput Fluids 38:1608-1628, 2009) is upgraded to fully MPI capabilities. The new code, PeliGRIFF (Parallel Efficient Library for GRains in Fluid Flow) is developed under the framework of the fully MPI open-source platform PELICANS. The parallel computing capabilities of PeliGRIFF offer new perspectives in the study of particulate flows and indeed increase the number of particles usually simulated in the literature. Solutions to address new issues raised by the parallelization of the DLM/FD method and assess the scalable properties of the code are proposed. Results on the 2D/3D sedimentation of a significant collection of isometric polygonal/polyhedral particles in a Newtonian fluid with collisions are presented as a validation test and an illustration of the class of particulate flows PeliGRIFF is able to investigate.

show abstract

Section: Resultsmentioning

confidence: 99%

Section: Computational Featuresmentioning

confidence: 99%

Section: Mpi Solution Of the Dlm/fd Saddle-pointmentioning

confidence: 99%

mentioning

confidence: 99%

See 2 more Smart Citations

PeliGRIFF, a parallel DEM-DLM/FD direct numerical simulation tool for 3D particulate flows

Wachs

2010

J Eng Math

Self Cite

View full text Add to dashboard Cite

show abstract

“…A number of numerical studies have been conducted investigating the hydrodynamic interaction between two spheres in Bingham fluids, formulating the problem with both fixed [26] and settling [17,27] spheres. The settling sphere studies , however, calculated Bi and Re based off the characteristic Stokes velocity of the settling particle, such that the terminal Bingham and Reynolds numbers varied between tests with the terminal settling velocity.…”

Section: Interaction Between Two Identical Spheresmentioning

confidence: 99%

Performance evaluation of an implicitly-regularised model for viscoplastic particle suspensions

Vaira¹

View full text Add to dashboard Cite

Firstly, I would like to show my appreciation towards my supervisor, Dr. Christopher Leonardi, for all of the support and advice over the last six months. Your thorough feedback and constant interest in and dedication to the project drove me to produce the highest quality work I could. My development throughout the project has been immeasurable, with the work opening my eyes to a world of future career possibilities, and for those things I am grateful.To Wojciech Regulski in Poland, thank you for the many code implementations, constant assistance, and the opportunity to share in your work. This was as much your project as it was mine. Good luck with your current work and future studies.Finally, to my family, thank you for keeping me grounded and providing an unwavering source of support that I could depend upon. To my friends, thank you for the much needed laughter (even if it did verge on hysterical at times) and for providing numerous outlets for me to blow off some steam. This year has been the most challenging of my life so far, but also the most enjoyable.ii Abstract Due to their unique properties, dense viscoplastic fluid-particle suspensions have significant potential to be used in a number of industrial applications, from hydraulic fracturing to pipeline particulate transport. However, the rheological properties of these suspensions are currently poorly understood, with no comprehensive modelling strategies existing to predict their behaviour. Numerous numerical difficulties arise when attempting to model viscoplastic suspensions, a key cause of which is the presence of a fluid yield-stress. Currently, a number of explicit regularisation techniques are used to approximate the yieldstress, the inherent numerical inaccuracies of which are quite often given little attention.A coupled LBM-DEM numerical approach presents an excellent solution to modelling the bulk movement of particles within suspensions, and has been successfully applied to the modelling of dense Newtonian suspensions. Of key importance, however, the LBM allows for an implicit regularisation of the yield-stress, in which the constitutive Bingham fluid is solved without the need for the approximations of current explicit solvers. Hence, it was hypothesised that the implicitly-regularised (IR) model is superior for the modelling of viscoplastic fluids, a claim which, if true, would lead to the eventual direct numerical simulation of dense viscoplastic particle suspensions. Consequently, the aim of the investigation was to evaluate the performance of the IR LBM-DEM in viscoplastic fluid-particle coupling, compared to that of the explicit Papanastasiou-regularised (PR) two-relaxation-time (TRT) method.To achieve this, five numerical testing models were developed, based off existing numerical and experimental studies of simplified benchmark problems in the literature. These 2-D and 3-D simulations provided a basis for the validation of the IR model, as well as a means for performance comparisons between the IR and TRT models.In initial 2-D simula...

show abstract

Motion and stability of cones in a yield stress fluid

2014

View full text Add to dashboard Cite

This experimental study focuses on the creeping flow of a shear thinning yield stress fluid around conical obstacles. The flow has been analyzed in steady state and with adherence conditions. Firstly, the influences of the cone apex angle and of the Oldroyd number, that is the ratio between plastic and viscous effects, on the drag coefficient have been analyzed. Correlations have been proposed to model the evolution of this coefficient as a function of these two parameters. The analysis provides a new alternative for measuring the yield stress. Then, the kinematic fields around the cones have been analyzed. These fields enable to describe the rigid zones and the sheared zone developing around the lateral edge of the cones as a function of the cone apex angle. Moreover, the wall shear stresses estimated from the particle image velocimetry measurements have enabled to quantify the contribution of the lateral drag force in the drag force. V C 2014 American Institute of Chemical Engineers AIChE J, 61: [709][710][711][712][713][714][715][716][717] 2015

show abstract

A fictitious domain method for dynamic simulation of particle sedimentation in Bingham fluids

Cited by 78 publications

References 34 publications

PeliGRIFF, a parallel DEM-DLM/FD direct numerical simulation tool for 3D particulate flows

PeliGRIFF, a parallel DEM-DLM/FD direct numerical simulation tool for 3D particulate flows

Performance evaluation of an implicitly-regularised model for viscoplastic particle suspensions

Motion and stability of cones in a yield stress fluid

Contact Info

Product

Resources

About