I think this is technically not a bug, but the title should explain the issue.
Expected behavior
A clear and concise description of what you expected to happen.
Minimal Reproducible Example 👇
It follows a simple stokes problem in 1D with different boundary conditions.
A) ramping up the inflow velocity linearly from t=0 to t=2 -> breaks adaptivity at t=2, but problem when adaptivity is turned off
B) ramping up the inflow velocity smoothly from t=0 to t=2 -> works adaptively and without adaptivity
Just comment in and out the first few lines in the combinations described above to reproduce the different behaviors.
...
[ Info: 1.9999991627202856
[ Info: 1.9999996552983934
[ Info: 1.9999998940463428
[ Info: 1.9999999317284627
[ Info: 1.9999999735509675
┌ Warning: dt(2.220446049250313e-16) <= eps(t)(1.9999999735509675) , and step error estimate = 15299.784213961508. Aborting. There is either an error in your model specification or the true solution is unstable (or the true solution can not be represented in the precision of Float64).
└ @ SciMLBase ~/.julia/packages/SciMLBase/JUp1I/src/integrator_interface.jl:632
[ Info: 1.9999999735509675
Environment (please complete the following information):
This is expected behavior. Setting d_discontinuities=[2.0] as a kwarg to init to tell the solver where the discontinuity occurs (or using a callback) solves it.
Describe the bug 🐞
I think this is technically not a bug, but the title should explain the issue.
Expected behavior
A clear and concise description of what you expected to happen.
Minimal Reproducible Example 👇
It follows a simple stokes problem in 1D with different boundary conditions. A) ramping up the inflow velocity linearly from t=0 to t=2 -> breaks adaptivity at t=2, but problem when adaptivity is turned off B) ramping up the inflow velocity smoothly from t=0 to t=2 -> works adaptively and without adaptivity
Just comment in and out the first few lines in the combinations described above to reproduce the different behaviors.
Error & Stacktrace ⚠️
Environment (please complete the following information):
using Pkg; Pkg.status()