Analysis of a BDF–DGFE scheme for nonlinear convection–diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection–diffusion equation. We employ a combination of the discontinuous Galerkin finite element method for the space semi-discretization and the k-step backward difference formula for the time discretization. The diffusive and stabilization terms are treated implicitly whereas the nonlinear convective term is treated by a higher order explicit extrapolation method, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the discrete L ^∞(L ²)-norm and the L ²(H ¹)-seminorm with respect to the mesh size h and time step τ for k = 2,3. Numerical examples verifying the theoretical results are presented. This work is a part of the research project MSM 0021620839 financed by the Ministry of Education of the Czech Republic and was partly supported by the Grant No. 316/2006/B-MAT/MFF of the Grant Agency of the Charles University Prague. The research of M. Vlasák was supported by the project LC06052 of the Ministry of Education of the Czech Republic (Jindřich Nečas Center for Mathematical Modelling).     

Backward Euler method for abstract time-dependent parabolic equations with variable domains

We consider semidiscretizations in time, based on the backward Euler method, of an abstract, non-autonomous parabolic initial value problem where , , is a family of sectorial operators in a Banach space X. The domains are allowed to depend on t. Our hypotheses are fulfilled for classical parabolic problems in the , , norms. We prove that the semidiscretization is stable in a suitable sense. We get optimal estimates for the error even when non-homogeneous boundary values are considered. In particular, the results are applicable to the analysis of the semidiscretizations of time-dependent parabolic problems under non-homogeneous Neumann boundary conditions. Received October 17, 1997 / Revised version received April 17, 1998     

Discrete approximations of BV solutions to doubly nonlinear degenerate parabolic equations

Summary. In this paper we present and analyse certain discrete approximations of solutions to scalar, doubly nonlinear degenerate, parabolic problems of the form under the very general structural condition . To mention only a few examples: the heat equation, the porous medium equation, the two-phase flow equation, hyperbolic conservation laws and equations arising from the theory of non-Newtonian fluids are all special cases of (P). Since the diffusion terms a(s) and b(s) are allowed to degenerate on intervals, shock waves will in general appear in the solutions of (P). Furthermore, weak solutions are not uniquely determined by their data. For these reasons we work within the framework of weak solutions that are of bounded variation (in space and time) and, in addition, satisfy an entropy condition. The well-posedness of the Cauchy problem (P) in this class of so-called BV entropy weak solutions follows from a work of Yin [18]. The discrete approximations are shown to converge to the unique BV entropy weak solution of (P). Received November 10, 1998 / Revised version received June 10, 1999 / Published online June 8, 2000     

Convergence estimates of a projection-difference method for an operator-differential equation

This article investigates the projection-difference method for a Cauchy problem for a linear operator-differential equation with a leading self-adjoint operator A(t) and a subordinate linear operator K(t) in Hilbert space. This method leads to the solution of a system of linear algebraic equations on each time level; moreover, the projection subspaces are linear spans of eigenvectors of an operator similar to A(t). The convergence estimates are obtained. The application of the developed method for solving the initial boundary value problem is given.     

Finite element analysis of exponentially fitted Lumped schemes for time-dependent convection-diffusion problems

Summary In the paper we consider a singularly perturbed linear parabolic initialboundary value problem in one space variable. Two exponential fitted schemes are derived for the problem using Petrov-Galerkin finite element methods with various choices of trial and test spaces. On rectangular meshes which are either arbitrary or slightly restricted, we derive global energy norm andL ² norm and localL error bounds which are uniform in the diffusion parameter. Numerical results are also persented.     

The Dirichlet problem of the quantum drift-diffusion model

We study the quantum drift-diffusion model, a fourth-order parabolic system, with Dirichlet boundary conditions. Using a semi-discretization approximate method with a compact argument and applying a new entropy estimate, we prove the existence of global regular weak solutions.     

Co-volume methods for degenerate parabolic problems

Summary A complementary volume (co-volume) technique is used to develop a physically appealing algorithm for the solution of degenerate parabolic problems, such as the Stefan problem. It is shown that, these algorithms give rise to a discrete semigroup theory that parallels the continuous problem. In particular, the discrete Stefan problem gives rise to nonlinear semigroups in both the discreteL ¹ andH ^–1 spaces.The first author was supported by a grant from the Hughes foundation, and the second author was supported by the National Science Foundation Grant No. DMS-9002768 while this work was undertaken. This work was supported by the Army Research Office and the National Science Foundation through the Center for Nonlinear Analysis.     

Galerkin and subspace decomposition methods in space and time for the Navier-Stokes equations

