[Nek5000-users] 2D axisymmetric problems

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

Hi Mike,


A little correction to the previous email -- for the current version of axisymmetric setup with azimuthal velocity one needs to set 

       ifaziv = .true.

somewhere in .usr routines (including nek5_svn/example/vortex2/v2d.usr) since we do not have a mechanism for reading this flag from .rea file yet -- we will have it in place in one of the next updates...

Best,
Aleks



----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Friday, December 17, 2010 9:27:55 PM
Subject: Re: [Nek5000-users] 2D axisymmetric problems



Hi Mike,

We've updated the repo and have added an example "vortex2" 
that shows axisymmetric plus swirl.   I confess I've not
tried it with .re2 format yet, but for small problems .rea
should be ok.

The key is to add a logical "ifaziv" to be true for azimuthal
velocity, and then to add the Coriolis force terms to the .usr
file, per the example.   As always with axisymmetric problems,
you must have "A  " for the boundary condition on the x-axis
if y=0 is part of your domain boundary.

Please let us know if this meets your needs, or not...!

Cheers,

Paul



>
> On Tue, 14 Dec 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hello All.
>> 
>> I have a couple of questions regarding simulation of 2D axisymmetric 
>> problems where there's a nonzero azimuthal velocity component.
>> 
>> First, does the new .rea + .re2 format work with such problems?  I tried 
>> this with the "axi" problem in the examples directory, and things halted 
>> with
>>
>>    ERROR READING INITIAL CONDITION/DRIVE FORCE DATA
>>  ABORTING IN ROUTINE RDICDF.
>> 
>> It worked fine if I stuck with the older .rea format.
>> 
>> Second, when doing such problems a while back, we would trick nek for such 
>> problems by using temperature to represent the azimuthal velocity.  This 
>> approach required one to add appropriate terms into the .usr file.  Is this 
>> still the case, or is Nek's third velocity component now the azimuthal 
>> velocity?
>> 
>> Thanks.
>> --Mike
>>
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users