An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

THEORETICAL RESULTS ON THE EXISTENCE,REGULARITY AND ASYMPTOTIC STABILITY OF ENHANCED PULLBACK ATTRACTORS:APPLICATIONS TO 3D PRIMITIVE EQUATIONS

Several new concepts of enhanced pullback attractors for nonautonomous dynamical systems are introduced here by uniformly enhancing the compactness and attraction of the usual pullback attractors over an infinite forward time-interval under strong and weak topologies.Then we provide some theoretical results for the existence,regularity and asymptotic stability of these enhanced pullback attractors under general theoretical frameworks which can be applied to a large class of PDEs.The existence of these enhanced attractors is harder to obtain than the backward case[33],since it is difficult to uniformly control the long-time pullback behavior of the systems over the forward time-interval.As applications of our theoretical results,we consider the famous 3D primitive equations modelling the large-scale ocean and atmosphere dynamics,and prove the existence,regularity and asymptotic stability of the enhanced pullback attractors in V×V and H^(2)×H^(2) for the time-dependent forces which satisfy some weak conditions.

Solution of Cauchy's Problem for Wave Equations in Higher Space Dimensions by Means of D'Alembert's Formula

A simple method for solving Cauchy’s problem of wave equations in higher space dimensions with initial condition of separated variables, has been given by using D’Alembert’s formula and some examples have been shown.

ON THE SUPERSTABILITY OF THE PEXIDER TYPE GENERALIZED TRIGONOMETRIC FUNCTIONAL EQUATIONS

The aim of this paper is to investigate the superstability problem for the pexiderized trigonometric functional equation∑ v∈Φ∫Kf（xkv（y）k^-1）dwK（k）= Φ g（x）h（y）, x, y ∈ G,where G is any topological group, K is a compact subgroup of G, ωK is the normalized Haar measure of K, Φ is a finite group of K-invariant morphisms of G and f, g, h are continuous complex-valued functions.Consequently, we have generalized the results of stability for d＇Alembert＇s and Wilson＇s equations by R. Badora, J. Baker, B. Bouikhalene, P. Gavruta, S. Kabbaj, Pl. Kannappan, G. H.Kim, J.M. Rassias, A. Roukbi, L. Sz′ekelyhidi, D. Zeglami, etc.

Lagrange equations of nonholonomic systems with fractional derivatives

This paper obtains Lagrange equations of nonholonomic systems with fractional derivatives. First, the exchanging relationships between the isochronous variation and the fractional derivatives are derived. Secondly, based on these exchanging relationships, the Hamilton＇s principle is presented for non-conservative systems with fractional derivatives. Thirdly, Lagrange equations of the systems are obtained. Furthermore, the d＇Alembert-Lagrange principle with fractional derivatives is presented, and the Lagrange equations of nonholonomic systems with fractional derivatives are studied. An example is designed to illustrate these results.

THE NAVIER-STOKES EQUATIONS IN STREAM LAYER AND ON STREAM SURFACE AND A DIMENSION SPLIT METHODS

In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.

A correction of the 2D KdV equation of Djordjevic & Redekopp in exponentially stratified fluid

To further study the nal solitons on the continental fission laws of initial intershelf/slope, we rederive and correct the 2D KdV equation of Djordjevic ＆ Redekopp for exponentially stratified fluid （or ocean） and with twodimensional topography. Through a combination of theoretical study and numerical experiments, we show that solitons in the odd vertical modes can fission. However, because of the corrections, the fission conditions are different from those of Djordjevic ＆ Redekopp. The even modes cannot fission unless the initial internal solitons propagate from shallow sea to deep sea. This conclusion is entirely opposite to that of Djordjevic ＆ Redekopp.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

Asymptotic behavior of 2D generalized stochastic Ginzburg-Landau equation with additive noise

The 2D generalized stochastic Ginzburg-Landau equation with additive noise is considered. The compactness of the random dynamical system is established with a priori estimate method, showing that the random dynamical system possesses a random attractor in H^1 0.

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.

One-dimensional dynamic equations of a piezoelectric semiconductor beam with a rectangular cross section and their application in static and dynamic characteristic analysis

