Cette thèse vise à étudier les propriétés thermo-poro-mécaniques de roches de faille, à partir de l'analyse structurale et microstructurale d'une faille aujourd'hui à l'affleurement et à partir d'expériences menées en laboratoire, en vue de déterminer les processus qui contrôlent l'efficacité de deux mécanismes responsables de d'affaiblissement cosismiques : la pressurisation thermique et le mécanisme d'affaiblissement par drainage. L'étude de terrain a été conduite sur deux affleurements appartenant à une faille décrochante potentiellement active appartenant au système de failles du Chugoku occidental (Japon) : la faille d'Usukidani. Le travail expérimental a quant à lui été mené dans le laboratoire de déformation des roches de l'Université de Kyoto. Les résultats majeurs de ce travail sont exposés ci-dessous. <br /> Les propriétés hydrologiques et poro-élastiques de la gouge et de la brèche de la faille d'Usukidani ont été determinées à partir d'échantillons prélevés sur le terrain. Ces données hydrauliques ont ensuite été utilisées dans un modèle numérique afin d'évaluer l'importance du phénomène de pressurisation thermique dans le cas d'un glissement cosismique le long de la zone de glissement principale et le long de zones de glissement secondaires. Les résultats de cette modélisation suggèrent que le mécanisme de pressurisation thermique n'est efficace que si la rupture reste localisée le long des zones de glissement contenant de la gouge, avec comme facteur de contrôle l'épaisseur de cette zone de glissement. <br /> Afin d'identifier les processus dynamiques particulaires responsables de l'affaiblissement cosismique dans la zone de glissement, plusieurs essais de friction ont été menés sur une machine à cisaillement annulaire. Ces expériences ont été conduites à des vitesses cosismiques (équivalentes à 0,09, 0,9 et 1,3 m/s) en conditions humides ou conditions sèches. Les données obtenues montrent que quelles que soient les conditions d'humidité initiales, les failles simulées montrent toutes un affaiblissement lors du déplacement. Un examen détaillé des microstructures des gouges cisaillées obtenues une fois l'équilibre frictionnel atteint permet de définir deux types de microstructures impliquant deux régimes de déformation : un régime de déformation par roulement avec la formation d'agrégats argileux, et un régime de déformation par glissement avec la formation d'une zone de cisaillement complexe localisée à l'interface gouge-éponte. L'affaiblissement observé lors des expériences semble être lié à une diminution de la proportion de grains roulants par rapport à celle de grains glissants, et semble être favorisé par le développement des agrégats argileux, lesquels sont contrôlés par la teneur en eau.<br /> A partir d'un modèle numérique (P2 FEM) et des données de contrainte cisaillante obtenues lors des essais de friction, il a été possible de calculer l'évolution de la température de la gouge en fonction du déplacement. Les résultats suggèrent que la distance dc pourrait représenter la distance nécessaire à une faille pour produire et diffuser assez de chaleur afin de casser les ponts d'eau capillaire (forces d'adhésion) et ainsi permettre à l'eau contenue dans la gouge d'être libérée. Ce mécanisme est appelé mécanisme d'affaiblissement par drainage.

<p>The elastic properties of oil/water/surfactant interfaces play an important role in the phase behaviour of microemulsions and for the stability of macroemulsions. The aim of this thesis is to obtain an understanding of the relationship between the structure of the surfactant molecules, the structure of the interface, and macroscopic interfacial properties. To achieve this aim, we performed molecular simulations of oil/water/surfactant systems. We made a quantitative comparison of various model surfactants to determine how structural changes affect interfacial properties and film rupture. The model consists of water, oil, head, and tail beads, and surfactants are constructed by coupling head and tail beads with harmonic springs. We used a hybrid dissipative particle dynamics-Monte Carlo scheme. The former was used to simulate particle trajectories and the Monte Carlo scheme was used to mimic experimental conditions: bulk-interface phase equilibrium, tensionless interfaces in microemulsions, and the surface force apparatus.</p><p>A detailed comparison of various non-ionic model surfactants showed how structural changes affect interfacial properties:</p><p>Comparison between linear and branched surfactants showed that the efficiency of adsorption is higher for linear surfactants, although branched surfactants are more efficient at a given surface density. Linear surfactants can be more efficient also at the same surface density if the head group is sufficiently soluble in oil, because low head-oil repulsion makes the branched isomers stagger at the interface. The bending rigidity is higher for linear surfactants. Furthermore, branched surfactants make oil droplets coalesce more easily than linear surfactants do, but linear and branched surfactants have roughly the same effect on water droplet coalescence. </p><p>Comparison of linear surfactants with varying chain lengths showed that longer surfactants have a lower surface tension and higher bending rigidity. The increase in rigidity with chain length follows a power law, but the exponent is higher for surfactant monolayers at a fixed density than at a fixed tension. Longer tails and/or denser monolayers influence the stability of water droplets in a positive direction, and the stability of oil droplets in a negative direction. </p><p>Addition of cosurfactant showed that mixed monolayers have a lower bending rigidity than pure monolayers at the same average chain length and tension. Cosurfactants have a negative effect on the stability of water droplets, and a positive effect on the stability of oil droplets.</p> / Paper I reprinted with kind permission of EDP Sciences. Paper III reprinted with kind permission of the American Institute of Physics. Paper IV reprinted with kind permission of the American Physics Society.

