É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 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 Oxford University Press (OUP)
Résumé : International audience
Estimating how topography is maintained provides insights into the different factors responsible for surface deformations and their relative roles. Here, we develop a new and simple approach to assess the degree of isostatic compensation of continental topography at regional scale from GOCE gravity gradients. We calculate the ratio between the radial gradient observed by GOCE and that calculated from topography only. From analytical and statistical formulations, simple relationships between this ratio and the degree of compensation are obtained under the Airy–Heiskanen isostasy hypothesis. Then, a value of degree of compensation at each point of study area can be easily deduced. We apply our method to the Alaska-Canada Cordillera and validate our results by comparison with a standard isostatic gravity anomaly model and additional geophysical information for this area. Both our GOCE-based results and the isostatic anomaly show that Airy–Heiskanen isostasy prevails for the Yukon Plateau whereas additional mechanisms are required to support topography below the Northwest Territories Craton and the Yakutat collision zone.
ISSN: 0956-540X

hal-01467540
https://hal.archives-ouvertes.fr/hal-01467540
https://hal.archives-ouvertes.fr/hal-01467540/document
https://hal.archives-ouvertes.fr/hal-01467540/file/cadioGJI2016.pdf
DOI : 10.1093/gji/ggw281

É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
Résumé : International audience
In this paper, via the algebra of generalized functions, we investigate the generalized Riemann's problem associated to conservation laws with analytical coefficients. This allows us to transform the problem into a system of ordinary differential equations. In some particular cases, such that Burgers' and conservative Richard's equation, approximated solutions are obtained by the truncation of the so called Hugoniot-Maslov's chain and numerical simulations are also presented in the case of equations with polynomial coefficients.
Novi Sad Journal of Mathematics

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

Éditeur(s) : HAL CCSD Society for Industrial and Applied Mathematics
Résumé : International audience
In this paper, we developed a theoretical study for nonconservative sytems in one dimension in order to construct numerical schemes for solving the Riemann problem. The nonconservative form of our model system required the use of a well-adapted theory in order to give us a sense of our problem. We chose a framework of generalized functions for solving a scalar hyperbolic equation with a discontinuous coefficient $sigma_t +usigma_x approx 0$, where u is the velocity solution of a Burgers's equation. After an explicit solution of the Riemann problem, we derived Godunov split schemes for computing an approximate solution of the Cauchy problem. We applied our study to a system modeling elasticity and a system modeling gas dynamics. Some stability properties of a scheme and its convergence to a generalized solution are proved for the first model. Numerical experiments confirmed this convergence result. For the second model, calculations of flows containing weak-to-moderate shocks showed that conservation errors are reduced when the mesh is refined but were not entirely eliminated.
ISSN: 0036-1429

hal-00770266
https://hal.univ-antilles.fr/hal-00770266
DOI : 10.1137/S0036142900378637

Éditeur(s) : HAL CCSD springer
Résumé : International audience
Although optimal control problems for linear systems have been profoundly investigated in the past 50-60 years, the issue of numerical approximations and precise error analyses remains challenging due the bang-bang structure of the optimal controls. Based on a recent paper by M. Quincampoix and V.M. Veliov on metric regularity of the optimality conditions for control problems of linear systems the paper presents new error estimates for the Euler discretization scheme applied to such problems. It turns out that the accuracy of the Euler method depends on the "controllability index" associated with the optimal solution, and a sharp error estimate is given in terms of this index. The result extends and strengthens in several directions some recently published ones.
Lecture notes in computer science
Sozopol, Bulgaria

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

