The Augmented Scattering Matrix and Exponentially Decaying Solutions of an Elliptic Problem in a Cylindrical Domain

The self-adjoint elliptic boundary-value problem in a domain with cylindrical outlets to infinity is considered. The notion of an augmented scattering matrix is introduced on the basis of artificial radiation conditions. Properties of the augmented scattering matrix are studied, and the relationship with the classical scattering matrix is demonstrated. The central point is the possibility of calculating the number of linearly independent solutions of a homogeneous problem with fixed rate of decrease at infinity by analyzing the spectrum of the augmented scattering matrix. This property is applied to the problem on diffraction on a periodic boundary as an example. Bibliography: 21 titles.     

A discontinuous Galerkin method for Stokes equation by divergence-free patch reconstruction

A discontinuous Galerkin method by patch reconstruction is proposed for Stokes flows. A locally divergence-free reconstruction space is employed as the approximation space, and the interior penalty method is adopted which imposes the normal component penalty terms to cancel out the pressure term. Consequently, the Stokes equation can be solved as an elliptic system instead of a saddle-point problem due to such weak form. The number of degree of freedoms of our method is the same as the number of elements in the mesh for different order of accuracy. The error estimations of the proposed method are given in a classical style, which are then verified by some numerical examples.     

On a Toeplitz Representation of the H1/2 Norm

Some representations of the H1/2 norm are used as Schur complement preconditioner in PCG based domain decomposition algorithms for elliptic problems. These norm representations are efficient preconditioners but the corresponding matrices are dense, so they need FFT algorithm for matrix-vector multiplications. Here we give a new matrix representation of this norm by a special Toeplitz matrix. It contains only O(log(n)) different entries at each row, where n is the number of rows and so a matrix-vector computation can be done by O(nlog(n)) arithmetic operation without using FFT algorithm. The special properties of this matrix assure that it can be used as preconditioner. This is proved by estimating spectral equivalence constants and this fact has also been verified by numerical tests.     

Optimal interior partial regularity¶for nonlinear elliptic systems: the method of A-harmonic approximation

We consider nonlinear elliptic systems of divergence type. We provide a new method for proving partial regularity for weak solutions, based on a generalization of the technique of harmonic approximation. This method is applied to both homogeneous and inhomogeneous systems, in the latter case with inhomogeneity obeying the natural growth condition. Our methods extend previous partial regularity results, directly establishing the optimal H?lder exponent for the derivative of a weak solution on its regular set. We also indicate how the technique can be applied to further simplify the proof of partial regularity for quasilinear elliptic systems. Received: 22 July 1999 / Revised version: 23 May 2000     

Discrete and continuous maximum principles for parabolic and elliptic operators

The major qualitative properties of linear parabolic and elliptic operators/PDEs are the different maximum principles (MPs). Another important property is the stabilization property (SP), which connects these two types of operators/PDEs. This means that under some assumptions the solution of the parabolic PDE tends to an equilibrium state when t→∞, which is the solution of the corresponding elliptic PDE. To solve PDEs we need to use some numerical methods, and it is a natural requirement that these qualitative properties are preserved on the discrete level. In this work we investigate this question when a two-level discrete mesh operator is used as the discrete model of the parabolic operator (which is a one-step numerical procedure for solving the parabolic PDE) and a matrix as a discrete elliptic operator (which is a linear algebraic system of equations for solving the elliptic PDE). We clarify the relation between the discrete parabolic maximum principle (DPMP), the discrete elliptic maximum principle (DEMP) and the discrete stabilization property (DSP). The main result is that the DPMP implies the DSP and the DEMP.     

On the convergence of a dual-primal substructuring method

Summary.   In the Dual-Primal FETI method, introduced by Farhat et al. [5], the domain is decomposed into non-overlapping subdomains, but the degrees of freedom on crosspoints remain common to all subdomains adjacent to the crosspoint. The continuity of the remaining degrees of freedom on subdomain interfaces is enforced by Lagrange multipliers and all degrees of freedom are eliminated. The resulting dual problem is solved by preconditioned conjugate gradients. We give an algebraic bound on the condition number, assuming only a single inequality in discrete norms, and use the algebraic bound to show that the condition number is bounded by for both second and fourth order elliptic selfadjoint problems discretized by conforming finite elements, as well as for a wide class of finite elements for the Reissner-Mindlin plate model. Received January 20, 2000 / Revised version received April 25, 2000 / Published online December 19, 2000     

