Wie kann ich die Gleichungen der klassischen Kontinuumsmechanik anpassen, um mit der Allgemeinen Relativitätstheorie übereinzustimmen?

Ich komme aus der Kontinuumsmechanik und mache numerische Simulationen von Flüssigkeiten/Festkörpern mit der Finite-Elemente-Methode. Die Grundgleichung, die wir dann lösen, ist Newtons Bewegungsgesetz, geschrieben in Bezug auf relevante Vektoren und Tensoren. Unter Verwendung einer Lagrange-Beschreibung haben wir die Gleichgewichtsgleichung:

ρ u ¨ _ = d ich v ( σ _ _ ) + f _ v

wo u ¨ _ ist der Beschleunigungs-3-Vektor (zweite zeitliche Ableitung des Verschiebungsvektors), σ _ _ ist der Cauchy-Spannungstensor zweiter Ordnung und f _ v ist der 3-Vektor der äußeren Kräfte.

Das Gleichungssystem wird durch eine konstitutive Gleichung vervollständigt, die den Spannungstensor mit einem Dehnungsmaß (typische Festkörper) oder einem Dehnungsgeschwindigkeitsmaß (viskoelastische Festkörper und Flüssigkeiten) verknüpft. Betrachten wir zum Beispiel eine übliche linearelastische Beziehung:

σ _ _ = C : ε _ _

mit C ist der Elastizitätstensor 4. Ordnung und ε _ _ der lineare Dehnungstensor, klassisch definiert in der Kompatibilitätsgleichung als der symmetrische Teil des Verschiebungsgradienten:

ε _ _ = 1 2 ( _ _ u + _ _ T u ) .

Nun, die angegebenen Gleichungen sind alles, was ich brauche, um das System zu diskretisieren und kleine Verformungen einer beliebigen Geometrie unter einem Kräftesystem im Kontext der Newtonschen Mechanik mit der Finite-Elemente-Methode zu simulieren.

Was ich wissen möchte, ist, wie man diese Gleichungen so anpasst, dass sie der Allgemeinen Relativitätstheorie genügen. Das heißt, ich möchte die Verformung einer mechanischen Struktur simulieren, wenn die beteiligten Geschwindigkeiten nahe an der Lichtgeschwindigkeit liegen und/oder wenn ein sehr massives Objekt in der Nähe ist. Ich bin bei Bedarf mit nichtlinearer Elastizität vertraut, aber was Tensoren angeht, bin ich mit der kovarianten / kontravarianten Notation nicht vertraut und bevorzuge die intrinsische Notation, obwohl ich in irgendeiner Weise ausgedrückte Antworten nehmen werde.

Welche Form haben die Gleichgewichts-, Konstitutiv- und Kompatibilitätsgleichungen? Ist die Simulation der Verformung von Körpern in einem relativistischen Kontext etwas, das bereits richtig gemacht wurde? Muss der Elastizitätstensor vielleicht in Bezug auf den metrischen Tensor neu definiert werden? Ich konnte keine gute Referenz finden, die dieses Problem anspricht, obwohl ich der Meinung bin, dass dies möglich ist. Ich wäre sehr dankbar für jede Hilfe in dieser Angelegenheit.

