[Nek5000-users] Constant mass flow rate in simulation

[nek5000-users at lists.mcs.anl.gov]

Hi Tony,

I would just set p55=1.0; p54=-1.0

This will fix < u_x > == 1.0

If you are using variable viscosity, this might be problematic, but for the
constant viscosity case this is the best way.

Be certain to also set ffx=ffy=ffz=0 in the .usr file, the the p54/p55 combination
sets the x-body force as needed to ensure the prescribed mass flux is realized.

Paul

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Sunday, August 31, 2014 5:57 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Constant mass flow rate in simulation

Dear Neks,

Recently I've been struggling with constant mass flow rate in my turbulent channel flow simulation. Based on the userchk subroutine in 'turbChannel' tutorial case, I also used torque_calc(1.0,x0,.false.,.false.) to calculate the wall shear and then got the mean wall shear. After that I set mean pressure gradient 'ffx' equal to the mean wall shear in order to keep constant mass flow rate. In this case, the flow rate started from 1 (nondimensionalised case, bulk mean equal to 1) and kept increasing. Although the increment was very small at each time step, it never stopped and the bulk mean velocity has increased to 1.002.

Is it appropriate to set constant mass flow rate by updating mean pressure gradient based on mean wall shear in Nek5000 or am I using the wrong calculation? Since the bulk mean is still increasing, the simulation is like acceleration channel flow.

Does anyone has any clue about how to set constant mass flow rate for turbulent channel flow in Nek5000? Any idea is appreciated and many thanks in advance.

Thanks and regards,
Tony

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140831/0fc955b0/attachment.html>