De baan van een baksteen bij een zwart gat

Bereken alle mogelijke banen van een baksteen die ik naar een niet-roterend zwart gat gooi.
Omdat dit baksteenprobleem nogal een uitgebreid verhaal is heb ik het onderverdeeld in vijftien hoofdstukken (als je zo veel mogelijk wiskunde wilt overslaan dan raad ik aan om gelijk naar hoofdstuk 5 te gaan):
  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting

Afleidingen van de vergelijkingen


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Schwarzschild
Schwarzschild

Dit is de Schwarzschild-oplossing:

Vergelijking

Hiermee kom je tot de vergelijking van een geodetische lijn rondom een puntmassa (voor de afleiding zie deze pagina, de vergelijkingen (28) en (30) aldaar), als differentiaalvergelijking en als integraal:
Vergelijking
Vergelijking
Hierin is u de reciproke van de radiële coördinaat r:
Vergelijking
Van deze pagina, waar ik de bovenstaande vergelijkingen (2) heb opgehaald, neem ik ook de vergelijkingen (15) en (21) over:
Vergelijking
Vergelijking
Hierin is J het impulsmoment, een constante, en K is ook een constante. Voor de eigentijd van een waarnemer geldt altijd:
Vergelijking
Dit ga ik invullen in vergelijking (4b):
Vergelijking
Ik stel:
Vergelijking
Omdat K een constante is, is E dat uiteraard ook. Om redenen die later duidelijk zullen worden noem ik E de energie van de baksteen. Daarmee wordt vergelijking (6):
Vergelijking

Onderzoek van de functie


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Ik ga verder werken met vergelijking (2b), de integraal. Wanneer ik die integraal oplos en daarna u vervang door 1/r dan heb ik de oplossing van dit vraagstuk in poolcoördinaten en klaar is Kees. Het venijn zit uiteraard in het oplossen van de integraal, want daarin staat een derdegraads vergelijking, en die zit ook nog verstopt onder een wortelteken, en tot overmaat van ramp staat de wortel in de noemer.

Om te beginnen breng ik de Schwarzschild-straal even in herinnering, de horizon van een zwart gat:
Vergelijking
Met behulp van de vergelijkingen (7) en (9) wordt vergelijking (2b):
Vergelijking
Nu ga ik mij eerst richten op die derdegraads vergelijking in de noemer. Een dergelijke vergelijking ziet er in het algemeen zo uit:
Vergelijking
Het is duidelijk in vergelijking (10) dat a > 0, dus die derdegraads vergelijking in de noemer ‘begint’ ergens linksonder (in het derde kwadrant) en ‘eindigt’ ergens rechtsboven (in het eerste kwadrant).
Grafiek
a < 0
Grafiek
a > 0
Tabel 1
Van een derdegraads vergelijking kan ik de discriminant uitrekenen:
Vergelijking
De discriminant vertelt mij hoeveel nulpunten er zijn (met nulpunten bedoel ik op deze pagina de reële nulpunten, oftewel de punten waar de functie de horizontale as snijdt of raakt, dus complexe nulpunten laat ik buiten beschouwing):
Grafiek
D < 0
Grafiek
D = 0
Grafiek
D > 0
Tabel 2
Laat ik daar dan maar eens mee beginnen, ik ga de discriminant uitrekenen:
Vergelijking
Ik heb nu een tweedegraads vergelijking in E2, een parabool, en daar maak ik even een plaatje van.
Grafiek
Figuur 1
De grafiek van D (E2) voor J = 1 (de rode lijn), J = √3 (de oranje lijn),
J = 3 (de groene lijn), J = 4 (de paarse lijn),
J = 5 (de blauwe lijn), J = 6 (de grijze lijn),
J = 7 (de bruine lijn) en J = 8 (de gele lijn),
Rs = 1
Zoals figuur 1 laat zien kan de discriminant negatief zijn, nul of positief. Oftewel, hij kan alle waarden aannemen. Het is natuurlijk wel interessant om te onderzoeken wanneer die discriminant dan gelijk aan nul wordt. Daarvoor gebruik ik uiteraard de abc-formule:
Vergelijking
We verkeren in de gelukkige omstandigheid dat het deel onder de wortel te ontbinden is in gelijke factoren:
Vergelijking
En zo vind ik de beide nulpunten:
Vergelijking
Hier maak ik ook een plaatje van.
Grafiek
Figuur 2
De grafiek van E2 (J) voor de plusoplossing (de rode lijn)
en de minoplossing (de groene lijn),
D = 0, Rs = 1
Om redenen die later duidelijk zullen worden maak ik de voorgaande grafiek nog een keer, maar dan met E in plaats van E2 langs de verticale as.
Vergelijking
Grafiek
Figuur 3
De grafiek van E (J) voor de plusoplossing (de rode lijn)
en de minoplossing (de groene lijn),
D = 0, Rs = 1
Ik kijk ook even wat er gebeurt voor grote waarden van J:
Vergelijking
De linkerterm is een parabool en komt overeen met de rode lijn in figuur 2, en de rechterterm komt overeen met de groene lijn in figuur 2.

Ik kan in vergelijking (16a) ook aflezen waar D een maximum bereikt, de top van de parabool:
Vergelijking
De waarde van D is dan:
Vergelijking
Er treedt een bijzonder geval op wanneer de top precies aan de horizontale as raakt. Deze top vind ik door vergelijking (19) nul te stellen:
Vergelijking
En inderdaad vallen in dit geval de nulpunten volgens vergelijking (16a) samen. Dan is E (vergelijking (18)):
Vergelijking
Als al dit rekenwerk klopt dan is voor deze waarden van J en E de discriminant D (vergelijking (13)) gelijk aan nul:
Vergelijking
Het punt (J, E2) = (√3, 8/9) ≈ (1.732, 0.889) in figuur 2 is inderdaad het punt waar de rode lijn en de groene lijn samenkomen. Links daarvan is de discriminant altijd negatief.

Ik ga vergelijking (13) anders opschrijven:
Vergelijking
Ik heb nu een tweedegraads vergelijking in J2 (het deel tussen de haakjes), een parabool, en daar maak ik twee plaatjes van, voor E < 1 en E > 1.
Grafiek
Figuur 4
De grafiek van D (J2) (alleen het deel tussen de haakjes)
voor E = 0.2 (de rode lijn), E = 0.3 (de oranje lijn),
E = 0.4 (de groene lijn), E = 0.5 (de paarse lijn),
E = 0.6 (de blauwe lijn), E = 0.7 (de grijze lijn),
E = 0.8 (de bruine lijn) en E = 0.9 (de gele lijn),
Rs = 1
Grafiek
Figuur 5
De grafiek van D (J2) (alleen het deel tussen de haakjes)
voor E = 1.1 (de rode lijn), E = 1.2 (de oranje lijn),
E = 1.3 (de groene lijn), E = 1.4 (de paarse lijn),
E = 1.5 (de blauwe lijn), E = 1.6 (de grijze lijn),
E = 1.7 (de bruine lijn) en E = 1.8 (de gele lijn),
Rs = 1
Zoals de figuren 4 en 5 laten zien kan de discriminant negatief zijn, nul of positief. Oftewel, hij kan alle waarden aannemen, maar voor E < 1 is de discriminant altijd negatief. Het is natuurlijk wederom interessant om te onderzoeken wanneer die discriminant dan gelijk aan nul wordt. Dus nogmaals de abc-formule in actie om de beide nulpunten te vinden:
Vergelijking
We verkeren weer in de gelukkige omstandigheid dat het deel onder de wortel te ontbinden is in gelijke factoren:
Vergelijking
En zo vind ik de beide nulpunten:
Vergelijking
Hier maak ik ook een plaatje van.
Grafiek
Figuur 6
De grafiek van J2 (E) voor de plusoplossing (de rode lijn)
en de minoplossing (de groene lijn),
D = 0, Rs = 1
Ook in dit geval maak ik de voorgaande grafiek nog een keer, maar dan met J in plaats van J2 langs de verticale as.
Vergelijking
Grafiek
Figuur 7
De grafiek van J (E) voor de plusoplossing (de rode lijn)
en de minoplossing (de groene lijn),
D = 0, Rs = 1
Ik kijk wat er gebeurt voor grote waarden van E:
Vergelijking
De linkerterm is een parabool en komt overeen met de rode lijn in figuur 6, en de rechterterm komt overeen met de groene lijn in figuur 6.

Ik kan in vergelijking (26a) ook aflezen waar D een maximum bereikt, de top van de parabool:
Vergelijking
De waarde van D is dan:
Vergelijking
Er treedt een bijzonder geval op wanneer de top precies aan de horizontale as raakt. Deze top vind ik door vergelijking (29) nul te stellen:
Vergelijking
En inderdaad vallen in dit geval de nulpunten volgens vergelijking (26a) samen. Dan is J (vergelijking (28)):
Vergelijking
Vergelijking
Vergelijking
Het moge duidelijk zijn dat het impulsmoment niet imaginair kan zijn noch oneindig en dat alleen de derde oplossing voldoet. Als al dit rekenwerk klopt dan is voor deze waarden van J en E de discriminant D (vergelijking (23)) gelijk aan nul:
Vergelijking
Het punt (E, J2) = (2/3√2, 3) ≈ (0.943, 3) in figuur 6 is inderdaad het punt waar de rode lijn en de groene lijn samenkomen. Links daarvan is de discriminant altijd negatief. Merk op dat de figuren 3 (E als functie van J) en 7 (J als functie van E) elkaars inverse zijn (door de ene grafiek te spiegelen in de lijn E = J onstaat de andere grafiek).
Grafiek
Figuur 8
De grafiek van E (J) = figuur 3 (links) en de grafiek van J (E) = figuur 7 (rechts),
D = 0, Rs = 1
We weten nu de criteria wanneer de derdegraads vergelijking één (D < 0), twee (D = 0) of drie (D > 0) nulpunten heeft.
Grafiek
Figuur 9
De grafiek van E (J),
D = 0, Rs = 1
Grafiek
Figuur 10
De grafiek van J (E),
D = 0, Rs = 1
Ik heb de discriminant eerst als functie van E2 geschreven (vergelijking (13)) en daarna als functie van J2 (vergelijking (23)):
Vergelijking
Vergelijking
Daarna heb ik in beide gevallen uitgebreid onderzocht waar zich de extrema bevinden, en in het tweede geval deed ik dat alleen voor het deel tussen de haakjes. Voor de volledigheid maak ik ook nog een plaatje van de gehele vergelijking (23) (en niet alleen het deel tussen de haakjes).
Grafiek
Figuur 11
De grafiek van D (J2)
voor E = 0 (de rode lijn), E = 0.25 (de oranje lijn),
E = 0.5 (de groene lijn), E = 0.75 (de paarse lijn),
E = 1 (de blauwe lijn), E = 1.25 (de grijze lijn),
E = 1.5 (de bruine lijn) en E = 1.75 (de gele lijn),
Rs = 1
Tenslotte wil ik proberen D als functie van E en J in een 3D-plot te vangen.
Grafiek
Figuur 12
De grafiek van D (E, J),
Rs = 1