Ich kann darauf keine kurze Antwort geben (und ehrlich gesagt nicht einmal eine bestimmte lange), aber ich kann ein paar vielleicht hilfreiche Kommentare geben: Der Schritt von der Newtonschen Mechanik zu GR und sogar SR ist groß ... sehr groß in Bezug auf die Komplexität . Es gibt GR-Ausdrücke für die Gleichungen der Hydrodynamik (rel. Euler-Gleichung, ...) und es gibt vollständige GR-4D-Simulationen von deformierten / kollidierenden Körpern: Neutronensternverschmelzungen, Supernovae, ... Ich würde Bücher und Literatur über Numerik empfehlen Relativitätstheorie und der 3+1/ADM-Formalismus. Zum Beispiel "Numerische Relativitätstheorie - TW Baumgarte und SL Shapiro".
Danke für die Leseempfehlungen. Ich bin neu im GR-Formalismus, also wird das sehr helfen!
Wenn Sie mit diesem Buch keinen Zugang zu einer Bibliothek haben; Dieses Papier arxiv.org/abs/gr-qc/0703035v1 gibt eine sehr schöne Einführung in die numerische Relativitätstheorie (es wurde ein Buch mit demselben Titel). Und dieses Papier adsabs.harvard.edu/abs/1993A%26A...278..421B ist ein schönes Beispiel für eine "einfache" Anwendung des 3+1-Formalismus auf rotierende Neutronensterne mit EM-Feldern. Es zeigt alle notwendigen Gleichungen (GR-Feld, Bewegungsgleichung, Maxwell) im Rahmen der numerischen Relativitätstheorie. Für mehr adv. Simulationen empfehle ich die Arbeit von L. Rezzolla et. al.: Whisky-Code, rel. Hydrodynamik ....

Antworten (2)

Die gesuchte Antwort scheint in
Rezzolla & Zanotti: Relativistic Hydrodynamics (Oxford UP 2013)
https://books.google.com/books/?id=KU2oAAAAQBAJ enthalten zu sein

aber es ist keine triviale Verallgemeinerung. Zitat von Disconzi's On the well-posedness of relativistic viscous fluids (Nonlinearity 27 (2014) 1915, arXiv:1310.1954 ):

