adding dispersion force term in UEqn leads to oscillations of velocity field
1 Attachment(s)
Hello Foamers,
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 Code:
#include "FAdh.H" Maybe someone of you has a clue about what part of my implementation is causing this physically incorrect velocity field? Kind regards, Fabian |
Solution
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 !):
Code:
#include "FAdh.H" |
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. http://www.cfd-online.com/Forums/ope...tml#post495673 |
All times are GMT -4. The time now is 05:57. |