HI Sandeep, I think Dave has implemented this type of feature but I do not know if it is in Master or a branch? I am cc’ing Dave for who is best able to comment. Cheers, Spencer. On 27 May 2018, at 14:25, Sandeep Pandey <spandey.ike@gmail.com<mailto:spandey.ike@gmail.com>> wrote: 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 _______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk<mailto:Nektar-users@imperial.ac.uk> https://mailman.ic.ac.uk/mailman/listinfo/nektar-users Spencer Sherwin FREng, FRAeS Head, Aerodynamics, Professor of Computational Fluid Mechanics, Department of Aeronautics, Imperial College London South Kensington Campus, London, SW7 2AZ, UK s.sherwin@imperial.ac.uk<mailto:s.sherwin@imperial.ac.uk> +44 (0)20 7594 5052 http://www.imperial.ac.uk/people/s.sherwin/