一类各向异性外问题的重叠型区域分解算法

本文以椭圆外调和问题的自然边界归化为基础，提出了求解各向异性常系数椭圆方程的一种重叠型区域分解算法，并分析了算法的收敛性及收敛速度．理论分析及数值实验表明，该方法对于求解各向异性外问题非常有效．     

Signorini传输问题的有限元-边界元耦合

对一类带Signorini接触条件的非线性传输问题, 提出了新的有限元-边界元耦合框架.为求解所得到的耦合变分不等式, 设计了一种区域分解型迭代方法, 并对其做了完全的收敛性分析.     

无穷凹角区域各向异性问题的重叠型区域分解算法

本文以凹角椭圆外区域上调和问题的自然边界归化为基础,提出了求解无穷凹角区域各向异性问题的重叠型区域分解算法,并分析了算法的收敛性及收敛速度.最后给出了数值例子,以示方法的可行性和有效性.     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

一类耦合的有限元-边界元变分不等式的预条件子

本文主要研究一类Signorini 接触条件的非线性传输问题. 这类问题可以用耦合的有限元- 边界元变分不等式来描述. 我们首先提出一种求解变分不等式的预处理梯度投影法. 然后对离散系统构造了有效的区域分解预条件子. 该预条件子能够使耦合的不等式问题分解成等式问题和小规模的不等式问题, 并且这些问题可以并行求解. 最后我们详细研究了该迭代方法的收敛性.     

箱约束变分不等式的一个简单光滑价值函数和阻尼牛顿法

通过引入中间值函数的一类光滑价值函数,构造了箱约束变分不等式的一种新的光滑价值函数,该函数形式简单且具有良好的微分性质．基于此给出了求解箱约束变分不等式的一种阻尼牛顿算法,在较弱的条件下,证明了算法的全局收敛性和局部超线性收敛率,以及对线性箱约束变分不等式的有限步收敛性．数值实验结果表明了算法可靠有效的实用性能．     

圆锥规划问题的光滑牛顿方法

给出求解圆锥规划问题的一种新光滑牛顿方法.基于圆锥互补函数的一个新光滑函数,将圆锥规划问题转化成一个非线性方程组,然后用光滑牛顿方法求解该方程组.该算法可从任意初始点开始,且不要求中间迭代点是内点.运用欧几里得代数理论,证明算法具有全局收敛性和局部超线性收敛速度.数值算例表明算法的有效性.     

一类耦合的有限元-边界元变分不等式的预条件子 谨以此文致《中国科学》创刊六十周年

本文主要研究一类Signorini接触条件的非线性传输问题.这类问题可以用耦合的有限元-边界元变分不等式来描述.我们首先提出一种求解变分不等式的预处理梯度投影法.然后对离散系统构造了有效的区域分解预条件子.该预条件子能够使耦合的不等式问题分解成等式问题和小规模的不等式问题,并且这些问题可以并行求解.最后我们详细研究了该迭代方法的收敛性.     

求解非光滑凸规划的一种混合束方法

提出一种求解非光滑凸规划问题的混合束方法. 该方法通过对目标函数增加迫近项, 且对可行域增加信赖域约束进行迭代, 做为迫近束方法与信赖域束方法的有机结合, 混合束方法自动在二者之间切换, 收敛性分析表明该方法具有全局收敛性. 最后的数值算例验证了算法的有效性．     

一类非线性传输问题的区域分解方法（英文）

针对一类非线性传输问题提出了有限元与边界元的耦合方法并设计了基于耦合法的区域分解算法.该算法避免了求解边界积分方程,从而计算量大大减少.算法的收敛性分析和数值算例验证了该算法的合理和有效性.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

二阶三点边值问题解的存在性与唯一性

利用单调迭代方法,研究了反序上、下解条件下的二阶微分方程三点边值问题解的存在性与唯一性,分别得到解存在与唯一的充分条件,在满足解的唯一性的条件下,给出了求解的迭代序列及误差估计式.     

数学物理中的总体和局部场GL方法

为了求解物理化学生物材料和金融中的微分方程,提出了一种总体(Global)和局部(Local)场方法.微分方程的求解区域可以是有限域,无限域,或具曲面边界的部分无限域.其无限域包括有限有界不均匀介质区域.其不均匀介质区域被分划为若干子区域之和.在这含非均匀介质的无限区域,将微分方程的解显式地表示为在若干非均匀介质子区域上和局部子曲面的积分的递归和.把正反算的非线性关系递归地显式化.在无限均匀区域,微分方程的解析解被称为初始总体场.微分方程解的总体场相继地被各个非均匀介质子区域的局部散射场所修正.这种修正过程是一个子域接着另个子域逐步相继地进行的.一旦所有非均匀介质子区域被散射扫描和有限步更新过程全部完成后,微分方程的解就获得了.称其为总体和局部场的方法,简称为GL方法.GL方法完全地不同于有限元及有限差方法,GL方法直接地逐子域地组装逆矩阵而获得解.GL方法无需求解大型矩阵方程,它克服了有限元大型矩阵解的困难.用有限元及有限差方法求解无限域上的微分方程时,人为边界及其上的吸收边界条件是必需的和困难的,人为边界上的吸收边界条件的不精确的反射会降低解的精确度和毁坏反算过程.GL方法又克服了有限元和有限差方法的人为边界的困难.GL方法既不需要任何人为边界又不需要任何吸收边界条件就可以子域接子域逐步精确地求解无限域上的微分方程.有限元和有限差方法都仅仅是数值的方法,GL方法将解析解和数值方法相容地结合起来.提出和证明了三角的格林函数积分方程公式.证明了当子域的直经趋于零时,波动方程的GL方法的数值解收敛于精确解.GL方法解波动方程的误差估计也获得了.求解椭圆型,抛物线型,双曲线型方程的GL模拟计算结果显示出我们的GL方法具有准确,快速,稳定的许多优点.GL方法可以是有网,无网和半网算法.GL方法可广泛应用在三维电磁场,三维弹塑性力学场,地震波场,声波场,流场,量子场等方面.上述三维电磁场等应用领域的GL方法的软件已经由作者研制和发展了。     

