Propagation des ondes acoustiques dans les milieux en mouvement 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 Propagation des ondes acoustiques dans les milieux en mouvement PDF full book. Access full book title Propagation des ondes acoustiques dans les milieux en mouvement by Laurent Dallois. Download full books in PDF and EPUB format.
Author: Laurent Dallois
Publisher:
ISBN:
Category :
Languages : fr
Pages : 152
Get Book Here
Book Description
De nos jours, la propagation du son dans les milieux inhomogènes est largement modélisée à partir d'une approximation paraxiale de l'équation d'onde appelée équation parabolique. Si dans le cadre de la propagation en milieu au repos il est possible d'obtenir facilement cette équation parabolique, l'extension au milieu en mouvement reste un problème ouvert. Jusqu'à maintenant, les équations paraboliques développées pour prendre en compte un champ de vitesse n'utilisaient que la seule composante de vitesse dans la direction de propagation de l'onde acoustique. Deux nouvelles équations paraboliques grand angle qui prennent en compte la nature vectorielle du champ de vitesse ont été dérivées, et chacune adaptée a un type d'écoulement particulier : l'équation mw-wape a été dérivée a partir d'une équation d'onde exacte quant à la prise en compte des effets de convection du milieu. Elle est utilisée pour traiter le cas de grandes structures d'écoulements (écoulement moyen). Elle est d'ordre deux par rapport au nombre de mach. L'équation tw-wape a été dérivée à partir d'une équation d'onde plus générale pour la propagation dans un milieu inhomogène en mouvement. Elle est principalement développée pour modéliser la propagation en écoulements turbulents. Elle est d'ordre un en nombre de mach et intègre les termes de gradient de l'écoulement. Nous avons développé des schémas numériques adaptés à des cas bidimensionnels de propagation. Une première étape, essentiellement de validation, consiste à traiter le cas de la propagation du son au dessus d'un sol plan pour un milieu non turbulent. On obtient un bon accord avec les solutions de référence. Dans une deuxième étape, on étudie des géométries plus complexes : la propagation d'une onde plane à travers un vortex, et la diffraction du son dans une zone d'ombre acoustique par une turbulence cinématique. On montre la validité de l'équation parabolique mw-wape pour des nombres de mach allant jusqu'a 0.5. Dans le cas de la zone d'ombre, nous obtenons des modifications du niveau acoustique moyen de quelques db en utilisant l'équation parabolique tw-wape par rapport aux résultats d'une équation parabolique standard. Enfin, nous appliquons a la résolution de l'équation parabolique tw-wape deux méthodes permettant d'augmenter son domaine d'application et de diminuer les besoins en ressources informatiques : la méthode de résolution de type split-step pade d'ordre (n,n) et la méthode des écrans de phases. On accroit ainsi l'angle de validité de l'équation parabolique et on a une meilleure résolution des effets des petites structures de notre champ turbulent. De plus, ces deux méthodes permettent d'envisager le passage à des résolutions tridimensionnelles pour prendre en compte les conditions d'une couche limite atmosphérique par exemple. Ces différents résultats prouvent que l'utilisation d'un indice effectif pour la modélisation des effets du champ de vitesse du milieu doit être limitée a des cas simples. Dans le cas général, on doit utiliser des équations paraboliques du type de celles développées dans cette thèse.
Author: Laurent Dallois
Publisher:
ISBN:
Category :
Languages : fr
Pages : 152
Get Book Here
Book Description
De nos jours, la propagation du son dans les milieux inhomogènes est largement modélisée à partir d'une approximation paraxiale de l'équation d'onde appelée équation parabolique. Si dans le cadre de la propagation en milieu au repos il est possible d'obtenir facilement cette équation parabolique, l'extension au milieu en mouvement reste un problème ouvert. Jusqu'à maintenant, les équations paraboliques développées pour prendre en compte un champ de vitesse n'utilisaient que la seule composante de vitesse dans la direction de propagation de l'onde acoustique. Deux nouvelles équations paraboliques grand angle qui prennent en compte la nature vectorielle du champ de vitesse ont été dérivées, et chacune adaptée a un type d'écoulement particulier : l'équation mw-wape a été dérivée a partir d'une équation d'onde exacte quant à la prise en compte des effets de convection du milieu. Elle est utilisée pour traiter le cas de grandes structures d'écoulements (écoulement moyen). Elle est d'ordre deux par rapport au nombre de mach. L'équation tw-wape a été dérivée à partir d'une équation d'onde plus générale pour la propagation dans un milieu inhomogène en mouvement. Elle est principalement développée pour modéliser la propagation en écoulements turbulents. Elle est d'ordre un en nombre de mach et intègre les termes de gradient de l'écoulement. Nous avons développé des schémas numériques adaptés à des cas bidimensionnels de propagation. Une première étape, essentiellement de validation, consiste à traiter le cas de la propagation du son au dessus d'un sol plan pour un milieu non turbulent. On obtient un bon accord avec les solutions de référence. Dans une deuxième étape, on étudie des géométries plus complexes : la propagation d'une onde plane à travers un vortex, et la diffraction du son dans une zone d'ombre acoustique par une turbulence cinématique. On montre la validité de l'équation parabolique mw-wape pour des nombres de mach allant jusqu'a 0.5. Dans le cas de la zone d'ombre, nous obtenons des modifications du niveau acoustique moyen de quelques db en utilisant l'équation parabolique tw-wape par rapport aux résultats d'une équation parabolique standard. Enfin, nous appliquons a la résolution de l'équation parabolique tw-wape deux méthodes permettant d'augmenter son domaine d'application et de diminuer les besoins en ressources informatiques : la méthode de résolution de type split-step pade d'ordre (n,n) et la méthode des écrans de phases. On accroit ainsi l'angle de validité de l'équation parabolique et on a une meilleure résolution des effets des petites structures de notre champ turbulent. De plus, ces deux méthodes permettent d'envisager le passage à des résolutions tridimensionnelles pour prendre en compte les conditions d'une couche limite atmosphérique par exemple. Ces différents résultats prouvent que l'utilisation d'un indice effectif pour la modélisation des effets du champ de vitesse du milieu doit être limitée a des cas simples. Dans le cas général, on doit utiliser des équations paraboliques du type de celles développées dans cette thèse.
Author: Iryna Malinouskaya
Publisher:
ISBN:
Category :
Languages : en
Pages : 0
Get Book Here
Book Description
Author: Patrick Rasolofosaon
Publisher:
ISBN:
Category : Fluid dynamics
Languages : fr
Pages : 602
Get Book Here
Book Description
UNE PREMIERE PARTIE THEORIQUE EXPOSE EN DETAIL LA PROPAGATION DES ONDES ACOUSTIQUES DANS LES MILIEUX POREUX A TRAVERS LE MODELE POROVISCOELASTIQUE. ON MONTRE QUE LES EFFETS POROELASTIQUES NE SONT SUPERIEURS AUX EFFETS VISCOELASTIQUES QUI EN PRESENCE DE DISCONTINUITES PHYSIQUES OU D'ECHANGE HYDRAULIQUE ENTRE DIFFERENTS MILIEUX. TROIS EXPERIENCES DANS LA DEUXIEME PARTIE, CELLE DE PLONA, RESONANCE DES MODES DE LAMB, PROPAGATION D'ONDES A L'INTERFACE ROCHE FLUIDE. DANS LA TROISIEME PARTIE, ON ENVISAGE L'INVERSION DES MESURES ACOUSTIQUES
Author: LIAN-JIE.. HUANG
Publisher:
ISBN:
Category :
Languages : fr
Pages : 206
Get Book Here
Book Description
LES EQUATIONS D'ONDES ACOUSTIQUES OU ELASTIQUES, C'EST A DIRE LES EQUATIONS AUX DERIVEES PARTIELLES QUI DECRIVENT LE COMPORTEMENT DES ONDES DANS LES MILIEUX CONTINUS SONT EN GENERAL RESOLUES NUMERIQUEMENT POUR UN CHAMP D'ONDES CONTINU A INTERVALLES DE TEMPS DISCRETS, DANS UN ESPACE DISCRET, EN UTILISANT DES EQUATIONS AUX DIFFERENCES FINIES, POUR ETUDIER LES PHENOMENES ASSOCIES AUX ONDES PAR UNE SIMULATION SUR ORDINATEUR. IL EST DIFFICILE POUR LES METHODES DE DIFFERENCES FINIES DE SIMULER LA PROPAGATION D'ONDES DANS UN MILIEU FORTEMENT HETEROGENE ET COMPLEXE, A CAUSE DE LA COMPLEXITE GEOMETRIQUE ET/OU A CAUSE DES INTERFACES A FORT CONTRASTES DANS UN MILIEU AUX PROPRIETES PHYSIQUES VARIEES. LE PRINCIPAL OBJECTIF DE CETTE THESE EST DE DEVELOPPER DES METHODES MICROSCOPIQUES DANS UN ESPACE DISCRET, AVEC UN TEMPS ET DES CHAMPS D'ONDE EUX-MEMES DISCRETISES. GRACE A CES METHODES, LA SIMULATION DES PHENOMENES D'ONDES DANS DES MILIEUX COMPLEXES SANS APPROXIMATIONS ESSENTIELLES EST POSSIBLE. COMME PREMIERE APPROCHE A LA PROPAGATION DES ONDES P DANS UN MILIEU HETEROGENE, LE SOLIDE PHONONIQUE SUR RESEAU PAR INTERPOLATION EST DEVELOPPE DANS LE CHAPITRE 2 DANS LE BUT DE MODELISER DIRECTEMENT LE COMPORTEMENT DES NOMBRES DE DENSITE DE QUASI-PARTICULES (QUI TRANSPORTENT LA PRESSION) SUR UN MAILLAGE DISCRET, PLUTOT QUE DE RESOUDRE L'EQUATION D'ONDE MACROSCOPIQUE CORRESPONDANTE PAR DIFFERENCES FINIES. L'ANALYSE THEORIQUE DEMONTRE QUE DANS LA LIMITE MACROSCOPIQUE, LA METHODE MODELISE L'EQUATION D'ONDES ACOUSTIQUES POUR UN MILIEU HETEROGENE A N-DIMENSIONS (N=1, 2, 3). ELLE EST COMPARABLE AUX GAZ SUR RESEAU DANS LESQUELS LES ERREURS DUES AUX DIFFERENCES FINIES SONT ABSENTES, MAIS MODELISE LES PHENOMENES D'ONDES DANS DES MILIEUX HETEROGENES PLUTOT QU'UN FLUX DE FLUIDE ET D'ONDES ACOUSTIQUES DANS UN GAZ A MODULE DE COMPRESSIBILITE CONSTANTE. UNE CONDITION DE BORD ABSORBANT ET QUATRE MODELES DE REFLEXION SUR SURFACE LIBRE UTILISES DANS LE SOLIDE PHONONIQUE SUR RESEAU PAR INTERPOLATION SONT P ROPOSES DANS LE CHAPITRE 3. PUISQUE LA METHODE PERMET DE TENIR COMPTE D'INTERFACES A FORT CONTRASTES EN UN POINT QUELCONQUE DE LA GRILLE, CELLE-CI A ETE UTILISEE POUR SIMULER LA PROPAGATION D'ONDES DANS UN MILIEU FINEMENT STRATIFIE ET FRACTURE AUSSI BIEN QUE DANS UN MILIEU AVEC DES PORES VIDES. LES PROPRIETES D'ANISOTROPIE D'UN MILIEU FINEMENT STRATIFIE ET D'UN MILIEU FRACTURE ISOTROPIQUE ONT AINSI ETE OBSERVEES. CES RESULTATS SONT DECRITS AU CHAPITRE 4. LES MOUVEMENTS DES SITES DE LA GRILLE (C'EST A DIRE LES DEFORMATIONS) CAUSES PAR LE PASSAGE D'UNE ONDE MACROSCOPIQUE SONT PARTICULIEREMENT IMPORTANTS POUR UN MILIEU DONT LES PORES SONT REMPLIS DE FLUIDE, ETANT DONNE QUE LES INTERACTIONS NON-LINEAIRES ENTRE SOLIDES ET FLUIDES SONT SUPPOSEES JOUER UN ROLE IMPORTANT DANS LES MECANISMES D'ATTENUATION. LES MOUVEMENTS DES SITES DE LA GRILLE ONT ETE PRIS EN COMPTE POUR LE SOLIDE PHONONIQUE SUR RESEAU PAR INTERPOLATION ET L'APPROCHE EST APPELEE LE SOLIDE PHONONIQUE SUR RESEAU AVEC FLUIDES, PARCE QU'IL POURRAIT APPORTER UNE MEILLEURE COMPREHENSION DE L'EFFET DE L'INTERACTION ENTRE FLUIDE ET SOLIDE DANS LES PROBLEMES DE PROPAGATION D'ONDES. LA LIMITE MACROSCOPIQUE DE L'EQUATION DE TRANSPORT DU TYPE BOLTZMANN POUR CETTE METHODE MICROSCOPIQUE ABOUTIT A L'EQUATION D'ONDES ACOUSTIQUES POUR UN MILIEU HETEROGENE, MODIFIEE PAR DES TERMES DE CISAILLEMENT, DE MODULE DE VISCOSITE AINSI QUE PAR LE TERME DE SECOND ORDRE DANS LA VITESSE MACROSCOPIQUE ET DES TERMES NON-LINEAIRES ADDITIONNELS DUS AUX MOUVEMENTS DES SITES DE LA GRILLE. LA METHODE A ETE UTILISEE POUR ETUDIER L'ATTENUATION DE LA PROPAGATION D'ONDES DANS UN MILIEU AVEC DES PORES REMPLIS DE FLUIDE, ATTENUATION PROVOQUEE PAR LA VISCOSITE DES FLUIDES, PAR LA DIFFUSION DES HETEROGENEITES ET LE MOUVEMENT DES PARTICULES DU MILIEU INDUIT PAR LE PASSAGE D'UNE ONDE MACROSCOPIQUE. LA THEORIE QUI CONCERNE CETTE METHODE ET LES RESULTATS NUMERIQUES SONT DECRITS AU CHAPITRE 5. LE CHAPITRE 6 ABORDE LES PROBLEMES DE LA LIMITATION DE LA RESOLUTION ET DE LA DISTORSION DES IMAGES OBTENUES PAR UNE TOMOGRAPHIE A DIFFRACTION MONO-FREQUENTIELLE. LA TOMOGRAPHIE A RETRO-DIFFUSION MULTIFREQUENTIELLE SUR UN FOND QUI VARIE VERTICALEMENT A ETE DEVELOPPEE POUR AMELIORER HAUTEMENT LA RESOLUTION ET LA QUALITE DES IMAGES. LA VITESSE DES CALCULS DE LA METHODE EST DE LOIN SUPERIEURE A CELLE DES AUTRES METHODES MULTIFREQUENTIELLES, COMME L'HOLOGRAPHIE MULTIFREQUENTIELLE OU LA MIGRATION AVANT-SOMMATION
Author: Amine Berbiche
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
La méthode de minimisation de l'intégrale d'action (principe variationnel) permet d'obtenir les équations de propagation des ondes. Cette méthode a été généralisée aux milieux poreux de dimensions fractales, pour étudier la propagation acoustique dans le domaine temporel, en se basant sur le modèle du fluide équivalent. L'équation obtenue réécrite dans le domaine fréquentiel représente une généralisation de l'équation d'Helmholtz. Dans le cadre du modèle d'Allard-Johnson, l'équation de propagation a été résolue de manière analytique dans le domaine temporel, dans les régimes des hautes et des basses fréquences. La résolution a été faite par la méthode de la transformée de Laplace, et a porté sur un milieu poreux semi-infini. Il a été trouvé que la vitesse de propagation dépend de la dimension fractale. Pour un matériau poreux fractal d'épaisseur finie qui reçoit une onde acoustique en incidence normale, les conditions d'Euler ont été utilisées pour déterminer les champs réfléchi et transmis. La résolution du problème direct a été faite dans le domaine temporel, par la méthode de la transformée de Laplace, et par l'usage des fonctions de Mittag-Leffler. Le problème inverse a été résolu par la méthode de minimisation aux sens des moindres carrés. Des tests ont été effectués avec succès sur des données expérimentales, en utilisant des programmes numériques développés à partir du formalisme établi dans cette thèse. La résolution du problème inverse a permis de retrouver les paramètres acoustiques de mousses poreuses, dans les régimes des hautes et des basses fréquences.
Author: Emmanuel Redon
Publisher:
ISBN:
Category :
Languages : fr
Pages : 299
Get Book Here
Book Description
LES EQUATIONS GENERALES DE LA MECANIQUE DES FLUIDES SONT RAPPELEES AFIN D'ETABLIR LES EQUATIONS QUI GOUVERNENT LA PROPAGATION ACOUSTIQUE DANS LES MILIEUX NON ISOTHERMES EN MOUVEMENT. L'ETUDE BIBLIOGRAPHIQUE DES TRAVAUX REALISES SUR LE SUJET NOUS CONDUIT A TRAITER LE PROBLEME DIRECTEMENT A PARTIR DES EQUATIONS DE L'ACOUSTIQUE LINEAIRE. LES CHAMPS MOYENS DE VITESSE ET DE TEMPERATURE, QUI SONT DES PARAMETRES DU PROBLEME ACOUSTIQUE, SONT ENSUITE DECRITS PAR DES PROFILS ANALYTIQUES SIMPLES OU DETERMINES NUMERIQUEMENT A L'AIDE D'UN LOGICIEL DE MECANIQUE DES FLUIDES. LA MODELISATION DU COMPORTEMENT DES ONDES ACOUSTIQUES DANS LES MILIEUX CONFINES EN PRESENCE D'ECOULEMENT NON ISOTHERME EST DEVELOPPEE EN UTILISANT LA METHODE DES ELEMENTS FINIS. ON PRESENTE LES PRINCIPALES ETAPES DE LA DISCRETISATION QUI PERMETTENT D'OBTENIR LA FORME MATRICIELLE DU PROBLEME. SUR CETTE BASE, UN CODE DE CALCUL A ETE DEVELOPPE A PARTIR DE LA FORMULATION VITESSE-PRESSION-MASSE VOLUMIQUE (UVPR) DES EQUATIONS PRIMITIVES DE L'ACOUSTIQUE. LA DETERMINATION DES MODES PROPRES D'UNE CAVITE SANS ECOULEMENT PERMET DE VALIDER LA METHODE. LES RESULTATS NUMERIQUES SONT COMPARES AVEC LA THEORIE POUR LE CAS ISOTHERME ET AVEC CEUX OBTENUS PAR LA RESOLUTION DE L'EQUATION DE HELMHOLTZ POUR LE CAS NON ISOTHERME. L'ETUDE DE MODES PROPRES DE CAVITES EN PRESENCE D'ECOULEMENT MONTRE ENSUITE LES MODIFICATIONS NOTABLES QUE SUBISSENT LES CHAMPS ACOUSTIQUES SOUS L'EFFET D'UNE COUCHE LIMITE LAMINAIRE OU TURBULENTE. LES PROPRIETES DE LA PROPAGATION DES ONDES HARMONIQUES SONT ETUDIEES DANS UN GUIDE D'ONDE UNIFORME ET DANS UN RESONATEUR A DEUX COLS. POUR CE DERNIER, LES CHAMPS MOYENS DE VITESSE ET DE TEMPERATURE PRESENTENT UNE COMPLEXITE PROCHE DES CONDITIONS REELLES. LES SPECTRES DE REPONSE EN AMPLITUDE ET LES COURBES DE PERTE PAR TRANSMISSION PERMETTENT D'EXAMINER LES CARACTERISTIQUES DU RESONATEUR EN FONCTION DES DIFFERENTS PARAMETRES MIS EN JEU
Author: BRAC Jean
Publisher: Lavoisier
ISBN: 2746217244
Category :
Languages : en
Pages : 274
Get Book Here
Book Description
Cet ouvrage propose une méthode de construction de schémas numériques de grande précision sur la base d'une analyse spectrale de l'erreur. Ces schémas sont appliqués à la propagation des ondes mais ils peuvent l'être à la résolution par différences finies de tout autre système d'équations aux dérivées partielles. Plusieurs formulations du problème continu sont exposées mais on retient la formulation en vitesses de déplacement et en contraintes. D'autre part, une analyse des caractéristiques des équations de la propagation conduit à faire une comparaison avec les caractéristiques des équations de la mécanique des fluides et d'indiquer les conditions de la filiation. La discrétisation des équations est basée sur les schémas en grilles décalées. On effectue des développements de Taylor à des ordres élevés et une analyse de l'erreur de discrétisation par transformée de Fourier. Puis, on introduit la notion d'approximation optimale en contraste avec les approximations basées sur l'erreur de troncature. Les schémas construits restent de type convolutif. Le calcul s'avère très efficace sur la base de l'erreur relative de discrétisation. L'algorithme de calcul des coefficients d'approximation optimale est fourni en Fortran dans une annexe. L'élévation de l'ordre en temps consiste à reporter le calcul des dérivées d'ordre élevé en temps sur des dérivées d'ordre élevé en espace. Enfin, l'analyse des conditions de stabilité et de dispersion est réalisée pour prendre en compte l'approximation optimale des dérivées pour des ordres élevés en espace et en temps.
Author:
Publisher: Ed. Techniques Ingénieur
ISBN:
Category :
Languages : fr
Pages : 23
Get Book Here
Book Description
Author: Benjamin Cotté
Publisher:
ISBN:
Category :
Languages : fr
Pages : 230
Get Book Here
Book Description
La propagation acoustique en milieu extérieur fait intervenir des phénomènes physiques complexes, liés essentiellement aux variations de température et de vent dans la couche limite atmosphérique et aux frontières du domaine (effet de l’impédance du sol, de la topographie, ...). De plus, dans le contexte des trains `a grande vitesse (TGV), les sources de bruit sont étendues, en mouvement `a une vitesse relativement élevée, et de nature diverse (bruit de roulement, bruit d’origine aérodynamique, ...). Ce travail de thèse a pour but de modéliser la propagation du bruit des TGV en milieu atmosphérique, et de comprendre les phénomènes physiques associés `a ce type de propagation. La première partie de ce travail s’intéresse à l’effet de diffusion des fluctuations turbulentes de température et de vent sur les ondes acoustiques en présence d’une zone d’ombre acoustique. L’effet de diffusion d’une taille de structure turbulente donnée dépend de la fréquence acoustique et de la géométrie de propagation. Ce couplage entre échelles de turbulence, fréquence acoustique et géométrie est étudié `a l’aide d’outils de la théorie de la propagation des ondes en milieu aléatoire et de simulations d’ ́equation parabolique, afin d’estimer les plus petites et les plus grandes tailles de structure turbulente à prendre en compte dans une configuration donnée. Dans une deuxième partie, une méthode de résolution des équations d’Euler linéarisées par différences finies dans le domaine temporel est décrite. L’utilisation de schémas numériques optimisés permet d’appliquer ce modèle de propagation `a des configurations de propagation acoustique longue distance. Une des principales difficultés rencontrées avec les modèles temporels de propagation est la prise en compte de la réflexion des acoustiques sur un sol d’impédance finie. Des conditions limites d’impédance performantes d’un point de vue numérique sont proposées pour des modèles d’impédance couramment employés dans les études de propagation en milieu extérieur. Ces conditions limites sont obtenues en approchant l’impédance par des fonctions-type particulières, ce qui permet d’utiliser la méthode de convolution récursive. Elles sont validées dans des configurations de propagation bi- et tridimensionnelle, en considérant une atmosphère homogène puis une atmosphère stratifiée. Enfin, les applications spécifiques au bruit des TGV sont présentées dans une troisième partie. Dans un premier temps, un modèle de propagation du bruit des TGV basée sur une décomposition du train en un ensemble de sources ponctuelles équivalentes est d ́écrit. Les résultats de ce modèle sont comparés à des mesures réalisées à différentes distances de la voie de circulation en supposant les conditions de propagation homogènes. Le modèle de sources équivalentes est également couplé à un code d’ ́equation parabolique afin de prendre en compte l’effet d’un profil vertical de température ou de vent. Dans un deuxième temps, les phases d’approche et d’ ́eloignement de passages de TGV sont analysées afin de caractériser un phénomène de “grondement”. Dans certaines circonstances, il est en effet possible d’entendre un bruit similaire à un passage d’avion une dizaine de secondes avant ou après un passage de TGV, bruit qualifié de “grondement”. L’analyse s’appuie principalement sur des résultats expérimentaux, qui permettent de déterminer les caractéristiques du “grondement” et les circonstances dans lesquelles il se produit. Cette étude est complétée par des simulations d’ ́equation parabolique qui montrent l’importance du vent dans l’apparition de ce phénomène.
Author: Nathalie Dabin
Publisher:
ISBN:
Category :
Languages : fr
Pages : 176
Get Book Here
Book Description
