Paramètres de simulation

 

    Afin de présenter les paramètres de nos simulations nous vous présentons un de type fichier .cas dans son intégralité, assorti de quelques explications.

 

 

/---------------------------------------------------------------------
/ ENTREES-SORTIES, FICHIERS
/---------------------------------------------------------------------

FICHIER DES RESULTATS              ='resultat_simulation' /fichier contenant les résultats.

FICHIER DES CONDITIONS AUX LIMITES ='conlim1' /fichier contenant les conditions limites entrées dans le mailleur Matisse.

FICHIER DE GEOMETRIE               ='geo1' /fichier contenant le maillage.

/FICHIER DU CALCUL PRECEDENT = 'res1' /si une suite de calculs est nécessaire il faut activer cette ligne.

FICHIER FORTRAN = 'ondes.f' /appelle le fichier fortran imposant la marée en frontière sud.

/---------------------------------------------------------------------
/ ENTREES-SORTIES, GRAPHIQUES ET LISTING
/---------------------------------------------------------------------

/SUITE DE CALCUL = OUI / si une suite de calculs est nécessaire il faut activer cette ligne.

PERIODE POUR LES SORTIES GRAPHIQUES =350 /une image sur 350 sera disponible en post simulation.

PERIODE POUR LES SORTIES LISTING    =20 /tous les 20 pas de temps un calcul est affiché dans le shell windows

VARIABLES POUR LES SORTIES GRAPHIQUES : 'U,V,TOB,S,B,H,Q,M,T1'

/U=vitesse selon x, V=vitesse selon y, S=côte surface libre, B=côte du fond, H=hauteur d'eau, /Q=débit, M=vitesse scalaire, T1=traceur.

VARIABLES A IMPRIMER = 'T1' /sert à afficher les valeurs prises de chaque points vis-à-vis du /traceur dans le shell.

/---------------------------------------------------------------------
/ EQUATIONS
/---------------------------------------------------------------------
DIFFUSION DES TRACEURS = OUI

NOMBRE DE TRACEURS = 1

NOMS DES TRACEURS = T1

COEFFICIENT DE DIFFUSION DES TRACEURS = 0.0005

MODELE DE TURBULENCE                  =3 /Modèle k-$\epsilon$

COEFFICIENT DE FROTTEMENT             =20.

REGIME DE TURBULENCE POUR LES PAROIS  =1 /régime turbulent lisse

COEFFICIENT DE DIFFUSION DES VITESSES =1.E-6

LOI DE FROTTEMENT SUR LE FOND         =3 /Modéle de Strickler

VENT = OUI /active l'option vent

COEFFICIENT D'INFLUENCE DU VENT = 2.16E-6 /coefficient dépendant de la vitesse

VITESSE DU VENT SUIVANT X = -5.7

VITESSE DU VENT SUIVANT Y = -15.7

/PROFONDEUR LIMITE POUR LE VENT = 2.0

/---------------------------------------------------------------------
/ EQUATIONS, CONDITIONS INITIALES
/---------------------------------------------------------------------

CONDITIONS INITIALES ='COTE CONSTANTE'

COTE INITIALE        =0.1

VALEURS INITIALES DES TRACEURS = 1. /on initialise le domaine à 1 en traceur à t=0.

VALEURS IMPOSEES DES TRACEURS = 0.;0. /on laisse évoluer ensuite librement le traceur dans le domaine.

/---------------------------------------------------------------------
/ EQUATIONS, CONDITIONS LIMITES
/---------------------------------------------------------------------

/OPTION POUR LES FRONTIERES LIQUIDES =1;1

DEBITS IMPOSES                      =0.;30.0

/PROFILS DE VITESSE                  =4;4

COTES IMPOSEES                      =0.1

/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES
/---------------------------------------------------------------------

REMISE A ZERO DU TEMPS =NON

BANCS DECOUVRANTS      =NON

PAS DE TEMPS           =0.4

DUREE DU CALCUL        =100000

/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES, SOLVEUR
/---------------------------------------------------------------------

OPTION DU SOLVEUR    =3

SOLVEUR              =7

PRECISION DU SOLVEUR =1.E-4

/---------------------------------------------------------------------
/ PARAMETRES NUMERIQUES, VITESSE-CELERITE-HAUTEUR
/---------------------------------------------------------------------

MASS-LUMPING SUR LA VITESSE   =1

IMPLICITATION POUR LA HAUTEUR =0.55

MASS-LUMPING SUR H            =1

IMPLICITATION POUR LA VITESSE =0.55