Extended convergence results for the method of multipliers for nonstrictly binding inequality constraints

In this paper, we extend the classical convergence and rate of convergence results for the method of multipliers for equality constrained problems to general inequality constrained problems, without assuming the strict complementarity hypothesis at the local optimal solution. Instead, we consider an alternative second-order sufficient condition for a strict local minimum, which coincides with the standard one in the case of strict complementary slackness. As a consequence, new stopping rules are derived in order to guarantee a local linear rate of convergence for the method, even if the current Lagrangian is only asymptotically minimized in this more general setting. These extended results allow us to broaden the scope of applicability of the method of multipliers, in order to cover all those problems admitting loosely binding constraints at some optimal solution. This fact is not meaningless, since in practice this kind of problem seems to be more the rule rather than the exception.In proving the different results, we follow the classical primaldual approach to the method of multipliers, considering the approximate minimizers for the original augmented Lagrangian as the exact solutions for some adequate approximate augmented Lagrangian. In particular, we prove a general uniform continuity property concerning both their primal and their dual optimal solution set maps, a property that could be useful beyond the scope of this paper. This approach leads to very simple proofs of the preliminary results and to a straight-forward proof of the main results.The author gratefully acknowledges the referees for their helpful comments and remarks. This research was supported by FONDECYT (Fondo Nacional de Desarrollo Científico y Technológico de Chile).     

Reflection of Surface Plasmon Polaritons by a Corner Reflector of Negative Index Materials

A novel type of corner reflector for the reflection of surface plasmon polariton (SPP) is proposed. The reflector consists of two layers of compensatory media whose permittivity and permeability take opposite signs to that of the corresponding layer of the SP waveguide. By rigorously solving Maxwell's equations, the reflected SPP is proved to exist. The transverse wave vector of the reflected SPP is always antiparallel to that of the incident SPP and no phase retardation is introduced.     

高阶半线性抛物型方程组的生命跨度

其中m,P,q>1．利用试验函数方法,首先推导一些积分不等式,然后对方程组爆破解的生命跨度 [0,T)给出估计．     

On parabolic equations with a small parameter

A singularly perturbed parabolic equation with a nonlinear right-hand side of a special form is examined. A numerical analytical study of such equations is performed.     

Mathematical modelling and analytical solution for workpiece temperature in grinding

This paper deals with modelling the workpiece temperature field produced during the grinding process. The proposed model is given in terms of a two-dimensional boundary-value problem where the interdependence among the grinding wheel, the workpiece and the coolant is described by two variable functions in the boundary condition. An explicit integral form solution is constructed using the Laplace and Fourier transforms and the Green’s function method.     

浅水中污染物扩散的自适应有限元分析法

介绍浅水中污染物扩散分析中的有限元法.分析包括两个部分:1)流场速度、水面高度的计算;2)根据扩散模型计算污染物浓度场.联合使用了自适应网格技术以期提高解的精度,同时减少计算时间和计算机内存的消耗.通过几个有已知解的实例验证了有限元公式和计算机程序.最后,使用这种联合方法分析泰国Chao Phraya河附近海湾中的污染物扩散.     

Effect of core-shell micelle formation on the redox properties of phenothiazine-labeled poly(ethyl glycidy ether)-block-poly(ethylene oxide)

Redox properties of phenothiazine-labeled poly(ethyl glycidy ether)-block-poly(ethylene oxide) (PT-EGE_n-b-EO_m) are reversibly changed by core-shell micelle formation. In the temperature range higher than the critical micellization temperature (cmt), the anodic potential of PT group positively shifts and concomitantly its anodic current decrease, or levels off compared to those of the reference polymer PT-EO_m without the thermo-responsive EGE_n segment. The former alteration is caused by incorporation of hydrophobic PT groups into a core of the micelle and the latter by the decrease in the diffusion coefficient of PT groups due to formation of the core-shell micelles. The cmt value and the temperature-dependent alteration in the redox properties strongly depend on the polymer structure, especially the length of thermo-responsive EGE_n segment. The electrochemically determined hydrodynamic radii of the polymer aggregates seem to be overestimated, compared to the values reported for the aggregates of other thermo-responsive polymers with similar molecular weights, implying the presence of electrochemically inactive PT groups in the copolymers having longer thermo-responsive segments.     

A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

等腰三角形Mindlin板的自由振动分析

提出了一种新方法来对基于 Mindlin剪切变形理论的等腰三角形板进行自由振动分析 .此方法采用了一种新的基函数并利用 pb-2 Rayleigh-Ritz边界函数得到了一种新型的 Ritz方法 .这种方法的有效性通过收敛性和对比性分析得到了证实 .数值结果表明此方法相当精确有效 .     

Finite element modal solutions of planar photonic crystal fibers with rectangular air-holes

Modal solutions of planar photonic crystal waveguides with rectangular air-holes are presented by using a rigorous full-vectorial finite element-based approach. The effective indices, mode field profiles, spot-sizes, power confinements, modal hybridness, beat lengths and group velocity dispersions are shown for the fundamental and higher order modes of the quasi-TE and TM polarizations.