In this Thesis, two problems were studied: a direct vacuum acceleration of electrons by a tightly focused ultrashort relativistic laser pulse and ion acceleration in the process of spherical laser-heated plasma explosion. The electromagnetic field of a tightly focused laser pulse was evaluated numerically by means of Stratton-Chu integrals. The properties of the focused field were analyzed in detail for a plane wave or a macroscopically large Gaussian beam incident onto the mirror. Free electrons moving in the tightly focused field were found to accelerate by two possible mechanisms: focal spot acceleration and capture-and-acceleration scenario. The two mechanisms were studied in detail. Comparison of the mirror-focused field with first- and fifth-order paraxial fields is performed. A 3D electromagnetic PIC code SCPIC was created for simulations of pulse interaction with targets having a finite number of particles interacting with each other by collective fields. Atto-second bunch formation was observed in the interaction with ultra-small or ultra-thin targets. Physical mechanism of bunch formation is explained. The problem of electrostatic explosion of a nano-scale spherical plasma with initially hot electrons and cold ions was solved numerically. Expansion in a wide regime of electron temperature $0 < T leq infty$ was studied in detail for different initial density profiles of plasma. Favorable conditions for obtaining mono-energetic ions resulting from the explosion were specified in single and two ionic species cases. In case of a two-species explosion, the number of mono-energetic, $deltavarepsilon/varepsilon < 10\%$, ions can be as high as 70-80\% of the total light ions for a wide range of electron temperatures.

The elastic properties of oil/water/surfactant interfaces play an important role in the phase behaviour of microemulsions and for the stability of macroemulsions. The aim of this thesis is to obtain an understanding of the relationship between the structure of the surfactant molecules, the structure of the interface, and macroscopic interfacial properties. To achieve this aim, we performed molecular simulations of oil/water/surfactant systems. We made a quantitative comparison of various model surfactants to determine how structural changes affect interfacial properties and film rupture. The model consists of water, oil, head, and tail beads, and surfactants are constructed by coupling head and tail beads with harmonic springs. We used a hybrid dissipative particle dynamics-Monte Carlo scheme. The former was used to simulate particle trajectories and the Monte Carlo scheme was used to mimic experimental conditions: bulk-interface phase equilibrium, tensionless interfaces in microemulsions, and the surface force apparatus. A detailed comparison of various non-ionic model surfactants showed how structural changes affect interfacial properties: Comparison between linear and branched surfactants showed that the efficiency of adsorption is higher for linear surfactants, although branched surfactants are more efficient at a given surface density. Linear surfactants can be more efficient also at the same surface density if the head group is sufficiently soluble in oil, because low head-oil repulsion makes the branched isomers stagger at the interface. The bending rigidity is higher for linear surfactants. Furthermore, branched surfactants make oil droplets coalesce more easily than linear surfactants do, but linear and branched surfactants have roughly the same effect on water droplet coalescence. Comparison of linear surfactants with varying chain lengths showed that longer surfactants have a lower surface tension and higher bending rigidity. The increase in rigidity with chain length follows a power law, but the exponent is higher for surfactant monolayers at a fixed density than at a fixed tension. Longer tails and/or denser monolayers influence the stability of water droplets in a positive direction, and the stability of oil droplets in a negative direction. Addition of cosurfactant showed that mixed monolayers have a lower bending rigidity than pure monolayers at the same average chain length and tension. Cosurfactants have a negative effect on the stability of water droplets, and a positive effect on the stability of oil droplets. / Paper I reprinted with kind permission of EDP Sciences. Paper III reprinted with kind permission of the American Institute of Physics. Paper IV reprinted with kind permission of the American Physics Society.