The Galerkin method and the subspace decomposition method in space and time for the two-dimensional incompressible Navier-Stokes equations with the H²-initial data are considered. The subspace decomposition method consists of splitting the approximate solution as the sum of a low frequency component discretized by the small time step Δt and a high frequency one discretized by the large time step pΔt with p>1. The H²-stability and L²-error analysis for the subspace decomposition method are obtained. Finally, some numerical tests to confirm the theoretical results are provided.     

Spectral element methods for parabolic problems

A spectral element method for solving parabolic initial boundary value problems on smooth domains using parallel computers is presented in this paper. The space domain is divided into a number of shape regular quadrilaterals of size h and the time step k   is proportional to h²$h^2$. At each time step we minimize a functional which is the sum of the squares of the residuals in the partial differential equation, initial condition and boundary condition in different Sobolev norms and a term which measures the jump in the function and its derivatives across inter-element boundaries in certain Sobolev norms. The Sobolev spaces used are of different orders in space and time. We can define a preconditioner for the minimization problem which allows the problem to decouple. Error estimates are obtained for both the h and p versions of this method.     

Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem

The Generalized Riemann Problem (GRP) for a nonlinear hyperbolic system of m balance laws (or alternatively “quasi-conservative” laws) in one space dimension is now well-known and can be formulated as follows: Given initial-data which are analytic on two sides of a discontinuity, determine the time evolution of the solution at the discontinuity. In particular, the GRP numerical scheme (second-order high resolution) is based on an analytical evaluation of the first time derivative. It turns out that this derivative depends only on the first-order spatial derivatives, hence the initial data can be taken as piecewise linear. The analytical solution is readily obtained for a single equation (m = 1) and, more generally, if the system is endowed with a complete (coordinate) set of Riemann invariants. In this case it can be “diagonalized” and reduced to the scalar case. However, most systems with m > 2 do not admit such a set of Riemann invariants. This paper introduces a generalization of this concept: weakly coupled systems (WCS). Such systems have only “partial set” of Riemann invariants, but these sets are weakly coupled in a way which enables a “diagonalized” treatment of the GRP. An important example of a WCS is the Euler system of compressible, nonisentropic fluid flow (m = 3). The solution of the GRP discussed here is based on a careful analysis of rarefaction waves. A “propagation of singularities” argument is applied to appropriate Riemann invariants across the rarefaction fan. It serves to “rotate” initial spatial slopes into “time derivative”. In particular, the case of a “sonic point” is incorporated easily into the general treatment. A GRP scheme based on this solution is derived, and several numerical examples are presented. Special attention is given to the “acoustic approximation” of the analytical solution. It can be viewed as a proper linearization (different from the approach of Roe) of the nonlinear system. The resulting numerical scheme is the simplest (second-order, high-resolution) generalization of the Godunov scheme.     

Semi-discretization in time and numerical convergence of solutions of a nonlinear cross-diffusion population model

Summary. A positivity-preserving numerical scheme for a strongly coupled cross-diffusion model for two competing species is presented, based on a semi-discretization in time. The variables are the population densities of the species. Existence of strictly positive weak solutions to the semidiscrete problem is proved. Moreover, it is shown that the semidiscrete solutions converge to a non-negative solution of the continuous system in one space dimension. The proofs are based on a symmetrization of the problem via an exponential transformation of variables and the use of an entropy functional. Received September 10, 2001 / Revised version received February 25, 2002 / Published online June 17, 2002     

A second-order accurate numerical method for a fractional wave equation

We study a generalized Crank–Nicolson scheme for the time discretization of a fractional wave equation, in combination with a space discretization by linear finite elements. The scheme uses a non-uniform grid in time to compensate for the singular behaviour of the exact solution at t = 0. With appropriate assumptions on the data and assuming that the spatial domain is convex or smooth, we show that the error is of order k ² + h ², where k and h are the parameters for the time and space meshes, respectively.     

Solving the generalized regularized long wave equation on the basis of a reproducing kernel space

On the basis of a reproducing kernel space, an iterative algorithm for solving the generalized regularized long wave equation is presented. The analytical solution in the reproducing kernel space is shown in a series form and the approximate solution u_n is constructed by truncating the series to n terms. The convergence of u_n to the analytical solution is also proved. Results obtained by the proposed method imply that it can be considered as a simple and accurate method for solving such evolution equations.     

Space–time discontinuous Galerkin finite element method for shallow water flows

A space–time discontinuous Galerkin (DG) finite element method is presented for the shallow water equations over varying bottom topography. The method results in nonlinear equations per element, which are solved locally by establishing the element communication with a numerical HLLC flux. To deal with spurious oscillations around discontinuities, we employ a dissipation operator only around discontinuities using Krivodonova's discontinuity detector. The numerical scheme is verified by comparing numerical and exact solutions, and validated against a laboratory experiment involving flow through a contraction. We conclude that the method is second order accurate in both space and time for linear polynomials.     

