A new method of studying the dynamical behavior of the sine-gordon equation

We try to connect the theory of infinite dimensional dynamical systems and nonlinear dynamical methods. The sine-Gordon equation is used to illustrate our method of discussing the dynamical behaviour of infinite dimensional systems. The results agree with those of Bishop and Flesch [SLAM J. Math. Anal. 21 (1990) 1511].

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

Elastodynamic stress intensity factor history for a semi-infinite crack under three-dimensional transient loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.

Spontaneous emission lifetime distribution of infinite line antennas in two-dimensional photonic crystals with finite size

The authors calculate the lifetime distribution functions of spontaneous emission from infinite line antennas embedded in two-dimensional disordered photonic crystals with finite size. The calculations indicate the coexistence of both accelerated and inhibited decay processes in disordered photonic crystals with finite size. The decay behavior of the spontaneous emission from infinite line antennas changes significantly by varying factors such as the line antennas' positions in the disordered photonic crystal, the shape of the crystal, the filling fraction, and the dielectric constant. Moreover, the authors analyze the effect of the degree of disorder on spontaneous emission. (c) 2007 American Institute of Physics.

Synthesis and Characterization of One-Dimensional and Molecular M(tren)InAsS4 (M = Mn, Co, and Zn) Compounds with a Noncondensed AsS33- Unit

Four transition-metal-amine complexes incorporating indium thioarsenates with the general formula M(tren)InAsS4 (M=Mn, Co, and Zn) and a noncondensed AsS33- unit have been prepared and characterized. Single-crystal X-ray diffraction analyses show that compound 1 (M=Mn) crystallizes in the triclinic crystal system (space group: P (1) over bar) and consists of a one-dimensional (1D) inorganic (1)(infinity){[InAsS4](2-)} chain and [Mn(tren)](2+) groups bonded to the opposite sides of an eight-membered In2As2S4 ring along the backbone of the infinite inorganic chains. Compounds 2 (M=Mn), 3 (M=Zn), and 4 (M=Co) are isomorphous molecular compounds. They all crystallize in the monoclinic crystal system (space group: P2(1)/c). The Mn2+ cation of [Mn(tren)](2+) in 1 has a distorted octahedral environment, while the transition-metal cations of [M(tren)](2+) in the other three compounds locate in trigonal-bipyramidal environments.

Poly[di-mu(3)-chloro-mu(2)-trans-1,2-di-4-pyridylethylene-dicopper(I)]: a two-dimensional organic-inorganic hybrid constructed from linear CuCl clusters and bridging ligands

The structure of the title compound, [Cu2Cl2(C12H10N2)](n), contains infinite CuCl staircase-like chains, which lie about inversion centres. The trans-1,2-di-4-pyrid-ylethyl-ene mol-ecules also lie about inversion centres and connect the CuCl chains through Cu-N coordination bonds into a two-dimensional organic-inorganic hybrid network. The planar sheets are stacked along the c axis and associated through weak C-H center dot center dot center dot Cl inter-actions. The results show a reliable structural motif with controllable separation of the CuCl chains by variation of the length of the ligand.

Novel two-dimensional sodium molybdenum phosphates containing {M(enMe)(2)}(2+) bridging groups (M = Ni, Cu; enMe=1,2-diaminopropane)

Two novel organic-inorganic hybrid compounds, (H(2)enMe)(4)(H3O)[Ni(enMe)(2)].[Na3Mo12O52P8(OH)(10)].5H(2)O (1) and (H(2)enMe)(4)(H3O)[Cu(enMe)(2)].[Na3Mo12O52P8(OH)(10)].5H(2)O (2) (enMe = 1,2-diaminopropane), have been hydrothermally synthesized and characterized by elemental analyses, IR, EPR, XPS, UV-Vis spectra and TG analyses. Single crystal X-ray diffraction shows that 1 and 2 are isostructural compounds. Both the compounds exhibit an unusual two-dimensional (2-D) window-like network consisting of one-dimensional (1-D) chains of sodium molybdenum phosphate anions connected by transition metal coordination complexes cations. Compound 1 and 2 represent the first 2-D molybdenum phosphate skeleton pillared by transition metal complex fragments.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Effective elastic moduli of two-dimensional solids with multiple cracks

An accurate method which directly accounts for the interactions between different microcracks is used for analyzing the elastic problem of multiple cracks solids. The effective elastic moduli for randomly oriented cracks and parallel cracks are evaluated for the representative volume element (RVE) with microcracks in infinite media. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method is simple and provides a direct and efficient approach to dealing with elastic solids containing multiple cracks.

Three-dimensional instability of an oscillating viscous flow past a circular cylinder

A systematically numerical study of the sinusoidally oscillating viscous flow around a circular cylinder was performed to investigate vortical instability by solving the three-dimensional incompressible Navier-Stokes equations. The transition from two- to three-dimensional flow structures along the axial direction due to the vortical instability appears, and the three-dimensional structures lie alternatively on the two sides of the cylinder. Numerical study has been taken for the Keulegan-Carpenter( KC) numbers from 1 to 3.2 and frequency parameters from 100 to 600. The force behaviors are also studied by solving the Morison equation. Calculated results agree well with experimental data and theoretical prediction.

Embedded cell model of three-dimensional two-phase composites

An embedded cell model is presented to obtain the effective elastic moduli and the elastic-plastic stress-strain relations of three-dimensional two-phase particulate composites. Each cell consists of an ellipsoidal inclusion surrounded by a finite ellipsoidal matrix that embedded in an infinite matrix. When both matrix and particle are elastic, the effective elastic moduli are derived which is an exact analytic formula without any simplified approximation that can be expressed in an explicit form. Further, the elastic-plastic stress-strain relations are obtained for spherical cells and oblate spheroid cells, in which the matrix is elastic and the particle is elastic-plastic. In addition, the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC) is investigated by a systematic approach [1] in which the matrix is elastic-plastic and the particle is elastic.

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Effects of lattice vibration on the effective mass of quasi-two-dimensional strong-coupling polaron

The effects of lattice vibration on the system in which the electron is weakly coupled with bulk longitudinal optical phonons and strongly coupled with interface optical phonons in an infinite quantum well were studied by using Tokuda' linear-combination operator and a modified LLP variational method. The expressions for the effective mass of the polaron in a quantum well QW as functions of the well's width and temperature were derived. In particular, the law of the change of the vibration frequency of the polaron changing with well' s width and temperature are obtained. Numerical results of the effective mass and the vibration frequency of the polaron for KI/AgCl/Kl QW show that the vibration frequency and the effective mass of the polaron decrease with increasing well's width and temperature, but the contribution of the interaction between the electron and the different branches of phonons to the effective mass and the vibration frequency and the change of their variation with the well's width and temperature are greatly different.