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