This testcase is based on: https://svn.oss.deltares.nl/repos/delft3d/branches/research/Deltares/20151209_34238_slib3D/examples/19_slib3d_McKay The information below belongs to the original McKay testcase: ============================================================ Simulatie met Gensig-model met aanpassingen op verzoek van Sitoni/Talmon (valsnelheid zand wordt beinvloed door slibconcentratie (fallsand).: Recent overleg tussen Rob En Arno hebben nog wat aanvullende aanpassingen in de subroutine BNGHAM en FALLSAND opgeleverd. RVW: Q_H model; waterstandsrand hard: alfa = 0.0 bin_cvisc = 0.252 Initiele conditie met schuine waterstand en put. 2 fracties: zandconcentratie (963, was 1000 g/l), let op dit moet stof 1 zijn!!! en uniforme slibconcentratie (fines) ( 193, was 180 g/l), let op dit moet stof lsed zijn (stof 2) Bij de uitstroomrand ontstaat er een zandhoop, zonder toevoeging van een put bij de uitstroomrand. Fallvelocity mud = 0.0 (in code); Csoil = 1193. Power in fallsand veranadert: nf_floc=2.15 power = 2./(3.-nf_floc) Deze waarde staat in de invoer, maar ook hard in de code om de parameter niet door te hoeven geven naar de BNGHAM-subroutine. FlcNf = 2.15 FlcK2 = 1. FlcK3 = 5. FlcK4 = .0070 FlcD0 = 10.e-06 FlcD = 0.44 FlWsmin = 0.3e-3 EroSchel = false EroSns = 0. EroScv = 3.e-9 EroScu = 400. EroSd50 = 20.e-6 cons_mud = true cons_kk = 2.0e-14 cons_ksig= 2.0e8 cons_eta = 1.0e5 bingham = true bin_cvisco = 0.252 bin_cnvisco = 2.35 bin_cyield = 30.0 bin_abingh = 1000 ero_intfc = true conc_intfc = 25.0 In de sourcecode wordt irobed op 1 gezet en de valsnelheid van de fines op 0.0. Let op consolidatie uitgeschakeld in deze .sed file door cons_kk op nul te zetten. Let op bij plotten settling velocity. Moet de grootheid op het generieke sigma rooster zijn, de andere is afkomstig uit FALLVE i.p.v. uit SETTOT. Valsnelheid moet het gedrag van dudz vertonen bij de bodem. Voor deze condities ontstaan er golven in de oplossing. Daarom simulatie herhaald met andere formulering voor de Binghmstress. Fractale dimensie aangepast van 2.15 naar 1.5 in subroutine Bngham en bin_cyield van 37.8 naar 30.0 in kk3_2frac.sed.. Simulatie gedaan op verzoek van Arno Talmon zonder turbulentuiemodel, dan verdwijnen slingeringen. Constant. Vicoww =10-3 en Dicoww=10-4.