Changes between Version 63 and Version 64 of SelfGravityDevel


Ignore:
Timestamp:
10/22/19 11:18:56 (5 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • SelfGravityDevel

    v63 v64  
    11= Self Gravity =
    2 
    3 (This is a skeleton while a new version of this page is under construction.)
    42
    53== Physics ==
    64=== Poisson Equation of Self Gravity ===
    75
    8   The electric field generated by a point charge[[latex($Q$)]]is
     6The inclusion of self gravity involves an additional force to the momentum equation
     7
     8[[latex($\frac{d \rho \mathbf{v}}{dt} = \mathbf{f}_g = -\rho \nabla \phi$)]]
     9
     10and energy equation
     11
     12[[latex($\frac{d E}{dt} = \mathbf{v} \cdot \mathbf{f}_g = -\rho \mathbf{v} \cdot \nabla \phi$)]]
     13
     14and requires the solution for the gas potential
     15
     16[[latex($\nabla^2 \phi = 4 \pi G \left (\rho - \rho_0\right)$)]]
     17
     18where
     19
     20* [[latex($\rho_0 = \bar{\rho}$)]] for periodic boundary conditions
     21* [[latex($\rho_0 = 0$)]] otherwise.
     22
     23
     24AstroBEAR uses hypre's linear solver to solve Poisson's equation for the potential.  However, there are two methods for including the source terms for momentum and energy.  The non-conservative approach simply includes the terms during a source update, while the conservative approach recasts the RHS of the momentum and energy equations as total divergences - so they can be differenced conservatively.
     25
     26
     27
     28
     29
     30The electric field generated by a point charge[[latex($Q$)]]is
    931
    1032 [[latex($E=\frac{1}{4\pi\epsilon_{0}}\frac{Q}{r^{2}}$)]]