Same behaviour with fluiddyn/fluidsim # To see the bugs ----------------- python no_energy_conservation_ns2dstrat.py """ # import numpy as np import matplotlib.pyplot as plt # from fluidsim.solvers.ns2d.strat.solver import Simul from fluidsim.solvers.ns2d.solver import Simul params = Simul.create_default_params() # params.oper.nx = nx = 128*2 # params.oper.ny = ny = nx // 4 params.oper.nx = nx = 64 params.oper.ny = ny = nx params.oper.Lx = 10. params.oper.Ly = params.oper.Lx * (ny / nx) # params.oper.coef_dealiasing = 0.5 # it is completely wrong! params.nu_8 = 1e-8 # there is clearly a problem # params.nu_8 = 1e-10
def _complete_params_with_default(params):
"""This static method is used to complete the *params* container.
"""
SimulNS2D._complete_params_with_default(params)
attribs = {"N": 1.0, "NO_SHEAR_MODES": False}
params._set_attribs(attribs)