Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
Lakes are major depositional systems for which the related depositional processes have long been considered relatively simple. Breaking this statement, this study presents a detailed analysis of deposits in Lake Saint-Jean, the third largest natural lake in Québec. In addition to postglacial deltaic and coastal depositional systems fringing the lake, current-controlled features such as a large subaqueous prograding wedge and three sediment drifts have been identifi ed in its central portion based on two-dimensional (2-D) acoustic high-resolution subbottom profiles. The large subaqueous prograding wedge is a 4-km-long and up to 15-m-thick heterolithic shelf-like construction in the southeastern part of the lake. The three sediment drifts are 0.1–0.5-km-long and 2–5-mthick mud mounds distributed on the lake floor in the central portion of the lake. Diatom analyses and radiocarbon dating show that the development of these current-controlled features occurred during the lacustrine phase, after the disconnection with the postglacial marine Laflamme Gulf at 8.5 cal. k.y. B.P. Depositional facies show evidence of recurrent bottom-current activity. Related deposits alternate with pelagic sedimentation stages characterized by the settling of mud and biogenic accumulations. We investigated the origin of bottom currents using a numerical simulation (SYMPHONIE, an oceanographic model), with the aim of modeling wind-induced lake-scale water circulation. Simulations suggest that the subaqueous prograding wedge and the three sediment drifts result from wind-induced bottom currents generated by storm events having wind speed greater than 10 m s–1. Such strong winds are able to significantly affect sedimentation in the central portion of Lake Saint-Jean. Theresulting wind-induced sedimentary features were integrated into a refi ned lacustrine depositional model that summarizes the evolution of a group of water bodies referred to as “wind-driven water bodies.” This study applies a new tool for lake strata characterization and highlights the potential diffi cultyin differentiating them from marine deposits in the geological record.
ISSN: 0016-7606

insu-01184724
https://hal-insu.archives-ouvertes.fr/insu-01184724
DOI : 10.1130/B31145.1

Éditeur(s) : HAL CCSD Wiley
Résumé :  
Aims: Astrocaryum sciophilum (Miq.) Pulle (Arecaceae) is an understorey palm, endemic to north-eastern South America with a patchy distribution. We tested the hypothesis that the spatial distribution of this palm species is not in equilibrium but is slowly colonizing the forest understorey. Location: Inventories and seed dispersal studies were conducted in the undisturbed tropical forest close to the Nouragues research station, French Guiana. Additional data were collected in the entire territory of French Guiana. Methods: We studied the demography of A. sciophilum on a 20-ha plot located at the edge of its distribution. The age of the palms was estimated by postulating an exponentially decreasing abundance by age class. Direct seed dispersal experiments were also conducted, to estimate dispersal parameters. The seeds of A. sciophilum were dispersed only by rodents. This information was used to parameterize a forest growth simulator, to study the spatial spread of this species. Results: Within the sampling plot, the density of A. sciophilum dropped sharply from about 500 individuals per hectare to zero. The maturation age was estimated to be 170 ± 70 years, and over 55 years with 95% confidence. Seed-dispersal experiments yielded an average seed dispersal distance of 11 m and a maximum estimated dispersal distance of 125 m across a generational span of 55 years to maturity. Therefore, the maximal estimated colonization speed is 2.3 m/y. Conclusions: Empirical results and numerical simulations suggest that the boundary of the A. sciophilum population is a colonization front, and that the range of this species is slowly expanding. The implications of this result in respect of palaeoenvironmental changes in this region are discussed.
ISSN: 1466-822X

hal-01032154
https://hal.archives-ouvertes.fr/hal-01032154
DOI : 10.1046/j.1466-822X.2003.00020.x

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
The association of experimental data showing that the plastic deformation of olivine, the main constituent of the upper mantle, is highly anisotropic and the ubiquitous seismic anisotropy in the upper mantle, which indicates that olivine crystals show coherent orientations over scales of tens to hundreds of kilometers, implies that the long-term deformation in the upper mantle is anisotropic. We propose a multiscale approach, based on a combination of finite element and homogenization techniques, to model the deformation of a lithospheric plate while fully considering the mechanical anisotropy stemming from a strain-induced orientation of olivine crystals in the mantle. This multiscale model explicitly takes into account the evolution of crystal preferred orientations (CPO) of olivine and of the mechanical anisotropy during the deformation. We performed a series of numerical experiments simulating the uniaxial extension of a homogeneous (100% olivine) but anisotropic plate to test the role of the olivine CPO on the plate mechanical behavior and the link between CPO and mechanical anisotropy evolution. Even for this simple solicitation, different orientations and intensity of the initial olivine CPO result in variable plate strengths and deformation regimes. A plate with an initial CPO where the olivine [100] and [010] axes are concentrated at 45 degrees to the extension direction has high resolved shear stresses on the easy (010)[100] and (001)[100] slip systems of olivine. This results in low strength and in deformation by transtension. Plates with an initial CPO where the maximum of [100] axes is parallel or normal to the extension direction show a high initial strength. Isotropic plates have an intermediate behavior. The progressive rotation of olivine [100] axes toward the imposed stretching direction results in hardening in all models, except in those characterized by an initial concentration of olivine [100] axes normal to the imposed extension, in which softening is followed by hardening.
ISSN: 1525-2027

