tc_kappa.c File Reference

#include "pluto.h"

Include dependency graph for tc_kappa.c:

Go to the source code of this file.

Functions
void	TC_kappa (double *v, double x1, double x2, double x3, double *kpar, double *knor, double *phi)

Function Documentation

void TC_kappa	(	double *	v,
		double	x1,
		double	x2,
		double	x3,
		double *	kpar,
		double *	knor,
		double *	phi
	)

Compute thermal conduction coefficients.

Parameters

[in]	v	array of primitive variables
[in]	x1	coordinate in the X1 direction
[in]	x2	coordinate in the X2 direction
[in]	x3	coordinate in the X3 direction
[out]	kpar	pointer to the conduction coefficient in the direction of magnetic field
[out]	knor	pointer to the conduction coefficient perpendicular to magnetic field
[out]	phi	pointer to the parameter controlling the magnitude of the saturated flux.

Definition at line 3 of file tc_kappa.c.

{
  double mu, nH, sqT, T, B2_cgs;
  double ctts_mu = 1.265060;

  mu   = ctts_mu/2.0;
  T    = v[PRS]/v[RHO]*(mu*CONST_mp*UNIT_VELOCITY*UNIT_VELOCITY/CONST_kB);
  sqT  = sqrt(T);

  *kpar = 5.6e-7*T*T*sqT;

  #if PHYSICS == MHD
   nH      = v[RHO]*UNIT_DENSITY/(mu*CONST_mp);
   B2_cgs  = EXPAND(v[BX1]*v[BX1], + v[BX2]*v[BX2], + v[BX3]*v[BX3]) + 1.e-12;
   B2_cgs *= 4.0*CONST_PI*UNIT_DENSITY*UNIT_VELOCITY*UNIT_VELOCITY;
   *knor   = 3.3e-16*nH*nH/sqT/B2_cgs;
  #else
   *knor = 0.0;
  #endif

  *kpar *= CONST_mp*mu/(UNIT_DENSITY*UNIT_VELOCITY*UNIT_LENGTH*CONST_kB);
  *knor *= CONST_mp*mu/(UNIT_DENSITY*UNIT_VELOCITY*UNIT_LENGTH*CONST_kB);

  *phi = 0.3;
}