Green's functions for elliptic and parabolic systems with Robin-type boundary conditions

The aim of this paper is to investigate Green's function for parabolic and elliptic systems satisfying a possibly nonlocal Robin-type boundary condition. We construct Green's function for parabolic systems with time-dependent coefficients satisfying a possibly nonlocal Robin-type boundary condition assuming that weak solutions of the system are locally Hölder continuous in the interior of the domain, and as a corollary we construct Green's function for elliptic system with a Robin-type condition. Also, we obtain Gaussian bound for Robin Green's function under an additional assumption that weak solutions of Robin problem are locally bounded up to the boundary. We provide some examples satisfying such a local boundedness property, and thus have Gaussian bounds for their Green's functions.     

Lyapunov diagonal semistability of acyclic matrices

It is shown that an acyclic matrix is Lyapunov diagonally semistable if and only if the matrix has the weak principal submatrix rank property. This result completes the solution of the problem of characterizing the various types of matrix stability for acyclic matrices. Also, those acyclic matrices which have a unique Lyapunov scaling factor are characterized.     

Regularity of Solutions to Problems for Nondiagonal Quasilinear Elliptic Systems with Conditions on Interface Surfaces

We consider a model problem for a quasilinear elliptic system with nondiagonal principal matrix and certain conditions on a surface dividing a medium. The partial regularity of a weak solution in a neighborhood of the interface is established. Bibliography: 11 titles.     

Convergence of the ghost fluid method for elliptic equations with interfaces

This paper proves the convergence of the ghost fluid method for second order elliptic partial differential equations with interfacial jumps. A weak formulation of the problem is first presented, which then yields the existence and uniqueness of a solution to the problem by classical methods. It is shown that the application of the ghost fluid method by Fedkiw, Kang, and Liu to this problem can be obtained in a natural way through discretization of the weak formulation. An abstract framework is given for proving the convergence of finite difference methods derived from a weak problem, and as a consequence, the ghost fluid method is proved to be convergent.

     

On the Non‐Hierarchical Preconditioner Constructions

In this paper a new cyclic matrix representation of the H^–1/2 norm is presented. Its application as Schur complement preconditioning matrix requires only matrix‐vector multiplication. The computational cost of this matrix‐vector multiplication is O (N · log(N)) arithmetic operations, where N is the number of unknowns. The efficiency of the construction to elliptic problems has been verified by numerical tests. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a class of degenerate elliptic operators arising from Fleming-Viot processes

We are dealing with the solvability of an elliptic problem related to a class of degenerate second order operators which arise from the theory of Fleming-Viot processes in population genetics. In the one dimensional case the problem is solved in the space of continuous functions. In higher dimension we study the problem in spaces with respect to an explicit measure which, under suitable assumptions, can be taken invariant and symmetrizing for the operators. We prove the existence and uniqueness of weak solutions and we show that the closure of the operator in such spaces generates an analytic -semigroup. Received December 4, 2000; accepted December 9, 2000.     

Supersonic flow onto a solid wedge

We consider the problem of two‐dimensional supersonic flow onto a solid wedge, or equivalently in a concave corner formed by two solid walls. For mild corners, there are two possible steady state solutions, one with a strong and one with a weak shock emanating from the corner. The weak shock is observed in supersonic flights. A longstanding natural conjecture is that the strong shock is unstable in some sense. We resolve this issue by showing that a sharp wedge will eventually produce weak shocks at the tip when accelerated to a supersonic speed. More precisely, we prove that for upstream state as initial data in the entire domain, the time‐dependent solution is self‐similar, with a weak shock at the tip of the wedge. We construct analytic solutions for self‐similar potential flow, both isothermal and isentropic with arbitrary γ ≥ 1. In the process of constructing the self‐similar solution, we develop a large number of theoretical tools for these elliptic regions. These tools allow us to establish large‐data results rather than a small perturbation. We show that the wave pattern persists as long as the weak shock is supersonic‐supersonic; when this is no longer true, numerics show a physical change of behavior. In addition, we obtain rather detailed information about the elliptic region, including analyticity as well as bounds for velocity components and shock tangents. © 2007 Wiley Periodicals, Inc.     

A two-well structure and intrinsic mountain pass points

