Poisson Kernel for a Degenerate Elliptic Equation and Application to Boundary Value Problems

In this paper we consider the Dirichlet problem to a degenerate elliptic equation in a domain whose interior contains a degenerate surface. By means of the method of expansion of Poisson kernel and applying the properties of special functions, we obtain the twice continuously differentiable solution of the problem on the entire space including infinity.     

Analysis of a Diagonal Form of the Fast Multipole Algorithm for Scattering Theory

Discretisation of the integral equations of acoustic scattering yields a system of linear equations with full coefficient matrices. In recent years a number of fast algorithms for the solution of this system have been proposed. In this paper we present a complete analysis for a fast multipole method for the Helmholtz equation. A one-level diagonal form of the multipole method is applied to a hypersingular integral equation arising from 2d scattering theory. The error of the approximation is analysed and the results used to establish the complexity of the method.     

A Posteriori Error Analysis of a Fully-Mixed Finite Element Method for a Two-Dimensional Fluid-Solid Interaction Problem

In this paper we develop an a posteriori error analysis of a fully-mixed finite element method for a fluid-solid interaction problem in 2D. The media are governed by the elastodynamic and acoustic equations in time-harmonic regime, respectively, the transmission conditions are given by the equilibrium of forces and the equality of the corresponding normal displacements, and the fluid is supposed to occupy an annular region surrounding the solid, so that a Robin boundary condition imitating the behavior of the Sommerfeld condition is imposed on its exterior boundary. Dual-mixed approaches are applied in both domains, and the governing equations are employed to eliminate the displacement u of the solid and the pressure $p$ of the fluid. In addition, since both transmission conditions become essential, they are enforced weakly by means of two suitable Lagrange multipliers. The unknowns of the solid and the fluid are then approximated by a conforming Galerkin scheme defined in terms of PEERS elements in the solid, Raviart-Thomas of lowest order in the fluid, and continuous piecewise linear functions on the boundary. As the main contribution of this work, we derive a reliable and efficient residual-based a posteriori error estimator for the aforedescribed coupled problem. Some numerical results confirming the properties of the estimator are also reported.     

Error Estimates for the Numerical Approximation of Boundary Semilinear Elliptic Control Problems

We study the numerical approximation of boundary optimal control problems governed by semilinear elliptic partial differential equations with pointwise constraints on the control. The analysis of the approximate control problems is carried out. The uniform convergence of discretized controls to optimal controls is proven under natural assumptions by taking piecewise constant controls. Finally, error estimates are established and some numerical experiments, which confirm the theoretical results, are performed.The first two authors were supported by Ministerio de Ciencia y Tecnología (Spain). The second author was also supported by the DFG research center “Mathematics for key technologies” (FZT86) in Berlin.     

一类非线性椭圆型方程奇摄动问题

A class of singularly perturbed problems for the nonlinear elliptic equations is considered. Under suitable conditions, using the theory of differential inequalities the asymptotic behavior of solution for the boundary value problems are studied, which reduced equations possess two intersecting solutions.     

一种基于LVI求解二次规划问题的数值算法

给出并研究了一种数值算法(简称94LVI算法),用于求解带等式和双端约束的二次规划问题. 这类带约束的二次规划问题首先被转换为线性变分不等式问题,该问题等价于分段线性投影等式.接着使用94LVI算法求解上述分段线性投影等式,从而得到QP问题的最优解. 进一步给出了94LVI算法的全局收敛性证明. 94LVI算法与经典有效集算法的对比实验结果证实了给出的94LVI算法在求解二次规划问题上的高效性与优越性.     

Algorithm Research Based on the PDE Sensitivity Filter北大核心CSCD

采用PDE灵敏度滤波器可以消除连续体结构拓扑优化结果存在的棋盘格现象、数值不稳定等问题,且PDE灵敏度滤波器的实质是具有Neumann边界条件的Helmholtz偏微分方程.针对大规模PDE灵敏度滤波器的求解问题,有限元分析得到其代数方程,分别采用共轭梯度算法、多重网格算法和多重网格预处理共轭梯度算法对代数方程进行求解,并且研究精度、过滤半径以及网格数量对拓扑优化效率的影响.结果表明：与共轭梯度算法和多重网格算法相比,多重网格预处理共轭梯度算法迭代次数最少,运行时间最短,极大地提高了拓扑优化效率.     

二维非稳态对流扩散边界控制问题的简化算法

