Energy Equation

Energy equation in nanoFluidX is implemented so that it accommodates for conduction and convection heat transfer with initial or Dirichlet boundary conditions.

The rate of temperature change of a specific material is given by:

\frac{d T}{d t} = \frac{1}{c_{p} ρ} \nabla (κ \nabla T)

Where,

$t$: Time
$T$: Temperature
$k$: Coefficient of thermal conductivity
$ρ$: Density
$c_{p}$: Specific heat capacity of the material

In SPH form the above analytic expression becomes:

{(\frac{d T}{d t})}_{i} = \frac{1}{c_{p, i}} \sum_{j} \frac{4 m_{j}}{ρ_{i} ρ_{j}} \frac{κ_{i} κ_{j}}{κ_{i} + κ_{j}} T_{i j} \frac{\nabla W_{i j}}{| r_{i j} |} \cdot n

Where,

$i$ and $j$ indices: Stand for owner and neighboring particles respectively
$ij$ index: Difference between the respective variables of particle $i$ and particle $j$
$m$: Mass of particle
$\nabla W$: Gradient of the kernel
$r$: Position vector
$n$: Unit coordinate vector

If you want to use the energy transport option, you need to turn the feature on by setting the energy_transport flag in Simulation parameters to true. Once this step is completed, you can opt for desired output, which can be either temperature or rate of temperature change (flux). The code will not work if an output is required and the energy transport flag is turned off.

Finally, set up the phase parameters, which are initial temperature (temp_0) in [K], evolve temperature flag (evolve_temp), specific heat capacity of the material (heat_cap) in [J/kg/K] and heat conduction coefficient of the material (heat_cond) in [W/m/K].

Note: Regarding the evolve_temp flag, this parameter is valid only for WALL or MOVINGWALL phases. If the evolve_temp flag is set to false, or if it is not defined (default is false), the set initial temperature will remain constant throughout the simulation. If set to true, the initial temperature will evolve in time and the temperature of the WALL or MOVINGWALL will be influenced by any surrounding phase. For FLUID phases, the temperature always evolves.