Under small dead-load perturbations, and the natural boundary value condition (Neumann problem), we establish the existence of an unstable critical point (mountain pass point) for a variational integral with a two-well structure. The integrands we consider are obtained by the quasiconvex relaxation [18] of the squared distance function and its quasiconvex lower bounds. The models are m otivated by the variational approach to material microstructure when the wells are incompatible. We show that these functions give quasimonotone gradient mappings. We introduce the weak Palais-Smale condition (weak PS) to deal with the lack of compactness in the borderline case where the integrand is . Received March 1, 1999 / Accepted March 29, 2000 / Published online December 8, 2000     

Verification of Second-Order Sufficient Optimality Conditions for Semilinear Elliptic and Parabolic Control Problems

We study optimal control problems for semilinear parabolic equations subject to control constraints and for semilinear elliptic equations subject to control and state constraints. We quote known second-order sufficient optimality conditions (SSC) from the literature. Both problem classes, the parabolic one with boundary control and the elliptic one with boundary or distributed control, are discretized by a finite difference method. The discrete SSC are stated and numerically verified in all cases providing an indication of optimality where only necessary conditions had been studied before.     

椭圆边界上的自然积分算子及各向异性外问题的耦合算法

1.引 言为求解微分方程的外边值问题常需要引进人工边界（见[1-4]）,对人工边界外部区域作自然边界归化得到的自然积分方程即Dirichlet-Neumann映射,正是人工边界上的准确的边界条件（见[2-6]）,这是一类非局部边界条件.自然积分算子即Dirichlet-Neumann算子,     

Discrete boundary element methods on general meshes in 3D

Abstract. This paper is concerned with the stability and convergence of fully discrete Galerkin methods for boundary integral equations on bounded piecewise smooth surfaces in . Our theory covers equations with very general operators, provided the associated weak form is bounded and elliptic on , for some . In contrast to other studies on this topic, we do not assume our meshes to be quasiuniform, and therefore the analysis admits locally refined meshes. To achieve such generality, standard inverse estimates for the quasiuniform case are replaced by appropriate generalised estimates which hold even in the locally refined case. Since the approximation of singular integrals on or near the diagonal of the Galerkin matrix has been well-analysed previously, this paper deals only with errors in the integration of the nearly singular and smooth Galerkin integrals which comprise the dominant part of the matrix. Our results show how accurate the quadrature rules must be in order that the resulting discrete Galerkin method enjoys the same stability properties and convergence rates as the true Galerkin method. Although this study considers only continuous piecewise linear basis functions on triangles, our approach is not restricted in principle to this case. As an example, the theory is applied here to conventional “triangle-based” quadrature rules which are commonly used in practice. A subsequent paper [14] introduces a new and much more efficient “node-based” approach and analyses it using the results of the present paper. Received December 10, 1997 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Superrelaxation and the rate of convergence in minimizing quadratic functions subject to bound constraints

The paper resolves the problem concerning the rate of convergence of the working set based MPRGP (modified proportioning with reduced gradient projection) algorithm with a long steplength of the reduced projected gradient step. The main results of this paper are the formula for the R-linear rate of convergence of MPRGP in terms of the spectral condition number of the Hessian matrix and the proof of the finite termination property for the problems whose solution does not satisfy the strict complementarity condition. The bound on the R-linear rate of convergence of the projected gradient is also included. For shorter steplengths these results were proved earlier by Dostál and Schöberl. The efficiency of the longer steplength is illustrated by numerical experiments. The result is an important ingredient in developming scalable algorithms for numerical solution of elliptic variational inequalities and substantiates the choice of parameters that turned out to be effective in numerical experiments.     

On one approach to a posteriori error estimates for evolution problems solved by the method of lines

Summary. In this paper, we describe a new technique for a posteriori error estimates suitable to parabolic and hyperbolic equations solved by the method of lines. One of our goals is to apply known estimates derived for elliptic problems to evolution equations. We apply the new technique to three distinct problems: a general nonlinear parabolic problem with a strongly monotonic elliptic operator, a linear nonstationary convection-diffusion problem, and a linear second order hyperbolic problem. The error is measured with the aid of the -norm in the space-time cylinder combined with a special time-weighted energy norm. Theory as well as computational results are presented. Received September 2, 1999 / Revised version received March 6, 2000 / Published online March 20, 2001