Jackie et le clan Kennedy en 160 photographies :: Angers-etc\nPersonnalité à la fabuleuse destinée, elle a incarné la grâce et l'élégance dans le monde au gré des déplacements de son époux, le président des États-Unis, John Fitzgerald Kennedy, assassiné le 22 novembre 1963 à Dallas. Jackie Kennedy-Onassis, née Jacqueline Bouvier, restera à jamais comme une icône devant l'éternité. Mais au-delà d'un destin autant romanesque qu'historique, elle deviendra un mythe absolu en restant la First Lady la plus élégante de l'histoire.\nDans la ville de l'impératrice Eugénie, l'exposition consacrée à Jackie Kennedy semble tout aussi impériale (Photo AngersEtc).\nSur le littoral basque, le temps d'un printemps-été enfin déconfiné et propice à la flânerie en plein air, la ville de Biarritz rend un bel hommage à Jackie Kennedy-Onassis par le biais d'une magnifique exposition intitulée Jackie et les Kennedy, qui se déroule sur plusieurs sites extérieurs (Grande plage, Golf du phare, square d'Ixelles, Grand jardin public, Côte des Basques). Cinq sites pour cinq thèmes (« Jacqueline Bouvier et John Fitzgerald Kennedy », « Dallas, mort d'un président », « Les rêves ne meurent jamais », « JFK, le meilleur golfeur de la Maison blanche », « 1007 jours à la Maison Blanche »).\nParce que l'ex-Première dame des États-Unis y aurait séjourné dans le passé, au même titre que son époux dans sa jeunesse, la cité balnéaire du Pays Basque s'est donc mise à l'heure de l'histoire américaine en organisant cette exposition d'envergure. À travers un cheminement à découvrir à travers la commune, la visite se savoure comme une balade,\nJusqu'au 5 septembre de cette saison estivale, ce sont 164 clichés au format conséquent qui attendent les vacanciers et potentiels visiteurs... Des textes - que nous restituons en regard de chaque photographie - accompagnent chacun de ces moments. Pour ceux qui n'auraient pas l'occasion de découvrir l'événement et souhaiteraient en savoir plus, l'organisateur met en vente les photographies de l'exposition, mais également le catalogue de l'exposition sur son site internet.\n" Ma vie à côté du président "\n(Bref aperçu en quelques exemples choisis)\nTout simplement, à propos de son mari de président, Jackie précise : « Sa vie est devenue la mienne ! »\nLe 15 septembre 1962, à bord du destroyer USS Joseph Patrick Kennedy Junior, baptisé ainsi en hommage à son frère défunt, John assiste à la première épreuve de l'America's Cup. Depuis leur mariage, Jackie ne cesse de se battre « contre ces fichues chaussettes de tennis » qu'il porte souvent avec ses costumes.\n« John lit de la manière la plus étrange qui soit. Il lit en marchant, en se rasant, à table, après le dîner. Même en nouant sa cravate ! Je le vois toujours lire, y compris dans son bain. Il aime bien les romans de Ian Flemming. En fait, il attrape le premier qui lui tombe sous la main. Tous les dimanches, il arrache les trois pages du supplément littéraire du Times et coche ce que je dois acheter. S'il pleut, il reste au lit pour lire encore et encore ».\nJackie ne rentrera pas dans l'ombre de John. Elle entreprend l'inventaire des 132 pièces de la Maison Blanche. Une partie du mobilier a disparu, dilapidée par d'anciennes Premières dames à court d'argent. Une autre est entreposée, en désordre, dans la cave de la résidence. L'un des rares portraits de George Washington trône dans des toilettes du rez-de-chaussée. Le décorateur français Stéphane Boudin est recruté de toute urgence.\n9 août 1963, Patrick Bouvier meurt deux jours après sa naissance. Il souffrait de détresse respiratoire. John se charge des obsèques. Jackie se sent coupable. Après une fausse couche, la disparition d'Arabella, ce fils était très attendu. Un mois après le drame, à Newport, ils célèbrent leurs dix ans de mariage. John lui présente le catalogue du célèbre antiquaire new-yorkais J.J. Kleijman : « Choisis ce que tu veux ! ». Jackie désigne un bracelet-serpent. De son côté, elle lui offre un album composé de photographies du jardin de la Maison Blanche et des roses qu'il aime tant.\nLe 10 juin 1963, John envoie un message historique au Kremlin : « Nous devons revoir notre attitude envers l'Union soviétique. Les Russes ont accompli de grandes choses dans l'espace et ont beaucoup souffert dans leur combat contre Hitler. (...) Car, en dernière analyse, notre lien commun fondamental, c'est que nous habitons tous sur cette planète. Nous respirons tous le même air. Nous chérissons tous l'avenir de nos enfants. Et nous sommes tous mortels. Touché, Krouchtchev exige que les quotidiens publient son discours dans son intégralité.\nLe 26 mai 1961, le prince souverain de Monaco Rainier III et la princesse Grace sont à la Maison blanche. Admirant le ravissant ensemble vert pomme porté par la princesse, John lui glisse à l'oreille : « Givenchy ? » Étonnée, elle lui sourit, il lui adresse un clin d'œil amical : « C'est moi qui paye les factures de Jackie ! »\nLe 21 mars 1962, le président du Pakistan, Ayub Kahn, est sous le charme de Jackie. Le lendemain, John lui fait savoir que son gouvernement est fort impressionné par son travail. « Ce que je souhaite le plus au monde, c'est être ta femme et la mère de tes enfants .»\nÀ des dizaines de milliers de kilomètres de Washington, en compagnie de sa sœur et de l'ambassadeur John Kenneth Galbraith, Jackie est en Inde avant d'être reçue au Pakistan. John l'a chargée d'une mission officieuse : consolider l'amitié des États-Unis avec les deux pays adversaires, sans froisser Nehru, ni Ayut Kahn.\n« Ces femmes sont littéralement folles de lui. Parfois, elles ne s'aperçoivent même pas de ma présence... Les hommes sont un savant mélange de divin et de malin. »\nExposition-événement : Jackie et les Kennedy. 64200 Biarritz. Jusqu'au 5 septembre 2021.\nPhotos : jackiebiarritz.com


Retour sur la fashionweek de Paris - Blog mode Caroline\nRetour sur la fashionweek de Paris\n26 décembre 2019 Émilie Fréchette Actualité 0\nL’univers de la mode a pris d’assaut les grandes villes de la planète au cours de la fashionweek 2019, et après Milan et New York, ce fut au tour de Paris d’apporter sa touche à cet évènement tant attendu dans la semaine du 1er Octobre, et les collections présentées n’ont pas manqué de faire tournemer les têtes.\nUne grande place accordée à l’écologie\nComme les thèmes de la préservation de la nature ont été largement mentionnés au cours des derniers mois, ceux-ci ont également occupé une large place dans cette fashionweek de Paris. Ainsi, la durabilité a été un élément récurrent dans les différentes tenues présentées, notamment par des créateurs tels que Marine Serre et Givenchy, qui ont proposé différentes créations utilisant des matériaux restaurés comme le denim. On a pu donc observer une large gamme de vêtements qui montrent qu’on peut être élégant sans pour autant puiser déraisonnablement dans les ressources de la planète.\nDes maisons comme Dior ont profité de l’occasion pour annoncer le fait qu’elles allaient procéder à des reboisements cette rentrée. Concevoir des pièces vestimentaires eco-friendly n’est pas une chose nouvelles pour les créateurs comme Stella McCarthney, et celle-ci a eu l’occasion de se mettre en avant autant que possible cette saison, et a d’ailleurs évoqué que sa collection de cette année sera parmi les plus respectueuses de l’environnement qu’elle n’ait jamais publiée.\nLes couleurs ont été mises à l’honneur\nComme ce fut le cas de bien des créateurs cette année, les pièces vestimentaires visibles sur les podiums à Paris au cours de la fashionweek ont été caractérisées par une grande richesse en couleurs. Cela se présente comme une rupture avec la tendance de l’année dernière, qui s’est plutôt orientée vers les costumes, autant pour femmes que pour hommes, donc des habits assez sobres. Balmain, reconnu pour ses créations en noir et blanc, n’a pas hésité à jouer avec les couleurs, tout en restant fidèle à ses palettes fétiches.\nCe fut également une occasion pour les vêtements du style années 1970 de refaire surface, dans la mesure où il n’y a rien de plus colorés que les vêtements de cette époque, et Léonard Paris n’a pas manqué de nous servir des pièces magnifiques dans cet esprit. En outre, les créateurs comme IsseyMiyake et JunyaWatanabe se sont tournées vers des tons de couleurs caractérisés par leur fraicheur.\nUn petit relent d’extravagance\nLa mode française est généralement caractérisée par le fait d’être sobre et élégante, mais beaucoup des créateurs ont décidé de casser cette tendance, et bien que l’élégance reste de mise, plusieurs des collections ont pris une tournure extravagante. Ainsi, on peut mentionner la collection un brin fantaisiste de ThornBrowne, « Marie Antoinette », qui s’inspire de la souveraine. De leur côté, les créateurs Vivienne Westwood et Manish Amora qui se sont tournés vers des tons de couleurs flashy, ainsi que des formes oversized et des motifs kaléidoscopiques un peu partout sur leurs créations présentées cette semaine.\nAstuces de beauté : quels avantages à réaliser un masque au visage ?\nPrendre soin de ses ongles avant une manucure : comment procéder ?\nQuels sont les avantages des savons au lait de chèvre et à l’huile d’olive ?\nJoindre le service client d’orange pour bénéficier d’une assistance de qualité\nUn peeling pour le visage est très agréable en hiver\nhdsex dans Choisir son maillot de bain selon sa morphologie, comment ?\nMan dans Choisir une robe de cocktail d’été, pour un mariage\nEden Fantasys dans Choisir une robe de cocktail d’été, pour un mariage\nClotilde Ranno dans Chemise blanche : Et si vous adoptez les looks de stars ?\nMode Eté 2015\nMode hiver 2015


Programme de rachat d’actions de bourse direct. Le fait que le contribuable ou l'un des membres de son foyer fiscal ait participé au financement de l'activité ne saurait à lui seul caractériser l'exercice à titre professionnel de l'activité. Ainsi, les frais de garde chez etoro pour avantages de Forex les actions à l’achat est de 0.22 dollar.\nAvantages de Forex, brokers régulé en Niger\n3 Les options sont des instruments financiers complexes présentant un risque de perte en capital. Les pertes peuvent être extrêmement rapides. Solde de la balance commerciale dans quelques pays de l'ue en 2016 (en milliards d'euros) - solde de la balance commerciale allemagne 255,1 pays-bas 60,0 italie 51,5 irlande 45,3 belgique 22,9 rép.tchèque 18,0 danemark 8,9 hongrie 7,5 pologne 4,7 slovénie 2,1 chypre -4,4 autriche -5,1 luxembourg -5,4 croatie -7,3 roumanie -10,0 portugal -10,8 grèce -18,7 espagne -20,4 france -64,8 royaume-uni -204,7. Cette activité ne demande pas obligatoirement un apport conséquent. Chaque investisseur peut investir en bourse en commençant par quelques centaines d’euros. Pour avantages de Forex acheter des actions, nul besoin d’une grosse somme d’argent. C’est pour cette raison que beaucoup se tournent vers cette option. Dans tous les cas, il faut toujours mieux se renseigner sur le terrain où vous voulez intégrer. La base doit être acquise avant de commencer une activité.\nActivités avantages de Forex du bureau du représentant résident.\nPlateforme d’échange de cryptomonnaies: le trading des actions boursières\nL’indicateur lui, est une pièce du puzzle. Il vient simplement donner une information pour un objectif bien précis: mesurer l’état d’une tendance, l’explosivité / volatilité d’un mouvement. En ce qui concerne l’outil financier, un pea dans une banque en ligne pour limiter les frais ou un compte-titre chez un courtier offrant des frais de courtage réduit. Ippi (avantages de Forex indice des prix des produits industriels) (mensuel) (déc).\nJust make sure you don’t use the rewards points as an excuse for unnecessary spending. Credit card interest almost always wipes out any rewards you earn, from a financial perspective, so it’s best to stick to what you would normally buy and pay off the card in full every month. En islande par exemple, il existe des fermes de «minage» de bitcoin ou autres cryptomonnaies, la température ambiante étant de – 5°, elle est idéale pour le matériel informatique. Ces fermes incluent des milliers de processeurs qui avantages de Forex fonctionnent en permanence. Je recommande vivement maxime. C'est quelqu'un de très compétent dans le domaine de la finance et très pédagogue pour transmettre son savoir.\nLe montant minimum d’achat/vente applicable aux crypto-monnaies équivaut à 25 eur/usd. Au vu de la avantages de Forex situation actuelle: une cassure des précédents plus hauts sur le btc (sur la zone des 20 000 donc) pourrait déclencher une nouvelle vague de fièvre acheteuse, et une nouvelle envolée spéculative (principalement motivée par le “fomo” comme disent les anglais, c’est à dire: la peur de rater le train). Dit voorbeeld toont een back-test van de dax index in een 5-minuten grafiek.\nForex; cfd sur matières premières, sur indices; cfd sur etf, sur obligations, sur actions; et cfd sur les crypto-monnaies. Volatilité: la volatilité est une mesure de l'ampleur des variatons du cours d'un actif financier. Elle sert de paramètre de quantification du risque avantages de Forex de rendement et de prix d'un actif financier. L’opérateur public, congo telecom qui a perdu une partie de son monopole st resté maître des services et infrastructures de télécommunication de base (stations terriennes, lignes fixes).


