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