In continuum one-dimensional space, a coupled directed continuous time random walk model is proposed, where the random walker jumps toward one direction and the waiting time between jumps affects the subsequent jump. ...
By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in...
详细信息
By the standard theory,the stable Qk+1,k−Qk,k+1/Qdck divergence-free element converges with the optimal order of approximation for the Stokes equations,but only order k for the velocity in H1-norm and the pressure in *** is due to one polynomial degree less in y direction for the first component of velocity,which is a Qk+1,k polynomial of x and *** this manuscript,we will show by supercloseness of the divergence free element that the order of convergence is truly k+1,for both velocity and *** special solutions(if the interpolation is also divergence-free),a two-order supercloseness is shown to *** tests are provided confirming the accuracy of the theory.
The multifractal properties of daily rainfall time series at the stations in Pearl River basin of China over periods of up to 45 years are examined using the universal multifractal approach based on the multiplicative...
详细信息
The multifractal properties of daily rainfall time series at the stations in Pearl River basin of China over periods of up to 45 years are examined using the universal multifractal approach based on the multiplicative cascade model and the multifractal detrended fluctuation analysis (MF-DFA). The results from these two kinds of multifractal analyses show that the daily rainfall time series in this basin have multifractal behavior in two different time scale ranges. It is found that the empirical multifractal moment function K ( q ) of the daily rainfall time series can be fitted very well by the universal multifractal model (UMM). The estimated values of the conservation parameter H from UMM for these daily rainfall data are close to zero indicating that they correspond to conserved fields. After removing the seasonal trend in the rainfall data, the estimated values of the exponent h ( 2 ) from MF-DFA indicate that the daily rainfall time series in Pearl River basin exhibit no long-term correlations. It is also found that K ( 2 ) and elevation series are negatively correlated. It shows a relationship between topography and rainfall variability.
Based on W -transformation, some parametric symplectic partitioned Runge–Kutta (PRK) methods depending on a real parameter α are developed. For α = 0 , the corresponding methods become the usual PRK methods, includ...
Based on W -transformation, some parametric symplectic partitioned Runge–Kutta (PRK) methods depending on a real parameter α are developed. For α = 0 , the corresponding methods become the usual PRK methods, including Radau IA – I A ¯ and Lobatto IIIA – IIIB methods as examples. For any α ≠ 0 , the corresponding methods are symplectic and there exists a value α ∗ such that energy is preserved in the numerical solution at each step. The existence of the parameter and the order of the numerical methods are discussed. Some numerical examples are presented to illustrate these results.
In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caput...
详细信息
In this paper, we apply the Legendre spectral-collocation method to obtain approximate solutions of nonlinear multi-order fractional differential equations (M-FDEs). The fractional derivative is described in the Caputo sense. The study is conducted through illustrative example to demonstrate the validity and applicability of the presented method. The results reveal that the proposed method is very effective and simple. Moreover, only a small number of shifted Legendre polynomials are needed to obtain a satisfactory result.
The Hagedorn wavepacket method is an important numerical method for solving the semiclassical time-dependent Schrödinger equation. In this paper, a new semi-discretization in space is obtained by wavepacket opera...
详细信息
The Hagedorn wavepacket method is an important numerical method for solving the semiclassical time-dependent Schrödinger equation. In this paper, a new semi-discretization in space is obtained by wavepacket operator. In a sense, such semi-discretization is equivalent to the Hagedorn wavepacket method, but this discretization is more intuitive to show the advantages of wavepacket methods. Moreover, we apply the multi-time-step method and the Magnus-expansion to obtain the improved algorithms in time-stepping computation. The improved algorithms are of the Gauss–Hermite spectral accuracy to approximate the analytical solution of the semiclassical Schrödinger equation. And for the given accuracy, the larger time stepsize can be used for the higher oscillation in the semiclassical Schrödinger equation. The superiority is shown by the error estimation and numerical experiments.
In this paper, the fractional variational integrators developed by Wang and Xiao (2012) [28] are extended to the fractional Euler–Lagrange (E–L) equations with holonomic constraints. The corresponding fractional dis...
详细信息
In this paper, the fractional variational integrators developed by Wang and Xiao (2012) [28] are extended to the fractional Euler–Lagrange (E–L) equations with holonomic constraints. The corresponding fractional discrete E–L equations are derived, and their local convergence is discussed. Some fractional variational integrators are presented. The suggested methods are shown to be efficient by some numerical examples.
Anisotropic meshes are known to be well-suited for problems which exhibit anisotropic solution features. Defining an appropriate metric tensor and designing an efficient algorithm for anisotropic mesh generation are t...
详细信息
Anisotropic meshes are known to be well-suited for problems which exhibit anisotropic solution features. Defining an appropriate metric tensor and designing an efficient algorithm for anisotropic mesh generation are two important aspects of the anisotropic mesh methodology. In this paper, we are concerned with the natural metric tensor for use in anisotropic mesh generation for anisotropic elliptic problems. We provide an algorithm to generate anisotropic meshes under the given metric tensor. We show that the inverse of the anisotropic diffusion matrix of the anisotropic elliptic problem is a natural metric tensor for the anisotropic mesh generation in three aspects: better discrete algebraic systems, more accurate finite element solution and superconvergence on the mesh nodes. Various numerical examples demonstrating the effectiveness are presented.
We will investigate the superconvergence for the semidiscrete finite element approximation of distributed convex optimal control problems governed by semilinear parabolic equations. The state and costate are approxima...
We will investigate the superconvergence for the semidiscrete finite element approximation of distributed convex optimal control problems governed by semilinear parabolic equations. The state and costate are approximated by the piecewise linear functions and the control is approximated by piecewise constant functions. We present the superconvergence analysis for both the control variable and the state variables.
In this paper,the time-dependent Maxwell’s equations used to modeling wave propagation in dispersive lossy bi-isotropic media are *** and uniqueness of the modeling equations are *** fully discrete finite element sch...
详细信息
In this paper,the time-dependent Maxwell’s equations used to modeling wave propagation in dispersive lossy bi-isotropic media are *** and uniqueness of the modeling equations are *** fully discrete finite element schemes are proposed,and their practical implementation and stability are discussed.
暂无评论