[Nek5000-users] How to add Coriolis force

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

Xianbel,

You can add the Coriolis term in the userf routine:

      omega = something

c add coriolis forces
      ffx = + 2.0*omega*uy
      ffy = - 2.0*omega*ux
      ffz = 0.

Mike


On May 27, 2014, at 8:45 PM, <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>> <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>> wrote:

Dear Neks;
        I want to add Coriolis force to do a rotating simulation. As known, the Coriolis term is defined as 2*omega(cross product)U, so how can I add it to the N-S equation? Is there any example of similar case?

Thank you
Xianbei


来自网易手机号码邮箱了解更多<http://shouji.163.com/>
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users