8 | | The electric field generated by a point charge[[latex($Q$)]]is |
| 6 | The 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 | |
| 10 | and energy equation |
| 11 | |
| 12 | [[latex($\frac{d E}{dt} = \mathbf{v} \cdot \mathbf{f}_g = -\rho \mathbf{v} \cdot \nabla \phi$)]] |
| 13 | |
| 14 | and requires the solution for the gas potential |
| 15 | |
| 16 | [[latex($\nabla^2 \phi = 4 \pi G \left (\rho - \rho_0\right)$)]] |
| 17 | |
| 18 | where |
| 19 | |
| 20 | * [[latex($\rho_0 = \bar{\rho}$)]] for periodic boundary conditions |
| 21 | * [[latex($\rho_0 = 0$)]] otherwise. |
| 22 | |
| 23 | |
| 24 | AstroBEAR 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 | |
| 30 | The electric field generated by a point charge[[latex($Q$)]]is |