geo-fluid-dynamics / sapphire

An engine for constructing PDE-based simulations using Firedrake
MIT License
19 stars 5 forks source link

NSB first-order time accuracy verification fails with latest Firedrake #98

Closed agzimmerman closed 4 years ago

agzimmerman commented 4 years ago

In tests/verification/test__unsteady_navier_stokes_boussinesq.py::test__verify_first_order_temporal_convergence_via_mms, after generating

Delta_t, errors, temporal_orders 0.5, [None, 0.011107126169912135, 0.025971296457037198], None 0.25, [None, 0.006439718841114489, 0.014279010781708205], [None, 0.7864159791056798, 0.8630219973726072] 0.125, [None, 0.003462654871468789, 0.007475212830499297], [None, 0.8951191025760915, 0.9337094754111044]

the nonlinear solver diverges with Delta_t = 0.06125 after solving at t = 0.8125.