A comprehensive study of the 26th December 2004 Sumatra earthquake: possible implications on Earth rotation and investigations on the coseismic and postseismic stress diffusion associated with the seismic rupture

In this work a multidisciplinary study of the December 26th, 2004 Sumatra earthquake has been carried out. We have investigated both the effect of the earthquake on the Earth rotation and the stress field variations associated with the seismic event. In the first part of the work we have quantified the effects of a water mass redistribution associated with the propagation of a tsunami wave on the Earth’s pole path and on the length-of-day (LOD) and applied our modeling results to the tsunami following the 2004 giant Sumatra earthquake. We compared the result of our simulations on the instantaneous rotational axis variations with some preliminary instrumental evidences on the pole path perturbation (which has not been confirmed yet) registered just after the occurrence of the earthquake, which showed a step-like discontinuity that cannot be attributed to the effect of a seismic dislocation. Our results show that the perturbation induced by the tsunami on the instantaneous rotational pole is characterized by a step-like discontinuity, which is compatible with the observations but its magnitude turns out to be almost one hundred times smaller than the detected one. The LOD variation induced by the water mass redistribution turns out to be not significant because the total effect is smaller than current measurements uncertainties. In the second part of this work of thesis we modeled the coseismic and postseismic stress evolution following the Sumatra earthquake. By means of a semi-analytical, viscoelastic, spherical model of global postseismic deformation and a numerical finite-element approach, we performed an analysis of the stress diffusion following the earthquake in the near and far field of the mainshock source. We evaluated the stress changes due to the Sumatra earthquake by projecting the Coulomb stress over the sequence of aftershocks taken from various catalogues in a time window spanning about two years and finally analyzed the spatio-temporal pattern. The analysis performed with the semi-analytical and the finite-element modeling gives a complex picture of the stress diffusion, in the area under study, after the Sumatra earthquake. We believe that the results obtained with the analytical method suffer heavily for the restrictions imposed, on the hypocentral depths of the aftershocks, in order to obtain the convergence of the harmonic series of the stress components. On the contrary we imposed no constraints on the numerical method so we expect that the results obtained give a more realistic description of the stress variations pattern.

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.

Integrating 40Ar-39Ar, 87Rb-87Sr and 147Sm-143Nd geochronology of authigenic illite to evaluate tectonic reactivation in an intraplate setting, central Australia

The Warburton-Cooper basins, central Australia, include a multitude of reactivated fracture-fault networks related to a complex, and poorly understood, tectonic evolution. We investigated authigenic illites from a granitic intrusion and sedimentary rocks associated with prominent structural features (Gidgealpa-Merrimelia-Innamincka Ridge and the Nappamerri Trough). These were analysed by 40Ar-39Ar, 87Rb-87Sr and 147Sm-143Nd geochronology to explore the thermal and tectonic histories of central Australian basins. The combined age data provide evidence for three major periods of fault reactivation throughout the Phanerozoic. While Carboniferous (323.3 ± 9.4 Ma) and Late Triassic ages (201.7 ± 9.3 Ma) derive from basin-wide hydrothermal circulation, Cretaceous ages (~128 to ~86 Ma) reflect episodic fluid flow events restricted to the synclinal Nappamerri Trough. Such events result from regional extensional tectonism derived from the transferral of far-field stresses to mechanically and thermally weakened regions of the Australian continent. Specifically, Cretaceous ages reflect continent-wide transmission of tensional stress from a > 2500 km long rifting event on the Eastern (and southern) Australian margin associated with break-up of Gondwana and opening of the Tasman Sea. By integrating 40Ar-39Ar, 87Rb-87Sr and 147Sm-143Nd dating, this study highlights the use of authigenic illite in temporally constraining the tectonic evolution of intracontinental basins that would otherwise remain unknown. Furthermore, combining Sr- and Ar-isotopic systems enables more accurate dating of authigenesis whilst significantly reducing geochemical pitfalls commonly associated with these radioisotopic dating methods.

Off-axis field approximations for ion traps with apertures

