Existence and stability of stellardynamic models [Elektronische Ressource] / von Achim Schulze

Description

Existence and stability ofstellardynamic modelsbei der Universit¨at Bayreuthzur Erlangung des Grades einesDoktors der Naturwissenschaften (Dr. rer. nat.)vorgelegte Abhandlung vonAchim Schulzegeb. 11.07.1979 in Amberg1. Gutachter:2. Gutachter:Tag der Einreichung:Tag des Kolloquiums:DanksagungIch m¨ochte mich bei Prof. Dr. Gerhard Rein recht herzlich fu¨r das Themaund die sehr gute Betreuung bedanken. Mein Dank gilt dabei besondersden vielen fruchtbaren Diskussionen und den Anregungen, die zum Gelingendieser Arbeit beigetragen haben.Mein herzlicher Dank gilt auch meinen Kollegen MSc. Roman Fiˇrt undDipl.-Math. Martin Seehafer fu¨r viele Impulse und fu¨r die unz¨ahligen, nichtnur fachlichen Gespr¨ache.Diese Arbeit wurde finanziell durch die Deutsche Forschungsgemeinschaftunterstu¨tzt.ContentsAbstract 1Preface 3Notation 91 Static shells 111.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111.2 Global existence . . . . . . . . . . . . . . . . . . . . . . . . . . . 151.3 A lower bound on H . . . . . . . . . . . . . . . . . . . . . . . 22C1.4 A scaling lemma . . . . . . . . . . . . . . . . . . . . . . . . . . 271.5 Existence and properties of minimizers . . . . . . . . . . . . . . 291.6 Dynamical stability . . . . . . . . . . . . . . . . . . . . . . . . . 341.7 The case M =0 . . . . . . . . . . . . . . . . . . . . . . . . . . . 36c2 Axially symmetric solutions 432.1 Introduction . . . . . . . . . . . . . . . . . . . . . .

Sujets

• Astronomie

Informations

• Publié par: universitat_bayreuth
• Publié le: 01 janvier 2008
• Nombre de lectures: 14
• Langue: Deutsch

Extrait

Existence and stability of stellardynamic models

beiderUniversit¨atBayreuth zur Erlangung des Grades eines Doktors der Naturwissenschaften (Dr. rer. nat.) vorgelegte Abhandlung von

1. Gutachter: 2. Gutachter:

Tag der Einreichung: Tag des Kolloquiums:

Achim Schulze

geb. 11.07.1979 in Amberg

Danksagung

Ichm¨ochtemichbeiProf.Dr.GerhardReinrechtherzlichfu¨rdasThema und die sehr gute Betreuung bedanken. Mein Dank gilt dabei besonders den vielen fruchtbaren Diskussionen und den Anregungen, die zum Gelingen dieser Arbeit beigetragen haben.

Mein herzlicher Dank gilt auch meinen Kollegen MSc. Roman Fiˇrt und Dipl.-Math.MartinSeehaferf¨urvieleImpulseundfu¨rdieunza¨hligen,nicht nurfachlichenGespra¨che.

Diese Arbeit unterst¨utzt.

wurde

finanziell

durch

die

Deutsche

Forschungsgemeinschaft

Contents

Abstract

Preface

Notation

1

2

Static shells 1.1 Introduction . . . . . . . . . . . . . . . 1.2 Global existence . . . . . . . . . . . . . 1.3 A lower bound onH C . . . . . . . . . 1.4 A scaling lemma . . . . . . . . . . . . 1.5 Existence and properties of minimizers 1.6 Dynamical stability . . . . . . . . . . . 1.7 The caseM c = 0 . . . . . . . . . . . . .

Axially symmetric solutions 2.1 Introduction . . . . . . . . . 2.2 The main result . . . . . . . 2.3Fre´chet-DifferentiabilityofT 2.4∂ ζ T(00) is an isomorphism 2.5 Appendix . . . . . . . . . .

Bibliography

. . . . .

i

. . . . .

. . . . .

. . . . .

. . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

. . . . . . .

. . . . .

1

3

9

11 11 15 22 27 29 34 36

43 43 46 56 70 75

77

ii

CONTENTS

Abstract

