- **Phoenics**
(*https://www.cfd-online.com/Forums/phoenics/*)

- - **How to calculate Velocity?
**
(*https://www.cfd-online.com/Forums/phoenics/52054-how-calculate-velocity.html*)

How to calculate Velocity?
Dear Friends! Kindly help me in implementing my probem of velocity calculation in PHOENICS.
I am solving a 1-Dimensional problem consisting of soil in a cylindrical tube fitted with porous electrodes at the end. Due to electroosmotic flow, there is flow of fluid in the pores from anode to cathode. How ever this flow is a function of potential which is non uniform along the length of the tube due to change in resistance of soil with time. But the net velocity since area is same will be uniform along the length of the soil, and the pressure will change to make the velocity uniform. So I have an equation dc/dt + v dc/dx = D d2c/dx2 This velocity v is due to electroosmosis and hydraulic conductivity in pores, then v=Ke*dphi/dx + Kh*dP/dx where v=velocity Ke=constant =Electroosmotic Permiability Kh=constant =Hydraulic conductivity dphi/dx = potential gradient dP/dx = Pressure Gradient From Continuity equation dv/dx=0 So I have an equation d(Ke*dphi/dx)/dx + d(Kh*dP/dx)/dx = 0 If I can solve this equation and find P, I can recalculate Velocity. My boundary conditions for P are atmospheric pressure at ends and I do not know the velocity at ends. Kindly suggest me how to calculate velocity by this pressure equation in PHOENICS Thanks in advance. Kali Sanjay |

Re: How to calculate Velocity?
In the Equation d(Ke*dphi/dx)/dx + d(Kh*dP/dx)/dx = 0 Ke, Kh and Phi as funtion of X are known, I need to find Pressure from this equation and activate conection term of solver for C. However, do I need to solve U1 and P1 in PHOENICS. If so how to implement my boundary conditions. The problem here at both ends the boundary condition for P is atmospheric pressure where as my d(Phi)/dx fuction is parabolic function.
Any help is greatly appreciated. I am struckup in implementation part in PHOENICS. |

All times are GMT -4. The time now is 10:47. |