A kinetic approach to comparison properties for degenerate parabolic-hyperbolic equations with boundary conditions

We study the comparison principle for degenerate parabolic-hyperbolic equations with initial and nonhomogeneous boundary conditions. We prove a comparison theorem for any entropy sub- and supersolution. The L¹ contractivity and, therefore, uniqueness of entropy solutions has been obtained so far by some authors, but it seems that any comparison theorem is not proven. The method used there is the doubling variable technique due to Kru?kov. Our method is based upon the kinetic formulation and the kinetic techniques. By developing the kinetic techniques for degenerate parabolic-hyperbolic equations with boundary conditions, we can obtain a comparison property which obviously extends the L¹ contractive property.     

Nonstationary boundary problem for model kinetic equations at critical parameters

Nonstationary solutions of the model kinetic equation at critical values of the motion of the wall (the boundary of the half-space occupied by gas) are studied. The characteristic equation is obtained by separating the variables. The eigenfunctions and the eigenvalue spectrum are found in the distribution space. A solution to the equation is expandable over the eigenfunction basis. The Rayleigh problem is considered as an application. The distribution function is continuous in the plane of the wall-motion parameters, including the closed curve of critical parameter values. Translated from Teoreticheskaya i Matematicheskaya Fizika. Vol. 116. No. 2, pp. 305–320. August. 1998.     

Fast parallel solvers for symmetric boundary element domain decomposition equations

Summary. The boundary element method (BEM) is of advantage in many applications including far-field computations in magnetostatics and solid mechanics as well as accurate computations of singularities. Since the numerical approximation is essentially reduced to the boundary of the domain under consideration, the mesh generation and handling is simpler than, for example, in a finite element discretization of the domain. In this paper, we discuss fast solution techniques for the linear systems of equations obtained by the BEM (BE-equations) utilizing the non-overlapping domain decomposition (DD). We study parallel algorithms for solving large scale Galerkin BE–equations approximating linear potential problems in plane, bounded domains with piecewise homogeneous material properties. We give an elementary spectral equivalence analysis of the BEM Schur complement that provides the tool for constructing and analysing appropriate preconditioners. Finally, we present numerical results obtained on a massively parallel machine using up to 128 processors, and we sketch further applications to elasticity problems and to the coupling of the finite element method (FEM) with the boundary element method. As shown theoretically and confirmed by the numerical experiments, the methods are of algebraic complexity and of high parallel efficiency, where denotes the usual discretization parameter. Received August 28, 1996 / Revised version received March 10, 1997     

Convergence of alternating domain decomposition schemes for kinetic and aerodynamic equations

A domain decomposition scheme linking linearized kinetic and aerodynamic equations is considered. Convergence of the alternating scheme is shown. Moreover the physical correctness of the obtained coupled solutions is discussed.     

A space decomposition method for parabolic equations

A convergence proof is given for an abstract parabolic equation using general space decomposition techniques. The space decomposition technique may be a domain decomposition method, a multilevel method, or a multigrid method. It is shown that if the Euler or Crank–Nicolson scheme is used for the parabolic equation, then by suitably choosing the space decomposition, only O(| log τ |) steps of iteration at each time level are needed, where τ is the time-step size. Applications to overlapping domain decomposition and to a two-level method are given for a second-order parabolic equation. The analysis shows that only a one-element overlap is needed. Discussions about iterative and noniterative methods for parabolic equations are presented. A method that combines the two approaches and utilizes some of the good properties of the two approaches is tested numerically. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 27–46, 1998     

A singularity-based eigenfunction decomposition for Kohn-Sham equations

We show that the eigenfunctions of Kohn-Sham equations can be decomposed as ■ = F ψ, where F depends on the Coulomb potential only and is locally Lipschitz, while ψ has better regularity than ■.     

A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations

In this paper we propose a new modified recursion scheme for the resolution of multi-order and multi-point boundary value problems for nonlinear ordinary and partial differential equations by the Adomian decomposition method (ADM). Our new approach, including Duan’s convergence parameter, provides a significant computational advantage by allowing for the acceleration of convergence and expansion of the interval of convergence during calculations of the solution components for nonlinear boundary value problems, in particular for such cases when one of the boundary points lies outside the interval of convergence of the usual decomposition series. We utilize the boundary conditions to derive an integral equation before establishing the recursion scheme for the solution components. Thus we can derive a modified recursion scheme without any undetermined coefficients when computing successive solution components, whereas several prior recursion schemes have done so. This modification also avoids solving a sequence of nonlinear algebraic equations for the undetermined coefficients fraught with multiple roots, which is required to complete calculation of the solution by several prior modified recursion schemes using the ADM.     

