
[Sponsors] 
January 23, 2013, 22:47 
I really need some help on this problem

#1 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
Hi everyone,
I am trying to simulate the heat transfer performance of supercritical water flow in vertical upward round tube with 3.2 mm diameter. Here is a problem that bothers me for long time. The distribution patterns of the calculated wall temperature and heat transfer coefficient are not reasonable. The wall temperature increases too fast near the inlet of the tube and then increases almost linearly along the length of the tube. Since the coolant temperature increases almost lnearly along the length of the tube, the heat transfer coefficint, which is defined as h=q/△T, decrease very fast near the inlet of the tube and then almost keeps constant. I can not find the reason why the wall temperature increases so fast near the inlet of the tube. It should increase step by step, and the difference between the wall temperature and coolant temperature should increase with the increasing of the coolant temperature. I tried several other ke models. I got similar unreasonable distribution patterns. I attached a figure of the distributions of the temperature and heat transfer coefficient. Would some one please tell me what possible reasons can cause such problem? Thank you all, thank you for your kind advise 

January 24, 2013, 03:41 

#2 
Senior Member
siamak rahimi ardkapan
Join Date: Jul 2010
Location: Copenhagen, Denmark
Posts: 218
Rep Power: 11 
Can you explain what is the coolant here?
__________________
Good luck Siamak 

January 24, 2013, 03:44 

#3 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
I am sorry that I didn't describe clearly. The coolant is supercritical water. The temperature of coolant is higher than 390 oC


January 24, 2013, 12:16 

#4 
Senior Member
Pauli
Join Date: Mar 2009
Posts: 190
Rep Power: 10 
Not knowing your geometry, consider this more an interpretation of the graph you shared.
The graph indicates the entrance (length < 0.1m) has a much higher heat transfer coefficient than the remainder of the tube. This is driving the computed temperatures. The question is why does the entrance have higher heat transfer coefficient. From a flow physics standpoint, my initial assumption would be entrance effects (separation & mixing) cause additional mixing (varying velocity profile). At length > 0.1 the flow is fully developed. This provides a constant velocity profile & the h=q/dt relationship you expect. From a numerical calculation standpoint, you should be asking questions regarding how the heat transfer coefficient is calculated. Are you using wall functions or resolving the boundary layer? Is grid resolution adequate? Etc, etc. Computing heat transfer coefficients under pipe entrance conditions is not wall functions strength. 

January 24, 2013, 13:26 

#5 
New Member
Prashanth Avireddi
Join Date: Sep 2012
Location: Farmington Hills, MI
Posts: 10
Rep Power: 7 
Could you put up a schematic of the model you are working with, I have worked with numerically modeling similar flows (multiphase flows, lot of empiricism). I have got results which are analogous to yours, I was working on MicroProcessor Cooling Techniques, Which is reversal of your application.


January 24, 2013, 21:06 
all informations of my simulation1

#6 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
The geometry is very simple. The length of tube is 500 mm. The diameter is 3.2mm. There is no wall thickness, and we only consider the fluid region. The uniform heat flux is 800kW/m2, the mass velocity is 1000kg/m2s , the temperature of inlet is 670 K and the pressure at outlet is 25 MPa. Flow direction is upward.


January 24, 2013, 21:07 
all informations of my simulation2

#7 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
1. The needed thermal properties of water are as follow:
density Number of interval 3 Interval ranges [663,700, 800, 1073] Numbers of coefficients [5,5,4] Coefficient [20888437.32897, 121466.279457929, 264.892650225018, 0.2567597311521, 0.000093333082540744, 79195.9972248006, 404.619712014594, 0.7784131025619, 0.00066713291581415, 2.1477224361164E07, 902.76678358608, 2.2681931455009, 0.002078145302, 6.5474813334491E07] Exponents [0.0,1.0,2.0,3.0,4.0, 0.0,1.0,2.0,3.0,4.0, 0.0,1.0,2.0,3.0]  Specific heat Number of interval 3 Interval ranges [663,683,783, 1073] Numbers of coefficients [5,5,5] Coefficient [72098255169.205, 426901983.90771, 947913.913893974, 935.475678541713, 0.3462040914552, 35740652.493319, 189447.825833259, 376.86407580743, 0.3333750897413, 0.00011063602503477, 223117.120300903, 870.188643358111, 1.2944485418665, 0.00085908519886809, 2.1442187362133E07] Exponents [0.0, 1.0, 2.0, 3.0, 4.0, 0.0, 1.0, 2.0, 3.0, 4.0, 0.0, 1.0, 2.0, 3.0, 4.0]  Dynamic viscosity ($Temperature<663.0)? 3.15716e5: (($Temperature<683.0)? 6.92242646247360.0408832952645*$Temperature +0.000090549079244015* pow($Temperature,2)8.913669797436E08* pow($Temperature,3)+3.2906108603975E11* pow($Temperature,4): (($Temperature<783.0)? 0.01089895126960.000057839685831513*$Temperature +1.1532390544032E07* pow($Temperature,2)1.0216997356426E10* pow($Temperature,3)+3.3955858119227E14* pow($Temperature,4): (($Temperature<1073.0)? 0.0000032427270621927+4.741942361337E08*$Temperature +4.78555217991E12* pow($Temperature,2): 4.21E5)))  Thermal conductivity ($Temperature<663.0)? 0.2264: (($Temperature<695.0)? 47661.033038555278.416789654012*$Temperature +0.6099315425136* pow($Temperature,2)0.00059388587645495* pow($Temperature,3)+2.1685650279812E07* pow($Temperature,4) : (($Temperature<775.0)? 185.8344102193190.9793944425447*$Temperature +0.0019389248748* pow($Temperature,2)0.0000017078378939546* pow($Temperature,3)+5.6465757723963E10* pow($Temperature,4): (($Temperature<1073.0)? 2.41995127661810.0094779285294*$Temperature +0.000014350317833471* pow($Temperature,2)9.5997271574745E09* pow($Temperature,3)+2.421812079001E12* pow($Temperature,4): 0.123)))  Prantl Number ($Temperature<663.0)? 3.9489: (($Temperature<700.0)? 712655.3702821414137.78650804068*$Temperature +9.0095886619281* pow($Temperature,2)0.0087191263008* pow($Temperature,3)+0.0000031643427014962* pow($Temperature,4): (($Temperature<800.0)? 1789.01131665129.1712038684109*$Temperature +0.017667884118* pow($Temperature,2) 0.00001514546110914* pow($Temperature,3)+4.8732160696751E09* pow($Temperature,4): (($Temperature<1073.0)? 35.19336111290590.132228857307*$Temperature +0.00019253019932984* pow($Temperature,2)1.2511500694886E07* pow($Temperature,3)+3.0554723673796E11* pow($Temperature,4): 0.91))) 

