We consider using a discrete network model in combination with continuous nonlinear optimization models to solve the problem of optimizing channels in nanoporous materials. The problem and the hierarchical optimization algorithm are described in [2]. A key feature of the model is the fact that we use the edges of the finite element grid as the locations of the channels. The focus here is on the use of the discrete model within that algorithm. We develop several approximations to the relevant flow and a greedy algorithm for quickly generating a “good” tree connecting all of the nodes in the finite-element mesh to a designated root node. We also consider Metropolis-Hastings (MH) improvements to the greedy result. We consider both a regular triangulation and a Delaunay triangulation of the region, and present some numerical results.
