[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