非协调区域分解的Lagrangian乘子法

§1.引言 近年来随着并行计算机的迅速发展,求解椭圆型方程的区域分解法愈来愈引起人们的兴趣和重视.但是,目前能够见到的有限元区域分解法几乎都要求有限元空间在跨过子区域的边界时是协调的,必然限制有限元区域分解算法的优越性. [3]提出了一种非协凋区域分解法——非协调区域分解的杂交法.采用简化杂交法处理各子区域交界处的非协调性,这种方法在子区域的内部和边界采用两套不同的变量,允许内部变量在跨过各子区域的边界时不连续.但是这种方法有它的局限性,即要求边界变量在各子区域的顶点处必须保持连续性,这对推广到三维空间的情形带来很大的困难.本文提出一种非协调区域分解的Lagrangian乘子法,引进Lagrangian乘子来处理各子区域交界处的非协调性.这种方法也在子区域内部和边界采用两套不同的变量,它不仅允许内部变量在越过各子区域边界时的非协调性,并且还允许边界变量在各子区域的顶点处可以不连续,这就弥补了[3]的不足.同时,这种算法具有[3]的优点,即在不     

基于几何非协调分解的区域分解方法误差分析

对三维可压线弹性问题采用一种基于几何非协调分解的区域分解方法进行求解,证明了数值解具有最优误差估计.     

基于Agent与分解协调的综合生产计划研究

以作业单元为局部决策Agent,车间管理者为全局协调Agent,引入生产节点间的内部结算价格,基于多Agent系统,建立了综合生产计划的分布式决策模型.通过将局部Agent决策目标的总和与全局Agent决策目标进行对比,证明了所引入的内部结算价格就是全局Agent目标函数关于物流平衡约束的Lagrange乘子.基于Lagrange分解协调原理,设计了局部作业单元Agent和全局协调Agent的迭代协调算法.该迭代算法以上次计算的中间结果作为对其它作业单元生产需求的估计,从而能将各个生产单元Agent的决策模型分离,实现了分布建模与求解.在算例研究中使用启发式规则来确定Lagrange乘子迭代的步长系数,保证了较好的收敛性,证明模型和算法是有效的.     

Lagrangian乘子区域分解法的一类预条件子

1．引言非重叠区域分解的Lagrangian乘子法已被许多作者讨论［1今它是一类非协调区域分解法（与通常的非协调元区域分解不同），特别适合于非匹配网格的情形（即相邻子域在公共边或公共面上的结点不重合，参见14］［6］）．这种方法的一个最大优点是不要求界面变量在内交点（或内交边）上的连续性，从而界面方程易于建立，程序易于实现，而又正因为这个特点，使得界面矩阵的预条件子不能按通常的方法构造，故目前还未见到理想的预条件子（或者条件数差，或者应用上不方便）．本文在很大程度上解决了这一问题．1）工作单位：湘潭大学数学系…     

一种新的局部时间积分的区域分解波形松弛算法

提出一种新的区域分解波形松弛算法, 使得可以在不同的子域采用不同的时间步长来并行求解线性抛物方程的初边值问题. 与传统的区域分解波形松弛算法相比, 该算法可以通过预条件子来加快收敛速度, 并且对内存的需求大大降低. 给出了局部时间步长一种具体的实现方法, 证明了离散解的存在唯一性, 并在时间连续水平分析了预条件系统. 数值实验显示了新算法的有效性.     

非协调区域分解Lagrange乘子法的超收敛

本文讨论分析非协调区域分解Ｌａｇｒａｎｇｅ乘子法对二阶椭圆型方程Ｄｉｒｉｃｈｌｅｔ问题的有限元超收敛现象。文中通过利用积分恒等式，适宜地引进Ｌ２投影过渡以及高次插值后处理等技巧，经过一系列误差分析及估计，得到了高出半阶的超收敛结果，实现了非协调区域分解法与高精度算法的结合。     

带有不等式约束的非线性规划问题的一个精确增广Lagrange函数

对求解带有不等式约束的非线性非凸规划问题的一个精确增广Lagrange函数进行了研究．在适当的假设下,给出了原约束问题的局部极小点与增广Lagrange函数,在原问题变量空间上的无约束局部极小点之间的对应关系．进一步地,在对全局解的一定假设下,还提供了原约束问题的全局最优解与增广Lagrange函数,在原问题变量空间的一个紧子集上的全局最优解之间的一些对应关系．因此,从理论上讲,采用该文给出的增广Lagrange函数作为辅助函数的乘子法,可以求得不等式约束非线性规划问题的最优解和对应的Lagrange乘子．     