Uns fehlt immer noch eine zufriedenstellende Formulierung viskoser Phänomene innerhalb von Einsteins allgemeiner Relativitätstheorie. [... T]hier gab es verschiedene Vorschläge, was das Richtige ist T a β sollte sein. [... Versuche, eine viskose relativistische Theorie zu formulieren, die auf einer einfachen kovarianten Verallgemeinerung des klassischen (dh nicht-relativistischen) Spannungs-Energie-Tensors für die Navier-Stokes-Gleichungen basiert, haben ebenfalls keine kausale Theorie hervorgebracht.

Ähnliche Probleme bestehen für elastische Materialien.

Sie können sich auch Kapitel 15 ("Relativistische Kontinuumsmechanik") von Maugin's Continuum Mechanics Through the Twentieth Century (Springer 2013) und seine Referenzen ansehen:
https://books.google.com/books?id=-QhAAAAAQBAJ

und bei Bressan's Relativistic Theories of Materials (Springer 1978):
https://books.google.com/books?id=kMTuCAAAQBAJ

Die allgemein-relativistische Kontinuumsmechanik hat leider noch keinen klaren mathematischen und begrifflichen Rahmen bekommen. Die Newtonsche Kontinuumsmechanik lässt sich leicht zusammenfassen:

  1. Wir wählen einen Referenzrahmen (vorzugsweise, aber nicht notwendigerweise träge).
  2. Wir haben eine Reihe von 11 raumzeitabhängigen Feldern mit klarer physikalischer Bedeutung: Masse, Impuls oder Verformung, Spannung, Körperkraft, innere Energie, Wärmefluss, Körpererwärmung, Temperatur, Entropie, Entropiefluss, Körperentropiezufuhr. Von diesen repräsentieren die "körperlichen" Eingriffe von außen.
  3. Wir haben 5 Bilanzgleichungen: Masse, Kraft-Impuls, Drehmoment-Drehimpuls, Energie, Entropie. Sie sind eindeutig in Bezug auf die obigen Felder geschrieben und gelten für alle Materialien.
  4. Wir wählen eine Reihe unabhängiger Felder (normalerweise Masse, Impuls oder Verformung, Temperatur).
  5. Wir wählen konstitutive Gleichungen (kompatibel mit den Ausgleichsgleichungen), die die verbleibenden Felder mit den unabhängigen in Beziehung setzen. Diese Gleichungen drücken die besonderen Eigenschaften (flüssig, fest, elastisch, plastisch, mit/ohne Gedächtnis...) des untersuchten Materials aus.

Und an diesem Punkt haben wir einen wohldefinierten Satz partieller Differentialgleichungen in einer Reihe unbekannter Felder, für die wir wohldefinierte Anfangs- und Randwertprobleme aufstellen können, die analytisch oder numerisch gelöst werden. (Ein erweiterter, aber analoger Rahmen nimmt Elektromagnetismus und Kontinua mit interner Struktur auf.)

Dieser Rahmen und diese Schritte sind sehr ordentlich – wir wissen genau, was die Felder sind, welche von ihnen abhängig und welche unabhängig sind; Welche Gleichungen gelten für alle Materialien und welche Gleichungen sind für jedes Material konstitutiv? Ich habe noch nie ein klar definiertes Verfahren wie das obige für die allgemeine Relativitätstheorie gesehen, obwohl ich glaube, dass es aus den Büchern von Rezzolla & Zanotti oder Bressan entnommen werden könnte. Darüber hinaus verwendet der Kern der allgemein-relativistischen Gemeinschaft einen anderen Jargon und eine andere Denkweise.

Die meisten Bücher über allgemeine Relativitätstheorie sagen Ihnen, dass die Einstein-Gleichungen alles bestimmen, aber sie sind nicht so klar darüber, welche Felder in ihnen unabhängig und welche abhängig sind; sogar Gravitation (Kap. 21) von Misner et al. enthält eine lange Diskussion und Erklärung zu diesem Punkt. Erst mit 3+1-Formulierungen und der Arbeit von Arnowitt, Deser, Misner, York und anderen in den 1970er Jahren wurde dieser Punkt geklärt. Dann sagen sie dir, dass wir „spezielle“ Zusatzgleichungen für den Spannungstensor brauchen – also konstitutive Gleichungen. Manchmal werden andere Erhaltungsgleichungen wie die baryonische Zahl (im Grunde Ruhemasse) ohne wirkliche Erklärung hinzugefügt. Dies ist eine Auswahl von Büchern, in denen "konstitutive Gleichungen" explizit erwähnt werden (in den meisten nur ein- oder zweimal):

  • Rezzolla & Zanotti oben (und sie erklären, was eine "konstitutive Gleichung" ist, als wäre es ein exotisches Konzept)
  • Choquet-Bruhat: Allgemeine Relativitätstheorie und Einsteins Gleichungen
  • Anile & Choquet-Bruhat: Relativistische Fluiddynamik
  • Bertotti et al.: Allgemeine Relativitätstheorie und Gravitation
  • Pützfeld et al.: Bewegungsgleichungen in der relativistischen Gravitation
  • Tonti: Die mathematische Struktur der klassischen und relativistischen Physik
  • Bini & Ferrarese: Einführung in die relativistische Kontinuumsmechanik
  • Tolman (offensichtlich): Relativitätstheorie, Thermodynamik und Kosmologie

aber sie bilden eine sehr kleine Minderheit in der riesigen relativistischen Literatur.

Dennoch kann man der allgemein-relativistischen Gemeinschaft nicht den wirren Begriffszustand und die irgendwie wirre Sprache des Subjekts vorwerfen. Die Newtonsche Kontinuumsmechanik kann heute sauber formuliert werden, weil sie über mehrere Jahrhunderte verfeinert wurde. Stattdessen ist die Allgemeine Relativitätstheorie noch sehr jung, und ihre konzeptionelle Verfeinerung ist noch im Gange. Einige der Schritte im Newtonschen Rahmen werden in der Allgemeinen Relativitätstheorie extrem kompliziert. Zum Beispiel: Schritt 1. (Wählen Sie ein Inertialsystem) kann nicht so einfach durchgeführt werden. Die aus Anfangsbedingungen entstandenen Einstein-Gleichungen konstruieren „nebenbei“ einen Bezugsrahmen , indem sie die Dynamik bestimmen. Dadurch entstehen eigentümliche Felder wie "lapse" und "shift", die nicht wirklich physisch sind,

Ein weiteres Beispiel: Die Metrik wird zu einer dynamischen Feldvariablen, und Sie stellen plötzlich fest, dass sie fast überall im Newtonschen Rahmen verborgen ist – Divergenzen, Locken, Vektoren/Kovektoren ... Daher kann ihre Entwicklung nicht einfach zwischen einem neuen Gleichgewicht und aufgeteilt werden konstitutive Gleichungen (wie wir es stattdessen mit Elektromagnetismus tun können). Sind alle seine Erscheinungen im Newtonschen System dynamisch signifikant? oder kann die Metrik an einigen Stellen eliminiert werden? Es gibt heute einige Forschungen zu dieser "Demetrisierung" der Newtonschen Gleichungen; siehe zum Beispiel Segev's Metric-independent analysis of the stress-energy tensor , J. Math. Phys. 43(2002) 3220. Diese Forschungslinie hat gezeigt, dass einige Newtonsche physikalische Objekte eigentlich keine Metrik benötigen: Sie werden durch Differentialformen und andere metrisch-freie differentiell-geometrische Objekte ausgedrückt (z. B. van Dantzigs Über die geometrische Darstellung elementarer Physical Objects and the Relations between Geometry and Physics , Nieuw Archief voor Wiskunde II (1954) 73; dazu gibt es eine umfangreiche Literatur, lassen Sie es mich wissen, wenn Sie weitere Referenzen wünschen). Dies ist noch in Arbeit – was bedeutet, dass offensichtlich noch nicht ganz klar ist, wie Masse-Energie-Impuls-Spannung und Metrik gekoppelt sind.

Abschließend denke ich, dass ein weiterer guter Ausgangspunkt, um zu verstehen, wie die Dinge in der allgemein-relativistischen Kontinuumsmechanik funktionieren, darin besteht, Bücher über numerische Formulierungen der allgemeinen Relativitätstheorie und Materiedynamik zu lesen. Der konzeptionelle Rahmen in ihnen ist etwas verwirrt, aber Sie können sehen, wie sie es tatsächlich tun . Wenn es Ihnen anhand der Praxis dieser Bücher gelingt, ein Framework wie das Newtonsche oben nachzukonstruieren, schreiben Sie bitte eine pädagogische Arbeit darüber!

Hier sind einige Bücher und Rezensionen zur numerischen Relativitätstheorie mit Kontinua:

  • Rezzolla & Zanotti oben
  • Gourgoulhon: 3+1-Formalismus in der Allgemeinen Relativitätstheorie (Springer 2012, arXiv:gr-qc/0703035 )
  • Baumgarte & Shapiro: Numerische Relativitätstheorie (Cambridge UP 2010)
  • Alcubierre: Einführung in die 3+1 Numerische Relativitätstheorie (Oxford UP 2008)
  • Palenzuela-Luque & Bona-Casas: Elemente der numerischen Relativitätstheorie und relativistischen Hydrodynamik (Springer 2009)
  • Lehner: Numerische Relativitätstheorie: eine Rezension , Klasse. Menge Grav. 18 (2001) R25, arXiv:gr-qc/0106072
  • Guzmán: Einführung in die numerische Relativitätstheorie durch Beispiele , Rev. Mex. Fis. S 53 (2007) 78

Gerne gebe ich weitere Referenzen an oder suche nach weiteren Referenzen.

Willkommen bei Physics SE und vielen Dank für Ihren Beitrag – dies muss der beeindruckendste erste Beitrag sein, den ich bisher auf dieser Website gelesen habe!
@Sanya Danke Sanya! Hast du mir ein Bronzeabzeichen gegeben? danke auch in diesem fall dafür! BrouHs Frage ist sehr wichtig, und ich habe sie noch nie so sauber und einfach formuliert gesehen wie er/sie. Ich wünschte, BrouH würde es an pädagogische Zeitschriften wie das American Journal of Physics oder das European Journal of Physics schicken und sehen, welche Antwort es dort erhalten würde ...
@pglm Du hast vollkommen Recht und eine gute pädagogische Behandlung würde viel Ärger machen. Ich habe selbst einen Cont-Mech-Hintergrund und bin kürzlich durch Bini & Ferrarese gesprungen, um nach einem Analogon von konstitutiven Gleichungen zu suchen, aber naja ... Daher war es sehr aufschlussreich, Ihren Beitrag zu lesen :)
@Sanya Ich verstehe, wir stehen vor ähnlichen Problemen. Sie werden feststellen, dass Bini & Ferrarese die Arbeit von Bressan, Maugin oder Toupin nicht erwähnt (nur einmal Truesdell & Toupin), und Maugin erwähnt Bini & Ferrarese nicht. Das gibt eine Vorstellung davon, wie verstreut und isoliert dieses Thema leider ist ... Aber ich kann Bressan wirklich empfehlen, auch wenn es etwas weitläufig ist. Über aufschlussreiche Texte, die auch Sie finden, würde ich mich freuen!
Um ehrlich zu sein, habe ich hauptsächlich aus Interesse nachgeschlagen und nicht, weil ich daran arbeite. Wenn ich Ihre Leseliste sehe, fürchte ich, dass ich nicht viel beitragen kann, aber ich werde es Sie wissen lassen. Sie können mich oft im Chat chat.stackexchange.com/rooms/71/the-h-bar finden, wenn Sie ein bisschen reden möchten :)
Ich werde mich dann mal dem Chat anschließen :) Ich habe eine neue Referenz hinzugefügt, sie scheint die vielversprechendste zu sein (muss sie selbst lesen).
Vielen Dank für eine so gründliche Antwort! Sie haben nicht nur Referenzen gegeben, die meine Fragen beantworten könnten, sondern Sie haben jedes Konzept der Newtonschen Kontinuumsmechanik, über das ich anfangs gesprochen habe, erfolgreich verstanden und Ihre Antwort darauf gegeben. Danke auch für den Hinweis, dass diese unterschiedlichen Teilgebiete der Physik möglicherweise (noch!) nicht in vollem Umfang zu einem produktiven Dialog fähig sind. Ich wechsle zu dieser Antwort als akzeptierte, weil sie die ganze Frage hervorragend abdeckt. Nochmals vielen Dank dafür und dafür, dass Sie mir einen neuen Stapel Bücher zum Lesen gegeben haben :-)
Danke, BroH! Es scheint, dass es viele nette Leute in dieser Community gibt :) Tut mir leid, dass ich keine konkrete Antwort gegeben habe. Bitte erwägen Sie, Ihre Frage als Brief an das American Journal of Physics oder das European Journal of Physics zu senden . Früher gab es in jedem Band einen eigenen Abschnitt für diese Art von Fragen, und ich denke, Ihrer ist sehr wichtig.

