A comparison between a collocation and weak implementation of the rigid‐body motion constraint on a particle surface

In this work, we implemented and compared two different methods to impose the rigid‐body motion constraint on a solid particle moving inside a fluid. We consider a fictitious domain method to easily manage the particle motion. As the solid as well as the fluid inertia are neglected, the particle can be discretized through its boundary only. The rigid‐body motion is imposed via Lagrange multipliers on the boundary. In the first method, such constraints are imposed in discrete points on the boundary (collocation), whereas in the second the constraint is imposed in a weak way on elements dividing the particle surface. Two test problems, that is, a spherical and an ellipsoidal particle in a sheared Newtonian fluid, are chosen to compare the methods. In both cases, the analysis is carried out in 2D as well as in 3D. The results show that for the collocation method an optimal number of collocation points exist leading to the smallest error. However, small variations in the optimal value can generate large deviations. In the weak implementation, the error is only mildly affected by the number of elements used to discretize the particle boundary and by the Lagrange multiplier's interpolation space. A further analysis is carried out to study the effect of an approximated integration of weak constraints. A comparison between the two methods showed that the same accuracy can be achieved by using less constraints if the weak discretization is used. Finally, the rigid‐body motion imposed via weak constraints leads to better conditioned linear systems. Copyright © 2009 John Wiley & Sons, Ltd.     

A parametric view on the Mangasarian–Fromovitz constraint qualification

The growth of the multipliers, when the parameter approaches such a critical parameter, is characterized by a parametric constraint qualification which is introduced here. It is equivalent to a bound on the growth of the multipliers. Received May 8, 1995 / Revised version received February 12, 1998 Published online February 25, 1999     

基于插值的全矢量有限差分法求解阶跃光波导本征模式

给出了一种适合求解阶跃折射率分布光波导本征模式的全矢量有限差分法.该方法充分考虑了电场分量在折射率分布跃变处的不连续性, 在拉格朗日插值的基础上得到了一种精度较高的差分格式.为了验证其有效性,对一种典型结构、折射率阶跃分布的波导作了计算,并与其它方法的结果作了比较.最后对不理想刻蚀情况做了模拟计算.     

大型动态规划的分解算法

An approach about large dynamic programming based on discrete linear system with a quadratic index function is proposed by importing two Lagrange multipliers.     

Criterion for Lagrange stability of a class of pendulum-like systems in Hilbert spaces

In this paper, a class of pendulum-like systems in infinite-dimensional Hilbert spaces is investigated. A frequency domain criterion guaranteeing the systems with the property of Lagrange stability is established.     

Constraint qualifications in multiobjective optimization problems: Differentiable case

In this paper, we are concerned with a multiobjective optimization problem with inequality constraints. We introduce a constraint qualification and derive the Kuhn-Tucker type necessary conditions for efficiency. Moreover, we give conditions which ensure the constraint qualification.This work was done while the author was visiting the University of California, Berkeley, California.     

真空光电探测器件的发展动向及其应用

真空光电探测器件的发展动向及其应用陈举忠（南京电子器件研究所，２１００１６）［摘要］介绍光电倍增管、高灵敏度高分辨率视像管、ＥＢＣＣＤ／像增强器、条纹管类光电探测器的发展现状及其应用。关键词：光电探测器，光电倍增管，视像管，条纹管ＴＨＥＤＥＶＥＬＯＰ...     

Asymptotic analysis of stochastic programs

In this paper we discuss a general approach to studying asymptotic properties of statistical estimators in stochastic programming. The approach is based on an extended delta method and appears to be particularly suitable for deriving asymptotics of the optimal value of stochastic programs. Asymptotic analysis of the optimal value will be presented in detail. Asymptotic properties of the corresponding optimal solutions are briefly discussed.     

Finite element and boundary element contact stress analysis with remeshing technique

The fundamental part of the contact stress problem solution using a finite element method is to locate possible contact areas reliably and efficiently. In this research, a remeshing technique is introduced to determine the contact region in a given accuracy. In the proposed iterative method, the meshes near the contact surface are modified so that the edge of the contact region is also an element’s edge. This approach overcomes the problem of surface representation at the transition point from contact to non-contact region. The remeshing technique is efficiently employed to adapt the mesh for more precise representation of the contact region. The method is applied to both finite element and boundary element methods. Overlapping of the meshes in the contact region is prevented by the inclusion of displacement and force constraints using the Lagrange multipliers technique. Since the method modifies the mesh only on the contacting and neighbouring region, the solution to the matrix system is very close to the previous one in each iteration. Both direct and iterative solver performances on BEM and FEM analyses are also investigated for the proposed incremental technique. The biconjugate gradient method and LU with Cholesky decomposition are used for solving the equation systems. Two numerical examples whose analytical solutions exist are used to illustrate the advantages of the proposed method. They show a significant improvement in accuracy compared to the solutions with fixed meshes.     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.