Generalized Bessel functions: Theory and their applications

This paper presents 2 new classes of the Bessel functions on a compact domain [0,T] as generalized‐tempered Bessel functions of the first‐ and second‐kind which are denoted by GTBFs‐1 and GTBFs‐2. Two special cases corresponding to the GTBFs‐1 and GTBFs‐2 are considered. We first prove that these functions are as the solutions of 2 linear differential operators and then show that these operators are self‐adjoint on suitable domains. Some interesting properties of these sets of functions such as orthogonality, completeness, fractional derivatives and integrals, recursive relations, asymptotic formulas, and so on are proved in detail. Finally, these functions are performed to approximate some functions and also to solve 3 practical differential equations of fractionalorders.     

高速移动荷载下黏弹性半空间体的动力响应

分别以移动荷载和黏弹性半空间体模拟运动列车荷载和地基,分析了地基在运动列车作用下的动力响应.首先采用Green函数法求解黏弹性半空间体在各种移动荷载模式作用下的动力响应的解析解,包括恒常和简谐移动点源、线源和面源荷载.然后采用IFFT算法和自适应数值积分算法计算解析解中的二维积分,得到了包括低音速、跨音速和超音速移动荷载作用下位移的数值结果.最后分析了速度对位移的分布和最大值的影响,发现当速度大于Rayleigh波速时,位移发生显著变化.     

Recursion relations for the extended Krylov subspace method

The evaluation of matrix functions of the form f(A)v, where A is a large sparse or structured symmetric matrix, f is a nonlinear function, and v is a vector, is frequently subdivided into two steps: first an orthonormal basis of an extended Krylov subspace of fairly small dimension is determined, and then a projection onto this subspace is evaluated by a method designed for small problems. This paper derives short recursion relations for orthonormal bases of extended Krylov subspaces of the type K^m,mi+1(A)=span{A^-m+1v,…,A^-1v,v,Av,…,A^miv}, m=1,2,3,…, with i a positive integer, and describes applications to the evaluation of matrix functions and the computation of rational Gauss quadrature rules.     

Numerical solution of two-dimensional reaction-diffusion Brusselator system

In this paper, polynomial based differential quadrature method (DQM) is applied for the numerical solution of a class of two-dimensional initial-boundary value problems governed by a non-linear system of partial differential equations. The system is known as the reaction-diffusion Brusselator system. The system arises in the modeling of certain chemical reaction-diffusion processes. In Brusselator system the reaction terms arise from the mathematical modeling of chemical systems such as in enzymatic reactions, and in plasma and laser physics in multiple coupling between modes. The numerical results reported for three specific problems. Convergence and stability of the method is also examined numerically.     

On a rational differential quadrature method in irregular domains for problems with boundary layers

A rational differential quadrature method in irregular domains (RDQMID) is investigated to deal with a kind of singularly perturbed problems with boundary layers. Through a transformation, the boundary layer, which may be not straight, is transformed into a segment of a line parallel to one of the Cartesian axes. The rational differential quadrature method (RDQM) is applied to discretize the governing equation. Finally, a direct expansion method of the boundary conditions (DEMBC) is raised to deal with the boundary conditions. Numerical experiments show that RDQMID is of high accuracy, efficiency and easy to programme.     

A note on one-point quadrature formula for Daubechies scale function with partial support

This note is concerned with an efficient computation of integrals of products of a smooth function and Daubechies scale function with partial support by using a one-point quadrature rule. The error estimate is obtained. The rule is illustrated by considering an example from the literature.     

Application of Gaussian quadrature method to characterize heavy ends of hydrocarbon fluids for modeling wax precipitation

The hydrocarbon plus fractions that comprise a significant portion of naturally occurring hydrocarbon fluids create major problems when determining the thermodynamic properties and the volumetric behavior of these fluids by equations of state. These problems arise due to the difficulty of properly characterizing the plus fractions (heavy ends). Proper characterization of the heavier components is important when cubic equations of state and/or solid formation thermodynamic models are used to describe complex phase behavior of reservoir fluids. The effect of heavy fractions characterization on thermodynamic modeling of wax precipitation has been investigated using different models including Won, Pan and proposed models. In order to characterize the plus fraction (heavier part) as a series of pseudocomponents, a probability model that expresses the mole fraction as a continuous function of the molecular weight has been used. The study has been conducted using several mixtures. Two different SCN (single carbon number), C₇₊${\text{C}}_{7^{+}}$ and C₁₀₊${\text{C}}_{{10}^{+}}$ were chosen. The Chosen SCN were distributed to multi-components of five, six, and/or ten using continuous function and Gaussian quadrature method. The results showed that the fractioning is required to be able to predict wax precipitation. Distribution of C₁₀₊${\text{C}}_{{10}^{+}}$ using a proper distribution function has shown improvement in predictions of WAT and the amount of wax deposited in comparison with the characterization of C₇₊${\text{C}}_{7^{+}}$ using semi-continuous approach. In predicting of WAT and the amount of wax build up the developed model showed superiority over the others.     

Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations

In this paper, we state and prove a new formula expressing explicitly the derivatives of shifted Chebyshev polynomials of any degree and for any fractional-order in terms of shifted Chebyshev polynomials themselves. We develop also a direct solution technique for solving the linear multi-order fractional differential equations (FDEs) with constant coefficients using a spectral tau method. The spatial approximation with its fractional-order derivatives (described in the Caputo sense) are based on shifted Chebyshev polynomials T_L,n(x) with x ∈ (0, L), L > 0 and n is the polynomial degree. We presented a shifted Chebyshev collocation method with shifted Chebyshev–Gauss points used as collocation nodes for solving nonlinear multi-order fractional initial value problems. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed techniques and to compare with the existing results.     

Efficient long-time computations of time-domain boundary integrals for 2D and dissipative wave equation

Linear hyperbolic partial differential equations in a homogeneous medium, e.g., the wave equation describing the propagation and scattering of acoustic waves, can be reformulated as time-domain boundary integral equations. We propose an efficient implementation of a numerical discretization of such equations when the strong Huygens’ principle does not hold.For the numerical discretization, we make use of convolution quadrature in time and standard Galerkin boundary element method in space. The quadrature in time results in a discrete convolution of weights W_j with the boundary density evaluated at equally spaced time points. If the strong Huygens’ principle holds, W_j converge to 0 exponentially quickly for large enough j. If the strong Huygens’ principle does not hold, e.g., in even space dimensions or when some damping is present, the weights are never zero, thereby presenting a difficulty for efficient numerical computation.In this paper we prove that the kernels of the convolution weights approximate in a certain sense the time domain fundamental solution and that the same holds if both are differentiated in space. The tails of the fundamental solution being very smooth, this implies that the tails of the weights are smooth and can efficiently be interpolated. Further, we hint on the possibility to apply the fast and oblivious convolution quadrature algorithm of Schädle et al. to further reduce memory requirements for long-time computation. We discuss the efficient implementation of the whole numerical scheme and present numerical experiments.