Re: [ESPResSo-devel] [ESPResSo-users] Simulation of fluid using LB-method

[Ulf Schiller]

On 01/30/2014 02:39 PM, Ulf Schiller wrote:
> On 01/30/2014 02:37 PM, Georg Rempfer wrote:
> > About the "pow(lbpar.__agrid,4)*tau*tau" issue. That has puzzled us
> > several times as well, but we are very confident that it is correct. I
> > can not remember the specific reason (too long ago) but I remember that
> > Owen checked that and Joost and I checked again, very thoroughly. I
> > think if you actually track down the values from input to application it
> > works out.
>
> Possible, but just to make sure you don't get the right result for the
> wrong reason: The units are incorrect here, so there must be a canceling
> factor elsewhere.

Which might be here...

    p_u[0] = j[0]/rho/tau/lbpar.agrid;
    p_u[1] = j[1]/rho/tau/lbpar.agrid;
    p_u[2] = j[2]/rho/tau/lbpar.agrid;

--
Dr. Ulf D. Schiller                        Building 04.16, Room 3006
Institute of Complex Systems (ICS-2)       Phone:   +49 2461 61-6144
Forschungszentrum Jülich, Germany          Fax:     +49 2461 61-3180

smime.p7s
Description: S/MIME Cryptographic Signature