eos.c File Reference

Ideal (constant gamma) EOS. 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

Ideal (constant gamma) EOS.

Implements essential functions for the constant-gamma law for hydro/MHD, relativistic hydro/relativistic 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 48 of file eos.c.

{
  int i;
  double gmmr, theta;

  gmmr = g_gamma/(g_gamma - 1.0);

/* ---------------------------------------------------------------
              Classical equations of state
   --------------------------------------------------------------- */

  #if PHYSICS == HD || PHYSICS == MHD
   for (i = beg; i <= end; i++) h[i] = gmmr*v[i][PRS]/v[i][RHO];
  #elif PHYSICS == RHD || PHYSICS == RMHD
   for (i = beg; i <= end; i++) {
     theta = v[i][PRS]/v[i][RHO];
     h[i] = 1.0 + gmmr*theta;
   }
  #endif
}

Here is the caller graph for this function:

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 80 of file eos.c.

{
  int i;
  double rho, th;

  #if PHYSICS == HD || PHYSICS == MHD
   for (i = beg; i <= end; i++){
     rho  = v[i][RHO];
     s[i] = v[i][PRS]/pow(rho,g_gamma);
   }
  #elif PHYSICS == RHD || PHYSICS == RMHD
   for (i = beg; i <= end; i++) {
     rho = v[i][RHO];
     s[i] = v[i][PRS]/pow(rho,g_gamma);
   }
  #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 16 of file eos.c.

{
  int  i;

  #if PHYSICS == HD || PHYSICS == MHD
   for (i = beg; i <= end; i++) cs2[i] = g_gamma*v[i][PRS]/v[i][RHO];
  #elif PHYSICS == RHD || PHYSICS == RMHD
   double theta;
   Enthalpy(v, h, beg, end);
   for (i = beg; i <= end; i++) {
     theta = v[i][PRS]/v[i][RHO];
     cs2[i] = g_gamma*theta/h[i];
   }
  #endif
}

Here is the caller graph for this function: