Re: [ESPResSo-users] lbboundary

[Axel Arnold]

On 22.07.14 14:58, roya moghaddasi wrote:

Hello dear everybody,

I have a cylinder in my script like this:

constraint cylinder center $c $c $c axis 0 0 1 radius $r length $box_length direction -1 type 4 penetrable 0

All of my particles are within this cylinder,

I want to have hydrodynamics, LB, in my system, so I've added the following to my script:

kill_particle_motion
thermostat off
thermostat lb 1.0
lbfluid dens 0.864  visc 3.0  agrid 1.0  tau 0.01
lbfluid friction 20.0
integrate 100000

Should I use lbboundary in my script or not?

If you also want the fluid to stay within the cylinder,  you need a lbboundary with the same parameters (except for the last two parameter that determine the interactions).

Axel

-- 
JP Dr. Axel Arnold
ICP, Universität Stuttgart
Allmandring 3
70569 Stuttgart, Germany
Email: address@hidden
Tel: +49 711 685 67609