MHD development - Revision history

Tbode at 21:33, 24 October 2012

2012-10-24T21:33:59Z

2012-10-24T21:19:08Z

2012-10-24T16:23:59Z

2012-10-24T16:06:31Z

2012-10-23T21:21:02Z

2012-10-23T20:11:35Z

2012-10-23T19:51:45Z

2012-10-23T19:39:08Z

2012-10-23T19:30:52Z

2012-10-23T19:23:33Z

@@ Line 32: / Line 32: @@
 References: [http://arxiv.org/abs/1110.4633 arXiv:1110.4633] for Lorenz gauge and [http://arxiv.org/abs/1007.2848 arXiv:1007.2848] for Algebraic gauge
 === Coding Tasks ===
@@ Line 41: / Line 45: @@
 * Riemannsolver -- Support for A^mu (Lorenz gauge) or verified no-flux for algebraic gauge
 * New prolongation operator for A_i

@@ Line 27: / Line 27: @@
 === 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.
+This was partially implemented in the last development workshop, but removed again within the last year. Initial implementation had been of cell-centered vector potential, though we should now be able to support staggered grid vector potential.
 === Coding Tasks ===
 As this workshop will take place during the pre-release phase, this development should take place in Zelmani.
-* Storage and interface (done)
+* Source terms for A -- Use B^i from appropriate locations for A_i to be located as expected
-* Source terms for A (done)
+* Conversion from (Aphi,A_i) -> Bvec
-* Convert (Aphi,A_i) -> Bvec, called with (or by) Con2Prim
+* Riemannsolver -- Support for A^mu (Lorenz gauge) or verified no-flux for algebraic gauge
 * New prolongation operator for A_i

@@ Line 34: / Line 34: @@
 * Storage and interface (done)
-* Source terms for A (just need to add source terms for A_i, Aphi)
+* Source terms for A (done)
 * Convert (Aphi,A_i) -> Bvec, called with (or by) Con2Prim
 * Riemannsolver support for A^mu (Lorenz gauge)

@@ Line 28: / Line 28: @@
 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 ===
@@ Line 35: / Line 46: @@
 * 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

@@ Line 41: / Line 41: @@
 * Storage and interface (done)
 * Evolution equations for Avec,Aphi
 * Bvec construction from A^mu (into HydroBase in the end, scheduled more often by GRHydro when Avec evolved)
@@ Line 46: / Line 48: @@
 * VecPot version of Con2Prim avoid Avec itself
 * Reconstruct A^mu

@@ Line 1: / Line 1: @@
-= MHD Development Agenda =
+= MHD Discussions/Agenda =
-Discussions
+* ''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)
-* '''Paper updates and progress'''
+** What place can we store SmoothedPuncture (Turduckening) function?
 * Issues with current MHD methods
 ** Con2Prim failures (GRHydro/GRMHD)
@@ Line 17: / Line 17: @@
 ** 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)
-= 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.
@@ Line 33: / Line 33: @@
 * 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
 === Coding Tasks ===

@@ Line 10: / Line 10: @@
 *** Save initial state in all Con2Prim routines, check results before overwriting
 *** Pass ideal & acceptable tolerances, have end C2P tolerance GF
 ** Divergence cleaning and mesh refinement
 ** Combining constraint transport and divergence cleaning