The practical Gauss type rules for Hadamard finite-part integrals using Puiseux expansions

A general framework is constructed for efficiently and stably evaluating the Hadamard finite-part integrals by composite quadrature rules. Firstly, the integrands are assumed to have the Puiseux expansions at the endpoints with arbitrary algebraic and logarithmic singularities. Secondly, the Euler-Maclaurin expansion of a general composite quadrature rule is obtained directly by using the asymptotic expansions of the partial sums of the Hurwitz zeta function and the generalized Stieltjes constant, which shows that the standard numerical integration formula is not convergent for computing the Hadamard finite-part integrals. Thirdly, the standard quadrature formula is recast in two steps. In step one, the singular part of the integrand is integrated analytically and in step two, the regular integral of the remaining part is evaluated using the standard composite quadrature rule. In this stage, a threshold is introduced such that the function evaluations in the vicinity of the singularity are intentionally excluded, where the threshold is determined by analyzing the roundoff errors caused by the singular nature of the integrand. Fourthly, two practical algorithms are designed for evaluating the Hadamard finite-part integrals by applying the Gauss-Legendre and Gauss-Kronrod rules to the proposed framework. Practical error indicator and implementation involved in the Gauss-Legendre rule are addressed. Finally, some typical examples are provided to show that the algorithms can be used to effectively evaluate the Hadamard finite-part integrals over finite or infinite intervals.     

第二类端点奇异Fredholm积分方程的分数阶退化核方法

本文针对第二类端点奇异Fredholm积分方程构造基于分数阶Taylor展开的退化核方法，设计了两种计算格式，一是在全区间上使用分数阶Taylor展开式近似核函数，二是在包含奇点的小区间上采用分数阶插值，在剩余区间上采用分段二次多项式插值逼近核函数.讨论了两种退化核方法收敛的条件，并给出了混合插值法的收敛阶估计.数值算例表明对于非光滑核函数分数阶退化核方法有着良好的计算效果，且混合二次插值法比全区间上的分数阶退化核方法有着更广泛的适用范围.     

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＦｉｎｉｔｅｖｏｌｕｍｅｅｌｅｍｅｎｔｍｅｔｈｏｄｏｒＦＶＥＭｕｓｅｓａｖｏｌｕｍｅｉｎｔｅｇｒａｌｆｏｒｍｕｌａｔｉｏｎｏｆｔｈｅｄｉｆ ｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｗｉｔｈａｆｉｎｉｔｅｐａｒｔｉｔｉｏｎｉｎｇｓｅｔｏｆｖｏｌｕｍｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｅｑｕａｔｉｏｎ ,ｔｈｅｎｒｅ ｓｔｒｉｃｔｓｔｈｅａｄｍｉｓｓｉｂｌｅｆｕｎｃｔｉｏｎｓｔｏａｌｉｎｅａｒｆｉｎｉｔｅｅｌｅｍｅｎｔｓｐａｃｅｔｏｄｉｓｃｒｅｔｉｚｅｔｈｅｓｏｌｕｔｉｏｎ ([1 -4 ] ) .…     

Cubic superconvergent finite volume element method for one-dimensional elliptic and parabolic equations

In this paper, a cubic superconvergent finite volume element method based on optimal stress points is presented for one-dimensional elliptic and parabolic equations. For elliptic problem, it is proved that the method has optimal third order accuracy with respect to H¹ norm and fourth order accuracy with respect to L² norm. We also obtain that the scheme has fourth order superconvergence for derivatives at optimal stress points. For parabolic problem, the scheme is given and error estimate is obtained with respect to L² norm. Finally, numerical examples are provided to show the effectiveness of the method.     

抛物型方程非齐次边值问题的推广型LOD有限差分及有限元格式

1引言本文考虑区域Ω=[0,1]~d(d=2,3)上的非齐次抛物型方程第一边值问题(?)-C_1△u C_2u=f(x,t),x∈Ω,t∈(0,T],(1．1) u(x,0)=u_0(x),x∈Ω,(1．2) u(x,t)=(?)(x,t),x∈(?)Ω,t∈(0,T],(1．3)其中C_1,C_2为常数且C_1＞0,C_2≥0．对于以上问题,可以使用有限差分方法及有限元方法进行离散,并采用交替方向方法求解．交替方向方法能够将高维问题转化为一系列的一维问题进行计算,具有计算量少,计算稳定且易于并行实现等优点,在大规模科学计算中起着非常重要的作用,一直是计算数     

Alternating direction finite volume element methods for 2D parabolic partial differential equations

On the basis of rectangular partition and bilinear interpolation, this article presents alternating direction finite volume element methods for two dimensional parabolic partial differential equations and gives three computational schemes, one is analogous to Douglas finite difference scheme with second order splitting error, the second has third order splitting error, and the third is an extended locally one dimensional scheme. Optimal L² norm or H¹ semi‐norm error estimates are obtained for these schemes. Finally, two numerical examples illustrate the effectiveness of the schemes. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

大鼠单侧睾丸扭转对未扭转侧睾丸的影响及红花的保护作用

为探讨大鼠单侧睾丸扭转对未扭转侧睾丸生精功能的影响及红花对未扭转侧睾丸早期损伤的保护作用,将21只健康雄性青春期前SD大鼠随机均分为假手术对照组(A组,n=7)、睾丸扭转复位组(B组,n=7)和睾丸扭转复位+红花组(C组,n=7).按Turner法建立睾丸扭转模型,喂养至术后7 d时处死,即刻采集各组大鼠未扭转侧睾丸标...     

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

<正>This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes,one is analogous to Douglas finite difference scheme with second-order splitting error,the other two schemes have third-order splitting error,and the last one is an extended LOD scheme.The L~2 norm and H~1 semi-norm error estimates are obtained for the first scheme and second one,respectively.Finally,two numerical examples are provided to illustrate the efficiency and accuracy of the methods.     

一维二阶椭圆型方程组的超收敛二次有限体积元方法

本文针对二阶椭圆型常微分方程组边值问题提出二次超收敛有限体积元方法,证明格式的H1和L2模误差估计,并给出应力佳点处的梯度超收敛估计.最后,编写计算格式的Fortran程序,用数值算例验证了理论分析的正确性和格式的有效性.     

铜中毒对SD大鼠肝的影响研究

将36只雄性大鼠80~100 g随机分为3组,分别喂以对照日粮(Cu质量分数13.15mg/kg)和高铜日粮(Cu质量分数650、850 mg/kg),实验期为9周,进行了铜中毒对大鼠肝脏影响的病理学研究。结果表明,与对照组比较,高铜Ⅰ、Ⅱ组大鼠肝脏出现不同程度的病理损害;大鼠肝脏脏器系数有增高趋势,但差异无显著(P0.05)。同时,乳酸脱氢酶、谷草转氨酶活性显著升高(P0.01),血清铜含量上升,SOD含量下降。提示日粮铜质量分数650~850 mg/kg即可引起大鼠肝脏的病理损伤,功能降低。