cfbb389fd3
This change adds representation of the shape of a dispersed phase. A
layer has been added to model the relationship between the
characteristic volume of a sizeGroup and its physical diameter.
Previously this relationship was represented by a constant form factor.
Currently, two shape models are available:
- spherical
- fractal (for modelling fractal agglomerates)
The latter introduces the average surface area to volume ratio, kappa,
of the entities in a size group as a secondary field-dependent internal
variable to the population balance equation, which makes the population
balance approach "quasi-"bivariate. From kappa and a constant mass
fractal dimension, a collisional diameter can be derived which affects
the coagulation rates computed by the following models:
- ballisticCollisions
- brownianCollisions
- DahnekeInterpolation
- turbulentShear
The fractal shape modelling also takes into account the effect of sintering
of primary particles on the surface area of the aggregate.
Further additions/changes:
- Time scale filtering for handling large drag and heat transfer
coefficients occurring for particles in the nanometre range
- Aerosol drag model based on Stokes drag with a Knudsen number based
correction (Cunningham correction)
- Reaction driven nucleation
- A complete redesign of the sizeDistribution functionObject
The functionality is demonstrated by a tutorial case simulating the
vapour phase synthesis of titania by titanium tetrachloride oxidation.
Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum Dresden -
Rossendorf (HZDR) and VTT Technical Research Centre of Finland Ltd.
40 lines
1.0 KiB
Bash
Executable File
40 lines
1.0 KiB
Bash
Executable File
#!/bin/sh
|
|
|
|
if ! which gnuplot > /dev/null 2>&1
|
|
then
|
|
echo 'gnuplot not found - skipping graph creation' >&2
|
|
exit 1
|
|
fi
|
|
|
|
gnuplot<<EOF
|
|
set terminal postscript eps monochrome
|
|
set output '../numberDensity.eps'
|
|
set decimalsign '.'
|
|
|
|
set format xy '%g'
|
|
set xtics 1e1 mirror
|
|
set xlabel 'v(m^{3})'
|
|
set ytics 1e1 mirror
|
|
set ylabel 'n(m^{-3}m^{-3})'
|
|
|
|
set logscale xy
|
|
set xrange [1e-5:1.5]
|
|
set yrange [1e-5:10]
|
|
set key at graph 0.8,0.5
|
|
|
|
N0 = 0.05
|
|
v0 = 1.0
|
|
|
|
n(x,t) = N0*exp(-t*x**2.0)*2.0*t*v0
|
|
|
|
numberDensity = '../postProcessing/numberDensity/10/numberDensity.dat'
|
|
|
|
plot '<printf "1 0.05\n1 0\n"' w l t 'Init. Cond.',\
|
|
n(x,10.0) ls 2 lc rgb 'black' t '[Ziff and McGrady, J. Phys. A: Math. Gen., 18: 3027-3037, 1985]',\
|
|
numberDensity every ::0::19 u 1:4 w p pt 1 t 'air1',\
|
|
numberDensity every ::20::24 u 1:4 w p pt 5 t 'air2',\
|
|
numberDensity every ::25::28 u 1:4 w p pt 9 t 'air3'
|
|
EOF
|
|
|
|
#------------------------------------------------------------------------------
|