[Nek5000-users] 2D axisymmetric problems

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

Axisymmetry with azimuthal velocity (swirl) should now work if you add

T IFAZIV

to the .rea file in addition to ifaxis.   You need to also set
IFHEAT to be true, since the swirl component is stored in the T()
array, and you need boundary condtions.

An example is now in    examples/vortex2

Hopefully I've properly checked everything into the repo.

Paul


On Sun, 19 Dec 2010, Aleksandr Obabko wrote:

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