[Nek5000-users] 2D axisymmetric problems
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Dec 17 09:30:11 CST 2010
Generally we store azimuthal velocity as temperature (which I think
you've done in the past...), there is a flag that needs to be set,
however, for the Helmholtz operator. The issue is how to deal with
this in a more general way, so that the user can control it ... and
then, what happens if the user also wants heat transfer???
Our current plan is that if you also want heat transfer then temperature
will be stored in PS1.
As I mentioned, Aleks will get you the details... I think we should be
able to have it as a flag from the .usr file so you don't need to
touch the source. The caveat will be that the user needs to recognize
that a separate set of bcs is required for "temperature", which will be
bcs on the azimuthal velocity.
Paul
On Fri, 17 Dec 2010, nek5000-users at lists.mcs.anl.gov wrote:
> Paul,
>
> Great! Katie mentioned a solution to the "ERROR" message -- I think the only remaining thing is the representation of the azimuthal velocity.
>
> --Mike
>
> On Dec 17, 2010, at 8:23 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>>
>> Hi Mike,
>>
>> Aleks is looking into this and will get you an answer today.
>>
>> 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
>
> _______________________________________________
> 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