Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

Approximate analytical description for fundamental-mode fields of graded-index fibers: Beyond the Gaussian approximation

An approximate analytical description for fundamental-mode fields of graded-index fibers is explicitly presented by use of the power-series expansion method, the maximum-value condition at the fiber axis, the decay properties of fundamental-mode fields at large distance from the fiber axis, and the approximate modal parameters U obtained from the Gaussian approximation. This analytical description is much more accurate than the Gaussian approximation and at the same time keep the simplicity of the latter. As two special examples, we present the approximate analytical formulas for the fundamental-mode fields of a step profile fiber and a Gaussian profile fiber, and we find that they are both highly accurate in the single-mode range by comparing them with the corresponding exact solutions.

Effects of the frequency chirp on the fields of a chirped Gaussian pulse passing through a hard-edged aperture

Based on the Huygens-Fresnel diffraction integral and Fourier transform, propagation expression of a chirped Gaussian pulse passing through a hard-edged aperture is derived. Intensity distributions of the pulse with different frequency chirp in the near-field and far-field are analyzed in detail by numerical calculations. In the near-field, amplitudes of the intensity peaks generated by the modulation of the hard-edged aperture decrease with increasing the frequency chirp, which results in the improving of the beam uniformity. A physical explanation for the smoothing effect brought by increasing the frequency chirp is given. The smoothing effect is achieved not only in the pulse with Gaussian transverse profile but also in the pulse with Hermite-Gaussian transverse profile when the frequency chirp increases. (C) 2005 Elsevier B.V. All rights reserved.

Modal fields and bending loss analyses of three-layer large flattened mode fibers

Theoretical method to analyze three-layer large flattened mode (LFM) fibers is presented. The modal fields, including the fundamental and higher order modes, and bending loss of the fiber are analyzed. The reason forming the different modal fields is explained and the feasibility to filter out the higher order modes via bending to realize high power, high beam quality fiber laser is given. Comparisons are made with the standard step-index fiber. (c) 2006 Elsevier B.V. All rights reserved.

Fields of apertured polychromatic laser beams with Gaussian and Hermite-Gaussian transverse modes

Starting from the Huygens-Fresnel diffraction integral, the field expressions of apertured polychromatic laser beams with Gaussian and Hermite-Gaussian transverse modes are derived. Influence of the bandwidth on the intensity distributions of the laser beams is analyzed. It is found that when the bandwidth increases, the amplitudes and numbers of the intensity spikes decrease and beam uniformity is improved in the near field and the width of transverse intensity distribution of the apertured beams decreases in the far field. Thus, the smoothing and narrowing effects can be achieved by increasing the bandwidth. Also, these effects are found in the laser beams with Hermite-Gaussian transverse modes as the bandwidth increases.(c) 2006 Elsevier Ltd. All rights reserved.

Electromagnetic fields and transmission properties in tapered hollow metallic waveguides

We analyze the electromagnetic spatital distributions and address an important issue of the transmission properties of spherical transverse-electric (TE) and transverse-magnetic (TM) eigenmodes within a tapered hollow metallic waveguide in detail. Explicit analytical expressions for the spatital distributions of electromagnetic field components, attenuation constant, phase constant and wave impedance are derived. Accurate eigenvalues obtained numerically are used to study the dependences of the transmission properties on the taper angle, the mode as well as the length of the waveguide. It is shown that all modes run continuously from a propagating through a transition to an evanescent region and the value of the attenuation increases as the distance from the cone vertex and the cone angle decrease. A strict distinction between pure propagating and pure evanescent modes cannot be achieved. One mode after the other reaches cutoff in the tapered hollow metallic waveguide as the distance from the cone vertex desreases. (C) 2008 Optical Society of America