| 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 |