Solution enclosures for scalar conservation laws [Elektronische Ressource] / von Duy Nguyen Vu Hoang

Sujets

• Mathematics

Informations

• Publié par: karlsruher_institut_fur_technologie
• Publié le: 01 janvier 2010
• Nombre de lectures: 205
• Langue: Deutsch

Extrait

Solution Enclosures for Scalar Conservation
Laws
Zur Erlangung des akademischen Grades eines
DOKTORS DER NATURWISSENSCHAFTEN
von der Fakultat¨ fur¨ Mathematik des
Karlsruher Instituts fur¨ Technologie (KIT)
genehmigte
DISSERTATION
von
Dipl. phys. Duy Nguyen Vu Hoang
aus Saarlouis
Tag der mundlichen¨ Prufung:¨ 30.6.2010
Referent: Prof. Dr. M. Plum
Korreferent: Prof. Dr. W. ReichelDanksagung
An dieser Stelle mochte¨ ich ganz herzlich Herrn Prof. Plum danken, der diese Arbeit betreut
hat und mir die Moglichk¨ eiten computergestutzter¨ Beweise in der Mathematik nahegebracht hat.
Ferner mochte¨ ich ihm fur¨ seine zahlreichen Verbesserungsvorschlage,¨ diese Arbeit betreffend,
danken.
Herrn Prof. Reichel danke ich herzlich fur¨ die freundliche Bereitschaft, das Korreferat zu
ubernehmen.¨
Maria Radosz mochte¨ ich fur¨ zahlreiche Diskussionen danken, die zum Gelingen der Arbeit
beigetragen haben, sowie fur¨ die genaue Durchsicht eines (schwierigen) Manuskripts auf Fehler.
Moralische Unterstutzung¨ ist mir von Giacomo K. zuteil geworden, weswegen ihm auch Dank
gebuhrt.¨Contents
1 Introduction 4
1.1 General Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.1 Conservation laws in continuum physics . . . . . . . . . . . . . . . . . . . 4
1.1.2 Hyperbolic conservation laws . . . . . . . . . . . . . . . . . . . . . . . . 7
1.1.3 Computer-assisted proofs for partial differential equations . . . . . . . . . 9
1.1.4 Concluding remarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.2 Formulation of the mathematical problem . . . . . . . . . . . . . . . . . . . . . . 11
1.2.1 Single-shock solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.2.2 Linearization of the problem . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.3 Some remarks on linear hyperbolic problems . . . . . . . . . . . . . . . . 15
1.3 Mathematical preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.3.1 Notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2 Fixed Point Formulation 20
2.1 Some operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.1.1 Fixed point formulation; preliminary considerations . . . . . . . . . . . . 20
2.1.2 Regularized operators. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.2 Compactness and Continuity Properties . . . . . . . . . . . . . . . . . . . . . . . 35
2.3 Main theorem on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3 Solution enclosure 42
3.1 The case of Burgers’ Equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.1.1 Estimation of T and K . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.1.2 Construction ofD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2 The general case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2CONTENTS 3
4 Numerical treatment 52
4.1 Computation of solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
4.1.1 A finite difference scheme . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.1.2 Construction of the smooth approximate solution! . . . . . . . . . . . . . 55
4.1.3 Evaluation of the defects . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.1.4 Verification procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
4.2 Numerical examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
5 Systems of Conservation Laws 73
5.1 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
5.2 Linearization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
5.3 Fixed-point formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
A Linear Hyperbolic Problems 80
A.1 Weak solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
A.2 Smoothing operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
A.3 Existence Theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
1A.3.1 A-priori estimates inL -norm . . . . . . . . . . . . . . . . . . . . . . . . 85
A.3.2 Existence of weak solutions. . . . . . . . . . . . . . . . . . . . . . . . . . 88
B Auxiliary estimations 97
B.1 Estimation ofr (v ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97Chapter 1
Introduction
1.1 General Introduction
1.1.1 Conservation laws in continuum physics
We typically find conservation laws in classical continuum physics, most notably in continuum
mechanics, continuum thermodynamics and electrodynamics. The basic unifying principle in these
disciplines consists in regarding physical systems as being continuously distributed in a region of
space. Thus the physical quantities we use to describe the electromagnetic field or a deformable
body of matter are represented by functions which assign values to every point in space. This point
of view has up to now been very successful in modeling physical processes. It must be noted that
even in the disciplines which take into account the quantified nature of matter (e.g. the quantum
theory of a single electron or the kinetic theory of gases) the basic quantity is a field (the wave
function or the molecular distribution function) in the sense mentioned before.
The basic physical laws are often expressed in the form of conservation laws. This is connected
to the concept of extensive quantities. Extensive physical quantities (like mass, charge, energy,
entropy) are distributed in space: to any region
in space at a certain time, a number Q(t; )
is assigned (the quantity of mass, charge etc. contained in
). A conservation law reflects the
idea that the change of a certain physical extensive quantity Q(t; ) , given by the integral over
its density(x;t) over the region
, is balanced by its fluxj through the boundary of
and the
productionr(x;t):
Z Z Z
d d
Q(t; ) = (x;t)dx = j(x;t)d + r(x;t)dx
dt dt
@

4CHAPTER1. INTRODUCTION 5
In generalQ can be scalar- or vector-valued. Given sufficient smoothness of the fields andj the
above integral form of the conservation law is equivalent to its local form
@
d(x;t) + divj(x;t) =r(x;t); (x;t)2R R:
@t
To arrive at the nonlinear hyperbolic PDE’s which are studied in this work, one has to combine
conservation laws and constitutive relations. In order to illustrate this, we take a look at the basic
structure of mechanics. Here, the motion of a deformable body is described by a function
3 3 :R [0;1)!R
(assumed sufficiently smooth) with the property (X; 0) = X, whose physical interpretation is
as follows: t! (X ;t) is the trajectory of the material point whose position at time t = 0 is0
3X 2 R . This is the ”Lagrangean” or ”referential” description of the continuum, in which we0
3think of X 2 R indexing the individual material points making up the material. Assume for
3 3simplicity that (;t) : R ! R is bijective for t fixed and that the inverse mapping (;t) is
sufficiently smooth.
The basic quantities are the velocity field
@
v(X;t) := (X;t);
@t
the mass density(X;t) (a scalar field) and the stressT (X;t) (a symmetric tensor field), thought
as functions of the particleX. The Eulerian description of the continuum is constructed by con-
^sidering the functionsv^; ;^ T defined by
^v^(x;t) :=v( (x;t);t); ^(x;t) :=( (x;t);t); T (x;t) :=T ( (x;t);t):
Henceforth we shall drop the hat from the notation (as is customary in mechanics) and write e.g.
v(X;t) for the velocity field in the Langrangean description andv(x;t) for the velocity field in the
Eulerian description. The basic laws of continuum mechanics, namely balance of mass, balance
1of momentum and the balance of kinetic energy in Eulerian description are given follows:
@
+ div( v ) = 0
@t
@( v )
+ div( v
v) = divT
@t
1For simplicity, we do not explain the balance of rotational momentum here.CHAPTER1. INTRODUCTION 6
1 @ 12 2jvj + divjvj v = div (Tv):
2@t 2
In the second line, the divergence of the tensor field T is given in cartesian components by
P3
(divT ) = @ T . Using the balance of mass, the second line can be rewritten asi j i;jj=1

@v
+ (vr)v = divT;
@t
a form familiar from fluid dynamics.
The balance laws themselves, which are assumed to hold for all material bodies in the domain
of classical mechanics, are not sufficient to determine the actual time evolution of the system. This
is of course to be expected, since we want to describe a variety of bodies, e.g. elastic and fluid
ones, gases, plasmas etc. within this framework. Therefore we supplement the above by so-called
constitutive relations, which characterize concrete materials (or types of materials). We do not
discuss the general principles on which the choice of the constitutive relations is based. It suffices
to say that a very general form is given symbolically by
dT (X;t) =F(f(Y;t s) :Y2R ; 0stg);
that is, the stress at some particlex at the present time depends through a particular functionalF
on the motion of the continuum at all earlier times.
Further specializing the constitutive relations we arrive at many special mech