Direct simulation of fluid particle motions. (English) Zbl 0754.76054

The authors develop a package that simulates the unsteady two-dimensional solid-liquid two-phase flows using the Navier-Stokes equations for the liquid and Newton’s equations of motion for the solid particles. The former are solved through a finite-element formulation and the latter through an explicit-implicit scheme. The authors show that the simplest fully explicit scheme to update the particle motion using Newton’s equations is unstable. To correct this instability they propose and implement the explicit-implicit scheme in which, at each time step, the positions of the particles are updated explicitly, the computational domain is remeshed, the solution at the previous time is mapped onto the new mesh and the Navier-Stokes equations and the discretized Newton’s equations are solved iteratively on the new mesh. The procedure is illustrated through its application to sedimentation of one and two discs (cylinders) in a vertical two-dimensional channel of finite width, and it is shown that it is able to reproduce successfully the effects of the vortex shedding on the motions of the cylinders and the drafting-kissing and tumbling scenario which is dominant in two-dimensional beds of spheres.
Reviewer: D.Jou (Ballaterra)


76M10 Finite element methods applied to problems in fluid mechanics
76T99 Multiphase and multicomponent flows
76A99 Foundations, constitutive equations, rheology, hydrodynamical models of non-fluid phenomena
76D05 Navier-Stokes equations for incompressible viscous fluids
Full Text: DOI


[1] Avalosse, Th., Couniot, A., and Crochet, M.J. (1990), Finite element simulation of compounding in domains of arbitrary shape, Private communication.
[2] Batchelor, G.K. (1967),An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge. · Zbl 0152.44402
[3] Couniot, A., Dheur, L., Hansen, O., and Dupret, F. (1989), A finite element method for simulating injection molding of thermoplastics,Proc. NUMIFORM ’89 Conf., Fort Collins, 26-30 June 1989, pp. 235-241, Balkema, Rotterdam.
[4] Crochet, M.J., Debaut, B., Keunings, R., and Marchal, J.M. (1991), A finite element program for modelling processes, inCAE for Polymer Processing: Application in Extrusion and Other Continuous Processes, K.T. O’Brian, ed., Hanser Verlag, Munchen.
[5] Cuthill, E., and McKee, J.M. (1969), Reducing the bandwidth of sparse symmetric matrices,Proc. 24th Nat. Conf. of the Association for Computing Machinery, p. 69
[6] Fortes, A.F., Joseph, D.D., and Lundgren, T.S. (1987), Nonlinear mechanics of fluidization of beds of spherical particles,J. Fluid Mech.,177 467-483.
[7] Hecht, F., and Saltel, E. (1989), Emc2 un logiciel d’√©dition de maillages et de contours bidimensionels, Rapport, INRIA.
[8] Milne-Thomson, L.M. (1960),Theoretical Hydrodynamics, Macmillan, New York. · Zbl 0089.42601
[9] Swanson, W.M. (1970),Fluid Mechanics, Holt, Rinehart and Winston, New York, Section 11.5.
[10] Tezduyar, T.E., Liou, J., and Behr, M. (1990a), A new strategy for finite element computations involving moving boundaries and interfaces?the DSD/ST procedure: I. The concept and the preliminary numerical tests. Supercomputer Institute Research Report UMSI 90/169, University of Minnesota, to appear inComput. Methods Appl. Mech. Engrg. · Zbl 0745.76044
[11] Tezduyar, T.E., Liou, J., Behr, M., and Mittal, S. (1990b), A new strategy for finite element computations involving moving boundaries and interfaces?the DSD/ST procedure: II. Computation of free-surface flows, two-liquid flows, and flows with drafting cylinders. Supercomputer Institute Research Report UMSI 90/170, University of Minnesota, to appear inComput. Methods Appl. Mech. Engrg. · Zbl 0745.76045
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.