Éditeur(s) : Styli 2003. Trente ans de crevetticulture en Nouvelle-Calédonie. Nouméa - Koné, 2-6 juin 2003. 2004. Gorarant C, Harache Y, Herbland A, Mugnier C. Ed. Ifremer, Actes de Colloque, n°38, pp.113-120
Résumé : Two oyster species are currently present along the French coasts : the indigenous European flat oyster (Ostrea edulis), and the Pacific cupped oyster (Crassostrea gigas), that has been introduced from Japan since the beginning of the 70ies. The flat oyster successively suffered from two protozoan diseases during the 60ies and its production decreased from 20 000 tons/year by that time to 1 500 tons/year nowadays. Consequently, the oyster production is principally (99%) based upon the Pacific oyster species with approximately 150 000 tons/year among which 90% are grown from the natural spat. However, the hatchery production of this species is developing and was estimated to 400 to 800 millions spat in 2002. Moreover, strengthened relationships between IFREMER and the 5 commercial hatcheries, that all joined the SYSAAF (Union of the French poultry, shellfish and fish farming selectors), allow to plan for new genetic breeding programs. At the end of the 80ies, IFREMER initiated a genetic breeding program for the resistance of the European flat oyster to the bonamiosis, and obtained strains more tolerant to this disease. After two generations of massal selection, molecular markers had identified a reduced genetic basis in this program. It was then reoriented to an intra-familial selection. However, we were confronted to a zootechnic problem to manage such a scheme and we compromised by an intra-cohorts of families selection scheme managed using molecular markers. The program has now reached the transfer level with experimentation at a professional scale. Concerning the Pacific cupped oyster, and in parallel with the obtaining and the study of polyploids, performance of different Asian cupped oyster strains were compared to the one introduced in France thirty years ago and currently suffering from summer mortalities. The local strain exhibited better performance, certainly based upon a good local adaptation. In other respects, although early growth is a relevant criteria for selection for growth to commercial stage, it is not to be privileged in the context of an oyster producing region with a limited food availability. Contrary, the spat summer mortality became a priority for numerous teams (genetic, physiology, pathology, ecology,...) joined in the MOREST program. The first results showed important survival differences between fullsib and halsib families. They indicate a genetic determinism to this character "survival" and promote for its selection.
Deux espèces d'huîtres sont aujourd'hui présentes sur les côtes françaises : l'huître plate européenne (Ostrea edulis), indigène, et l'huître creuse du Pacifique (Crassostrea gigas), introduite du Japon dans le début des années 1970. L'huître plate a subi coup sur coup, dans les années 1960, deux maladies dues à des protozoaires parasites et sa production est ainsi passée de 20 000 t/an dans les années 1960 à 1500 t/an aujourd'hui en France. Il en résulte que la filière ostréicole française repose quasiment exclusivement (99%) sur la production de l'huître creuse, C. gigas, avec 150 000 tonnes en 2000 dont environ 90% est issu du captage naturel. Cependant, la reproduction de cette espèce en écloserie est actuellement en expansion avec une estimation de 400 à 800 millions de naissains commercialisés en 2002. Des rapports renforcés entre l'IFREMER et les cinq écloseries françaises de la filière ostréicole, qui adhèrent toutes au SYSAAF (SYndicat des Sélectionneurs Avicoles et Aquacoles Français), permettent désormais d'envisager le développement de nouveaux programmes d'amélioration génétique. A la fin des années 1980, un programme de sélection pour la résistance de l'huître plate à la bonamiose a été initié au sein de l'IFREMER, et a permis d'obtenir des souches plus tolérantes à cette maladie. Après deux générations en sélection massale, les marqueurs moléculaires avaient mis en évidence une base génétique réduite au sein de ce programme qui s'était alors orienté vers une sélection intra-familiale. Devant la lourdeur zootechnique d'un tel schéma, un compromis a été trouvé en réalisant une sélection au sein de cohortes de familles biparentales, la gestion de ces cohortes étant permise par les analyses moléculaires. Le programme est maintenant dans sa phase de transfert, avec des expérimentations à plus grande échelle en partenariat avec la profession. En ce qui concerne l'huître creuse, et en parallèle de l'obtention et l'étude de polyploïdes, les performances de différentes souches d'huîtres creuses provenant d'Asie ont été comparées à celles de la souche introduite en France il y a trente ans, sujette à des mortalités estivales. La souche locale présente de meilleures performances, vraisemblablement grâce à une adaptation bien réalisée. Par ailleurs, bien que l'amélioration de la croissance jusqu'à la taille commerciale s'avère possible par une sélection précoce au stade larvaire, le critère de croissance n'est pas nécessairement à privilégier dans le contexte d'un bassin ostréicole disposant de ressources trophiques limitées. En revanche, les importantes mortalités rencontrées par le naissain en périodes estivales constituent le sujet d'étude prioritaire sur lequel se sont penchées de nombreuses équipes (génétique, physiologie, pathologie, écologie, ...) dans le cadre du grand défi MOREST. Les premiers résultats ont mis en évidence de fortes différences de survie entre des familles biparentales d'animaux indiquant un déterminisme génétique du caractère « survie », et encourageant dans la voie d'une sélection pour ce caractère.