IndervorliegendenArbeitwirddieExistenzundStabilit¨atvonstationa¨ren L¨osungendesVlasov–Poisson–Systemsuntersucht.DiesesSystemwirdu.a.in der Stellardynamik zur Beschreibung von Galaxien verwendet. Dabei werden Kollisionenvernachla¨ssigtunddieInteraktionzwischendenTeilchen,d.hder Sterne wird nur durch das von ihnen erzeugte Gravitationspotential bestimmt. Das System besteht aus folgenden Gleichungen:

∂ t f+v ∇ x f− ∇ x U ∇ v f= 0

ΔU= 4πρlimU(tx) = 0 |x|→∞ ρ(tx) = Z f(txv)dv

Dabei bezeichnetf:R×R 6 →R f=f(txv)≥r0f¨ut∈R xv∈R 3 die Phasenraumdichte der Teilchen,U=U(tx) das Gravitationspotential, und ρ=ρ(txcithhcDemuilrea¨)die. ImerstenTeilderArbeitwerdenExistenzundStabilit¨tstation¨arer a Lo¨sungenderForm f 0 (xv) = (E 0 −E) k + (L−L 0 ) l + unter dem Einfluss einer im Ursprung fixierten Punktmasse mittels einer Vari- ationsmethode bewiesen. Dabei ist

E=E(xv)= 1 2 |v| 2 +U(x)− | M x c | 

2 L=|x×v|

undklmit 0< k≤l, sowieE 0 undL 0 sind Konstanten.M c ist die Masse derPunktmasse.DieseLo¨sungenk¨onnenalseinfacheModellef¨urGalaxien verwendet werden, die in ihrem Zentrum ein massives schwarzes Loch besitzen, dasdannann¨aherndalsfixePunktmasseangesehenwerdenkann.Fu¨rdas Vlasov–Poisson System unter der Wirkung eines solchen externen Potentials wirdaucheinglobalerExistenzsatzfu¨rgeeigneteAnfangsdatenbewiesen.

1

2

ABSTRACT

ImzweitenTeilwirddieExistenzaxialsymmetrischerLo¨sungenbehandelt. Eswerdenstation¨areZust¨andekonstruiert,dieineinemlangsamrotierenden SystemdurchDeformationeinergegebenenspha¨rischsymmetrischenLo¨sung entstehen.DazuwirdeinSatz¨uberimpliziteFunktionenaufeinenmodi- fiziert n L¨ tor der Poisson–Gleichung angewendet. Im rotierenden e osungsopera SystemhabendieseZust¨andedieForm f 0 (xv) =ϕ  1 2 v 2 +U(x)− 1 2 ω 2 (x 1 2 +x 2 2 )  

wobeiωdie Rotationsgeschwindigkeit bezeichnet undϕ∈C 1 (R) noch gewissen Zusatzannahmen unterliegt. Der Ausdruck

E 1 2 |v| 2 +U(x)− 2 1 ω 2 (x 1 2 +x 2 2 ) J := wirdauchalsJakobisIntegralbezeichnet.Diehierpra¨sentiertenErgebnisse ¨ sind in Ubereinstimmung mit den numerischen Beobachtungen des Astro- physikers P.O. Vandervoort.

Preface

In the present work, we analyse the time evolution of large particle ensembles, e.g. stars, under the influence of a self-consistent force, which will be the gravitational force here. If relativistic effects are neglected and if the distances between the particles are large enough, an appropiate way to describe such a system are the equations of theN-body-problem:

m j X ¨ j =−γ k6= X N m j m | X X j − X X k k j k j − | 3 

j= 1 N

(0.1)

whereX 1 (s)X 2 (s) X N (s) are threedimensional time-dependent vectors denoting the position of theNpoint masses andm 1 m 2  m N are their masses. ˙ The gravitational constant is denoted byγand we prescribeX j (0)X j (0) j= 1 NwithX j (0)6=X i (0) fori6=j. Many important developments in astro- physics originate from theN-body-problem such as the analysis of motions of the solar system and with the growing computing power of modern processors, the problem is numerically tractable even for largeN. However, for very largeN, such as in galaxies or star clusters, whereN > 10 10 , cf. Figure 1, one is in a very unpromising situation concerning analytical results for this system and it is indeed not reasonable to observe each particle if the focus is for example on stability or the distribution of the particles. In fact, an “averaged” view seems to be more adequate and for many questions it is sufficient to know only how the particles are distributed in phase spaceR 6 = R 3 ×R 3  that purpose, one introduces a function. Forf:R×R 6 →R, which is the phase space densityf=f(txv)≥0, wheret∈Rdenotes time andxv∈R 3 denote position and velocity. Its value is - up to a multiplicative constant - the number of particles per phase space volume at the timet. More precisely, if we take a phase space volume Π⊂R 6 , the number of particlesN(Πt) at time tin the volume Π is N(Π Z Π t) =f(txv)dxdv We have to ascertain the equations which determinefand we start by analysing the motion of a particleXa fixed force field, induced by a givenin

