Difference between revisions of "MHD development"
(Add/reshuffle topics) |
|||
| Line 3: | Line 3: | ||
List of Potential Topics | List of Potential Topics | ||
* Paper updates and progress | * Paper updates and progress | ||
| − | * Divergence cleaning and mesh refinement | + | * Issues with current MHD methods |
| − | * Combining constraint transport and divergence cleaning | + | ** Divergence cleaning and mesh refinement |
| − | * Vector Potential implementation | + | ** Combining constraint transport and divergence cleaning |
| + | * New methods | ||
| + | ** Con2Prim Inversion -- Force-free | ||
| + | ** Vector Potential implementation | ||
| + | |||
| + | = Con2Prim Inverter = | ||
| + | |||
| + | * add code to transition into force-free inversion (no force free evolution) if magnetic pressure starts to dominate gas pressure (done in SpEC, seems robust) | ||
= Vector Potential = | = Vector Potential = | ||
| − | Evolving the vector potential allows for | + | Evolving the vector potential allows for better control of the B-field divergence and better behavior around a puncture. Evolving the vector potential is also better behaved surrounding puncture BHs. |
=== Current State === | === Current State === | ||
| − | This was partially implemented in the last development workshop, | + | This was partially implemented in the last development workshop, but removed again within the last year. We had agreed to implement the algebraic gauge first with a view towards the Lorenz gauge. Initial implementation had been of cell-centered vector potential, though we should now be able to support staggered grid vector potential. |
=== Matters to Discuss === | === Matters to Discuss === | ||
| Line 31: | Line 38: | ||
* Riemannsolver support for A^mu | * Riemannsolver support for A^mu | ||
* Recovery point-wise support for A^mu | * Recovery point-wise support for A^mu | ||
| − | |||
| − | |||
| − | |||
| − | |||
Revision as of 16:34, 23 October 2012
Contents
MHD Development Agenda
List of Potential Topics
- Paper updates and progress
- Issues with current MHD methods
- Divergence cleaning and mesh refinement
- Combining constraint transport and divergence cleaning
- New methods
- Con2Prim Inversion -- Force-free
- Vector Potential implementation
Con2Prim Inverter
- add code to transition into force-free inversion (no force free evolution) if magnetic pressure starts to dominate gas pressure (done in SpEC, seems robust)
Vector Potential
Evolving the vector potential allows for better control of the B-field divergence and better behavior around a puncture. Evolving the vector potential is also better behaved surrounding puncture BHs.
Current State
This was partially implemented in the last development workshop, but removed again within the last year. We had agreed to implement the algebraic gauge first with a view towards the Lorenz gauge. Initial implementation had been of cell-centered vector potential, though we should now be able to support staggered grid vector potential.
Matters to Discuss
- Staggered vs unstaggered A^i
- Prolongation operators given the location on the cell ( for staggered A only )
- Scheduling for reconstructing B from the vector potential
Coding Tasks
As this workshop will take place during the pre-release phase, this development should take place in Zelmani.
- Storage and evolution (simpler now that E^i has been implemented and tested for the constraint transport)
- B-field construction from A^mu
- Initial Prim2Con support for A^mu
- Reconstruction support for A^mu
- Riemannsolver support for A^mu
- Recovery point-wise support for A^mu