hal-00420063
https://hal.archives-ouvertes.fr/hal-00420063
DOI : 10.1029/2009GC002423

Éditeur(s) : Université Paris 7 – Denis Diderot
Résumé : Sea surface salinity (SSS) influences numerous oceanic phenomena, for instance surface water ventilation, deep water formation and thermohaline circulation. SSS also controls some ocean-atmosphere coupled processes, such as the intensity of freshwater flux and the penetration of heat flux and turbulence. Salinity is more difficult to measure than temperature from in situ surveys, which results in 20 times less data being currently available. Moreover, sea surface temperature (SST) is routinely estimated from satellites, which is not possible yet for SSS. Two space missions will fill this gap soon: SMOS from the European Space Agency and Aquarius/SAC-D from NASA and CONAE. To contribute to the SMOS project, we propose a method for estimating SSS from current satellite observations and for studying the mechanisms governing its variability. We developed a simplified model of the ocean mixed layer, based on the "slab mixed layer" formulation (Frankignoul et Hasselmann, 1977). This 2D horizontal model is implemented over the global ocean, using a near 100 km resolution, and integrated during a climatological year. Air-sea fluxes are taken from the ECMWF meteorological model (ERA40 reanalysis) and the surface currents are provided by altimeter data (SSALTO-DUACS analysis). The mixed layer depth (MLD) is derived from SST observations, using an original inversion technique. The MLD fields obtained from this inversion are well correlated to in situ estimates. This effective depth represents the air-sea fluxes penetration and ensures consistency between fluxes, les currents and SST. We first validate the simulation through examination of the heat balance in north-eastern Atlantic, by comparing to measurements and models from the POMME experiment. Then we study the salinity balance in the global domain, in terms of its geographical distribution and seasonal evolution. Equilibrium between the various processes appears generally more complex than for temperature. Noteworthy, the role of atmospheric flux is less predominant (22%), while geostrophic advection (33%) and diapycnal mixing (22%) contribute more strongly. Our results indicate this model succeeds in reconstructing SSS variability over most of the oceans. Daily SSS variations are also simulated, whereas they are not represented in current observed data at a global scale. Owing to its simplicity and fast computation, the model will be useful to the SMOS mission. It can help for the measurement calibration/validation and provide a first guess estimate to the sophisticated algorithm required for SSS restitution.
La salinité de surface des océans (SSS) influence de nombreux phénomènes océaniques, parmi lesquels la ventilation des eaux de surface, la formation d'eaux profondes et la circulation thermohaline. Elle détermine aussi certains processus couplés océan-atmosphère, notamment l'intensité du flux d'eau douce, la pénétration du flux de chaleur et de la turbulence. La mesure in situ de salinité est plus compliquée que celle de température, si bien qu'on dispose actuellement de 20 fois moins de données pour cette propriété. De plus, la température de surface (SST) est couramment estimée par satellite, ce qui n'est pas encore le cas de la SSS. Deux missions spatiales vont prochainement combler cette lacune : le satellite SMOS de l'Agence Spatiale Européenne et le satellite Aquarius/SAC-D de la NASA et de la CONAE. En préparation du projet SMOS, nous proposons une méthode pour estimer la SSS à partir des observations satellitaires actuelles et étudier les mécanismes de sa variabilité. Nous avons développé un modèle simplifié de couche mélangée océanique, basé sur la formulation "slab mixed layer" (Frankignoul et Hasselmann, 1977). Ce modèle 2D horizontal est implémenté sur l'océan global, avec une résolution proche de 100 km, et intégré au cours d'une année climatologique. Les flux air-mer proviennent du modèle météorologique ECMWF (réanalyse ERA40) et les courants de surface sont issus de l'altimétrie (analyse SSALTO-DUACS). La profondeur de la couche mélangée (MLD) est dérivée des observations de SST, grâce à une technique d'inversion originale. La MLD obtenue par inversion est bien corrélée aux estimations basées sur des données in situ. Cette profondeur effective représente la pénétration des flux air-mer et assure la cohérence entre les flux, les courants et la SST. La simulation est d'abord validée en examinant le bilan de chaleur dans l'Atlantique Nord-Est, par comparaison aux mesures et aux modèles de l'expérience POMME. Puis le bilan de salinité est étudié dans le domaine global, en termes de distribution géographique et d'évolution saisonnière. L'équilibre entre les différents processus est généralement plus complexe que pour la température. Notamment, le rôle du flux atmosphérique est moins prépondérant (22%), tandis que l'advection géostrophique (33%) et le mélange diapycnal (22%) contribuent fortement. Nos résultats montrent que ce modèle parvient à restituer la variabilité de la SSS sur la majeure partie des océans. Le modèle simule aussi les variations journalières de SSS, qui ne sont pas représentées à l'échelle globale dans les observations actuelles. Grâce à sa simplicité et à sa rapidité, le modèle pourra être utile dans le cadre de SMOS. Il pourra aider à la calibration/validation de la mesure et fournir une estimation a priori pour l'algorithme complexe nécessaire à la restitution de la SSS.

