Hello all,
I am a newbie to Nektar++ and I validated turbulent channel flow DNS with the standard DNS database. For now, the flow was driven by the constant pressure gradient. As a next step, I want to drive flow with such forcing so that I have constant
mass flux at every time step. The forcing should take the following form to achieve the constant mass flux:
dpdx(time)= dpdx(time-1)
-[ 0.3*(massFlowRef
(time=0)
- 2*massFlow
(time)
+ massFlow(time-1)
) / (Ly*Lz*dt)]
where; massflow should be volume integrated value.
Does anyone has already implemented such forcing? Or does anyone any idea? I am familiar with ./library/SolverUtils/Forcing, but any hint in the specific direction will be very useful to me.
Thanks and have a nice weekend!
Sandeep