tc_kappa.c File Reference

Define the thermal conduction coefficients. More...

#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)

Detailed Description

Define the thermal conduction coefficients.

Use this function to supply the thermal conduction coefficients and along and across magnetic field lines and the parameter used to control the magnitude of the saturated flux . To exclude saturation, simply set to a very large number.

Authors

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

Date

Oct 31, 2012

Definition in file tc_kappa.c.

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 20 of file tc_kappa.c.

{
  double mu, T, B2_cgs, nH;

  mu   = 0.5;
  T    = v[PRS]/v[RHO]*mu*KELVIN;

  *kpar = 5.6e-7*T*T*sqrt(T);
  #if PHYSICS == MHD
   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;
   nH    = v[RHO]*UNIT_DENSITY/CONST_mp;
   *knor = 3.3e-16*nH*nH/(sqrt(T)*B2_cgs);
  #endif
  
  *phi = 0.3;
}