|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Join Date: Dec 2009
Posts: 49
Rep Power: 15 ![]() |
I have written finite difference SIMPLE algorithm to solve the Navier Stokes equation in Matlab. My problem is that I don't know how to specify boundary condition for Poisseuille flow. If i specify inlet velocity, the solution tend to blowup. I suspect that somehow I don't specify the pressure boundary correctly.(In the matlab code if I specify the pressure boundary for inlet and outlet but no velocity boundary at the inlet, the solution does not blow up but the problem is that I cannot control the inlet velocity) I managed to solve Couette flow and Driven cavity problem with no problem though..Care to advice me on what I have to do. I attached also the code for the SIMPLE algorithm that I write..
Last edited by katakgoreng; December 27, 2009 at 08:48. Reason: need to give mo info |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
jk
Join Date: Jun 2009
Posts: 64
Rep Power: 15 ![]() |
You are solving flow in a duct basically. Use velocity inlet and pressure outlet bc. Which grid r u using ? colocated or staggered
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Join Date: Dec 2009
Posts: 49
Rep Power: 15 ![]() |
I'm using staggered grid..if i set the inlet velocity to v=1 and pressure outlet p=0, how about inlet pressure, do i need to extrapolate it with inner pressure?
|
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
jk
Join Date: Jun 2009
Posts: 64
Rep Power: 15 ![]() |
yes exactly. p(1).value = p(2).value where p is the pressure and numbers indicate the node location
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
Member
Join Date: Dec 2009
Posts: 49
Rep Power: 15 ![]() |
||
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
SIMPLE algorithm in 3D cylindrical coordinates | zouchu | Main CFD Forum | 1 | January 20, 2014 17:02 |
A simple heat transfer problem ? | iggor | FLUENT | 1 | June 15, 2009 04:42 |
SIMPLE algorithm for collocated storage | alex | Main CFD Forum | 4 | April 8, 2008 23:48 |
DPM - simple problem (I think) | Maciej | FLUENT | 0 | November 9, 2005 16:28 |
? Simple CFD problem | Davy | Main CFD Forum | 0 | August 21, 1998 06:19 |