Difference between revisions of "MHD development"
| Line 1: | Line 1: | ||
| − | = MHD | + | = MHD Discussions/Agenda = |
| − | + | * ''Pre-Release Points: Need a publishable version of MHD used for the paper'' | |
| − | * ''Pre-Release Points: Need a publishable version of MHD'' | ||
** Need some addressing of Con2Prim for GRMHD pushed | ** Need some addressing of Con2Prim for GRMHD pushed | ||
** Test cases from MHD paper examples (~0.5-1.5 minutes tops) | ** Test cases from MHD paper examples (~0.5-1.5 minutes tops) | ||
| − | * | + | ** What place can we store SmoothedPuncture (Turduckening) function? |
| + | * ''Paper updates and progress'' | ||
* Issues with current MHD methods | * Issues with current MHD methods | ||
** Con2Prim failures (GRHydro/GRMHD) | ** Con2Prim failures (GRHydro/GRMHD) | ||
| Line 17: | Line 17: | ||
** Vector Potential implementation | ** Vector Potential implementation | ||
| − | = Con2Prim Inverter = | + | == 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) | * 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 better control of the B-field divergence and better behavior around a puncture. Evolving the vector potential is also better behaved surrounding puncture BHs. | 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. | ||
| Line 33: | Line 33: | ||
* A^\mu should live in HydroBase so initial data doesn't have to depend on GRHydro | * A^\mu should live in HydroBase so initial data doesn't have to depend on GRHydro | ||
* Want to apply proper staggered A^i so need to extend appropriate prolongation operators | * Want to apply proper staggered A^i so need to extend appropriate prolongation operators | ||
| + | * Start with this at heart once release version of GRMHD ready | ||
=== Coding Tasks === | === Coding Tasks === | ||
Revision as of 19:51, 23 October 2012
Contents
MHD Discussions/Agenda
- Pre-Release Points: Need a publishable version of MHD used for the paper
- Need some addressing of Con2Prim for GRMHD pushed
- Test cases from MHD paper examples (~0.5-1.5 minutes tops)
- What place can we store SmoothedPuncture (Turduckening) function?
- Paper updates and progress
- Issues with current MHD methods
- Con2Prim failures (GRHydro/GRMHD)
- Save initial state in all Con2Prim routines, check results before overwriting
- Pass ideal & acceptable tolerances, have end C2P tolerance GF
- Hacks: if failures within horizon, just use initial prims and run prim2con
- Divergence cleaning and mesh refinement
- Combining constraint transport and divergence cleaning
- Con2Prim failures (GRHydro/GRMHD)
- 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.
Discussion
- A^\mu should live in HydroBase so initial data doesn't have to depend on GRHydro
- Want to apply proper staggered A^i so need to extend appropriate prolongation operators
- Start with this at heart once release version of GRMHD ready
Coding Tasks
As this workshop will take place during the pre-release phase, this development should take place in Zelmani.
- Storage and interface (done)
- Evolution equations
- Bvec construction from A^mu
- Prim2Con/Con2Prim support for A^mu
- Reconstruction support for A^mu
- Riemannsolver support for A^mu
- Recovery point-wise support for A^mu
