tc_flux.c File Reference

Compute the thermal conduction flux. More...

#include "pluto.h"

Macros
#define	HYPERBOLIC_SAT_FLUX   YES

Functions
void	TC_Flux (double ***T, const State_1D *state, double **dcoeff, int beg, int end, Grid *grid)

Detailed Description

Compute the thermal conduction flux along one row of computational zones for the HD and MHD modules according to Spitzer (1962):

where is the classical thermal conduction flux, is the saturated flux. Since the first term is purely parabolic, it is discretized using standard finite difference. The saturated flux is treated in an upwind manner following the guidelines given in Balsara (2008) (see also Mignone et al. 2012)

The classical MHD flux further splits into 2 components, along and across the magnetic field lines. This function also computes the inverse of the time step and return its maximum over the current sweep.

References

• "The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics"
  Mignone et al, ApJS (2012) 198, 7M

Authors

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

Date

Sep 13, 2012

Macro Definition Documentation

#define HYPERBOLIC_SAT_FLUX   YES

When set to YES, saturated flux is computed using an upwind selection rule. When set to NO, staurated flux is treated in the same manner as the conduction flux.

Function Documentation

void TC_Flux	(	double ***	T,
		const State_1D *	state,
		double **	dcoeff,
		int	beg,
		int	end,
		Grid *	grid
	)

Compute the thermal conduction flux, state->par_flx.

Parameters

[in]	T	3D array containing the dimensionless temperature
[in,out]	state	pointer to a State_1D structure
[out]	dcoeff	the diffusion coefficient needed for computing the time step.
[in]	beg	initial index of computation
[in]	end	final index of computation
[in]	grid	pointer to an array of Grid structures

Returns

This function has no return value.