Hi Riadh,
unfortunately decreasing the time step to 1/10 of the employed time step when using another scheme doesn't help. I'm using v7p1r1.
The problem arises only when using MOC for the advection of k-e model. For example the combination MOC for velocities and MURD for k-e works well.
Here my steering file with a time step of 0.005s.
N.B.: all the other schemes work well with a time step of 0.05s!
________________________________________________________________
PARALLEL PROCESSORS =4
/
/ EQUATIONS
/
TURBULENCE REGIME FOR THE BOTTOM =2
TURBULENCE REGIME FOR LATERAL SOLID BOUNDARIES =2
COEFFICIENT FOR HORIZONTAL DIFFUSION OF VELOCITIES =1.E-6
COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES =1.E-6
HORIZONTAL TURBULENCE MODEL =3
VERTICAL TURBULENCE MODEL =3
LAW OF BOTTOM FRICTION =5
LAW OF FRICTION ON LATERAL BOUNDARIES =5
FRICTION COEFFICIENT FOR THE BOTTOM =0.0005
FRICTION COEFFICIENT FOR LATERAL SOLID BOUNDARIES =0.0005
NON-HYDROSTATIC VERSION =YES
OPTION FOR THE BOUNDARY CONDITIONS OF K-EPSILON =1
/
/ EQUATIONS, ADVECTION
/
SCHEME FOR ADVECTION OF VELOCITIES =1
SCHEME FOR ADVECTION OF K-EPSILON =1
/
/ EQUATIONS, BOUNDARY CONDITIONS
/
VELOCITY VERTICAL PROFILES =2;2 /2:logarithmic
VELOCITY PROFILES =1;1
/
/ EQUATIONS, INITIAL CONDITIONS
/
INITIAL ELEVATION =0.189
INITIAL CONDITIONS ='CONSTANT ELEVATION'
PRESCRIBED FLOWRATES =0.0;0.0453
PRESCRIBED ELEVATIONS =0.189;0.0
/
/ INPUT-OUTPUT, FILES
/
/ LIQUID BOUNDARIES FILE ='Hydrograph_to_steady.txt'
GEOMETRY FILE ='mesh.slf'
BOUNDARY CONDITIONS FILE ='BOTTOM_BCmesh.cli'
3D RESULT FILE ='r3d_MOC_ts0.005.slf'
2D RESULT FILE ='r2d_MOC_ts0.005.slf'
/
/ INPUT-OUTPUT, GRAPHICS AND LISTING
/
LISTING PRINTOUT PERIOD =120
VARIABLES FOR 3D GRAPHIC PRINTOUTS ='Z,U,V,W,K,EPS'
GRAPHIC PRINTOUT PERIOD =3000
VARIABLES FOR 2D GRAPHIC PRINTOUTS ='U,V,H,S,B,US'
MASS-BALANCE =true
/
/ NUMERIC PARAMETERS
/
NUMBER OF HORIZONTAL LEVELS =10
/
/ NUMERICAL PARAMETERS
/
FREE SURFACE GRADIENT COMPATIBILITY =0.9
TIME STEP =0.005
NUMBER OF TIME STEPS =120000
/ IMPLICITATION FOR DEPTH =1.0
/ IMPLICITATION FOR VELOCITIES =1.0
________________________________________________________________
Best regards,
Clemens
