diff --git a/examples/telemac3d/stratif_wind/doc/stratif_wind.tex b/examples/telemac3d/stratif_wind/doc/stratif_wind.tex index 692cb7fee0a8d52519c61c804c3f1122613e2ebe..c9cc0bf00de75833f961cee2d0804e734c9e10a8 100644 --- a/examples/telemac3d/stratif_wind/doc/stratif_wind.tex +++ b/examples/telemac3d/stratif_wind/doc/stratif_wind.tex @@ -132,7 +132,7 @@ and \telkey{ACCURACY FOR PPE} was set to 10$^{-4}$ (default value until release \subsubsection{Improved numerical parameters since 2015} To greatly improve the results (less diffusion for temperature field), -the \telkey{ACCURACY FOR PPE} has to be decreased (from 10$^{-4}$ to 10$^{-7}$ +the \telkey{ACCURACY FOR PPE} has to be decreased (from 10$^{-4}$ to 10$^{-8}$ e.g.). It also enables to decrease the differences between sequential and parallel runs, in particular for temperature. @@ -148,6 +148,9 @@ temperature field: The extra cost for CPU time is very moderate for this test-case compared to the PSI-type MURD scheme. +Using \telkey{TIDAL FLATS} = YES (default choice) improves mass conservation +(water depth up to machine accuracy, but only a little bit for temperature). + Setting the 3 keywords \telkey{IMPLICITATION FOR DEPTH} and \telkey{IMPLICITATION FOR VELOCITIES}, \telkey{MASS-LUMPING FOR DEPTH} to their default values (0.55 ; 0.55 ; 0.) does not significantly change the results and diff --git a/examples/telemac3d/stratif_wind/f3d_stratif_wind.slf b/examples/telemac3d/stratif_wind/f3d_stratif_wind.slf index 72e944ccc0c7241fa1a62a949d75f2df063324b0..329071f57411b43de4fcf8b6b5c723a4d1371123 100644 --- a/examples/telemac3d/stratif_wind/f3d_stratif_wind.slf +++ b/examples/telemac3d/stratif_wind/f3d_stratif_wind.slf @@ -1,3 +1,3 @@ version https://git-lfs.github.com/spec/v1 -oid sha256:1012400168ee988d2694d4c60eb2fdbcc4c29880e0ec2482e67c8366c11f2dc5 +oid sha256:52db0bbf74c7d7ac577d703af061958b9818da3ec6c8a05eecaae69d83840b9f size 7340940 diff --git a/examples/telemac3d/stratif_wind/t3d_stratif_wind.cas b/examples/telemac3d/stratif_wind/t3d_stratif_wind.cas index 877ab7246bb2f1f7897b8020fd0e8c4e2461c45f..0c474b2a82d9d6b279540a316eb03f22e5aafe7b 100644 --- a/examples/telemac3d/stratif_wind/t3d_stratif_wind.cas +++ b/examples/telemac3d/stratif_wind/t3d_stratif_wind.cas @@ -24,8 +24,6 @@ NUMBER OF HORIZONTAL LEVELS : 40 /---------------------------------------------------------------------- / NUMERICAL OPTIONS /---------------------------------------------------------------------- -TIDAL FLATS : NO - MASS-BALANCE : YES /CONSTANT VISCOSITY VERTICAL TURBULENCE MODEL : 1 @@ -35,9 +33,12 @@ COEFFICIENT FOR VERTICAL DIFFUSION OF VELOCITIES : 0.01 SCHEME FOR ADVECTION OF VELOCITIES : 1 ACCURACY FOR DIFFUSION OF VELOCITIES : 1.E-8 -ACCURACY FOR PPE = 1.E-7 PRECONDITIONING FOR DIFFUSION OF VELOCITIES = 34 +/ +SOLVER FOR PPE = 1 PRECONDITIONING FOR PPE = 34 +/ +MASS-LUMPING FOR DEPTH = 1. /---------------------------------------------------------------------- / PHYSICAL OPTIONS /---------------------------------------------------------------------- @@ -65,10 +66,5 @@ WIND VELOCITY ALONG X : 10. WIND VELOCITY ALONG Y : 0. /---------------------------------------------------------------------- / -/ DEFAULT VALUES UNTIL V8P0 KEPT FOR NON REGRESSION -SOLVER FOR PROPAGATION = 1 -ACCURACY FOR PROPAGATION = 1.E-6 -SOLVER FOR PPE = 1 -/ / DEFAULT VALUE UNTIL V8P1 KEPT FOR NON REGRESSION COEFFICIENT OF WIND INFLUENCE VARYING WITH WIND SPEED = NO diff --git a/examples/telemac3d/stratif_wind/vnv_stratif_wind.py b/examples/telemac3d/stratif_wind/vnv_stratif_wind.py index 54dd0abc61050b91353d04dd2f1bb59a8a58e68e..b2ea508fbfd47a0d9685f261c43a4d675e389473 100644 --- a/examples/telemac3d/stratif_wind/vnv_stratif_wind.py +++ b/examples/telemac3d/stratif_wind/vnv_stratif_wind.py @@ -50,17 +50,17 @@ class VnvStudy(AbstractVnvStudy): # Comparison with the last time frame of the reference file. self.check_epsilons('vnv_1:T3DRES', 'f3d_stratif_wind.slf', - eps=[1.E-5, 1.E-5, 1.E-5, 1.E-5, 0.04]) + eps=[1.E-6, 1.E-5, 1.E-5, 1.E-5, 2.E-3]) # Comparison with the last time frame of the reference file. self.check_epsilons('vnv_2:T3DRES', 'f3d_stratif_wind.slf', - eps=[1.E-5, 2.E-5, 1.E-5, 1.E-5, 0.05]) + eps=[1.E-6, 1.E-5, 1.E-5, 1.E-5, 6.E-3]) # Comparison between sequential and parallel run. self.check_epsilons('vnv_1:T3DRES', 'vnv_2:T3DRES', - eps=[1.E-5, 2.E-5, 1.E-5, 1.E-5, 0.05]) + eps=[1.E-6, 1.E-5, 1.E-5, 1.E-5, 6.E-3]) def _post(self):