B-spline finite-element method in polar coordinates

A numerical solution method for two-dimensional electromagnetic field problems is presented using the B-spline finite-element expression based on polar coordinates. The technique has two main advantages: (1) to avoid the truncation errors at some curved boundaries and (2) to improve the accuracy of singular boundary-value problems with a sharp corner. The B-spline finite-element formulation in polar coordinates is derived and its numerical applications are illustrated by an eddy current problem and several waveguide eigenvalue problems.     

对称荷载作用下半圆盘的应力分析

本文应用双极坐标系中的弹性力学相容方程,获得了受对称荷载作用且τξη=0情况下半圆盘的应力解答。运用所得公式,分析了预应力钢丝缠绕式压机中半圆盘形梁等实例。     

GEOMETRIC ANALYSIS OF FINITE STRAIN TENSORS IN TECTONODEFORMATION MECHANICS

According to the peculiarities of tectonodeformation this paper systematically researches the graphic method of finite strain theory in tectonomechanics, proposes a geometrical analytic method describing the strain state of a point with metric tensor, expounds and proves the structure of deformation gradient circle, and clearly illustrates the relations of various quantities in polar decomposition theory. Finally, a simple shear deformation which often occurs in tectonodeformation was calculated through various different methods in order to explain and to compare the peculiarities among various graphic methods of finite strain.     

Combined perturbation bounds:Ⅱ.Polar decompositions

In this paper,we study the perturbation bounds for the polar decomposition A=QH where Q is unitary and H is Hermitian.The optimal (asymptotic) bounds obtained in previous works for the unitary factor,the Hermitian factor and singular values of A areσ_r~2||ΔQ||_F~2≤||ΔA||_F~2, 1/2||ΔH||_F~2≤||ΔA||_F~2 and ||Δ∑||_F~2≤||ΔA||_F~2,respectively,where∑=diag(σ_1,σ_2,...,σ_r,0,...,0) is the singular value matrix of A andσ_r denotes the smallest nonzero singular value.Here we present some new combined (asymptotic) perturbation boundsσ_r~2||ΔQ||_F~2 1/2||ΔH||_F~2≤||ΔA||_F~2 andσ_r~2||ΔQ||_F~2 ||Δ∑||_F~2≤||ΔA||_F~2 which are optimal for each factor.Some corresponding absolute perturbation bounds are also given.     

A formally fourth-order accurate compact scheme for 3D Poisson equation in cylindrical coordinates

In this paper, we extend our previous work (M.-C. Lai, A simple compact fourth-order Poisson solver on polar geometry, J. Comput. Phys. 182 (2002) 337–345) to 3D cases. More precisely, we present a spectral/finite difference scheme for Poisson equation in cylindrical coordinates. The scheme relies on the truncated Fourier series expansion, where the partial differential equations of Fourier coefficients are solved by a formally fourth-order accurate compact difference discretization. Here the formal fourth-order accuracy means that the scheme is exactly fourth-order accurate while the poles are excluded and is third-order accurate otherwise. Despite the degradation of one order of accuracy due to the presence of poles, the scheme handles the poles naturally; thus, no pole condition is needed. The resulting linear system is then solved by the Bi-CGSTAB method with the preconditioner arising from the second-order discretization which shows the scalability with the problem size.     

Compressible flow equations based on moving coordinates determined by the wave speed

In this paper, we consider moving coordinates with a speed determined by the wave speed. This allows us to move the coordinates at the speed of contact discontinuities, expansion fans, or other wave structures. While recognizing that the Lagrangian coordinates, which can well resolve the sliplines, behave badly for expansion waves, we demonstrate that a coordinate system moving at the characteristic speed of the expansion fan behaves much better for expansion waves. Moreover, the new coordinate system allows one to well capture the shock wave, the slipline and the expansion waves at the same time. Copyright © 2006 John Wiley & Sons, Ltd.     

Multiple ground-based observations at Zhongshan Station during the April/May 1998 solar events

Simultaneous observations at Zhongshan Station, Antarctica, during May 1-7, 1998 are presented to show the responses of the polar ionosphere to the April/May 1998 solar events. One of the main geo-effects of the solar events resulted in the major magnetic storm on May 4. At the storm onset on May 2 the ionosphere F2 layer abruptly increased in altitude, the geomagnetic H-component started negative deviation and the spectral amplitude of the ULF wave intensified. Both large isolated riometer absorption and large negative deviation of the geomagnetic H-component occurred at about 0639UT. There was a time lag of about one hour and ten minutes between the storm onset and the IMF southward turning, as measured by the WIND satellite. The polar ionosphere was highly disturbed, as shown by frequent large deviations of the geomagnetic H-component, large riometer absorption events and strong ULF waves in all the courses of the storm. The absorption increased greatly causing the digisonde to be blackout most of the time. However, the data still showed a substantial decrease in the F2 electron density and oscillation of the F2 layer peak height with an amplitude exceeding 200 km.     

超球坐标下Li原子基态能量的快速收敛计算

将Ｌｉ原子２Ｓ态的波函数分解为两项之积，ψ＝ＸΦ，Ｘ＝ｅｘｐ［－α（ｒ１＋ｒ２＋ｒ３）］，而Φ向以九维超球谐为投影函数而构造的Ｓ３置换群的二维不可约表示的基函数展开。Ｓｃｈｒｏｄｉｎｇｅｒ方程转化为一起球径耦合微分方程。再利用广义Ｌａｇｕｅｒｒｅ函数展开法对其求解，从而得到一广义能量本征方程，给出了库仑势及有关超球角函数在非对称九维超球谐之间矩阵元的简单解析式，可由一起几何函数４Ｆ３：（Ａ，Ｂ，Ｃ，Ｄ；Ｅ，Ｆ，Ｇ；１）的线性组合表示。取超球谐函数的子集（ｌ１＝ｌ２＝ｌ３＝０），基态１２Ｓ的本征能量在超球角方向的收敛速度加快，首次用非变分法得到了优于Ｈａｒｔｒｅｅ－Ｆｏｃｋ精度的基态能量。结果表明本文给出的方法是对直接超球谐广义Ｌａｇｕｅｒｒｅ函数展开法（ＨＨ－ＧＬＦ）的一个非常有效的改进。     

完整系统三阶Lagrange方程的一种推导与讨论

从牛顿运动方程出发，推导了完整系统关于广义加速度的Lagrange方程.讨论了该方程与传统分析力学中的Lagrange方程的相容性问题.结果显示，三阶Lagrange方程可以通过对Lagrange方程求一阶时间导数得到，表明它们是相容的.因此三阶Lagrange方程提供了一种不同于传统Lagrange方程方法的求解物体运动方程的途径. 关键词： Lagrange方程 加速度能量 广义坐标     

Extended affine Weyl groups and Frobenius manifolds

We define certain extensions of affine Weyl groups (distinct from these considered by K. Saito [S1] in the theory of extended affine root systems), prove an analogue of Chevalley Theorem for their invariants, and construct a Frobenius structure on their orbit spaces. This produces solutions F(t₁, ..., t_n) of WDVV equations of associativity polynomial in t₁, ..., t_n-1, exp t_n.