Droits : info:eu-repo/semantics/openAccess
http://archimer.ifremer.fr/doc/2006/these-2302.pdf
http://archimer.ifremer.fr/doc/00000/2302/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
To evaluate the performance of downhole and surface geophysical monitoring methods, a series of shallow gas injection-monitoring experiments has been performed in a coastal saline aquifer at Maguelone, France. The recorded data include pressure measurements with a Westbay multilevel completion and CO2 saturation at an observation well derived from electrical resistivity with a modified Waxman-Smits (MWS) model. In this work, the aim is to develop a simulation model capturing the gas transport behavior and consistent with field data. For this purpose, the simulation of the CO2 injection experiment is carried out with two conceptual models, a homogeneous model and a heterogeneous model treated with multiple realization Monte Carlo simulations. Numerical simulator TOUGH2 with the equation of state module EOS7C is used for the simulations. Comparison of the model results with field data suggests that the pressure responses are captured with relatively good accuracy. Similarly, the model also provides an overall reasonable agreement and correct order of magnitude for predicted gas saturation values. However, as the heterogeneity pattern in the field data remains largely unknown, the model predictions can only be used to capture the mean behavior as well as to provide insights into how heterogeneity can influence the system behavior, by means of sensitivity analyses of the influence of heterogeneities on individual realizations.
ISSN: 1750-5836

hal-01467240
https://hal.archives-ouvertes.fr/hal-01467240
DOI : 10.1016/j.ijggc.2016.09.006

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
The way faults control upward fluid flow in nonmagmatic hydrothermal systems in extensional context is still unclear. In the Eastern Pyrénées, an alignment of twenty-nine hot springs (29 ∘ C to 73 ∘ C), along the normal Têt fault, offers the opportunity to study this process. Using an integrated multiscale geological approach including mapping, remote sensing, and macro-and microscopic analyses of fault zones, we show that emergence is always located in crystalline rocks at gneiss-metasediments contacts, mostly in the Têt fault footwall. The hot springs distribution is related to high topographic reliefs, which are associated with fault throw and segmentation. In more detail, emergence localizes either (1) in brittle fault damage zones at the intersection between the Têt fault and subsidiary faults or (2) in ductile faults where dissolution cavities are observed along foliations, allowing juxtaposition of metasediments. Using these observations and 2D simple numerical simulation, we propose a hydrogeological model of upward hydrothermal flow. Meteoric fluids, infiltrated at high elevation in the fault footwall relief, get warmer at depth because of the geothermal gradient. Topography-related hydraulic gradient and buoyancy forces cause hot fluid rise along permeability anisotropies associated with lithological juxtapositions, fracture, and fault zone compositions.
ISSN: 1468-8115

hal-01582968
https://hal-brgm.archives-ouvertes.fr/hal-01582968
https://hal-brgm.archives-ouvertes.fr/hal-01582968/document
https://hal-brgm.archives-ouvertes.fr/hal-01582968/file/Taillefer_et_al_Geofluids_2017.pdf
DOI : 10.1155/2017/8190109