针对二维非稳态对流扩散边界控制问题计算量大的问题,提出了基于降阶模型的最优实时控制方法.利用POD(the Proper Orthogonal Decomposition)和奇异值分解以及Galerkin投影方法得到了具有高精度离散形式的状态空间降阶模型.在所得的降阶状态空间模型中,利用离散时间线性二次调节器方法设计出了最优控制器.对流-扩散过程的控制模拟结果说明了所提方法的有效性和准确性.     

基于凝聚函数的互补问题的光滑化算法

对于不可微的"极大值"形式的函数,可以利用凝聚函数对其进行光滑逼近.借助这个技术,给出了求解线性互补问题的光滑方程组算法.首先是将互补问题转化为等价的非光滑方程组,再利用凝聚函数进行光滑逼近,从而转化为光滑方程组的求解问题.通过一些考题对这个算法进行了数值试验,结果显示了该算法的有效性和稳定性.     

Fast Computational Procedure for Solving Multi-Item Single-Machine Lot Scheduling Optimization Problems

In this paper, we deal with the numerical solution of the optimal scheduling problem in a multi-item single machine. We develop a method of discretization and a computational procedure which allows us to compute the solution in a short time and with a precision of order k, where k is the discretization size. In our method, the nodes of the triangulation mesh are joined by segments of trajectories of the original system. This special feature allows us to obtain precision of order k, which is in general impossible to achieve by usual methods. Also, we develop a highly efficient algorithm which converges in a finite number of steps.     

A Dual Dynamic Programming for Multidimensional Elliptic Optimal Control Problems

This paper deals with the optimal control problems with multiple integrals and an elliptic partial differential equation. The sufficient conditions for optimality in these problems are proved through a dual dynamic programming. The concept of an optimal dual feedback is introduced, and the theorem guaranteeing its existence is established. For the purposes of numerical methods, the ε-version of the verification theorem provided appears to be very useful.     

线性规划的目标函数最速递减算法

在对偶单纯形方法的基础上，提出了线性规划的目标函数最速递减算法。它避开求初始可行基或初始基，以目标函数全局快速递减作为选基准则，将选基过程与换基迭代合二为一，从而大大减少了迭代次数。数值算例显示了该算法的有效性和优越性。     

New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

Boundary-value Problems for Integro-differential Equations of Elliptic Type

In this paper we study the existence of solutions to the Dirichlet problem for a class of integro-differential equations of elliptic type by using the weakly continuous method.     

Three Dimensional Interface Problems for Elliptic Equations

The author studies the structure of solutions to the interface problems for second order linear elliptic partial differential equations in three space dimension.The set of singular points consists of some singular lines and some isolated singular points.It is proved that near a singular line or a singular point,each weak solution can be decomposed into two parts,a singular part and a regular part.The singular parts are some finite sum of particular solutions to some simpler equations,and the regular parts are bounded in some norms,which are slightly weaker than that in the Sobolev space H~2.     

一类半线性椭圆型方程奇摄动广义边值问题

讨论了一类非线性椭圆方程奇摄动广义边值问题。在适当的条件下,研究了边值问题广义解的存在、唯一性及其渐近性态。     

基于Burton-Miller边界积分方程的二维声学波动问题对角形式快速多极子边界元及其应用

论述了二维声学问题的快速多极子边界元（FMBEM）方程及实现步骤．概述了核函数展开理论,并对FMBEM的4个重要组成部分:源点矩计算、源点矩转移、源点矩至本地展开转移、本地展开转移进行了详细的描述．提出了一种有利于四叉树建立的数据结构．推导了一种比直接数值计算更精确、稳定和高效的解析源点矩计算公式．数值算例验证了FMBEM的正确性和高效性．最后,使用FMBEM对轨道二维声学辐射模型进行了模拟计算．     

椭圆外区域上Helmholtz问题的耦合法

以椭圆外区域上Helmholtz方程为例,研究一种带有椭圆人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及误差分析并给出数值例子.理论分析及数值结果表明,用方法求解椭圆外问题是十分有效的.为求解具有长条型内边界外Helmholtz问题提供了一种很好的数值方法.     

求解不等式约束优化问题的一个改进算法

本文改进了一个求解不等式约束优化问题的对偶算法，建立了一个相应的算法，进一步证明了该算法的收敛性．最后，给出数值结果以验证该算法的有效性。