Dirchilet boundary condition implementation in FVM
Hi,
I am solving 1D non-linear hyperbolic equation - the Buckley Leverett equation using Finite volume methods. Background information: Here is the equation:Buckley Leverett http://www.cfd-online.com/Forums/dat...AASUVORK5CYII= I am solving this system using upwinded (for simplicity purposes right now) using fully implicit solver (I have to use fully implicit method due to other reason). I solve it implicitly at each time step using Newton's method - a standard way in reservoir simulation. Question: I am having some troubles in implementing Dirichlet boundary condition i.e. how to impose a boundary condition on the inlet face in the numerical framework. 1) Should I just add another cell on the left (which is the boundary) of same length "dx" as other cells or should I add a cell of very large length? 2) Further, do I consider this cell - ghost cell as part of my system of equations and solve it? Thanks in advance for your help. |
All times are GMT -4. The time now is 07:58. |