The present invention discloses a method, device and computer readable storage medium (CRSM) for generating a volume foliation. The method comprises: constructing a pants decomposition graph on a smooth closed surface S of genus g>1, wherein the surface S has g handles, and the pants are a genus-zero triangular mesh surface with g boundaries; constructing an initial mapping f0:(T,g)→(Gr, h) from a tetrahedral mesh T to a graph Gr by a process of extending a topological disk Dk into a solid cylinder Ck; evolving the initial mapping into a generalized harmonic map by the nth iterative mapping; and obtaining a volume foliation of the tetrahedral mesh T from preimages of all nodes of a metric graph (Gr, h) and points on edges of the metric graph (Gr, h) under the generalized harmonic map f.