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.

An Adaptive Spectral Conjugate Gradient Method with Restart Strategy

As a generalization of the two-term conjugate gradient method(CGM),the spectral CGM is one of the effective methods for solving unconstrained optimization.In this paper,we enhance the JJSL conjugate parameter,initially proposed by Jiang et al.(Computational and Applied Mathematics,2021,40:174),through the utilization of a convex combination technique.And this improvement allows for an adaptive search direction by integrating a newly constructed spectral gradient-type restart strategy.Then,we develop a new spectral CGM by employing an inexact line search to determine the step size.With the application of the weak Wolfe line search,we establish the sufficient descent property of the proposed search direction.Moreover,under general assumptions,including the employment of the strong Wolfe line search for step size calculation,we demonstrate the global convergence of our new algorithm.Finally,the given unconstrained optimization test results show that the new algorithm is effective.

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.

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.

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.

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.

High-precision solution to the moving load problem using an improved spectral element method

In this paper, the spectral element method(SEM)is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem.In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases.Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.

Re-study on Recurrence Period of Stokes Wave Train with High Order Spectral Method

Owing to the Benjamin-Feir instability, the Stokes wave train experiences a modulation-demodulation process, and presents a recurrence characteristics. Stiassnie and Shemer researched the unstable evolution process and provided a theoretical formulation for the recurrence period in 1985 on the basis of the nonlinear cubic Schrodinger equation （NLS）. However, NLS has limitations on the narrow band and the weak nonlinearity. The recurrence period is re-investigated in this paper by using a highly efficient High Order Spectral （HOS） method, which can be applied for the direct phase- resolved simulation of the nonlinear wave train evolution. It is found that the Stiassnie and Shemer＇s formula should be modified in the cases with most unstable initial conditions, which is important for such topics as the generation mechanisms of freak waves. A new recurrence period formula is presented and some new evolution characteristics of the Stokes wave train are also discussed in details.

An efficient source wavefield reconstruction scheme using single boundary layer values for the spectral element method

In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.

THE DYNAMICAL BEHAVIOR OF FULLY DISCRETE SPECTRAL METHOD FOR NONLINEAR SCHRODINGER EQUATION WITH WEAKLY DAMPED

Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to

Dynamic analysis of beam-cable coupled systems using Chebyshev spectral element method

The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of highorder accuracy. To attain high accuracy in time discretization, the fourth-order AdamsBashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Hermite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation （nonlinear convection- diffusion problem） and KdV equation（single solitary and 2-solitary wave problems）, where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.

Cross-spectral recognition method of bridge deck aerodynamic admittance function

This study proposes a new identification algorithm about the admittance function, which can estimate the full set of six aerodynamic admittance functions considering cross power spectral density functions about the forces and the turbulence components. The method was first numerically validated through Monte Carlo simulations, and then adopted to estimate the aerodynamic admittance of a streamlined bridge deck. The identification method was further validated through a comparison between the numerical calculation and wind tunnel tests on a moving bridge section. © 2015, Institute of Engineering Mechanics, China Earthquake Administration and Springer-Verlag Berlin Heidelberg.

Seismic wave modeling in viscoelastic VTI media using spectral element method

Spectral element method（SEM） for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity of earth media,SEM for elastic media is no longer appropriate.On fundamental of the second-order elastic SEM,this work takes the viscoelastic wave equations and the vertical transversely isotropic（VTI） media into consideration,and establishes the second-order SEM for wave modeling in viscoelastic VTI media.The second-order perfectly matched layer for viscoelastic VTI media is also introduced.The problem of handling the overlapped absorbed corners is solved.A comparison with the analytical solution in a twodimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling.Furtherly,numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media.This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.

SPECTRAL METHOD IN TIME FOR KdV EQUATIONS

This paper presents a fully spectral discretization method for solving KdV equations with periodic boundary conditions.Chebyshev pseudospectral approximation in the time direction and Fourier Galerkin approximation in the spatial direction.The expansion coefficients are determined by minimizing an object funictional.Rapid convergence of the method is proved.

A lumped mass Chebyshev spectral element method and its application to structural dynamic problems

A diagonal or lumped mass matrix is of great value for time-domain analysis of structural dynamic and wave propagation problems,as the computational efforts can be greatly reduced in the process of mass matrix inversion.In this study,the nodal quadrature method is employed to construct a lumped mass matrix for the Chebyshev spectral element method(CSEM).A Gauss-Lobatto type quadrature,based on Gauss-Lobatto-Chebyshev points with a weighting function of unity,is thus derived.With the aid of this quadrature,the CSEM can take advantage of explicit time-marching schemes and provide an efficient new tool for solving structural dynamic problems.Several types of lumped mass Chebyshev spectral elements are designed,including rod,beam and plate elements.The performance of the developed method is examined via some numerical examples of natural vibration and elastic wave propagation,accompanied by their comparison to that of traditional consistent-mass CSEM or the classical finite element method(FEM).Numerical results indicate that the proposed method displays comparable accuracy as its consistent-mass counterpart,and is more accurate than classical FEM.For the simulation of elastic wave propagation in structures induced by high-frequency loading,this method achieves satisfactory performance in accuracy and efficiency.