Einstein Toolkit Documentation - User contributions [en]

MHD implementation details and discussions

2011-04-26T18:25:03Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber <strike>will commit</strike> has committed the routines mentioned in these notes <strike>soon</strike>.
<li> Scott Noble and Josh Faber worked on Con2PrimM and have committed.
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:<br>
[[Image:Rel_cleaning_new.pdf]]<br>
[[Image:Rel_cleaning_new.tex]]<br>

File:Rel cleaning new.tex

2011-04-26T18:24:36Z

Snoble: Simplified and corrected derivation of the general relativistic divergence cleaning equations in 3+1 form.

Simplified and corrected derivation of the general relativistic divergence cleaning equations in 3+1 form.

File:Rel cleaning new.pdf

2011-04-26T18:24:14Z

Snoble: Simplified and corrected derivation of the general relativistic divergence cleaning equations in 3+1 form.

Simplified and corrected derivation of the general relativistic divergence cleaning equations in 3+1 form.

MHD implementation details and discussions

2011-04-25T20:17:16Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber <strike>will commit</strike> has committed the routines mentioned in these notes <strike>soon</strike>.
<li> Scott Noble and Josh Faber worked on Con2PrimM and have committed.
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:<br>
[[Image:rel_cleaning.pdf]]<br>
[[Image:rel_cleaning.tex]]<br>

MHD implementation details and discussions

2011-04-25T20:15:34Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber <strike>will commit</strike> has committed the routines mentioned in these notes <strike>soon</strike>.
<li> Scott Noble is working Con2PrimM, will hopefully commit by the end of this week (or next week at least).
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:<br>
[[Image:rel_cleaning.pdf]]<br>
[[Image:rel_cleaning.tex]]<br>

MHD implementation details and discussions

2011-04-25T20:14:40Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber will commit the routines mentioned in these notes soon.
<li> Scott Noble is working Con2PrimM, will hopefully commit by the end of this week (or next week at least).
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:<br>
[[Image:rel_cleaning.pdf]]<br>
[[Image:rel_cleaning.tex]]<br>

MHD implementation details and discussions

2011-04-25T20:14:17Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber will commit the routines mentioned in these notes soon.
<li> Scott Noble is working Con2PrimM, will hopefully commit by the end of this week (or next week at least).
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:
[[Image:rel_cleaning.pdf]]
[[Image:rel_cleaning.tex]]

MHD implementation details and discussions

2011-04-25T20:09:40Z

Snoble:

<h2>Scientific Goals</h2>

<h2>Random bits and pieces</h2>
<ul>
<li><strike>Reuse parts of Harm3D</strike></li>
<li>Used Harm3d as a reference, but the only parts of Harm3d directly included were for the Con2PrimM routines. </li>
</ul>

<h2>Conventions / Variable names etc.</h2>
<ul>
<li> Principally follow the Valencia formulation as out-lined in [http://adsabs.harvard.edu/abs/2006ApJ...637..296A Anton et al. (2006)] </li>
<li> HydroBase's Bvec[1-3] GF will be identified with Valencia's B^i vector </li>
<li> Bvec[] will be used both as an evolved quantity and a primitive variable</li>
<li> During the Carpet workshop, we concluded it is easier to define B^i to be the undensitized magnetic field; this requires some care, since psi^6 B^i is the quantity that appears in the evolution equation, but simplifies the rest of the code
<li> Only use locally calculated "comoving magnetic tensor" b^a when needed </li>
<li> In general, but especially for MHD, it seems worthwhile to promote sqrt(gamma)=psi^6 to an independent gridfunction
<li> See Josh Faber's notes here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD.tex], or check out the code (still a work in progress!!!) here [http://ccrg.rit.edu/~carpet/wiki/File:GRHydro_MHD_routines.tgz].
<li> Josh Faber will commit the routines mentioned in these notes soon.
<li> Scott Noble is working Con2PrimM, will hopefully commit by the end of this week (or next week at least).
</ul>

<h2>Subprojects</h2>
<table border="1">
<tr><th>Task </th><th>Main people </th><th>Notes</th><th>Status</th></tr>
<tr><td>Riemann Problem </td><td>Faber </td><td></td><td>HLLE done</td></tr>
<tr><td>Con2Prim </td><td>Noble, Faber</td><td>from Harm3D</td><td>done, except for general EOS</td></tr>
<tr><td>Constraint transport </td><td>? </td><td>from Harm3D</td><td>still on drawing board</td></tr>
<tr><td>Reconstruction </td><td>Faber, Mundim</td><td></td><td>done</td></tr>
<tr><td>Failure Recovery Procedure </td><td>? </td><td></td><td>just using old tools for now</td></tr>
<tr><td>Carpet Refluxing </td><td>Erik Schnetter </td><td></td><td>Finished?</td></tr>
<tr><td>EMF Restriction/Prolongation</td><td>?</td><td>Needed for CT only</td><td></td></tr>
<tr><td>Constraint Damping?</td><td>Faber, Mundim, Noble</td><td>SR case done</td><td>GR derived, not implemented yet</td></tr>
</table>

During the Cactus workshop, some progress was made on implementing MHD within GRHydro by reimplementing various subroutines with alternate versions for MHD and calling them out with the MHD flag in the scheduler. Basically, instead of GRHydro_filename.F90, we will have GRHydro_filenameM.F90.

Most of these *M.F90 files have generated and are part of ET already.

Notes on the derivation of the GR divergence cleaning equations:
https://docs.einsteintoolkit.org/et-docs/images/6/68/rel_cleaning.pdf
https://docs.einsteintoolkit.org/et-docs/images/6/68/rel_cleaning.tex

File:rel cleaning.pdf

2011-04-25T20:08:09Z

Snoble: PDF version of notes.

PDF version of notes.

File:rel cleaning.tex

2011-04-25T20:07:37Z

Snoble: Latex file containing notes on the derivation of the general relativistic divergence cleaning equations in 3+1 equations.

Latex file containing notes on the derivation of the general relativistic divergence cleaning equations in 3+1 equations.