带非局部边界条件的非线性抛物方程的隐式差分解法

在准静态弹性力学中常遇到求解带有非局部边界条件的抛物方程初边值问题.本文构造了一个数值求解带有非局部边界条件的非线性抛物方程的隐式差分格式,利用离散泛函分析的知识和不动点定理证明了差分解是存在的,且在离散最大模意义下关于时间步长一阶收敛,关于空间步长二阶收敛,并给出了数值算例.     

求解界面方程的Chebyshev加速算法

非重叠区域分解算法在于建立和求解相关的界面方程．建立界面方程在理论上虽。然容易推导，例如某些问题可用Gauss块消去法，但在实际计算时并不可行，所以界面方程在一些算法中是陷式的．而求解界面方程一般要进行预处理，本提出一种区域分解算法，可得出界面方程的显式表达．算法是完全并行的，所得出的界面方程的系数矩阵的条件数已与网参数无关，事实上就是(Sh^（1）)^-1Sh，进而可直接用收敛速度较快的Chebyshev加速算法求解该界面方程，在充分应用并行计算方法的条件下，本算法与[4]中的算法相比计算效率提高．     

另一类非完整力学系统的Lagrange方程

用文[1]的方法,导出另一类一阶非完整力学系统不带乘子的Lagrange方程.这种形式的方程也是新的.     

A mortar element method for hyperbolic problems

A non-conforming finite element method based on non-overlapping domain decomposition is extended to linear hyperbolic problems. The method is based on streamline-diffusion/discontinuous Galerkin methods and the mortar element method. A weak flux continuity condition at the inflow interface is enforced by means of Lagrange multipliers. This weak flux continuity condition replaces the usual mortar condition for elliptic problems, and allows non-matching grids at the subdomain interfaces. To cite this article: Y. Bourgault, A. El Boukili, C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

A regularized domain decomposition method with Lagrange multiplier

In this paper, we are concerned with the nonoverlapping domain decomposition method with Lagrange multiplier for three-dimensional second-order elliptic problems with no zeroth-order term. It is known that the methods result in a singular subproblem on each internal (floating) subdomain. To handle the singularity, we propose a regularization technique which transforms the corresponding singular problems into approximate positive definite problems. For the regularized method, one can build the interface equation of the multiplier directly. We first derive an optimal error estimate of the regularized approximation, and then develop a cheap preconditioned iterative method for solving the interface equation. For the new method, the cost of computation will not be increased comparing the case without any floating subdomain. The effectiveness of the new method will be confirmed by both theoretical analyzes and numerical experiments. The work is supported by Natural Science Foundation of China G10371129.     

从几何角度给予拉格朗日乘数法新的推导思路

拉格朗日乘数法是求条件极值的重要方法,该文通过数形结合给出定理推导的新路径,相比教材上纯代数推导更直观,体现了"几何意义"的重要性.     

高阶拉氏乘子法和弹性理论中更一般的广义变分原理

作者曾指出[1],弹性理论的最小位能原理和最小余能原理都是有约束条件限制下的变分原理采用拉格朗日乘子法,我们可以把这些约束条件乘上待定的拉氏乘子,计入有关变分原理的泛函内,从而将这些有约束条件的极值变分原理,化为无条件的驻值变分原理.如果把这些待定拉氏乘子和原来的变量都看作是独立变量而进行变分,则从有关泛函的驻值条件就可以求得这些拉氏乘子用原有物理变量表示的表达式.把这些表达式代入待定的拉氏乘子中,即可求所谓广义变分原理的驻值变分泛函.但是某些情况下,待定的拉氏乘子在变分中证明恒等于零.这是一种临界的变分状态.在这种临界状态中,我们无法用待定拉氏乘子法把变分约束条件吸收入泛函,从而解除这个约束条件.从最小余能原理出发,利用待定拉氏乘子法,企图把应力应变关系这个约束条件吸收入有关泛函时,就发生这种临界状态,用拉氏乘子法,从余能原理只能导出Hellinger-Reissner变分原理[2],[3],这个原理中只有应力和位移两类独立变量,而应力应变关系则仍是变分约束条件,人们利用这个条件,从变分求得的应力中求应变.所以Hellinger-Reissner变分原理仍是一种有条件的变分原理.     