Résumé : Cette thèse à pour objet l'étude théorique et numérique de solutions dans les équations aux dérivées partielles non linéaires de la physique, en particulier en dynamique des fluides. La présence de discontinuités dans les solutions de ces équations complique la compréhension mathématique des phénomènes mis enjeu et leur traitement numérique, notamment en vue de simulations informatiques . Nous étudions ces équations par une méthode de régularisation dans un espace fonctionnel approprié. Lorsque des schémas numériques construits par des méthodes différentes conduisent à des résultats identiques, ceci jusque dans leurs moindres détails, il semble alors naturel de s'interroger dans quelle mesure ces suites de solutions numériques constituent une approximation d'une solution des équations étudiées. Nous construisons des suites de solutions approchées à partir d'un schéma numérique original,stable et suffisamment simple pour démontrer que ses suites constituent une méthode asymptotique de Maslov au sens des distributions en dimension trois d'espèce. La technique de régularisation employée consiste à étendre les variables réelles du problème ne des variables complexes, ce qui nous permet de construire des familles de solutions particulières que l'on ramène au cas réel en faisant tendre un petit paramètre vers O. Les solutions physiques recherchées apparaissent alors comme valeurs au bord de fonction holomorphes. Nous illustrons les résultats obtenus par des applications en cosmologie dans les cadres Newtoniens et relativistes pour des systèmes sans pression, puis avec pression et auto-gravitation, ainsi que pour le système des gaz parfaits.
This thesis is devoted to the theoretical and numerical study of singular solutions appearing in nonlinear partial differential complicates the mathematical understanding of the phenomena under concem as well as their numerical treatment, in particular in view of computation. These equations are studied by a regularization method in an appropriate functional space. When completely different numerical methods give the same results up to the smallest details one can reasonably expect that these numerical results suggest the existence of a mathematical solution of theses equations. We construct sequences of approximate solutions from an original numerical scheme, which is stable and simple enough to prove that these sequences constitute a Maslov asymptotic method in three space dimension. The regularization technique in use consits in extending the real variables of the problem into complex ones, which perrnits to construct families of particular equations that we bring back to the real case by letting a small paramater tend to zero. The expected physical solutions appear as boundary values of holomorphie functions . Illustrations are given by applications to cosmology in the Newtorian and re1ativistic settings for pressure1ess fluid dynamics, then in presence of self-gravitation and pressure as weil as for the systemof ideal gases

http://www.theses.fr/2011AGUY0478/document

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This paper aims to present an overview on the influence of rheological heterogeneity and mechanical anisotropy on the deformation of continents. After briefly recapping the concept of rheological stratification of the lithosphere, we discuss two specific issues: (1) as supported by a growing body of geophysical and geological observations, crust=mantle mechanical coupling is usually efficient, especially beneath major transcurrent faults which probably crosscut the lithosphere and root within the sublithospheric mantle; and (2) in most geodynamic environments, mechanical properties of the mantle govern the tectonic behaviour of the lithosphere. Lateral rheological heterogeneity of the continental lithosphere may result from various sources, with variations in geothermal gradient being the principal one. The oldest domains of continents, the cratonic nuclei, are characterized by a relatively cold, thick, and consequently stiff lithosphere. On the other hand, rifting may also modify the thermal structure of the lithosphere. Depending on the relative stretching of the crust and upper mantle, a stiff or a weak heterogeneity may develop. Observations from rift domains suggest that rifting usually results in a larger thinning of the lithospheric mantle than of the crust, and therefore tends to generate a weak heterogeneity. Numerical models show that during continental collision, the presence of both stiff and weak rheological heterogeneities significantly influences the large-scale deformation of the continental lithosphere. They especially favour the development of lithospheric-scale strike-slip faults, which allow strain to be transferred between the heterogeneities. An heterogeneous strain partition occurs: cratons largely escape deformation, and strain tends to localize within or at the boundary of the rift basins provided compressional deformation starts before the thermal heterogeneity induced by rifting are compensated. Seismic and electrical conductivity anisotropies consistently point towards the existence of a coherent fabric in the lithospheric mantle beneath continental domains. Analysis of naturally deformed peridotites, experimental deformations and numerical simulations suggest that this fabric is developed during orogenic events and subsequently frozen in the lithospheric mantle. Because the mechanical properties of single-crystal olivine are anisotropic, i.e. dependent on the orientation of the applied forces relative to the dominant slip systems, a pervasive fabric frozen in the mantle may induce a significant mechanical anisotropy of the whole lithospheric mantle. It is suggested that this mechanical anisotropy is the source of the so-called tectonic inheritance, i.e. the systematic reactivation of ancient tectonic directions; it may especially explain preferential rift propagation and continental break-up along pre-existing orogenic belts. Thus, the deformation of continents during orogenic events results from a trade-off between tectonic forces applied at plate boundaries, plate geometry, and the intrinsic properties (rheological heterogeneity and mechanical anisotropy) of the continental plates.
ISSN: 0040-1951

