Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Predicting the impact of underground engineering on the environment requires the knowledge of natural media at different scales. In particular, understanding basic phenomena controlling the properties of rocks in the presence of complex fluids necessitates a detailed atomic description of the solid/fluid/solid contacts, the subject of Part I of the present study. First, building the solid interspace between two different crystals in a non-periodic situation is achieved using the ab initio and molecular mechanics code GenMol (TM). A description of the fluid confined within the interspace is then derived from the original genetic iterative multi-species (GIMS) algorithm implemented in the same code. This approach consists of equilibrating chemical potentials, cycle after cycle and species after species, between the confined fluid and the free natural fluid. An elementary iteration for a species k consists of different steps incrementing the number N-k of particles k, with other numbers N-k remaining constant. At each step, an optimum fluid composition is obtained by a genetic process distributing the fluid particles on a grid by stochastic shots, followed in fine by a refining process. The effectiveness of the GIMS approach is demonstrated in the case study of a fluid confined between two (0 0 1) kaolinite faces, with apertures h varying between 4 and 10 A. connected to a 9-species external solution [H2O, Cl-, Na+ CO2(aq), NaCl(aq), Ca2+, Mg2+, HCO3-, H3O+] where concentrations are ranging from 55 to 10(-4) mol/L. The results show a drastic variation in the solute/solvent and cations/ions ratios in the confined fluid when aperture h is lowered to less than 1 nm. These results obtained with a very rapid convergence of the iterative algorithm combined with a very competitive genetic optimizer are validated with high precision on a free solution. This description of contacts between crystals is original and unattainable by standard crystal interface approaches. It opens the way to understanding and/or predicting phenomena between crystals in a complex natural environment, such as adhesion or repulsion investigated in the Part II of this study.
ISSN: 0022-0248

hal-00544569
https://hal.archives-ouvertes.fr/hal-00544569
DOI : 10.1016/j.jcrysgro.2010.08.009

Résumé : Cette thèse est consacrée à la résolution d'un problème fondamental de l'analyse variationnelle qu'est la recherchede zéros d'opérateurs maximaux monotones dans un espace de Hilbert. Nous nous sommes tout d'abord intéressés au cas de l'opérateur somme étendue de deux opérateurs maximaux monotones; la recherche d'un zéro de cet opérateur est un problème dont la bibliographie est peu fournie: nous proposons une version modifiée de l'algorithme d'éclatement forward-backward utilisant à chaque itération, l'epsilon-élargissement d'un opérateur maximal monotone,afin de construire une solution. Nous avons ensuite étudié la convergence d'un nouvel algorithme de faisceaux pour construire ID zéro d'un opérateur maximal monotone quelconque en dimension finie. Cet algorithme fait intervenir une double approximation polyédrale de l'epsilon-élargissement de l'opérateur considéré
This thesis is devoted to solving a basic problem of variational analysis which is the search of zeros of maximal monotone operators in a Hilbert space. First of aIl, we concentrate on the case of the extended som of two maximal monotone operators; the search of a zero of this operator is a problem for which the bibliography is not abondant: we purpose a modified version of the forward-backward splitting algorithm using at each iteration, the epsilon-enlargement of a maximal monotone operator, in order to construet a solution. Secondly, we study the convergence of a new bondie algorithm to construet a zero of an arbitrary maximal monotone operator in a finite dimensional space. In this algorithm, intervenes a double polyhedral approximation of the epsilon-enlargement of the considered operator

http://www.theses.fr/2012AGUY0548/document

Éditeur(s) : HAL CCSD
Résumé : International audience
It is established that the multivalued quasi variational inequalities in uniformly smooth Banach spaces are equivalent to the fixed-point problem. We use this equivalence to suggest and analyze some iterative algorithms for quasi variational inequalities with noncompact sets in Banach spaces. Our results are new and represent a significant improvement of the previously known results.;
JIPAM

hal-00778168
https://hal.univ-antilles.fr/hal-00778168

