## MATLAB code for discretization of regular sphere packings

Here you can find the routines used for discretization of regular sphere packings (of BCC, FCC, and SC types) in the work I co-authored. Their purpose is to reduce the scatter of output transport coefficients, drag force or permeability. The main routine is pack_discretization.m with documentation located in its comments. This routine calls binary2vector.m for initialization of sphere coordinates and calc_rad_scal_factor.m calculating radius for a given porosity value.

An optional routine isBoxBad.m helps to check the ratio between the edge *E* of a lattice and the number of unit cells *U* to be replicated along each Cartesian direction — both are input parameters for pack_discretization.m. If isBoxBad.m returns one for a given pair (*E*, *U*), other value of *E* should be taken, otherwise results of a flow simulation for the given lattice edge *E* can be obtained using smaller number of unit cells.