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.

A Numerical Study of Riemann Problem Solutions for the Homogeneous One-Dimensional Shallow Water Equations

The solution of the Riemann Problem (RP) for the one-dimensional (1D) non-linear Shallow Water Equations (SWEs) is known to produce four potential wave patterns for the scenario where the water depth is always positive. In this paper, we choose four test problems with exact solutions for the 1D SWEs. Each test problem is a RP with one of the four possible wave patterns as its solution. These problems are numerically solved using schemes from the family of Weighted Essentially Non-Oscillatory (WENO) methods. For comparison purposes, we also include results obtained from the Random Choice Method (RCM). This study has three main objectives. Firstly, we outline the procedures for the implementation of the methods employed in this paper. Secondly, we assess the performance of the schemes in conjunction with a second-order Total Variation Diminishing (TVD) flux on a variety of RPs for the 1D SWEs (for both short- and long-time simulations). Thirdly, we investigate if a single method yields optimal outcomes for all test problems. Optimal outcomes refer to numerical solutions devoid of spurious oscillations, exhibiting high resolution of discontinuities, and attaining high-order accuracy in the smooth parts of the solution.

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.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

The Global Uniqueness of Solutions for a Class of Inverse Problem in 1-D Wave Equations of Hyperbolic Type

This paper has given the global uniquene ss theory of solutions for a class of inverse problem in 1_D Wave equation of hype rbolic type.

New algorithm for solving 3D incompressible viscous equations based on projection method

A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.

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.

A novel(2+1)-dimensional integrable KdV equation with peculiar solution structures

The celebrated(1+1)-dimensional Korteweg de-Vries(KdV)equation and its(2+1)-dimensional extension,the Kadomtsev-Petviashvili(KP)equation,are two of the most important models in physical science.The KP hierarchy is explicitly written out by means of the linearized operator of the KP equation.A novel(2+1)-dimensional KdV extension,the cKP3-4 equation,is obtained by combining the third member(KP3,the usual KP equation)and the fourth member(KP4)of the KP hierarchy.The integrability of the cKP3-4 equation is guaranteed by the existence of the Lax pair and dual Lax pair.The cKP3-4 system can be bilinearized by using Hirota's bilinear operators after introducing an additional auxiliary variable.Exact solutions of the cKP3-4 equation possess some peculiar and interesting properties which are not valid for the KP3 and KP4 equations.For instance,the soliton molecules and the missing D'Alembert type solutions(the arbitrary travelling waves moving in one direction with a fixed model dependent velocity)including periodic kink molecules,periodic kink-antikink molecules,few-cycle solitons,and envelope solitons exist for the cKP3-4 equation but not for the separated KP3 equation and the KP4 equation.

Generalized Wronskian Solutions to Differential-Difference KP Equation

A new method for constructing the Wronskian entries is proposed and applied to the differential-difference Kadomtsev-Petviashvilli （DΔKP） equation. The generalized Wronskian solutions to it are obtained, including rational solutions and Matveev solutions.

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.

Direct Reduction and Exact Solutions for Generalized Variable Coefficients 2D KdV Equation under Some Integrability Conditions

Based on the closed connections among the homogeneous balance （HB） method and Clarkson-KruSkal （CK） method, we study the similarity reductions of the generalized variable coefficients 2D KdV equation. In the meantime it is shown that this leads to a direct reduction in the form of ordinary differential equation under some integrability conditions between the variable coefficients. Two different cases have been discussed, the search for solutions of those ordinary differential equations yielded many exact travelling and solitonic wave solutions in the form of hyperbolic and trigonometric functions under some constraints between the variable coefficients.

Two-Dimensional Rossby Waves： Exact Solutions to Petviashvili Equation

The two-dimensional （2D） nonlinear Rossby waves described by the Petviashvili equation, which has been invoked as an ageostrophic extension of the barotropic quasi-geostrophic potential vorticity equation, can be investigated through the exact periodic-wave solutions for the Petviashvili equation, while the exact analytical periodic-wave solutions to the Petviashvili equation are obtained by using the Jacobi elliptic function expansion method. It is shown that periodicwave 2D Rossby solutions can be obtained by this method, and in the limit cases, the 213 Rossby soliton solutions are also obtained.

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.

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.

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.