Sub(L)中的五边形格特征及其数量不等式

五边形结构在刻画格的特征方面具有十分重要作用，应用格论及组合数学的方法讨论了格L与其子格格Sub(L)中所含五边形格之间的数量关系，给出了有限格L的子格格中三个元生成五边形格的充要条件，同时给出了Sub(L)所含不同五边形格数量的一个下界。     

Localization effects and measure source terms in numerical schemes for balance laws

This paper investigates the behavior of numerical schemes for nonlinear conservation laws with source terms. We concentrate on two significant examples: relaxation approximations and genuinely nonhomogeneous scalar laws. The main tool in our analysis is the extensive use of weak limits and nonconservative products which allow us to describe accurately the operations achieved in practice when using Riemann-based numerical schemes. Some illustrative and relevant computational results are provided.

     

一维光学格子孤子的传输特性及控制研究

利用解析和数值方法研究了在具有横向折射率周期性调制的克尔型非线性介质中光学格子孤子的传输，得到了孤子参数的演化方程以及格子孤子的形成和稳定传输的条件.结果表明：当光束的入射角小于某临界角度时，光束可被类似波导形式的路径俘获而稳定传输，该临界角随折射率调制周期、调制深度的增加而增大，且光束越窄临界值越大.此外，线性空间啁啾虽然对光束传输的中心位置没有任何影响，但会导致光束发散从而破坏格子孤子的形成和稳定传输，对此提出了采用特定功率取值来补偿啁啾作用从而形成格子孤子的方案. 关键词： 光孤子 光学格子 光传输 矩方法     

宽带膜厚实时监控过程中膜层折射率的确定方法

为了准确计算出镀膜过程中每层膜的折射率，介绍了实时监控过程中确定膜层折射率的2种方法：一种是由实测的透射比光谱直接反算出膜层的折射率；另一种是用最小二乘法的优化算法实时拟合折射率。试验结果表明：在线反算适合单点监控，所得折射率误差小于2%。然而在实际镀膜过程中,由于宽带内膜层参数误差较大,一般大于25%。为此，采用最小二乘法拟合，即在整个宽光谱范围内采集每个波长点的信息，所得结果误差很小，一般都在2％～5%之间,有时可达到10%,在很大程度上提高了实际镀膜时膜厚监控的精度。     

A Modified Quasi-Newton Method for Structured Optimization with Partial Information on the Hessian

This paper develops a modified quasi-Newton method for structured unconstrained optimization with partial information on the Hessian, based on a better approximation to the Hessian in current search direction. The new approximation is decided by both function values and gradients at the last two iterations unlike the original one which only uses the gradients at the last two iterations. The modified method owns local and superlinear convergence. Numerical experiments show that the proposed method is encouraging comparing with the methods proposed in [4] for structured unconstrained optimization Presented at the 6th International Conference on Optimization: Techniques and Applications, Ballarat, Australia, December 9–11, 2004     

Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Acoustic reflection from a baffled elastic/piezoelectric plate

A finite difference/boundary integral procedure to determine the acoustic reflected pressure from a fluid-loaded bi-laminated plate is described. The bi-laminate is composed of a piezoelectric layer and an elastic layer in contact with the fluid, and is held by an acoustically hard baffle. In the numerical model, the fluid pressure at fluid/solid interface is replaced by a continuum of point sources weighted by the normal acceleration of the elastic plate, and the governing equation system is solved in the solid domain. With the normal acceleration found, the reflected pressure in the fluid is determined by an integral expression involving the Green's function. It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling either the reflected or scattered pressure of the plate at any chosen field points in the fluid. Project supported by the National Natural Science Foundation of China (No. 10172039).     

烷基硫酸盐表面张力的量子化学研究

用AMl量子化学计算法优化了十三烷基硫酸根阴离子和——CH3在不同取代位的十二烷基硫酸根阴离子的几何构型，得到最优构型时最高占据分子轨道能级EHOMO、最低空轨道能级ELUMO、电子能量Eele和偶极矩μ等数据．将这些电子结构数据分别与表面张力相拟合，得到很好的相关性。文中讨论了——CH3在不同位置取代对表面张力的影响。     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

解奇异非光滑方程组的牛顿法和不精确牛顿法

本文主要解决奇异非光滑方程组的解法。应用一种新的次微分的外逆，我们提出了牛顿法和不精确牛顿法，它们的收敛性同时也得到了证明。这种方法能更容易在一引起实际应用中实现。这种方法可以看作是已存在的解非光滑方程组的方法的延伸。