hal-01389719
http://hal.univ-reunion.fr/hal-01389719
http://hal.univ-reunion.fr/hal-01389719/document
http://hal.univ-reunion.fr/hal-01389719/file/vauchez_etal_rheology_1998_hal.pdf
DOI : 10.1016/S0040-1951(98)00137-1

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
We propose a new thermomechanical numerical model of mid-ocean ridge accretion aimed at investigating asymmetric spreading rates, diverse configurations of lens and sill magma injections, crystallization and depth, and on- and off-axis patterns of hydrothermal cooling. The numerical algorithm iteratively resolves temperature and motion equations until it reaches a stationary solution. The motion equation was written in a vorticity-stream function formalism, with boundary and internal conditions applied to the stream function to impose the style of magma injection. Unlike in previous models, our model does not assume an a priori shape for the temperature field, which is initiated by an initial half-space cooling according to the left and right spreading rates. Complex patterns of hydrothermal cooling are simulated by enhanced thermal diffusivity. The model succeeds in describing the dynamic and thermal effects of spreading rates, the style of magma intrusion, and the hydrothermal cooling. Accurate descriptions of these are essential to study the cooling histories of crustal rocks and geophysical observables.
ISSN: 1525-2027

hal-00413118
https://hal.archives-ouvertes.fr/hal-00413118
DOI : 10.1029/2008GC002270

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Experimentally Assisted Molecular Modelling (EAMM) is an original approach for predicting the structure of organic micro-crystals and deducing their habit in the presence of various solvents and additives. It is applied here in the case of n-octylamino-nitrobenzoxadiazole (nOA-NBD), a fluorescent compound. This general approach is first of all described and validated by its blind application to three known and closely-related crystals. Then, the whole process is applied to predict the molecular crystal generated by n-octylamino-NBD, including the molecule conformation, its structure (cell parameters, molecular packing, X-ray powder diagrams) and its theoretical plus actual habits in the presence of a solvent (p-xylene) and two additives (acetic acid and n-dodecane). The conformation of the predicted molecule is validated by comparison with two similar molecules embedded in observable crystals. Then, the choice between two proposed structures (with the same energy and two equiprobable packings) is based on the comparison between experimental and re-computed X-ray powder diffraction diagrams, and also on the interpretation of an actual TEM image in the light of the proposed cell parameters. This choice is confirmed by the comparison between the vapour-grown face development (proposed by the BFDH or the attachment energy models) and the actual face development in a pure solvent. Finally, the actual habit deduced from comparing the crystal attachment energy with the adsorption energy of solvent or additive molecules is validated by SEM images of NBD micro-crystals in presence of the solvent alone or in presence of both additives. In conclusion, the EAMM approach appears to be a decisive tool for quickly simulating structural and habit properties of the molecular crystals, accessible or not by experimentation. Thus, a numerical selection of molecules becomes possible in view of deriving organic materials with predetermined properties, such as the fluorescence properties of the compound considered in the present study.
ISSN: 0022-0248

hal-00857894
https://hal.archives-ouvertes.fr/hal-00857894
DOI : 10.1016/j.jcrysgro.2011.07.035

Éditeur(s) : HAL CCSD Institute for Scientific Computing and Information
Résumé : International audience
In this paper, we solve a problem describing the mechanical deformations of a porous medium in the presence of a monophasic linear flow or a two phase nonlinear flow with the purpose of modelizing subsidence of hydrocarbon reservoirs. An essential characteristics of this problem is that the mechanical deformation and the flow have different time scales. In petroleum industry, one uses different very efficient simulators for the flow problem and the mechanical deformations, which enables to handle complex models. Therefore it is necessary to be able to combine as efficiently as possible the exploitation of these simulators. We propose two alternative splitting approaches. The first one is the staggered algorithm used by engineers, which amounts to a Gauss-Seidel method in the one phase linear case. The second approach is based upon the preconditioned conjugate gradient method. We use a numerical multi-scale method in both of these algorithms. We compare these two approaches and we show that the preconditioned conjugate gradient algorithm is faster and more robust than the staggered algorithm. Applying the preconditioned conjugate gradient algorithm therefore seems to compensate for the fact that the inf-sup condition for the mixed discretization method is not satisfied when combining the simulators for the mechanical deformations and for the flow computations.
ISSN: 1705-5105

