eos.c File Reference

Taub-Matthews (TM) EOS for relativistic hydro and MHD. More...

#include "pluto.h"

Include dependency graph for eos.c:

Go to the source code of this file.

Functions
void	SoundSpeed2 (double **v, double *cs2, double *h, int beg, int end, int pos, Grid *grid)
void	Enthalpy (double **v, real *h, int beg, int end)
void	Entropy (double **v, double *s, int beg, int end)

Detailed Description

Taub-Matthews (TM) EOS for relativistic hydro and MHD.

Author

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

Date

April 14, 2014

Definition in file eos.c.

Function Documentation

void Enthalpy	(	double **	v,
		real *	h,
		int	beg,
		int	end
	)

Compute the enthalpy.

Parameters

[in]	v	1D array of primitive quantities
[in]	h	1D array of enthalpy values
[in]	beg	initial index of computation
[in]	end	final index of computation

Returns

This function has no return value.

Definition at line 50 of file eos.c.

{
  int i;
  double gmmr, theta;

  #if PHYSICS == RHD || PHYSICS == RMHD
   for (i = beg; i <= end; i++) {
     theta = v[i][PRS]/v[i][RHO];
     #if EOS == IDEAL
      h[i] = 1.0 + gmmr*theta;
     #elif EOS == TAUB
      h[i] = 2.5*theta + sqrt(2.25*theta*theta + 1.0);
    #endif
   }
  #endif
}

void Entropy	(	double **	v,
		double *	s,
		int	beg,
		int	end
	)

Compute the entropy.

Parameters

[in]	v	1D array of primitive quantities
[in]	s	1D array of entropy values
[in]	is	initial index of computation
[in]	ie	final index of computation

Returns

This function has no return value.

Definition at line 78 of file eos.c.

{
  int i;
  double rho, th;

  #if PHYSICS == RHD || PHYSICS == RMHD
   for (i = beg; i <= end; i++) {
     rho = v[i][RHO];
     #if EOS == IDEAL
      s[i] = v[i][PRS]/pow(rho,g_gamma);
     #elif EOS == TAUB
      th   = v[i][PRS]/rho;
      s[i] = v[i][PRS]/pow(rho,5./3.)*(1.5*th + sqrt(2.25*th*th + 1.0));
     #endif
   }
  #endif
}

void SoundSpeed2	(	double **	v,
		double *	cs2,
		double *	h,
		int	beg,
		int	end,
		int	pos,
		Grid *	grid
	)

Define the square of the sound speed.

Parameters

[in]	v	1D array of primitive quantities
[out]	cs2	1D array containing the square of the sound speed
[in]	h	1D array of enthalpy values
[in]	beg	initial index of computation
[in]	end	final index of computation
[in]	pos	an integer specifying the spatial position inside the cell (only for spatially-dependent EOS)
[in]	grid	pointer to an array of Grid structures

Returns

This function has no return value.

Definition at line 13 of file eos.c.

{
  int  i;
  double theta;

  #if PHYSICS == RHD || PHYSICS == RMHD
   Enthalpy(v, h, beg, end);
   for (i = beg; i <= end; i++) {
     theta = v[i][PRS]/v[i][RHO];
     #if EOS == IDEAL
      cs2[i] = g_gamma*theta/h[i];
     #elif EOS == TAUB
      cs2[i] = theta/(3.0*h[i])*(5.0*h[i] - 8.0*theta)/(h[i] - theta);
     #endif
   }
  #else
   print ("! SoundSpeed2: Taub EOS not defined for this physics module.\n");
   QUIT_PLUTO(1);
  #endif
}

Here is the call graph for this function: