[Nek5000-users] How to add Coriolis force

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

Hi,praveen:
    Thank you for your explaination, it's really helpful. By the way, I have two more questions(sorry, I'm really newer to Nek):
     1/ Is there any way to map  a coarse mesh solution to a finer mesh solution?
     2/ When the "0 RESTART OPTIONS FOLLOW" is changed to '1 RESTART OPTIONS FOLLOW' with a specified field, is the average field been loaded in the starting? Because I do have some variables in the subroutine userchk, such as 'common /plane/  uavg_pl(ly1*lely)'







At 2014-05-31 10:04:11, nek5000-users at lists.mcs.anl.gov wrote:

vx is the solution array, see the include file trunk/nek/SOLN whereas ux is just a real number which is filled with the velocity at the current evaluation point. So you should use ux in userf.


Best
praveen



On Sat, May 31, 2014 at 3:47 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

Hi,Mike:
    Is there any difference by using ux and vx? As in the help, most frequently used is vx. However, when I use vx in the formulation, the compile turn out error.

Xianbei







At 2014-05-29 09:09:16, nek5000-users at lists.mcs.anl.gov wrote:

Hi,Mike
   Thank you very much for you help. Now I have a question about the built-in variables. As the help document is not so sufficient, there are so many variables not introduced, so how can I find what's the meaning of the variables? I'm sorry to trouble you for this simple question, but it's really important to me.


Thank you
Xianbei








在 2014-05-28 11:53:58，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
>
>_______________________________________________
>Nek5000-users mailing list
>Nek5000-users at lists.mcs.anl.gov
>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users



来自网易手机号码邮箱了解更多



来自网易手机号码邮箱了解更多

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users


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