adding dispersion force term in UEqn leads to oscillations of velocity field
I'm trying to implement an additional force in my solver, to represent dispersion forces between a wall and the fluid particles near the wall. The modelling approach for the magnitude of this force is as follows:
f = -3A / (delta^4)
where delta is the normal distance to the wall.
For my implementation I use a steady state SIMPLE procedure
Maybe someone of you has a clue about what part of my implementation is causing this physically incorrect velocity field?
It seems, that I forgot about calculating an additional flux. Now my implementation is based on the modelling of the gravity term in this wiki-entry. In case someone is interested in the final implementation (comments are welcome !):
In case the force is a vectorField, How could that be implemented?
I am trying to include, OF2.3.0, a vector force calculated over cell centers.
|All times are GMT -4. The time now is 13:50.|