fluxes.c File Reference

Compute the RMHD flux. More...

#include "pluto.h"

Include dependency graph for fluxes.c:

Go to the source code of this file.

Functions
void	Flux (double **ucons, double **uprim, double *h, double **fx, real *pr, int beg, int end)

Detailed Description

Compute the RMHD flux.

Compute the flux of the conservative RMHD equations in the direction given by g_dir. This function defines the component of the hyperbolic flux tensor of the standard MHD equations.
In what follows:

• VXn, MXn, BXn are the velocity, momentum and magnetic field components in the direction given by g_dir (normal, "n")
• VXt, MXt, BXt and VXb, MXb, BXb are the transverse components (tangent "t" and bi-tangent "b").

Author

A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)

Date

Jun 10, 2015

Definition in file fluxes.c.

Function Documentation

void Flux	(	double **	ucons,
		double **	uprim,
		double *	h,
		double **	fx,
		real *	pr,
		int	beg,
		int	end
	)

Parameters

[in]	u	1D array of conserved quantities
[in]	w	1D array of primitive quantities
[in]	a2	1D array of sound speeds
[out]	fx	1D array of fluxes (total pressure excluded)
[out]	p	1D array of pressure values
[in]	beg	initial index of computation
[in]	end	final index of computation

Returns

This function has no return value.

Definition at line 23 of file fluxes.c.

{
  int    nv, i;
  double vB, uB, b2, wt, wtg2, Bmag2, pt; 
  double *u, *v, b[4], g, g2, g_2;

  for (i = beg; i <= end; i++) {

    u = ucons[i];
    v = uprim[i];

    g     = u[RHO]/v[RHO];
    g2    = g*g;
    g_2   = 1.0/g2;

    Bmag2 = EXPAND(v[BX1]*v[BX1], + v[BX2]*v[BX2], + v[BX3]*v[BX3]);    

    vB    = EXPAND(v[VX1]*v[BX1], + v[VX2]*v[BX2], + v[VX3]*v[BX3]);
    EXPAND(b[IDIR] = g*(v[BX1]*g_2 + vB*v[VX1]); , 
           b[JDIR] = g*(v[BX2]*g_2 + vB*v[VX2]); ,
           b[KDIR] = g*(v[BX3]*g_2 + vB*v[VX3]);)
    b2 = Bmag2*g_2 + vB*vB;
   
    pt   = v[PRS] + 0.5*b2;
    wt   = v[RHO]*h[i] + b2;
    wtg2 = wt*g2;

    fx[i][RHO]  = u[RHO]*v[VXn];
    EXPAND(fx[i][MX1] = wtg2*v[VX1]*v[VXn] - b[IDIR]*b[g_dir];  ,
           fx[i][MX2] = wtg2*v[VX2]*v[VXn] - b[JDIR]*b[g_dir];  ,
           fx[i][MX3] = wtg2*v[VX3]*v[VXn] - b[KDIR]*b[g_dir];)
    EXPAND(fx[i][BXn] = 0.0;             ,
           fx[i][BXt] = v[VXn]*v[BXt] - v[BXn]*v[VXt];   ,
           fx[i][BXb] = v[VXn]*v[BXb] - v[BXn]*v[VXb]; )

    fx[i][ENG] = u[MXn];
  #if RMHD_REDUCED_ENERGY
    fx[i][ENG] -= fx[i][RHO];
  #endif

    pr[i] = pt;

#ifdef GLM_MHD
    fx[i][BXn]     = v[PSI_GLM];
    fx[i][PSI_GLM] = glm_ch*glm_ch*v[BXn];
#endif

  }
}