Droits : http://creativecommons.org/licenses/by/
hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039/document
https://hal.archives-ouvertes.fr/hal-01449039/file/DHLE.pdf

Éditeur(s) : HAL CCSD Global Science Press
Résumé : International audience
In this paper, the authors present two different time-discretization schemes combined with a finite-volume space discretization for the incompressible Navier-Stokes equations. They show that the collocated space discretization, usually used with a projection method, can also be extended to other time discretizations over previous projection methods, like the splitting methods. The main advantage of this splitting method is that the pressure approximation is truly of second order. In the first and second parts of this article, two different time-discretization schemes are considered with a collocated space discretization and it is explained how the unknown can be correctly coupled. Numerical simulations are presented in the last part. Tests of spatial and time accuracies are made. The numerical results when the Navier-Stokes problem is solved are shown, in the model problem of the driven cavity.
ISSN: 1815-2406

hal-00730191
https://hal.archives-ouvertes.fr/hal-00730191

Éditeur(s) : HAL CCSD EDP Sciences
Résumé : International audience
We show first results of 3D simulation of sea currents in a realistic context. We use the full Navier-Stokes equations for incompressible viscous fluid. The problem is solved using a second order incremental projection method associated with the finite volume of the staggered (MAC) scheme for the spatial discretization. After validation on classical cases, it is used in a numerical simulation of the Pointe à Pitre harbour area. The use of the fictious domain method permits us to take into account the complexity of bathymetric data and allows us to work with regular meshes and thus preserves the efficiency essential for a 3D code.
ISSN: 1270-900X
Guidel, France

hal-00728894
https://hal.archives-ouvertes.fr/hal-00728894
DOI : 10.1051/proc/201235025

Éditeur(s) : HAL CCSD
Résumé : To understand the effect of inherited lithospheric structures on the continental plates deformation process, we have developed a multi scale model that explicitly takes into account an anisotropic viscosity due to the Crystal Preferred Orientation (CPO) of olivine in the mantle by coupling a viscoplastic self-consistent simulation of the deformation at the polycrystalline aggregate to a 3D Finite Element code. Using this tool, we calculate the mechanical behavior of a continental lithosphere submitted to extension for various initial structures. The first set of experiments investigates the deformation of a plate with a homogeneous initial olivine CPO. Results show that the macroscopic anisotropy is strongly dependent on the initial orientation of the olivine CPO relative to the solicitation and to its evolution. The second set of experiments aims to analyze the deformation of a continental plate containing pre-existing lithospheric scale shear zones. These multi-domain models show that the olivine mechanical anisotropy may induce the reactivation of pre-existing lithospheric scale shear zones in a plate submitted to rifting. In both cases, the results shows that the macroscopic mechanic anisotropy induced by the evolution of olivine CPO in the upper mantle is a key parameter to understand the deformation of the continental plates and to explain the observed reactivation of inherited structure from past tectonics events
Afin de comprendre l'effet des structures préexistantes sur la déformation des plaques continentales, nous avons développé une méthode de modélisation mécanique multi-échelle utilisant une viscosité anisotrope induite par les Orientations Préférentielles du Réseau (OPR) des cristaux d'olivine présents dans le manteau. Cette méthode en couple le calcul viscoplastique autocohérent de la déformation d'un agrégat d'olivine à un code éléments finis 3D. Elle a été utilisée pour le calcul du comportement mécanique d'une plaque continentale soumise à une extension. Une première série d'expériences numériques sur l'extension de plaque lithosphérique présentant des OPR initiales d'olivine homogènes a permis de montrer que l'anisotropie mécanique macroscopique est fortement dépendante de l'orientation initiale de l'OPR. La seconde série de simulations analyse l'extension d'une plaque continentale présentant une zone présentant des OPR d'olivine héritées d'un événement tectonique précédant. Les résultats de ces calculs multi-domaines que l'anisotropie mécanique de l'olivine peut expliquer la réactivation des zones de cisaillement héritées lors du rifting. Les résultats obtenus dans les deux séries de modèles montrent que l'anisotropie mécanique macroscopique induite par l'évolution des OPR d'olivine dans le manteau supérieur est un paramètre clé pour la compréhension de la déformation des plaques continentales et pour expliquer la réactivation observée de structures tectoniques, tels que les grands décrochements, lors des épisodes tectoniques ultérieurs.
https://tel.archives-ouvertes.fr/tel-00369463

tel-00369463
https://tel.archives-ouvertes.fr/tel-00369463
https://tel.archives-ouvertes.fr/tel-00369463/document
https://tel.archives-ouvertes.fr/tel-00369463/file/KnollThese09.pdf

Résumé : La thèse est consacrée à la recherche des catastrophes marines (tsunamis, ondes de tempête) dans les Antilles françaises, en utilisant des modèles analytiques et numériques de mécanique des fluides. L'accent est mis sur le développement de modèles de mouvement de glissements de terrain et des vagues causées par ces glissements. Le modèle le plus simple du glissement de terrain «solide block» est utilisé pour évaluer les caractéristiques des flux pyroclastiques du volcan Soufrière-Hills (Montserrat). Un modèle plus complexe de modélisation de glissement de terrain (modèle de Savage-Hutter) a été étudié analytiquement, donnant ainsi une nouvelle famille de solutions exactes décrivant le mouvement de l'écoulement par gravité non linéaire dans une vallée ou dans un canyon sous-marin. Le modèle comprend des ondes simples (Riemann waves),le cas d'un barrage qui cède (dam-break problem), des solutions auto-similaires dans la M-fonne et « chapeaux paraboliques ». Grâce à la théorie linéaire de l'eau peu profonde nous avons étudié le processus de génération de tsunamis par des glissements de terrain de volume variable, se déplaçant à une vitesse variable dans un bassin de profondeur variable. Dans le cas d'un fond marin particulier (cas sans "réflexion"), les phénomènes de résonnance ont été étudiés dans un bassin à profondeur variable. Nous avons utilisé des méthodes numériques pour la résolution non-linéaire des équations des eaux peu profondes afin d'analyser des catastrophes marines réelles (les ondes de tempête causées par le cyclone Lili en 2002, le tsunami volcanique de 2003 à Montserrat) et probables (un tsunami prés des côtes de la Martinique). Des données sur les catastrophes
The dissertation is devoted to research in the field of marine natural hazards (tsunamis, storm surges) in the French West Indies, using analytical and numerical models of fluid mechanics. Emphasis is placed on the development of models of landslide motion and generated tsunami waves. The simple "solid block" model is used to evaluate the characteristics of pyroclastic flow Soufriere Hills volcano (Montserrat). The "fluid model" of a landslide (so called Savage-Hutter model) is studied analytically; within this model a new family of exact solutions that describe the motion of nonlinear gravity flow in a valley or underwater canyon is found: nonlinear Riemann wave, dan break problem, self-similar solutions (M - wave and parabolic cap). In the framework of the linear shallow water theory the process of generation of tsunami waves by landslides of variable volume moving with variable velocity above the basin of variable depth is studied. For the specific bottom profile ("reflectionless" one) the resonant phenomena is investigated in the basin of variable depth. Numerical methods are used to analyze marine hazards: historical (storm surges, caused by Cyclone LILI in 2002; volcanic tsunami 2003 on Montserrat) and possible events (possible tsunami of the coast of Martinique). Various data on marine natural disasters are obtained during field surveys (volcanic tsunami in 2003, stonn surges caused by Hurricane Dean in 2007). Designed catalogs of tsunamis and storm surges are created based on results of numerical modeling and field studies; some statistical analysis is perfomed

http://www.theses.fr/2011AGUY0435/document

Éditeur(s) : HAL CCSD
Résumé : International audience
Proceedings - 2014 International Conference on Electromag- netics in Advanced Applications, ICEAA 2014
Palm Beach, Aruba, Aruba

hal-01116596
https://hal.inria.fr/hal-01116596

Éditeur(s) : HAL CCSD EDP Sciences
Résumé : International audience
In this work, we point out and apply some ideas for the numerical simulation by multi- level schemes of several wave equations including dispersive (Korteweg de Vries) and dissipative ones (Kuramoto-Sivashinski). We start from the approach developed for dissipative equations decomposing the solution into large (Y) and small (Z) scales, and we show that it is necessary to have a regularization in time to allow stable numerical treatment of Y and Z which are here generated using hierarchical processes such as Wavelets, Incremental Unknowns or Fourier.
ISSN: 1270-900X