Within the framework of continuum mechanics, the double power series ex- pansion technique is proposed, and a series of reduced one-dimensional （1D） equations for a piezoelectric semiconductor beam are obtained. These derived equations are universal, in which extension, flexure, and shear deformations are all included, and can be degen- erated to a number of special cases, e.g., extensional motion, coupled extensional and flexural motion with shear deformations, and elementary flexural motion without shear deformations. As a typical application, the extensional motion of a ZnO beam is analyzed sequentially. It is revealed that semi-conduction has a great effect on the performance of the piezoelectric semiconductor beam, including static deformations and dynamic be- haviors. A larger initial carrier density will evidently lead to a lower resonant frequency and a smaller displacement response, which is a little similar to the dissipative effect. Both the derived approximate equations and the corresponding qualitative analysis are general and widely applicable, which can clearly interpret the inner physical mechanism of the semiconductor in the piezoelectrics and provide theoretical guidance for further experimental design.

A Variational Approach for Numerically Solving the Two-Component Radial Dirac Equation for One-Particle Systems

In this paper we propose a numerical approach to solve the relativistic Dirac equation suitable for computational calculations of one-electron systems. A variational procedure is carried out similar to the well-known Hylleraas computational method. An application of the method to hydrogen isoelectronic atoms is presented, showing its consistency and high accuracy, relative to the exact analytical eigenvalues.

SOLUTIONS AND STABILITY OF A GENERALIZATION OF WILSON'S EQUATION

In this paper we study the solutions and stability of the generalized Wilson＇s functional equation fc f（xty）dtt（t） ＋ fc f（xtσ（y））dtt（t） =2f（x）g（y）, x,y C G, where G is a locally compact group, a is a continuous involution of G and # is an idempotent complex measure with compact support and which is a-invariant. We show that ∫Gg（xty）dp（t） ＋ fcg（xta（y））dp（t） = 2g（x）g（y） if f = 0 and fcf（t.）dp（t） =0, where [fcf（t.）dp（t）]（x） = fc f（tx）dμ（t）. We also study some stability theorems of that equation and we establish the stability on noncommutative groups of the classical Wilson＇s functional equation f（xy） ＋ X（y）f（xa（y）） = 2f（x）g（y） x, y C G, where X is a unitary character of G.

Exact Solutions of D-Dimensional Schr(o|¨)dinger Equation for an Energy-Dependent Potential by NU Method

We study the D-dimensional Schrodinger equation for an energy-dependent Hamiltonian that linearlydepends on energy and quadraticJy on the relative distance.Next,via the Nikiforov-Uvarov (NU) method,we calculatethe corresponding eigenfunctions and eigenvalues.

Computation of the Eigenvalues of 3<i>D</i>“Charged” Integral Equations

The Rayleigh-Ritz and the inverse iteration methods are used in order to compute the eigenvalues of 3D Fredholm-Stieltjes integral equations, i.e. 3D Fredholm equations with respect to suitable Stieltjes-type measures. Some applications are shown, relevant to the problem of computing the eigenvalues of a body charged by a finite number of masses concentrated on points, curves or surfaces lying in.

Two-dimensional equations for thin-films of ionic conductors

A theoretical model is developed for predicting both conduction and diffusion in thin-film ionic conductors or cables. With the linearized Poisson-Nernst-Planck(PNP)theory, the two-dimensional(2D) equations for thin ionic conductor films are obtained from the three-dimensional(3D) equations by power series expansions in the film thickness coordinate, retaining the lower-order equations. The thin-film equations for ionic conductors are combined with similar equations for one thin dielectric film to derive the 2D equations of thin sandwich films composed of a dielectric layer and two ionic conductor layers. A sandwich film in the literature, as an ionic cable, is analyzed as an example of the equations obtained in this paper. The numerical results show the effect of diffusion in addition to the conduction treated in the literature. The obtained theoretical model including both conduction and diffusion phenomena can be used to investigate the performance of ionic-conductor devices with any frequency.

2. 5-D RESISTIVITY TOMOGRAPHY USINGBOUNDARY INTEGRAL EQUATIONS

DC Resistivity Tomography is a non-linear inversion problem. So far there are mainly two kinds of inversion methods, based on the finite-element method and alpha centers method. In this paper, the disadvantages of these two kinds of methods were analysed,and a new method of forward modeling and inversion (Tomography) based on boundary integral equations was proposed. This scheme successfuly overcomes the difficulties of the two formarly methods. It isn’t necessary to use the linearization approximation and calculate the Jacobi matrix. Numerical modeling results given in this paper showed that the computation speed of our method is fast, and there is not any special requirement for initial model, and satisfying results of tomography can be obtained in the case of great contrast of conductivity. So it has wide applications.