The existing implementation of the ver-1d trapping cases has numerical instabilities due to the nature of the problem and boundary condition. It is is also missing documentation detailing what equations are actually being solved.
Impact
Numerical instabilities make the model predictions not match well with the analytical solution, and also make the solution depend on the choice of initial time step size and mesh element size. A lack of documentation makes it harder for users to understand what the point of this verification case is.
The existing implementation of the ver-1d trapping cases has numerical instabilities due to the nature of the problem and boundary condition. It is is also missing documentation detailing what equations are actually being solved.
Impact
Numerical instabilities make the model predictions not match well with the analytical solution, and also make the solution depend on the choice of initial time step size and mesh element size. A lack of documentation makes it harder for users to understand what the point of this verification case is.