[Nek5000-users] turbulent heat transfer

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Feb 13 09:01:40 CST 2013


PS - sorry - I forgot to answer your question about
how to set p8.

I guess I would set p2=-750 and p8 = -(750*0.7)



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

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