In recent work (Int. J. Mass Spec., vol. 282, pp. 112–122) we have considered the effect of apertures on the fields inside rf traps at points on the trap axis. We now complement and complete that work by considering off-axis fields in axially symmetric (referred to as “3D”) and in two dimensional (“2D”) ion traps whose electrodes have apertures, i.e., holes in 3D and slits in 2D. Our approximation has two parts. The first, EnoAperture, is the field obtained numerically for the trap under study with apertures artificially closed. We have used the boundary element method (BEM) for obtaining this field. The second part, EdueToAperture, is an analytical expression for the field contribution of the aperture. In EdueToAperture, aperture size is a free parameter. A key element in our approximation is the electrostatic field near an infinite thin plate with an aperture, and with different constant-valued far field intensities on either side. Compact expressions for this field can be found using separation of variables, wherein the choice of coordinate system is crucial. This field is, in turn, used four times within our trap-specific approximation. The off-axis field expressions for the 3D geometries were tested on the quadrupole ion trap (QIT) and the cylindrical ion trap (CIT), and the corresponding expressions for the 2D geometries were tested on the linear ion trap (LIT) and the rectilinear ion trap (RIT). For each geometry, we have considered apertures which are 10%, 30%, and 50% of the trap dimension. We have found that our analytical correction term EdueToAperture, though based on a classical small-aperture approximation, gives good results even for relatively large apertures.

Electric Field Modelling of Composite High Voltage Insulators

This paper deals with the two-dimensional electric field modelling and electric field stress calculations of different types of composite insulators used in high voltage distribution and transmission systems. The computer simulations are carried out by using a commercially available software package. The potential and electric filed results obtained for the actual insulator profiles for three types of composite/polymeric insulators are discussed and presented.

Thermomechanical behaviour of pin joints subjected to sheet loads

An important problem regarding pin joints in a thermal environment is addressed. The motivation emerges from structural safety requirements in nuclear and aerospace engineering. A two-dimensional model of a smooth, rigid misfit pin in a large isotropic sheet is considered as an abstraction. The sheet is subjected to a biaxial stress system and far-field unidirectional heat flow. The thermoelastic analysis is complex due to non-linear load-dependent contact and separation conditions at the pin-hole interface and the absence of existence and uniqueness theorems for the class of frictionless thermoelastic contact problems. Identification of relevant parameters and appropriate synthesis of thermal and mechanical variables enables the thermomechanical generalization of pin-joint behaviour. This paper then proceeds to explore the possibility of multiple solutions in such problems, especially interface contact configuration.

Estimation of electric stress and surface potential for traction insulators

Traction insulators are solid core insulators widely used for railway electrification. Constant exposure to detrimental effects of vandalism, and mechanical vibrations begets certain faults like shorting of sheds or cracks in the sheds. Due to fault in one/two sheds, stress on the remaining healthy sheds increases, owing to atmospheric pollution the stress may lead to a flashover of the insulator. Presently due to non availability of the electric stress data for the insulators, simulation study is carried out to find the potential and electric field for most widely used traction insulators in the country. The results of potential and electric field stress obtained for normal and faulty imposed insulators are presented.

An exact solution for the three-phase piezoelectric cylinder model under antiplane shear and its applications to piezoelectric composites

A three-phase piezoelectric cylinder model is proposed and an exact solution is obtained for the model under a farfield antiplane mechanical load and a far-field inplane electrical load. The three-phase model can serve as a fiber/interphase layer/matrix model, in terms of which a lot of interesting mechanical and electrical coupling phenomena induced by the interphase layer are revealed. It is found that much more serious stress and electrical field concentrations occur in the model with the interphase layer than those without any interphase layer. The three-phase model can also serve as a fiber/matrix/composite model, in terms of which a generalized self-consistent approach is developed for predicting the effective electroelastic moduli of piezoelectric composites. Numerical examples are given and discussed in detail.

The strength distribution of brittle materials with a high-concentration of cavities

The maximum stress concentration factor in brittle materials with a high concentration of cavities is obtained. The interaction between the nearest cavities, in addition to the far field interactions, is taken into account to evaluate the strength distribution based on the statistical analysis of the nearest distance distribution. Through this investigation, it is found that the interaction between the nearest neighbors is much more important than the far field interactions, and one has to consider it in calculating the strength of brittle materials even if the volume fraction of cavities it contains is small. The other important conclusion is that the maximum stress concentration factor has a wide scattered distribution.