De effectieve potentiaal


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dan is er nog een ander onderwerp dat besproken moet worden en daarvoor ga ik uit van de snelheid van de baksteen. Die snelheid heb ik afgeleid op deze pagina, dat is de snelheid van een baksteen die een zwart gat nadert (vergelijking (6) op die pagina) bezien vanuit een waarnemer ergens ver verwijderd van het zwarte gat:
Vergelijking
Ik ga links en rechts kwadrateren en nog wat knutselen:
Vergelijking
Dit combineer ik met de vergelijkingen (5), (7), (8) en (9):
Vergelijking
De rechterterm noemen we de effectieve potentiaal (in het kwadraat) V2:
Vergelijking
Oftewel:
Vergelijking
Hier maak ik ook een plaatje van.
Grafiek
Figuur 13
De grafiek van V (r) voor J = 0.0 (de rode lijn),
J = 0.5 (de lichtgroene lijn), J = 1.0 (de oranje lijn),
J = 1.5 (de paarse lijn), J = √3 (de blauwe lijn),
J = 2.0 (de grijze lijn), J = 2.5 (de bruine lijn),
J = 3.0 (de lichtbruine lijn), J = 3.5 (de donkergroene lijn),
J = 4.0 (de gele lijn) en J = 4.5 (de lichtblauwe lijn),
Rs = 1
Ook hier is het weer interessant om de extrema op te zoeken en daarom ga ik de afgeleide naar r bepalen:
Vergelijking
Met behulp van de abc-formule (toegepast op de teller) vind ik de extrema:
Vergelijking
Deze punten geef ik aan in figuur 13.
Grafiek
Figuur 14
De grafiek van V (r) voor J = 0.0 (de rode lijn),
J = 0.5 (de lichtgroene lijn), J = 1.0 (de oranje lijn),
J = 1.5 (de paarse lijn), J = √3 (de blauwe lijn),
J = 2.0 (de grijze lijn), J = 2.5 (de bruine lijn),
J = 3.0 (de lichtbruine lijn), J = 3.5 (de donkergroene lijn),
J = 4.0 (de gele lijn) en J = 4.5 (de lichtblauwe lijn),
de zwarte punten zijn de extrema volgens vergelijking (39),
Rs = 1
Merk op dat de onderste vier lijnen (rood, lichtgroen, oranje en paars) helemaal geen nulpunten hebben in de afgeleide van de potentiaal en dat de potentiaalfunctie voor die waarden van J (J < √3) dus nergens horizontaal loopt. Merk ook op dat de linker zes punten allemaal op een top liggen en dat dat dus geen stabiele punten zijn, of in andere woorden: dat zijn geen stabiele banen om om een zwart gat te draaien. Deze zes punten zijn trouwens allemaal minoplossingen van vergelijking (39), de plusoplossingen liggen rechts buiten beeld. Alleen de blauwe lijn heeft een punt (op drie Schwarzschild-stralen van het centrum van het zwarte gat) die op de rand van stabiliteit verkeert (de plusoplossing en de minoplossing vallen samen). Conclusie: binnen twee Schwarzschild-stralen van de horizon bevinden zich geen stabiele banen en verdwijnt uiteindelijk alles in het zwarte gat (voor meer details zie verderop op deze pagina en ook deze pagina). Ik zal figuur 14 nogmaals laten zien, maar dan met een andere horizontale schaalverdeling zodat de plusoplossingen ook zichtbaar worden.
Grafiek
Figuur 15
De grafiek van V (r) voor J = 0.0 (de rode lijn),
J = 0.5 (de lichtgroene lijn), J = 1.0 (de oranje lijn),
J = 1.5 (de paarse lijn), J = √3 (de blauwe lijn),
J = 2.0 (de grijze lijn), J = 2.5 (de bruine lijn),
J = 3.0 (de lichtbruine lijn), J = 3.5 (de donkergroene lijn),
J = 4.0 (de gele lijn) en J = 4.5 (de lichtblauwe lijn),
de zwarte punten zijn de nulpunten volgens vergelijking (39),
Rs = 1
De combinatie van de vergelijkingen (35) en (36) levert aldus op:
Vergelijking

Nulpunten in de noemer


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dan kunnen we ons nu vol op het oplossen van de derdegraads vergelijking storten. Die ziet er zo uit:
Vergelijking
Ik ga vier hulpvariabelen in het leven roepen:
Vergelijking
Vergelijking
Vergelijking
Vergelijking
Dit stelt mij in staat in alle situaties de nulpunten uit te rekenen (voor alle details zie de pagina over de derdegraads vergelijking).
Grafiek
D < 0
Eén nulpunt:
Vergelijking
Grafiek
D = 0
Twee nulpunten, een snijpunt en een raakpunt
(het raakpunt is feitelijk twee samenvallende nulpunten):
Vergelijking
Vergelijking
Indien q'' > 0 dan ligt het raakpunt rechts van het snijpunt
en indien q'' < 0 dan ligt het raakpunt links van het snijpunt.
Grafiek
D > 0
Drie nulpunten, van links naar rechts:
Vergelijking
Vergelijking
Vergelijking
Tabel 3
Een logische volgende stap is om die vier hulpvariabelen uit te rekenen:
Vergelijking
Vergelijking
Vergelijking
Vergelijking

Radiële inval


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Vervolgens begin ik met het allersimpelste geval: radiële inval. Dan is het impulsmoment J = 0 en wordt de integraal:
Vergelijking
Zoals verwacht is in het geval van radiële inval de hoek φ een constante, en of ik de baksteen gooi of laat vallen maakt niets uit voor (de vorm van) de baan die de baksteen volgt. In een plaatje ziet dat er zo uit.
Grafiek
Figuur 16
De grafiek van φ (x, y)
voor φ0 = 0,
J = 0,
Rs = 1
Grafiek
Figuur 17
De grafiek van φ (x, y)
voor φ0 = π/6,
J = 0,
Rs = 1
In het geval dat J = 0 is de waarde van E niet relevant en eindigt de baksteen onherroepelijk via de kortst mogelijke weg in het zwarte gat.
Condities Baan van de baksteen
J = 0 De baksteen valt radieel in het zwarte gat
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
Tabel 4

