|
From: | Wolfgang Riefler |
Subject: | [ESPResSo-users] Lattice Boltzmann External Force |
Date: | Thu, 9 Feb 2012 10:49:51 +0000 |
j[i] += 0.5*lbpar.ext_force[i], where i=0,1,2 and j is the local momentum density.
Comparing this to the theory, I think it should be
j[i] += 0.5*lbpar.ext_force[i]*lbpar.tau,
i.e. the force should be scaled by the time step tau.
Implementing this in my source code indeed gave me the right solution for a grid size of a=1. Changing this to a=0.5 again caused an offset. Could it be possible that the force has to be scaled not only by time, but also by the grid size? I have noticed similar scalings in older versions of Espresso (Version 2.1.2).
On another note, is there any progress, or is it planned to allow the user to set the node populations manually (like e.g. the velocity)? I have noticed that the declaration of the corresponding function is already implemented, but the function is empty.
Thank you very much for your troubles in advance,
best,
Wolfgang
[Prev in Thread] | Current Thread | [Next in Thread] |