On Thu, Jul 10, 2008 at 1:02 PM, Tony S Yu <[EMAIL PROTECTED]> wrote: > > > Thanks, this solution works as you suggest. However, I don't understand the > math.
There is not really much to it. It just fixed one cell to have a permanent value.. > It looks like the last term is what I would normally consider to be a > source/sink, with the magnitude of the term being proportional to the flux > into/out-of the system. This is similar. Essentially, the cell is being removed from the calculation and acting as a fixed value boundary condition for its neighboring cells. The boundary condition is at the cell center rather than the faces. > But, I don't understand how ImplicitSourceTerm works > (I looked at the docstring, but I'm still confused). ImplicitSourceTerm(coeff) adds "coeff * cellVolume" to the diagonal of the matrix. If you have a source "a * phi" for example, where phi is the variable that you are solving for, then you would write ImplicitSourceTerm(a). The phi is included implicitly. > I don't want to waste your time, but if you could point me to a good > reference on this type of source term, I would be very grateful. Take a look at Patanker "Numerical Heat Transfer and Fluid Flow" for a simple reference. Cheers -- Daniel Wheeler