Analysis and computation of a least-squares method for consistent mesh tying

In the finite element method, a standard approach to mesh tying is to apply Lagrange multipliers. If the interface is curved, however, discretization generally leads to adjoining surfaces that do not coincide spatially. Straightforward Lagrange multiplier methods lead to discrete formulations failing a first-order patch test [T.A. Laursen, M.W. Heinstein, Consistent mesh-tying methods for topologically distinct discretized surfaces in non-linear solid mechanics, Internat. J. Numer. Methods Eng. 57 (2003) 1197–1242].     

一个求解条件极值问题的极值点的新方法

利用齐次线性方程组理论,建立了一个求解条件极值问题的极值点的新方法.该方法的优点是:能有效地避免在运用Lagrange乘数法求解条件极值时,因引进了参数而给解方程组带来的困扰.也可以说,对于有些问题我们仅从已知条件入手,不必引进参数就可以直接求得极值点.     

A Robin-type non-overlapping domain decomposition procedure for second order elliptic problems

This article deals with the analysis of an iterative non-overlapping domain decomposition (DD) method for elliptic problems, using Robin-type boundary condition on the inter-subdomain boundaries, which can be solved in parallel with local communications. The proposed iterative method allows us to relax the continuity condition for Lagrange multipliers on the inter-subdomain boundaries. In order to derive the corresponding discrete problem, we apply a non-conforming Galerkin method using lowest order Crouzeix–Raviart elements. The convergence of the iterative scheme is obtained by proving that the spectral radius of the matrix associated with the fixed point iterations is less than 1. Parallel computations have been carried out and the numerical experiments confirm the theoretical results established in this paper.     

Applications of the method of partial inverses to convex programming: Decomposition

A primal–dual decomposition method is presented to solve the separable convex programming problem. Convergence to a solution and Lagrange multiplier vector occurs from an arbitrary starting point. The method is equivalent to the proximal point algorithm applied to a certain maximal monotone multifunction. In the nonseparable case, it specializes to a known method, the proximal method of multipliers. Conditions are provided which guarantee linear convergence.This research was sponsored, in part, by the Air Force Office of Scientific Research under grant 80-0195.     

A Dual-Primal FETI method for incompressible Stokes equations

In this paper, a dual-primal FETI method is developed for incompressible Stokes equations approximated by mixed finite elements with discontinuous pressures. The domain of the problem is decomposed into nonoverlapping subdomains, and the continuity of the velocity across the subdomain interface is enforced by introducing Lagrange multipliers. By a Schur complement procedure, the solution of an indefinite Stokes problem is reduced to solving a symmetric positive definite problem for the dual variables, i.e., the Lagrange multipliers. This dual problem is solved by the conjugate gradient method with a Dirichlet preconditioner. In each iteration step, both subdomain problems and a coarse level problem are solved by a direct method. It is proved that the condition number of this preconditioned dual problem is independent of the number of subdomains and bounded from above by the square of the product of the inverse of the inf-sup constant of the discrete problem and the logarithm of the number of unknowns in the individual subdomains. Numerical experiments demonstrate the scalability of this new method. This work is based on a doctoral dissertation completed at Courant Institute of Mathematical Sciences, New York University. This work was supported in part by the National Science Foundation under Grants NSF-CCR-9732208, and in part by the U.S. Department of Energy under contract DE-FG02-92ER25127.     

A mixed hybrid formulation based on oscillated finite element polynomials for solving Helmholtz problems

A mixed-hybrid-type formulation is proposed for solving Helmholtz problems. This method is based on (a) a local approximation of the solution by oscillated finite element polynomials and (b) the use of Lagrange multipliers to “weakly” enforce the continuity across element boundaries. The computational complexity of the proposed discretization method is determined mainly by the total number of Lagrange multiplier degrees of freedom introduced at the interior edges of the finite element mesh, and the sparsity pattern of the corresponding system matrix. Preliminary numerical results are reported to illustrate the potential of the proposed solution methodology for solving efficiently Helmholtz problems in the mid- and high-frequency regimes.