Éditeur(s) : HAL CCSD Seismological Society of America
Résumé : International audience
We present an improvement of the maximum intersection (MAXI) method in which absolute earthquake location is defined by the maximum intersection number of hyperbolic equal differential time (EDT) volumes (one EDT volume being described as all grid nodes satisfying the arrival‐time differences between two stations, ± a tolerance value known as TERR). This 3D technique is well adapted to a strongly heterogeneous environment, avoids the depth versus origin‐time trade‐off (even using P arrivals alone), and objectively filters possible erroneous arrival times. Improvements consist of introducing an iterative multiscale approach on the TERR parameter and conducting a hierarchy of grid discretization to refine the preliminary solution (multigrid algorithm). The outlier cleaning process is strengthened by conducting statistical analyses based on the number of EDT volume intersections per station, rather than on travel‐time residues. Processing two series of synthetic arrival times computed in the 3D seismological configuration of the southernmost Ryukyu subduction zone, offshore Taiwan, allows us to (1) confirm MAXI robustness, even when arrival times are highly disturbed or azimuthal coverage largely reduced, and (2) show that MAXI, by limiting trade‐offs between origin time and depth and between epicentral position and depth, is less affected by a wrong velocity model than a technique that minimizes travel‐time residues. To conclude, we apply the MAXI method to locate earthquakes recorded by the combined Central Weather Bureau of Taiwan (CWB) and the Japanese Meteorological Agency (JMA) networks (1991-2008) within a georealistic a priori 3D P‐velocity model. Hypocenter determinations are selected based on MAXI confidence factors, and earthquake clustering along well‐known tectonic features witnesses hypocenter determination accuracy.
ISSN: 0037-1106

hal-00725227
https://hal.archives-ouvertes.fr/hal-00725227
DOI : 10.1785/0120100311

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : This paper presents an inversion method using the interseismic velocity field to determine effective rigidity of the lithosphere. The method is based on the minimization of a cost function defined as the quadratic measure of the difference between measured and modelled velocity fields on a discrete set of points. The continuous mapping of the rigidity is fulfilled with a limited set of parameters and the forward solution is achieved using a plane stress finite element code. The computation of the cost function gradient in the parameters' space allows one to iteratively find the best parameters set through a suitable optimization algorithm.;We first design a benchmark including an abrupt rigidity variation that cannot be described by a continuous function. For such a case, we show that increasing the number of parameters is a way to accurately describe sharp variations of the rigidity map. Then, we use a dense GPS velocity field over the southwestern United States to estimate the corresponding rigidity variations for different spatial resolutions of the parameters' grid. We analyse the conceptual and practical difficulties associated with our methodology. Finally, rigidity maps obtained by our inversion method in southwestern United States and particularly across the San Andreas Fault System are reviewed and compared to current plate rigidity estimates and geophysical data over this area.
ISSN: 0956-540X

hal-00669897
https://hal.archives-ouvertes.fr/hal-00669897
DOI : 10.1111/j.1365-246X.2011.05192.x

É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 Springer Verlag
Résumé : International audience
We study the convergence of the Mann Iteration applied to the partial complement of a firmly nonexpansive operator with respect to a linear subspace of a Hilbert space. A new concept considered here. A regularized version is also proposed. Furthermore, to motivate this concept, some applications to robust regression procedures and location problems are proposed.
ISSN: 1862-4472

hal-00776655
https://hal.univ-antilles.fr/hal-00776655
DOI : 10.1007/s11590-009-0172-3

Éditeur(s) : HAL CCSD
Résumé : The convex feasibility problem (CFP) is to find a member of the intersection of finitely many closed convex sets in Euclidean space. When the intersection is empty, one can minimize a proximity function to obtain an approximate solution to the problem. The split feasibility problem (SFP) and the split equality problem (SEP) are generalizations of the CFP. The approximate SFP (ASFP) and approximate SEP (ASEP) involve finding only approximate solutions to the SFP and SEP, respectively. We present here the SSEA, a simultaneous iterative algorithm for solving the ASEP. When this algorithm is applied to the ASFP it resembles closely, but is not equivalent to, the CQ algorithm. The SSEA involves orthogonal projection onto the given closed convex sets. The relaxed SSEA (RSSEA) is an easily implementable variant of the SSEA that uses orthogonal projection onto half-spaces at each step to solve the SEP. The perturbed version of the SSEA (PSSEA) is similar to the RSSEA, but uses orthogonal projection onto a sequence of epi-convergent closed convex sets.
https://hal.univ-antilles.fr/hal-00776640

Droits : info:eu-repo/semantics/OpenAccess
hal-00776640
https://hal.univ-antilles.fr/hal-00776640
https://hal.univ-antilles.fr/hal-00776640/document
https://hal.univ-antilles.fr/hal-00776640/file/CBAM10.pdf