Boundary Integral Equation Solution of Non-linear Plane Potential Problems

This paper presents the boundary integral equation (BIE) formulation,and numerical solution procedure for two-dimensional problemsgoverned by Laplac's equation and subject to non-linear boundaryconditions. The introduction of non-linear terms constitutesa fundamental extension of the BIE method, as previous applicationshave been restricted entirely to linear problems. Furthermore,non-linearities necessitate the use of iterative solution techniqueswhich present the conceptual disadvantage that a solution isnot guaranteed. However, such difficulties were not encounteredwith the Newton—Raphson method employed in this study.The various features of the BIE technique are illustrated bythe application to a physical problem which is of significancein heat exchanger design.     

Two dimensional solitary waves in shear flows

In this paper we study existence and asymptotic behavior of solitary-wave solutions for the generalized Shrira equation, a two-dimensional model appearing in shear flows. The method used to show the existence of such special solutions is based on the mountain pass theorem. One of the main difficulties consists in showing the compact embedding of the energy space in the Lebesgue spaces; this is dealt with interpolation theory. Regularity and decay properties of the solitary waves are also established.     

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

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

Solution of Burgers' equation for large Reynolds number using finite elements with moving nodes

Burgers' equation often arises in the mathematical modelling used to solve problems in fluid dynamics involving turbulence. Numerical difficulties arise in the solution for the case of large Reynolds number. To obtain high accuracy, finite element methods are important. The aim of this paper is to summarize relevant past work and to use a moving node finite element method to obtain a solution of Burgers' equation under certain prescribed conditions. The results for high Reynolds number are compared with accurate results obtained by other authors.     

SOLUTION OF SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD

The aim of this paper is to apply the relatively new Adomian decomposition method to solving the system of linear fractional, in the sense of Riemann-Liouville and Caputo respectively, differential equations. The solutions are expressed in terms of Mittag-Leffler functions of matric argument. The Adomian decomposition method is straightforward, applicable for broader problems and avoids the difficulties in applying integral transforms. As the order is 1, the result here is simplified to that of first order differential equation.     

On an algorithm for finding the electric potential distribution in the DG-MOSFET transistor

An efficient numerical algorithm for finding the electric potential distribution in the DG-MOSFET transistor is proposed and discussed in detail. The class of hydrodynamic models describing the charge transport in semiconductors includes the Poisson equation for the electric potential. Since the equations of hydrodynamic models are nonlinear and involve small parameters and specific conditions on the boundary of the DG-MOSFET transistor domain, the numerical solution of the Poisson equation meets significant difficulties. An original algorithm is proposed that is based on the stabilization method and the idea of schemes without saturation and helps to cope with these difficulties.     

M.E.Shvez迭代解法及其在力学中的应用

当微分方程中含有微量项时,可用M.E.Shvez迭代法求解。但当微量项出现奇性,或在某一区间内微量项并非微量时。用此法求解将遇到困难。本文针对这类问题,把原M.E.Shvez迭代解法稍加改变。算例表明,用改进了的M.E.Shvez法求解上述问题。其精度比原M.E.Shvez法的有所提高。     

An improved boundary element formulation for the motion of spherical bubbles

We suggest a modified boundary element method for modeling the potential flow caused by the motion of many spherical bubbles. The problem for the fluid velocity potential is reduced to a Fredholm integral equation of the second kind. The kernel of that integral equation has no singularity; as a result, its numerical solution does not encounter any difficulties. The matrix representation of the integral equation is diagonally dominant and is well suited to handle multiple bubble system.     

基于不动点方法求解非线性Falkner-Skan流动方程

Falkner-Skan流动方程描述绕楔面的流动,该方程具有很强的非线性.首先通过引入变换式,将原半无限大区域上的流动问题转化为有限区间上的两点边值问题.接着基于泛函分析中的不动点理论,采用不动点方法求解两点边值问题从而得到Falkner Skan流动方程的解.最后将不动点方法给出的结果和文献中的数值结果相比较,发现不动点方法得到的结果具有很高的精度,并且解的精度很容易通过迭代而不断得到提高.表明不动点方法是一种求解非线性微分方程行之有效的方法.     

边界元的分区解法

边界元法是一种比‘区域’型解法(如有限单元法,有限差分法)有着更多优点的方法。由于它能使所考虑的问题在维数上降阶,只需要对结构的边界进行离散,离散误差仅来自边界[1],因而具有输入量小,精度高等优点,在工程界受到广泛的重视.随着边界元网格的增加,其非对称满值的系数矩阵呈自由度平方级的增大,系数矩阵的总量可高达几千上万个。如此巨大的数组空间,一次送入计算机内存进行整体计算显然是办不到的,本文所提出的分区解法,使用了静力凝聚的概念,将一个整体的定解问题转化为若干个彼此有联系的分区定解问题,这样不仅解决了容量问题,而且提高了计算效率。     

Percentage problems in bridging courses