hal-00768332
https://hal.archives-ouvertes.fr/hal-00768332
https://hal.archives-ouvertes.fr/hal-00768332/document
https://hal.archives-ouvertes.fr/hal-00768332/file/esaim_proc_2009_v27.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Numerical programs for simulating flow and reactive transport in porous media is essential for predicting reservoir properties related to CO2 sequestration performance, subsurface storage and risk assessment. In this paper we solve the Navier Stokes' equation using finite difference method, on a simulated porous rock structure, to study the velocity distribution of fluid flowing through it under a constant pressure gradient. A reactive solute carried through the fluid is allowed to interact with the minerals in the rock. This chemical reaction dissolves the mineral which changes the rock structure thus affecting its flow properties. These changes of flow properties are studied with variation in reactive solute concentration and pressure field. The different mechanisms of dissolution responsible for the variation of flow properties for the different parameters is predicted. Before the onset of homogeneous dissolution, variation in porosity follows a power-law behaviour with change in permeability when the latter is scaled by the concentration of the reactive species. The simulation results are compared with available experimental data and found to give a reasonable match.
ISSN: 0022-1694

hal-00757760
https://hal.archives-ouvertes.fr/hal-00757760
DOI : 10.1016/j.jhydrol.2012.05.024

Éditeur(s) : American Geophysical Union
Résumé : The origins and fates of the Atlantic salinity maximum water (SMW), formed through excess evaporation in the tropics and subtropics of both hemispheres, are studied using monthly mean outputs of a numerical simulation of the world ocean climatological circulation. After defining formation domains from the surface salinity field and the vertical stratification, a Lagrangian technique is used to estimate the formation rates and main pathways in each hemisphere and the role of this water in the framework of the warm water return flow of the meridional overturning cell. Formation rates around 9 and 11 Sv are found in the Southern and Northern Hemispheres, respectively. While the export of the southern SMW from its formation area is realized by the western boundary currents, that of the northern SMW mainly results from interior subduction. Equatorward of the formation regions, a fraction of each SMW variety is entrained in the subtropical cells that connect the subtropics to the equatorial region. Poleward of them, both varieties are seen to feed the regions of subtropical mode water (STMW) formation around 35degrees of latitude in both hemispheres. The bulk of the transport associated with each variety eventually turns northward: This amounts to similar to6 Sv of southern SMW gathered in the North Brazil Undercurrent (NBUC), and similar to10 Sv of northern SMW found in the Gulf Stream at 35degreesN, of which 8 Sv have gone through the Caribbean Sea. Of the 13.4-Sv northward transport of the meridional overturning cell estimated by the model at 47degreesN, more than 50% (6.9 Sv) is found to have transited through at least one of the SMW regions. This gives an indication of the likely important role of SMW formation in the observed northward salinity increase of the upper Atlantic Ocean.
Journal of Geophysical Union - Research C - Oceans (0148-0227) (American Geophysical Union), 2002-09 , Vol. 107 , N. C10 , P. NIL_610-NIL_624

Droits : 2002 by the American Geophysical Union
http://archimer.ifremer.fr/doc/2002/publication-769.pdf
DOI:10.1029/2002JC001318
http://archimer.ifremer.fr/doc/00000/769/

Éditeur(s) : HAL CCSD
Résumé : International audience
In this paper, we compare the effect of rolling resistance at the contacts in granular systems composed of disks with the effect of angularity in granular systems composed of regular polygonal particles. For this purpose, we use contact dynamics simulations. By means of a simple shear numerical device, we investigate the mechanical behavior of these materials in the steady state in terms of shear strength, solid fraction, force and fabric anisotropies, and probability distribution of contact forces. We find that, based on the energy dissipation associated with relative rotation between two particles in contact, the effect of rolling resistance can explicitly be identified with that of the number of sides in a regular polygonal particle. This finding supports the use of rolling resistance as a shape parameter accounting for particle angularity and shows unambiguously that one of the main influencing factors behind the mechanical behavior of granular systems composed of noncircular particles is the partial hindrance of rotations as a result of angular particle shape.
AIP Conference Proceedings
Sydney, Australia

hal-00842799
https://hal.archives-ouvertes.fr/hal-00842799
https://hal.archives-ouvertes.fr/hal-00842799/document
https://hal.archives-ouvertes.fr/hal-00842799/file/Comparison_effects_rolling_Estrada_al.pdf