twoFluidEMHD.h

Go to the documentation of this file.

#ifndef TWOFLUIDEMHD_H
#define TWOFLUIDEMHD_H

#include "model.h"

// Two-Fluid ElectroMagnetoHydroDynamics
/*  Human readable description
                                                            <br>
  The model has eighteen conserved variables:                                                         <br>
(0)   D, Sx, Sy, Sz, tau,                                                                             <br>
(5)   Dbar, Sbarx, Sbary, Sbarz, taubar,                                                              <br>
(10)  Bx, By, Bz,                                                                                     <br>
(13)  Ex, Ey, Ez,                                                                                     <br>
(16)  psi, phi                                                                                        <br>
  Sixteen primitive variables:                                                                        <br>
(0)   rho1, vx1, vy1, vz1, p1,                                                                        <br>
(5)   rho2, vx2, vy2, vz2, p2,                                                                        <br>
(10)  Bx, By, Bz,                                                                                     <br>
(13)  Ex, Ey, Ez                                                                                      <br>
  Thirty-five auxiliary variables:                                                                   <br>
(0)   h1, W1, e1, vsq1, Z1, D1, Stildex1, Stildey1, Stildez1, tauTilde1,                              <br>
(10)  h2, W2, e2, vsq2, Z2, D2, Stildex2, Stildey2, Stildez2, tauTilde2,                              <br>
(20)  Jx, Jy, Jz,                                                                                     <br>
(23)  Stildex, Stildey, Stildez, tauTilde,                                                            <br>
(27)  Bsq, Esq                                                                                        <br>
(29)  rhoCh0, rhoCh,                                                                                  <br>
(31)  ux, uy, uz, W                                                                                   <br>
*/




class TwoFluidEMHD : public Model
{

  public:

    TwoFluidEMHD();     


    TwoFluidEMHD(Data * data);

    virtual ~TwoFluidEMHD() { }     



    void sourceTermSingleCell(double *cons, double *prims, double *aux, double *source, int i=-1, int j=-1, int k=-1);


    void sourceTerm(double *cons, double *prims, double *aux, double *source);


    void getPrimitiveVarsSingleCell(double *cons, double *prims, double *aux, int i=-1, int j=-1, int k=-1);


    void getPrimitiveVars(double *cons, double *prims, double *aux);


    void primsToAll(double *cons, double *prims, double *aux);


    void fluxVector(double *cons, double *prims, double *aux, double *f, const int dir);


    void finalise(double *cons, double *prims, double *aux) { };
};


#endif