CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Phoenics

Velocity in Porous medium : HELP! HELP! HELP!

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 6, 2006, 06:10
Default Velocity in Porous medium : HELP! HELP! HELP!
  #1
Kali Sanjay
Guest
 
Posts: n/a
Dear Friends!

I have a query regarding solving of velocity in a porous domain. I have a sort of tube which is filled with soil in a section which causes the area to decrease. The tube is left to open at both ends. I am solving for a velocity profile and the source of velocity is due to electroosmosis in the soil which is a function of x(sort of pumps in series of different capacities ).

--Non-porous--|---Soil-Porous---|--Non-porous--

So soil acts like a pump. I am solving for 1-D Cartesian coordinates and I can write the velocity as

d((Ke*dV/dx)+(Kh*dp/dx)/dx=0

Actually the part Ke*dV/dx is source of velocity which I can calculate. The second part of the above equation will be calculated for satisfying continuity equation. As the velocity is different , the pressure actually will compensate for making it uniform. I hope I can solve this problem using PHOENICS.

I gave the porosity as

CONPOR(CAREA,PRST,CELL,#2,#2,1,1,1,1)

Where my region 2 is packed with soil.

I calculated Ke*dv/dx at every point in the soil and specified as velocity source in ground.

PATCH(VSPD,CELL,1,NX,1,1,1,1,1,LSTEP) COVAL(VSPD,U1,FIXVAL,GRND)

In GROUP 13 for GRND (Value), I specified Ke* dv/dx........v is potential (volts)

I specified pressure boundary conditions as P=0 at ends with initial Pressure P specified as Zero at all points.

These are my inital and boundary conditions for U1 and P1

FIINIT(U1)=1E-10 FIINIT(P1)=0

PATCH(ANODE,WEST,1,1,1,1,1,1,1,LSTEP) COVAL(ANODE,P1,FIXP,0.0) COVAL(ANODE,U1,ONLYMS,0.0)

PATCH(CATHODE,EAST,NX,NX,1,1,1,1,1,LSTEP) COVAL(CATHODE,P1,FIXP,0.0) COVAL(CATHODE,U1,ONLYMS,0.0)

Is it the correct way? I have found that when CONPOR is specified as CELL, the AEAST shows the same area as non-porous region and I think the volume decreases, but area remains same (Internally propably PHOENICS takes volume as DXU2D *AEAST*Porosity). In that case the velocity that I need to specify should be superficial velocity. What I know Ke * dv/dx is actual velocity in pore. Do I need to multiply this velocity with porosity and give as area patch for velocity?

I am not sure whether I had pose my problem correctly or not, in summary, 1. Is CONPOR fixing the way I did is OK? 2. Can I specify known velocity field in the cells and set PHOENICS to solve for pressure to satisfy the continuity equation and get uniform velocity? I just know pressure at ends (atmosphere) and velocity source Ke*dv/dx in soil. 3. Do I need to activate DARCY=T and fix DARCON for #2 region?

Please please help me. I am really struck up.

Thank you in advance,

Regards,

Sanjay

  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
porous medium J Siemens 2 January 24, 2007 06:09
porous medium rati FLUENT 1 June 29, 2006 18:18
Velocity vs Pressure in Porous Medium K Sanjay Phoenics 0 February 23, 2006 06:39
porous medium antonietta Main CFD Forum 0 March 15, 2005 04:37
porous medium-velocity inlet a qn? Shankar FLUENT 2 July 25, 2002 13:30


All times are GMT -4. The time now is 09:53.