Difference between revisions of "MHD implementation details and discussions"

From Einstein Toolkit Documentation
Jump to: navigation, search
Line 20: Line 20:
 
<table border="1">
 
<table border="1">
 
  <tr><th>Task                        </th><th>Main people  </th><th>Notes</th></tr>
 
  <tr><th>Task                        </th><th>Main people  </th><th>Notes</th></tr>
  <tr><td>Riemann Problem            </td><td>             </td><td></td></tr>
+
  <tr><td>Riemann Problem            </td><td>Scott Noble  </td><td></td></tr>
 
  <tr><td>Con2Prim                    </td><td>Christian Ott</td><td>from Harm3D</td></tr>
 
  <tr><td>Con2Prim                    </td><td>Christian Ott</td><td>from Harm3D</td></tr>
  <tr><td>Constraint transport        </td><td>             </td><td>from Harm3D</td></tr>
+
  <tr><td>Constraint transport        </td><td>Scott Noble  </td><td>from Harm3D</td></tr>
 
  <tr><td>Reconstruction              </td><td>Frank Löffler</td><td></td></tr>
 
  <tr><td>Reconstruction              </td><td>Frank Löffler</td><td></td></tr>
  <tr><td>Failure Recovery Procedure  </td><td>             </td><td></td></tr>
+
  <tr><td>Failure Recovery Procedure  </td><td>Scott Noble  </td><td></td></tr>
 
  <tr><td>Carpet Refluxing            </td><td>            </td><td></td></tr>
 
  <tr><td>Carpet Refluxing            </td><td>            </td><td></td></tr>
 
  <tr><td>EMF Restriction/Prolongation</td><td>            </td><td></td></tr>
 
  <tr><td>EMF Restriction/Prolongation</td><td>            </td><td></td></tr>
 
</table>
 
</table>

Revision as of 15:53, 14 June 2010

Scientific Goals

Random bits and pieces

  • Reuse parts of Harm3D

Conventions / Variable names ect.

  • Principally follow the Valencia formulation as out-lined in Anton et al. (2006)
  • HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector
  • Bvec[] will be used both as an evolved quantity and a primitive variable
  • Introduce "densitized_state" as a scalar that controls interpretation of Bvec[]'s meaning;
  • if( densitized_state) then Bvec[] = B^i \sqrt(\gamma) else Bvec[] = B^i
  • Only use locally calculated "comoving magnetic tensor" b^a when needed


Subprojects

Task Main people Notes
Riemann Problem Scott Noble
Con2Prim Christian Ottfrom Harm3D
Constraint transport Scott Noble from Harm3D
Reconstruction Frank Löffler
Failure Recovery Procedure Scott Noble
Carpet Refluxing
EMF Restriction/Prolongation