[Nek5000-users] 2D axisymmetric problems

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Dec 20 13:41:05 CST 2010


Paul,  Aleks,

Thanks!  All seems to be working great.  Haven't tried with .re2, but .rea works fine.  Validated the code against the numerical closed-cell Taylor-Couette results of Czarny et al (Physics of Fluids 15, 2003).

--Mike


On Dec 19, 2010, at 3:44 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> 
> 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
>> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users




More information about the Nek5000-users mailing list