Purpose

Solver for coupled CFD-DEM simulations involving incompressible fluids. Supports transport of temperature as a passive scalar.

Warning

GPU support for this model or solver has not been tested and may not work as expected.

cfdemSolverPisoScalar

Description

cfdemSolverPisoScalar is a coupled CFD-DEM solver using CFDEMcoupling, an open source parallel coupled CFD-DEM framework. Based on pisoFoam®(*), a finite volume based solver for turbulent Navier-Stokes equations applying PISO algorithm, cfdemSolverPisoScalar has additional functionality for a coupling to the DEM code “Aspherix” as well as a scalar transport equation. The volume averaged Navier-Stokes Equations are solved accounting for momentum exchange and volume displacement of discrete particles, whose trajectories are calculated in the DEM code Aspherix.

The scalar transport equation using the voidfraction \alpha, the thermal diffusivities \alpha_\mathrm{t} and \alpha_\mathrm{eff} and the volumetric flux \phi

\alpha_\mathrm{t} &= \frac{\nu_\mathrm{t}}{\mathit{Pr}_\mathrm{t}}\\
\alpha_\mathrm{eff} &= \frac{\nu}{\mathit{Pr}} + \alpha_\mathrm{t}

\frac{\partial \alpha T}{\partial t} &+ \nabla \cdot \left( \phi T \right)
    - \nabla \cdot \left( \alpha_\mathrm{eff} \alpha \nabla T \right)
    = S_\mathrm{particles} + T_\mathrm{source}

models heat transfer withing the fluid and, in addition, heat transfer between the fluid and particles by the additional source term S_\mathrm{particles}. Thus, convective heat transfer in a fluid granular system can be modeled with cfdemSolverPisoScalar.

The transport equation uses a field alphat to calculate local kinematic turbulent thermal conductivities based on the laminar and turbulent Prandtl numbers. Additional sources are treated via the source term TSource.

In addition to the coupling fields listed in cfdemSolverPiso, this solver requires the following additional input:

name

location

type

unit

default

T

0

scalar field

K

required

TSource

0

scalar field

K/s

0

alphat

0

scalar field

m²/s

0

Pr

transportProperties

scalar

1

required

Prt

transportProperties

scalar

1

required

Literature

GONIVA, C., KLOSS, C., HAGER, A. and PIRKER, S. (2010): “An Open Source CFD-DEM Perspective”, Proc. of OpenFOAM Workshop, Göteborg, June 22.-24.

The heat transfer equation is implemented according to Nield & Bejan (2013), Convection in Porous Media, DOI 10.1007/978-1-4614-5541-7_2, Springer


(*) This offering is not approved or endorsed by OpenCFD Limited, the producer of the OpenFOAM software and owner of the OPENFOAM® and OpenCFD® trade marks.