Numerical Mathematics Theory,Methods and Applications

Aims and Scope: Numerical Mathematics:Theory, Methods and Applications (NM-TMA) publishes high-quality original research papers on the construction,analysis and application of numerical methods for solving scientific problems.Important research and expository papers devoted to the numerical solution of mathematical problems arising in all areas of science and technology are expected.The journal originates from the journal Numerical Mathematics:A Journal of Chinese Universities (English Edition).

Numerical Mathematics Theory,Methods and Applications

Aims and Scope: Numerical Mathematics:Theory, Methods and Applications(NM-TMA) publishes high-quality original

Real-time 3-D space numerical shake prediction for earthquake early warning

In earthquake early warning systems, real-time shake prediction through wave propagation simulation is a promising approach. Compared with traditional methods, it does not suffer from the inaccurate estimation of source parameters. For computation efficiency, wave direction is assumed to propagate on the 2-D surface of the earth in these methods. In fact, since the seismic wave propagates in the 3-D sphere of the earth, the 2-D space modeling of wave direction results in inaccurate wave estimation. In this paper, we propose a 3-D space numerical shake pre- diction method, which simulates the wave propagation in 3-D space using radiative transfer theory, and incorporate data assimilation technique to estimate the distribution of wave energy. 2011 Tohoku earthquake is studied as an example to show the validity of the proposed model. 2-D space model and 3-D space model are compared in this article, and the prediction results show that numerical shake prediction based on 3-D space model can estimate the real-time ground motion precisely, and overprediction is alleviated when using 3-D space model.

Focusing of elliptically polarized Gaussian beams through an annular high numerical aperture

Based on the vectorial Debye theory, the focusing properties of the Gaussian beam through an annular high numerical aperture are studied numerically, including the intensity, the phase and the orbital angular momentum properties. Then the influence of certain parameters on the focusing properties is also investigated. It is shown that sub-wavelength elliptical light spots can be obtained. And there exists a vortex in the longitudinal component of the focused field even though the incident beam is Gaussian beam, indicating that the spin angular momentum of the elliptically polarized Gaussian beam is converted into the orbital angular momentum by the focusing.

Real-time numerical shake prediction and updating for earthquake early warning

Ground motion prediction is important for earthquake early warning systems, because the region＇s peak ground motion indicates the potential disaster. In order to predict the peak ground motion quickly and pre- cisely with limited station wave records, we propose a real- time numerical shake prediction and updating method. Our method first predicts the ground motion based on the ground motion prediction equation after P waves detection of several stations, denoted as the initial prediction. In order to correct the prediction error of the initial prediction, an updating scheme based on real-time simulation of wave propagation is designed. Data assimilation technique is incorporated to predict the distribution of seismic wave energy precisely. Radiative transfer theory and Monte Carlo simulation are used for modeling wave propagation in 2-D space, and the peak ground motion is calculated as quickly as possible. Our method has potential to predict shakemap, making the potential disaster be predicted before the real disaster happens. 2008 Ms8.0 Wenchuan earthquake is studied as an example to show the validity of the proposed method.

Fracture initiation and propagation in intact rock—A review

The initiation and propagation of failure in intact rock are a matter of fundamental importance in rock engineering. At low confining pressures, tensile fracturing initiates in samples at 40%-60% of the uniaxial compressive strength and as loading continues, and these tensile fractures increase in density, ultimately coalescing and leading to strain localization and macro-scale shear failure of the samples. The Griffith theory of brittle failure provides a simplified model and a useful basis for discussion of this process. The Hoek-Brown failure criterion provides an acceptable estimate of the peak strength for shear failure but a cutoff has been added for tensile conditions. However, neither of these criteria adequately explains the progressive coalition of tensile cracks and the final shearing of the specimens at higher confining stresses. Grain-based numerical models, in which the grain size distributions as well as the physical properties of the component grains of the rock are incorporated, have proved to be very useful in studying these more complex fracture processes.

Distinguishing Oil and Water Layers in a Porous Cracked Medium by Interpreting Acoustic Logging Data on the Basis of Hudson Theory

During surveys, water layers may interfere with the detection of oil layers. In order to distinguish between oil and water layers in a porous cracked medium, research on the properties of cracks and oil and water layers and their relation to acoustic logging rules is essential. On the basis of Hudson＇s crack theory, we simulated oil and water layers in crack-porous medium with different crack parameters corresponding to the well-field response. We found that in a cracked medium with high crack angle or low number density of cracks, compressional and shear wave velocities are sensitive to crack characteristics; further, these velocities are more sensitive to crack characteristics when the waves propagate through the water layer than when they propagate through the oil layer. Compressional and shear wave velocities increase with an increase in crack angle： in the water layer, the increase is approximately linear. On comparing the full waveforms observed in the oil and water layers, we find that the amplitudes of most waves are higher in the water layer. Among the considered waves, the Stoneley wave suffers maximum amplitude attenuation in the oil layer. The maximum excitation intensity for oil layer is greater than that for the water layer. These results can guide further cracked media logging field exploration work.

ON BLOCK MATRICES ASSOCIATED WITH DISCRETE TRIGONOMETRIC TRANSFORMS AND THEIR USE IN THE THEORY OF WAVE PROPAGATION

Block matrices associated with discrete Trigonometric transforms （DTT＇s） arise in the mathematical modelling of several applications of wave propagation theory including discretizations of scatterers and radiators with the Method of Moments, the Boundary Element Method, and the Method of Auxiliary Sources. The DTT＇s are represented by the Fourier, Hartley, Cosine, and Sine matrices, which are unitary and offer simultaneous diagonalizations of specific matrix algebras. The main tool for the investigation of the aforementioned wave applications is the efficient inversion of such types of block matrices. To this direction, in this paper we develop an efficient algorithm for the inversion of matrices with U-diagonalizable blocks （U a fixed unitary matrix） by utilizing the U- diagonalization of each block and subsequently a similarity transformation procedure. We determine the developed method＇s computational complexity and point out its high efficiency compared to standard inversion techniques. An implementation of the algorithm in Matlab is given. Several numerical results are presented demonstrating the CPU-time efficiency and accuracy for ill-conditioned matrices of the method. The investigated matrices stem from real-world wave propagation applications.