Some theorems in classical elastodynamic

This investigation is concerned with various fundamental aspects of the linearized dynamical theory for mechanically homogeneous and isotropic elastic solids. First, the uniqueness and reciprocal theorems of dynamic elasticity are extended to unbounded domains with the aid of a generalized energy identity and a lemma on the prolonged quiescence of the far field, which are established for this purpose. Next, the basic singular solutions of elastodynamics are studied and used to generate systematically Love's integral identity for the displacement field, as well as an associated identity for the field of stress. These results, in conjunction with suitably defined Green's functions, are applied to the construction of integral representations for the solution of the first and second boundary-initial value problem. Finally, a uniqueness theorem for dynamic concentrated-load problems is obtained.

Numerical analysis and comparison of three metal-oxide-type super-resolution near field structures

Based on the Fresnel-Kirchkoff diffraction theory, we build up a Gaussian diffraction model of metal-oxide-type super-resolution near field structure (super-RENS), which can describe far field optical properties. The spectral contrast induced by refractive index and the structural changes in AgOx, PtOx and PdOx thin films, which are the key functional layers in super-RENS, are studied by using this model. Comparison results indicate that the spectral contrast intensively on laser-induced distribution and change of the refractive index in the metal-oxide films. The readout mechanism of the metal-oxide-type super-RENS optical disc is further clarified. This Gaussian diffraction model can be used as a simple and effective method for choosing proper active materials in super-RENS.

Structure-change-dependent transmittance of PtOx thin film in super-resolution near-field structure: Diffraction effect analysis

To study working mechanism of super-resolution near-field structure (super-RENS) optical disk from a far-field optics view is very necessary because of the actual far-field writing/readout process in the optical disk system. A Gaussian diffraction model based on Fresnel-Kirchhoff diffraction theory of PtOx-type super-RENS has been set up in this Letter. The relationship between micro-structural deformation (change of bubble structure and refractive index profile) with far-field optical response of PtOx thin film has been studied with it in detail. The simulation results are in good agreement with the experimental results reported in literatures with a designed configuration. These results may provide more quantitative information for better understanding of the working mechanism of metal-oxide-type super-RENS. (c) 2007 Elsevier B.V. All rights reserved.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Scalar and its dissipation in the near field of turbulent lifted jet flame

In this study various scalar dissipation rates and their modelling in the context of partially premixed flame are investigated. A DNS dataset of the near field of a turbulent hydrogen lifted jet flame is processed to analyse the mixture fraction and progress variable dissipation rates and their cross dissipation rate at several axial positions. It is found that the classical model for the passive scalar dissipation rate ε{lunate}̃ZZ gives good agreement with the DNS, while models developed based on premixed flames for the reactive scalar dissipation rate ε{lunate}̃cc only qualitatively capture the correct trend. The cross dissipation rate ε{lunate}̃cZ is mostly negative and can be reasonably approximated at downstream positions once ε{lunate}̃ZZ and ε{lunate}̃cc are known, although the sign cannot be determined. This approach gives better results than one employing a constant ratio of turbulent timescale and the scalar covariance c'Z'̃. The statistics of scalar gradients are further examined and lognormal distributions are shown to be very good approximations for the passive scalar and acceptable for the reactive scalar. The correlation between the two gradients increases downstream as the partially premixed flame in the near field evolves ultimately to a diffusion flame in the far field. A bivariate lognormal distribution is tested and found to be a reasonable approximation for the joint PDF of the two scalar gradients. © 2011 The Combustion Institute.

High-output very small aperture laser and its near-field distribution properties - art. no. 60200L

A Very-Small-Aperture Laser with a 250 X 500 nm(2) aperture has been created on a 650nm edge emitting LD. The highest far-field output power is 1.9mW and the power per unit emission area is about 15 MW/mu m(2). The special fabrication process and high output power mechanism are demonstrated respectively. The near-field distribution properties are also analyzed theoretically and experimentally.