Video X salope Mature – Sexe Femme Mure » 2013 » juillet\nqu’est devenu l’actrice porno Amber Lynn?\nAmber Lynn , égérie du cinema x dans les années 80, nous accueille chez elle en bikini vintage qu’elle retire vite pour nous montrer ses nouveaux gros seins ainsi que ses petites fesses hautes! Puis l’ancienne porno star americaine caresse sa vulve bien mouillée, aux lèvres tombantes et peu poilue. Elle stimule son clito de […]\nune mature parisienne se fait baiser par deux jeunes\nune cougar francaise de 55 ans baize avec un jeune puceau qui voulait une longue baise avec sa femme de menage mature et ce 20 juillet il est enfin arrivé à ses fins.belle italienne mature se caresse les tetons en le regardant droit dans les yeux pour l’exciter puis cette cougare qui prends sont pied […]\nle gendarme baise avec l’inculpée mature\nVoici un extrait porno montrant un gendarme au beau milieu d’un interrogatoire qui vire rapidement au fist de chatte pour faire parler la présumée coupable de vol. Sans plus attendre, il retire son uniforme de policier pour etre plus à l’aise et lui broutte le vagin avec une agilité lingual incroyable, on sent l’homme qui […]\net si je te faisais une ejac bucale?\nBlack baise vieille, ejaculation bucale, mature à gros seins\ncette femme de 45 ans salope est une bourgeoise délurée qui en attendant de lecher plusieurs vielles chattes de 70 ans dans son prochain plan cul lesbien, va se rassasier avec la teub d’un antillais qui va savoir profiter de ses grosses mamelles matures comme il se doit. En effet, le black bien membré ne […]\nPatricia se fait baiser dans une voiture par 2 amis\nOn dit que les matures de 60 ans sont des exemple à prendre, mais Patricia est l’exception qui infirme la regle. Dans une voiture qui roule à pleine vitesse sur la nationale, la salope du nord pas de calais ne se refuse aucun plaisir et quand elle a envie de baiser, c’est tout de suite […]


