[论文解读] Optimal packing of spheres in $\mathbb R^d$ and extremal effective conductivity

Vladimir Mityushev|arXiv (Cornell University)|Dec 23, 2014

Advanced Mathematical Modeling in Engineering参考文献 6被引用 3

一句话总结

本文建立了 $ℝ^d$ 中最优球体密堆积与具有理想导电球形夹杂物的复合材料中有效导电率最小化之间的联系。通过将问题表述为基于Voronoi-Delaunay结构的离散能量最小化，推导出球体密堆积密度的上界，并提出一种基于结构近似和 $p$-Laplacian 同伦的构造性算法，用于确定最优球体位置。

ABSTRACT

Optimal packing of spheres in $\mathbb R^d$ is studied by optimization of the energy $E$ (effective conductivity) of composites with ideally conducting spherical inclusions. It is demonstrated that the minimum of $E$ over locations of spheres is attained at the optimal packing. The energy is estimated in the framework of structural approximations. This method yields upper bounds and sometimes exact values for the maximal concentrations of spheres in $\mathbb R^d$. A constructive algorithm for the optimal locations of spheres associated to the classes of the Delaunay graphs is constructed.

研究动机与目标

• 建立 $ℝ^d$ 中最优球体密堆积与周期性复合材料中有效导电率最小化之间的对应关系。
• 为 $ℝ^d$ 中 $p$-Laplacian 同伦构建一个结构近似框架，以实现基于能量的密堆积效率估计。
• 基于Delaunay图类和晶格周期性，构建一种计算算法，用于确定最优球体构型。
• 通过离散能量最小化推导出 $ℝ^d$ 中最大球体密堆积密度的上界，并在层状结构中加以验证。
• 提出一种构造性方法，通过求解从能量泛函导出的解耦线性系统，计算最优球体位置。

提出的方法

• 将有效导电率问题表述为离散能量最小化：$\min_{t_j} \sum_{k,j} g_{kj}^{(0)} |t_k - t_j|^2$，其中 $g_{kj}^{(0)}$ 表示粒子间的一阶通量。
• 将结构近似理论应用于 $ℝ^d$ 中的 $p$-Laplacian 系统，假设非线性导电率满足 $p > \frac{d+1}{2}$，以确保当粒子间距趋于零时通量奇点的出现。
• 在周期性胞元 $Q_0$ 上采用环面拓扑，利用基本向量 $\boldsymbol{\nu}_j$ 定义通过周期性识别的距离，以模拟无限周期性复合材料。
• 通过在准周期性势场和周期性通量条件下求解能量最小化，推导出最优球体位置 $\mathbf{a}_k^*$ 的线性系统（4.8）。
• 通过分离坐标，将系统分解为 $d$ 个独立的 $n$-维线性系统，从而实现高效的数值求解。
• 通过 $\mathbf{a}_k = \sum_{\ell=1}^d \mathbf{a}_k^{(\ell)} \boldsymbol{\nu}_\ell$ 构造最优构型，其中 $\mathbf{a}_k^{(\ell)}$ 是坐标分离系统（5.2）的解，保持晶格对称性。

实验结果

研究问题

• RQ1 能否通过具有理想夹杂物的复合材料中有效导电率的最小化来求解 $\u211d^d$ 中的最优球体密堆积问题？
• RQ2 基于 $p$-Laplacian 同伦的离散能量泛函与几何密堆积构型之间有何关系？
• RQ3 Delaunay 图类在确定最优球体排列及密堆积密度上界中起什么作用？
• RQ4 在何种条件下，能量最小化问题的解能产生宏观各向同性的结构？
• RQ5 该结构近似框架能否扩展到非周期性或非球形夹杂物的情况，其计算影响如何？

主要发现

• 在球体位置上对离散能量泛函的最小值恰好对应于 $\u211d^d$ 中的最优球体密堆积构型。
• 对于层状结构，推导出的球体密堆积密度上界是精确的，为这些类别提供了分析解。
• 能量最小化问题可简化为求解 $d$ 个大小为 $n$ 的独立线性系统，从而实现最优构型的高效计算。
• 最优球体位置 $\mathbf{a}_k^*$ 作为晶格向量 $\boldsymbol{\nu}_\ell$ 的线性组合获得，其系数通过坐标分离解求得。
• 该方法为 $\u211d^d$ 中最大球体密堆积密度提供了上界，且在如层状排列等结构类别中可达到精确值。
• 宏观各向同性的必要条件是周期胞元 $Q_0$ 的相对面之间存在 $d$ 条连通链，该条件可通过数值方法验证。