Cheviakov Alexei F, Zawada Daniel
Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada S7N 5E6.
Phys Rev E Stat Nonlin Soft Matter Phys. 2013 Apr;87(4):042118. doi: 10.1103/PhysRevE.87.042118. Epub 2013 Apr 22.
A narrow-escape problem is considered to calculate the mean first passage time (MFPT) needed for a Brownian particle to leave a unit sphere through one of its N small boundary windows (traps). A procedure is established to calculate optimal arrangements of N>>1 equal small boundary traps that minimize the asymptotic MFPT. Based on observed characteristics of such arrangements, a remarkable property is discovered, that is, the sum of squared pairwise distances between optimally arranged N traps on a unit sphere is integer, equal to N(2). It is observed for 2≤N≤1004 with high precision. It is conjectured that this is the case for such optimal arrangements for all N. A dilute trap limit of homogenization theory when N→∞ can be used to replace the strongly heterogeneous Dirichlet-Neumann MFPT problem with a spherically symmetric Robin problem for which an exact solution is readily found. Parameters of the Robin homogenization problem are computed that capture the first four terms of the asymptotic MFPT. Close agreement of asymptotic and homogenization MFPT values is demonstrated. The homogenization approach provides a radically faster way to estimate the MFPT since it is given by a simple formula and does not involve expensive global optimization to determine locations of N>>1 boundary traps.
考虑一个九死一生的问题,以计算布朗粒子通过其N个小边界窗口(陷阱)之一离开单位球体所需的平均首次通过时间(MFPT)。建立了一种程序来计算N>>1个相等小边界陷阱的最优排列,以使渐近MFPT最小化。基于此类排列的观测特征,发现了一个显著特性,即在单位球体上最优排列的N个陷阱之间的成对距离平方和为整数,等于N(2)。在2≤N≤1004范围内进行了高精度观测。推测对于所有N的此类最优排列都是如此。当N→∞时,均匀化理论的稀陷阱极限可用于将强非均匀狄利克雷 - 诺伊曼MFPT问题替换为一个球对称罗宾问题,该问题的精确解很容易找到。计算了罗宾均匀化问题的参数,这些参数捕获了渐近MFPT的前四项。证明了渐近MFPT值和均匀化MFPT值的紧密一致性。均匀化方法提供了一种从根本上更快的方式来估计MFPT,因为它由一个简单公式给出,并且不涉及确定N>>1个边界陷阱位置的昂贵全局优化。