A boundary only procedure for time-dependent diffusion equations

In the recent literature, the boundary element method (BEM) is extensively used to solve time-dependent partial differential equations. However, most of these formulations yield algorithms where one has to include all interior points in the computation process if finite difference procedures are used to approximate the temporal derivative. This obviously restricts the advantages of the BEM, which is mainly considered to be a boundary only algorithm for time-independent problems. A new algorithm is demonstrated here, which extends the boundary only nature of the method to time-dependent partial differential equations. Using this procedure, one can reduce the finite difference time integration algorithm, generated in a standard manner, to a boundary only process. The proposed method is demonstrated with considerable success for diffusion problems. Results obtained in these applications are presented comparatively with analytical and other boundary element time integration procedures. The algorithm proposed may utilize several coordinate functions in the secondary reduction phase of the formulation. A summary of such functions is described here and performances of these functions are tested and compared in three applications. It is shown that some coordinate functions perform better than others under certain conditions. Using these results, we propose a general coordinate function, which may be used with satisfactory results in all parabolic partial differential equation applications.     

A boundary multiplier/fictitious domain method for the steady incompressible Navier-Stokes equations

Summary. We analyze the error of a fictitious-domain method with boundary Lagrange multiplier. It is applied to solve a non-homogeneous steady incompressible Navier-Stokes problem in a domain with a multiply-connected boundary. The interior mesh in the fictitious domain and the boundary mesh are independent, up to a mesh-length ratio. Received February 24, 1999 / Revised version received January 30, 2000 / Published online October 16, 2000     

A kinetic scheme for the Savage–Hutter equations

The Savage–Hutter (SH) equations describe the motion of granular material under the influence of friction. Based on the kinetic formulation of the SH equations, we present a kinetic scheme in one dimension, which describes the deformation of the mass profile and allows it to start and to stop. Moreover the method is able to preserve the steady states of granular masses at rest. The method is tested on several numerical examples. Copyright © 2008 John Wiley & Sons, Ltd.     

A kinetic decomposition for singular limits of non-local conservation laws

We consider a non-local regularization of nonlinear hyperbolic conservation laws in several space variables. The regularization is motivated by the theory of phase dynamics and is based on a convolution operator. We formulate the initial value problem and begin by deriving a priori estimates which are independent of the regularization parameter. Following Hwang and Tzavaras we establish a kinetic decomposition associated with the problem under consideration, and we conclude that the sequence of solutions generated by the non-local model converges to a weak solution of the corresponding hyperbolic problem. Depending on the scaling introduced in the non-local dispersive term, this weak limit is either a classical Kruzkov solution satisfying all entropy inequalities or, more interestingly, a nonclassical entropy solution in the sense defined by LeFloch, that is, a weak solution satisfying a single entropy inequality and containing undercompressive shock waves possibly selected by a kinetic relation. Finally, we illustrate our analytical conclusions with numerical experiments in one spatial variable.     

A Nitsche-based domain decomposition method for hypersingular integral equations

We introduce and analyze a Nitsche-based domain decomposition method for the solution of hypersingular integral equations. This method allows for discretizations with non-matching grids without the necessity of a Lagrangian multiplier, as opposed to the traditional mortar method. We prove its almost quasi-optimal convergence and underline the theory by a numerical experiment.     

A quasi-wavelet algorithm for second kind boundary integral equations

In solving integral equations with a logarithmic kernel, we combine the Galerkin approximation with periodic quasi-wavelet (PQW) [4]. We develop an algorithm for solving the integral equations with only O(N log N) arithmetic operations, where N is the number of knots. We also prove that the Galerkin approximation has a polynomial rate of convergence.     

A boundary value problem for fourth-order elastic beam equations

Multiplicity results for a fourth-order nonlinear boundary value problem are presented. The proof of the main result is based on the critical point theory.