fluxes.c File Reference

Compute the hydro (HD) flux. More...

#include "pluto.h"

Include dependency graph for fluxes.c:

Go to the source code of this file.

Functions
void	Flux (double **u, double **w, double *a2, double **fx, double *p, int beg, int end)

Detailed Description

Compute the hydro (HD) flux.

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

• VXn, MXn are the velocity, momentum components in the direction given by g_dir (normal, "n")
• VXt, MXt and VXb, MXb 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

Aug 16, 2012

Definition in file fluxes.c.

Function Documentation

void Flux	(	double **	u,
		double **	w,
		double *	a2,
		double **	fx,
		double *	p,
		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, ii;

  for (ii = beg; ii <= end; ii++) {
    fx[ii][RHO] = u[ii][MXn];
    EXPAND(fx[ii][MX1] = u[ii][MX1]*w[ii][VXn]; ,
           fx[ii][MX2] = u[ii][MX2]*w[ii][VXn]; ,
           fx[ii][MX3] = u[ii][MX3]*w[ii][VXn];)
#if HAVE_ENERGY
    p[ii] = w[ii][PRS];
    fx[ii][ENG] = (u[ii][ENG] + w[ii][PRS])*w[ii][VXn];
#elif EOS == ISOTHERMAL
    p[ii] = a2[ii]*w[ii][RHO];
#endif
/*
#if DUST == YES
    fx[ii][RHO_D] = u[ii][MXn_D];
    EXPAND(fx[ii][MX1_D] = u[ii][MX1_D]*w[ii][VXn_D]; ,
           fx[ii][MX2_D] = u[ii][MX2_D]*w[ii][VXn_D]; ,
           fx[ii][MX3_D] = u[ii][MX3_D]*w[ii][VXn_D];)
#endif
*/
  }
}

Here is the caller graph for this function: