[Nek5000-users] 2D axisymmetric problems

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

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
>>