I meet a question to deal with boundary conditions. It just needs to implement inhomogeneous Dirichlet boundary conditions for saturation, not for pressure and velocity. But I do not know how to do it.
Apart from the question below, you are aware that you may only impose Dirichlet conditions on the *inflow* part of the boundary, but not on the *outflow part*. This is nontrivial to do and you'll have to think how to implement this.
The violated condition was: n_components == i->second->n_components The name and call sequence of the exception was: ExcDimensionMismatch(n_components, i->second->n_components) Additional Information: Dimension 4 not equal to 1
Like the documentation says, the Function object you provide needs to have as many components as the finite element has (even if three of the components will be ignored). In this case, that's four, but the function you provide only has one.
Best W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
