CFL condition is based on water flux
The CFL condition currently estimates the maximum time step based on the water flux j_w
Summary
The water flux \mathbf{j}_w = \theta_w \mathbf{v}
is not the velocity of the solute. Instead, the water flux has to be divided by the water content to yield the water velocity \mathbf{v} = \mathbf{j}_w / \theta_w
.
What is the current bug behaviour?
CFL condition is based on water flux.
What is the expected correct behaviour?
CFL condition is based on convective velocity.
Ideas how to fix this?
Pass the water content grid function to the CFL condition function and divide the local water flux by the local water content to receive the local velocity.
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information