Éditeur(s) : American Geophysical Union
Résumé : As part of the French CLIPPER project, an eddy permitting model of the Atlantic circulation has been run for 22 years. The domain has open boundaries at Drake passage and at 30 degreesE, from Africa to Antarctica. The simulated mean circulation, as well as the eddy activity, is satisfactory for a 1/3 degrees model resolution, and the meridional heat transport at 30 degreesS is within the range estimated from observations. We use the "mixed" open boundary algorithm of Barnier et al. [1998], which has both a radiation condition and a relaxation to climatology. The climatological boundary forcing strongly constrains the solution in the whole domain. The model heat balance adjusts through the surface (heat flux retroaction term) more than the open boundaries. The radiation phase velocities calculated within the algorithm are analyzed. This shows, quite surprisingly, that both the eastern and western boundaries have a similar behavior, regardless of the preferred directions for advection (mainly eastward) and wave propagation (mainly westward). Our results confirm that open boundary algorithms behave differently according to the dynamics of the region considered. The passive boundary condition that Penduff et al. [2000] applied successfully in the north eastern Atlantic does not work in the present South Atlantic model. We emphasize the need for a careful prescription of the climatology at the open boundary, for which a new approach based on synoptic sections is implemented.
Journal of Geophysical Union - Research C - Oceans (0148-0227) (American Geophysical Union), 2001-10 , Vol. 106 , N. C10 , P. 22115-22129

Droits : 2001 American Geophysical Union
http://archimer.ifremer.fr/doc/2001/publication-454.pdf
DOI:10.1029/2000JC000376
http://archimer.ifremer.fr/doc/00000/454/

