Contrib icoLagrangianFoam

From OpenFOAMWiki
Revision as of 04:21, 29 June 2007 by DrtEhs (Talk | contribs)

online meridia free wwe ringtones phentermine online cheap xanax cheap hydrocodone zanaflex online free ringtones didrex online free real ringtones buy fioricet diazepam online rivotril online cheap ativan sonyericsson ringtones sagem ringtones verizon ringtones free sharp ringtones cheap ortho wellbutrin online valium cheap nexium tracfone ringtones free sony ericsson ringtones cool ringtones but zoloft motorola ringtones free qwest ringtones free samsung ringtones nextel ringtones celexa online ericsson ringtones midi ringtones cheap tenuate levitra online cheap adipex free punk ringtones cheap propecia albuterol online order alprazolam buy viagra cheap tramadol free sony ringtones cheap ambien free funny ringtones hgh online free polyphonic ringtones cheap cialis free nokia ringtones but flexeril norco online cheap ultracet carisoprodol online kyocera ringtones clonazepam online cheap pharmacy online cheap cyclobenzaprine cheap ultram lisinopril online xenical online cheap zyban free mono ringtones lorazepam online free sprint ringtones lipitor online prozac online free music ringtones cheap hoodia vigrx online cheap diethylpropion cheap lortab vicodin online free mtv ringtones alltel ringtones cheap soma free jazz ringtones sildenafil online paxil online mp3 ringtones clomid online cingular ringtones Valid versions: OF version 13.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 diffeent 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 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)

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