
[Sponsors] 
October 2, 2013, 13:09 
COMSOL  turbulent flow over a plate

#1 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi guys,
I calculated the cooling of a plate by a turbulent forced convection analytically. Now I want to check my results with COMSOL 4.3 but I have some troubles, maybe someone would be so kind and help me. My plate is 0.075 m long, to keep the plate at +4 °C i need a power of ~130 W. The cooling fluid is air and has a temperature of 30 °C with a "inflow"  velocity of 45 m/s. I tried to model my plate in COMSOL with a rectangle ( 0.075 m x 0.005 m ), this plate is in the upper right corner of an "larger" rectangle. The inlet is on the left. The material of my plate is copper. My first problem is how to fit the power, because for my analytic approach i have just a surface (~0.017 m^2) but no volume, which is necessary for a heat source with a unit of W/m^3. When I try to run the simulation with "Total Power" equal my required 130 W, the result is far away from my analytical one. I'm thankful for any advice, I'm pretty new at COMSOL and an electrical engineer so I'm not really familiar with themal stuff ;). BR Last edited by stenzer; October 7, 2013 at 12:33. 

October 2, 2013, 13:26 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,915
Rep Power: 34 
If I understand correctly your question, you should work by setting the heat flux in W/m^2:
q =  k grad T so that, by fixing q, you set a nonhomogeneous Neumann BC.s dT/dnwall = q/k 

October 2, 2013, 14:34 

#3 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi,
I don't know if I understand you correctly. I calculated the power which is needed to keep my plate at a certain temperature (+4°C). So in my opinion I have to model the plate as a heat source and not as heat flux ?!?!???? BR 

October 2, 2013, 17:32 

#4  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,915
Rep Power: 34 
Quote:
Could you explain your problem by means of a sketch? 

October 3, 2013, 04:08 

#5 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi,
here is the link for a sketch, I hope you can read everything. https://www.dropbox.com/s/ggs1e7wrfvk2y5c/Cooling.jpg thx 

October 3, 2013, 04:24 

#6 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,915
Rep Power: 34 
it is a 2D model? therefore, your power for m^3 must be considered for 1m of unit lenght in the third direction (outward the sketch plane) in such a way you have a volume.
I think you could model a mixed fluidcopper problem, so that the heat exchange at wall is due to the heat flux. In the copper region you can solve a standard diffusion equation for the temperature which is governed by a parabolic equation like dT/dt = k Lap T 

October 7, 2013, 12:10 

#7 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi,
I did it like you suggested. But it doesn't fit with my "analytical" result. I already checked my calculations, but in my opinoin they a are right. At +4°C the kinematic viscosity (nu) of airis v =13.85e6 m^2/s, Pr = 0.7116 and the heat conductivity k = 24.6e3 W/mK. Thus my Reynolds number is: Re = 243682 Because I want to investigate a plate with a plunt edge I assumed a complete turbulent flow, thus I can calculate my Nusselt number by Nu = [0.037 * Re^(0.8) * Pr]/[1 + 12.7 * sqrt{ 0.037 * Re^0.2} * (Pr^(2/3)  1)] = 626.795 Now I can determine my heat transfere coefficient h and the needed power P to keep the surface at +4°C. h = 205.5 W/(m^2 K) P = 0.176 m^2 * h * (4°C + 30°C) = 133 W I think I made no mistake, but maybe someone can find a mistake or give me any advise how to model that in COMSOL. I'm useing Conjugate Heat Transfer and k  epsilon Turbulent Model. thx 

October 7, 2013, 12:19 

#8  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,915
Rep Power: 34 
Quote:
well, standard results of a flow over a flat plate established that a fully turbulent flow should be achieved for Re number greater than yours ... I think you are still in a transitional region and I am doubtful a keps model can manage succesfull such a condition 

October 7, 2013, 12:32 

#9 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi,
thank you for your fast reply. I already thougt about that, because the turbulent flow can be assumed for Reynolds numbers between 5*10^5 and 10^7. To be honest I thought that would be modeled by setting the boundary/wall condition as Wall function . I will calculate the model for a laminar flow and check if the result is correct. But my problem is still that I have to know how a plate with a plunt edge would act. Is there a possibility to get results for such a model??? 

October 7, 2013, 12:35 

#10 
Member
Serge A. Suchkov
Join Date: Oct 2011
Location: Moscow, Russia
Posts: 74
Blog Entries: 5
Rep Power: 7 
The use of the wall function (as shown on your sketch) for keps turbulence model assumes that y+ on the plate is in the range of 30200, but for heat transfer problems the recommended value of y+ ~ 1 and no wall function (for reasonable results)
__________________
OpenHyperFLOW2D Project 

October 8, 2013, 05:34 

#11 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Hi SergeAS,
I use conjugate heat transfer, so I think that should be considered by the model automatically!?!?!??? My idea is to calculate and simulate the whole model for a longer plate, thus I get a larger Reeynolds number. And because of the Similitude the result for my smaller plate should be correct. Can I do it that way??? BR 

October 8, 2013, 05:59 

#12  
Member
Serge A. Suchkov
Join Date: Oct 2011
Location: Moscow, Russia
Posts: 74
Blog Entries: 5
Rep Power: 7 
Quote:
PS: Take a look http://www.comsol.com/blogs/whichtu...dapplication/ Quote:
__________________
OpenHyperFLOW2D Project Last edited by SergeAS; October 8, 2013 at 10:11. 

October 10, 2013, 10:19 

#13 
New Member
Join Date: Oct 2013
Posts: 7
Rep Power: 5 
Thank you for the link, I will try it with Low Reynolds number kepsilon.
thx 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
problem in turbulent flow over flat plate with stagnation point  mb.pejvak  Main CFD Forum  0  September 12, 2012 21:25 
Turbulent Flow in Straight Square Duct  kep Validation  mcclud  STARCCM+  2  August 7, 2012 16:56 
laminar and turbulent flow in one simulation  msna  FLUENT  0  January 27, 2007 18:35 
Flow passing a porous plate.  Hao Yu  Main CFD Forum  2  October 19, 2000 01:52 
Turbulent flow over flat plate  Narayanan  Main CFD Forum  1  May 8, 2000 07:03 