Diffusion Processes

In general there are two types of diffusion processes: one for hydro, such as viscosity and thermal conduction; and one for the field, such as Ohmic dissipation and Ambipolar diffusion. We have included these four diffusion processes in current implementation.

Configuration

We do not need to manually switch on the diffusion in the configuration stage. They are determined by run-time parameters.

However some diffusion processes depend on specific physics and will be ignored if the code is not configured that way. For example, in order to calculate Ohmic dissipation and Ambipolar diffusion, magnetic field has to be enabled during configuration with -b. Similarly, thermal conduction is calculated when non-barotropic equation of state is enabled.

Input File

The additional microphysics are enrolled in automatically through the following input parameters:

<problem>
nuiso     = 0.01       # isotropic viscosity coefficient
nuani     = 0.0        # anisotropic viscosity coefficient
kiso      = 0.01       # isotropic thermal conduction coefficient
kani      = 0.0        # anisotropic thermal conduction coefficient
coef_o    = 0.01       # Ohmic resistivity coefficient
coef_a    = 0.0        # Ambipolar diffusion coefficient

Positive coefficients invoke the calculation of the corresponding physical processes.

Algorithms

Currently all diffusion processes are calculated with the conservative form, the diffusion fluxes are directly added to the hydro and field flux in a non-operator-split way.