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.

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.

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.

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.

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.

Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Abundant solutions of Wick-type stochastic fractional 2D KdV equations

A modified fractional sub-equation method is applied to Wick-type stochastic fractional two-dimensional （2D） KdV equations. With the help of a Hermit transform, we obtain a new set of exact stochastic solutions to Wick-type stochastic fractional 2D KdV equations in the white noise space. These solutions include exponential decay wave solutions, soliton wave solutions, and periodic wave solutions. Two examples are explicitly given to illustrate our approach.

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.

Solving Cauchy Issues of Highly Nonlinear Elliptic Equations Using a Meshless Method

In this paper,we address 3D inverse Cauchy issues of highly nonlinear elliptic equations in large cuboids by utilizing the new 3D homogenization functions of different orders to adapt all the specified boundary data.We also add the average classification as an approximate solution to the nonlinear operator part,without requiring to cope with nonlinear equations to resolve the weighting coefficients because these constructions are owned many conditions about the true solution.The unknown boundary conditions and the result can be retrieved straightway by coping with a small-scale linear system when the outcome is described by a new 3D homogenization function,which is right to find the numerical solutions with the errors smaller than the level of noise being put on the over-specified Neumann conditions on the bottom of the cuboid.Besides,note that the new homogenization functions method(HFM)does not require dealing with the regularization and highly nonlinear equations.The robustness and accuracy of the HFM are verified by comparing the recovered results of several numerical experiments to the exact solutions in the entire region,even though a very large level of noise 50%is imposed on the over specified Neumann conditions.The numerical errors of our scheme are in the order of O(10^(−1))-O(10^(−4)).

Modeling Microbial Decomposition in Real 3D Soil Structures Using Partial Differential Equations

Partial Differential Equations (PDEs) have been already widely used to simulate various complex phenomena in porous media. This paper is one of the first attempts to apply PDEs for simulating in real 3D structures. We apply this scheme to the specific case study of the microbial decomposition of organic matter in soil pore space. We got a 3D geometrical representation of the pore space relating to a network of volume primitives. A mesh of the pore space is then created by using the network. PDEs system is solved by free finite elements solver Freefem3d in the particular mesh. We validate our PDEs model to experimental data with 3D Computed Tomography (CT) images of soil samples. Regarding the current state of art on soil organic matter decay models, our approach allows taking into account precise 3D spatialization of the decomposition process by a pore space geometry description.

Higher-order differential variational principle and differential equations of motion for mechanical systems in event space

In this paper we study the higher-order differential variational principle and differential equations of motion for mechanical systems in event space. Based on the higher-order d＇Alembert principle of the system, the higher-order velocity energy and the higher-order acceleration energy of the system in event space are defined, the higher-order d＇Alembert- Lagrange principle of the system in event space is established, and the parametric forms of Euler-Lagrange, Nielsen and Appell for this principle are given. Finally, the higher-order differential equations of motion for holonomic systems in event space are obtained.

A Local Discontinuous Galerkin Method for Two-Dimensional Time Fractional Diffusion Equations

For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.