|
[Sponsors] |
![]() |
![]() |
#1 |
Guest
Posts: n/a
|
Hi there,
Im not sure how to implement boundary conditions in a finite volume method. This is what I have so far: Discretized transport equation (pure convection) is of the form ap*phi_P -aw*phi_W -ae*phi_E = S and aP=aE+aW I'm using upwind for face values ( and U is positive) so aE=0 On a boundary, I set the coefficient of the boundary face to zero and add the flux to the source term. So on the west boundary, aW is set to zero and the source term is nonzero. However, this will make aP=(aW+aE)=0 so the TDMA solver will not converge. Any help greatly appreciated Thanks! |
|
![]() |
![]() |
![]() |
#2 |
Guest
Posts: n/a
|
At the west boundary, calculate aW (non-zero) by the usual method (from the prescribed flux, likely). Set aP=(aW+aE) = aW+0 = aW. Next, set aW to zero (without changing aP), and add the prescribed flux to the source term S (making the latter non-zero). Solve using the TDMA. Allow to cool, stir, and enjoy.
|
|
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Editing code | Joseph | Main CFD Forum | 0 | May 28, 2007 17:18 |
What is the Better Way to Do CFD? | John C. Chien | Main CFD Forum | 54 | April 23, 2001 09:10 |
Pseudo compressibility implementation free code | vijay pargaonkar | Main CFD Forum | 0 | March 18, 2001 02:49 |
own Code vs. commercial code | Bernhard Mueck | Main CFD Forum | 10 | February 16, 2000 11:07 |
Use of an hydrodynamic model integrated in the vertical by the PHOENICS code | Hafsia Zouhaier | Main CFD Forum | 2 | March 15, 1999 09:23 |