A granular material is a collection of discrete, solid particles of macroscopic size dispersed in an interstitial fluid, in which the fluid has an insignificant effect on the particle dynamics. Because they exhibit fascinating properties because of dissipative interactions, due to their importance in geophysical and industrial processes, flows of granular materials have been the focus of large amount of research involving physicists and engineers. A good understanding of the physics of granular materials is desired in order to design efficient processing and handling systems. Granular materials can be heaped like a solid, and can flow like a fluid. Though the two distinct regimes of granular flows are well described by kinetic theory (rapid flows) and plasticity theories (quasi-static), the intermediate dense flow regime, where collisional and frictional interactions are important, is not yet described successfully. In this thesis, we examine the applicability of kinetic theory for dense granular flows, the structure and dynamics in sheared inelastic hard disks systems and dynamics of sheared non-spherical particles. Two complementary simulation techniques, the discrete element (DE) technique for soft particles and the event driven (ED) simulation technique for hard particles, are used to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model when the particle stiffness becomes large. First, we examine the average co-ordination number for the particles in the flow down an inclined plane using the DE technique using both linear and Hertzian contact models. The simulations show that the average co-ordination number decreases below 1 for values of the spring stiffness corresponding to real materials such as sand and glass, even when the angle of inclination is only 1olarger than the angle of repose. The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are found to be in quantitative agreement. In addition, we also conduct the comparison of the pre-collisional relative velocities of particles in contact. Since momentum is transported primarily by particle contacts in a dense flow, the relative velocity distribution is a sensitive comparison of the dynamics in the two simulation techniques. It is found that the relative velocity distribution in both simulation techniques are well approximated by an exponential distribution for small coefficients of restitution, indicating that the dynamics of a dense granular flow can be adequately described by a hard particle model. The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analysed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analysed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analysed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated. The orientational ordering and dynamical properties of the shear flow of inelastic dumbbells in two dimensions are studied, as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are smooth fused disks characterised by the ratio of the distance between centers (L) and the disk diameter (D), and the ratio (L/D)varies between 0 and 1 in our simulations. Area fractions studied are in the range 0.1 to 0.7, while coefficients of normal restitution from 0.99 to 0.6 are considered. The simulations are similar to the event driven simulations for circular disks, but the procedure for predicting collisions is much more complicated due to the non-circular shape of the particles and due to particle rotation. The average orientation is measured using an orientational order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axis of all dumbbells in the same direction). It is found that there is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased, and the order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. However, there is no discontinuous nematic transition for all the parameters studied here. The axis of the dumbbells are preferentially oriented along the extensional axis (at an angle of 45ofrom the flow direction) at low area fraction, but the orientation is closer to the flow direction as the area fraction is increased. The orientation distribution is calculated, and it is found that the orientation distribution is well described by a function of the form P(θ) =(1/π)+ (2S/π)cos(2(θ−θp)), where θis the angle from the flow direction and θpis the principal orientation direction. The mean energy of the velocity fluctuations in the flow direction is found to be higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocity are found to be Gaussian distributions to a very good approximation. The equation of state for the pressure is calculated, and it is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation even at the lowest area fractions studied here. The mean angular velocity of the particles is examined in some detail. It is found that the mean angular velocity is equal to half the vorticity at low area fractions, but the magnitude of the mean angular velocity systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

