[Nek5000-users] LES and heat transfer in channel

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Sun Feb 21 22:12:41 CST 2010


Hi,

thanks for the answers.
Before considering DNS calculations, I'd like to look into the 
Smagorinski model. Extending the turbulent channel example in the nekton 
source, is the following the correct way to set user defined properties 
for the energy equation?

       subroutine uservp (ix,iy,iz,ieg)
       include 'SIZE'
       include 'TOTAL'
       include 'NEKUSE'

       common /cdsmag/ ediff(lx1,ly1,lz1,lelv)
      $              , thermdiff(lx1,ly1,lz1,lelv)

       ie     = gllel(ieg)
       udiff  = ediff(ix,iy,iz,ie)
       utrans = 1.
       if (ifield.eq.2) then
         udiff  = thermdiff(ix,iy,iz,ie)
         utrans = 1.
       endif

       return
       end

To answer some previous questions, the excess heat is subtracted from 
the inflow temperature (via an energy balance between the periodic 
sections).
The flow is driven with a constant flow rate, much like in the jet example.

Markus


nek5000-users at lists.mcs.anl.gov wrote:
> Hi Markus,
> 
> it's hard if something is wrong with your setup or if you're running into other issues.
> 
> I would start with a simple DNS of turbulent heat transfer in a channel e.g. similar to Hawamura et. al Journal of Head and Fluid Flow, 1998.
> 
> Let's say you do the following simulation:
> Re_tau = 180
> Pr = 0.71
> Domain: (12.8/2/6.4) where you're reference length is the channel half-height
> Resolution: 24x12x24 elements using N=9 (this should be pretty well resolved)
> Thermal BC: uniform heat-flux heating
> 
> - How do the Nek results compare to the results of Hawamura?
> 
> I think that's a good way of testing your setup and after doing this exercise you're sure that everything works correct (assuming you get the right answers). 
> Then the next step is to use a much coarser resolution together with a SGS-model (e.g. the simple filtering and/or dynamic Smagorinski).
> 
> - How do the LES results compare to the DNS results?   
> 
> I would be not too surprised if you cannot get good results in your LES. It's well known that turbulent passive scalar mixing behaves quite different and a flow SGS model will not do a very good job in this case.
> 
> 
> Cheers,
> Stefan
> 
> 
> On Feb 9, 2010, at 4:34 AM, nek5000-users at lists.mcs.anl.gov wrote:
> 
>> Hi,
>>
>> I am running a channel LES simulation (periodic in streamwise and spanwise, length=width=9*channel height) in nek with the following parameters:
>> -LES model: Filter last 3 modes, 5% filter weight; lx1=14; ld1=20
>> -Piece of the rea file
>> "
>>   1.00000         DENSITY
>>   0.31250E-04     VISCOS
>>   1.00000         RHOCP
>>   1.00000         CONDUCT
>>   4.00000           p99=3 ==> dealiasing turned on
>> 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
>> T      IFCHAR
>> "
>> -Reynolds number based on 2*channel height: 12800
>> -Constant heat flux on walls
>> -The streamwise boundary is a recycling one, I set it up similar to the turbJet example, but I added temperature recycling. To avoid "over-heating", I subtract the added thermal energy from the temperature when recycling
>> -Some material parameters are set in the .usr file:
>> "
>> param(8)   = param(2)/0.71 ! Prandt=0.71
>> cpfld(2,1) = param(8)      ! conductivity
>> "
>>
>> When I compare the spanwise averaged and time averaged Nusselt number on a wall, it is fairly constant with streamwise direction (which is good), but about 50% higher then the experimental one from the correlation
>> Nu=0.022*Re^0.8*Pr^0.5
>> This might be because I am only 2.5 flow throughs away from the initial  condition, but I wanted to make sure that everything is right before putting in more computing hours.
>> Does the above look OK? Is there another explanation for the heat transfer coefficient being too high?
>>
>> Thanks,
>> Markus
>> _______________________________________________
>> 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