Moment Tensor Analysis of the Acoustic Emission Source in the Rock Damage Process

To further investigate the mechanism of acoustic emission (AE) in the rock fracture experiment, moment tensor analysis was carried out. The AE sources characterized by crack sizes, orientations and fracture modes, are represented by a time-dependent momen

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

Flux closure vortexlike domain structures in ferroelectric thin films

Enhanced piezoresponse force microscopy was used to study flux closure vortexlike structures of 90° ferroelastic domains at the nanoscale in thin ferroelectric lead zirconium titanate (PZT) films. Using an external electric field, a vortexlike structure was induced far away from a grain boundary, indicating that physical edges are not necessary for nucleation contrary to previous suggestions. We demonstrate two different configurations of vortexlike structures, one of which has not been observed before. The stability of these structures is found to be size dependent, supporting previous predictions. © 2010 The American Physical Society.

Uniform stability of a particle approximation of the optimal filter derivative

Sequential Monte Carlo methods, also known as particle methods, are a widely used set of computational tools for inference in non-linear non-Gaussian state-space models. In many applications it may be necessary to compute the sensitivity, or derivative, of the optimal filter with respect to the static parameters of the state-space model; for instance, in order to obtain maximum likelihood model parameters of interest, or to compute the optimal controller in an optimal control problem. In Poyiadjis et al. [2011] an original particle algorithm to compute the filter derivative was proposed and it was shown using numerical examples that the particle estimate was numerically stable in the sense that it did not deteriorate over time. In this paper we substantiate this claim with a detailed theoretical study. Lp bounds and a central limit theorem for this particle approximation of the filter derivative are presented. It is further shown that under mixing conditions these Lp bounds and the asymptotic variance characterized by the central limit theorem are uniformly bounded with respect to the time index. We demon- strate the performance predicted by theory with several numerical examples. We also use the particle approximation of the filter derivative to perform online maximum likelihood parameter estimation for a stochastic volatility model.

Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Influence of thermal residual stresses on the composite macroscopic behavior

The influence of the thermal residual stress on the deformation behavior of a composite has been analyzed with a new micromechanical method. The method is based on secant moduli approximation and a new homogenized effective stress to characterize the plastic state of the matrix. It is found that the generated thermal residual stresses after cooling and their influence on the subsequent deformation behavior depends significantly on the aspect ratio of the inclusions. With prolate inclusions, the presence of thermal residual stresses generate a higher compressive hardening curves of the composite, but it is reversed with oblate inclusions. For particle reinforced composite, thermal residual stresses induce a tensile hardening curve higher than the compressive one and this is in agreement with experimental observations. (C) 1998 Elsevier Science Ltd.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Fatigue-crack closure measurements on 2024-t3 sheet specimens

A computer-controlled procedure has been developed for automatic measurement of the crack opening stress S-op during fatigue tests. A crack opening displacement gauge (GOD meter) is used to obtain digital data on the load versus COD curves. Three methods for deriving S-op from the data sets are compared: (1) a slope method, (2) a tangent lines intersecting method, and (3) a tangent point method. The effect of the position of the COD meter with respect to the crack tip on S-op is studied in tests of 2024-T3 specimens. Results of crack growth and S-op are presented for CA loading with an overload, and with an overload followed by an underload.

Statistical-theory of multiphoton excitation of polyatomic-molecules based on the wigner function

This paper is aimed at establishing a statistical theory of rotational and vibrational excitation of polyatomic molecules by an intense IR laser. Starting from the Wigner function of quantum statistical mechanics, we treat the rotational motion in the classical approximation; the vibrational modes are classified into active ones which are coupled directly with the laser and the background modes which are not coupled with the laser. The reduced Wigner function, i.e., the Wigner function integrated over all background coordinates should satisfy an integro-differential equation. We introduce the idea of ``viscous damping'' to handle the interaction between the active modes and the background. The damping coefficient can be calculated with the aid of the well-known Schwartz–Slawsky–Herzfeld theory. The resulting equation is solved by the method of moment equations. There is only one adjustable parameter in our scheme; it is introduced due to the lack of precise knowledge about the molecular potential. The theory developed in this paper explains satisfactorily the recent absorption experiments of SF6 irradiated by a short pulse CO2 laser, which are in sharp contradiction with the prevailing quasi-continuum theory. We also refined the density of energy levels which is responsible for the muliphoton excitation of polyatomic molecules.