
[Sponsors] 
August 16, 2010, 15:00 

#21  
Member
Join Date: Dec 2009
Location: Kanpur, India
Posts: 54
Rep Power: 16 
Quote:
For confirming that this would give me the variable heat loss. Though I had arrived at what I wanted by following changes. Here what I did that I wrote Tamb=(internalField(T)); which is variable. so as you said (TTamb) changes then wall should loss. so I am getting the variations. type groovyBC; gradientExpression "gradT"; fractionExpression "0"; variables "hcof=0.88;Cp0=520.3;rho0=0.8947;Tamb=(internalField(T)); alphaEff=1.2956*pow(10,4);heatFlux=hcof *(TTamb);gradT=heatFlux/(alphaEff * Cp0 * rho0);" timelines ( ); value nonuniform List<scalar> 3300 ( ); As you had a question about 0.5.That I was trying to use temperature BC also; which was 393 K (known). So I thought to make a Robin type of BC [ i.e T and heat loss (i.e. gradient )], so I used 0.5. Again a small question if you could answer please, should I use T and heat loss together or not. If I use them together then does it similar to just using heat loss without T ?. How these two changes makes difference in the system. thank you dinesh 

September 6, 2010, 01:02 

#22  
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 16 
Quote:
I am doing a similar heat transfer problem through a series of concentric ' annular walls'. May I ask you several things regarding your groovyBC? 1.) You specify your heat trf coeff as 0.88. How did you calculate this value? In a convection problem, the hc is a function of local Nusselt number; hence varies along the wall. 2.) Are the rho0 and cp0 properties of the solid or the fluid? 3.) What is alphaEff ? Is alphaEff the heat conductivity of the fluid? Thank you very much for your help. Best Regards, Stefano
__________________
Stefano Wahono Defence Science and Technology Organisation Propulsion Systems 

September 6, 2010, 16:42 

#23  
Member
Join Date: Dec 2009
Location: Kanpur, India
Posts: 54
Rep Power: 16 
Quote:
I use Nussult number Nu=5.62*10^3*(Ra^(0.443)); (i) ///this correlation is specific to my problem I found from  Martin B W and cohen H British, Jounal of applied Physics, Vol5 1954. These type of corelation are obtained form experimetnal studies. your problem may also have this if not I suggest that it can be calculated using Darcy formula. which is like as h=0.023 (K/De)*( (De*U*rho/mu)^0.8 ) *(Pr^0.4) (ii) //here K=thermal conductivity, D=hydrolic radius, rho =density, mu=kinamatic viscosity. Pr=prandtl number, U =max valocity in the sytem (which you shall have to guess). (ref. M. M ElWakil page no. 252 nuclear Power Engineering ) and then I calculate h_coff=Nu*K/De (iii)//well known I set Rayleigh number in my system by changing temperature difference (TT_wall) and characteristic length (L). and then Nu is calculated, which gives h_coff using eqn (iii). Sorry to say that I do not have answer for if local hcoff varies. My knowledge of CFD is limited. I did not think about this aspect, I had searched in this forum; I found above type of setting only. I shall try this aspect too. If I find then let you know soon, meanwhile you find in your way. (2) rho0 and cp0 are the Argon gas values. (3) alphaEff is total thermal diffusivity of Argon alphaEff= thermal diffusivity + turbulent thermal diffusivity. further I would like to add this I have made following changes which I have tested properly. ******************************* type groovyBC; valueExpression "T_wall"; gradientExpression "gradT"; fractionExpression "0"; variables "hcof=20;Cp0=520.3;T_wall=460;rho0=0.9764;alphaEff =0.0020;heatFlux=hcof*(TT_wall);gradT=heatFlux/(alphaEff * Cp0 * rho0);" timelines (); **************************** here if heat tranfer is considered then fractionExpression "0"; if temp. + heat transfer is considered then fractionExpression "0.5"; hope this helps I would be happy to answer you questions if it falls in domain of knowledge. thank you dinesh 

September 6, 2010, 21:35 

#24 
Member
Stefano Wahono
Join Date: Aug 2010
Location: Melbourne, Australia
Posts: 42
Rep Power: 16 
Hi Dinesh,
Wow thanks for the quick reply! Yes, the Nusselt number is dependant of both the momentum and thermal boundary layers. But, most people indeed just use correlations to obtain a 'lump' Nusselt number. However, from your hceqn, U, rho and mu are local informations in the flow field. Do you think it'll be better to read it off the solution (i.e. use internalField(U), internalField(rho)) at every time step (or iteration)? If the wall is 'twosided' as follows: > hot flow, T_hot= 600 C, MACH =0.6 ==== THIN WALL ================= > cold flow, T_cold=50 C, MACH 0.1 Do you think your BC is applicable to this scenario? In this case, what is the T_wall (since the wall is 'twosided'; and varies along the wall depending on the temp of the gas on both sides) ? Thank you very much again for your discussion and help. It seems like OF is missing an important (and seemingly trivial) feature here (which has been taken for granted when using commercial package, i.e. FLUENT) 

September 9, 2010, 16:54 

#25  
Member
Join Date: Dec 2009
Location: Kanpur, India
Posts: 54
Rep Power: 16 
Quote:
I think that the h_coff written in the groovyBC is just initialization. As the solution advances; the wall heat transfer changes. If one need to h_coff; he can calculate from later time solutions. But I am firm that the h_coff written under the groovyBC does not remain fixed; it has to be calculated form different temperature values in later solutions. I my opinion how would I proceed 1) use various values of (TT_wall) in any time dir. 2) use this value in formula of Rayleigh number Ra =g*beta*(TT_wall)*L^3*Pr/nu^2; this gives various Rayleigh numbers. 3) then use these Ra's to calculate various Nusslt numbers, say, using Nu=5.62*10^3*(Ra^(0.443)) or whatever applicable for ones case. Note that these would be local Nusselt numbers. 4) then these Nusselt numbers can give various h_coffs at that time dir. Still I am wondering how to decide the physical locations of these values. So as soon as I'll find, post it. So in my opinion that groovyBC is useful. I again repeat that the confusion was that we were thinking that it is fixed throughout, which is not true. That was just an initialization, h_coffs has to be calculated from new solutions. please Let me know if it convince you. thank you 

May 10, 2011, 04:26 

#26  
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 16 
Quote:
how did you manage to connect the alphaEff with nuSgs? i'm desperately trying to get a modified pisoFoam solver (energy equation with alphaEff included) to work, that uses dynSmagorinsky turbulence model. can you show me the c++ source of wallHeatFluxLES? best regards 

May 24, 2011, 23:11 

#27 
New Member
Oky Andytya
Join Date: Nov 2010
Posts: 26
Rep Power: 15 
Hi everyone,
I need help, how to get the value of convection coefficient [h] from Fluent directly? Thank you. Regrads, Oky 

Tags 
wallheatflux 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
equilibrium constant in fluent  ashishpande  FLUENT  1  May 30, 2014 05:43 
Chemical properties that are used in Fluent  kimej  FLUENT  4  August 19, 2011 11:29 
Restart for Solidification  atulverma  FLOW3D  3  May 6, 2009 17:40 
v4 restart bug  optima prime  Siemens  3  February 2, 2009 22:51 
Really big problems when creating own materials da  Jop  FLUENT  0  June 28, 2007 11:15 