Droits : 2004 Ifremer
http://archimer.ifremer.fr/doc/2003/acte-3491.pdf
http://archimer.ifremer.fr/doc/00000/3491/

Éditeur(s) : HAL CCSD
Résumé : In this work, we develop a variant of a bundle method in order to find a zero of a maximal monotone operator. This algorithm relies on two polyhedral approximations of the epsilon-enlargement of the considered operator, via a systematic use of the transportation formula. Moreover, the use of a double polyhedral approximation in our algorithm could inspire other bundle methods for the case where the given operator can be split as the sum of two other maximal monotone operators.
https://hal.archives-ouvertes.fr/hal-00828930

hal-00828930
https://hal.archives-ouvertes.fr/hal-00828930
ARXIV : 1305.5810

Éditeur(s) : HAL CCSD Elsevier
Résumé : The Gulf of Aigues-Mortes (NW Mediterranean Sea) is a midshelf zone whose scale is an intermediate between the nearshore scale (0–10 m depth) and the coastal scale (including the whole continental shelf). Its hydrodynamics is investigated for the first time. ADCP, CTD and thermosalinograph data were collected during three short cruises (HYGAM; March 6–7, 20–21, April 5–6, 2005). They were scheduled approximately every 15 days to sample the gulf circulation under different weather conditions. Moreover, the cruise data were used to validate the Symphonie model, a 3D primitive equations circulation model. The circulation features displayed by in situ data were well reproduced by Symphonie. A downscaling modelling approach was implemented, the largest scale being obtained by the replay of the MFSTEP regional model of the North-Western Mediterranean Sea. The analysis, closely coupling in situ measurements and model results, provides information that would not have been obtained using data separately. The great variability of the oceanic circulation at this scale is well evidenced. Winds are the main forces, not only locally but also at a larger scale. North winds generate coherent structures interacting with the general Mediterranean circulation. The explanation of the induced currents is then not straightforward, some of them being, for example, northward. South and southeast winds reinforce the sloping surface, the latter allowing the geostrophic equilibrium with longshore currents. This study, focused on the end of winter 2004–2005, also enhances the heat fluxes from offshore to nearshore, as well as the initiation of shelf dense waters. The knowledge of the midshelf oceanic circulation is then applied to run a numerical experiment which describes both the dispersion of a fictive contaminant from a point S (at the centre of the GAM), and evaluates the risk of contamination of the beaches along the GAM shores.
ISSN: 0278-4343

hal-00196620
https://hal.archives-ouvertes.fr/hal-00196620
DOI : 10.1016/j.csr.2007.06.006