-12. Dispersion atmosphérique (Mécanismes et outils de calcul) - PDF Free Download\nDownload "-12. Dispersion atmosphérique (Mécanismes et outils de calcul)"\n1 Méthodes pour l évaluation et la prévention des risques accidentels (DRA-006) -1 Dispersion atmosphérique (Mécanismes et outils de calcul) Ministère de l Ecologie et du Développement Durable Direction des Risques Accidentels D é c e m b r e 0 0\n2 -1 Dispersion atmosphérique (Mécanismes et outils de calcul) INERIS Ce document comporte 61 pages DIRECTION DES RISQUES ACCIDENTELS D é c e m b r e 0 0\n3 PAGE DE VALIDATION DISPERSION ATMOSPHERIQUE Rédaction initiale Auteurs Qualité Date Emargement J.C Couillet Ingénieur Unité PHEN 00 Signé Dans le cadre de la procédure générale qualité de l INERIS et en respect du paragraphe 14. du manuel qualité, ce document a fait l objet de relectures et d un contrôle par des vérificateurs Vérificateur final Qualité Date Emargement O. Salvi Délégué scientifique 00 Signé Approbateur Qualité Date Emargement B. Faucher Directeur DRA 00 Signé Autres personnes ayant participé à la relecture du document : Y. MOUILLEAU S. DUPLANTIER J.M. LACOME E. BERNUCHON 1/61\n4 PREAMBULE Le présent document a été établi : au vu des données scientifiques et techniques disponibles ayant fait l objet d une publication reconnue ou d un consensus entre experts, au vu du cadre légal, réglementaire ou normatif applicable. Il s agit de données et informations en vigueur à la date de l édition du document, le 30 Septembre 00. Le présent document comprend des propositions ou recommandations. Il n a en aucun cas pour objectif de se substituer au pouvoir de décision du ou des gestionnaire(s) du risque ou d être partie prenante. /61\n5 TABLE DES MATIERES 1 OBJECTIF ET DOMAINE D APPLICATION OBJECTIF ET CONTEXTE DOMAINE D APPLICATION Rappels succincts sur le phénomène Démarche observée... 5 MÉCANISMES PHYSIQUES CONDITIONS DE REJET Etat physique initial du polluant Réactivité avec l atmosphère (stabilité chimique) Cas particulier des particules Mode d émission Déplacement du nuage Dilution du nuage - Entraînement d air Echanges thermiques CONDITIONS MÉTÉOROLOGIQUES Structure de l atmosphère Echelle des mouvements atmosphériques Structure de la turbulence atmosphérique à petite échelle Stabilité atmosphérique et conditions météorologiques ENVIRONNEMENT Effets mécaniques Effets thermiques MODÉLISATION CONTEXTE Essais en grandeur réelle Essais à échelle réduite Outils de calcul Détermination des effets sur la santé humaine MODÈLES NUMÉRIQUES Modèles gaussiens Modèles intégraux Modèles tri-dimensionnels Avantages et inconvénients Autres types de modèles : Modèles simplifiés RÉFÉRENCES /61\n6 1 OBJECTIF ET DOMAINE D APPLICATION 1.1 OBJECTIF ET CONTEXTE Depuis l année 000, le Ministère en charge de l Environnement (anciennement Ministère de l Aménagement du Territoire et de l Environnement devenu Ministère de l Ecologie et du Développement Durable) finance un programme d études et de recherches, intitulé «Recueil des méthodes utilisées à l INERIS dans le domaine des risques accidentels» (DRA-006). L objet de ce programme est de réaliser un recueil global formalisant l expertise de l INERIS dans le domaine des risques accidentels. Ce recueil sera constitué de différents rapports consacrés aux thèmes suivants : les phénomènes physiques impliqués en situation accidentelle (incendie, explosion, BLEVE ) l analyse et la maîtrise des risques, les aspects méthodologiques pour la réalisation de prestations réglementaires (étude de dangers, analyse critique..) Chacun de ces documents reçoit un identifiant propre du type «-X» afin de faciliter le suivi des différentes versions éventuelles du document. In fine, ces documents décrivant les méthodes pour l évaluation et la prévention des risques accidentels, constitueront un recueil des méthodes de travail de l INERIS dans le domaine des risques accidentels. 1. DOMAINE D APPLICATION Le présent rapport, baptisé "-1" présente la démarche adoptée par l INERIS pour l étude de la dispersion atmosphérique, un des sujets retenus dans le thème «phénomènes physiques» cité ci-dessus. Il s inscrit dans une démarche de dissémination et de valorisation du savoir-faire de l INERIS auprès des pouvoirs publics, des industriels et du public. 4/61\n7 1..1 RAPPELS SUCCINCTS SUR LE PHENOMENE La dispersion atmosphérique caractérise le devenir dans le temps et dans l espace d un ensemble de particules (aérosols, gaz, poussières) rejetées dans l atmosphère. L émission d un produit à l atmosphère peut revêtir un caractère : soit chronique, avec des émissions à l atmosphère plus ou moins continues ou périodiques dans le temps. Les rejets sont les sous-produits indésirables de toute activité humaine, par exemple, les gaz d'échappement des voitures, les fumées d'usines, de chauffages urbains ou toute forme de pollution diffuse (émission de gaz issu d'une décharge...). soit accidentel, avec des émissions à l atmosphère ponctuelles dans le temps, non désirées comme la fuite d'une cuve ou un dégagement de fumées dû à un incendie. Ce rapport traite des mécanismes de dispersion atmosphérique en situation accidentelle. Dans ce cadre, la connaissance des mécanismes de dispersion atmosphérique peut s appliquer dans plusieurs contextes : lors d une étude des dangers, la réalisation d une étude de dispersion atmosphérique permet d'envisager a priori, à titre prévisionnel, les risques potentiels d une installation industrielle ; lors d'un rejet accidentel, la modélisation de la dispersion atmosphérique permet d'évaluer les mesures à prendre en temps réel ; en situation post-accidentelle, l analyse des conditions de dispersion dans l atmosphère peut permettre de mieux comprendre le déroulement et les conséquences de ces accidents. Il s agit ici à la fois d identifier les mécanismes physiques intervenant dans la dispersion atmosphérique d une substance, et de présenter les principales méthodes d estimation de la dispersion atmosphérique d un produit et en particulier celles relatives aux outils de résolution numérique. En amont à la modélisation de la dispersion atmosphérique d une substance, il convient généralement de quantifier un «terme source», c est-à-dire de caractériser le rejet de la substance étudiée vers l air en termes de débit, température, état physique... L étude des mécanismes physiques relatifs aux termes sources et l étude des méthodes d évaluation associées n est pas traitée dans ce document et doit faire l objet d un document spécifique. Ce rapport traite de la phase aval, c est-à-dire des mécanismes de dispersion atmosphérique et des outils de modélisation associés. 1.. DEMARCHE OBSERVEE La démarche adoptée pour le présent rapport consacré à l étude du phénomène de la «dispersion atmosphérique» est décrite ci-après : chapitre : Mécanismes physiques ; chapitre 3 : Modélisation. 5/61\n8 MECANISMES PHYSIQUES Les conditions de dispersion atmosphérique d un produit vont dépendre de plusieurs paramètres : les conditions de rejet (nature du nuage de produit, mode d émission...) ; les conditions météorologiques (champ de vent, de température...) ; l environnement (nature du sol, présence d obstacles, topographie...). Différents processus dont l importance relative dépend des conditions de rejet, entrent en jeu et agissent simultanément ou successivement. Ces différents mécanismes physiques sont présentés dans cette partie..1 CONDITIONS DE REJET.1.1 ETAT PHYSIQUE INITIAL DU POLLUANT Avant d être rejetées dans l atmosphère, les substances sont stockées dans des citernes, des sphères, des bouteilles, des containers, des fûts... Les substances peuvent se trouver sous la forme : d un gaz (sous pression ou non), d un liquide (réfrigéré ou non), d un gaz liquéfié. Pour les deux derniers cas, les rejets vont donner lieu à une émission diphasique pouvant conduire à la formation d une nappe. Pour diverses raisons, ces produits peuvent être relâchés dans l environnement. La rupture du confinement va entraîner des phénomènes variables selon le produit considéré, le type de stockage, l importance et la position du point de rejet..1. REACTIVITE AVEC L ATMOSPHERE (STABILITE CHIMIQUE) Dans la plupart des cas, les produits rejetés sont des produits inertes ou stables (le monoxyde de carbone, par exemple), c est-à-dire ne subissant pas ou peu de transformations ou d altération d origine chimique en contact avec l atmosphère. Lorsque des produits réactifs primaires (typiquement le dioxyde de soufre SO, les oxydes d azote NO x ) sont émis, les transformations chimiques sont généralement négligées pour des périodes d étude de l ordre de quelques minutes, voire quelques heures, comme cela peut être le cas en risque accidentel 1. Ces produits sont alors traités comme des produits inertes ou stables. 1 Dans le domaine du risque chronique, les mécanismes de dispersion sont étudiés sur des périodes de temps plus grandes, il devient alors nécessaire de prendre en compte ces transformations chimiques. De la même façon, pour les polluants secondaires issus de la pollution primaire comme l ozone et le benzène, il convient de faire appel à des modèles munis de modules chimiques. 6/61\n9 .1.3 CAS PARTICULIER DES PARTICULES Les produits émis peuvent aussi se trouver sous forme de particules (poussières...). Dans ce cas, il convient d emblée de tenir compte de l influence des effets de la pesanteur pour analyser les possibilités de dépôt éventuel. A cet égard, il est généralement admis que lorsque les particules ont un rayon moyen inférieur à 5 m, elles se comportent comme un gaz. En revanche, les particules ayant un rayon supérieur à 5 m sont soumises à la force gravitationnelle non négligeable [Hanna et. al., 198b]. La force gravitationnelle agit sur ces particules en même temps que les autres forces de transport et de dispersion et lorsqu'elle devient prépondérante, la particule peut se déposer sur le sol..1.4 MODE D EMISSION Parmi les types de rejet à l'atmosphère, on distingue les rejets instantanés des rejets continus. a- Rejet instantané La ruine consiste en la rupture totale du confinement entraînant un relâchement immédiat et instantané de la totalité du contenu. Les rejets de type instantané sont principalement caractérisés par la quantité émise à l atmosphère et la pression de la capacité au moment de la rupture (éclatement). Cette dernière peut être différente de la pression de stockage dans certains cas. Lors d'un rejet instantané de gaz, en l'absence d'obstacle de taille significative à proximité, l'extension initiale du nuage est relativement isotrope. Il en résulte un volume gazeux sphérique, ou semi-sphérique si l'émission a lieu au niveau du sol (voir la Figure 1). Figure 1 : Rejet instantané d un produit Dont les propriétés physiques sont identiques dans toutes les directions de l espace. 7/61\n10 b- Rejet continu La fuite continue se produit à partir d une conduite ou d une cuve de stockage ou de l émission par une cheminée d extraction. La fuite va principalement dépendre des conditions de stockage au moment de la fuite mais aussi de la taille, de la géométrie et de la hauteur de l orifice de fuite. A partir de ces données, il est possible d estimer un terme-source (débit, vitesse, température... du rejet). Pour un rejet à partir d un réservoir, le débit varie généralement dans le temps puisque durant la vidange du réservoir, la quantité de produit stocké diminuant, la pression de fuite s'atténue dans le temps. Le panache est dans la plupart des cas de forme allongée (voir la Figure ). c- Détermination du type de rejet Figure : Rejet continu d un produit Plusieurs critères permettent de différencier les rejets instantanés des rejets continus. Hanna [Hanna, 198a] propose de comparer la durée d'émission du rejet à la durée de transfert. La durée de transfert est définie comme la durée moyenne nécessaire au produit pour migrer de la source d'émission vers le point considéré. Cette durée de transfert, notée t, peut être estimée par le rapport : avec : x t u x : distance entre la source et le point considéré ; u : vitesse du nuage. Si la durée du rejet est quatre fois supérieure au temps de transfert au point d observation, le rejet est à admettre comme continu [UIC, 1995]. D autres modèles retiennent comme critère le ratio entre la largeur et la longueur du nuage. Si ce ratio est plus grand qu une valeur limite, le rejet est qualifié d instantané [Witlox, 000]. 8/61\n11 .1.5 DEPLACEMENT DU NUAGE Plusieurs mécanismes participent au déplacement et à la dilution du nuage avec une intensité différente et variable selon l évolution des caractéristiques du nuage. pour des rejets de grande vitesse, la quantité de mouvement va être prédominante dans les premiers instants, puis va être réduite avec les effets d entraînement de l air et du frottement du sol. Tout rejet est émis avec une certaine quantité de mouvement dépendant de ses conditions de stockage. Ainsi, plus la pression de stockage est grande, plus la quantité de mouvement du nuage est initialement importante. Pour un rejet continu non dirigé directement vers le haut ou vers le bas, cette quantité de mouvement va transporter le nuage de gaz vers l avant. Pour un rejet instantané, cette quantité de mouvement est supposée être distribuée tout autour du centre du nuage. Elle participe donc à son expansion, mais ne contribue pas de ce fait au déplacement du nuage. Quand le nuage est en contact avec le sol ou atteint une couche d inversion de température (voir le chapitre..4.5 Conditions d inversion de température), la composante verticale de la quantité de mouvement est transformée en composante horizontale, ce qui se traduit par une augmentation de la largeur du nuage. si le nuage est plus dense que l air, l effet des forces de gravité est l effet le plus fort ; Si le nuage est plus lourd que l air, il va s effondrer sur le sol ce qui va générer de la turbulence, entraîner de l air et changer sa section transversale circulaire en une ellipse tronquée. Le sol exerce alors une force de frottement sur les zones du nuage qui sont en contact, ce qui diminue la quantité de mouvement du nuage (voir le chapitre Densité du gaz). si le nuage au sol a une densité suffisamment petite par rapport à celle de l air alors le nuage peut décoller du sol ; enfin, la quantité de mouvement va être perdue progressivement du fait de l entraînement de l air et de l éventuel frottement du nuage avec le sol, et à mesure que la densité devient proche de celle de l air, l effet lié à la turbulence atmosphérique devient le plus important, entraînant en permanence de l air dans le nuage. Le nuage est entraîné par le vent. 9/61\n12 .1.6 DILUTION DU NUAGE - ENTRAINEMENT D AIR Plusieurs sources de turbulence vont entraîner de l air dans le nuage et le diluer QUANTITE DE MOUVEMENT INITIALE JET TURBULENT En fonction des conditions de stockage, le nuage peut avoir une quantité de mouvement importante au début du rejet. La différence entre la quantité de mouvement du nuage et celle de l air environnant provoque une turbulence importante, à l origine d un taux important d entraînement d air favorisant la dilution du rejet. On parle alors de jet turbulent. Les caractéristiques du jet dépendent de la quantité de mouvement du gaz émis. L absence de prise en compte de ce phénomène peut conduire à une surestimation des valeurs de concentration dans l atmosphère CISAILLEMENT DU CHAMP DE VENT Si le nuage monte ou descend, alors il cisaille le champ de vent horizontal. La différence entre la quantité de mouvement du nuage et de l air produit des tourbillons favorisant le mélange de l air dans le nuage DENSITE DU GAZ Gaz dense Lorsque le gaz est plus dense que l air, il est qualifié de «gaz lourd». Ce type de rejet est rencontré lors des fuites : d un produit rejeté à température ambiante de masse volumique plus importante que celle de l'air (propane, chlore, phosgène...) ; d un produit plus léger que l'air à température ambiante, mais rejeté à une température assez faible pour que sa masse volumique soit supérieure à celle de l'air. Lors de sa dispersion, le nuage de gaz voit sa masse volumique diminuer du fait de sa dilution avec l'air ambiant et de son réchauffement ; d un aérosol. La phase condensée provoque par sa seule présence une augmentation considérable de la masse volumique globale du mélange. En outre, le changement de phase intervenant lors de la dispersion avec l évaporation des gouttelettes de liquide maintient le nuage à une température voisine de la température d'ébullition du produit. La densité du rejet reste de ce fait supérieure à celle de l air. A titre d illustration, c est le cas de rejets de gaz liquéfiés tels l ammoniac, le chlore... de produits qui en réaction avec l air vont engendrer par réactions chimiques de nouveaux produits qui augmentent plus ou moins localement la densité du nuage (exemple du monoxyde d'azote réagissant avec l'oxygène de l'air pour donner du dioxyde d'azote plus dense que l'air). Lorsque la densité du gaz rejeté est plus grande que celle de l atmosphère, les forces de gravité influencent de façon importante la dispersion du nuage. A proximité du rejet, la force de gravité est à l origine d un mouvement d ensemble vers le sol et l expansion verticale du nuage est réduite. En contact avec le sol, le nuage se répand, poursuit son effondrement ce qui favorise son développement transversal. 10/61\n13 La dilution du nuage s opère sur la face supérieure du nuage et sur les faces latérales. Sur la face supérieure, les échanges sont assurés par les mouvements de turbulence de l atmosphère et grâce à la turbulence induite par la différence de vitesse entre l atmosphère et le nuage. Enfin, il s établit sur les faces latérales du nuage une recirculation interne favorisée par l étalement du nuage au sol à mesure de son effondrement. Ce comportement est illustré par la Figure 3 ci-après avec la représentation de coupes d isoconcentrations à 1%, 6%, 11% et 16% relatives à la modélisation de l essai n 18 de Thorney Island avec le modèle Mercure - Gaz lourds [Riou, 1989]. 11/61\n14 Figure 3 : Dispersion d une bouffée de gaz dense [Riou, 1989] 1/61\n15 Gaz léger Lorsque le gaz est moins dense que l air, il est qualifié de «gaz léger». Ce type de rejet concerne les fuites : d un produit rejeté à température ambiante de masse volumique moins importante que celle de l'air (hydrogène...) ; d un produit moins léger que l'air à température ambiante, mais rejeté à une température assez importante pour que sa masse volumique soit plus petite que celle de l'air (exemple des fumées d incendie). Lors de sa dispersion, le nuage de gaz voit sa masse volumique augmenter du fait de sa dilution avec l'air ambiant et de son refroidissement. Si la densité du gaz est suffisamment faible, l effet des forces d Archimède 3 est susceptible de favoriser l ascension du nuage de gaz. Si la vitesse ascensionnelle est suffisante, le frottement des bords du nuage avec l atmosphère lors de l ascension peut favoriser la dilution (apparition de turbulence sur les bords du nuage). Enfin, le cœur du nuage plus léger que la périphérie peut s élever plus rapidement avec pour conséquence un accroissement de la surface de contact entre l atmosphère et le gaz. Tous ces phénomènes sont susceptibles de favoriser la dilution du nuage [Proust, 1999]. Gaz neutre (ou passif) Le gaz est qualifié de «gaz neutre», lorsque le gaz : n a pas de quantité de mouvement, a la même densité que l air, a la même température que l air. Ce type de rejet est rencontré lors des fuites d un gaz de même densité que l air ou très dilué quelle que soit sa densité à l état pur. Le gaz neutre n'apporte aucune perturbation mécanique à l'écoulement atmosphérique et va se disperser du fait de la seule action du fluide porteur, l air. Le déplacement et la dilution du gaz vont alors dépendre du vent et de la turbulence atmosphérique d origine mécanique ou thermique..1.7 ECHANGES THERMIQUES Outre les effets mécaniques, les effets thermiques jouent un rôle significatif sur la dilution. Pour des gaz dont la température est très différente de celle de l air et du sol, les échanges thermiques vont accélérer la dilution par convection thermique. 3 Les forces d Archimède sont les forces de pression hydrostatiques qui s exercent sur un corps placé dans un fluide et soumis au champ de pesanteur terrestre. 13/61\n16 . CONDITIONS METEOROLOGIQUES Les phénomènes météorologiques à prendre en compte dépendent de la dimension du domaine étudié. Dans la majorité des cas traités dans le domaine du risque industriel, l échelle spatiale ne dépasse pas quelques kilomètres (au plus une dizaine) et celle de temps quelques minutes, à quelques heures. Le suivi du nuage de polluant se fera alors dans la couche limite atmosphérique...1 STRUCTURE DE L ATMOSPHERE Les propriétés de l atmosphère vont jouer un rôle plus ou moins moteur dans la dispersion d un produit. Dans le cadre du risque accidentel, les phénomènes liés à la dispersion d un produit ne se produisent que dans la couche inférieure de l atmosphère, dite couche limite atmosphérique, de quelques centaines de mètres d épaisseur (au plus 1 à km) dans la plupart des cas. Dans ce domaine, il importe donc de connaître, en particulier, la structure de l écoulement du vent et du champ des températures qui vont être caractérisés par deux paramètres fondamentaux : le gradient vertical de vitesse ; le gradient vertical de température. Ces gradients vont particulièrement être influencés principalement par deux mécanismes physiques : d une part, les effets de frottement de l air sur le sol. Ces effets dépendent des conditions orographiques 4 locales ; d autre part, les échanges de chaleur entre le sol et l atmosphère. Ces échanges vont varier avec le cycle diurne, les conditions météorologiques et la nature du sol. 4 Relatif à l agencement des reliefs terrestres. 14/61\n17 Selon la localisation géographique, ces mécanismes vont être plus ou moins actifs et vont jouer sur l épaisseur de la couche limite atmosphérique. Celle-ci va dépendre de la rugosité du sol mais aussi de la météorologie (ensoleillement, nébulosité...). Ainsi, la couche limite atmosphérique peut être décomposée en deux parties avec (voir la Figure 4) [Turbelin G., 000] : Figure 4 : Décomposition de la couche limite atmosphérique Une couche de surface dont l épaisseur est de l'ordre de quelques dixièmes de l'épaisseur de la couche limite atmosphérique (soit entre une centaine de mètres les jours de vent fort et quelques mètres les nuits claires et par vent faible). La force de cisaillement et les flux de chaleur y sont supposés quasiment constants en fonction de l altitude. Le principal moteur du flux radiatif est le rayonnement solaire. La distribution verticale de température (stratification thermique de l air) induit aussi des mouvements verticaux des masses d air chauffées et refroidies à proximité du sol. Le vent résulte de l équilibre entre les forces de pression 5, et les forces de frottement (les effets de la force de Coriolis sont, dans cette zone, négligeables) et de la stratification thermique de l air. 5 Les forces de pression proviennent des inhomogénéités de pression dues à l échauffement inégal de l air. Elles sont perpendiculaires aux isobares et dirigées des hautes pressions vers les basses pressions. 15/61\n18 Le profil vertical de la vitesse du vent est perturbé du fait du frottement de l air sur le sol. La présence d obstacles ou de discontinuités va modifier la trajectoire du vent. La couche de surface est donc une zone à fort gradient de vitesse de vent. Au sein de la couche de surface, il est possible de distinguer dans la partie inférieure, la sous-couche rugueuse. L épaisseur de cette zone varie de moins d un mètre (en mer) à quelques dizaines de mètres (dans les zones fortement urbanisées). Du fait de la présence des éléments de rugosité dans cette zone, l écoulement de l air est fortement turbulent, non homogène et instationnaire. Cette zone est souvent caractérisée par une rugosité globale homogène à une longueur reliée à la nature du recouvrement. Une couche dite de transition surmontant la précédente et s'étendant jusqu'à la limite supérieure de la couche limite atmosphérique. La force de cisaillement et flux de chaleur y décroissent avec l altitude et deviennent négligeables au sommet de la couche. Le champ de vent résulte de l équilibre entre les forces de pression 6, les forces de frottement et la force de Coriolis 7 ainsi que de la stratification thermique de l air. Cette dernière provoque une rotation de la direction du vent à travers la couche (spirale d Ekman) vers la direction du vent géostrophique propre à l atmosphère libre. Plus haut en altitude, dans l atmosphère libre, les forces de frottement n ont plus d effet, le vent résulte uniquement de l équilibre entre les forces de pression et la force de Coriolis. Le vent dit géostrophique, est uniforme, horizontal et très peu turbulent. Au niveau des réseaux d'observation météorologique, ce sont les caractéristiques du vent synoptique qui sont disponibles. Le vent synoptique est observé sur une zone considérée comme non influencée par la topographie ou encore l'occupation du sol. Il caractérise le flux atmosphérique de grande échelle. La vitesse du vent synoptique est mesurée à 10 m audessus de la surface du sol... ECHELLE DES MOUVEMENTS ATMOSPHERIQUES La vitesse moyenne du vent joue deux rôles dans le processus de diffusion : l un dans le transport du produit, l autre dans sa diffusion. La vitesse instantanée du vent reflète la nature turbulente du vent. Elle se décompose en une partie moyenne et une partie fluctuante (ou turbulente) aléatoire : la vitesse moyenne est une quantité qui traduit la «force» du vent en un point de l espace sur des périodes de temps qui varient entre 10 minutes et une heure pour que le résultat ne soit sensible ni aux fluctuations lentes, ni aux variations rapides du vent ; 6 Les forces de pression proviennent des inhomogénéités de pression dues à l échauffement inégal de l air. Elles sont perpendiculaires aux isobares et dirigées des hautes pressions vers les basses pressions. 7 La force de Coriolis est perpendiculaire à la vitesse du vent, orientée vers sa droite dans l hémisphère Nord, du fait de la rotation de la terre. 16/61\n19 la partie fluctuante du vent est une fonction aléatoire qui représente les fluctuations turbulentes de la vitesse autour de sa moyenne. Ces fluctuations sont décrites statistiquement (écarts types, densités spectrales, par exemple). Le vent est un écoulement turbulent d air, constitué d une multitude de tourbillons de tailles diverses, imbriqués les uns dans les autres, les petits étant transportés par les plus gros, euxmêmes transportés par le mouvement d ensemble. L écoulement de l air atmosphérique est donc constitué d une large gamme d écoulements interdépendants caractérisés par des tailles allant du millimètre au millier de kilomètres [Crabol et al., 1999]. Chaque phénomène atmosphérique est identifié par son extension horizontale et par sa durée de vie (ces deux paramètres sont fortement corrélés). De ce fait, la vitesse du vent en un point donné de l espace, présente de fortes variations plus ou moins irrégulières en termes d amplitudes et de fréquences différentes. Les mouvements à grandes échelles (échelles synoptiques) ont une taille supérieure à la centaine de kilomètres. Ces mouvements ont des périodes variant : de 1 an, du fait des variations saisonnières ; à 4 jours avec les variations associées aux perturbations qui traversent une région donnée ; et 4 heures, étant données les variations journalières, et les phénomènes thermiques jour-nuit. Les mouvements à petites échelles ont une taille inférieure au kilomètre et ont une durée de vie de quelques minutes au maximum (micro-échelles). Ils sont liés à la turbulence et sont générés, dans la couche limite atmosphérique, par la présence d obstacles ou par la rugosité des sols. A cette échelle, l écoulement d air turbulent est constitué d une multitude de tourbillons de tailles différentes, emportés par le mouvement d ensemble 8. _ Des mouvements de tailles intermédiaires (meso-échelles) assurent la transition entre les précédents. Entre les petites et les meso-échelles, on distingue parfois des échelles dites sub-meso. 8 [Hunt (199)] précise qu il existe des «vortical eddies» qui sont des zones à fort rotationnel transportées de façon aléatoire dans l écoulement, et des «structural eddies», qui sont des tourbillons plus localisés liés aux propriétés particulières de l écoulement. En anglais, le vocable «eddies» désigne les structures tourbillonnaires liées aux écoulements turbulents, tandis qu un autre vocable, «vortex», désigne de façon plus générale tous les mouvements de rotation de fluides autour d un axe central. Cette distinction n a pas d équivalent en français ou, dans les deux cas, le vocable «tourbillon» est utilisé. 17/61\n20 Le Tableau 1 donne une correspondance approximative entre la taille moyenne des mouvements atmosphériques, les périodes qui leur sont associées. Taille [km] Période 1 s 10 s 1 min 10 min h 15 h 7 j 10 j Echelles Micro-échelles Méso-échelles Grandes échelles Tableau 1 : Echelles des mouvements atmosphériques [Atkinson, 1995] Parallèlement, une analyse spectrale permet de retrouver l échelle des fluctuations de l écoulement de l air (voir la Figure 5) [Crabol,1999]. Figure 5 : Spectre énergétique de la vitesse du vent Les pics du spectre montrent dans quelle mesure les tourbillons participent effectivement à l énergie turbulente. On retrouve les tendances explicitées précédemment avec : le pic relatif à une période d environ 100 heures traduisant les variations de vitesse du vent associées aux passages successifs de cyclones et d anticyclones ; le pic relatif à une période d environ 1 heures traduisant l augmentation de la vitesse du vent dans la journée et sa diminution durant la nuit ; le pic relatif à une période comprise entre 10 s et 10 min traduisant la turbulence de petite échelle. Pour l étude de la dispersion de produits, dans la couche limite atmosphérique, seules les informations concernant la structure des mouvements de petites échelles (micro-échelles) sont pertinentes. A cette échelle, la forme du spectre dépend fortement de la stratification thermique de l atmosphère. 18/61\n21 ..3 STRUCTURE DE LA TURBULENCE ATMOSPHERIQUE A PETITE ECHELLE La turbulence atmosphérique à petite échelle a deux origines distinctes : la première d origine «thermique», due à la distribution de température dans l atmosphère ; la seconde d origine «mécanique» générée par le cisaillement du vent en contact avec le sol. Les conditions météorologiques ne vont jouer que sur l aspect thermique de la turbulence en influençant le gradient vertical de température TURBULENCE D ORIGINE THERMIQUE Le mécanisme sous-jacent à la turbulence d origine thermique est le suivant. Un volume d air élémentaire déplacé adiabatiquement 9 vers le haut se détend du fait de la baisse de pression et par conséquent se refroidit. L air sec se refroidit de 0,98 C pour une élévation de 100 m (gradient de l adiabatique sèche). L air saturé en eau se refroidit de 0,55 C pour une élévation de 100 m (gradient de l adiabatique saturée). L évolution de la température de l atmosphère avec l altitude peut s écarter de ces conditions (comme le montre la Figure 6) : si le gradient de température de l atmosphère est inférieur au gradient thermique de l adiabatique alors tout volume d air déplacé vers le haut a, avant équilibre thermique, une température plus élevée que l air qui l entoure. La masse volumique du volume élémentaire est plus petite que l air qui l environne et continue de ce fait de se déplacer vers le haut. L atmosphère favorise alors les déplacements de gaz et peut être qualifiée d instable ; si le gradient de température de l atmosphère est supérieur au gradient thermique de l adiabatique alors tout volume d air déplacé vers le haut a, avant équilibre thermique, une température plus petite que l air qui l entoure. La masse volumique du volume élémentaire est plus importante que l air qui l entoure et tend à se déplacer vers la bas à sa position initiale. L atmosphère ne favorise alors pas les déplacements de gaz et peut être qualifiée de stable. 9 Sans échanges de chaleur avec l extérieur. 19/61\n22 Altitude [m] Adiabatique saturée Adiabatique séche Atmosphère instable 100 Atmosphère stable Gradient thermique [ C] Figure 6 : Comparaison entre le gradient de température de l atmosphère et le gradient thermique de l adiabatique Des phénomènes d inversion de température peuvent aussi survenir se traduisant par une augmentation de la température avec l altitude puis une baisse...3. TURBULENCE D ORIGINE MECANIQUE Lorsque le nuage est en contact avec le sol, les effets de rugosité dus à la nature de la surface créent une zone de turbulence à l interface d origine mécanique. La notion de rugosité sous-entend que les éléments (végétation...) qui la définissent, sont beaucoup plus petits que la taille du nuage. Si ces éléments sont de taille comparable à celle du nuage, ils doivent être considérés comme des obstacles qui génèrent des perturbations beaucoup plus importantes que la turbulence de petite échelle (tourbillons, zones de recirculations, brise). Cet aspect est étudié au chapitre.3 Environnement DEGRE DE STABILITE DE L ATMOSPHERE Nombre de Richardson [Richardson, 190] Comme indiqué précédemment, selon l importance d un type de turbulence par rapport à l autre, l atmosphère peut être stable, neutre ou instable. Le degré de stabilité atmosphérique peut être caractérisé en comparant le gradient vertical de température gradient de température adiabatique. z au 0/61\n23 Cette comparaison s effectue généralement au moyen du nombre de Richardson de gradient, qui intègre à la fois le paramètre lié au gradient de température et celui lié à la vitesse du vent : R i g z U z avec : R i 0 traduisant une atmosphère instable ; R i 0, une atmosphère neutre ; R i 0, une atmosphère stable. Longueur de Monin-Obukhov [Monin-Obukhov, 1954] La longueur de Monin-Obukhov permet aussi de qualifier l importance relative de la turbulence mécanique et de la turbulence convective, à une altitude donnée. Elle est donnée par : U L avec : 3 * C g k q p T : masse volumique de l air [kg/m 3 ] ; U * : vitesse de frottement [m/s]. Cette vitesse est égale à cisaillement de vent au sol ; C p : chaleur massique de l air à pression constante [J/kg.K] ; T : température de l air [K] ; g : constante de gravité [9,81 m/s ] ; k : constante de Von Karman [0,4] ; q : flux moyen de chaleur dans la direction verticale [W/m]. Ainsi, lorsque : avec, force de q=0, L, le gradient vertical de température est adiabatique ; q<0, L>0, le gradient vertical de température est positif et l atmosphère stable ; q>0, L<0, le gradient vertical de température est négatif et l atmosphère instable. La signification géométrique de cette longueur est la suivante : à des altitudes très petites devant L, la turbulence mécanique prédomine ; à des altitudes de l ordre de L ou plus grandes, la turbulence convective prédomine. Le rapport L z indique l importance relative de la turbulence mécanique et de la turbulence convective à l altitude z. 1/61\n24 Classification de Pasquill Une classification de la stabilité atmosphérique peut être effectuée au travers des classes de Pasquill [Pasquill, 1974] qui varient de A à F, de l atmosphère la plus instable à la plus stable. Les classes de stabilité peuvent être associées à des conditions météorologiques données (voir le Tableau ). Vitesse du vent à 10 m JOUR Rayonnement solaire incident NUIT Nébulosité [m/s] Fort Modéré Faible entre 4/8 et 7/8 <3/8 < A A-B B F F -3 A-B B C E F 3-5 B B-C C D E 5-6 C C-D D D D >6 C D D D D Tableau : Classes de stabilité de Pasquill Il est aussi possible de relier la rugosité moyenne du terrain z 0, et la longueur L aux classes de stabilité de Pasquill, (voir la Figure 7) [Lees]. Figure 7 : Relation entre la longueur de Monin-Obukhov L et la rugosité z 0 pour différentes classes de stabilité de Pasquill. /61\n25 Enfin, Doury fait appel à seulement deux classes de stabilité : les conditions de diffusion normale (DN). Cette situation est celle existant le jour, ou la nuit avec des vitesses de vent supérieures à 3 m/s ; les conditions de diffusion faible (DF). Cette situation est celle existant la nuit avec des vitesses de vent inférieures à 3 m/s...4 STABILITE ATMOSPHERIQUE ET CONDITIONS METEOROLOGIQUES..4.1 MECANISMES GENERAUX Il est intéressant d examiner les phénomènes météorologiques qui peuvent amener l'atmosphère à s'éloigner des conditions d'adiabatiques. Durant la journée, par temps clair, la répartition verticale des températures est proche de l'adiabatique ou supérieure en raison du fort réchauffement du sol par le soleil. Après le coucher du soleil, le sol se refroidit rapidement et refroidit par rayonnement les très basses couches de l'atmosphère. Par vents faibles, le refroidissement peut gagner un niveau de plus en plus élevé mais toutefois limité à quelques centaines de mètres au-dessus du sol, ce qui détermine une stabilité considérable des basses couches de l'atmosphère. Deux heures après l'aube environ, la chaleur du soleil brise l'inversion. La capacité qu'a le sol de perdre ou d'emmagasiner de la chaleur par rayonnement peut donc agir de façon très nette sur la répartition verticale des températures. Or, les propriétés de rayonnement du sol peuvent varier dans de larges limites suivant sa nature et sa couverture végétale. Les sols nus absorbent et perdent la chaleur rapidement. Par contre, les gains et pertes de chaleur se font de façon modérée pour les sols où la végétation est dense et humide. Un sol couvert de neige par exemple est plus favorable à la formation et à la persistance d'inversions thermiques qu'une région verdoyante ou boisée. La couverture nuageuse peut également contribuer à la stabilité des basses couches de l'atmosphère en interceptant le rayonnement solaire direct et diffus et en émettant son rayonnement propre de grandes longueurs d'ondes. En particulier, dans le cas d'une couche nuageuse continue et basse, le rayonnement de celle-ci peut neutraliser les effets de la turbulence dus au rayonnement du sol. La Figure 8 illustre les mécanismes d absorption des rayonnements solaires dans la couche limite atmosphérique. Figure 8 : Absorption du rayonnement solaire dans la couche limite atmosphérique 3/61\n26 La Figure 8 montre que les nuages réfléchissent 5 à 7% des rayons du soleil, tandis que l atmosphère absorbe à elle seule 16%. Environ 7% des radiations sont réfléchies par l atmosphère même sans nuage. Le reste (environ 50%) va directement au sol, ce dernier ne retournant que 3 à 5 % de ce qu il reçoit...4. CONDITIONS DE STABILITE NEUTRE Par vent fort ou par vent modéré et ciel couvert : des transferts importants de quantité de mouvement ont lieu dans la couche limite atmosphérique ; le gradient vertical de température ; z tend vers le gradient de température adiabatique le nombre de Richardson est nul. La stabilité atmosphérique est qualifiée de «neutre». Les rugosités de surface telles que les arbres, les bâtiments, les irrégularités topographiques conditionnent alors notablement l écoulement au voisinage du sol. Dans ces conditions, la turbulence produite est essentiellement d origine mécanique. Pour des vents faibles, la turbulence d origine thermique prédomine va influencer de façon appréciable la stabilité atmosphérique CONDITIONS DE STABILITE INSTABLE Au cours d une journée ensoleillée, avec des vents faibles, le sol s échauffe par rayonnement beaucoup plus rapidement que l air. Le gradient de température vertical négatif est inférieur au gradient de température adiabatique. Le nombre de Richardson est alors négatif. Des tourbillons convectifs apparaissent introduisant une instabilité qui a tendance à amplifier les mouvements verticaux. La stabilité atmosphérique est qualifiée «d instable» CONDITIONS DE STABILITE STABLE Au cours d une nuit claire, avec des vents faibles, le sol se refroidit par rayonnement plus rapidement que les couches d air adjacentes. Le gradient de température vertical est supérieur au gradient de température adiabatique. Les transferts turbulents sont faibles, particulièrement dans la direction verticale. Le nombre de Richardson est alors positif. Dans ces conditions, l atmosphère est stable CONDITIONS D INVERSION DE TEMPERATURE En fin de journée, à mesure que le soleil se couche, le sol se refroidit plus rapidement que les couches d air adjacentes. A proximité du sol, le gradient de température vertical est supérieur au gradient de température adiabatique. Par contre, plus en altitude, les couches d air sont plus froides à mesure que l on s élève (gradient de température vertical négatif inférieur au gradient de température adiabatique) du fait des conditions rencontrées au cours de la journée. L atmosphère est donc stable près du sol et instable plus haut en altitude. En début de journée, à mesure que le soleil se lève, le sol se réchauffe davantage et plus rapidement que les couches d air adjacentes. A proximité du sol, le gradient de température 4/61\n27 vertical est inférieur au gradient de température adiabatique. Par contre, plus en altitude, les couches d air sont plus chaudes à mesure que l on s élève (gradient de température vertical négatif supérieur au gradient de température adiabatique) du fait des conditions rencontrées au cours de la nuit. L atmosphère est donc instable près du sol et stable plus haut en altitude. Ces deux dernières situations sont qualifiées «d inversion de température»..3 ENVIRONNEMENT La dispersion atmosphérique d un produit est influencée par l occupation du sol, à la fois dans le champ proche avec la présence d obstacles isolés (bâtiments industriels, lotissements...), mais aussi dans le champ lointain avec l existence d accidents topographiques (vallées, falaises, collines, buttes...). Ces éléments perturbent le champ de vent de façon mécanique. Par ailleurs, les éventuelles discontinuités de recouvrement du sol vont être à l origine d effets thermiques propres à modifier le champ de vent avec la formation de mécanismes de brise..3.1 EFFETS MECANIQUES Les obstacles et les accidents topographiques perturbent la trajectoire du vent et modifient les caractéristiques moyennes et turbulentes de l écoulement de l air. Ces modifications sont, naturellement, fonction de la taille et de la forme des obstacles topographiques rencontrés par le vent. Lorsque les irrégularités au sol sont de faible taille par rapport à celle du nuage, les perturbations qu elles engendrent, affectent de façon globale la dispersion du nuage. L hypothèse d un terrain idéalement plat et de rugosité uniforme est généralement bien adaptée à cette problématique. Dans le tableau ci-après, des valeurs de rugosité sont données pour quelques types caractéristiques de surface. 5/61\n28 Désignation courante Exemple Rugosité (m) Terrain plat Aéroport 0,03 Terrain agricole Plaine avec arbres dispersés 0,1 Zone d'habitat dispersé Maisons éparpillées 0,3 Zone résidentielle Zone urbaine Sites industriels sans obstacles importants, forêts Ville, sites industriels ave bâtiments importants 1 3 Tableau 3 : Valeurs typiques de rugosité Ceci étant, lorsque des obstacles de la taille du nuage existent, que le recouvrement du sol n est pas uniforme ou que la topographie est accidentée, les mécanismes de dispersion sont plus complexes et réellement propres aux caractéristiques du site. A titre d illustration, il est intéressant d analyser la perturbation induite par ces éléments sur le champ de vent au travers de trois cas généraux avec : un site présentant un changement de rugosité ; un site avec une topographie accentuée ; un site avec un obstacle CHANGEMENT DE RUGOSITE Si la nature du sol présente une discontinuité et que la hauteur de rugosité passe brusquement de la valeur z 1 à z, alors la structure du vent se trouve profondément perturbée en aval du changement de rugosité (voir la Figure 9) [Turbelin G., 000]. La perturbation engendrée est fonction de l importance relative des longueurs z 1 et z, avant et après la discontinuité. 6/61\n29 Figure 9 : Ecoulement sur une surface à rugosité variable [Turbelin, 000] Ainsi, à l aval immédiat de la discontinuité, on va distinguer (voir la Figure 9) : une zone située loin de la surface, dans laquelle les propriétés du vent sont toujours reliées aux paramètres caractéristiques de la surface amont u 1 et z 1 ; une zone d épaisseur, appelée couche limite interne, dans laquelle le profil des vitesses est affecté par la discontinuité. Cette couche limite interne est constituée : d une zone d épaisseur située près du sol, appelée sous-couche interne, dans laquelle l écoulement est en équilibre vis a vis des nouvelles conditions aux limites ; d une zone de transition dans laquelle les vitesses passent du profil amont (loin du sol) au profil de la sous-couche interne (près du sol). et sont des fonctions croissantes. Enfin, loin en aval de la discontinuité, l écoulement retrouve un équilibre et les vitesses s alignent sur les caractéristiques de la surface en aval dépendant de u et z TOPOGRAPHIE ACCENTUEE Un écoulement qui aborde transversalement une colline subit une ascendance le long du versant au vent puis une descente le long du versant sous le vent. Au niveau du sommet, une zone de dépression et de survitesse se forme, tandis qu au pied du versant au vent, l air est ralenti (voir la Figure 10) [Turbelin G., 000]. 7/61\n30 Le rapport entre la vitesse moyenne incidente et la perturbation de vitesse au sommet est appelé rapport fractionnaire de survitesse. Ce paramètre est proportionnel à la pente moyenne de l obstacle H L. Quand la pente est faible, on distingue : une couche interne, dans laquelle les perturbations sont principalement liées aux transferts turbulents ; une couche externe, dans laquelle ils sont négligeables ; et une zone de sillage, située en aval de la topographie, dans laquelle la turbulence est importante. Si la pente ou la rugosité augmente fortement, une zone de recirculation se forme au pied du versant au vent et modifie les valeurs des survitesses. Pour certaines valeurs critiques du rapport H L, l écoulement décolle immédiatement après le sommet puis se recolle en aval. La structure de l écoulement dans la poche de décollement est très complexe et influence l ensemble des caractéristiques du vent. Figure 10 : Ecoulement au-dessus d une colline [Turbelin, 000] Plus particulièrement, s agissant de la dispersion d un panache de produit, l expérience permet d admettre que : en atmosphère instable ou neutre, le panache suit les variations du relief. Cela revient à considérer qu en tout point de l axe du panache, la distance verticale est constante, égale à la hauteur effective initiale du rejet ; en atmosphère stable, le panache reste horizontal et contourne le relief si la hauteur de celui-ci est supérieure à la hauteur effective du rejet. 8/61\n31 OBSTACLE La perturbation du champ de vent par la présence d obstacles isolés (constructions, haies...) va dépendre : de la géométrie et des dimensions de l obstacle ; de son état de surface ; des caractéristiques du vent incident ; de l environnement proche. Figure 11 : Ecoulement en présence d un obstacle isolé [Turbelin, 000] La complexité des écoulements autour d obstacles isolés peut être illustrée par l écoulement autour d un bâtiment de forme cubique (voir la Figure 11) [Turbelin G., 000] : au niveau des parois de l obstacle, se forment des couches limites qui se séparent et se rattachent, créant des zones de recirculation ; au niveau de la rencontre de l écoulement décollé et du sillage, se forme une couche fortement cisaillée ; au niveau de la face au vent, se forme un écoulement descendant qui, à la rencontre du sol, forme un rouleau tourbillonnaire. Ce rouleau se déplace en aval en contournant latéralement l obstacle, formant des tourbillons dit en fer à cheval. 9/61\n32 Relativement à la dispersion d un produit, on peut observer : 1) en amont de l obstacle : l accroissement de la concentration ; l augmentation du temps d exposition ; la présence de zones de recirculation ; ) en aval de l obstacle : la décroissance de la concentration ; l augmentation du temps d exposition ; la diminution de la largeur du nuage ; l homogénéité verticale de la concentration. La vitesse du vent et la densité du gaz rejeté jouent des rôles particulièrement importants. Ainsi, plus le vent est fort et le gaz léger, moins la présence de l obstacle se fait sentir. Les mécanismes de dispersion peuvent être considérés comme dégagés de l influence de tout obstacle au sol à des distances supérieures à 10 fois la dimension caractéristique de cet obstacle dans la direction du vent, ou à,5 fois perpendiculairement à cette direction à la fois dans le plan horizontal et dans le plan vertical (voir la Figure 1) [Hug, 1975]. Figure 1 : Influence d un obstacle sur l écoulement du vent [Hug, 1975] Pour un observateur éloigné, les turbulences induites par un bâtiment proche de la source engendrent une dilution initiale du rejet à peu près homogène dans le sillage de l obstacle. Le rejet peut alors être considéré comme ayant des dimensions initiales non nulles, de l ordre de celle de l obstacle. En termes de conséquences, si le gaz rejeté avait été toxique, en amont de l obstacle, les concentrations sont plus importantes et le temps d exposition plus grand, la situation se serait bien évidemment aggravée. En aval de l obstacle, la concentration est moins importante mais il faut souligner que l augmentation du temps d exposition ne permet pas de conclure à des effets de dose plus faibles dans tous les cas. Par ailleurs, si le nuage déborde sur les côtés de l obstacle, il faut parfois redouter la formation de filets de nuage à forte toxicité (la concentration en amont étant plus importante) de chaque côté de la zone du sillage de l obstacle. 30/61\n33 Dans le cas d un gaz explosible, la présence de l obstacle dans le champ très proche fait chuter la concentration en aval de l obstacle ce qui risque d augmenter pour certains gaz sensiblement le volume explosible à cette distance. De plus le passage de l obstacle génère une turbulence non négligeable dans le nuage, ce qui revient à le brasser. Ces deux conséquences peuvent dans de nombreux cas générer une explosion pouvant être encore amplifiée par la présence de l obstacle..3. EFFETS THERMIQUES Les mécanismes mis en évidence dans ce chapitre sont relatifs aux phénomènes dits de brise. La définition qui en est donnée ci-après est issue du "Dictionnaire du climat" Editions Larousse - Troisième édition, Brise n.f. Vent à caractère local ou régional et à alternance diurne, qui s'établit à proximité des lacs, des mers, dans les régions montagneuses, près des villes, etc., et qui résulte des différences de température dans les basses couches de l'atmosphère. Le terme «brise» est utilisé à la fois comme un synonyme de «vent faible» et pour désigner un courant d'origine thermique. Dans le premier cas, il est surtout employé dans la navigation maritime, où l'on distingue, en fonction de la force du vent, cinq types de brise correspondant aux niveaux à 6 de l'échelle de Beaufort. Dans la deuxième acception, la brise (ou brise thermique) est un vent généré par un gradient de pression horizontal. Ce gradient résulte des réchauffements et des refroidissements différenciés des diverses surfaces (terre, eau, prairie, forêt, immeubles...). Ces contrastes thermiques dépendent de la nature des surfaces et/ou de l'exposition des reliefs au soleil. Les littoraux. les campagnes, les montagnes et même les villes peuvent être à l'origine de ces courants thermiques dont les deux plus importants sont les brises de terre et de mer (ou de lac), dites «brises littorales», et les brises de vallée et de montagne, dites «brises orographiques» INFLUENCE DE LA PRESENCE D'UNE ETENDUE D'EAU PHENOMENES DE BRISE LITTORALE Au cours d'une journée, la terre se réchauffe plus rapidement que la mer. Sur la terre une basse pression s'établit localement propice à un écoulement d'air de l'étendue d'eau (mer, lac) vers la terre. Ce phénomène est appelé brise de mer (voir la Figure 13). Air chaud montant Air froid descendant Brise de mer Terre Figure 13 : Mécanisme conduisant à une brise de mer 31/61\n34 La nuit, c'est le phénomène inverse qui se produit. La terre perd sa chaleur plus rapidement que la mer. Un gradient de pression et de température s'établit lorsque la terre devient plus froide que la mer. Une zone de basse pression se crée alors sur la mer à l'origine d'un écoulement d'air de la terre vers la mer appelé brise de terre (voir la Figure 14). Air froid descendant Air chaud montant Terre Brise de terre Figure 14 : Mécanisme conduisant à une brise de terre Au crépuscule, il se produit souvent une période calme, les températures sur terre et sur mer étant plus ou moins égales. La brise de terre est en général moins forte que la brise de mer du fait d'une différence de température entre la terre et la mer moins importante la nuit. L'effet des brises de mer et de terre s'étend sur quelques centaines de mètres à quelques kilomètres près des côtes. Les lacs engendrent parfois un phénomène analogue, mais avec une extension et une intensité généralement plus limitées. Sous les latitudes tempérées, les brises apparaissent surtout en été lorsque le vent à l'échelle synoptique ne les contrarie pas. Lorsque ce vent souffle dans la même direction que la brise, on peut enregistrer une accélération locale du vent (le vent résultant étant la somme des deux flux)..3.. INFLUENCE D'UNE PENTE PHENOMENES DE BRISE DE PENTE Le relief peut être à l'origine d'un régime de vent particulier, dit brise de pente. La brise de pente est due à la différence de température entre le sol incliné d'un relief et l'air libre au même niveau. Durant le jour, sous l'influence du rayonnement solaire, le sol se réchauffe plus vite que l'air libre. L'air qui se trouve directement au contact du sol tend donc à s'élever le long de la pente tandis que l'air libre s'affaisse : c'est la brise montante ou anabatique (voir la Figure 15). Figure 15 : Brise montante 3/61\n35 La nuit, le sol refroidit l'air qui se trouve en contact avec lui. cet air, devenu plus lourd s'écoule le long de la pente : c'est la brise descendante ou catabatique (voir la Figure 16). Figure 16 : Brise descendante.3..3 INFLUENCE D'UNE VALLEE L'influence des vallées sur la circulation du vent est assez caractéristique. Plusieurs mécanismes interviennent sur les circulations locales du vent. Tout d'abord, les phénomènes de brise de pente se retrouvent ici naturellement du fait de la nature même des vallées (pentes de chaque côté). Ainsi, dans la journée, le fond des vallées et les versants directement exposés au soleil (adrets) se réchauffent plus que les versants à l'ombre (ubacs). L'air proche du sol se réchauffe au contact des adrets et tend à s'élever le long de ces pentes, donnant naissance à un vent de pente ascendant ou anabatique. La nuit, les pertes du sol par rayonnement infrarouge sont importantes. L'air de surface devient plus froid et donc plus lourd que l'air ambiant et il s'écoule par gravité vers la vallée donnant naissance à un vent de pente descendante ou catabatique (Voir la Figure 17). Figure 17 : Brises dans une vallée 33/61\n36 Les conditions locales (topographie, forme des vallées, exposition...) jouent de façon évidente un rôle majeur dans l'intensité des brises orographiques. Ces mécanismes interagissent avec le champ de vent général. A l'échelle régionale et locale, les reliefs canalisent les flux et sont la cause d'accélérations ou de ralentissements des vents. L'influence du champ de vent général est maximale lorsqu'il souffle dans l'axe de la vallée, ce qui accroît d'ailleurs sa vitesse du fait de la réduction de la section de passage offerte à l'écoulement, par les limites de la vallée. Lorsque le champ de vent général a une direction perpendiculaire à l'axe de la vallée, son influence sur la brise est réduite car sa vitesse dans la vallée est fortement réduite. 34/61\n37 3 MODELISATION 3.1 CONTEXTE La dispersion atmosphérique d un produit peut être étudiée selon différentes approches avec [Mouilleau, 1991b] : la réalisation d essais en grandeur réelle ; la simulation sur maquette (hydraulique ou aéraulique) ; l utilisation de codes de calcul mathématiques. L emploi de modèles numériques pour l évaluation de la dispersion atmosphérique d un produit présente de nombreux avantages par rapport à des expérimentations : la rapidité de l'étude et la possibilité d'envisager un grand nombre de cas. Toutefois, les essais à grande échelle ou sur maquette permettent aussi de connaître les phénomènes à modéliser et par conséquent de participer au développement des outils de calcul. Ce chapitre présente tout d abord de façon synthétique l ensemble de ces méthodes pour se concentrer ensuite sur l utilisation des outils de calcul destinés à évaluer la dispersion atmosphérique ESSAIS EN GRANDEUR REELLE Les essais en grandeur réelle peuvent être réalisés pour estimer la dispersion sur un site donné ou approfondir les connaissances sur certains mécanismes de dispersion. D une manière générale, ces essais sont relativement coûteux du fait de l importance des moyens dont il faut disposer (systèmes d acquisition de données, systèmes de rejet, quantités de gaz à rejeter, main d œuvre, site disponible pour la réalisation de ces essais). A ces contraintes, il convient d ajouter l impossibilité de maîtriser les situations météorologiques. Des campagnes à caractère plus fondamental ont été faites pour l essentiel en Grande Bretagne (Porton Down, Thorney Island, Maplin Sands...) et aux Etats Unis (China Lake, Desert Tortoise...). A titre d illustration, la campagne de Thorney Island représente 9 rejets réalisés pendant un an et demi. Cette campagne a permis de collecter plus de mesures par essai recueillies par environ 00 capteurs de concentrations, de vitesse, de température, de pression [CCPS, 1996 ; Lees]. Pour un certain nombre de ces sites de production, EDF a réalisé des campagnes de traçage à l hexafluorure de soufre (SF 6 ) pour évaluer l influence des conditions météorologiques et du relief sur la dispersion des gaz passifs. Enfin, l INERIS a mené entre 1996 et 1997, des essais expérimentaux destinés à mieux appréhender la phénoménologie de la dispersion de l ammoniac [Bouet, 1999]. 35/61\n38 3.1. ESSAIS A ECHELLE REDUITE Les essais à échelle réduite se font généralement sur une maquette placée dans une soufflerie ou une veine hydraulique. L échelle de travail est en général de l ordre du 1/50 e. Pour représenter un phénomène identique à ce que l on observait dans la réalité, il faut s assurer de la conservation de grandeurs adimensionnelles de façon à garder le même système d équations pour les essais en grandeur réelle et ceux à échelle réduite OUTILS DE CALCUL La modélisation numérique présente a priori des avantages certains par rapport aux techniques physiques que sont les essais in situ ou à échelle réduite : rapidité de l étude ; possibilité d envisager un grand nombre de situations. L évolution des outils numériques de dispersion s est faite parallèlement à la capacité des ordinateurs, même si certaines techniques de calcul sophistiquées ont des origines assez anciennes. Les modèles numériques de dispersion peuvent se ranger par ordre de complexité croissante en trois principales familles : les modèles gaussiens ; les modèles intégraux ; les modèles CFD (Computational Fluid Dynamics). Globalement, les deux premiers types d outils s attachent à modéliser la dispersion à partir d équations paramétrées et simplifiées. Les temps de calcul sont courts, de l ordre de la minute. Ils peuvent être mis en œuvre à partir d un matériel informatique de type PC. La troisième famille résout directement le système d équation décrivant les mécanismes physiques de la dispersion. Les temps de calcul sont plus longs, de l ordre de l heure, voire de la journée. Ce type d outil fonctionne généralement sous station de calcul. Les premiers modèles de dispersion utilisés sur ordinateur sont de type Gaussien. Leur origine remonte aux travaux de Sutton (193), de Pasquill (1961, 1974), de Gifford (1961, 1968) et en France, de Le Quinio (1964, 1975) et de Doury (1977, 1986). Ces modèles permettent d estimer la dispersion des gaz neutres ou passifs (de masse volumique très proche de celle de l air). Cette limitation s est vite avérée inacceptable pour bon nombre d études faisant intervenir des gaz dont la masse volumique était sensiblement plus importante que celle de l air. Van Ulden reprit en 1974, une méthodologie déjà existante et créa le premier modèle de type intégral appliqué à la dispersion des gaz lourds. L augmentation et la démocratisation croissante des moyens de calcul ont permis d envisager de résoudre directement le système d équations qui régit réellement la dispersion du nuage. C est ainsi que sont apparus depuis la fin des années 70, les modèles 36/61\n39 tridimensionnels permettant de simuler le comportement des nuages de gaz plus lourds que l air. Un domaine d application maintenant en développement est l intégration de ces outils dans des systèmes de surveillance et de prévision en temps réel DETERMINATION DES EFFETS SUR LA SANTE HUMAINE Les modèles déterminent les lieux de l'espace où une concentration donnée est observée. Il est donc possible, en première approximation, de ne s'intéresser qu'à la distance maximale sous le vent de la source où une concentration donnée est atteinte. Toutefois, d une manière générale, quel que soit le phénomène considéré (dispersion d'un gaz toxique ou inflammable), les conséquences d un accident industriel sont entre autres quantifiées en termes de distances limites en deçà desquelles il pourrait être observé des effets sur la santé humaine irréversibles et létaux 10. Pour cela, dans le cas d un rejet de gaz inflammable, il faut apprécier : la distance où est atteinte la limite inférieure d explosibilité (LIE) ; celle où est atteinte la limite supérieure d explosibilité (LSE) ; ainsi que la masse de gaz inflammable (concentration comprise entre la LIE et LSE).. Lorsqu'une personne respire une atmosphère polluée par un produit toxique, les effets auxquels on s'intéresse habituellement sont définis comme étant l'apparition (à faible probabilité) de la létalité, des malaises, de la toux... Ces effets sont directement fonction de la concentration C et du temps t pendant lequel le sujet est exposé à cette concentration. Les effets varient bien évidemment selon que la personne est un enfant ou un vieillard, ainsi qu'en fonction de son état de santé ou de ses facultés d'accoutumance. Aussi les courbes dans un plan (C, t) correspondant à un effet donné sont-elles établies pour une population représentative de l'ensemble des situations susceptibles de se produire. Dans le plan Log (C), Log (t), ces courbes sont assimilables à des droites et donc, à effet E constant, la concentration et le temps se trouvent pratiquement liés par une relation C n.t = E. Sur ces courbes, les coordonnées d'un point (C, t) représentent : l'échelon de concentration C, et le temps d'application t de cet échelon, nécessaires pour que l'effet E se produise. Cet effet se produira dès que C n.t E. En pratique lors d'un accident, un observateur n'est jamais soumis à un échelon de concentration constant dans le temps et il y a donc lieu d'intégrer les apports de chacun des pas de temps pendant lesquels la concentration est supposée constante en calculant : 10 Les seuils de référence à retenir pour les installations classées sont donnés dans l arrêté du /10/004 paru au J.O. du 19/1/ /61\n40 L'effet se produira si I E. n I C ( t) dt Les valeurs des variables E et n dépendent à la fois du type d effet considéré (apparition à faible probabilité de la létalité par exemple) et de la nature du produit toxique en cause. 3. MODELES NUMERIQUES 3..1 MODELES GAUSSIENS GENERALITES Un gaz passif va se disperser du fait de la seule action du fluide porteur, l air (voir le chapitre Densité du gaz). Le transport et la diffusion du gaz vont alors dépendre du vent et de la turbulence atmosphérique d origine mécanique ou thermique. En considérant la diffusion moléculaire négligeable 11 la diffusion turbulente homogène et isotrope et un champ de vent uniforme dans l'espace, la concentration de produit est supposée suivre une distribution gaussienne le long des plans perpendiculaires à la direction du rejet (voir la Figure 18). Figure 18 : Répartition gaussienne de la concentration dans un panache de gaz passif [Turner, 1970] DOMAINE DE VALIDITE Le modèle gaussien s applique aux rejets de gaz passifs, le produit rejeté doit donc avoir : une densité à peu près égale à celle de l'air (ou bien il est très dilué) ; 11 La diffusion moléculaire de 10 3 à 10 4 inférieure à la diffusion turbulente. 38/61\n41 une température identique à celle de l'air ; une vitesse initiale relative nulle. De plus, les modèles gaussiens s appuient sur un champ de concentration tridimensionnel généré par une source ponctuelle. La diffusion moléculaire est négligée devant la diffusion turbulente, la vitesse du vent doit donc être d au moins 1 à m/s. Dans la plupart des cas, le champ de vent est uniforme avec un profil vertical constant. Des outils plus élaborés peuvent considérer la variation de la vitesse du vent, de la température de l air et de la densité de l atmosphère avec l altitude. Une couche d inversion de température peut parfois aussi être prise en compte. La turbulence atmosphérique est prise en compte par l'intermédiaire de classes de stabilité (voir le chapitre..3.3 Degré de stabilité de l atmosphère). Le terrain doit être homogène et plat puisque la présence de reliefs, d'obstacles (murs, bâtiments ) introduirait des perturbations de l'écoulement de l'air importantes. Aussi, de façon pratique, les résultats sont valables au-delà d au moins 100 m depuis le point de rejet. Par ailleurs, au-delà de distances de dispersion de l ordre de la dizaine de kilomètres, les résultats ne sont plus valables car d autres phénomènes de turbulence et de diffusion doivent être considérés METHODES DE RESOLUTION Rejet instantané Dans le cas du rejet ponctuel et instantané d une masse de gaz M, la concentration C du gaz dans l atmosphère en un point (x,y,z) est de la forme : C(x, y,z,t) () 3 M x y z (x x 0 ut) exp x (y y 0 ) y (z z exp 0 ) z (z z exp 0 ) z - C : concentration [kg/m 3 ] ; - M : masse de produit libéré [kg] ; - x y z :coordonnées de la source de produit [m] ; x - y z : coordonnées du point où l on calcule la concentration [m] ; - u : vitesse moyenne du vent [m/s] ; - t : temps depuis l émission du gaz ; x - y z : écarts-types de la distribution gaussienne de la quantité M de gaz par rapport à sa localisation à l instant t [m] ; 39/61\n42 - : coefficient de réflexion au sol 1. Rejet continu Modèle «à bouffées» Pour modéliser un rejet continu, il est possible de considérer l émission d une succession de rejets instantanés qui évolueront de façon gaussienne. Le rejet de débit en fonction du temps Q(t) peut se décomposer en n rejets instantanés de masse M tels que : i ti1 ti Mi Q (ti ti1) Le i ème rejet instantané, indicé i, est émis à l instant t i et a une masse M i. La concentration C du gaz dans l atmosphère en un point (x,y,z) est alors de la forme : C(x, y,z,t) i n i 1 () 3 M i xi yi zi exp x x 0 u (t t ) x i i (y y 0 ) y i (z z exp 0 ) z i (z z exp 0 ) z i - M i : Masse du i ème rejet instantané [kg] ; - u : vitesse moyenne du vent [m/s] ; - n : nombre de rejets instantanés considérés ; - t i : instant de fin d émission du i ème rejet [s] ; - t i-1 : instant de fin d émission du (i-1) ème rejet et début d émission du i ème rejet [s] ; - xi yi zi : écarts-types de la distribution gaussienne du i ème rejet instantané de masse M i par rapport à sa localisation à l instant t [m] ; - : coefficient de réflexion au sol. 1 Un terme de réflexion au sol peut être pris en compte pour modéliser la capacité de réflexion ou d absorption du produit sur le sol, l eau ou les végétaux. De façon évidente, ce phénomène est plus ou moins accentué suivant la nature du produit et l environnement. 0 : absorption totale et 1 : réflexion totale pour un sol non poreux avec un produit ne pouvant réagir avec ce sol ou la végétation (un gaz comme l air sur du béton). 40/61\n43 Rejet continu Modèle «panache» Si le rejet est continu (voir le chapitre.1.4 Mode d émission) et si la vitesse de vent est suffisamment importante pour que le phénomène de diffusion dans la direction du vent soit négligeable devant le phénomène de convection, la concentration C du gaz dans l atmosphère en un point (x,y,z) peut s écrire : Q C(x, y,z) u Calcul des écarts types (y y exp (z z exp (z z exp 0 ) 0 ) y z y z L utilisation des modèles gaussiens impose donc la détermination des écarts-types [Doury ; Hanna, 198b ; Pasquill, 1974 ; Turner, 1969]. A partir de campagnes de mesures sur des rejets de produits, des formules ont été établies, donnant la valeur des écarts-types en fonction : de la distance d éloignement de la source de rejet (Pasquill, Turner, Briggs, Hosker) ; ou encore du temps de transfert (Doury). Les écarts types ont été calés pour des rejets relativement importants à la source (au minimum de l'ordre de la tonne), et pour des distances de dispersion de l'ordre de la dizaine de kilomètres au maximum. Les écarts types de la loi gaussienne dépendent : de la distance à la source ou de la durée de transfert ; des caractéristiques de la structure de l atmosphère ; et de la rugosité du site. Corrélation de Pasquill-Turner Les écart-types de Pasquill ont été déterminés pour les modèles «panache» à partir des mesures de «Prairies Grass» menées en Grande Bretagne, en Ces mesures ont été faites sur un terrain plat, bien dégagé, peu rugueux (de l ordre de 3 cm). En fonction des classes de stabilité atmosphérique de Pasquill (voir le chapitre..3.3 Degré de stabilité de l atmosphère), les valeurs des écarts-types ont été formulées en fonction à la distance de la source (valable pour des distances supérieures à 100 m et inférieur à 10 km). Les valeurs des écarts types y et z 0 ) z proposées par Pasquill correspondent à des durées d échantillonnage de 10 minutes, et une hauteur de source qui n excède pas les premières centaines de mètres. La relation pour calculer les écarts type est de la forme : a x b c Les valeurs de a, b et c sont reportées dans les tableaux suivants pour respectivement y et z (x, y et z sont exprimés en km). 41/61\n44 Stabilité atmosphérique (Pasquill) a b c A 0,15 0,858 B 0,155 0,889 C 0,105 0,903 D 0,068 0,908 0 E 0,05 0,914 F 0,034 0,908 Tableau 4 : Coefficients relatifs à y Stabilité atmosphérique (Pasquill) a b c A 0,467 1,89 0,01 B 0,103 1,11 0 C 0,066 0,915 0 D 0,0315 0,8 0 E si x < 1 km 0,03 0,745 0 E si x > 1 km 0,148 0,15-0,16 F si x < 1 km 0,0144 0,77 0 F si x > 1 km 0,031 0,306-0,017 Corrélation de Davidson Tableau 5 : Coefficients relatifs à z Pasquill est par ailleurs associé à d autres auteurs (Turner, Gifford...) qui ont établi des formules plus ou moins proches. L une de ces formulations est celle de Davidson (90). La relation de Davidson (90) pour estimer les écarts types est de la forme : a x bc.log x Les valeurs de a, b et c sont reportées dans les tableaux suivants pour respectivement y et z (x est exprimé en km et y, z en m). 4/61\n45 Stabilité atmosphérique (Pasquill) a b c A 09,6 0,8804-0,00690 B 154,7 0,893-0,00671 C 103,3 0,911-0, D 68,8 0,911-0, E 51,05 0,911-0, F 33,96 0,911-0, Tableau 6 : Coefficients relatifs à y Stabilité atmosphérique (Pasquill) a b c A 417,9,058 0,499 Pour plus de précisions avec la classe A si 0,10 < x < 0,45 km Pour plus de précisions avec la classe A si 0,45 < x < 3,1 km 310,4 1,773 0, ,9,117 0 B 109,8 1,064 0,01163 C 61,14 0, D 30,38 0,7309-0,03 E 1,14 0,680-0,045 F 13,7 0,6584-0,05367 Tableau 7 : Coefficients relatifs à z Corrélation de Briggs Une autre corrélation des écarts types basée sur les essais de Pasquill est celle de Briggs. La relation pour estimer les écarts types est de la forme : a x bc.log x Les valeurs de a, b et c sont reportées dans les tableaux suivants pour respectivement y et z (x est exprimé en km et y, z en m). 43/61\n46 Stabilité atmosphérique (Pasquill) y (m) z (m) 1 A 0.x x 1 B 0.16x x 0.0x 0.1x C x x x x D x x x x E x x 0.03x x 1 F x x 0.016x x 1 Corrélation de Doury Tableau 8 : Coefficients relatifs à y et z Les écart-types de Doury ont été déterminés à partir d un grand nombre de données expérimentales recueillies par Sandia Corporation, en En fonction des classes de stabilité atmosphérique de Doury (voir le chapitre..3.3 Degré de stabilité de l atmosphère), les valeurs des écarts-types sont fonction du temps de transfert et de la stabilité atmosphérique. La relation pour calculer les écarts type est de la forme : Kh y (A ht) Kz z (A zt) Les valeurs de A h, A z, Kh et Kz sont reportées dans les tableaux suivants pour respectivement y et z (t est exprimé en s et y, z en m). Diffusion normale (DN) Temps de transfert [s] A h A z [m 1/k /s] k= Kh ou Kz Kh Kz 0 < t < 40 0,405 0,4 0,859 0, < t < ,135 1,00 1,130 0, < t < ,135 0,0 1,130 0, < t < ,463 0,0 1,000 0, < t < ,50 0,0 0,84 0,500 t > ,0 0,500 0,500 Tableau 9 : Coefficients relatifs aux écarts-types 44/61\n47 Diffusion faible (DF) Temps de transfert [s] A h A z [m 1/k /s] k= Kh ou Kz Kh Kz 0 < t < 40 0,405 0,0 0,859 0, < t < ,135 0,0 1,130 0, < t < ,463 0,0 1,000 0, < t < ,50 0,0 0,84 0,500 t > ,0 0,500 0,500 Tableau 10 : Coefficients relatifs aux écarts-types LIMITATIONS ET CORRECTIONS POSSIBLES Définition des écarts-types relativement à la structure turbulente de l atmosphère La pertinence d un modèle gaussien dépend principalement d une estimation correcte des écarts-types. Les modèles de type gaussien ne donnent pas nécessairement tous les mêmes résultats du fait des différentes approches retenues pour la détermination des écarts-types. Aussi, certaines considérations relatives à la caractérisation de la turbulence atmosphérique doivent être gardées en mémoire pour le choix de la paramétrisation des écarts-types vis-àvis du problème étudié. Ainsi, le spectre de turbulence est décrit par des formulations empiriques qui s appuient sur des considérations de similitude. Les lois de similitude prévoient que l énergie turbulente augmente avec la vitesse du vent. Lorsque la vitesse du vent tend vers zéro, la théorie dit que la turbulence tend vers zéro. Néanmoins, les lois de similitude ne s appliquent pas à des turbulences de grande échelle. Même si la vitesse du vent tend vers zéro, une turbulence résiduelle constituée de basses fréquences persiste dans l atmosphère. Cette turbulence ne dépend pas de paramètres locaux. De ce fait, pour la turbulence de petite échelle, toutes choses égales par ailleurs, lorsque la vitesse du vent augmente, l énergie de mélange est plus importante, mais en contre-partie il faut moins de temps pour atteindre un point donné. Ces deux effets se compensent, sur une distance de transfert donnée, la diffusion reste la même quelle que soit la vitesse du vent. Autrement dit, l écart-type de la répartition du produit dépend de la distance de transfert et non pas de la vitesse à laquelle est parcourue cette distance. Lorsque la turbulence de grande échelle prédomine, c est-à-dire dans le cas de distances de transfert importantes (la diffusion de nuage de grandes dimensions est régie par des structures importantes) ou des situations de vent faible (la turbulence de petite échelle est faible), les écarts-types dépendent de ce fait de façon évidente du temps de transfert et non de la distance. 45/61\n48 Autrement dit, pour les temps de transfert courts, la turbulence de petite échelle joue un rôle non négligeable. Par contre pour les temps de transfert importants, la turbulence de grande échelle, indépendante de la hauteur de rugosité, devient dominante. Ainsi, dans la direction horizontale, le processus de diffusion comprend deux phénomènes [Bakkum, 1997] : la diffusion relative du produit autour du centre de masse de la bouffée (cas a) ; la diffusion des centres de masse des bouffées de produit (cas b). Selon la dimension de la bouffée, les hautes et basses fréquences du spectre de turbulence seront ou non actives sur la dispersion de la bouffée. Par ailleurs, lorsque la bouffée est dans un champ turbulent de dimension comparable à celle de la bouffée, les deux processus sont actifs (cas c). Les écarts-types sont fonction de la distance près de la source et dépendent uniquement du temps de transfert loin de la source. Pour des temps de transfert importants, l écart-type horizontal devient indépendant de la stabilité atmosphérique car la turbulence de grande échelle prévaut. Figure 19 : Dispersion d une bouffée de produits transportée par l air, dans trois conditions de turbulence [Bakkum et. al ;, 1997] 46/61


Vrai/faux : toutes les idées reçues sur l'eczéma - Le Parisien\nVrai/faux : toutes les idées reçues sur l'eczéma\nSi l'eczéma entraîne des démangeaisons bien souvent interminables, cette maladie est également source de nombreuses idées reçues. Décryptage avec des professionnels de santé pour enfin reconnaître le vrai du faux.\nLe 2 juillet 2017 à 13h09\n1. L'eczéma est une maladie fréquente\nVrai. La dermatite atopique ou eczéma est la deuxième maladie dermatologique chronique la plus fréquente en France après l'acné, avec une prévalence de 4,65%. Elle touche 20% des enfants de moins de 7 ans et 18% des enfants âgés de 7 à 16 ans.\n2. L'eczéma est une maladie contagieuse\nFaux. L'eczéma ne peut pas se transmettre par le toucher. Il s'agit d'une prédisposition génétique.\n3. L'eczéma est dû à un manque d'hygiène.\nFaux. Au contraire, l'eczéma peut être aggravé par un excès de soins d'hygiène. Le contact prolongé avec la peau et l'utilisation de produits non adaptés peuvent être responsables de plaques d'eczéma.\n4. Les poussées d'eczéma disparaissent souvent avec le temps\nVrai et Faux. L'eczéma du bébé s'améliore en général avec le temps. La plupart guérissent spontanément en grandissant. L'eczéma atopique de l'adulte peut persister depuis l'enfance, ressurgir après une période d'accalmie ou apparaître à l'âge adulte.\n5. Les crèmes à la cortisone soignent l'inflammation de la peau\nVrai. Les dermocorticoïdes ont une activité anti-inflammatoire dès l'apparition des rougeurs. L'inflammation de la peau lors d'une poussée d'eczéma doit être traitée sous peine de s'aggraver ou de se surinfecter.\n6. L'eczéma peut entraîner de l'asthme\nFaux. Le terrain atopique, qui se manifeste par une hypersensibilité aux allergènes, peut s'exprimer sous plusieurs formes : l'eczéma, la rhinoconjonctivite allergique ou l'asthme. Ce n'est donc pas l'eczéma qui provoque l'asthme. La maladie est également associée à des allergies alimentaires.\n7. Il y a forcément un risque de rebond à l'arrêt des dermocorticoïdes\nFaux. On parle d'effet rebond quand les patients observent que l'arrêt de l'application déclenche une poussée de l'Eczéma. Cet effet rebond est causé par une mauvaise utilisation du traitement. Pour savoir si l'on peut arrêter les dermocorticoïdes, il suffit de regarder la peau. Si celle-ci est revenue à sa couleur d'origine et qu'elle est douce, c'est bon.\n8. La bonne astuce naturelle pour soulager l'eczéma\nTrès efficaces, l'huile d'argan et l'huile de bourrache seraient des alliées exemplaires pour lutter contre les poussées d'eczéma. En effet, appliquées directement sur la peau (ou prises sous formes de gélules) elles atténuent le prurit et la rougeur des lésions.


Reprise tendue des discussions sur le nucléaire iranien à Moscou\nMonde Reprise tendue des discussions sur le nucléaire iranien à Moscou\nMis à jour le 18/06/12 à 05h39\nL'Iran et six puissances internationales se retrouvent ce lundi à Moscou pour tenter de maintenir en vie les efforts de règlement par la voie diplomatique de la crise liée au programme nucléaire iranien. La menace d'enlisement, voire de rupture, des discussions alimente l'hypothèse d'une intervention militaire d'Israël, qui menace de bombarder les sites nucléaires iraniens en cas d'échec des efforts diplomatiques.\nAprès 15 mois d'interruption, des pourparlers ont été relancés en avril entre l'Iran, d'une part, et le groupe P5+1 (Etats-Unis, Russie, Chine, France, Grande-Bretagne et Allemagne), d'autre part. Lors de ce troisième cycle de discussions à Moscou, le groupe P5+1 va insister sur son principal motif de préoccupation: l'enrichissement d'uranium à 20% par l'Iran, qui constitue un pas supplémentaire vers une possible utilisation militaire de ce combustible.\nAprès deux séances de discussions guère fructueuses, à Istanbul en avril puis à Bagdad en mai, les six puissances internationales sont aussi soucieuses d'obtenir des résultats afin de dissiper l'impression que l'Iran cherche avant tout à gagner du temps avec ce processus diplomatique.\n«Nous sommes très déterminés à poursuivre ce processus tant qu'existe un élan et tant qu'il y a un engagement (de l'Iran) à discuter du fond de la question nucléaire», a dit un diplomate européen de haut rang.\nReprésentés par Catherine Ashton, haute représentante de la diplomatie de l'Union européenne, les Six espèrent au moins obtenir l'assurance que Téhéran est disposé à aborder des questions concrètes, ce qui permettra d'entretenir les discussions.\n08:05 23 commentaires 9 partages\nTour de France 2019 EN DIRECT: Clarke en tête à l'entame de la dernière montée... Tout est ouvert pour la victoire... Suivez la 12e étape avec nous\n13:15 15 commentaires 11 partages


C4 carbon remover - Club 50 60 - 22hunter\n> C4 carbon remover\nBore Tech C4 Remover Carbon Bore solvant de nettoyage, pour éliminer le carbone de la combustion de la poudre. Nettoie efficacement, les endroits difficiles comme les tuners. Idéal pour les canons, les tuners, les silencieux, les gorges de chambre, les cylindres de revolver, et tout endroit où l'encrassement au carbone est présent.\nRéférence : jar-002\nFormulation biodégradable non-inflammable, non dangereux et 100% exclusive de la C4 Carbon Remover repose sur un savant mélange de tensioactifs, des agents de pénétration, des agents de suspension et sur l'action chimique des produits pour enlever les dépôts de carbone tenaces.\nIdéal pour les canons, les tuners, les silencieux, gorges de la chambre, les cylindres de revolver, et tout endroit où l'encrassement au carbone est présent.\nPar JOSÉ V. le 04 Avr. 2014 :\nCommentaire : M'a permis de venir à bout d'un anneau de carbone sans m'énerver ! Il a éliminé les petites traces noires toujours présentes aux angles des rayures qui me font toujours penser que je n'ai pas fini !! J'en suis content, j'alterne avec du Rimfire Blend pour les traces de plomb.\n(C4 carbon remover)\nC4 carbon remover - Bore Tech Inc - Bore Tech C4 Remover Carbon Bore solvant de nettoyage,... - Code GTIN : 0667739980050 - Référence produit : jar-002 - Prix : 21,45 € - Etat : Neuf - Stock En stock\nNote(s) et avis : C4 carbon remover - Moyenne: 5 / 5 - Basé sur 1 note(s)