Error growth analysis via stability regions for discretizations of initial value problems

This paper deals with numerical methods for the solution of linear initial value problems. Two main theorems are presented on the stability of these methods. Both theorems give conditions guaranteeing a mild error growth, for one-step methods characterized by a rational function ϕ(z). The conditions are related to the stability regionS={z:z∈ℂ with |ϕ(z)|≤1}, and can be viewed as variants to the resolvent condition occurring in the reputed Kreiss matrix theorem. Stability estimates are presented in terms of the number of time stepsn and the dimensions of the space. The first theorem gives a stability estimate which implies that errors in the numerical process cannot grow faster than linearly withs orn. It improves previous results in the literature where various restrictions were imposed onS and ϕ(z), including ϕ′(z)≠0 forz∈σS andS be bounded. The new theorem is not subject to any of these restrictions. The second theorem gives a sharper stability result under additional assumptions regarding the differential equation. This result implies that errors cannot grow faster thann ^β, with fixed β<1. The theory is illustrated in the numerical solution of an initial-boundary value problem for a partial differential equation, where the error growth is measured in the maximum norm.     

A numerical method for the Benjamin-Ono equation

This paper is concerned with the numerical solution of the Cauchy problem for the Benjamin-Ono equationu _t +uu _x −Hu _xx =0, whereH denotes the Hilbert transform. Our numerical method first approximates this Cauchy problem by an initial-value problem for a corresponding 2L-periodic problem in the spatial variable, withL large. This periodic problem is then solved using the Crank-Nicolson approximation in time and finite difference approximations in space, treating the nonlinear term in a standard conservative fashion, and the Hilbert transform by a quadrature formula which may be computed efficiently using the Fast Fourier Transform.     

Convergence rates of iterated Tikhonov regularized solutions of nonlinear III — posed problems

Summary In this paper we investigate iterated Tikhonov regularization for the solution of nonlinear ill-posed problems. In the case of linear ill-posed problems it is well-known that (under appropriate assumptions) then-th iterated regularized solutions can converge likeO(²² /(2n+1)), where denotes the noise level of the data perturbation. We give conditions that guarantee this convergence rate also for nonlinear ill-posed problems, and motivate these conditions by the mapping degree. The results are derived by a comparison of the iterated regularized solutions of the nonlinear problem with the iterated regularized solutions of its linearization. Numerical examples are presented.Supported by the Austrian Fonds zur Förderung der wissenschaftlichen Forschung,project P-7869 PHY, and by the Christian Doppler Society     

The Fujita exponent for the Cauchy problem in the hyperbolic space

It is well known that the heat kernel in the hyperbolic space has a different behavior for large times than the one in the Euclidean space. The main purpose of this paper is to study its effect on the positive solutions of Cauchy problems with power nonlinearities. Existence and non-existence results for local solutions are derived. Emphasis is put on their long time behavior and on Fujita?s phenomenon. To have the same situation as for the Cauchy problem in RN, namely finite time blow up for all solutions if the exponent is smaller than a critical value and existence of global solutions only for powers above the critical exponent, we must introduce a weight depending exponentially on the time. In this respect the situation is similar to problems in bounded domains with Dirichlet boundary conditions. Important tools are estimates for the heat kernel in the hyperbolic space and comparison principles.     

Multi-dimensional scalar balance laws with discontinuous flux

We consider scalar balance laws with a dissipative source term. The flux function may be discontinuous with respect to both the space variable x and the unknown quantity u. We formulate the definition of entropy weak solutions and provide existence and uniqueness to the considered problem. The problem is formulated in the framework of multi-valued mappings. The notion of entropy measure-valued solutions is used to prove the so-called contraction principle and comparison principle.     

Conservation law with the flux function discontinuous in the space variable—II: Convex–concave type fluxes and generalized entropy solutions

We deal with a single conservation law with discontinuous convex–concave type fluxes which arise while considering sign changing flux coefficients. The main difficulty is that a weak solution may not exist as the Rankine–Hugoniot condition at the interface may not be satisfied for certain choice of the initial data. We develop the concept of generalized entropy solutions for such equations by replacing the Rankine–Hugoniot condition by a generalized Rankine–Hugoniot condition. The uniqueness of solutions is shown by proving that the generalized entropy solutions form a contractive semi-group in L¹$L^1$. Existence follows by showing that a Godunov type finite difference scheme converges to the generalized entropy solution. The scheme is based on solutions of the associated Riemann problem and is neither consistent nor conservative. The analysis developed here enables to treat the cases of fluxes having at most one extrema in the domain of definition completely. Numerical results reporting the performance of the scheme are presented.