[Nek5000-users] Fixed mass flux

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

Hi Fred,

in general there are two different ways to do that. 
Anyway just have a look how we do it in the turbChannel example.

The general idea is the following:
- call set_forcing() in userchk()  (e.g. call set_forcing(ffx_,xv,1)
- add the acceleration term ffx_ to the momentum eqn. (e.g. set ffx = ffx_ in userf())

Please note that the variable ffx_ has to be visible for userf() and userchk() (e.g. use a common block).
The initial acceleration can be specified by param(71) .

#Stefan 


On Dec 11, 2009, at 2:12 PM, nek5000-users at lists.mcs.anl.gov wrote:

> Hello,
> 
> currently I'm trying to run some simple cases to getting used to the code. But unfortunately I didn't understand how it's possible to predefine a fixed mass flux. Do I just have to enter a fixed value for Parameter 54 resp. 55 in the rea-File? Maybe this question is a little stupid, but whatever I tried, the velocity in the result file .fld* was more or less equal zero..
> 
> Best regards,
> Fred
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users