[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [ESPResSo-users] Tube Flow Lattice Boltzmann
From: |
Georg Rempfer |
Subject: |
Re: [ESPResSo-users] Tube Flow Lattice Boltzmann |
Date: |
Thu, 18 Apr 2013 20:15:26 +0200 |
Hello Wolfgang,
here is a script that simulates EOF in a cylinder. It definitely
worked 1.5 years ago. You should be able to extract the necessary
parts from it to compare to yours.
Greetings,
Georg
2013/4/18 Wolfgang Riefler <address@hidden>:
> Hello all,
>
> I am trying to simulate a flow within a cylindrical tube, using the Lattice
> Boltzmann method.
>
> I am using the following script:
>
>
> setmd box_l 54 18 18
> setmd time_step 0.005
> setmd skin 0.1
>
> lbfluid agrid 1 dens 1 visc 10 tau 0.005 ext_force 0.04
> thermostat lb 0.0
>
> lbboundary cylinder center 27.0 9.0 9.0 axis 1.0 0.0 0.0 radius 8.26 length
> 27 direction -1
>
> integrate x
>
> According to the user guide this should define a cylindircal surface with
> symmetry axis parallel to the x-axis, where all nodes outside the cylinder
> are defined as boundary nodes.
>
> However the fluid stays inert. If I change the direction to 1, I can see the
> fluid flow outside the cylinder in x-direction as expected.
>
> What am I doing wrong?
>
> Thanks in advance,
>
> Wolfgang
>
run.tcl
Description: Tcl script