Éditeur(s) : HAL CCSD
Résumé : The present PhD thesis is devoted to the numerical simulation and analysis of shear on cohesive granular materials. In the framework of the Contact Dynamics method, we developed a model of cohesive bonding, which mimics cohesion by cementation. This model involves three contact strength parameters (a tensile strength, a coefficient of sliding friction, and a coefficient of rolling friction) as well as a local criterion, which controls the creation and the breaking of cohesive bonds. All contacts are cohesive at the beginning of the numerical simulations, and contact rupture leads to an irreversible loss of tensile strength (debonding). We varied systematically the three contact strength parameters and studied their influence on the shear strength, the dilatancy, and the microstructure of the granular media. Our results show that internal friction and macroscopic tensile strength have different origins at the contact scale. We also find that the microstructure of the granular media, when sheared at the steady state (also known as the `critical state"), is strongly affected by the local coefficients of sliding and rolling friction. In particular, the mechanical role of the weak network changes with these two contact parameters, which results from the formation of a new structure of force transmission. On the other hand, we also applied the same numerical simulation method to the simulation of landslides. We developed a methodology, which takes into account the initial morphology of the hillslope, the strength of geomaterials, and a triggering factor, integrating in the same model the initiation and propagation phases of the landslide. To illustrate this methodology, we applied it to the simulation of a hypothetical example of avalanche triggered by the progressive rise of a water table.
Ce travail de thèse porte sur la modélisation numérique et l'analyse du cisaillement dans les milieux granulaires cohésifs. Dans le cadre de la méthode de la Dynamique des Contacts, nous avons développé un modèle de cohésion qui cherche à imiter la cohésion par cimentation. Ce modèle fait intervenir trois paramètres de résistance au niveau des contacts (une résistance à la traction, un coefficient de frottement de glissement, et un coefficient de frottement de roulement) ainsi qu'un critère de création et de rupture des liens cohésifs. Les contacts sont cohésifs au début des simulations et leur rupture entraîne la perte irréversible de la résistance à la traction. Nous avons fait varier les trois paramètres microscopiques et nous avons étudié leur influence dans la résistance au cisaillement, la dilatance, et la microstructure du milieu. Nos résultats montrent que le frottement interne et l'adhésion macroscopiques ont des origines microscopiques différentes. Nous observons aussi que la microstructure du milieu cisaillé à l'état stationnaire (aussi connu comme ``l'état critique") est fortement affectée par l'angle de frottement et la résistance au roulement au niveau des contacts. En particulier, le rôle mécanique du réseau faible change avec ces deux paramètres microscopiques, ce qui résulte de la formation d'une nouvelle structure de transmission de forces. D'autre part, nous avons appliqué cette méthode de simulation numérique à l'étude d'un processus naturel : les glissements de terrain. Nous avons développé une méthodologie qui tient compte de la morphologie initiale du système, de la résistance des géomatériaux, et d'un mécanisme de déclanchement, ce qui permet d'intégrer dans le même modèle les phases d'initiation et de propagation du glissement. Pour illustrer cette méthodologie, nous l'avons appliquée à la modélisation d'un exemple hypothétique d'avalanche déclenchée par la montée progressive du niveau d'une nappe d'eau.
https://tel.archives-ouvertes.fr/tel-00311303

tel-00311303
https://tel.archives-ouvertes.fr/tel-00311303
https://tel.archives-ouvertes.fr/tel-00311303/document
https://tel.archives-ouvertes.fr/tel-00311303/file/these_estrada.pdf

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Direct numerical simulations of flow through two millimeter-scale rock samples of limestone and sandstone are performed using three diverse fluid dynamic simulators. The resulting steady-state velocity fields are compared in terms of the associated empirical probability density functions (PDFs) and key statistics of the velocity fields. The pore space geometry of each sample is imaged at 5.06−μm voxel size resolution using X-ray microtomography. The samples offer contrasting characteristics in terms of total connected porosity (about 0.31 for the limestone and 0.07 for the sandstone) and are typical of several applications in hydrogeology and petroleum engineering. The three-dimensional fluid velocity fields within the explicit pore spaces are simulated using ANSYS® FLUENT® ANSYS Inc. (2009), EULAG Prusa et al. (Comput. Fluids 37, 1193–1207 2008), and SSTOKES Sarkar et al. (2002). These computational approaches are highly disperse in terms of algorithmic complexity, differ in terms of their governing equations, the adopted numerical methodologies, the enforcement of internal no-slip boundary conditions at the fluid-solid interface, and the computational mesh structure. As metrics of comparison to probe in a statistical sense the internal similarities/differences across sample populations of velocities obtained through the computational systems, we consider (i) integral quantities, such as the Darcy flux and (ii) main statistical moments of local velocity distributions including local correlations between velocity fields. Comparison of simulation results indicates that mutually consistent estimates of the state of flow are obtained in the analyzed samples of natural pore spaces despite the considerable differences associated with the three computational approaches. We note that in the higher porosity limestone sample, the structures of the velocity fields obtained using ANSYS FLUENT and EULAG are more alike than either compared against the results obtained using SSTOKES. In the low-porosity sample, the structures of the velocity fields obtained by EULAG and SSTOKES are more similar than either is to the fields obtained using ANSYS FLUENT. With respect to macroscopic quantities, ANSYS FLUENT and SSTOKES provide similar results in terms of the average vertical velocity for both of the complex microscale geometries considered, while EULAG tends to render the largest velocity values. The influence of the pore space structure on fluid velocity field characteristics is also discussed.
ISSN: 1420-0597

hal-01174152
https://hal.archives-ouvertes.fr/hal-01174152
DOI : 10.1007/s10596-015-9486-7

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
1] Recent statistical analysis by Lallemand et al. (2008) of subduction zone parameters revealed that the back-arc deformation mode depends on the combination between the subducting (nu(sub)) and upper (nu(up)) plate velocities. No significant strain is recorded in the arc area if plate kinematics verifies nu(up) = 0.5 vsub - 2.3 (cm/a) in the HS3 reference frame. Arc spreading ( shortening) occurs if nu(up) is greater ( lower) than the preceding relationship. We test this statistical law with numerical models of subduction, by applying constant plate velocities far away from the subduction zone. The subducting lithosphere is free to deform at all depths. We quantify the force applied on the two converging plates to sustain constant surface velocities. The simulated rheology combined viscous (non-Newtonian) and brittle behaviors, and depends on water content. The influence of subduction rate vs is first studied for a fixed upper plate. After 950 km of convergence ( steady state slab pull), the transition from extensional to compressive stresses in the upper plate occurs for vs similar to 1.4 cm/a. The effect of upper plate velocity is then tested at constant subduction rate. Upper plate retreat ( advance) with respect to the trench increases extension ( compression) in the arc lithosphere and increases ( decreases) the subducting plate dip. Our modeling confirms the statistical kinematic relationship between vsub and nu(up) that describes the transition from extensional to compressive stresses in the arc lithosphere, even if the modeled law is shifted toward higher rates of upper plate retreat, using our set of physical parameters ( e. g., 100 km thick subducting oceanic plate) and short- term simulations. Our results make valid the choice of the HS3 reference frame for assessing plate velocity influence on arc tectonic regime. The subduction model suggests that friction along the interplate contact and the mantle Stokes reaction could be the two main forces competing against slab pull for upper mantle subductions. Besides, our simulations show that the arc deformation mode is strongly time dependent.
ISSN: 1525-2027

hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328/document
https://hal.archives-ouvertes.fr/hal-00412328/file/arcayG32008.pdf
DOI : 10.1029/2007GC001875

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Previous work has demonstrated that geoelectrical measurements, acquired either along the Earth’s surface or in boreholes, can be sensitive to the presence of fractures. However, a lack of numerical approaches that are well suited to modelling electric current flow in fractured media prevents us from systematically exploring the links between geoelectrical measurements and fractured rock properties. To address this issue, we present a highly computationally efficient methodology for the numerical simulation of geoelectrical data in 2.5-D in complex fractured domains. Our approach is based upon a discrete-dual-porosity formulation, whereby the fractures and rock matrix are treated separately and coupled through the exchange of electric current between them. We first validate our methodology against standard analytical and finite-element solutions. Subsequent use of the approach to simulate geoelectrical data for a variety of different fracture configurations demonstrates the sensitivity of these data to important parameters such as the fracture density, depth, and orientation.
ISSN: 0956-540X

hal-01622151
https://hal.archives-ouvertes.fr/hal-01622151
https://hal.archives-ouvertes.fr/hal-01622151/document
https://hal.archives-ouvertes.fr/hal-01622151/file/sanzGJI2017.pdf
DOI : 10.1093/gji/ggx080

Résumé : Ce travail est l'expression d'une volonté de chercheurs, de la Carabe d'améliorer les connaissances scientifiques méso- et micro-météorologiques appliquées aux milieux insulaires sous influence des alizés et de développer la recherche dans ces domaines.On sait que le phénomène météorologique le plus remarquable impactant les îles de la Caraïbe reste le cyclone tropical. Mais d'autres phénomènes, a des échelles inferieures, telles que les pluies intenses, les houles, la dégradation de la qualité de l'air ont une importance aigüe en termes de risques naturels ou de risques sur la sante. Ces exemples attestent la nécessité d'utiliser des méthodes de descente d'échelle pour exploiter l'information météorologique et/ou climatique de grande échelle et dériver des scenarios locaux et régionaux appliques aux territoires insulaires. Ce défi est important car l'attente d'analyses scientifiques pertinentes dans ces domaines est grande.Les travaux présentés dans ce mémoire ont pour principaux objectifs la simulation numérique puis l'analyse des mécanismes méso- et micro-échelles qui induisent des circulations locales diurnes et nocturnes sur l'archipel de la Guadeloupe à l'aide de codes numériques météorologiques éprouves car largement utilises en recherche et en prévision opérationnelle.Ils constituent donc la première étude de modélisation numérique à haute résolution en basse atmosphère, par descente d'échelle dynamique, pour des intervalles d'espace compris entre 1 km et 111 m sur cet archipel.Le modelé atmosphérique Weather Research and Forecasting ARW V3.4 (WRF ARW), non-hydrostatique, a été utilise pour l'ensemble des simulations pour modéliser la troposphère depuis l'échelle globale à l'échelle de la turbulence.Trois situations météorologiques classiques réelles d'une durée de 48 heures, correspondant à 80 % des situations météorologiques observées dans la zone, alizes soutenus (AS), alizes moyens (AM), alizes faibles (AF) ont été examinées. Ces situations sont caractérisées par les valeurs du nombre de Froude local suivantes : 0,82 (AS), 0,41 (AM) et 0,21 (AF). Six domaines de modélisation ont été sélectionnés pour effectuer les descentes d'échelle dynamiques : D01 (maille de 27 km), D02 (maille de 9 km), D03 (maille de 3 km), D04 (maille de 1 km), D05 (maille de 333 m) et D06 (maille de 111 m) avec soixante-dix niveaux verticaux. Les quatre premiers domaines (D01 à D04) couvrent l'archipel de la Guadeloupe et sont utilisés en mode méso-échelle à l'aide d'un schéma d'ensemble de couche limite planétaire YSU. Les domaines D05 (couvrant l'île de la Basse-Terre et le centre de l'archipel) et D06 (couvrant la zone littorale et rurale du Petit Cul-de-Sac Marin et la zone urbaine de l'agglomération pointoise) sont utilisés en mode Large Eddy Simulation avec une fermeture de la turbulence 1,5 TKE 3D. Le modèle WRF a été forcé toutes les six heures par l'assimilation des champs d'analyses globales du modèle NCEP FNL (1° de résolution). Les simulations effectuées ont permis d'obtenir des champs de variables météorologiques 10-minutes à très hautes résolutions spatiales.Les résultats des simulations méso- et micro-échelles ont été confrontés aux valeurs expérimentales obtenues à l'aide de capteurs places sur des mâts météorologiques (campagne Atmo-Mangrov et réseau de mesure Météo-France). Il s'agissait d'optimiser l'utilisation des couplages de codes numériques tout en conservant la possibilité de les confronter aux observations expérimentales.Les résultats des simulations numériques micro-échelles, des différents cas (AS, AM, AF) sont utilisées pour forcer (c'est-à-dire définir les conditions limites) un modèle lagrangien de dispersion de particules : FLEXPART. Le système couple FLEXPART-WRF a été employé dans le domaine D06 pour étudier la dispersion du panache d'oxydes d'azote émis par la principale usine de production d'électricité de l'île.
The present work expresses the will of Caribbean researchers to improve the meso- and micro-meteorological scientific knowledge of the trade winds influenced island areas, and to develop research in these domains.It is well known that tropical hurricanes remain the most remarkable meteorological phenomena that affect the Caribbean islands. But some other phenomena, of smaller scale, such as intense rainfall events, swells, or air quality degradation, are of extreme importance for natural or health hazards. These examples show the need to use downscaling methods to exploit large scale meteorological or climatic information, and to obtain local and regional scenarios for the island areas. This is an important challenge, as sound scientific studies in these matters are eagerly expected. The aim of the research works exposed in the present dissertation is numerical simulation and analysis of the meso- and microscale mechanisms that induce diurnal and nocturnal local circulations in the Guadeloupe archipelago, using numerical meteorological models that are widely used in research and in operational forecasting.These works represent the first high-resolution (1 km to 111 m) numerical study of the lower atmosphere over the Guadeloupe archipelago. The Weather Research and Forecasting ARW 3.4 (WRF-ARW) model is used to simulate the troposphere from global scale to turbulence scale. Real cases of three typical weather types (80% of cases during a year) are examined during 48 hours : strong trade winds (STW), medium trade winds (MTW) and weak trade winds (WTW). These weather types are characterized by values of the local Froude number : 0.82 (STW), 0.41 (MTW) and 0.21 (WTW). Six domains have been selected for the dynamical downscaling : D01 (grid spacing of 27 km), D02 (grid spacing of 9 km), D03 (grid spacing of 3 km), D04 (grid spacing of 1 km), D05 (grid spacing of 333 m) and D06 (grid spacing of 111 m), including 70 vertical levels. The first four domains (D01 to D04) cover the Guadeloupe archipelago and are used in the meso-scale simulations with the planetary boundary layer scheme YSU (ensemble mean). Domain D05 (covering the Basse-Terre island and the middle of the archipelago) and domain D06 (covering the coastal and rural area of Le Petit Cul-de-Sac Marin and the urban area of Pointe-à-Pitre), are employed in the micro-scale simulation (LES) with the 3D TKE 1.5 order closure scheme. WRF has been 6 hourly reinitialized with the NCEP FNL global analyses (resolution of 1°). These simulations permitted to obtain 10-minutes meteorological variable fields with a very high resolution (111 m).Meso-scale and micro-scale model results have been evaluated with observational data from meteorological stations (field campaign Atmo-Mangrov, French Met Office).Once validated, the micro-scale model outputs have been used for the assimilation of the lagrangian particle dispersion model : FLEXPART. The coupling FLEXPART-WRF has been employed in domain D06, in order to analyze the dispersion of the nitrogen oxide plume emitted by the main power plant of the archipelago.

http://www.theses.fr/2014AGUY0735/document

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The full numerical simulation of turbulent flows in the context of industrial applications remains a very challenging problem. One of the difficulties is the presence of a large number of interacting scales of different orders of magnitude, ranging from the macroscopic scale to the Kolmogorov dissipation scale. In order to better understand and simulate these interactions, new algorithms of incremental type have recently been introduced, stemming from the theory of infinite-dimensional dynamical systems (see, e.g., the algorithms of Foias, Jolly et al. (1988), Foias, Manley and Temam (1988), Laminie et al. (1993, 1994), Marion and Temam (1989, 1990), and Temam (1990)). These algorithms are based on decompositions of the unknown functions into a large- and a small-scale component, one of the underlying ideas being to approximate the corresponding attractor. In the context of finite element methods, the decomposition of a solution into small- and large-scale components appears when we consider hierarchical bases (Yserentant, 1986; Zienkiewicz et al., 1982). In the present article we derive several estimates of the size of the structures for the linear and nonlinear terms which correspond to interactions of different hierarchical components of the velocity field, and also their time variation. The one-step time variation of these terms can be smaller than the expected accuracy of the computation. Using this remark, we implement an adaptive spatial and temporal multilevel algorithm which treats differently the small- and large-scale components of the flow. We derive several a priori estimates in order to study the perturbation introduced into the approximated equations. All the interaction terms between small- and large-scale structures are frozen during several time steps. Finally, we implement the multilevel method in order to simulate a bidimensional flow described by the Burgers equations. We perform a parametric study of the procedure and its efficiency. The gain in CPU time is significant and the trajectories computed by our multi-scale method remain close to the trajectories obtained with the classical Galerkin method.
ISSN: 0168-9274

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

Éditeur(s) : HAL CCSD
Résumé : International audience
Structure and evolution of networks have been areas of growing interest in recent years, especially with the emergence of Social Network Analysis (SNA) and its application in numerous fields. Researches on diffusion are focusing on network modeling for studying spreading phenomena. While the impact of network properties on spreading is now widely studied, involvement of network dynamicity is very little known. In this paper, we address the epidemiology context and study the consequences of network evolutions on spread of diseases. Experiments are conducted by comparing incidence curves obtained by evolution strategies applied on two generated and two real networks. Results are then analyzed by investigating network properties and discussed in order to explain how network evolution influences the spread. We present the MIDEN framework, an approach to measure impact of basic changes in network structure, and DynSpread, a 2D simulation tool designed to replay infections scenarios on evolving networks.
Database and Expert Systems Applications - 22nd International Conference, DEXA 2011, Toulouse, France, August 29 - September 2, 2011, Proceedings, Part II
Toulouse, France

hal-00614925
https://hal.archives-ouvertes.fr/hal-00614925
DOI : 10.1007/978-3-642-23091-2_49

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Numerical models show that small-scale convection (SSC) occurring atop a mantle plume is a plausible mechanism to rejuvenate the lithosphere. The triggering of SSC depends on the density contrast and on the rheology of the unstable layer underlying the stagnant upper part of the thermal boundary layer (TBL). Partial melting may change both properties. We analyze, using 2-D numerical simulations, how partial melting influences the dynamics of time-dependent SSC instabilities and the resulting thermo-mechanical rejuvenation of an oceanic plate moving atop of a plume. Our simulations show a complex behavior, with acceleration, no change, or delay of the SSC onset, due to competing effects of the latent heat of partial melting, which cools the plume material, and of the buoyancy increase associated with both melt retention and depletion of residue following melt extraction. The melt-induced viscosity reduction is too localized to affect significantly SSC dynamics. Faster SSC triggering is promoted for low melting degrees (low plume temperature anomalies, thick lithosphere, or fast moving plates), which limit both the temperature reduction due to latent heat of melting and the accumulation of depleted buoyant residue in the upper part of the unstable layer. In contrast, high partial melting degrees lead to a strong temperate decrease due to latent heat of melting and development of a thick depleted layer within the sublithospheric convecting layer, which delay the development of gravitational instabilities. Despite differences in SSC dynamics, the thinning of the lithosphere is not significantly enhanced relatively to simulations that neglect partial melting.
ISSN: 1525-2027

hal-01277175
https://hal.archives-ouvertes.fr/hal-01277175
DOI : 10.1002/2015GC005967

Éditeur(s) : HAL CCSD
Résumé : Meteorological phenomena with scales lower than hurricanes, from extreme rainfalls to pollution peaks, represent important environmental and health risks in the Caribbean. These examples show the need to downscale large-scale meteorological data to island-scale fields.The main goal of the present numerical study is the analysis of the meso- and micro-scale mechanisms inducing diurnal and nocturnal circulations in the Guadeloupe archipelago. The Weather Research and Forecasting model is employed. The micro-scale results are obtained by the coupling of the 1-km simulated fields and the Large Eddy Simulation method. Then, in order to analyze the dispersion of a pollutants plume, the micro-scale outputs with 111 m of resolution are used to force the lagrangian particles dispersion model FLEXPART. The study has been conducted for three real cases of weather types linked with strong, medium and weak trade winds.The results of simulations, validated with available observed data, permitted to simulate all local circulations, including a nocturnal counter-current flow which may induce a pollution peak in the urban area of Pointe-à-Pitre city.These works, which may be applied to other Lesser Antilles islands, suggest many environmental and energy applications.
Les phénomènes météorologiques d'échelles inférieures aux cyclones, depuis les pluies intenses, jusqu'à la dégradation de la qualité de l'air, ont pour les Caraïbes une importance aiguë en termes de risques naturels ou de risques sur la santé. Ces exemples attestent de la nécessité d'utiliser des méthodes de descentes d'échelle pour exploiter l'information météorologique de grande échelle et dériver des scénarios locaux appliqués aux territoires insulaires. Les travaux présentés ont pour principaux objectifs la simulation numérique et l'analyse des mécanismes méso- et micro-échelles qui induisent des circulations locales diurnes et nocturnes sur l'archipel de la Guadeloupe, à l'aide du modèle numérique météorologique Weather Research and Forecasting. Les résultats des simulations micro-échelles sont obtenus par le couplage des champs modélisés à 1 km avec la méthode Large Eddy Simulation. Les sorties micro-échelles, à 111 m, sont ensuite utilisées pour forcer le modèle lagrangien de dispersion de particules FLEXPART afin d'étudier la dispersion d'un panache de polluants. L'étude a été menée pour trois cas réels de situations météorologiques correspondant à des alizés soutenus, moyens et faibles. Les résultats des simulations, après validation, et comparaison aux données d'observation disponibles, permettent de restituer l'ensemble des circulations locales et soulignent l'existence d'un écoulement nocturne à contre courant susceptible d'induire une forte pollution sur la zone urbaine de Pointe-à-Pitre. Cette étude, transposable aux autres îles des petites Antilles, laisse envisager de nombreuses applications environnementales et énergétiques.
https://hal.archives-ouvertes.fr/tel-01074600

Droits : info:eu-repo/semantics/OpenAccess
tel-01074600
https://hal.archives-ouvertes.fr/tel-01074600
https://hal.archives-ouvertes.fr/tel-01074600/document
https://hal.archives-ouvertes.fr/tel-01074600/file/rcece_thesis.pdf

Éditeur(s) : HAL CCSD
Résumé : Today, one of the main challenges in the oil industry is the exploitation of new resources in naturally fractured reservoirs often located in structurally complex areas such as plate boundaries, mountain ranges or near salt structures. While numerical models of rock deformation based on continuum mechanics are becoming industry standard in providing efficient means for modelling natural fractures in reservoirs, the knowledge of past perturbed stress fields, which is referred as the type (normal, wrench or reverse), orientation and magnitude of the regional or local tectonic stresses that needs to be applied as boundary condition in numerical simulations, is often unknown or approximated.The essence of this thesis is to improve the recovery of the paleostress field needed to constrain numerical simulations. I first described and discuss through a comprehensive mechanically-based multi-parametric study how far we can rely on the commonly used paleostress inversion methods based on fault slip data and the Wallace and Bott assumptions. I then compared this method with a new generation technique based on geomechanics using field observations. Since slip markers on faults are hardly observed in core or image log, I used observed natural fracture data as main drivers for the inversion of the paleostress using geomechanically-based method. I demonstrated, through various outcrop and subsurface examples, how this can efficiently be done. Finally, I aimed at addressing the problem of polyphase fracture data sets with unknown mechanical type (joints, faults, stylolites ...) and expanded the mechanical stress inversion to the separation of tectonic phases.
Aujourd'hui, l'un des principaux défis dans l'industrie pétrolière est l'exploitation de nouvelles ressources dans les réservoirs naturellement fracturés souvent situés dans des zones structurellement complexes tels que les frontières de plaques, les chaînes de montagnes ou les structures de proches des dômes de sel. Alors que les modèles numériques de déformation des roches basée sur la mécanique des milieux continus fournissant des moyens efficaces pour la modélisation de fractures naturelles dans les réservoirs, sont devenus la norme dans l'industrie, les paléo-contraintes régionales ou locales appliquées comme condition aux limites dans les simulations numériques, définies par le régime (normal, décrochant ou inverse), l'orientation, et la magnitude des contraintes tectoniques, sont souvent inconnues ou approximées.L'essence de cette thèse est d'améliorer la prédiction du champ de paléo-contrainte nécessaire pour contraindre les simulations numériques. J’ai décrit dans un premier temps à travers une étude comparative multi-paramétrique complète basée sur la mécanique, dans quelle mesure nous pouvions compter sur les méthodes d'inversion des paléo-contraintes basées sur les hypothèses de Wallace et Bott et les données de glissement sur les failles. J’ai ensuite comparé ces méthodes avec une technique d’inversion de nouvelle génération basée sur la géomécanique en utilisant les observations de terrain. Puisque les marqueurs de glissement sur les failles sont rarement observés le long du forage, j’ai utilisé l’inversion géomécanique des contraintes pour les données de fractures naturelles observées. J’ai alors démontré, l’efficacité de cette technique à travers divers exemples d'affleurements et de subsurface. Enfin, j’ai proposé une solution afin de résoudre le problème des données de fractures polyphasées à type cinématique inconnu (joints, failles, stylolites ...) et étendu l’inversion de contrainte basée sur la géomécanique à la séparation des phases tectoniques.
https://hal.archives-ouvertes.fr/tel-01324110

tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110
https://hal.archives-ouvertes.fr/tel-01324110/document
https://hal.archives-ouvertes.fr/tel-01324110/file/theselejri.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Holocene Lake Mega-Chad (LMC) was the largest late Quaternary water-body in Africa. The development of this giant paleo-lake is related to a northward shift of the isohyetes interpreted as evidence for an enhanced Monsoon (African Humid Period). Numerous preserved coastal features have been described all around the LMC shore. Such features reveal the main paleo-hydrodynamical tendencies. In the context of a closed water-body like LMC, hydrodynamics are forced mainly by winds. We use a three-dimensional numerical model (SYMPHONIE) to simulate the mean hydrodynamics in LMC under both Harmattan-like (northeasterly trade winds) and Monsoon-like (southwesterly winds) forcings. The northern part of LMC displays coastal features, such as sand spits, that are consistent with the simulations forced by Harmattan-like winds. Geomorphic features related to Monsoon-driven hydrodynamics are not clearly expressed. They could have developed during the early stage of LMC but subsequently reworked. At the time of sand-spit building, Harmattan-like driven hydrodynamics prevailed and related coastal features were preferentially preserved in the sedimentary record.
ISSN: 0033-5894

hal-00475633
https://hal.archives-ouvertes.fr/hal-00475633
DOI : 10.1016/j.yqres.2009.10.010

Éditeur(s) : HAL CCSD American Chemical Society
Résumé : Ants have evolved venoms rich in peptides and proteins used for predation, defense, and communication. However, they remain extremely understudied due to the minimal amount of venom secreted by each ant. The present study investigated the differences in the proteome and peptidome of the venom from the bullet ant, Paraponera clavata. Venom samples were collected from a single colony either by manual venom gland dissection or by electrical stimulation and were compared using proteomic methods. Venom proteins were separated by 2D-PAGE and identified by nanoLC-ESI-QTOF MS/MS. Venom peptides were initially separated using C18 reversed-phase high-performance liquid chromatography, then analyzed by MALDI-TOF MS. The proteomic analysis revealed numerous proteins that could be assigned a biological function (total 94), mainly as toxins, or roles in cell regulation and transport. This investigation found that ca. 73% of the proteins were common to venoms collected by the two methods. The peptidomic analysis revealed a large number of peptides (total 309) but with <20% shared by the two collection methods. There was also a marked difference between venoms obtained by venom gland dissection from different ant colonies. These findings demonstrate the rich composition and variability of P. clavata venom.
ISSN: 1535-3893

hal-01603415
https://hal.archives-ouvertes.fr/hal-01603415
DOI : 10.1021/acs.jproteome.6b00948
PRODINRA : 406775
PUBMED : 28118015

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the deformation mechanisms and the seismic properties of 10 eclogite samples from different localities (Alps, Norway, Mali and eastern China) through the analysis of their microstructures and lattice-preferred orientations (LPO). These samples are representative of various types and intensity of deformation under eclogitic metamorphic conditions. Omphacite and garnet LPO were determined from electron backscatter diffraction (EBSD) technique. Garnet appears to be almost randomly oriented whereas omphacite develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the (010)-poles concentrated sub-perpendicular to the foliation. In order to analyze the deformation mechanisms that produced such omphacite LPO, we compare our observations to LPO simulated by viscoplastic self-consistent numerical models. A good fit to the measured LPO is obtained for models in which the dominant slip systems are 1/2h110i{11 ¯ 0}, [001] {110} and [001] (100). Dominant activation of these slip systems is in agreement with TEM studies of naturally deformed omphacite. Seismic properties of eclogite are calculated by combining the measured LPO and the single crystal elastic constants of omphacite and garnet. Although eclogite seismic anisotropies are very weak (less than 3% for both P-and S-wave), they are generally characterized by a maximum P-wave velocity sub-parallel to the lineation and by a minimum velocity approximately normal to foliation. The mean P-and S-wave velocities are high (respectively, 8.6 and 4.9 km/s). The S-wave anisotropy pattern displays complex relationships with the structural frame but the fast polarization plane generally tends to be parallel to the foliation. Calculated reflection coefficients show that an eclogite/crust interface is generally a good reflector (Rc > 0.1), whereas an eclogite body embedded in the upper mantle would be hardly detectable.
ISSN: 0040-1951

hal-01389715
http://hal.univ-reunion.fr/hal-01389715
http://hal.univ-reunion.fr/hal-01389715/document
http://hal.univ-reunion.fr/hal-01389715/file/Bascou_eclo-sismo_tphy2001_hal.pdf
DOI : 10.1016/S0040-1951(01)00156-1

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The outstanding exposure of deformation bands in a multilayer sequence of the Orange quarry, Provence, France, motivates a study investigating the development and controls of deformation band geometries by field mapping combined with finite element modeling. Field mapping yields new insights into angular relations between, and evolution of, the deformation bands as well as provides input parameters for the setup and boundary conditions of the numerical simulations using ADELI 2D. In particular, reverse-sense deformation bands are found to have developed within three of the four exposed layers of the multilayer stack with similar orientations to bedding, indicating layer-parallel contraction. Numerical simulations, carried out as a parametric study to investigate the influence on deformation band development and geometry, successfully reproduce the sequential development of deformation bands in the three layers as seen in the field as well as recreate the observed angular geometries of similar to 39 degrees to the maximum compression. This angular relationship is put into context with observations from other field studies and theoretical values to explore controls of deformation band orientations. Furthermore, band localization is impeded in the model by increased layer cohesive strength, demonstrating the importance of host-rock property variations on where localization of bands within the multilayer sequence. Results from both field investigation and numerical simulations are consistent with regional tectonics, where the strata underwent gentle folding due to thrust faulting at depth, and bring further insight into the development and geometry of deformation bands in both extensional and contractional tectonic regimes.
ISSN: 0148-0227

hal-00639452
https://hal.archives-ouvertes.fr/hal-00639452
DOI : 10.1029/2011JB008365

Éditeur(s) : HAL CCSD Elsevier
Résumé : This paper presents a numerical investigation on the mechanical behavior of cemented granular materials by means of discrete element simulations. In particular, two conceptual elements of this behavior were explored: the yield surfaces and the plastic potentials, as defined in the framework of elastoplastic constitutive models. For this purpose, we implemented a contact model that mimics the mechanical behavior of cemented bonds, we constructed a set of polydisperse samples with different solid fractions, and we tested these samples by means of a biaxial device with four rigid walls. The yield surfaces and plastic potentials were thus obtained and analyzed as functions of the solid fraction of the material, drawing interesting conclusions, confirming previous experimental findings, and, consequently, evidencing the great potential of discrete element simulations in exploring the behavior of real cemented geomaterials.
ISSN: 0266-352X

hal-00853748
https://hal.archives-ouvertes.fr/hal-00853748
DOI : 10.1016/j.compgeo.2012.11.001

Éditeur(s) : Blackwell science
Résumé : This paper deals with the design of modelling tools suitable for investigating the consequences of alternative policies on the dynamics of resources and fisheries, such as the evaluation of marine protected areas (MPA). We first review the numerous models that have been developed for this purpose, and compare them from several standpoints: population modelling, exploitation modelling and management measure modelling. We then present a generic fisheries simulation model, Integration of Spatial Information for FISHeries simulation (ISIS-Fish). This spatially explicit model allows quantitative policy screening for fisheries with mixed-species harvests. It may be used to investigate the effects of combined management scenarios including a variety of policies: total allowable catch (TAC), licenses, gear restrictions, MPA, etc. Fisher's response to management may be accounted for by means of decision rules conditioned on population and exploitation parameters. An application to a simple example illustrates the relevance of this kind of tool for policy screening, particularly in the case of mixed fisheries. Finally, the reviewed models and ISIS-Fish are discussed and confronted in the light of the underlying assumptions and model objectives. In the light of this discussion, we identify desirable features for fisheries simulation models aimed at policy evaluation, and particularly MPA evaluation.
Fish and Fisheries (1467-2960) (Blackwell science), 2005-12 , Vol. 6 , N. 4 , P. 307-349

Droits : 2005 Blackwell Publishing, Inc
http://archimer.ifremer.fr/doc/2005/publication-1031.pdf
DOI:10.1111/j.1467-2979.2005.00199.x
http://archimer.ifremer.fr/doc/00000/1031/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Properties of composite materials depend on interatomic phenomena occurring between binder crystals. Experimental information of Atomic Force Microscopy (A.F.M.) is of prime importance; however understanding is helped by molecular modeling. As underlined in Section 1, the present study is able to simulate crystal interfaces in presence of a solution within apertures less than 1 Nanometer at a full atomic scale. Section 2 presents the case study of a gypsum solution between (120) gypsum faces, with related boundary conditions and atomic interactions. Section 3 deals with the mass equilibrium of the solution within interfaces <5 angstrom, using the original Semi Analytical Stochastic Perturbations (SASP) approach. This information becomes in Section 4 the key for explaining the peak of adhesion obtained in A.F.M. around an aperture of 3 A and gives enlightenments on gypsum cohesion. In conclusion, this illustration shows the potentialities of full atomic modeling which could not be attained by any numerical approach at a mesoscopic scale.
ISSN: 0008-8846

hal-00411868
https://hal.archives-ouvertes.fr/hal-00411868
DOI : 10.1016/j.cemconres.2007.09.022

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : The properties of the subduction interplate domain are likely to affect not only the seismogenic potential of the subduction area but also the overall subduction process, as it influences its viability. Numerical simulations are performed to model the long-term equilibrium state of the subduction interplate when the diving lithosphere interacts with both the overriding plate and the surrounding convective mantle. The thermomechanical model combines a non-Newtonian viscous rheology and a pseudo-brittle rheology. Rock strength here depends on depth, temperature and stress, for both oceanic crust and mantle rocks. I study the evolution through time of, on one hand, the brittle-ductile transition (BDT) depth, zBDT, and, on the other hand, of the kinematic decoupling depth, zdec, simulated along the subduction interplate. The results show that both a high friction and a low ductile strength at the asthenospheric wedge tip shallow zBDT. The influence of the weak material activation energy is of second order but not negligible. zBDT becomes dependent on the ductile strength increase with depth (activation volume) if the BDT occurs at the interplate decoupling depth. Regarding the interplate decoupling depth, it is shallowed (1) significantly if mantle viscosity at asthenospheric wedge tip is low, (2) if the difference in mantle and interplate activation energy is weak, and (3) if the activation volume is increased. Very low friction coefficients and/or low asthenospheric viscosities promote zBDT = zdec. I then present how the subducting lithosphere age affects the brittle-ductile transition depth and the kinematic decoupling depth in this model. Simulations show that a rheological model in which the respective activation energies of mantle and interplate material are too close hinders the mechanical decoupling at the down-dip extent of the interplate, and eventually jams the subduction process during incipient subduction of a young (20-Myr-old) and soft lithosphere under a thick upper plate. Finally, both the BDT depth and the decoupling depth are a function of the subducting plate age, but are not influenced in the same fashion: cool and old subducting plates deepen the BDT but shallow the interplate decoupling depth. Even if BDT and kinematic decoupling are intrinsically related to different mechanisms of deformation, this work shows that they are able to interact closely. Comparison between modelling results and observations suggests a minimum friction coefficient of 0.045 for the interplate plane, even 0.069 in some cases, to model realistic BDT depths. The modelled zdec is a bit deeper than suggested by geophysical observations. Eventually, the better way to improve the adjustment to observations may rely on a moderate to strong asthenosphere viscosity reduction in the metasomatised mantle wedge.
ISSN: 1869-9510

hal-00795559
https://hal.archives-ouvertes.fr/hal-00795559
DOI : 10.5194/se-3-467-2012, 2012

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Tropical islands are characterized by thermal and orographical forcings which may generate microscale air mass circulations. The Lesser Antilles Arc includes small tropical islands (width lower than 50 km) where a total of one-and-a-half million people live. Air quality over this region is affected by anthropogenic and volcanic emissions, or saharan dust. To reduce risks for the population health, the atmospheric dispersion of emitted pollutants must be predicted. In this study, the dispersion of anthropogenic nitrogen oxides (NOx) is numerically modelled over the densely populated area of the Guadeloupe archipelago under weak trade winds, during a typical case of severe pollution. The main goal is to analyze how microscale resolutions affect air pollution in a small tropical island. Three resolutions of domain grid are selected: 1 km, 333 m and 111 m. The Weather Research and Forecasting model (WRF) is used to produce real nested microscale meteorological fields. Then the weather outputs initialize the Lagrangian Particle Dispersion Model (FLEXPART). The forward simulations of a power plant plume showed good ability to reproduce nocturnal peaks recorded by an urban air quality station. The increase in resolution resulted in an improvement of model sensitivity. The nesting to subkilometer grids helped to reduce an overestimation bias mainly because the LES domains better simulate the turbulent motions governing nocturnal flows. For peaks observed at two air quality stations, the backward sensitivity outputs identified realistic sources of NOx in the area. The increase in resolution produced a sharper inverse plume with a more accurate source area. This study showed the first application of the FLEXPART-WRF model to microscale resolutions. Overall, the coupling model WRF-LES-FLEXPART is useful to simulate the pollutant dispersion during a real case of calm wind regime over a complex terrain area. The forward and backward simulation results showed clearly that the subkilometer resolution of 333 m is necessary to reproduce realistic air pollution patterns in this case of short-range transport over a complex terrain area. Globally, this work contributes to enrich the sparsely documented domain of real nested microscale air pollution modelling. This study dealing with the determination of the proper resolution grid and proper turbulence scheme, is of significant interest to the near-source and complex terrain air quality research community.
ISSN: 1352-2310

hal-01349208
https://hal.archives-ouvertes.fr/hal-01349208
DOI : 10.1016/j.atmosenv.2016.05.028