# porous inertial resistance

 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 31, 2009, 06:57 porous inertial resistance #1 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Hallo, I want to model an tube bundle heat exchanger. The tube bundles are replaced by a porous media. My problem is, that the porous inertial resistance is a function of to the fluid velocity. Is it possible to create a table, where the porous inertial resistance is shown for x,y,z -components with the according velocity? Star-CCM+ should compute with the porous inertial resistance for the equal velocity from the table. Lukas

 August 31, 2009, 09:03 #2 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 Your inertial resistance is in function of your velocity??? OK I have never heard of this but for everything a first time. When i derive my inertial and viscous factor, i usually have experimental pressure drops vs volume flow rates. Convert these in a plot where you have velocity on the X axis and pressure loss at the Y axis, then use excel or matlab to derive this 2nd degree polyline. Convert according to porous media expression and voila you have your inertial and viscous resistance. dp/dl=Av2+Bv use A and B to find inertial and viscous If now these factors are again dependant of velocity, this would mean a third degree polyline, just derive 3rd degree factors and make sure to use the correct ones to form your inertial and viscous terms.

 August 31, 2009, 10:44 #3 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 I calculated the pressure drop with equations from the VDI heat atlas. The porous inertial resistance can be calculated like this: Pi = dp/dl * 1/vē Is it allowed to remove the porous viscous resistance because of high Re-numbers (Re about 10^5)? dp is a function of velocity, so Pi is a function of velocity, too?!

 August 31, 2009, 15:45 #4 Senior Member   Aroon Join Date: Apr 2009 Location: Racine WI Posts: 148 Rep Power: 15 You can model your inertial (as well as viscous) porous media co-efficient as a function of anything you want (Temp, velocity). Define new Field Function using Tools > Field Function > User Defined Field Function. Name it appropriately and then in Porous Inertial Resistance select Field Function under Method and you can set it to the defined field function. Hope this helps

 September 1, 2009, 04:30 #5 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Thanks for your answer, but how I have to declare the field function, that it is according to my table? Lukas

 September 1, 2009, 09:07 #6 Senior Member   Aroon Join Date: Apr 2009 Location: Racine WI Posts: 148 Rep Power: 15 Hello Lukas, I guess you have a table or a function of Pi = f(V). You can enter this function in the new field function. If you have a table you can always fit a curve and define the new field function based on the defined curve. Regards Aroon

 September 1, 2009, 09:17 #7 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Ok, it sounds plausible, but Iīm not a expert in using Star-CCM+ Can you give me a step by step manual to do that. I just managed to create a table where the porous inertial resistance is shown for x,y,z -components with the according velocity. But I canīt combine this table with my field function. Lukas

 September 2, 2009, 09:21 #8 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 Visharoon, I dont think he needs resistance in function of velocity, Lukas Pi = dp/dl * 1/vē --> this is your formula, so I have no idea where you get that Pi is in function of v??? Pi is a constant value, the second derative from your 2nd order poly that is dp=f(v) Last edited by bramv101; September 2, 2009 at 10:30.

 September 2, 2009, 12:15 #9 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 The problem is ,that my velocity isnīt constant, so my Pi changes, too. Lukas

 September 8, 2009, 08:17 #10 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 Thats the whole idea of a porous media coefficient, your Pi and Pv stay constant while your dp and v change. for dp/dl=Pi*v^2+Pv*v --> dl, Pi and Pv remain constant Just plot your pressure drop (Y axis) and your velocity (X axis) in excel, you will see that this is a perfect 2nd degree polyline. If you have a Pi dependant on your velocity, that is third degree and this is highly unlikely. Cant imagine such a media.... Inertial = mass*acceleration

 September 8, 2009, 10:09 #11 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 This is fact. Iīm aware of this. But the velocity of my fluid, before it permeates the porous media, is between 0 and 2,5 m/s. And for this velocity I have to define my coefficients. If i have a constant velocity before the fluid permeates the porous media, I will have constant coefficients. Because of the various velocities I want to make the coefficients a function of the velocity.

 September 8, 2009, 10:35 #12 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 I understand, but when you say that your velocity goes from 0 to 2.5 [m/s]. Do you mean distributed over the face of the porous media? or is varying in time? Either way, I still can not find a good enough reason to model your resistance values in function of velocity. In terms of physics this doesn't make sense... Say you can model your dp=f(v) with a third degree poly, what is your goal?

 September 10, 2009, 06:17 #13 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Yes, the velocity is distributed over the face of the porous media. Because of this I think, that I will make a failure, if I calculate one coefficient for this different velocities.

 September 10, 2009, 06:46 #14 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 Unless your velocities range from 0 to 5000 m/s (unlikely), go for 1 constant value. If your resistance is constant and your velocity fluctuates then your pressure is a direct result of this. Adapting your resistance to your velocity is like adapting the young modulus of a material to the local force applied -> very wrong So in porous materials inertial and viscous resistance are material properties, these are inherent to the material itself and therefore a material constant. If you have experimental pressure drops vs approach velocity for this type of material just use Darcy to derive Pi and Pv

 September 11, 2009, 04:24 #15 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 I want to make the coefficients a function of the velocity outside of the porous media, because this velocity isnīt constant.

 September 11, 2009, 05:15 #16 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 Could you tell me what your case is a representation off: * filter applications * geological * solid/gas/liquid * ......

 September 11, 2009, 09:54 #17 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Ok, I will go for one constant Pi and Pv. Iīve got a tube bundle heat exchanger with liquid media. The tube bundle should be replaced by a porous media. If I plot the dp/L (y- axis) and v (x-axis) I got a equation like this: y= ax^2+ bx +c a is Pi, b is Pv. Do I have to respect the coefficient c?

 September 11, 2009, 11:23 #18 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 you can ignore this c value, usually this is neglible

 September 11, 2009, 11:24 #19 New Member   Join Date: Aug 2009 Posts: 12 Rep Power: 14 Ok, I have the solution. Bramv101, thanks for your help.

 September 11, 2009, 11:24 #20 Member   TonyD Join Date: Apr 2009 Posts: 85 Rep Power: 15 And to simulate u tube heat exchanger I would go for the heat exchanger tutorial from StarCCM+