Efficient simulation of spatially correlated non-stationary ground motions by wavelet-packet algorithm and spectral representation method

Although the classical spectral representation method(SRM)has been widely used in the generation of spatially varying ground motions,there are still challenges in efficient simulation of the non-stationary stochastic vector process in practice.The first problem is the inherent limitation and inflexibility of the deterministic time/frequency modulation function.Another difficulty is the estimation of evolutionary power spectral density(EPSD)with quite a few samples.To tackle these problems,the wavelet packet transform(WPT)algorithm is utilized to build a time-varying spectrum of seed recording which describes the energy distribution in the time-frequency domain.The time-varying spectrum is proven to preserve the time and frequency marginal property as theoretical EPSD will do for the stationary process.For the simulation of spatially varying ground motions,the auto-EPSD for all locations is directly estimated using the time-varying spectrum of seed recording rather than matching predefined EPSD models.Then the constructed spectral matrix is incorporated in SRM to simulate spatially varying non-stationary ground motions using efficient Cholesky decomposition techniques.In addition to a good match with the target coherency model,two numerical examples indicate that the generated time histories retain the physical properties of the prescribed seed recording,including waveform,temporal/spectral non-stationarity,normalized energy buildup,and significant duration.

THE SPACE-TIME FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS

Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.

Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.

THE QUASI-BOUNDARY VALUE METHOD FOR IDENTIFYING THE INITIAL VALUE OF THE SPACE-TIME FRACTIONAL DIFFUSION EQUATION

In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.

An accurate and efficient space-time Galerkin spectral method for the subdiffusion equation

In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the initial time. Consequently, traditional high-order numerical methods in time are inefficient. The second obstacle is that the resulting system of the space-time spectral approach is usually large and time-consuming to solve. We aim at overcoming the first difficulty by proposing a novel approach in time, which is based on variable transformation techniques. Suitable ψ-fractional Sobolev spaces and a new variational framework are introduced to establish the well-posedness of the associated variational problem. This allows us to construct our space-time spectral method using a combination of temporal generalized Jacobi polynomials(GJPs) and spatial Legendre polynomials. For the second difficulty, we propose a fast algorithm to effectively solve the resulting linear system. The fast algorithm makes use of a matrix diagonalization in space and QZ decomposition in time. Our analysis and numerical experiments show that the proposed method is exponentially convergent with respect to the polynomial degrees in both space and time directions, even though the exact solution has very limited regularity.

Generation of endurance time excitation functions using spectral representation method

In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

H^1 space-time discontinuous finite element method for convection-diffusion equations

An H1 space-time discontinuous Galerkin （STDG） scheme for convection- diffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L∞ （H1） norm is derived. The numerical exper- iments are presented to verify the theoretical results.

An Eigenspace Method for Detecting Space-Time Disease Clusters with Unknown Population-Data

Space-time disease cluster detection assists in conducting disease surveillance and implementing control strategies.The state-of-the-art method for this kind of problem is the Space-time Scan Statistics(SaTScan)which has limitations for non-traditional/non-clinical data sources due to its parametric model assumptions such as Poisson orGaussian counts.Addressing this problem,an Eigenspace-based method called Multi-EigenSpot has recently been proposed as a nonparametric solution.However,it is based on the population counts data which are not always available in the least developed countries.In addition,the population counts are difficult to approximate for some surveillance data such as emergency department visits and over-the-counter drug sales,where the catchment area for each hospital/pharmacy is undefined.We extend the population-based Multi-EigenSpot method to approximate the potential disease clusters from the observed/reported disease counts only with no need for the population counts.The proposed adaptation uses an estimator of expected disease count that does not depend on the population counts.The proposed method was evaluated on the real-world dataset and the results were compared with the population-based methods:Multi-EigenSpot and SaTScan.The result shows that the proposed adaptation is effective in approximating the important outputs of the population-based methods.

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

SPACE-TIME DISCONTINUOUS GALERKIN METHOD FOR MAXWELL EQUATIONS IN DISPERSIVE MEDIA

In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability and error estimate of order O？τr＋1＋hk＋1/2？ are obtained when polynomials of degree at most r and k are used for the temporal dis-cretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r＋1 in temporal variable t.

