Flow through porous media

marango · June 28, 2010, 12:40

Hello everybody,

I'm a newbie in OpenFOAM and have some problems with the implementation in OpenFOAM. As you can see in the title I'm examing flow through porous media using a new approach.

Similar to the Darcy-Forchheimer-Equation I want to add a new term to rhs of the incompressible Navier-Stokes-Equation. As a basis I use the simple-solver in OpenFOAM. Now I'd like to change the implemented differential equation. The rhs looks in x-direction as follows:

... = dp/dx - a* u_x - b* (u_x)² .

( a,b - scalar values)
I tried a lot to implement it. But certainly I don't know the right commands. Can anybody help me?

Many thanks in advance!

herbert · June 28, 2010, 12:46

Hi marango,

you don't have to do anything on your own. Please have a look on rhoPorousSimpleFoam and its tutorials. This is a compressible solver, but you can use nearly the same commands for incompressible solvers.

Regards,
Stefan

marango · June 28, 2010, 14:15

Thanks, Stefan, for your instant answer.

I already used the rhoPorousSimpleFoam-solver and implemented an incompressible version of this solver. Unfortunately the results were completely wrong.

So I use a new approach, which also changes the lhs of the NSE. This is working right, but only the rhs makes problems.

I also considered to use some commands from rhoPorousSimpleFoam, but either I couldn't find out what this command actually does or the mathematical background doesn't suite to my approach.

June 28, 2010, 12:40	Flow through porous media	#1
marango New Member Join Date: Jun 2010 Location: Germany Posts: 13 Rep Power: 16	Hello everybody, I'm a newbie in OpenFOAM and have some problems with the implementation in OpenFOAM. As you can see in the title I'm examing flow through porous media using a new approach. Similar to the Darcy-Forchheimer-Equation I want to add a new term to rhs of the incompressible Navier-Stokes-Equation. As a basis I use the simple-solver in OpenFOAM. Now I'd like to change the implemented differential equation. The rhs looks in x-direction as follows: ... = dp/dx - a* u_x - b* (u_x)² . ( a,b - scalar values) I tried a lot to implement it. But certainly I don't know the right commands. Can anybody help me? Many thanks in advance!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
no converge - for Low Re Porous media seepage flow	lee	CFX	3	August 31, 2007 02:12
Modeling porous media flow in prostar	Nirmal Kumar Pathak	Siemens	8	November 23, 2005 08:18
Flow Through Porous Media + Biochemical Reactions	Paresh Jain	CFX	0	September 20, 2003 06:06
How to simulate flow in porous media?	may song	Phoenics	2	August 29, 2002 00:10
Flow resistance of porous media	Desheng Kong	Main CFD Forum	3	December 19, 2000 12:07

June 28, 2010, 12:46		#2
herbert Senior Member Stefan Herbert Join Date: Dec 2009 Location: Darmstadt, Germany Posts: 129 Rep Power: 18	Hi marango, you don't have to do anything on your own. Please have a look on rhoPorousSimpleFoam and its tutorials. This is a compressible solver, but you can use nearly the same commands for incompressible solvers. Regards, Stefan

June 28, 2010, 14:15		#3
marango New Member Join Date: Jun 2010 Location: Germany Posts: 13 Rep Power: 16	Thanks, Stefan, for your instant answer. I already used the rhoPorousSimpleFoam-solver and implemented an incompressible version of this solver. Unfortunately the results were completely wrong. So I use a new approach, which also changes the lhs of the NSE. This is working right, but only the rhs makes problems. I also considered to use some commands from rhoPorousSimpleFoam, but either I couldn't find out what this command actually does or the mathematical background doesn't suite to my approach.