Wie Sie bereits erwähnt haben, ist der Schlüssel zu allen Modifikationen die Metrik, denn die Metrik ist alles, was Sie brauchen, um die Hintergrundraumzeit zu charakterisieren, sei sie flach oder gekrümmt.

Wenn Sie die abstrakte Indexnotation übernehmen, kann die Metrik geschrieben werden als g a b , ein symmetrischer Tensor zweiten Ranges, der zwei Vektoren aufnimmt und einen Skalar ausgibt. Schauen wir uns nun einige allgemeine Beispiele für die Umwandlung von flacher in gekrümmte Raumzeit an, bevor wir uns mit der speziellen Bewegungsgleichung befassen, die Sie interessiert.

Die Erhaltung des Energieimpulses in der flachen Raumzeit für einen beliebigen Spannungsenergietensor T kann ausgedrückt werden als:

a T a b = 0
wobei die Einstein-Summierungskonvention durch die Kontraktion über den Index impliziert wird a . Wie sieht diese Gleichung in gekrümmter Raumzeit aus? Beachten Sie nur, dass der partielle Ableitungsoperator nicht koordinatenunabhängig ist! Um das zu beheben, definieren wir den kovarianten Ableitungsoperator a die so funktioniert:
a u b = a u b + Γ c a b u c
Γ c a b wird das Christoffel-Symbol genannt und Sie können sehen, dass es eine Korrektur für den partiellen Ableitungsoperator bereitstellt, um ihn koordinatenunabhängig zu machen. Beachten Sie, dass in der flachen Raumzeit Γ c a b verschwindet und die kovariante Ableitung gleich der partiellen Ableitung ist. So, jetzt ist die Modifikation klar. Sie konvertieren einfach a zu a und die Erhaltungsgleichung in GR lautet:
a T a b = 0
Sie können diesen Prozess als Tensorisierung bezeichnen . Wann immer Sie eine Gleichung in flacher Raumzeit haben, wandeln Sie sie einfach in eine Tensorgleichung um, sodass sie sich ohne Krümmung auf Ihre flache Raumzeitgleichung reduziert. Es gibt einige zusätzliche Feinheiten in diesem Prozess. Es funktioniert nicht immer. (Siehe Kapitel 4 von Robert Walds Allgemeine Relativitätstheorie für eine Diskussion) Aber ich denke, diese Technik funktioniert für Ihre Gleichung, also wenden wir sie an:

ρ u ¨ = d ich v ( σ ) + f

Sehen wir uns diese Gleichung Term für Term an. Der erste Term beinhaltet eine zweite zeitliche Ableitung, die von der Definition einer Koordinatenzeit abhängt. Um es zu straffen, ändern Sie es einfach in die richtige Zeit. Im zweiten Term ist die Divergenz in flacher Raumzeit definiert als:

d ich v σ a b = a σ a b
In der gekrümmten Raumzeit wird es also einfach:
d ich v σ a b = a σ a b = 1 g a ( g σ a b )
wo g ist die Determinante des metrischen Tensors g a b .

Im dritten Semester f ist bereits ein Vektor, der koordinatenunabhängig ist. Also geht es dir gut. Somit lautet die endgültige Form:

ρ d 2 d τ 2 u b = a σ a b + f b

PS Als ich diese Antwort schrieb, wurde mir das klar ρ ist keine koordinatenunabhängige Größe. Die endgültige Form ist also wahrscheinlich falsch. Ich weiß nicht, wie die von Ihnen aufgestellte Gleichung hergeleitet wurde. Wenn Sie mir eine grundlegendere Gleichung geben, kann ich sie vielleicht korrigieren. Aber ich hoffe, die oben genannten Techniken können Ihnen helfen, es selbst herauszufinden!