Invallende banen vanaf ‘ver weg’ (D < 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dit was een eenvoudig inkoppertje, maar in alle andere situaties zal toch echt de integraal helemaal uitgewerkt moeten worden. Ik begin met de situatie waarbij D < 0, dus er is één nulpunt en dat nulpunt kan ik uitrekenen volgens vergelijking (43):
Vergelijking
De integraal zoek ik op in de tabel met integralen:
Vergelijking
Taylor
Taylor

Deze oplossing is bereikt door de integraal om te schrijven naar een elliptische integraal van de eerste soort. Je ziet in de oplossing dat de oorspronkelijke variabelen onherkenbaar veranderd zijn doordat er onderweg talrijke transformaties hebben plaatsgevonden, maar het ergste is dat de oplossing in dit geval (ruimschoots) onvoldoende praktisch nut heeft omdat de oplossing tot stand gekomen is middels een Taylor-reeks. Willen we daarmee verder werken dan zullen we vele duizenden termen mee moeten nemen en dat wordt zelfs voor de computer (mijn computer) een onhaalbare rekenexercitie.

Daarom ga ik het heel anders aanpakken en ga ik voor het uitwerken van de integraal gebruik maken van de trapeziummethode:
Vergelijking
Bij ieder volgend plaatje vermeld ik daarom het aantal intervallen dat ik gebruikt heb in de berekening, en tevens is vanaf nu φ0 = 0.

In een simpel geval dat het inkomende object nauwelijks impulsmoment heeft dan ziet het resultaat er als volgt uit, vergelijkbaar met de radiële inval.
Grafiek
Figuur 18
De grafiek van φ (x, y)
voor J = 0.1,
D = −7.044 ∙ 106,
E = 2.0117769801886,
n1 = −2.761,
Rs = 1,
999.999 intervallen
Ik ga het impulsmoment opvoeren, in het volgende plaatje is J = 2.
Grafiek
Figuur 19
De grafiek van φ (x, y)
voor J = 2,
D = −16.051,
E = 2.0117769801886,
n1 = −0.614,
Rs = 1,
999.999 intervallen
Langzaam maar zeker wordt er afbuiging zichtbaar, in het volgende plaatje is J = 3.
Grafiek
Figuur 20
De grafiek van φ (x, y)
voor J = 3,
D = −2.411,
E = 2.0117769801886,
n1 = −0.447,
Rs = 1,
999.999 intervallen
Door naar J = 4.
Grafiek
Figuur 21
De grafiek van φ (x, y)
voor J = 4,
D = −0.429,
E = 2.0117769801886,
n1 = −0.353,
Rs = 1,
999.999 intervallen
En ook nog J = 4.9
Grafiek
Figuur 22
De grafiek van φ (x, y)
voor J = 4.9,
D = −0.021,
E = 2.0117769801886,
n1 = −0.297,
Rs = 1,
9.999.999 intervallen
Voor E = 2.01177698018854097538... en J = 5 is D precies gelijk aan nul en krijgt g (u) een raakpunt (of preciezer gezegd: twee samenvallende nulpunten):
Vergelijking
De onderstaande grafiek laat zien hoe dat raakpunt zich vormt.
Grafiek
Figuur 23
De grafiek van g (u) voor J = 3 (de rode lijn)
in stappen van 0.2 oplopend tot J = 5 (de lichtblauwe lijn),
E = 2.01177698018854, Rs = 1
De vorming van dat raakpunt zorgt voor een piek in de integrand f (u):
Vergelijking
Ook dat laat ik zien in een grafiek.
Grafiek
Figuur 24
De grafiek van f (u) voor J = 3 (de rode lijn)
in stappen van 0.2 oplopend tot J = 5 (de lichtblauwe lijn),
E = 2.01177698018854, Rs = 1
Die piek in de integrand maakt het integreren lastig (dat is heel zacht uitgedrukt), maar zorgt wel voor de mooiste en interessantste plaatjes. Daarom ga ik proberen zo dicht mogelijk bij D = 0 (waar het raakpunt ontstaat) te komen.
Grafiek
Figuur 25
De grafiek van φ (x, y)
voor J = 5,
D = −3.125 ∙ 10−14,
E = 2.0117769801886,
n1 = −0.292,
Rs = 1,
5 rondjes om het zwarte gat,
40.959.718 intervallen
Het is wel mooi om even in te zoomen op de omgeving van het zwarte gat.
Grafiek
Figuur 26
De grafiek van φ (x, y)
voor J = 5,
D = −3.125 ∙ 10−14,
E = 2.0117769801886,
n1 = −0.292,
Rs = 1,
5 rondjes om het zwarte gat,
40.959.718 intervallen
Voor alle duidelijkheid: de baksteen cirkelt in dit geval vijfmaal om het zwarte gat heen en zelfs bij een inzoomcompositie wordt dit niet duidelijk zichtbaar (het blijkt uit de datafile).
Grafiek
Figuur 27
De grafiek van φ (x, y)
voor J = 5,
D = −3.125 ∙ 10−14,
E = 2.0117769801886,
n1 = −0.292,
Rs = 1,
5 rondjes om het zwarte gat,
40.959.718 intervallen
Ik ga nog een stap verder door nog dichter tegen de lijn D = 0 aan te kruipen.
Grafiek
Figuur 28
De grafiek van φ (x, y)
voor J = 5,
D = −4.428 ∙ 10−21,
E = 2.01177698018854097539,
n1 = −0.292,
Rs = 1,
7.5 rondjes om het zwarte gat,
192.772.784 intervallen
Wat ik in de voorgaande stappen heb gedaan is de energie constant gehouden en het impulsmoment steeds verder opgevoerd tot aan de rand van D = 0, zie het onderstaande overzicht.
Grafiek
Figuur 29
De grafiek van E (J), D = 0, Rs = 1,
met daarin de punten (de groene bolletjes) van de voorgaande figuren

Invallende banen vanaf ‘dichtbij’ (D < 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dan is er nog een ander punt van aandacht en daarvoor pak ik nogmaals het polynoom g (u) erbij:
Vergelijking
De laatste term bepaalt de waarde van g (0), het snijpunt met de verticale as. Voor E > 1 is g (0) > 0 en voor E < 1 is g (0) < 0. Echter, wanneer g < 0 dan bestaat de integrand f (u) niet:
Vergelijking
Grafiek
Figuur 30
De grafiek van g (u) voor E = 0 (de rode lijn)
in stappen van 0.2 oplopend tot E = 2 (de lichtblauwe lijn),
de grijze lijn is E = 1,
J = 1, Rs = 1
Grafiek
Figuur 31
De grafiek van f (u) voor E = 0 (de rode lijn)
in stappen van 0.2 oplopend tot E = 2 (de lichtblauwe lijn),
de grijze lijn is E = 1,
J = 1, Rs = 1
Voor E < 1 bestaat de integrand dus niet vanaf u = 0 tot aan het nulpunt van de derdegraads vergelijking. En omdat r = 1/u betekent dit dat de integrand niet bestaat voor grote waarden van r. De waarden van r < Rs zijn niet interessant, omdat dat binnen de horizon van het zwarte gat is. Dit komt overeen met u > 1 en dan blijft er het volgende over van figuur 31.
Grafiek
Figuur 32
De grafiek van f (u) voor E = 0 (de rode lijn)
in stappen van 0.2 oplopend tot E = 2 (de lichtblauwe lijn),
de grijze lijn is E = 1,
J = 1, Rs = 1
Ik zal dit even schematisch laten zien. Voor D < 0 en E > 1 bestaat de integrand voor het gehele traject van r = ∞ tot r = Rs.
Grafiek
Figuur 33
Het traject waar de integrand f (u) bestaat (de oranje balk),
D < 0, E > 1
Voor E < 1 bestaat de integrand vanaf het nulpunt van de derdegraads vergelijking.
Grafiek
Figuur 34
Het traject waar de integrand f (u) bestaat (de oranje balk),
D < 0, E < 1
Hoe ziet de baan van de baksteen er dan uit? Ik neem daarvoor de blauwe lijn uit figuur 32.
Grafiek
Figuur 35
De grafiek van φ (x, y)
voor E = 0.8,
D = −1.459,
J = 1,
n1 = 0.480,
1/n1 = 2.084,
Rs = 1,
9.999.999 intervallen
Ik ga E in stapjes verhogen.
Grafiek
Figuur 36
De grafiek van φ (x, y)
voor E = 0.9,
D = −1.315,
J = 1,
n1 = 0.231,
1/n1 = 4.328,
Rs = 1,
9.999.999 intervallen
Grafiek
Figuur 37
De grafiek van φ (x, y)
voor E = 0.95,
D = −1.892,
J = 1,
n1 = 0.108,
1/n1 = 9.269,
Rs = 1,
9.999.999 intervallen
Grafiek
Figuur 38
De grafiek van φ (x, y)
voor E = 0.99,
D = −2.732,
J = 1,
n1 = 0.020,
1/n1 = 49.252,
Rs = 1,
9.999.999 intervallen
Grafiek
Figuur 39
De grafiek van φ (x, y)
voor E = 0.999,
D = −2.972,
J = 1,
n1 = 0.002,
1/n1 = 499.250,
Rs = 1,
9.999.999 intervallen
Wellicht denk je “voer het impulsmoment dan even flink op”, maar zo simpel ligt dat niet zoals uit onderstaande grafiek blijkt (het is dezelfde grafiek als figuur 32, maar dan met J = 5 in plaats van J = 1).
Grafiek
Figuur 40
De grafiek van f (u) voor E = 0.1 (de oranje lijn)
in stappen van 0.1 oplopend tot E = 1 (de lichtblauwe lijn),
J = 5, Rs = 1
Ik zal even inzoomen op het rechterdeel.
Grafiek
Figuur 41
De grafiek van f (u) voor E = 0.1 (de oranje lijn)
in stappen van 0.1 oplopend tot E = 1 (de lichtblauwe lijn),
J = 5, Rs = 1
De baan van de baksteen ziet er dan bijvoorbeeld zo uit.
Grafiek
Figuur 42
De grafiek van φ (x, y)
voor E = 0.9,
D = −0.025,
J = 5,
n1 = 0.967,
1/n1 = 1.034,
Rs = 1,
9.999.999 intervallen
Voor E < 1 ontstaan er dus banen die niet ergens ver weg (ver van het zwarte gat) beginnen, maar banen die dichtbij het zwarte gat ontstaan (omdat ik vanaf die plek de baksteen gooi). Toch zit daar ook nog weer een kanttekening aan en daarvoor grijp ik helemaal terug op figuur 9 en het gaat mij om het deel dat ik hieronder rood gemaakt heb.
Grafiek
Figuur 43
De grafiek van E (J), D = 0, Rs = 1
Voor dat rode vlakje geldt dat E < 1, maar D > 0 en dus heeft de functie g (u) daar drie nulpunten.
Grafiek
Figuur 44
De grafiek van g (u) voor J = 1.6 de rode lijn)
in stappen van 0.06 oplopend tot J = 2.2 (de lichtblauwe lijn),
E = 0.999, Rs = 1
En daar waar g (u) < 0 bestaat de integrand f (u) niet. In dit geval wordt de integrand in tweeën gesplitst, een linkerdeel en een rechterdeel.
Grafiek
Figuur 45
De grafiek van f (u) voor J = 1.6 de rode lijn)
in stappen van 0.06 oplopend tot J = 2.2 (de lichtblauwe lijn),
E = 0.999, Rs = 1
Deze situatie parkeer ik even (daar kom ik later op terug, zie hoofdstuk 13). Op dit moment kunnen we dus ook de conclusie trekken dat voor D < 0 de baksteen gegarandeerd in het zwarte gat eindigt.
Condities Baan van de baksteen
J = 0 De baksteen valt radieel in het zwarte gat
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
D < 0 E ≥ 1 De baksteen eindigt in het zwarte gat vanaf ‘ver weg’
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n1 < u < 1/Rs oftewel Rs < r < 1/n1)
Grafiek
Grafiek
Grafiek
Tabel 5

Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
De volgende categorie is die waarbij de discriminant gelijk is aan nul. In dat geval hebben we twee nulpunten: een snijpunt en een raakpunt. We kunnen dan twee gevallen onderscheiden.
Grafiek
q'' < 0
Grafiek
q'' > 0
Tabel 6
De waarde van q'' had ik hiervoor al uitgerekend:
Vergelijking
En ik heb uitdrukkingen gevonden voor E2 en J2 in het geval dat D = 0:
Vergelijking
Vergelijking
Door de vergelijkingen (16a) en (26a) bij toerbeurt in te vullen in vergelijking (46b) vind ik uitdrukkingen voor q'' waarbij q'' enkel en alleen afhangt van E respectievelijk J. Eerst vul ik vergelijking (16a) in:
Vergelijking
Hier maak ik weer een plaatje van.
Grafiek
Figuur 46
De grafiek van q'' (J) voor de plusoplossing (de rode lijn)
en de minoplossing (de groene lijn),
D = 0, Rs = 1
Vervolgens vul ik vergelijking (26a) in:
Vergelijking
Ik ga de noemer even apart onder handen nemen:
Vergelijking
Dan komt vergelijking (54) er zo uit te zien:
Vergelijking
Nu ga ik de teller even apart onder handen nemen:
Vergelijking
Hiermee wordt vergelijking (56):
Vergelijking
Ook dit vraagt natuurlijk om een plaatje.
Grafiek
Figuur 47
De grafiek van q'' (E) voor de minoplossing (de rode lijn)
en de plusoplossing (de groene lijn),
D = 0, Rs = 1
In de figuren 9 en 10 kan ik nu informatie toevoegen betreffende q''.
Grafiek
Figuur 48
De grafiek van E (J),
D = 0, Rs = 1
Grafiek
Figuur 49
De grafiek van J (E),
D = 0, Rs = 1
Het moge duidelijk zijn dat het een stuk prettiger is om verder te werken met vergelijking (53) dan met vergelijking (58). Vergelijking (53) ga ik dan ook gebruiken om de nulpunten uit te rekenen (volgens de vergelijkingen (44)):
Vergelijking
Vergelijking
Deze nulpunten vragen natuurlijk ook om een plaatje.
Grafiek
Figuur 50
De grafiek van n1 (J) (de minoplossing, de rode lijn),
n1 (J) (de plusoplossing, de oranje lijn),
n2/3 (J) (de plusoplossing, de groene lijn)
en n2/3 (J) (de minoplossing, de paarse lijn),
Rs = 1
Ik begin met de situatie dat q'' < 0, dus het raakpunt ligt links van het snijpunt, oftewel, n2/3 < n1. n2/3 volgt dan de paarse lijn in figuur 50 en n1 de oranje lijn. Voor q'' < 0 geldt sowieso dat E < 1 zoals de figuren 48 en 49 laten zien. Schematisch ziet dat er als volgt uit, de integrand bestaat voor het traject van r = 1/n1 tot r = Rs.
Grafiek
Figuur 51
Het traject waar de integrand f (u) bestaat (de oranje balk),
D = 0, E < 1, q'' < 0
Dat we een raakpunt hebben (twee samenvallende nulpunten) is goed nieuws, want daardoor is de integraal exact op te lossen. Dus de oplossing van de integraal zoek ik simpelweg op in de tabel met integralen:
Vergelijking
Voor de variabele u geldt:
Vergelijking
Hiermee wordt vergelijking (60):
Vergelijking
Vervolgens reken ik het verschil van de nulpunten uit (en in figuur 50 kan ik aflezen dat ik bij n1 de plusoplossing moet gebruiken en bij n2/3 de minoplossing):
Vergelijking
Vergelijking
Vergelijking
Hiermee wordt vergelijking (61):
Vergelijking
Het is weer de hoogste tijd voor plaatjes. Ik begin met J = 3 en ik ga het impulsmoment in stappen verlagen.
Grafiek
Figuur 52
De grafiek van φ (x, y) voor
J = 3,
E = 0.9851121580978086,
D = 0,
n1 = 0.878,
1/n1 = 1.139,
q'' = −2.016 ∙ 10−2,
Rs = 1,
1.000.002 datapunten
Vervolgens J = 2 en lager.
Grafiek
Figuur 53
De grafiek van φ (x, y) voor
J = 2,
E = 0.9622504486493763,
D = 0,
n1 = 0.667,
1/n1 = 1.500,
q'' = −4.630 ∙ 10−3,
Rs = 1,
1.000.002 datapunten
Grafiek
Figuur 54
De grafiek van φ (x, y)
voor J = 1.9,
E = 0.9568176372873567,
D = 0,
n1 = 0.607,
1/n1 = 1.646,
q'' = −2.573 ∙ 10−3,
Rs = 1,
1.000.001 datapunten
Grafiek
Figuur 55
De grafiek van φ (x, y)
voor J = 1.8,
E = 0.9496463385962782,
D = 0,
n1 = 0.515,
1/n1 = 1.943,
q'' = −7.467 ∙ 10−4,
Rs = 1,
1.000.002 datapunten
Grafiek
Figuur 56
De grafiek van φ (x, y)
voor J = 1.7321,
E = 0.9428156854432492,
D = 0,
n1 = 0.338,
1/n1 = 2.955,
q'' = −1.585 ∙ 10−8,
Rs = 1,
5 rondjes om het zwarte gat,
249.999.998 datapunten
Wat ik in de voorgaande stappen heb gedaan is het impulsmoment steeds verder verlaagd, waarbij ik E zo gekozen heb dat D = 0 en q'' < 0. Zie het onderstaande plaatje.
Grafiek
Figuur 57
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D = 0, Rs = 1
Zoals je ziet eindigt ook in al deze gevallen de baksteen in het zwarte gat, maar naarmate ik het omslagpunt van q'' < 0 naar q'' > 0 (het punt q'' = 0) nader duurt het wel steeds langer voordat de baksteen geconfronteerd wordt met het onvermijdelijke. In figuur 56 draait de baksteen vijf rondjes om het zwarte gat voordat hij erin verdwijnt. Door nog dichter tegen het punt q'' = 0 aan te kruipen zal de baksteen nog meer rondjes om het zwarte gat draaien, in het onderstaande plaatje doet ie dat om precies te zijn maar liefst 173 keer!
Grafiek
Figuur 58
De grafiek van φ (x, y)
voor J = 1.7320508076,
E = 0.9428090415862986,
D = 0,
n1 = 0.333337,
1/n1 = 2.999964,
q'' = −7.979 ∙ 10−18,
Rs = 1,
173.5 rondjes om het zwarte gat,
400.000.000 datapunten

Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dan ga ik nu de situatie onderzoeken waarbij q'' > 0, dus het raakpunt ligt nu rechts van het snijpunt, oftewel, n2/3 > n1. n2/3 volgt dan de groene lijn in figuur 50 en n1 de rode lijn. Wat natuurlijk direct opvalt is dat de rode lijn onder de horizontale as duikt en n1 dus negatief kan worden (daar waar E ≥ 1 overgaat naar E < 1). Daarnaast is de integrand f (u) bij het raakpunt niet gedefinieerd (omdat de noemer daar nul is) en wordt dus in twee delen gesplitst. Het linkerdeel (de groene lijn in figuur 60) ligt tussen het snijpunt en (links van) het raakpunt, het rechterdeel (de blauwe lijn in figuur 60) ligt rechts van het raakpunt.
Grafiek
Figuur 59
De grafiek van g (u) met een snijpunt n1 = 0.5
en een raakpunt n2/3 = 3
Grafiek
Figuur 60
De grafiek van f (u) met een snijpunt n1 = 0.5
en een raakpunt n2/3 = 3
Schematisch gezien voorspelt dat drie categorieën banen.
Grafiek
Figuur 61
Het traject waar de integrand f (u) bestaat (de oranje balk),
D = 0, E ≥ 1, q'' > 0
Grafiek
Figuur 62
Het traject waar de integrand f (u) bestaat (de oranje balk),
D = 0, E ≥ 1, q'' > 0
Grafiek
Figuur 63
Het traject waar de integrand f (u) bestaat (de oranje balk),
D = 0, E < 1, q'' > 0
Dat we een raakpunt hebben (twee samenvallende nulpunten) is nog steeds goed nieuws, want daardoor is de integraal exact op te lossen. Dus de oplossing van de integraal zoek ik weer op in de tabel met integralen (en ik werk het eerst uit voor de situatie links van het raakpunt, dus u < n2/3, oftewel het groene deel in figuur 60, oftewel de oranje balk in de figuren 62 en 63):
Vergelijking
Voor de variabele u geldt:
Vergelijking
Hiermee wordt vergelijking (64):
Vergelijking
Vervolgens reken ik het verschil van de nulpunten uit (en in figuur 50 kan ik aflezen dat ik bij n1 de minoplossing moet gebruiken en bij n2/3 de plusoplossing):
Vergelijking
Vergelijking
Vergelijking
Hiermee wordt vergelijking (65):
Vergelijking
En zo wordt het weer de hoogste tijd voor plaatjes. Bij de volgende reeks grafieken vermeld ik het aantal rondjes dat de baksteen om het zwarte gat draait, maar dat is puur voor de vorm. De baksteen komt in een stabiele baan, dus indien ik de simulatie langer laat doorlopen dan neemt het aantal rondjes evenredig toe. Ik begin met J = 4 en ik ga het impulsmoment in stappen verlagen.
Grafiek
Figuur 64
De grafiek van φ (x, y)
voor J = 4,
E = 1.6491982708836967,
D = 0,
n1 = −0.268,
n2/3 = 0.634,
1/n2/3 = 1.578,
q'' = 2.713 ∙ 10−2,
Rs = 1,
3 rondjes om het zwarte gat,
27.365.598 datapunten
Vervolgens J = 2.5 en lager.
Grafiek
Figuur 65
De grafiek van φ (x, y)
voor J = 2.5,
E = 1.141591103064672,
D = 0,
n1 = −0.147,
n2/3 = 0.574,
1/n2/3 = 1.743,
q'' = 1.389 ∙ 10−2,
Rs = 1,
3.5 rondjes om het zwarte gat,
29.187.870 datapunten
Grafiek
Figuur 66
De grafiek van φ (x, y)
voor J = 2.1,
E = 1.025861056739398,
D = 0,
n1 = −0.044,
n2/3 = 0.522,
1/n2/3 = 1.916,
q'' = 6.696 ∙ 10−3,
Rs = 1,
4 rondjes om het zwarte gat,
31.099.521 datapunten
Merk op dat er nu een overgang plaatsvindt van E ≥ 1 naar E < 1 en dat de integrand daarom niet meer verloopt volgens figuur 62, maar volgens figuur 63.
Grafiek
Figuur 67
De grafiek van φ (x, y)
voor J = 1.9,
E = 0.9760369650397409,
D = 0,
n1 = 0.059,
1/n1 = 16.866,
n2/3 = 0.470,
1/n2/3 = 2.126,
q'' = 2.573 ∙ 10−3,
Rs = 1,
5 rondjes om het zwarte gat,
159.499.431 datapunten
Grafiek
Figuur 68
De grafiek van φ (x, y)
voor J = 1.8,
E = 0.954727796245619,
D = 0,
n1 = 0.152,
1/n1 = 6.584,
n2/3 = 0.424,
1/n2/3 = 2.358,
q'' = 7.467 ∙ 10−4,
Rs = 1,
6 rondjes om het zwarte gat,
61.426.903 datapunten
Grafiek
Figuur 69
De grafiek van φ (x, y)
voor J = 1.7321,
E = 0.9428157863474071,
D = 0,
n1 = 0.328,
1/n1 = 3.046,
n2/3 = 0.336,
1/n2/3 = 2.978,
q'' = 1.585 ∙ 10−8,
Rs = 1,
30.5 rondjes om het zwarte gat,
475.377.539 datapunten
Wat ik in de voorgaande stappen heb gedaan is het impulsmoment steeds verder verlaagd, waarbij ik E zo gekozen heb dat D = 0 en q'' > 0. Zie het onderstaande plaatje.
Grafiek
Figuur 70
De grafiek van E (J), D = 0, Rs = 1,
met daarin de punten (de groene bolletjes) van de voorgaande figuren

Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dan krijgen we nu nog een categorie banen ‘aan de andere kant’ van het raakpunt n2/3, oftewel de situatie rechts van het raakpunt, oftewel u > n2/3, oftewel het blauwe deel in figuur 60, oftewel de oranje balk in figuur 61. Dit zijn weer banen waarbij de baksteen in het zwarte gat verdwijnt. Ik kan nog steeds de integraal exact oplossen (want we hebben twee samenvallende nulpunten), maar waar ik in vergelijking (64) de plustekens gebruikte gebruik ik nu de mintekens:
Vergelijking
Voor u vul ik weer 1/r in:
Vergelijking
Hiermee wordt vergelijking (68):
Vergelijking
Voor het verschil van de nulpunten geldt nog steeds:
Vergelijking
Vergelijking
Vergelijking
Hiermee wordt vergelijking (69):
Vergelijking
Die oneindig in vergelijking (70) is slechts een beginwaarde en daar ga ik gewoon een groot getal voor invullen. Dan komen nu uiteraard weer de plaatjes. Ik begin met J = 4 en ik ga het impulsmoment in stappen verlagen.
Grafiek
Figuur 71
De grafiek van φ (x, y)
voor J = 4,
E = 1.6491982708836967,
D = 0,
n1 = −0.268,
n2/3 = 0.634,
1/n2/3 = 1.578,
q'' = 2.713 ∙ 10−2,
Rs = 1,
2 rondjes om het zwarte gat,
400.000.030 datapunten
Grafiek
Figuur 72
De grafiek van φ (x, y)
voor J = 2.5,
E = 1.141591103064672,
D = 0,
n1 = −0.147,
n2/3 = 0.574,
1/n2/3 = 1.743,
q'' = 1.389 ∙ 10−2,
Rs = 1,
2.5 rondjes om het zwarte gat,
400.000.022 datapunten
Grafiek
Figuur 73
De grafiek van φ (x, y)
voor J = 2.1,
E = 1.025861056739398,
D = 0,
n1 = −0.044,
n2/3 = 0.522,
1/n2/3 = 1.916,
q'' = 6.696 ∙ 10−3,
Rs = 1,
2.5 rondjes om het zwarte gat,
399.999.997 datapunten
Grafiek
Figuur 74
De grafiek van φ (x, y)
voor J = 1.9,
E = 0.9760369650397409,
D = 0,
n1 = 0.059,
1/n1 = 16.866,
n2/3 = 0.470,
1/n2/3 = 2.126,
q'' = 2.573 ∙ 10−3,
Rs = 1,
3.5 rondjes om het zwarte gat,
399.999.999 datapunten
Grafiek
Figuur 75
De grafiek van φ (x, y)
voor J = 1.8,
E = 0.954727796245619,
D = 0,
n1 = 0.152,
1/n1 = 6.584,
n2/3 = 0.424,
1/n2/3 = 2.358,
q'' = 7.467 ∙ 10−4,
Rs = 1,
4 rondjes om het zwarte gat,
400.000.012 datapunten
Grafiek
Figuur 76
De grafiek van φ (x, y)
voor J = 1.7321,
E = 0.9428157863474071,
D = 0,
n1 = 0.328,
1/n1 = 3.046,
n2/3 = 0.336,
1/n2/3 = 2.978,
q'' = 1.585 ∙ 10−8,
Rs = 1,
20.5 rondjes om het zwarte gat,
399.999.995 datapunten
De serie plaatjes hiervoor had dezelfde parameters als deze laatste serie, maar eerder was dat voor u < n2/3 (de figuren 64 tot en met 69) en nu voor u > n2/3 (de figuren 71 tot en met 76).
Grafiek
Figuur 77
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D = 0, Rs = 1
Het is wel instructief om de figuren die paarsgewijs vermeld staan in figuur 77 samen te brengen (per paar) in één grafiek.
Grafiek
Figuur 78
(= figuur 64 (de blauwe lijn) + figuur 71 (de rode lijn))
De grafiek van φ (x, y)
voor J = 4,
E = 1.6491982708836967,
D = 0,
n1 = −0.268,
n2/3 = 0.634,
1/n2/3 = 1.578,
q'' = 2.713 ∙ 10−2,
Rs = 1,
427.365.628 datapunten
Grafiek
Figuur 79
(= figuur 65 (de blauwe lijn) + figuur 72 (de rode lijn))
De grafiek van φ (x, y)
voor J = 2.5,
E = 1.141591103064672,
D = 0,
n1 = −0.147,
n2/3 = 0.574,
1/n2/3 = 1.743,
q'' = 1.389 ∙ 10−2,
Rs = 1,
429.187.892 datapunten
Grafiek
Figuur 80
(= figuur 66 (de blauwe lijn) + figuur 73 (de rode lijn))
De grafiek van φ (x, y)
voor J = 2.1,
E = 1.025861056739398,
D = 0,
n1 = −0.044,
n2/3 = 0.522,
1/n2/3 = 1.916,
q'' = 6.696 ∙ 10−3,
Rs = 1,
431.099.518 datapunten
Grafiek
Figuur 81
(= figuur 67 (de blauwe lijn) + figuur 74 (de rode lijn))
De grafiek van φ (x, y)
voor J = 1.9,
E = 0.9760369650397409,
D = 0,
n1 = 0.059,
1/n1 = 16.866,
n2/3 = 0.470,
1/n2/3 = 2.126,
q'' = 2.573 ∙ 10−3,
Rs = 1,
559.499.430 datapunten
Grafiek
Figuur 82
(= figuur 68 (de blauwe lijn) + figuur 75 (de rode lijn))
De grafiek van φ (x, y)
voor J = 1.8,
E = 0.954727796245619,
D = 0,
n1 = 0.152,
1/n1 = 6.584,
n2/3 = 0.424,
1/n2/3 = 2.358,
q'' = 7.467 ∙ 10−4,
Rs = 1,
461.426.915 datapunten
Grafiek
Figuur 83
(= figuur 69 (de blauwe lijn) + figuur 76 (de rode lijn))
De grafiek van φ (x, y)
voor J = 1.7321,
E = 0.9428157863474071,
D = 0,
n1 = 0.328,
1/n1 = 3.046,
n2/3 = 0.336,
1/n2/3 = 2.978,
q'' = 1.585 ∙ 10−8,
Rs = 1,
875.377.534 datapunten
De blauwe banen zijn weliswaar stabiel, maar daar is ook gelijk alles mee gezegd. De integrand bestaat maar op één punt niet, het raakpunt n2/3, en dat is de oneindig dunne scheidslijn tussen de stabiele banen (blauw) en de banen die in het zwarte gat eindigen (rood). Dit is het raakpunt n2/3:
Vergelijking
De reciproke hiervan is de afstand tot het middelpunt van het zwarte gat:
Vergelijking
Merk op dat dit exact gelijk is aan de minoplossing van de extrema van de effectieve potentiaal:
Vergelijking
Daarom toon ik nogmaals figuur 14.
Grafiek
Figuur 14
De grafiek van V (r) voor J = 0.0 (de rode lijn) oplopend tot J = 4.5 (de lichtblauwe lijn),
de zwarte punten zijn de extrema volgens vergelijking (39),
Rs = 1
Daar waar het raakpunt optreedt bevindt zich precies een potentiaaltop. Dat de baksteen volgens een blauwe baan om het zwarte gat cirkelt is net zo stabiel als een potlood dat op zijn punt balanceert, en de minste of geringste verstoring zal de baksteen alsnog in het zwarte gat doen storten.