Éditeur(s) : HAL CCSD Springer-Verlag
Résumé : International audience
High-resolution borehole strainmeters are usually installed in tectonically active regions in order to detect slow-slip events, and to estimate slow transients related to earthquake swarms. However, they are also sensitive to other numerous influences, internal or external. Furthermore, the quality of their coupling to the rock through cementation, and the mechanical properties of the rock mass around them, have a critical influence on their records. Many of the existing strainmeters present such problems, and the correction for these effects often remains a challenge. In this paper, we present the analysis of the records of a high-resolution borehole dilatometer (Sacks–Evertson), located in the seismically active rift of Corinth (Greece) (station TRZ in the Trizonia island). We show that the instrument suffers from an only partial-solid coupling, and that the nearby sea tides have a direct (through elastic response) and indirect (through pore-pressure diffusion) effect on the dilatation signal, which adds up to the solid tidal strain source. We propose a methodology that allows, in a first step, to better separate the internal (solid tide) from the external (air pressure, sea level) influences, by calculating a frequency-dependent transfer function outside the range of the tidal periods. We then extrapolate this function, in particular at the tidal periods. In a second step, the resulting variation with frequency of the coupling coefficients with sea level led us to estimate the proportion of instrument not solidly cemented to rock (thus in contact with water pore pressure), which is about 90 % of the total height. Despite the small proportion of solid coupling, the sensor resolution remains very good up to a few tens of hours of a time period, thanks to the confining effects of the rocks on the local pore pressure. These results allow us to correct for the external effects, and reduce the associated variance by 80–90 % (in the period range of minutes to days). The empirical correction of the sea level effect could be explained using a simple Boussinesq’s approximation and 1D pore-pressure diffusion model, which contributed to better constraint of some of the poro-elastic parameters in the vicinity of the instrument. After correction, the solid tidal signal at the 24-h period is almost anti-correlated with those of the theoretical solid tide. This surprising result is consistent with a similar anti-correlation observed for the longest period surface waves (200 s) comparing the TRZ dilatometer signals to the strain measured by a nearby borehole strainmeter (MOK, 15 km). This could be related to the presence of a shallow fault close to the instrument, which would creep in response to seismic wave-related stress.
ISSN: 0033-4533

hal-01122853
https://hal.archives-ouvertes.fr/hal-01122853
DOI : 10.1007/s00024-013-0742-2

Éditeur(s) : HAL CCSD
Résumé : International audience
Extreme rainfall events can trigger numerous landslides in mountainous areas and a prolonged increase of river sediment load. The resulting mass transfer at the Earth surface in turn induces stress changes at depth, which could be sufficient to trigger shallow earthquakes. The 2009 Morakot typhoon represents a good case study as it delivered 3 m of precipitation in 3 days and caused some of the most intense erosion ever recorded. Analysis of seismicity time-series before and after the Morakot typhoon reveals a systematic increase of shallow (i.e. 0-15 km of depth) earthquake frequency in the vicinity of the areas displaying a high spatial density of landslides. This step-like increase in frequency lasts for at least 2-3 years and does not follow an Omori-type aftershock sequence. Rather, it is associated to a step change of the Gutenberg-Richter b-value of the earthquake catalog. Both changes occurred in mountainous areas of southwest Taiwan, where typhoon Morakot caused extensive landsliding. These spatial and temporal correlations strongly suggest a causal relationship between the Morakot-triggered landslides and the increase of earthquake frequency and their associated b-value. We propose that the progressive removal of landslide materials from the steep mountain landscape by river sediment transport acts as an approximately constant increase of the stress rate with respect to pre-typhoon conditions, and that this in turn causes a step-wise increase in earthquake frequency. To test this hypothesis, we investigate the response of a rate-and-state fault to stress changes using a 2-D continuum elasto-dynamic model. Consistent with the results of Ader et al. (2013), our preliminary results show a step-like increase of earthquake frequency in response to a step-like decrease of the fault normal stress. We also investigate the sensitivity of the amplitude and time-scale of the earthquake frequency increase to the amplitude of the normal stress change and to rheological parameters. Our study offers new insights on the potential influence of extreme erosional events on the short-time scale dynamics of faults and earthquakes.
American Geophysical Union Fall Meeting 2016
San Francisco, United States

insu-01416987
https://hal-insu.archives-ouvertes.fr/insu-01416987