Danke für die ausführliche Antwort, das wollte ich und mehr. Ich verstehe jedoch nicht den Schritt, die kovariante Ableitung des Tensors auszudrücken σ in gekrümmter Raumzeit in Bezug auf die Determinante der Metrik. Diese Determinante ist eine Invariante des Tensors g , warum ist es also wichtig, es im Ableitungsoperator zu belassen? Ich glaube auch, dass ich die Antwort kenne, aber für die anderen Gleichungen, die ich erwähnt habe (konstitutive Gleichung und Kompatibilität), muss ich wohl die gleiche Argumentation auf ihre räumlichen Ableitungen anwenden, und die Christoffel-Symbole werden die Metrik berücksichtigen.
Auch der Begriff ρ ist die Massendichte ( k g . m 3 ) und ist im klassischen Kontext generell positions- und zeitabhängig. Die lokale Bewegungsgleichung lautet also: ρ ( x , t ) u ¨ ( x , t ) = d ich v ( σ ( u ( x , t ) ) ) + f v ( x , t )
Ich habe versucht, die kovariante Ableitung des Tensors auszudrücken g a b in Bezug auf partielle Ableitungen (die NICHT invariant sind). Es sieht also etwas seltsam aus. Aber Sie können dieses Dokument nachschlagen: web.mit.edu/edbert/GR/gr1.pdf und herausfinden, wie der Ausdruck abgeleitet wurde.
Die Massendichte ist eine Komponente des Spannungsenergietensors (nämlich der T 00 Komponente) und ist somit koordinatenabhängig. Meine Ableitung war also wahrscheinlich falsch ... Können Sie mir eine Referenz geben, wie das lokale EoM abgeleitet wurde?
1/2. Ich kann derzeit keine einfache Referenz für die Ableitung des EoM finden, aber die Ableitung ist nicht allzu schwierig. Über den ganzen Körper ist Masse x Beschleunigung gleich der Summe der Kräfte. Diese bestehen aus Volumenkräften f v (Nm-3) und Oberflächenkräfte f s (Nm-2). Über den gesamten Bereich haben wir also: Ω ρ u ¨ d v = Ω f s d S + Ω f v d v .
2/2. Der Flächenkraftvektor heißt Traktionsvektor und es wird angenommen, dass er linear vom nach außen gerichteten Normalenvektor abhängt (Satz von Cauchy), also f s = σ . n wo σ heißt Cauchy-Spannungstensor. Übertragen Sie das Oberflächenintegral zurück in ein Volumenintegral, indem Sie eine partielle Integration verwenden, die Folgendes ergibt: Ω σ . n d S = Ω d ich v ( σ ) d v . Mit Ihrer gesamten Gleichung integriert über eine beliebige Ω , können Sie das Integral entfernen, um das lokale Gesetz zu erhalten, das ich ursprünglich angegeben habe, aktualisiert mit der x- und t-Abhängigkeit in den obigen Kommentaren.
Danke für die Klarstellung! Ich verstehe es nicht ganz, aber ich sehe, wo die σ Begriff kommt von jetzt. Übrigens, ρ könnte als Invariante interpretiert werden, wenn Sie es wie folgt definieren: ρ = T a b u a u b . Also könnte meine Umrechnung stimmen. Steils Verweis auf die Numerische Relativitätstheorie ist jedoch definitiv viel präziser und umfassender. Also lies das :)
Ja, ich muss mich mit vielen Konzepten vertraut machen, bevor ich etwas vollständig begreife. Vielen Dank für die tolle Hilfe. Ich werde eine Weile warten, um zu sehen, ob andere darauf antworten, dann werde ich dies zur akzeptierten Antwort machen.
Nehmen Sie sich Zeit, Mann :) Ich bin auch gespannt auf andere bessere Antworten!
Suchen Sie nach Artikeln und Büchern über allgemeine relativistische Hydrodynamik und kinetische Theorie. Es funktioniert gut für Kosmologie und relativistische Materie und "Flüssigkeiten" und "Gase", und es ist sehr gut definiert und ich bin mir sicher, dass es auch einige Codes gibt. Keine Notwendigkeit, neu zu erfinden.
@ZhengyanShi Ich glaube nicht, dass Ihre Idee, einen Newtonschen Gl. funktioniert oder ist in diesem Fall sogar notwendig. Wie Sie bereits erwähnt haben, die Massendichte T 00 , somit a T a b = 0 enthält zwangsläufig eine Gleichung für die Massendichte. Warum nicht rückwärts von der explizit kovarianten Form zu einer "Divergenzform" in Bezug auf reguläre partielle Ableitungen arbeiten, die dann diskretisiert werden können usw., so etwas wie a T a b + = 0 ? Ich bin mit numerischem GR nicht vertraut, aber ich denke, dies sollte zumindest ein Teil des Ausgangspunkts sein.
@udrv Du hast wahrscheinlich recht. Obwohl ich denke, dass die resultierende Gleichung dieselbe wäre, nachdem Sie die Gleichung erweitert haben a T a b = 0 und interpretieren Sie Komponenten des Spannungs-Energie-Tensors (ein Teil davon stellt beispielsweise etwas dar, das mit dem Cauchy-Spannungstensor verwandt ist, und ein Teil davon die Energiedichte).
@BrouH Um es klar auszudrücken, der Schritt, durch den Sie in Ihrem ersten Kommentar verwirrt waren (Umwandeln der Divergenz in einen Ausdruck mit "expliziten" Ableitungen der Determinante der Metrik), ist nicht trivial und verdient eine Rechtfertigung. Siehe Seite 34 des folgenden Dokuments. people.physics.tamu.edu/pope/grav-phys.pdf