A method for calculating invariant subspaces of symmetric hyperbolic equations

An algorithm is constructed for calculating invariant subspaces of symmetric hyperbolic systems arising in electromagnetic, acoustic, and elasticity problems. Discrete approximations are calculated for subspaces that correspond to minimal eigenvalues and smooth eigenfunctions. Difficulties related to the presence of an infinite-dimensional kernel in the differential operator are successfully handled. The efficiency of the algorithm is demonstrated using acoustics equations.     

Comparison of growth stress measurements with modelling in thin iron oxide films

High temperature oxidation of metals leads to residual stresses both in the metal and in the growing oxide. In this work, the evolution of this residual stresses is theoretically predicted in the growing oxide layers. The origin of these stresses is based on a microstructural model. Using experimental results providing from the oxidation kinetics, and an analysis proposed to describe the growth strain occurring in the thin layers, a set of equations is established allowing determining the stresses evolution with oxidation time. Then, the model is compared with experimental results obtained on both α-Fe and phosphated α-Fe, oxidised at different temperatures. Numerical data are extracted from experiments either with an asymptotic formulation or with an inverse method. These two methods give good agreement with experiments and allow extracting the model parameters.     

有限时滞微分系统依照双测度的实际稳定性

本文对具有有限时滞的泛函微分方程建立了关于依照两种测度的实际稳定性的Razumikhin型判定定理，其中未采用通常的辅助函数，且可运用多个含有状态变量x的部分变元的Lyapunov函数，得出部分变远实际稳定性的判定定理，从而改进了已有的结果。     

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).