hugomvale/fixedpivot2d

Solution of population balance equations for 2D aggregation by fixed pivot method

FixedPivot2D

Description

This code is a Fortran 90 implementation of the extended fixed pivot method developed by Vale & McKenna to solve population balance equations (PBE) for two-component aggregation processes (also known as bivariate aggregation).

The method is described in detail in Vale2005b.pdf.

How to use the code

1. File INPUTS.f90:

• Function COAGKERNEL
  • Define the desired aggregation/coagulation kernel.
• Function PSDINI
  • Define initial particle size distribution (PSD).

2. File MAIN.f90:

• Program FIXEDPIVOT2D
  • Variable time_final: define end integration time.
  • Variables xmin, xmax, ymin, ymax: define grid range.
  • Variables meshx, meshy: define grid type. See subroutine GRID in file GRID.f90 for more information.

3. File GLOBAL.f90

• Module GLOBAL
  • Variables Mx, My: define number of grid cells.

4. File OUTPUT.f90

• Subroutine OUTPUT
  • Update folder path.

5. Build and run

• Do not forget to include DLSODE or ODEPACK, e.g. libDLSODE.a