Simulation numérique et stabilité des écoulements de convection mixte en conduite rectangulaire chauffée par le bas PDF Download
Are you looking for read ebook online? Search for your book and save it on your Kindle device, PC, phones or tablets. Download Simulation numérique et stabilité des écoulements de convection mixte en conduite rectangulaire chauffée par le bas PDF full book. Access full book title Simulation numérique et stabilité des écoulements de convection mixte en conduite rectangulaire chauffée par le bas by Xavier Nicolas. Download full books in PDF and EPUB format.
Author: Xavier Nicolas Publisher: ISBN: Category : Languages : en Pages : 0
Book Description
ON PROPOSE UNE ETUDE NUMERIQUE 2D DES ECOULEMENTS LAMINAIRES DE CONVECTION MIXTE EN CONDUITE RECTANGULAIRE HORIZONTALE CHAUFFEE PAR LE BAS (ECOULEMENTS DE POISEUILLE-BENARD). LA METHODE DE RESOLUTION EST DE TYPE VOLUMES FINIS ; ELLE EST BASEE SUR LA METHODE DU LAGRANGIEN AUGMENTE. ON PRESENTE UNE REVUE EXHAUSTIVE DES TRAVAUX EFFECTUES SUR CE SUJET. ON RAPPORTE DE MANIERE DETAILLEE LES RESULTATS DE TROIS THESES, ANALYSANT LA STABILITE LINEAIRE DE L'ECOULEMENT DE POISEUILLE ET LA STABILITE FAIBLEMENT NON-LINEAIRE DES ROULEAUX THERMOCONVECTIFS TRANSVERSES A L'AXE DE LA CONDUITE (RT), ET FOURNISSANT DES MESURES EXPERIMENTALES TRES FINES DES CHAMPS DE VITESSE DANS LES DOMAINES LINEAIRES ET NON-LINEAIRES. L'OBJECTIF DE CE TRAVAIL EST D'ETUDIER DANS QUELLE MESURE LES SIMULATIONS NUMERIQUES S'ACCORDENT AVEC THEORIES ET EXPERIENCES EN EFFECTUANT DES COMPARAISONS QUANTITATIVES SYSTEMATIQUES AVEC CES TROIS THESES. ON ETUDIE L'INFLUENCE DES CONDITIONS AUX LIMITES OUVERTES (CLO) ET PERIODIQUES (CLP) SUR LA STABILITE ET SUR LE DEVELOPPEMENT SPATIAL ET TEMPOREL DES RT, EN PRENANT EN COMPTE LE CARACTERE CONVECTIF OU ABSOLU DE L'INSTABILITE. ON TESTE CINQ SORTES DE CLO. ON MONTRE QUE L'AMPLITUDE DE LA PERTURBATION EN AMONT DE LA FRONTIERE DE SORTIE PEUT CONSIDERABLEMENT VARIER D'UNE CLO A L'AUTRE, ALORS QUE LA LONGUEUR DE LA ZONE PERTURBEE VARIE TRES PEU, ET QUE LA STRUCTURE DE L'ECOULEMENT EN DEHORS DE LA ZONE PERTURBEE SE CONSERVE. AVEC LES CLP, LES RT TRANSITENT VERS L'ECOULEMENT DE POISEUILLE AU PASSAGE DU SEUIL ENTRE INSTABILITE CONVECTIVE ET STABILITE LINEAIRE. AVEC LES CLO, CETTE TRANSITION A LIEU PRES DU SEUIL ENTRE INSTABILITE ABSOLUE ET INSTABILITE CONVECTIVE. TOUS LES RESULTATS DE LA STABILITE FAIBLEMENT NON-LINEAIRE SONT REPRODUITS PAR LA SIMULATION 2D ET CERTAINS RESULTATS EXPERIMENTAUX 3D SONT RETROUVES AVEC UNE TRES GRANDE PRECISION. ENFIN, ON MET EN EVIDENCE LES POINTS NECESSITANT DES APPROFONDISSEMENTS OU NE POUVANT ETRE TRAITES QUE PAR LA SIMULATION NUMERIQUE 3D.
Author: Xavier Nicolas Publisher: ISBN: Category : Languages : en Pages : 0
Book Description
ON PROPOSE UNE ETUDE NUMERIQUE 2D DES ECOULEMENTS LAMINAIRES DE CONVECTION MIXTE EN CONDUITE RECTANGULAIRE HORIZONTALE CHAUFFEE PAR LE BAS (ECOULEMENTS DE POISEUILLE-BENARD). LA METHODE DE RESOLUTION EST DE TYPE VOLUMES FINIS ; ELLE EST BASEE SUR LA METHODE DU LAGRANGIEN AUGMENTE. ON PRESENTE UNE REVUE EXHAUSTIVE DES TRAVAUX EFFECTUES SUR CE SUJET. ON RAPPORTE DE MANIERE DETAILLEE LES RESULTATS DE TROIS THESES, ANALYSANT LA STABILITE LINEAIRE DE L'ECOULEMENT DE POISEUILLE ET LA STABILITE FAIBLEMENT NON-LINEAIRE DES ROULEAUX THERMOCONVECTIFS TRANSVERSES A L'AXE DE LA CONDUITE (RT), ET FOURNISSANT DES MESURES EXPERIMENTALES TRES FINES DES CHAMPS DE VITESSE DANS LES DOMAINES LINEAIRES ET NON-LINEAIRES. L'OBJECTIF DE CE TRAVAIL EST D'ETUDIER DANS QUELLE MESURE LES SIMULATIONS NUMERIQUES S'ACCORDENT AVEC THEORIES ET EXPERIENCES EN EFFECTUANT DES COMPARAISONS QUANTITATIVES SYSTEMATIQUES AVEC CES TROIS THESES. ON ETUDIE L'INFLUENCE DES CONDITIONS AUX LIMITES OUVERTES (CLO) ET PERIODIQUES (CLP) SUR LA STABILITE ET SUR LE DEVELOPPEMENT SPATIAL ET TEMPOREL DES RT, EN PRENANT EN COMPTE LE CARACTERE CONVECTIF OU ABSOLU DE L'INSTABILITE. ON TESTE CINQ SORTES DE CLO. ON MONTRE QUE L'AMPLITUDE DE LA PERTURBATION EN AMONT DE LA FRONTIERE DE SORTIE PEUT CONSIDERABLEMENT VARIER D'UNE CLO A L'AUTRE, ALORS QUE LA LONGUEUR DE LA ZONE PERTURBEE VARIE TRES PEU, ET QUE LA STRUCTURE DE L'ECOULEMENT EN DEHORS DE LA ZONE PERTURBEE SE CONSERVE. AVEC LES CLP, LES RT TRANSITENT VERS L'ECOULEMENT DE POISEUILLE AU PASSAGE DU SEUIL ENTRE INSTABILITE CONVECTIVE ET STABILITE LINEAIRE. AVEC LES CLO, CETTE TRANSITION A LIEU PRES DU SEUIL ENTRE INSTABILITE ABSOLUE ET INSTABILITE CONVECTIVE. TOUS LES RESULTATS DE LA STABILITE FAIBLEMENT NON-LINEAIRE SONT REPRODUITS PAR LA SIMULATION 2D ET CERTAINS RESULTATS EXPERIMENTAUX 3D SONT RETROUVES AVEC UNE TRES GRANDE PRECISION. ENFIN, ON MET EN EVIDENCE LES POINTS NECESSITANT DES APPROFONDISSEMENTS OU NE POUVANT ETRE TRAITES QUE PAR LA SIMULATION NUMERIQUE 3D.
Book Description
Le modèle numérique que nous avons développé au cours de cette thèse présente deux caractéristiques principales : un modèle dilatable pour l'eau et la prise en compte de domaines ouverts. Les difficultés associées au premier aspect concernent l'adaptation de la loi d'état de l'eau au modèle dilatable sous l'approximation à faibles nombres de Mach, tandis que celles associées au second sont relatives à la mise en œuvre de conditions aux limites numériques de sortie compatibles avec l'algorithme de projection utilisé. Les résultats de simulations d'écoulement de convection mixte en canal horizontal chauffé par le bas ont été confrontés à celles utilisant l'approximation de Boussinesq et aux expériences.
Author: JEAN.. EL HAJAL Publisher: ISBN: Category : Languages : fr Pages : 187
Book Description
CE TRAVAIL CONCERNE L'ETUDE NUMERIQUE ET EXPERIMENTALE DU COUPLAGE THERMO-ELECTROHYDRODYNAMIQUE DANS UN SYSTEME DE CHAUFFAGE VOLUMIQUE DE LIQUIDES IONIQUES PAR CONDUCTION ELECTRIQUE DIRECTE CHAUFFAGE OHMIQUE. IL S'AGIT DE CHAUFFER UN FLUIDE PAR EFFET JOULE EN UTILISANT LE PRODUIT MEME COMME RESISTANCE ELECTRIQUE. LE PRINCIPE AINSI QUE LES NOMBREUX AVANTAGES ET APPLICATIONS DE CE PROCEDE SONT PRESENTES. LA CONFIGURATION DE L'ETUDE CONSISTE EN UN ECOULEMENT LAMINAIRE DANS UN CANAL VERTICAL A SECTION RECTANGULAIRE AVEC DISSIPATION VOLUMIQUE D'ENERGIE. UN MONTAGE EXPERIMENTAL A PERMIS L'EXPLORATION DES CHAMPS DE VITESSE ET DE TEMPERATURE. DIFFERENTES METHODES EXPERIMENTALES ONT ETE MISES EN UVRE : LA VELOCIMETRIE LASER A EFFET DOPPLER, LA VISUALISATION PAR FLUORESCENCE INDUITE PAR LASER ET LA MESURE DE TEMPERATURE PAR THERMOCOUPLES. LES RESULTATS ONT PERMIS D'ANALYSER LE COUPLAGE DES PHENOMENES PHYSIQUES MIS EN JEU. CES RESULTATS ONT EGALEMENT SERVI A LA VALIDATION D'UNE VERSION SPECIFIQUE D'UN CODE DE CALCUL 3D EN VOLUMES FINIS, ESTET. L'ETUDE NUMERIQUE A PERMIS D'ABORDER DIFFERENTS ASPECTS DU PROBLEME. L'INFLUENCE DES PARAMETRES, TELS QUE LA VARIATION DES PROPRIETES PHYSIQUES, LA DENSITE DE CHAUFFAGE VOLUMIQUE ET LES CONDITIONS AUX LIMITES, A ETE ANALYSEE. L'UNIFORMITE DU CHAUFFAGE A ETE CARACTERISEE. UNE COMPARAISON AVEC UN CHAUFFAGE PARIETAL A FLUX CONSTANT A PERMIS DE METTRE EN EVIDENCE CETTE HOMOGENEITE DE LA TEMPERATURE PLUS IMPORTANTE DANS LE CAS DU CHAUFFAGE VOLUMIQUE. UNE PREMIERE PHASE D'UNE ETUDE CONCERNANT LE CONTROLE ACTIF D'UN TEL SYSTEME A ETE REALISEE. LE BUT FINAL CONSISTE A DETERMINER UNE BOUCLE DE CONTROLE PERMETTANT D'OPTIMISER LE SYSTEME SELON DES CRITERES PROPRES A L'APPLICATION VISEE. UNE ANALYSE LINEAIRE DE STABILITE DANS LE CAS D'ECOULEMENT FORCE EN REGIME ETABLI A ETE REALISEE. LA REPONSE D'UN TEL ECOULEMENT A UN ECHELON DE CHAMP ELECTRIQUE A ETE ETUDIEE GRACE A UNE SIMULATION NUMERIQUE. LA CONVECTION NATURELLE EN CAVITE SOUS L'EFFET DU CHAMP ELECTRIQUE A ETE ETUDIEE. IL S'AGIT D'ETUDIER LA CONVECTION NATURELLE AVEC CHAUFFAGE VOLUMIQUE. AINSI QUE L'EFFET EVENTUEL DE LA FORCE ELECTRIQUE DANS CE SYSTEME. UN MONTAGE EXPERIMENTAL ET UN MODELE NUMERIQUE ONT ETE DEVELOPPES. LA STRUCTURE DE L'ECOULEMENT ET DES REGIMES QUI APPARAISSENT DANS CETTE CONFIGURATION SONT CARACTERISES.
Book Description
L'essentiel de mon travail de recherche porte sur l'étude des instabilités thermoconvectives en convection mixte. Il comporte aussi bien des aspects expérimentaux que numériques.Dès qu'un fluide est chauffé, en présence ou non d'un écoulement forcé, des instabilités apparaissent en fonction des différents paramètres (thermiques, opératoires, géométrie,..). Pour ce qui est des dispositifs parcourus par un écoulement de convection mixte, la bibliographie est très riche concernant les écoulements de fluides tel que l'air mais reste limitée pour tout ce qui est écoulement d'eau et surtout dans les configurations avec une condition aux limites de flux imposé. Récemment, ces études ont suscité un intérêt pour les applications industrielles telles que la CVD et les échangeurs de chaleur en général et certains phénomènes naturels tels que la météo. Faisant suite à une série d'études d'écoulements de convection mixte, notamment en conduite cylindrique soumise à un flux uniforme à la paroi (coexistence de gradients thermiques horizontal et vertical), rectangulaire avec un gradient thermique horizontal et rectangulaire avec un gradient thermique vertical ; dans ce dernier cas, le rapport d'aspect considéré était égal à 2 ce qui a introduit des effets de bords non négligeables sur le développement des instabilités. L'objectif de mon travail était de comprendre le déclenchement des instabilités dans le cas du gradient thermique vertical sans subir les effets de bords. Ainsi, j'ai étudié l'écoulement de convection mixte dans un canal horizontal, à section rectangulaire, de grand rapport d'aspect transversal (B=10), chauffé uniformément par le bas à flux constant pour les paramètres de fonctionnement suivants : 20 ≤ Re ≤ 200, 0 ≤ Ra ≤ 106, Pr = 7. Les techniques de visualisation que j'ai mises en place m'ont permis de mettre en évidence l'existence de deux mécanismes d'initiation des rouleaux longitudinaux dans ce type de configuration. Pour les faibles valeurs du nombre de Rayleigh, les rouleaux sont initiés le long des parois latérales et se propagent progressivement jusqu'à envahir toute la section droite par un effet de cascade visqueuse. Au-delà d'une valeur critique du nombre de Rayleigh, de nombreuses paires de rouleaux apparaissent simultanément à une même abscisse, dans la zone centrale de la section droite, indépendamment du premier mécanisme grâce à une instabilité du type Rayleigh Bénard. Ces deux mécanismes d'initiation ont été observés expérimentalement et numériquement. La condition de flux imposé à la paroi basse fait que l'établissement thermique n'est jamais atteint ; nous assistons dans ce cas à une évolution longitudinale du régime d'écoulement qui peut s'étaler le long du canal. Ainsi, j'ai mis en évidence quatre zones d'écoulement : une zone I d'entrée ou de formation des rouleaux ; une zone II établie où les rouleaux longitudinaux ont envahi toute la section droite, ces rouleaux sont stationnaires ; une zone III correspondant au régime déstabilisé où les rouleaux adoptent un comportement chaotique et enfin la zone IV qui correspond au régime turbulent, les rouleaux ne sont plus identifiables. J'ai déterminé en fonction des paramètres de contrôle les longueurs caractéristiques pour les différentes zones et proposé des corrélations permettant de les calculer en fonction des nombres de Reynolds et de Rayleigh. Pour l'analyse thermique de ce type d'écoulement, j'ai développé et exploité une technique de fluorescence induite par plan laser (PLIF) pour déterminer le champ thermique dans les différentes sections transversales du canal pour les différentes gammes des paramètres. Cette caractérisation a également été effectuée numériquement. Les champs de température obtenus dans les différentes sections droites fluide ont permis la détermination du transfert de chaleur, notamment le nombre de Nusselt. J'ai déterminé le nombre de Nusselt à la paroi chauffée, moyenné dans la direction transversale Y et représenté sa variation selon la direction longitudinale X pour diverses valeurs du nombre de Rayleigh et de Reynolds. On retrouve la signature des deux mécanismes d'initiation au travers du profil longitudinal du nombre de Nusselt. On peut observer que le transfert convectif est fortement accru en convection mixte par rapport à la convection forcée pure, d'un rapport de 5 à 10 selon les valeurs des nombres de Reynolds et de Rayleigh. L'ensemble des résultats expérimentaux que j'ai obtenus ont permis, en fonction des nombres de Reynolds et de Rayleigh, l'établissement de corrélations du nombre de Nusselt, valables en convection mixte pour ce type d'écoulement. Ces corrélations ont été établies pour le régime d'entrée et le régime établi. La PLIF semble être une technique performante pour la mesure de la température du fluide. L'ensemble de ces résultats ont été rédigés en deux articles et ont été soumis pour publication.
Book Description
Le travail présenté porte sur l'étude de la convection mixte en fluide binaire avec effet Soret (gradient de concentration induit par un gradient de température). La géométrie considérée est une conduite parallélépipédique, horizontale, chauffée par le bas. L'écoulement moyen est généré par un gradient de pression horizontal. Nous nous intéresserons à la naissance de structures thermoconvectives sous la forme de rouleaux transversaux 2D, perpendiculaires à l'axe de la conduite. Après avoir introduit les différentes notions de base et effectué une revue bibliographique, nous présentons le formalisme mathématique adopté. L'étude de la convection mixte en fluide binaire avec effet Soret est ensuite décomposée en trois parties : théorique, numérique et expérimentale. La partie théorique consiste en l'étude de la stabilité linéaire de l'écoulement de Poiseuille, linéairement stratifié en température et en concentration. Pour un mélange caractérisé par un nombre de Lewis élevé, un facteur de séparation négatif et pour de faibles valeurs du nombre de Reynols, nous retrouvons conformément à la littérature, à la naissance de la convection, la présence de deux branches de solutions correspondant à des rouleaux transversaux 2D se propageant dans le sens de l'écoulement moyen ou en sens inverse. D'un point de vue numérique, nous avons utilisé un code spectral pour détecter la naissance des rouleaux transversaux 2D. Nous avons également effectué des simulations numériques 2D et 3D avec le logiciel FEMLAB. Ces simulations nous ont permis de mettre en évidence les deux branches de solutions distinctes à la naissance de la convection mixte avec effet Soret, déterminées par l'analyse de stabilité linéaire. Enfin, nous avons mis au point un dispositif expérimental de PIV permettant la visualisation et la détermination du champ de vitesses des différentes structures thermoconvectives observables à la naissance de la convection dans des mélanges à facteur de séparation négatif.
Book Description
On présente les résultats d’une simulation numérique de la convection mixte conjuguée transitoire dans un tube vertical soumis à un flux DED chaleur constant et uniforme. Le fluide pénètre au haut du tube pour se diriger vers le bas ; par conséquent on est en présence d’un écoulement de convection mixte contrariée. Les équations de conservation dans le fluide et dans la paroi sont résolues numériquement en utilisant la méthode des volumes finis. Le couplage pression-vitesse est assuré en utilisant l’algorithme simple. On étudie l’influence des propriétés physiques et géométriques sur l’évolution transitoire des grandeurs thermiques (flux de chaleur à l’interface paroi-fluide et distribution radiale de température) et sur les grandeurs dynamiques (coefficient de frottement et champs vecteurs vitesses).
Book Description
LES TRANSFERTS THERMO-CONVECTIFS ASSOCIES A LA CLIMATISATION D'UN LOCAL DE GRAND VOLUME SONT ETUDIES NUMERIQUEMENT. LE MODELE DE TURBULENCE - A BAS NOMBRE DE REYNOLDS DE LAUNDER ET SHARMA A ETE INTEGRE DANS LE CODE DE CALCUL PHOENICS. UNE ETUDE PHENOMENOLOGIQUE DE L'INFLUENCE DES DIFFERENTS PARAMETRES CARACTERISTIQUES DE LA VENTILATION D'UN LOCAL SIMPLE DE GRAND VOLUME EN CONDITIONNEMENT D'AIR CHAUFFE OU REFROIDI A ETE MENEE. L'ETUDE DES PRINCIPAUX ECOULEMENTS QUI INTERAGISSENT DANS UN LOCAL DE GRAND VOLUME A ETE DECOUPLEE POUR EVALUER LES PERFORMANCES DU MODELE DE TURBULENCE A FAIBLE NOMBRE DE REYNOLDS: D'UNE PART, DIFFERENTS ECOULEMENTS DE CONVECTION MIXTE DE JET ET, D'AUTRE PART, UNE CONFIGURATION DE CONVECTION NATURELLE REPRESENTEE PAR UNE CAVITE DIFFERENTIELLEMENT CHAUFFEE. DANS LE CAS DE LA CAVITE DIFFERENTIELLEMENT CHAUFFEE DE RAPPORT DE FORME UNITE ET UN NOMBRE DE RAYLEIGH DE 10#1#0, LA COMPARAISON DES DIFFERENTS MODELES DE TURBULENCE - MONTRE UN ACCORD GLOBALEMENT SATISFAISANT. LA VERSION A BAS NOMBRE DE REYNOLDS AMELIORE LES RESULTATS LORSQUE LA SIMULATION DIRECTE DE LA TURBULENCE EST PRISE POUR REFERENCE. L'ETUDE DES JETS MONTRE LE BON COMPORTEMENT DES MODELES DE TYPE - POUR LES ECOULEMENTS DE JETS ISOTHERMES ET LES PANACHES EN AMBIANCE NEUTRE. DANS LE CAS DES JETS SOUFFLES VERTICALEMENT DANS UNE AMBIANCE STABLE LINEAIREMENT STRATIFIEE, UN BON ACCORD EST TROUVE AVEC DES RESULTATS EXPERIMENTAUX SUR LA HAUTEUR DE CULMINATION DANS LA GAMME DE GRAND NOMBRE DE FROUDE ET FORTE STRATIFICATION. L'APPLICATION DU MODELE DE TURBULENCE A LA CLIMATISATION DES LOCAUX DE GRAND VOLUME DONNE DES STRUCTURES D'ECOULEMENT TRES DIFFERENTES SELON LES VALEURS DES PARAMETRES CARACTERISTIQUES TELS QUE LE DEBIT DE SOUFFLAGE, L'ANGLE SOLIDE D'INJECTION DE L'AIR, LE POSITIONNEMENT ET L'ESPACEMENT D'UNE DISTRIBUTION REGULIERE DES AEROTHERMES AINSI QUE LA TEMPERATURE DE L'AIR VENTILE
Book Description
Ce travail traite le problμeme de la convection mixte dans un canal rectangulaire horizontal soumis a un gradient horizontal de temperature. Ce gradient cree un ecoulement de convection naturelle qui se superpose μa l'ecoulement principal ; ce qui modifie sa structure initiale. A faibles Reynolds et pour des Rayleigh moderes, on a un rouleau longitudinal montrant une trajectoire helicoidale du fluide. Pour des Rayleigh plus eleves, on observe un ecoulement de retour qui vient se superposer au rouleau longitudinal. Leur mise en evidence est faite μa l'aide d'outils experimental et numerique et la determination du coefficient de transfert montre un meilleur echange convectif. Pour les Reynolds grands, on observe un regime d'instabilites qui tend vers une turbulence hydrodynamique ou thermique. L'absence d'un regime intermittent prouve que le gradient horizontal de temperature est stabilisant vis-a-vis de l'instabilite thermoconvective observee dans le cas cylindrique.
Book Description
Dans de nombreux dispositifs industriels comme les conduites de refroidissement des moteurs de fusées, un écoulement turbulent est en contact avec une paroi solide chaude.L'instabilité due aux courbures ajoutée au confinement créent des courants orthogonaux au flux principal d'une intensité pouvant aller jusqu'à 20% de celui-ci. Pour comprendre les interactions entre ces courants secondaires et les transferts thermiques, des résultats de simulations numériques des grandes échelles sont obtenus dans des conduites pouvant présenter des courbures et des parois chauffées. Un écoulement pleinement développé est injecté à l'entrée du domaine par la simulation synchrone d'une conduite périodique. Un modèle de viscosité turbulente a été implanté dans le code de calcul et aucune loi de paroi n'est utilisée. Une température supérieure est d'abord appliquée sur une des quatre parois d'une conduite carrée rectiligne pour une validation du code de calcul. Les mécanismes d'éjections de fluide chaud de la paroi chauffée vers le coeur de la conduite sont bien représentés et on peut observer le développement progressif de la couche limite thermique. L'évolution spatiale de l'écoulement a ensuite permis d'appliquer un flux constant.Une conduite de section rectangulaire en forme de S a ensuite été simulée. A chaque courbure, apparaît un fort flux secondaire inhérent au gradient de pression normal aux parois courbes. Ceci implique la naissance d'une paire de tourbillons longitudinaux contra-rotatifs de type "Dean" ou "Ekman" évoluant au voisinage de la paroi convexe.L'interaction entre ces courants secondaires et le mélange turbulent est étudiée suivant les modes de chauffage.
Author: Xavier Nicolas
Author: Xavier Nicolas
Author: 
Author: Adel Haddad
Author: JEAN.. EL HAJAL
Author: Rani Taher
Author: Estelle Piquer
Author: Abdeslam Omara
Author: PHILIPPE.. BUCHMANN
Author: Kouadio Firmin Koffi
Author: Jérôme Hébrard