
[Sponsors] 
June 17, 2014, 12:15 
Comsol CFD with porous media

#1 
New Member
Carreira
Join Date: Jun 2014
Posts: 4
Rep Power: 3 
Hi everyone,
I would like to ask your for help in COMSOL. I'm trying to join three consecutive tubes, the first with turbulent flow, the second with a porous medium, and the last turbulent again. I'd like to solve all the parameters continuously, in particular the speed and pressure. I discovered that it is enough to put the same name of parameters in both physicals ( u,v,w,p), but the speed is still not continuous. That is, at the end of the first tube velocity falls to zero, the second starts at zero and so forth. Does anybody know which is the better way to solve this problem? Any suggestion? Thanks, Carreira 

June 17, 2014, 14:19 
Interesting

#2 
New Member
Иван Сташко
Join Date: Dec 2013
Posts: 27
Rep Power: 4 
Carreira,
I haven't seen this happen. Do you have a sample .mph file that you could post? I would look at it for you. 

June 18, 2014, 07:34 
mph file

#3 
New Member
Carreira
Join Date: Jun 2014
Posts: 4
Rep Power: 3 
Hi Ivan, thank you for your interest!!
Here is the sample file: https://www.dropbox.com/s/kax6kztmdc...ous_3tubes.mph you can clearly see (in the second chart 1D plot) that the velocity drops to zero. I may have to make a merge of the mesh, but I don't know how to do this. Thank you! 

June 18, 2014, 16:19 
Inlet and outlet defined

#4 
New Member
Иван Сташко
Join Date: Dec 2013
Posts: 27
Rep Power: 4 
Carreira,
At least in the model file which I received, there is no inlet and outlet boundary condition. The only boundary conditions are the default ones, Fluid Properties, Wall, and Initial Values. You need to define an Inlet and an Outlet. Without these, there will be no flow. Do you have these defined in your model on your system? Ivan 

June 18, 2014, 16:45 
wrong file

#5 
New Member
Carreira
Join Date: Jun 2014
Posts: 4
Rep Power: 3 
Yes i have,
but... I think you didn't open the correct file because I did't put any materials and in this image I can see two materials and....no porous media. https://www.dropbox.com/s/kax6kztmdc...ous_3tubes.mph Ty Ivan 

June 18, 2014, 19:30 
Interesting

#6 
New Member
Иван Сташко
Join Date: Dec 2013
Posts: 27
Rep Power: 4 
OK, now I see the behaviour.
It looks to me like it is set up correctly. Normally, however, I do assign a material to every region. But I see you have the viscosity and density specified directly. I don't know if it makes a difference, maybe you should assign a material to all regions just to check. Also it appears the continuity equation is violated. If V is high going into the porous section and low coming out of it, then there must be mass buildup in the porous section. One other interesting thing, you have the pressure as p in the porous section and 0 in the turbulent section. Normally, you can open the Model Inputs section of the Fluid Properties tab, and there will be a reference pressure. But this model does not offer this ability. I have seen this before, as well, and I don't know how to consistently get to and set the reference pressure. 

June 19, 2014, 06:01 
continuity equation

#7 
New Member
Carreira
Join Date: Jun 2014
Posts: 4
Rep Power: 3 
Do you have any ideia how to solve the problem with continuity equation?
I added one material to my geometry and this made no difference...in fact i can't converge the solution. Last edited by Carreira; June 19, 2014 at 08:32. 

June 19, 2014, 16:21 

#8  
New Member
Иван Сташко
Join Date: Dec 2013
Posts: 27
Rep Power: 4 
Quote:
The continuity equation, which is commonly "\div J = 0" appears to be part of that equation set, written as "\rho \div u = 0". It just *looks* like it's violated, you would have to integrate across one of the crosssectional planes to actually check. They should all have the same amount af flux. (Expressed in your system, it would be the integral at one crosssectional plane minus the integral at another crosssectional plane should be zero.) And that brings up another point. Make sure you integrate "\rho u" across the crosssections. Since "\rho" is involved, this means looking at a map of u may be misleading. On the other hand, the density shouldn't change that much. By the way, I'm still curious how to get to the reference pressure, which is under Fluid Properties > Model Inputs. As I have said, I've seen it in some models, but not it all, even when they employ the exact same physics, and I don't know how to make it just appear. Let me know if you find out. Your pressure results look correct, but I do not see them for the outlet. Could you generate them? Why are you using the komega model instead of kepsilon? Have you tried to run it using kepsilon? 

June 20, 2014, 11:45 
ke v komega

#9 
New Member
Иван Сташко
Join Date: Dec 2013
Posts: 27
Rep Power: 4 
Carreira,
I have been using ke with great success in models without curl, such as long pipe sections similar to your model. But you have selected komega. I ran across this thread: Use of kepsilon and komega Models and I recommend you read it, particularly near the bottom, vivek**'s comment about the applicability of ke. It seems, after reading this thread, I would chose ke as well for your model. I'm still curious how this turns out. 

Tags 
comsol, multiphysics, porous media, turbulent kw 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Porous media setup issues in Fluent  Bernard Van  FLUENT  26  July 24, 2015 18:24 
about porous media  major  FLUENT  5  March 6, 2013 11:00 
Porous Media coupled with internal flow  Samuel Andrade  FLUENT  2  August 26, 2012 09:43 
species mass source in porous media ?  PK  FLUENT  0  February 16, 2007 12:12 
porous media: Fluent or StarCD?  Igor  Main CFD Forum  0  December 5, 2002 16:16 