January 24, 2013, 21:10 
all informations of my simulation3

#8 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
2. The generated mesh is shown in Fig.1.


January 24, 2013, 21:12 
all informations of my simulation4

#9 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
3. The selected physical model


January 24, 2013, 21:16 
all informations of my simulation5

#10 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
4. the wall y+ is from 16 to 33


January 24, 2013, 21:21 
all informations of my simulation6

#11 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
There are totally three boundaries: Inlet Outlet and wall


January 24, 2013, 21:24 
all informations of my simulation7

#12 
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
A lineprobe is generated on the wall boundary to be used to get the wall temperature.
Nine plane sections are generated along Z direction to get coolant temperature by using Surface Average of the Reports 

January 24, 2013, 21:37 

#13  
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
Quote:


January 24, 2013, 22:14 

#14  
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
Quote:
I can understand the entrance effect at length < 0.1 m. But the wall temperature distribution pattern is more like a simulation result with constant thermal properities, not with changed thermal properities. In fact, in the high enthalpy region of supercritical water, all thermal properties of water are functions of coolant temperature. I defined all these thermal properties using field function of STAR CCM, or polynomial density and specific heat. I attached all informations used in my simulation. Are there any important steps that I missed during the simulation process? Thank you very much anyway! 

January 25, 2013, 01:03 

#15 
Senior Member
Pauli
Join Date: Mar 2009
Posts: 190
Rep Power: 10 
Did you use a uniform inlet velocity profile? Or did you apply a parabolic profile representing fully developed flow? From the information you provided, it appears you specified a uniform profile. I believe that is the primary source of your high h at the inlet.
Your setup panel shows you chose a 2layer all y+ wall treatment. The y+ plot shows your yplus is in the range where wall functions are used. With wall functions, the heat transfer function coefficient is strongly driven by the near wall velocity. Your mesh & y+ values are nice. I'm guessing the red spot in the y+ plot is the inlet. That observation coupled with your setup panel led me to believe you have a uniform inlet velocity profile. 

January 25, 2013, 01:31 

#16  
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
Quote:
I did use uniform inlet velocity profile in my simulation. The problem of my simulation is that the simulation results are more like we use constant thermal properties, rather than changed thermal properties. I have defined all thermal properities as the functions of the coolant temperature. But it seems all these functions are not rightly used. I don't know what causes this problem 

January 27, 2013, 10:54 

#17 
New Member
Prashanth Avireddi
Join Date: Sep 2012
Location: Farmington Hills, MI
Posts: 10
Rep Power: 7 
In the Physics model, change the fluid setting from Liquid (H2o) to Multiphase Mixture with inlet being only Liquid H2o. That should solve the problem if there was any.


January 30, 2013, 04:12 

#18  
New Member
Zhu
Join Date: Jan 2013
Posts: 18
Rep Power: 6 
Quote:
I met a problem here, as shown in attached figure. If I want to use only Liquid H2O, the vomume fraction of velocity inlet should be defined as [1.0, 0.0]. But I can not define it in STAR CCM. I was informed that I gave an Invalid array input. How to solve this problem? Would you please give me a bit more detailed information about how to set multiphase mixture model? I am looking forward to your help. Thank you very much!! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
UDF compiling problem  Wouter  Fluent UDF and Scheme Programming  6  June 6, 2012 04:43 
Gambit  meshing over airfoil wrapping (?) problem  JFDC  FLUENT  1  July 11, 2011 05:59 
natural convection problem for a CHT problem  SeHee  CFX  2  June 10, 2007 06:29 
Adiabatic and Rotating wall (Convection problem)  ParodDav  CFX  5  April 29, 2007 19:13 
Is this problem well posed?  Thomas P. Abraham  Main CFD Forum  5  September 8, 1999 14:52 