Research on teaching high school mathematics shows that the topic of percentages often causes learning difficulties. This article describes a method of teaching percentages that the authors used in university bridging courses. In this method, the information from a word problem about percentages is presented in a two-way table. Such a table gives a logical structure to the problem and provides an algorithm for finding a simple equation for the unknown value of interest. The use of this procedure is illustrated by several examples of different levels of difficulty. The method can be applied to many types of percentage problems, so it is quite universal.     

A starting method for the numerical solution of Volterra's integral equation of the second kind

The acquisition of starting values is one of the chief difficulties encountered in computing a numerical solution of Volterra's integral equation of the second kind by a multi-step method. The object of this note is to present a procedure which is derived from certain quadrature formulas and which provides these starting values, to provide a sufficient condition for the approximate solution to be unique, to bound the approximate solution and the error, and to give a numerical example.     

Symmetry of the blow-up set of a porous medium type equation

We study the blow-up set of a porous medium type equation with source. Under some technical conditions, we prove that if the blow-up set is a bounded smooth region, then it must be a ball with a certain radius. This problem can be reduced to a sublinear elliptic equation coupled with an overdetermined boundary condition. Roughly speaking, the overdetermined boundary condition forces the domain to be a ball. Because the nonlinear term is sublinear and then non-Lipschitz, many difficulties arise if one wants to use the moving plane method to reach the goal. In particular, the Hopf boundary lemma is not applicable to this problem. Instead, we investigate various related problems in a half space and a problem in the first quadrant of the entire space, and then use the symmetry results obtained for these problems to overcome the obstacles encountered. ©1995 John Wiley & Sons, Inc.     

Symmetry and uniqueness of minimizers of Hartree type equations with external Coulomb potential

In the present article we study the radial symmetry and uniqueness of minimizers of the energy functional, corresponding to the repulsive Hartree equation in external Coulomb potential. To overcome the difficulties, resulting from the “bad” sign of the nonlocal term, we modify the reflection method and obtain symmetry and uniqueness results.     

Local Existence of Analytical Solutions to an Incompressible Lagrangian Stochastic Model in a Periodic Domain

We consider an incompressible kinetic Fokker Planck equation in the flat torus, which is a simplified version of the Lagrangian stochastic models for turbulent flows introduced by S.B. Pope in the context of computational fluid dynamics. The main difficulties in its treatment arise from a pressure type force that couples the Fokker Planck equation with a Poisson equation which strongly depends on the second order moments of the fluid velocity. In this paper we prove short time existence of analytic solutions in the one-dimensional case, for which we are able to use techniques and functional norms that have been recently introduced in the study of a related singular model.     

On a boundary variational inequality of the second kind modelling a friction problem

A scalar contact problem with friction governed by the Yukawa equation is reduced to a boundary variational inequality. The presence of the non‐differentiable friction functional causes some difficulties when approximated. We present two methods to overcome this difficulty. The first one is a regularization leading to a non‐linear boundary variational equation, for which we propose an iterative procedure, whereas the second method is based on the boundary mixed variational formulation involving Lagrange multipliers. We propose Uzawa's algorithm to compute the saddle point of the corresponding boundary Lagrangian and investigate the discretization of various formulations by the boundary element Galerkin method. Convergence of the boundary element solution is proved and a convergence order is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Detection and construction of an elliptic solution of the complex cubic-quintic Ginzburg-Landau equation

In evolution equations for a complex amplitude, the equation for the phase is much more intricate than for the amplitude. Nevertheless, general methods should be applicable to both variables. In the example of the traveling-wave reduction of the complex cubic-quintic Ginzburg-Landau (CGL5) equation, we explain how to overcome the difficulties arising in two methods: (1) the criterion that the sum of residues of an elliptic solution is zero and (2) the construction of a first-order differential equation admitting a given equation as a differential consequence (subequation method).     

A note on the wave propagation in water of variable depth

In the present work, utilizing the two dimensional equations of an incompressible inviscid fluid and the reductive perturbation method we studied the propagation of weakly nonlinear waves in water of variable depth. For the case of slowly varying depth, the evolution equation is obtained as the variable coefficient Korteweg-de Vries (KdV) equation. Due to the difficulties for the analytical solutions, a numerical technics so called “the method of integrating factor” is used and the evolution equation is solved under a given initial condition and the bottom topography. It is observed the parameters of bottom topography causes to the changes in wave amplitude, wave profile and the wave speed.     

A NEW FINITE ELEMENT APPROXIMATION OF A STATE-CONSTRAINED OPTIMAL CONTROL PROBLEM

In this paper, we study numerical methods for an optimal control problem with pointwise state constraints. The traditional approaches often need to deal with the deltasingularity in the dual equation, which causes many difficulties in its theoretical analysis and numerical approximation. In our new approach we reformulate the state-constrained optimal control as a constrained minimization problems only involving the state, whose optimality condition is characterized by a fourth order elliptic variational inequality. Then direct numerical algorithms （nonconforming finite element approximation） are proposed for the inequality, and error estimates of the finite element approximation are derived. Numerical experiments illustrate the effectiveness of the new approach.