Die Multiphase Particle in Cell Methode in Matlab

Translated title of the contribution: The Multiphase Particle in Cell Method in Matlab

Christian Witz

Research output: ThesisMaster's Thesis


This paper tries to implement the Multiphase Particle-In-Cell method in Matlab. The method was developed by Andrews and O'Rourke. It has the advantage of the Euler-Lagrange method, that different particle sizes can be simulated relatively simple. But at the same time this method reduces the computational effort for higher solid fractions because the solid pressure is not calculated for each individual particle, but for each cell. The created program solves one-dimensional problems and deals with different particle sizes and densities. It is possible to define the viscosity, density and flow velocity of the surrounding fluid. The particles are grouped into packages and the number of particles in each package can be specified. Wall collisions are simulated with a factor of restitution. It determines the fraction of the initial speed the particle has after the collision. The first test case deals with the free fall of a single particle. The difference in falling time between the analytical solution and the solution calculated by the program is only half a percent. The sedimentation of a pile of particles is the topic of the second test case. The result shows, that numerical diffusion appears. Particles with different diameters are used in the third test case. It can be concluded that the first results are positive, however a lot of effort is needed to implement correction algorithms for the numerical diffusion. The behaviour of the program when dealing with high solid fractions will be a subject to further research.
Translated title of the contributionThe Multiphase Particle in Cell Method in Matlab
Original languageGerman
Awarding Institution
  • Montanuniversität
  • Weiß, Christian, Supervisor (internal)
Award date26 Mar 2010
Publication statusPublished - 2010

Bibliographical note

embargoed until null


  • Multiphase Particle in Cell Matlab

Cite this