Re: [ESPResSo-users] Setting up Lattice Boltzmann

[Vincent Ustach]

Stefan,

>>What is your mass and your Diffusion coefficient? Most likely you'll find out that if you put in the numbers, both your viscosity and your friction will be incredibly high.<<

Do you mean what are the empirical mass and Diffusion coefficients of the particles I am simulating? Indeed for the parameters

density 1 viscosity 0.8 friction 10 agrid 1 tau 0.01

the spheres moved extremely slowly.

Is there a reference I can find for what the LB units are?

--Vincent Ustach
  University of California, Davis

On Tue, May 7, 2013 at 12:14 AM, Stefan Kesselheim <address@hidden> wrote:

Hi,

On May 7, 2013, at 5:43 AM, Vincent Ustach <address@hidden> wrote:

> Stefan,
>
> Thanks for your reply! I have the code running but of course I still have a lot of work to do. I did not want the discussion to get stale, however, so I am replying now with the info that I have.
>
> I am abandoning the DPD thermostat and using LB for a thermostat and for the HIs. I was not using inter_dpd. Sorry that I was not clear.

That means that you also plan to throw out the DPD particles? OK. That sounds fair. Have you tried the GPU version of LB? It is incredibly fast!

> I will look into the pore constraint, but I don't want to restrict myself to cylindrical and conical pores.

That makes sense. However it is probably better to use the bounce back method for walls. One thing that we always wanted to do, is making it possible to flag single nodes as boundary nodes, but this is not implemented at the moment. Though quite easy, both for GPU and CPU. However I'm rather busy with my thesis right now, and thus would prefer to do it later :-). You can do it quite easily.

> This is indeed a concentration gradient driven diffusion, however electrostatics will be involved later. Viscosity seems to be (0.8 Pa * s) *10^-3. What are the units for friction?

What is your mass and your Diffusion coefficient? Most likely you'll find out that if you put in the numbers, both your viscosity and your friction will be incredibly high.

> I am working on printing the fluid velocity field to an output file and then hopefully it will shed some light on if the wall of particles serves as an adequate LB boundary.

You can try. It depends on your system, if it is. This I can tell you for sure. And on you definition of "adequate" :-).

Please don't think that you can get LB for free (not computationally, but in terms of carefully thinking about your system). Just as in DPD you have to worry seriously about getting parameters right. I have spent quite some time learning about that and I'm quite sure I can get helpful advice now. But you need to know: Will you have "macroscopic" fluid motion in your system, or is just correlations that matter? Probably in both cases it is necessary to fit the "effective particle radius" 1/(6 pi eta D) to the other lengths scales in you system. But in the first case it is OK if the radius is small, while in the second case it is important to have a high resolution to get the near field right.

Don't worry about the discussion getting stale, we can revive it later :-).  My desk is also quite full.

Cheers and good luck
Stefan