Implementation
As already seen, we obtain by starting from the integral
conservation law the cell average update rule
This ODE is implemented in FVMIntegrator block. The use of
this block is similar to the use of the integrator block in the MOL
package.
We must specify the initial and boundary conditions and give them
to the inputs IC and gcl, gcr respectively.
The only difference here is that
instead of writing special formulas for the boundary conditions
like in the MOL package, we extend the domain with additional
cells, called ghost cells,
and fill them with values. This way, the first and the last cell of
the domain for example, can use the same formula as all other cells
in the domain.
Finally, we must pass the \emph{flux} vector F to the
F input of the integrator, which can now use the cell
average update rule with the information
just provided and compute the averages for the next time step. See
Examples implemented in the package for better understanding.
Release Notes:
.PDE.UsersGuide.FiniteVolumeMethod.Implementation