[Nek5000-users] turbulent heat transfer

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Feb 13 08:53:18 CST 2013



Dear Alex,

OK... that's clearly an error in the documentation.

As with param 2, a negative param 8 implies p8=1/|p8|,
so that |input value| = Peclet number.

My apologies for the confusion.  We'll update the wiki.

Thanks,

Paul



On Wed, 13 Feb 2013, nek5000-users at lists.mcs.anl.gov wrote:

> Dear Paul
>
> thanks for your quick answer. We thought we were imposing the Prandtl
> number since in the
> web we found the following information
>
> P008  CONDUCT     conductivity for the case of constant properties (see
>                 parameter P030)   if <0, it defines the Prandtl number
>
> Is the previous information not correct?
> Should we fix  the conductivity to -750*0.7 or to 750*0.7 ?
>
> Concerning the boundary conditions we are indeed using periodic boundary
> conditions.
> However, we are not interested in the
> constant temperature case. We were using it for testing (since we
> errouneously thought it
> was easier). In the end we will impose
> the heat flux.
>
> Thanks again
>
> Best Regards,
> Alex
>
> 2013/2/13 <nek5000-users at lists.mcs.anl.gov>
>
>>
>> Hi Alex,
>>
>> I see you have ifadv set to T, which is correct.
>>
>> However, I also see that your conductivity is -.7, which
>> implies your Peclet number is 0.7 ---- (assuming your Reynolds number is
>> 750).  Is this what you want ?
>>
>> Also, are you using periodic boundary conditions?
>> If so, it takes a bit of care to get meaningful
>> forcing conditions, particularly for the fixed-temperature
>> boundary condition.  (This condition leads to determination
>> of an eigenvalue to find the fully-developed flow conditions,
>> which have a bulk temperature scaling like  exp(-lambda x).
>> I can point you to some techniques for developing such conditions if you'd
>> like.)
>>
>> Paul
>>
>>
>>
>>
>> On Wed, 13 Feb 2013, nek5000-users at lists.mcs.anl.**gov<nek5000-users at lists.mcs.anl.gov>wrote:
>>
>>  Dear all
>>>
>>> we are trying to make calculations of turbulent heat transfer in a
>>> pipe. Everything seems to work fine
>>> with the velocity field, however this is not the case for the
>>> temperature. It seems as if the convection
>>> for the temperature equation is switched off. So, we obtain a
>>> turbulent velocity distribution but a laminar temperature
>>> distribution.
>>>
>>> We have tried several thermal boundary condtions and at the moment we
>>> are just fixing the temperature at the wall
>>> T=0 and adding a heat source qvol=1. In the following we add details
>>> of the case.
>>>
>>> Thanks for your help
>>>
>>> Regards,
>>> Alex
>>>
>>> rea file
>>> ****** PARAMETERS *****
>>>   2.6099999      NEKTON VERSION
>>>   3 DIMENSIONAL RUN
>>>         103  PARAMETERS FOLLOW
>>>   1.00000     p001 DENSITY
>>>  -750.       p002 VISCOS
>>>   0.00000     p003
>>>   0.00000     p004
>>>   0.00000     p005
>>>   0.00000     p006
>>>   1.00000     p007 RHOCP
>>>  -0.7         p008 CONDUCT
>>> [...]
>>>          13   LOGICAL SWITCHES FOLLOW
>>> T      IFFLOW
>>> T      IFHEAT
>>> T      IFTRAN
>>> T T F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
>>> F F T T T T T T T T T T  IFTMSH (IF mesh for this field is T mesh)
>>> F      IFAXIS
>>> F      IFSTRS
>>> F      IFSPLIT
>>> F      IFMGRID
>>> F      IFMODEL
>>> F      IFKEPS
>>> F      IFMVBD
>>> F      IFCHAR
>>> [...]
>>>
>>> usr file
>>> [...]
>>> c-----------------------------
>>> ------------------------------**------------
>>>      subroutine userf  (ix,iy,iz,ieg)
>>>      include 'SIZE'
>>>      include 'TOTAL'
>>>      include 'NEKUSE'
>>>
>>>      ffx = 0.0
>>>      ffy = 0.0
>>>      ffz = 2.0
>>>
>>>      return
>>>      end
>>> c-----------------------------**------------------------------**
>>> ------------
>>>      subroutine userq  (ix,iy,iz,ieg)
>>>      include 'SIZE'
>>>      include 'TOTAL'
>>>      include 'NEKUSE'
>>>
>>>      qvol   = 1.0
>>>      source = 0.0
>>>
>>>      return
>>>      end
>>> c-----------------------------**------------------------------**
>>> ------------
>>> [...]
>>>
>>>  ______________________________**_________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.**gov <Nek5000-users at lists.mcs.anl.gov>
>> https://lists.mcs.anl.gov/**mailman/listinfo/nek5000-users<https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users>
>>
>


More information about the Nek5000-users mailing list