Ik kan mijn tabel met conclusies weer uitbreiden.
Condities Baan van de baksteen
J = 0 De baksteen valt radieel in het zwarte gat
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
D < 0 E ≥ 1 De baksteen eindigt in het zwarte gat vanaf ‘ver weg’
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n1 < u < 1/Rs oftewel Rs < r < 1/n1)
Grafiek
Grafiek
Grafiek
D = 0 q'' < 0 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n1 < u < 1/Rs oftewel Rs < r < 1/n1)
Grafiek
Grafiek
Grafiek
q'' > 0 E ≥ 1 u > n2/3 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n2/3 < u < 1/Rs oftewel Rs < r < 1/n2/3)
Grafiek
Grafiek
Grafiek
u < n2/3 De baksteen eindigt in een baan om het zwarte gat vanaf ‘dichtbij’
(n1 < u < n2/3 oftewel 1/n2/3 < r < 1/n1)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen eindigt in een baan om het zwarte gat vanaf ‘ver weg’
(0 < u < n2/3 oftewel 1/n2/3 < r < ∞)
Grafiek
Grafiek
Grafiek
Tabel 7
Voor een echte stabiele baan zal de scheidslijn tussen stabiele banen en invallende banen breder moeten zijn dan één enkel punt (het raakpunt n2/3). Het raakpunt moet als het ware uit elkaar getrokken worden tot twee aparte nulpunten en dat gebeurt wanneer de discriminant groter is dan nul.