Uniform regularity and vanishing viscosity limit for the chemotaxis‐Navier‐Stokes system in a 3D bounded domain

We investigate the uniform regularity and vanishing viscosity limit for the incompressible chemotaxis‐Navier‐Stokes system with Navier boundary condition for velocity field and Neumann boundary condition for cell density and chemical concentration in a 3D bounded domain. It is shown that there exists a unique strong solution of the incompressible chemotaxis‐Navier‐Stokes system in a finite time interval, which is independent of the viscosity coefficient. Moreover, this solution is uniformly bounded in a conormal Sobolev space, which allows us to take the vanishing viscosity limit to obtain the incompressible chemotaxis‐Euler system.     

上下解方法在四阶非局部边值问题中的应用

考虑具有p-Laplacian算子的四阶微分方程非局部边值问题.通过构造Green函数,利用上下解方法,结合单调迭代技巧得到了四阶非局部边值问题迭代解的存在性.最后给出实例验证主要结论.     

GL method for solving equations in math-physics and engineering

In this paper, we propose a GL method for solving the ordinary and the partial differential equation in mathematical physics and chemics and engineering. These equations govern the acustic, heat, electromagnetic, elastic, plastic, flow, and quantum etc. macro and micro wave field in time domain and frequency domain. The space domain of the differential equation is infinite domain which includes a finite inhomogeneous domain. The inhomogeneous domain is divided into finite sub domains. We present the solution of the differential equation as an explicit recursive sum of the integrals in the inhomogeneous sub domains. Actualy, we propose an explicit representation of the inhomogeneous parameter nonlinear inversion. The analytical solution of the equation in the infinite homogeneous domain is called as an initial global field. The global field is updated by local scattering field successively subdomaln by subdomain. Once all subdomains are scattered and the updating process is finished in all the sub domains, the solution of the equation is obtained. We call our method as Global and Local field method, in short , GL method. It is different from FEM method, the GL method directly assemble inverse matrix and gets solution. There is no big matrix equation needs to solve in the GL method. There is no needed artificial boundary and no absorption boundary condition for infinite domain in the GL method. We proved several theorems on relationships between the field solution and Green＇s function that is the theoretical base of our GL method. The numerical discretization of the GL method is presented. We proved that the numerical solution of the GL method convergence to the exact solution when the size of the sub domain is going to zero. The error estimation of the GL method for solving wave equation is presented. The simulations show that the GL method is accurate, fast, and stable for solving elliptic, parabolic, and hyperbolic equations. The GL method has advantages and wide applications in the 3D electromagnetic （EM）     

THE ARTIFICIAL BOUNDARY CONDITION FOR EXTERIOR OSEEN EQUATION IN 2-D SPACE

A finite element method for the solution of Oseen equation in exterior domain is proposed. In this method, a circular artificial boundary is introduced to make the computational domain finite. Then, the exact relation between the normal stress and the prescribed velocity field on the artificial boundary can be obtained analytically. This relation can serve as an boundary condition for the boundary value problem defined on the finite domain bounded by the artificial boundary. Numerical experiment is presented to demonstrate the performance of the method.     

Steady solutions of the Navier–Stokes equations with threshold slip boundary conditions

We establish the wellposedness of the time‐independent Navier–Stokes equations with threshold slip boundary conditions in bounded domains. The boundary condition is a generalization of Navier's slip condition and a restricted Coulomb‐type friction condition: for wall slip to occur the magnitude of the tangential traction must exceed a prescribed threshold, independent of the normal stress, and where slip occurs the tangential traction is equal to a prescribed, possibly nonlinear, function of the slip velocity. In addition, a Dirichlet condition is imposed on a component of the boundary if the domain is rotationally symmetric. We formulate the boundary‐value problem as a variational inequality and then use the Galerkin method and fixed point arguments to prove the existence of a weak solution under suitable regularity assumptions and restrictions on the size of the data. We also prove the uniqueness of the solution and its continuous dependence on the data. Copyright © 2006 John Wiley & Sons, Ltd.     

Monotone iterative method for the second-order three-point boundary value problem with upper and lower solutions in the reversed order

This article introduces a new concept of upper and lower solutions and studies the existence and uniqueness of solutions of second-order three-point boundary value problems with upper and lower solutions in the reversed order. The sufficient conditions for the existence and uniqueness of solutions are obtained by using the monotone iterative method, meantime, the iterative sequence for solving a solution and its error estimate formula under the condition of unique solution are given. Some results of previous literature are extended and improved. A numerical example is also included to illustrate the effectiveness of the proposed results.     

Existence and uniqueness of solutions of second-order three-point boundary value problems with upper and lower solutions in the reversed order

This paper studies the existence and uniqueness of solutions of second-order three-point boundary value problems with lower and upper solutions in the reversed order, obtains the sufficient conditions for the existence and uniqueness of solutions by use of the monotone iterative method, and gives the iterative sequence for solving a solution and its error estimate formula under the condition of unique solution.