In this thesis we studied theoretically and numerically dynamics, order and fluctuations in two dimensional active matter with specific reference to the nematic phase in collections of self-driven particles.The aim is to study the ways in which a nonequilibrium steady state with nematic order differs from a thermal equilibrium system of the same spatial symmetry. The models we study are closely related to “flocking”[1], as well as to equations written down to describe the interaction of molecular motors and filaments in a living cell[2,3] and granular nematics [4]. We look at (i) orientational and density fluctuations in the ordered phase, (ii) the way in which density fluctuations evolve in a nematic background, and finally (iii) the coarsening of nematic order and the density field starting from a statistically homogeneous and isotropic initial state. Our work establishes several striking differences between active nematics and their thermal equilibrium counterparts. We studied two-dimensional nonequilibrium active nematics. Two-dimensional nonequilibrium nematic steady states, as found in agitated granular-rod monolayers or films of orientable amoeboid cells, were predicted [5] to have giant number fluctuations, with the standard deviation proportional to the mean. We studied this problem more closely, asking in particular whether the active nematic steady state is intrinsically phase-separated. Our work has close analogy to the work of Das and Barma[6] on particles sliding downhill on fluctuating surfaces, so we looked at a model in which particles were advected passively by the broken-symmetry modes of a nematic, via a rule proposed in [5]. We found that an initially homogeneous distribution of particles on a well-ordered nematic background clumped spontaneously, with domains growing as t1/2, and an apparently finite phase-separation order parameter in the limit of large system size. The density correlation function shows a cusp, indicating that Porod’s Law does not hold here and that the phase-separation is fluctuation-dominated[7]. Dynamics of active particles can be implemented either through microscopic rules as in[8,9]or in a long-wavelength phenomenological approach as in[5]It is important to understand how the two methods are related. The purely phenomenological approach introduces the simplest possible (and generally additive)noise consistent with conservation laws and symmetries. Deriving the long-wavelength equation by explicit coarse-graining of the microscopic rule will in general give additive and multiplicative noise terms, as seen in e.g., in [10]. We carry out such a derivation and obtain coupled fluctuating hydrodynamic equations for the orientational order parameter (polar as well as apolar) and density fields. The nonequilibrium “curvature-induced” current term postulated on symmetry grounds in[5]emerges naturally from this approach. In addition, we find a multiplicative contribution to the noise whose presence should be of importance during coarsening[11]. We studied nonequilibrium phenomena in detail by solving stochastic partial differential equations for apolar objects as obtained from microscopic rules in[8]. As a result of “curvature-induced” currents, the growth of nematic order from an initially isotropic, homogeneous state is shown to be accompanied by a remarkable clumping of the number density around topological defects. The consequent coarsening of both density and nematic order are characterised by cusps in the short-distance behaviour of the correlation functions, a breakdown of Porod’s Law. We identify the origins of this breakdown; in particular, the nature of the noise terms in the equations of motion is shown to play a key role[12]. Lastly we studied an active nematic steady-state, in two space dimensions, keeping track of only the orientational order parameter, and not the density. We apply the Dynamic Renormalization Group to the equations of motion of the order parameter. Our aim is to check whether certain characteristic nonlinearities entering these equations lead to singular renormalizations of the director stiffness coefficients, which would stabilize true long-range order in a two-dimensional active nematic, unlike in its thermal equilibrium counterpart. The nonlinearities are related to those in[13]but free of a constraint that applies at thermal equilibrium. We explore, in particular, the intriguing but ultimately deceptive similarity between a limiting case of our model and the fluctuating Burgers/KPZequation. By contrast with that case, we find that the nonlinearities are marginally irrelevant. This implies in particular that 2-dactive nematics too have only quasi-long-range order[14].

No description available.

We consider fast computation methods for simulation of dynamics of a collection of particles dispersed in an unbounded Stokesian suspension. Stokesian suspensions are of great practical interest in the manufacturing and processing of various commercial products. The most popular dynamic simulation method for these kind of suspensions was developed by Brady and Bossis (Brady and Bossis [1988]). This method uses a truncated multipole expansion to represent the fluid traction on particle surfaces. The hydrodynamic interactions in Stoke-sian suspension are long ranged in nature, resulting in strong coupled motion of all particles. For an N particle system, this method imposes an O(N3) computational cost, thus posing limitations to the number of particles that may be simulated. More recent methods (Sierou and Brady [2001], Scintilla, Darve and Shaqfeh [2005]) have attempted to solve this problem using Particle Mesh Ewald summation techniques by distributing the moments on a grid and using Fast Fourier Transform algorithms, resulting in an O(N log N) computational cost. We review these methods and propose a version that we believe is some-what superior. In the course of this study, we have identified and corrected errors in previous studies that maybe of some importance in determining the bulk properties of suspensions. Finally, we show the utility of our method in determining certain properties of suspensions and compare them to existing analytical results for the same.

No description available.

No description available.