Invallende banen vanaf ‘dichtbij’ (D > 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Indien D > 0 dan zijn er drie nulpunten als volgt (en vanaf nu ook telkens in deze volgorde, dus n1 is het linkernulpunt, n2 het middelste nulpunt en n3 het rechternulpunt).
Grafiek
Figuur 84
De grafiek van g (u) voor D > 0
De integrand bestaat tussen n1 en n2 en rechts van n3. Schematisch gezien voorspelt dit weer drie categorieën banen.
Grafiek
Figuur 85
Het traject waar de integrand f (u) bestaat (de oranje balk),
D > 0, E ≥ 1, u > n3
Grafiek
Figuur 86
Het traject waar de integrand f (u) bestaat (de oranje balk),
D > 0, E ≥ 1, u < n2
Grafiek
Figuur 87
Het traject waar de integrand f (u) bestaat (de oranje balk),
D > 0, E < 1, n1 < u < n2
De drie nulpunten kan ik uitrekenen volgens de vergelijkingen (45) en de hulpvariabelen die hierin voorkomen had ik al uitgerekend (zie de vergelijkingen (46)). Dus dat ga ik doen, ik ga de nulpunten uitrekenen:
Vergelijking
Vergelijking
Vergelijking
Ik begin met de eerste categorie, u > n3, en de integraal zoek ik op in de tabel met integralen:
Vergelijking
Waarbij voor h geldt:
Vergelijking
Net als in de situatie D < 0 is deze oplossing bereikt door de integraal om te schrijven naar een elliptische integraal van de eerste soort. Om al die termen mee te nemen wordt weer een kolossale klus (als het al überhaupt lukt) en daarom ga ik ook hier gebruik maken van de trapeziummethode:
Vergelijking
Bij ieder volgend plaatje vermeld ik daarom het aantal intervallen dat ik gebruikt heb in de berekening, en ik stel φ0 = 0.

Binnen het domein “D > 0” ga ik een aantal punten eruit pikken en daarvoor de baan van de baksteen uitrekenen. Omdat u > n3 en dus 1/n3 < r < Rs zijn dit allemaal trajecten waarbij de baksteen in het zwarte gat eindigt (dit project kost een hoop bakstenen...).
Grafiek
Figuur 88
De grafiek van φ (x, y)
voor E = 1.0,
J = 5.0,
D = 1.344 ∙ 10−3,
n1 = 5.551 ∙ 10−17,
1/n1 = 1.801 ∙ 1016,
n2 = 0.0417,
1/n2 = 23.956,
n3 = 0.958,
1/n3 = 1.044,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 89
De grafiek van φ (x, y)
voor E = 1.5,
J = 5.0,
D = 9.784 ∙ 10−2,
n1 = −0.189,
n2 = 0.297,
1/n2 = 3.372,
n3 = 0.892,
1/n3 = 1.121,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 90
De grafiek van φ (x, y)
voor E = 2.0,
J = 5.0,
D = 6.144 ∙ 10−3,
n1 = −0.290,
n2 = 0.600,
1/n2 = 1.667,
n3 = 0.690,
1/n3 = 1.449,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 91
De grafiek van φ (x, y)
voor E = 1.0,
J = 4.0,
D = 2.930 ∙ 10−3,
n1 = 1.665 ∙ 10−16,
1/n1 = 6.005 ∙ 1015,
n2 = 0.0670,
1/n2 = 14.928,
n3 = 0.933,
1/n3 = 1.072,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 92
De grafiek van φ (x, y)
voor E = 1.5,
J = 4.0,
D = 6.274 ∙ 10−2,
n1 = −0.228,
n2 = 0.428,
1/n2 = 2.334,
n3 = 0.800,
1/n3 = 1.251,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 93
De grafiek van φ (x, y)
voor E = 1.1,
J = 3.0,
D = 3.883 ∙ 10−2,
n1 = −0.104,
n2 = 0.271,
1/n2 = 3.696,
n3 = 0.833,
1/n3 = 1.201,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 94
De grafiek van φ (x, y)
voor E = 1.0,
J = 2.1,
D = 4.780 ∙ 10−3,
n1 = 0,
1/n1 = ∞,
n2 = 0.348,
1/n2 = 2.877,
n3 = 0.652,
1/n3 = 1.533,
Rs = 1,
999.999 intervallen
Grafiek
Figuur 95
De grafiek van φ (x, y)
voor E = 0.98,
J = 2.0,
D = 2.304 ∙ 10−3,
n1 = 0.0486,
1/n1 = 20.584;,
n2 = 0.326,
1/n2 = 3.071,
n3 = 0.626,
1/n3 = 1.598,
Rs = 1,
999.999 intervallen
Zo zijn we al heel wat bakstenen kwijtgeraakt. Anderzijds, door de voorgaande reeks plaatjes samen te brengen in één grafiek (met verschillende kleuren) wordt het bijna kunst.
Grafiek
Figuur 96
De grafieken van φ (x, y)
voor D > 0,
u > n3,
Rs = 1,
999.999 intervallen (per lijn)
Ik zal in figuur 9 aangeven voor welke punten binnen het domein “D > 0” ik de baan van de baksteen heb uitgerekend.
Grafiek
Figuur 97
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D > 0, u > n3, Rs = 1
Aan het verlies van bakstenen komt nu (eindelijk) een eind, want in de volgende hoofdstukken gaan we de stabiele banen onderzoeken.

Stabiele open banen (D > 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
De laatste twee categorieën die aan bod komen zijn stabiele banen waarbij we eerst kijken naar de zogenaamde open banen, zie nogmaals figuur 86.
Grafiek
Figuur 86
Het traject waar de integrand f (u) bestaat (de oranje balk),
D > 0, E ≥ 1, u < n2
Een open baan wil zeggen dat de baksteen aan komt vliegen vanuit ver weg, eenmalig om het zwarte gat heen beweegt (in meer of mindere mate, want de baksteen kan alleen een beetje worden afgebogen of wellicht zelfs een heel rondje maken), en vervolgens weer vertrekt naar ver weg om nooit meer terug te keren. Omdat E ≥ 1 is n1 ≤ 0 en beweegt de baksteen zich van r = ∞ naar r = 1/n2 en daarna weer terug naar r = ∞. Ik begin met J = 5 en ik varieer de energie.
Grafiek
Figuur 98
De grafiek van φ (x, y)
voor E = 1.01 (de rode lijn),
E = 1.1 (de blauwe lijn),
E = 1.2 (de groene lijn),
E = 1.3 (de oranje lijn),
E = 1.4 (de roze lijn),
E = 1.5 (de lichtblauwe lijn),
E = 1.6 (de gele lijn),
E = 1.7 (de paarse lijn),
E = 1.8 (de bruine lijn)
en E = 1.9 (de grijze lijn),
J = 5,
Rs = 1,
1.999.997 intervallen (per lijn)
Door de energie E verder op te voeren draait de baksteen zelfs om het zwarte gat heen voordat ie weer vertrekt.
Grafiek
Figuur 99
De grafiek van φ (x, y)
voor E = 2.0 (de rode lijn),
E = 2.01 (de blauwe lijn),
E = 2.011 (de groene lijn),
E = 2.0117 (de oranje lijn)
en E = 2.01177 (de grijze lijn),
J = 5,
Rs = 1,
1.999.997 intervallen (per lijn)
Ik ga even inzoomen op de omgeving van het zwarte gat en ik vermeld het aantal rondjes die de baksteen om het zwarte gat draait (bezien vanaf r = +∞ tot r = −∞).
Grafiek
Figuur 100
De grafiek van φ (x, y)
voor E = 2.0 (de rode lijn, 0.5 rondje),
E = 2.01 (de blauwe lijn, 1 rondje),
E = 2.011 (de groene lijn, 1 rondje),
E = 2.0117 (de oranje lijn, 1.5 rondjes)
en E = 2.01177 (de grijze lijn, 2 rondjes),
J = 5,
Rs = 1,
1.999.997 intervallen (per lijn)
Zoals je in figuur 98 ziet komt de baksteen niet telkens op dezelfde ‘hoogte’ aanvliegen, maar verschilt dat afhankelijk van de waarde van E. Deze ‘hoogte’ is de y-waarde voor r = ∞ en die ga ik uitrekenen. Voor y geldt:
Vergelijking
Voor r = ∞ wordt dit:
Vergelijking
Die φ0 kan ik uitrekenen door het eerste trapezium uit te rekenen van de trapeziummethode.
Grafiek
Figuur 101
Trapeziummethode
Dat doe ik als volgt:
Vergelijking
Hierin is ∆ de breedte van één trapezium en die wil ik uiteraard minimaal hebben:
Vergelijking
Dit vul ik in in vergelijking (76):
Vergelijking
De functiewaarden kan ik uitrekenen volgens vergelijking (52) en zo vind ik y:
Vergelijking
De y-waarde die ik zojuist gevonden heb heet de impactparameter (of in beter Nederlands: inslagparameter, in de literatuur vaak aangegeven met de letter b) en het is de loodrechte afstand van ‘iets’ (in dit geval de baksteen) tot de as die door het middelpunt van het aantrekkende (of afstotende) object gaat (in dit geval het zwarte gat) op het moment dat dat ‘iets’ nog op grote afstand is van het object. Indien ik de bakstenen telkens met dezelfde impactparameter wil laten aanvliegen dan moet voor twee banen gelden:
Vergelijking
Ik maak figuur 98 nog een keer, maar dan met een gecorrigeerd impulsmoment volgens vergelijking (81). De baksteen komt nu telkens met de dezelfde impactparameter aanvliegen, maar zowel energie als impulsmoment nemen toe.
Grafiek
Figuur 102
De grafiek van φ (x, y)
voor E = 1.1, J = 5.00 (de blauwe lijn),
E = 1.2, J = 7.24 (de groene lijn),
E = 1.3, J = 9.06 (de oranje lijn),
E = 1.4, J = 10.69 (de roze lijn),
E = 1.5, J = 12.20 (de lichtblauwe lijn),
E = 1.6, J = 13.63 (de gele lijn),
E = 1.7, J = 15.00 (de paarse lijn),
E = 1.8, J = 16.33 (de bruine lijn)
en E = 1.9, J = 17.63 (de grijze lijn),
Rs = 1,
1.999.997 intervallen (per lijn)
Omdat de baksteen steeds meer energie (lees: snelheid) én impulsmoment (lees: snelheid) krijgt ontstaat steeds meer het beeld dat de baksteen rechtdoor vliegt. Logisch, de zwaartekracht van het zwarte gat krijgt steeds minder vat op de baksteen.

De volgende plaatjes zijn met een constant impulsmoment en een variërende energie. Ik begin met J = 4 en laat die vervolgens stapsgewijs dalen (en om de plaatjes overzichtelijk te houden laat ik de bakstenen tijdelijk weg).
Grafiek
Figuur 103
De grafiek van φ (x, y)
voor E = 1.01 (de rode lijn),
E = 1.1 (de blauwe lijn),
E = 1.2 (de groene lijn),
E = 1.3 (de oranje lijn),
E = 1.4 (de roze lijn),
E = 1.5 (de lichtblauwe lijn),
E = 1.6 (de gele lijn),
E = 1.64 (de paarse lijn),
E = 1.649 (de bruine lijn),
E = 1.6491 (de grijze lijn)
en E = 1.64919 (de lichtgroene lijn),
J = 4,
Rs = 1,
1.999.997 intervallen (per lijn)
Grafiek
Figuur 104
De grafiek van φ (x, y)
voor E = 1.01 (de rode lijn),
E = 1.1 (de blauwe lijn),
E = 1.2 (de groene lijn),
E = 1.3 (de oranje lijn),
E = 1.302 (de roze lijn),
E = 1.3023 (de lichtblauwe lijn),
en E = 1.30239 (de gele lijn),
J = 3,
Rs = 1,
1.999.997 intervallen (per lijn)
Grafiek
Figuur 105
De grafiek van φ (x, y)
voor E = 1.01 (de rode lijn),
E = 1.05 (de blauwe lijn),
E = 1.10 (de groene lijn),
E = 1.14 (de oranje lijn),
E = 1.141 (de roze lijn),
E = 1.1415 (de lichtblauwe lijn),
en E = 1.14159 (de gele lijn),
J = 2.5,
Rs = 1,
1.999.997 intervallen (per lijn)
Ik zal in figuur 9 weer aangeven voor welke punten binnen het domein “D > 0” ik de baan van de baksteen heb uitgerekend.
Grafiek
Figuur 106
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D > 0, E ≥ 1, u < n2, Rs = 1
Dan kunnen we ons nu opmaken voor de laatste etappe van deze lange pagina.

Stabiele gesloten banen (D > 0)


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Tot slot komen we bij de banen volgens figuur 87, de stabiele gesloten banen.
Grafiek
Figuur 87
Het traject waar de integrand f (u) bestaat (de oranje balk),
D > 0, E < 1, n1 < u < n2
Het feit dat n1 > 0 impliceert dat E < 1. Voor de duidelijkheid haal ik figuur 9 er weer even bij.
Grafiek
Figuur 9
De grafiek van E (J),
D = 0, Rs = 1
In dit hoofdstuk, de stabiele gesloten banen, zijn we alleen geïnteresseerd in het deel dat E < 1, dus ik ga even inzoomen op het onderste deel van de grafiek.
Grafiek
Figuur 107
De grafiek van E (J),
D = 0, Rs = 1
Ik begin met J = 5.
Grafiek
Figuur 108
De grafiek van φ (x, y)
voor E = 0.995 (de rode lijn),
E = 0.996 (de blauwe lijn),
E = 0.997 (de groene lijn),
E = 0.998 (de oranje lijn),
en E = 0.999 (de paarse lijn),
J = 5,
Rs = 1,
1.999.997 intervallen (per lijn)
Het eerste dat opvalt is dat naarmate de waarde van E stijgt dat de baan dan naar een open baan neigt. Logisch, zie het vorige hoofdstuk, want voor E ≥ 1 is de baan open. Het tweede dat opvalt is dat naarmate de waarde van E daalt dat de baan dan steeds cirkelvormiger wordt, zie het volgende plaatje.
Grafiek
Figuur 109
De grafiek van φ (x, y)
voor E = 0.994882114,
J = 5,
Rs = 1,
1.999.997 intervallen
Maar ook hier is een grens aan, want indien ik E nog een fractie verlaag dan wordt D ≤ 0 en kan de baan niet meer bestaan. In bovenstaande figuur is E = 0.994882114, de discriminant is nog slechts 2.294 ∙ 10−10 en de nulpunten n1 en n2 schuren al tegen elkaar aan (n1 = 0.0206303 en n2 = 0.0206476) en gaan bijna samen een raakpunt vormen. Wanneer de nulpunten n1 en n2 zouden kunnen samenvallen dan werd de baan perfect cirkelvormig, maar dat punt is dus onbereikbaar. En de baan is dan weliswaar niet cirkelvormig, de baan is wel stabiel.

Het derde dat opvalt in figuur 108 is duidelijk zichtbaar bij de blauwe lijn: de baan sluit niet na één omwenteling om het zwarte gat, de baan precesseert! Ik pik die blauwe baan er even uit en laat de baksteen tien rondjes om het zwarte gat maken hetgeen de precessie overduidelijk zichtbaar maakt.
Grafiek
Figuur 110
De grafiek van φ (x, y)
voor E = 0.996,
J = 5,
Rs = 1,
10 rondjes om het zwarte gat,
19.999.979 intervallen
En met vijftig rondjes ontstaat deze prachtige figuur.
Grafiek
Figuur 111
De grafiek van φ (x, y)
voor E = 0.996,
J = 5,
Rs = 1,
50 rondjes om het zwarte gat,
99.999.899 intervallen
Nu met J = 4.
Grafiek
Figuur 112
De grafiek van φ (x, y)
voor E = 0.992 (de rode lijn),
E = 0.993 (de blauwe lijn),
E = 0.994 (de groene lijn),
E = 0.995 (de oranje lijn),
E = 0.996 (de paarse lijn),
E = 0.997 (de roze lijn),
E = 0.998 (de lichtblauwe lijn)
en E = 0.999 (de gele lijn),
J = 4,
Rs = 1,
1.999.997 intervallen (per lijn)
In alle gevallen treedt ook hier precessie op. Ik pik nu de oranje baan er uit en laat de baksteen tien rondjes om het zwarte gat maken.
Grafiek
Figuur 113
De grafiek van φ (x, y)
voor E = 0.995,
J = 4,
Rs = 1,
10 rondjes om het zwarte gat,
19.999.979 intervallen
En met vijftig rondjes ontstaat weer een prachtige figuur.
Grafiek
Figuur 114
De grafiek van φ (x, y)
voor E = 0.995,
J = 4,
Rs = 1,
50 rondjes om het zwarte gat,
99.999.899 intervallen
Ik ga het impulsmoment verder verlagen:
Grafiek
Figuur 115
De grafiek van φ (x, y)
voor E = 0.986 (de rode lijn),
E = 0.988 (de blauwe lijn),
E = 0.990 (de groene lijn),
E = 0.992 (de oranje lijn),
E = 0.994 (de paarse lijn),
E = 0.996 (de roze lijn)
en E = 0.998 (de lichtblauwe lijn),
J = 3,
Rs = 1,
1.999.997 intervallen (per lijn)
Grafiek
Figuur 116
De grafiek van φ (x, y)
voor E = 0.978 (de rode lijn),
E = 0.980 (de blauwe lijn),
E = 0.982 (de groene lijn),
E = 0.984 (de oranje lijn),
E = 0.986 (de paarse lijn),
E = 0.988 (de roze lijn),
E = 0.990 (de lichtblauwe lijn),
E = 0.992 (de gele lijn),
E = 0.994 (de lichtgroene lijn),
E = 0.996 (de bruine lijn)
en E = 0.998 (de grijze lijn),
J = 2.5,
Rs = 1,
1.999.997 intervallen (per lijn)
Grafiek
Figuur 117
De grafiek van φ (x, y)
voor E = 0.963 (de rode lijn),
E = 0.969 (de blauwe lijn),
E = 0.975 (de groene lijn),
E = 0.981 (de oranje lijn),
E = 0.987 (de paarse lijn),
E = 0.993 (de roze lijn)
en E = 0.999 (de lichtblauwe lijn),
J = 2,
Rs = 1,
1.999.997 intervallen (per lijn)
Ik zal in figuur 107 aangeven voor welke punten binnen het domein “D > 0” ik de baan van de baksteen heb uitgerekend.
Grafiek
Figuur 118
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D > 0, E < 1, u < n2, Rs = 1
Ter lering (en vermaak, voor de mooie plaatjes) is het wel leuk om het domein “D > 0” nog iets verder te verkennen voor lagere waarden van J. Merk op dat de baksteen steeds meer rondjes om het zwarte gat draait naarmate D naar nul nadert. En merk ook op dat er in dit domein geen open banen meer bestaan, of de neiging daartoe, omdat de waarde E = 1 niet bereikt kan worden (onder de voorwaarde dat D > 0 wel te verstaan).
Grafiek
Figuur 119
De grafiek van φ (x, y)
voor E = 0.97,
J = 1.9,
Rs = 1,
1.8 rondjes om het zwarte gat,
1.999.997 intervallen
Grafiek
Figuur 120
De grafiek van φ (x, y)
voor E = 0.97,
J = 1.9,
Rs = 1,
10 × 1.8 = 18 rondjes om het zwarte gat,
19.999.979 intervallen
Grafiek
Figuur 121
De grafiek van φ (x, y)
voor E = 0.976,
J = 1.9,
Rs = 1,
3 rondjes om het zwarte gat,
1.999.997 intervallen
Grafiek
Figuur 122
De grafiek van φ (x, y)
voor E = 0.976,
J = 1.9,
Rs = 1,
10 × 3 = 30 rondjes om het zwarte gat,
19.999.979 intervallen
Grafiek
Figuur 123
De grafiek van φ (x, y)
voor E = 0.97603696,
J = 1.9,
Rs = 1,
5 rondjes om het zwarte gat,
1.999.997 intervallen
Grafiek
Figuur 124
De grafiek van φ (x, y)
voor E = 0.97603696,
J = 1.9,
Rs = 1,
10 × 5 = 50 rondjes om het zwarte gat,
19.999.979 intervallen
Als uitsmijter heb ik een baan uitgezocht waarbij de precessie minimaal is, dat levert deze plaatjes op.
Grafiek
Figuur 125
De grafiek van φ (x, y)
voor E = 0.9704,
J = 1.888,
Rs = 1,
2 rondjes om het zwarte gat,
1.999.997 intervallen
Grafiek
Figuur 126
De grafiek van φ (x, y)
voor E = 0.9704,
J = 1.888,
Rs = 1,
10 × 2 = 20 rondjes om het zwarte gat,
19.999.979 intervallen
Grafiek
Figuur 127
De grafiek van φ (x, y)
voor E = 0.9704,
J = 1.888,
Rs = 1,
50 × 2 = 100 rondjes om het zwarte gat,
99.999.899 intervallen
Voor de laatste reeks plaatjes heb ik deze punten gebruikt binnen het domein “D > 0”.
Grafiek
Figuur 128
De grafiek van E (J) met daarin de punten (de groene bolletjes)
van de voorgaande figuren,
D > 0, E < 1, u < n2, Rs = 1

Belangrijke vergelijkingen


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Dit is een overzicht met de belangrijkste vergelijkingen die de revue gepasseerd zijn.

De Schwarzschild-oplossing:
Vergelijking
Het impulsmoment J (van de baksteen):
Vergelijking
De energie E (van de baksteen):
Vergelijking
De Schwarzschild-straal:
Vergelijking
De variabele u is de reciproke van de radiële coördinaat r:
Vergelijking
Het kernprobleem van deze pagina, de integraal die φ als functie van u beschrijft:
Vergelijking
De integrand f (u):
Vergelijking
De derdegraads vergelijking g (u) in de noemer van f (u):
Vergelijking
De discriminant D van g (u) als functie van E2 respectievelijk J2:
Vergelijking
Vergelijking
E als functie van J en J als functie van E voor D = 0:
Vergelijking
Vergelijking
De snelheid van de baksteen (voor een waarnemer die zich ver van het zwarte gat bevindt):
Vergelijking
De effectieve potentiaal:
Vergelijking
De extrema van de effectieve potentiaal:
Vergelijking
Vier hulpvariabelen om de nulpunten van de derdegraads vergelijking uit te rekenen:
Vergelijking
Vergelijking
Vergelijking
Vergelijking
De waarden van deze hulpvariabelen voor dit specifieke probleem (deze pagina) zijn:
Vergelijking
Vergelijking
Vergelijking
Vergelijking
Het nulpunt n1 indien D < 0:
Vergelijking
Het nulpunt n1 indien D = 0:
Vergelijking
Het raakpunt n2/3 indien D = 0:
Vergelijking
De oplossing van de integraal voor D = 0 én n2/3 < n1:
Vergelijking
De oplossing van de integraal voor D = 0 én n2/3 > n1 én u < n2/3:
Vergelijking
De oplossing van de integraal voor D = 0 én n2/3 > n1 én u > n2/3:
Vergelijking
De drie nulpunten indien D > 0 (van links naar rechts):
Vergelijking
Vergelijking
Vergelijking
De impactparameter, de y-coördinaat waarop de baksteen aan komt vliegen vanuit r = ∞:
Vergelijking

Samenvatting


  1. Afleidingen van de vergelijkingen
  2. Onderzoek van de functie
  3. De effectieve potentiaal
  4. Nulpunten in de noemer
  5. Radiële inval
  6. Invallende banen vanaf ‘ver weg’ (D < 0)
  7. Invallende banen vanaf ‘dichtbij’ (D < 0)
  8. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' < 0)
  9. Stabiele banen vanaf ‘ver weg’ en ‘dichtbij’ (D = 0)
  10. Invallende banen vanaf ‘dichtbij’ (D = 0, q'' > 0)
  11. Invallende banen vanaf ‘dichtbij’ (D > 0)
  12. Stabiele open banen (D > 0)
  13. Stabiele gesloten banen (D > 0)
  14. Belangrijke vergelijkingen
  15. Samenvatting
Het kernprobleem van deze lange pagina is de integraal die φ als functie van u beschrijft:
Vergelijking
Nog afgezien van het feit dat het oplossen van deze integraal zeer problematisch is kent de noemer één of meerdere nulpunten die allemaal om een specifieke aanpak van het probleem vragen. De discriminant D van de derdegraads vergelijking in de noemer bepaalt het aantal nulpunten van die derdegraads vergelijking:
Grafiek
D < 0
Grafiek
D = 0
Grafiek
D > 0
Tabel 2
Er zijn twee parameters in het spel, de energie E (van de baksteen) en het impulsmoment J (van de baksteen). In het geval dat D = 0 kan ik E als functie van J schrijven of J als functie van E:
Vergelijking
Vergelijking
Grafisch ziet dat er zo uit, voor de blauwe lijn geldt D = 0:
Grafiek
Figuur 9
De grafiek van E (J),
D = 0, Rs = 1
Grafiek
Figuur 10
De grafiek van J (E),
D = 0, Rs = 1
Ik introduceerde een effectieve potentiaal:
Vergelijking
Grafiek
Figuur 13
De grafiek van V (r) voor J = 0.0 (de rode lijn),
J = 0.5 (de lichtgroene lijn), J = 1.0 (de oranje lijn),
J = 1.5 (de paarse lijn), J = √3 (de blauwe lijn),
J = 2.0 (de grijze lijn), J = 2.5 (de bruine lijn),
J = 3.0 (de lichtbruine lijn), J = 3.5 (de donkergroene lijn),
J = 4.0 (de gele lijn) en J = 4.5 (de lichtblauwe lijn),
Rs = 1
Ik vond voor de extrema van de effectieve potentiaal:
Vergelijking
Vergelijking (39) geeft r als functie van J, maar dat kan ik omschrijven naar J als functie van r:
Vergelijking
Door dit in te vullen in vergelijking (37) ontstaat:
Vergelijking
Dit voeg ik toe aan de grafiek van figuur 13.
Grafiek
Figuur 129
De grafiek van V (r) voor J = 0.0 (de rode lijn),
J = 0.5 (de lichtgroene lijn), J = 1.0 (de oranje lijn),
J = 1.5 (de paarse lijn), J = √3 (de blauwe lijn),
J = 2.0 (de grijze lijn), J = 2.5 (de bruine lijn),
J = 3.0 (de lichtbruine lijn), J = 3.5 (de donkergroene lijn),
J = 4.0 (de gele lijn) en J = 4.5 (de lichtblauwe lijn),
de roze lijn wordt gevormd door de extrema volgens vergelijking (83),
Rs = 1
We leerden onderweg ook dat vergelijking (39), de extrema van de effectieve potentiaal, exact overeenkomt met (de reciproke van) het raakpunt n2/3 (waar de integrand niet bestaat):
Vergelijking
Oftewel, de top van de potentiaalfunctie (de roze lijn in figuur 129), het ontstaan van een raakpunt (zie figuur 23) en het gelijk aan nul worden van de discriminant zijn allemaal rechtstreeks aan elkaar gekoppeld.

Hoe werkt dit uit in de praktijk? Voor D < 0 en E ≥ 1 komt de baksteen van ver weg en verdwijnt onherroepelijk in het zwarte gat, zie het rode gebied in onderstaande figuur (zie hoofdstuk 6 voor alle details).
Grafiek
Figuur 130
De grafiek van E (J) (de blauwe lijn),
D = 0, Rs = 1
Voor D < 0 en E < 1 komt de baksteen van dichtbij en verdwijnt eveneens onherroepelijk in het zwarte gat, zie het rode gebied in onderstaande figuur (zie hoofdstuk 7 voor alle details).
Grafiek
Figuur 131
De grafiek van E (J) (de blauwe lijn),
D = 0, Rs = 1
Indien D = 0, de blauwe lijn in bovenstaande figuur, dan kan de baksteen in theorie in een circulaire baan om het zwarte gat komen maar balanceert daarbij wel op de rand van stabiliteit. De minste of geringste verstoring, bijvoorbeeld een botsing met een zandkorrel, zal het lot van de baksteen bezegelen. De circulaire baan bevindt zich net wel of net niet aan de goede kant van de potentiaaltop. In onderstaande figuur komt de baksteen aanvliegen vanuit ver weg en eindigt in een circulaire baan, de blauwe lijn, maar de rode lijn is een circulaire baan van waaruit de baksteen in het zwarte gat stort. Zoals je ziet is de scheidslijn tussen beide banen minimaal (zie hoofdstuk 8, hoofdstuk 9 en hoofdstuk 10 voor alle details).
Grafiek
Figuur 78
(= figuur 64 (de blauwe lijn) + figuur 71 (de rode lijn))
De grafiek van φ (x, y)
voor J = 4,
E = 1.6491982708836967,
D = 0,
n1 = −0.268,
n2/3 = 0.634,
1/n2/3 = 1.578,
q'' = 2.713 ∙ 10−2,
Rs = 1,
427.365.628 datapunten
De stabiele banen treden op voor D > 0. Bevindt de baksteen zich te dicht bij het zwarte gat dan gaat het alsnog mis (zie hoofdstuk 11), maar verder weg van het zwarte gat is de baan stabiel. Indien E ≥ 1 dan is de baan open, de baksteen komt van ver weg en verdwijnt daar uiteindelijk ook weer naar toe, zie het gele gebied in onderstaande figuur (zie hoofdstuk 12 voor alle details).
Grafiek
Figuur 132
De grafiek van E (J) (de blauwe lijn),
D = 0, Rs = 1
In de buurt van het zwarte gat wordt de baksteen (sterk) afgebogen of maakt zelfs één of meerdere rondjes om het zwarte gat.
Grafiek
Figuur 133
De grafiek van φ (x, y)
voor E = 2.01177,
J = 5,
Rs = 1,
2 rondjes om het zwarte gat,
1.999.997 intervallen
Indien E < 1, zie het gele gebied in onderstaande figuur, dan is de baan gesloten (zie hoofdstuk 13 voor alle details).
Grafiek
Figuur 134
De grafiek van E (J) (de blauwe lijn),
D = 0, Rs = 1
Ik zal even inzoomen.
Grafiek
Figuur 135
De grafiek van E (J) (de blauwe lijn),
D = 0, Rs = 1
Deze gesloten banen precesseren, een relativistisch verschijnsel dat onbekend is in de klassieke mechanica.
Grafiek
Figuur 110
De grafiek van φ (x, y)
voor E = 0.996,
J = 5,
Rs = 1,
10 rondjes om het zwarte gat,
19.999.979 intervallen
Dit alles brengt ons bij de afsluitende overzichtstabel.
Condities Baan van de baksteen
J = 0 De baksteen valt radieel in het zwarte gat
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
D < 0 E ≥ 1 De baksteen eindigt in het zwarte gat vanaf ‘ver weg’
(0 < u < 1/Rs oftewel Rs < r < ∞)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n1 < u < 1/Rs oftewel Rs < r < 1/n1)
Grafiek
Grafiek
Grafiek
D = 0 q'' < 0 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n1 < u < 1/Rs oftewel Rs < r < 1/n1)
Grafiek
Grafiek
Grafiek
q'' > 0 E ≥ 1 u > n2/3 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n2/3 < u < 1/Rs oftewel Rs < r < 1/n2/3)
Grafiek
Grafiek
Grafiek
u < n2/3 De baksteen eindigt in een baan om het zwarte gat vanaf ‘dichtbij’
(n1 < u < n2/3 oftewel 1/n2/3 < r < 1/n1)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen eindigt in een baan om het zwarte gat vanaf ‘ver weg’
(0 < u < n2/3 oftewel 1/n2/3 < r < ∞)
Grafiek
Grafiek
Grafiek
D > 0 u > n3 De baksteen eindigt in het zwarte gat vanaf ‘dichtbij’
(n3 < u < 1/Rs oftewel Rs < r < 1/n3)
Grafiek
Grafiek
Grafiek
u < n2 E ≥ 1 De baksteen beweegt in een stabiele open baan om het zwarte gat
(0 < u < n2 oftewel 1/n2 < r < ∞)
Grafiek
Grafiek
Grafiek
E < 1 De baksteen beweegt in een stabiele gesloten baan om het zwarte gat
(n1 < u < n2 oftewel 1/n2 < r < 1/n1)
Grafiek
Grafiek
Grafiek
Tabel 8