!****************************************************
!**                                                **
!**      Auteur: Ambroise BROU                     ** 
!**      Cree en Aout 2005                         **
!**      Derniere mise a jour Novembre 2006        **
!**                                                **
!****************************************************



/title, Ecoulement autour d'un cylindre 

/PREP7 					! entree dans le mode preprocessing
/UNITS,SI				! systeme international des unites
R=0.05					! Rayon du cylindre en m
MU=1.8000e-5				! viscosite dynamique en kg/(m.s)
RHO=1.230				! masse volumique en kg/m3
RE=150					! Nombre de Reynolds
VE=RE*MU/(RHO*2*R)			! Vitesse a l'infini m/s
NR1=30					! nombre de division 
NR2=30					!
NX1=5*NR1				!
NX2=6
NX3=2
RX1=0.2


ET,1,FLUID141				! elements finis de type FLUID141

!coordonnees des points definissant le domaine 
K,1,-10*R,-20*R  		
K,2,-5*R,-20*R
K,3,5*R,-20*R 
K,4,40*R,-20*R
K,5,40*R,-5*R
K,6,40*R,5*R
K,7,40*R,20*R
K,8,5*R,20*R
K,9,-5*R,20*R
K,10,-10*R,20*R
K,11,-10*R,5*R
K,12,-10*R,-5*R
K,13,-5*R,-5*R
K,14,5*R,-5*R
K,15,5*R,5*R
K,16,-5*R,5*R
K,17,-R*(2**0.5)/2,R*(2**0.5)/2
K,18,-R*(2**0.5)/2,-R*(2**0.5)/2
K,19,R*(2**0.5)/2,-R*(2**0.5)/2
K,20,R*(2**0.5)/2,R*(2**0.5)/2
K,21,0,0				! centre du cercle

!creation des aires delimitees par les differents points cles
A,1,2,13,12
A,2,3,14,13
A,3,4,5,14
A,14,5,6,15
A,15,6,7,8
A,16,15,8,9
A,11,16,9,10
A,12,13,16,11

LARC,18,19,21,R				! creation d'un arc de cercle ayant pour centre le point 21, de rayon R reliant les points 18 et 19
LARC,19,20,21,R
LARC,20,17,21,R
LARC,17,18,21,R
L,16,17
L,13,18
L,14,19
L,15,20

! Creations des aires delimitees par les differentes lignes
AL,7,31,25,30
AL,13,26,31,32
AL,17,27,29,32
AL,23,28,29,30
SAVE

!maillage du domaine
LSEL,S,,,25,28
LESIZE,ALL,,,NR1,1
LSEL,S,,,7,17,10
LSEL,A,,,13,23,10
LSEL,A,,,11
LESIZE,ALL,,,NR1,1
LSEL,S,,,29,32,1
LESIZE,ALL,,,NR2,RX1
LSEL,S,,,10,12,2
LESIZE,ALL,,,NX1,1
LSEL,S,,,8,15,7
LESIZE,ALL,,,NX2,1
LSEL,S,,,5,18,13
LESIZE,ALL,,,NX3,1
LSEL,S,,,1,3,2
LSEL,A,,,20,21,1
LESIZE,ALL,,,NX3,1
SAVE					! sauvegarde de la base de donnees avant maillage

MSHKEY,1				
MSHAPE,0,2D				! maillage 2D mappe
AMESH,9,12,1				! maillage des surfaces 9 a 12
AMESH,4					! maillage de la surface 4
MSHKEY,0				! maillage libre
MSHAPE,1,2D				! choix d'un maillage triangulaire libre
AMESH,ALL				! maillage de toutes les surfaces non encore maillees

!Conditions limites a l'entree
LSEL,S,,,4				! selection de la ligne 4
LSEL,A,,,22,24,2			! selection additionnelle des lignes 22 et 24
NSLL,S,1				! selection des noeuds contenus dans les lignes selectionnees			
D,ALL,VX,VE				! vitesse dans la direction x egale a VE
D,ALL,VY,0				! vitesse dans la direction y nulle

!Conditions de non glissement sur le cylindre
LSEL,S,,,25,28				! selection des lignes 25 a 28
NSLL,S,1				! selection des noeuds contenus dans les lignes
D,ALL,VX,0				! vitesse dans la direction x nulle
D,ALL,VY,0				! vitesse dans la direction y nulle

!Conditions limites a la sortie
LSEL,S,,,9,11,2				! selection des lignes 9 et 12
LSEL,A,,,14				! selection additionnelle de la ligne 14
NSLL,S,1				! selection des noeuds contenus dans les lignes selectionnees
D,ALL,PRES,0				! pression relative nulle sur les noeuds selectionnes

!Conditions limites aux frontieres
LSEL,S,,,1,5,4				! selection des lignes 1 et 5
LSEL,A,,,8				! selection additionnelle de la ligne 8
LSEL,A,,,15,21,3			! selection additionnelle des lignes 15, 18 et 21
NSLL,S,1				! selection des noeuds contenus dans les lignes 
D,ALL,PRES,0				! pression relative nulle sur les noeuds selectionnes

ALLS					! selection de tout
SAVE                     		! sauvegarde de la base de donnees avant resolution
!Proprietes du fluide
FLDA,PROT,DENS,CONSTANT			! masse volumique constante
FLDA,NOMI,DENS,RHO			! masse volumique nominale
!FLDATA,PROT,DENS,AIR-SI
FLDATA,PROT,VISC,CONSTANT		! viscosite constante
FLDATA,NOMI,VISC,MU			! viscosite du fluide
FINISH					! sortie du mode preprocessing

/SOLU					! entree dans le mode solution
FLDATA,SOLU,TRAN,1                      ! ecoulement transitoire
FLDATA,SOLU,TURB,0			! laminaire
FLDATA,TIME,METH,NEWM			! Methode d'integration pour le pas de temps
FLDATA,ALGR,SEGR,SIMPLEN		! choix du type d'algorithme
FLDATA,TIME,STEP,-1			! pas de temps gouverne par l'advection
FLDATA,TIME,GLOB,30			! nombre d'iterations par pas de temps
FLDATA,TIME,NUMB,4000			! duree totale en seconde 
FLDATA,TIME,APPE,1			! ecriture des resultats dans le fichier .rfl toutes les secondes
FLDATA,RELX,Vx,0.95			! taux de relaxation
FLDATA,RELX,Vy,0.95
FLDATA,OUTP,PCOE,1 			! parametres supplementaies pour le post-traitement
FLDATA,OUTP,TAUW,1
SOLVE					! demarrage de la resolution
FINISH					! sortie du mode solution