Contrib icoLagrangianFoam

From OpenFOAMWiki

Valid versions: OF version 13.png OF version 141.png OF version 15.png

1 Short description

A solver for incompressible flow and hard particles

2 Usage

Used like icoFoam but has an additional file cloudProperties for the specification of the particles and their injection.

2.1 Problems with 1.2

Note: the code of the solver has to modified to work with 1.2

This code works with the vanilla 1.2 executables. The problem is that not all the data of the particles is output. Applying this fix and recompiling OpenFOAM makes everything work.

2.2 Post-Processing

paraFoam can't display the particles, but paraview can. Look at Postprocessing of Lagrangian particles in the FAQ

3 Physics

The particle code is a heavily lobotomized version of stuff found in the dieselSpray classes. It features:

  • a simple random injector
  • a drag force model that is horrible and not very physical
  • particles can bounce from walls or die (switchable)
  • particles leave the system at in or outlet. All other boundary types are not treat correctly
  • the particles can add a source term to the moment equation of the gas (switchable)
  • there is no particle-particle interaction

This solver is not to be used for simulations that resemble the real world. It's just a demo.

4 Example cases

4.1 Driven cavity

Particles in a driven cavity

This case is a three-dimensional driven cavity. Particles are injected in the middle. They don't perturbate the gas velocity.

The state at the end of the simulation can be seen on the right:

  • Gas velocity is shown by a plane of glyphs
  • Particles are shown by balls scaled according to their diameter. The velocity of the particle is shown by an arrow

Heavy particles have a tendency to settle at the bottom.

4.2 Open box

Particles interacting with the velocity field
Particles falling with different speeds

In this case a box with walls on five sides and a pressure outlet at the top is calculated. There is no initial velocity of the gas. Particles are eject upward and they perturbate the flow field. Some of the particles (the heavier ones) leave the box at the top. The others fall to the bottom of the box.

On the left a picture from the start of the simulation: it shows

  • the particles and their current velocity
  • isosurfaces of the regions where the gas has a non-zero velocity

On the right a picture from the end of the simulation. Particles are falling with different terminal velocities depending on their size.

5 Download

The most up-to-date (but not very extensively tested) version of the sources (and the one working on 1.5) can be downloaded via subversion:

svn checkout svn://

A version working on 1.4.1 can be found at:

svn checkout svn://

The packages below only work with 1.3:

The solver

Cavity case

Ejector case

6 History

  • 2005-11-11: Initial upload
  • 2006-04-18: Changed to work with 1.3. Now also works in parallel (Thanks to Jeff Allen)
  • 2010-07-23: A replacement for this solver can now be found in the tutorials of the 1.5-dev. It uses the new lagrangian/intermediate-classes that were introduced in 1.5. Therefore maintaining the solver described here doesn't make sense anymore

--Bgschaid 10:35, 11 Nov 2005 (CET)