[論文レビュー] 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$ における最適な球体充填と、理想導電性球状不純物を有する複合材料における有効導電率の最小化との間の関係を確立する。問題をボロノイ＝デローラン図構造上の離散的エネルギー最小化として定式化することにより、球体充填密度の上界を導出し、構造的近似と $p$-ラプラシアン均質化を用いて最適な球体位置を求める構成的アルゴリズムを構築する。

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$ における最適な球体充填と周期的複合材料における有効導電率の最小化との間の対応関係を確立すること。
• $p$-ラプラシアン均質化における構造的近似フレームワークを $ℝ^d$ で開発し、充填効率のエネルギーに基づく推定を可能にすること。
• デローラングラフクラスと格子周期性に基づいて、最適な球体配置を特定する計算アルゴリズムを構築すること。
• 離散的エネルギー最小化を用いて $ℝ^d$ における最大球体充填密度の上界を導出し、積層構造に対して検証すること。
• エネルギー関数から導かれる分離可能な線形系を解くことにより、最適な球体位置を計算する構成的手法を提供すること。

提案手法

• 有効導電率問題を離散的エネルギー最小化として定式化：$\min_{t_j} \sum_{k,j} g_{kj}^{(0)} |t_k - t_j|^2$、ここで $g_{kj}^{(0)}$ は一次の粒子間フラックスを表す。
• $p$-ラプラシアン系における構造的近似理論を $ℝ^d$ に適用し、非線形導電率を仮定し、$p > \frac{d+1}{2}$ とすることで、粒子間距離 $\to 0$ のときのフラックス特異性を保証する。
• 周期セル $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$-ラプラシアン均質化に基づく離散的エネルギー関数は、幾何的充填配置とどのように関係するか？
• RQ3 デローラングラフクラスは、最適な球体配置と充填密度の上界を決定する上で果たす役割は何か？
• RQ4 エネルギー最小化問題の解が、巨視的に等方的構造をもたらす条件は何か？
• RQ5 構造的近似フレームワークは、非周期的または非球形不純物へと拡張可能か？その計算的影響は何か？

主な発見

• $\u211d^d$ における球体位置の離散的エネルギー関数の最小値は、まさに最適な球体充填配置に一致する。
• 積層構造では、導出された充填密度の上界が正確であり、これらのクラス内での解析的解を提供する。
• エネルギー最小化問題は、$n$ 個の独立した $d$ 次元線形系の解に帰着され、最適配置の効率的計算が可能になる。
• 最適な球体位置 $\mathbf{a}_k^*$ は、格子ベクトル $\boldsymbol{\nu}_\ell$ の線形結合として得られ、係数は座標別に解かれた系から計算される。
• 本手法により、$\u211d^d$ における最大球体充填密度の上界が得られ、積層配置などの構造的クラスでは正確な値が達成可能である。
• 巨視的等方性のための必要条件は、周期セル $Q_0$ の向かい合う面を接続する $d$ 個の透過経路（パーコレーションチェーン）の存在であり、これは数値的に確認可能である。