3

4

Figure 1: The spiral galaxy M74, rgtth//:pbbuhiselo.et

potentialU=U(tx). ton’s Law,

PREFACE

The behaviour ofX=X(s) is again governed by New-  V X ˙ = = ˙ V −∇ x U(sX) (0.2) with some initial conditionX(0)V(0). If collissions are neglected, it is natural to require thatfis constant along such particle paths, ifUin (0.2) is the self-consistent potential which is collectively created by all particles. More precisely, we require f(sX(s)V(s)) = const. for all (XV) solving (0.2) and differentiating with respect tosimplies

∂ t f(sZ(s)) +V(s) ∇ x f(sZ(s))− ∇ x U(sX(s)) ∇ v f(sZ(s)) = 0

(0.3)

where we abbreviated (X(s)V(s)) =:Z(s). If we take for given (txv) a suit- ableZ=Z(s) withZ(t) = (X(t)V(t)) = (xv), we arrive at the so-called Vlasov equation, ∂ t f+v ∇ x f− ∇ x U ∇ v f= 0(0.4) and we easily verify thatf=f(txv) solves (0.4), ifffis constant along every solution of (0.2). This statement is just the fact that (0.2) is the characteristic system of the first order PDE (0.4). To reflect the requirement, that the potentialUis created by the particles, we have to couple Eq. (0.4) with Poisson’s equation,

Δ x U(tx) = 4πρ(tx)

limU(tx) = 0 |x|→∞

(0.5)

PREFACE

where we define

ρ(tx) := Z 3 f(txv)dv R sity. Equation (0.5) just says that U(tx) =−  | 1 | ∗ x ρ  (tx) =− Z R 3 | ρ x (t − y y ) | dy

as the spatial den

5

and we also have Z R x−y y)dy ∇ x U(tx) = 3 |x−y| 3 ρ(t which is the averaged form of the right-hand side of equation (0.1), and we assumedρ(t)∈C c (R 3 ). Equations (0.4) and (0.5) form a closed system, the so-called Vlasov–Poisson system. At the beginning of the last century, this system was used by Sir J. Jeans to model stellar clusters and galaxies [14] and to investigate their stability properties. In this context, it is also used by today’s astrophysicists, cf. [4]. Meanwhile, the existence theory for this system is well-understood: Local existence and uniqueness was established by R. Kurth in 1952, cf. [15] and global existence for spherically symmetric solutions was proved by J. Batt in 1977, cf. [3]. Afterwards, global solutions with small initial data were obtained by C. Bardos and P. Degond in 1985, cf. [2]. The final breakthrough was achieved in 1989, when global existence for general classical solutions was independently from each other proved by Pfaffelmoser, cf. [21] and P. L. Lions/B. Perthame, cf. [19]. They proved that for every inital datum ˚ f∈C c 1 (R 6 ), the corresponding solution to the Vlasov–Poisson system exists for all time. We are interested in existence and stability of stationary soltutions to (0.4)– (0.5). For the construction of stationary solutions, we recall thatfis a solution of the Vlasov equation, iff it is constant along characteristics, i.e., the solutions of the charateristic system (0.2). Thus natural candidates for stationary solu- tions are functions of conserved quantities of the characteristic system. One immediate expression for a conserved quantity is the particle energy E(xv):= 1 2 |v| 2 +U(x)(0.6)

so that the ansatz f 0 (xv) :=ϕ(E) (0.7) automatically satisfies the Vlasov equation forϕ∈C 1 . Of course, we still have to solve Poisson’s equation. Plugging (0.7) into (0.5) yields ΔU(x) = 4π Z ϕ  2 1 |v| 2 +U(x)  dv(0.8)