Program combining non-Newtonian pipe flow and aggregation-orientation (AO) model. A fluid flowing
in a pipe is thixotropic -- in this case nanocellulose suspension. Particles in suspension are
rigid and ellipsoidal. In order to find out rheological properties of suspension, behavior of
particles must be calculated first. To do this, program is using particle distributions instead of
single particles in order to calculate distribution function for particle aggregations and
orientations. From AO distribution, local stress of fluid is calculated. Basically, in order to get
stress, program is calculating
where the fourth moment of orientation tensor is calculated from
See e.g. paper Rheological modeling of carbon nanotube aggregate suspensions
(W. K. A. Ma et al. 2008)
for more detailed description of AO model.
As a local stress of fluid is calculated, it is used to calculate a local viscosity from which
fluid velocity is calculated in discretization points in a following way
where λn are roots of Bessel function of the first kind of the first order, i.e.
J0(λn)=0.
Detailed description of source code functions (those written by me and few others) can be found
under pipeaggrot/doc/ directory.
Program requires C++11 and any compiler which supports it. Currently in the makefile Intel C++
compiler is used.
Also Boost C++ Libraries and SUNDIALS/CVODE ODE solver
library are required. Installation
instructions of those libraries can be found in their documentation. However, both libraries must
be installed under pipeaggrot/cppsrc/
directory. So, if using cmake in SUNDIALS installation, use the following commands
cmake -DCMAKE_INSTALL_PREFIX=$SOME_PATH/pipeaggrot/cppsrc \
-DEXAMPLES_INSTALL_PATH=$SOME_PATH/pipeaggrot/cppsrc/examples \
-DOPENMP_ENABLE=ON \
../sundials-VERSION