You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 18 Next »

Unable to render {include} The included page could not be found.

Pre-Analysis & Start-Up

This site is currently under construction. Please come back after it is fully built. Thank you!

In the Pre-Analysis & Start-Up step, we'll review the following:

  • Theory for Fluid Phase
  • Theory for Particle Phase
  • Choosing the Cases

Pre-Analysis:

A particle laden flow is a multiphase flow where one phase is the fluid and the other is dispersed particles. Governing equations for both phases are implemented in Fluent. To run a meaningful simulation, a review of the theory is necessary.

Fluid Phase:

In the simulations considered for this tutorial, the fluid flow is a 2D perturbed periodic double shear layer as described in the first section. The geometry is Lx = 59.15m, Ly = 59.15m, and the mesh size is chosen as

Unknown macro: {latex} $\Delta x = L_x / n_x$

in order to resolve the smallest vorticies. As a rule of thumb. One typically needs about 20 grid points across the shear layers, where the vorticies are going to develop. The boundary conditions are periodic in the x and y  directions. The fluid phase satisfies the Navier-Stokes Equations:
-Momentum Equations

Unknown macro: {latex}

\begin

Unknown macro: {eqnarray*}

 
\rho_f (\frac{d \textbf

Unknown macro: {u}

_f}

Unknown macro: {dt}

+\textbf

_f \cdot \nabla \textbf

Unknown macro: {u} _f)=- \nabla p + \mu \nabla ^2 \textbf

_f + \textbf

Unknown macro: {f}

 
\end

 

-Continuity Equation

Unknown macro: {latex}

 
\begin

Unknown macro: {align*}

 
\frac

Unknown macro: {partial rho_f}
Unknown macro: {partial t}

+ \nabla \cdot (\rho \textbf

Unknown macro: {u}

_f)=0
\end

 

where

Unknown macro: {latex}

$\textbf

Unknown macro: {u}

$

is the fluid velocity,

Unknown macro: {latex} $p$

the pressure,

Unknown macro: {latex} $\rho_f$

the fluid density and

Unknown macro: {latex}

$\textbf

Unknown macro: {f}

$

is a momentum exchange term due to the presence of particles. When the particle volume fraction and the particle mass loading

Unknown macro: {latex} $M=\phi \rho_p/\rho_f$

are very small, it is legitimate to neglect the effects of the particles on the fluid:

Unknown macro: {latex}

$\textbf

Unknown macro: {f}

$

can be set to zero. This type of coupling is called one-way. In these simulations the fluid phase is air, while the dispersed phase is constituted of about 400 glass beads of diameter a few dozens of micron. This satises both conditions

Unknown macro: {latex} $\phi \ll 1$

and

Unknown macro: {latex} $M \ll 1$

One way-coupling is legitimate here. See ANSYS documentation (16.2) for further details about the momentum exchange term.

Particle Phase:

 

The suspended particles are considered as rigid spheres of same diameter d, and density

Unknown macro: {latex} $\rho_p$

. Newton's second law written for the particle i stipulates:

Unknown macro: {latex}

$m_p \frac{d \textbf

Unknown macro: {u}

_p^i}

Unknown macro: {dt}

$

where uip is the velocity of particle i, Fiex the forces exerted on it, and mp its mass.
In order to know accurately the hydrodynamic forces exerted on a particle one needs to resolve the flow to a scale significantly smaller than the particle diameter. This is computationnaly pro- hibitive. Instead, the hydrodynamic forces can be approximated roughly to be proportional to the drift velocity:
d u ip = v − u ip dt τp
where τp = XX is known as the particle response time. This equation needs to be solved for all particles present in the domain. This is done in Fluent via the module: Discrete phase ?.

 

Expected Results

Under Construction


Go to Step 2: Geometry

Go to all FLUENT Learning Modules


  • No labels