A Space-Time Interior Penalty Discontinuous Galerkin Method for the Wave Equation

A new higher-order accurate space-time discontinuous Galerkin(DG)method using the interior penalty flux and discontinuous basis functions,both in space and in time,is pre-sented and fully analyzed for the second-order scalar wave equation.Special attention is given to the definition of the numerical fluxes since they are crucial for the stability and accuracy of the space-time DG method.The theoretical analysis shows that the DG discre-tization is stable and converges in a DG-norm on general unstructured and locally refined meshes,including local refinement in time.The space-time interior penalty DG discre-tization does not have a CFL-type restriction for stability.Optimal order of accuracy is obtained in the DG-norm if the mesh size h and the time stepΔt satisfy h≅CΔt,with C a positive constant.The optimal order of accuracy of the space-time DG discretization in the DG-norm is confirmed by calculations on several model problems.These calculations also show that for pth-order tensor product basis functions the convergence rate in the L∞and L2-norms is order p+1 for polynomial orders p=1 and p=3 and order p for polynomial order p=2.

PRECONDITIONED METHODS FOR SPACE-TIME ADAPTIVE PROCESSING

This paper introduces the preconditioned methods for Space-Time Adaptive Processing(STAP).Using the Block-Toeplitz-Toeplitz-Block(BTTB)structure of the clutter-plus-noise covari-ance matrix,a Block-Circulant-Circulant-Block(BCCB)preconditioner is constructed.Based on thepreconditioner,a Preconditioned Multistage Wiener Filter(PMWF)which can be implemented by thePreconditioned Conjugate Gradient(PCG)method is proposed.Simulation results show that thePMWF has faster convergence rate and lower processing rank compared with the MWF.

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.The staggered DG scheme defines the discrete pressure on the primal triangular mesh,while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh.In this paper,a new pair of equal-order-interpolation velocity-pressure finite elements is proposed.On the primary triangular mesh(the pressure elements),the basis functions are piecewise polynomials of degree N and are allowed to jump on the boundaries of each triangle.On the dual mesh instead(the velocity elements),the basis functions consist in the union of piecewise polynomials of degree N on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries,while they are continuous inside.In other words,the basis functions on the dual mesh arc built by continuous finite elements on the subtriangles.This choice allows the construction of an efficient,quadrature-free and memory saving algorithm.In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations,the arbitrary high order of accuracy in time is achieved through tire use of time-dependent test and basis functions,in combination with simple and efficient Picard iterations.Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes,not only from a computer memory point of view,but also concerning the computational time.

Focused Wave Properties Based on A High Order Spectral Method with A Non-Periodic Boundary

In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.

Three-dimensional simulations of strong ground motion in the Shidian basin based upon the spectral-element method

The strong motion of a small long and narrow basin caused by a moderate scenario earthquake is simulated by using the spectral-element method and the parallel computing technique.A total of five different geometrical profiles within the basin are used to analyze the generation and propagation of surface waves and their relation to the basin structures in both the time and frequency domain.The amplification effects are analyzed by the distribution of peak ground velocity(PGV)and cumulative kinetic energy(Ek) in the basin.The results show that in the 3D basin,the excitation of the fundamental and higher surface wave modes are similar to that of the 2D model.Small bowls in the basin have great influence on the amplification and distribution of strong ground motion,due to their lateral resonances when the wavelengths of the lateral surface waves are comparable to the size of the bowls.Obvious basin edge effects can be seen at the basin edge closer to the source for constructive interference between direct body waves and the basin-induced surface waves.The Ek distribution maps show very large values in small bowls and some corners in the basin due to the interference of waves propagating in different directions.A high impedance contrast model can excite more surface wave modes,resulting in longer shaking durations as well as more complex seismograms and PGV and Ek distributions.

CONVERGENCE ANALYSIS OF THE JACOBI SPECTRAL-COLLOCATION METHOD FOR FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L^2-norm. The numerical examples are given to illustrate the theoretical results.

COMPACT FINITE DIFFERENCE-FOURIER SPECTRAL METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Q-factor estimation in CMP gather and the continuous spectral ratio slope method

The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.