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