動作環境

GeForce GTX 1070 (8GB)
ASRock Z170M Pro4S [Intel Z170chipset]
Ubuntu 16.04 LTS desktop amd64
TensorFlow v1.2.1
cuDNN v5.1 for Linux
CUDA v8.0
Python 3.5.2
IPython 6.0.0 -- An enhanced Interactive Python.
gcc (Ubuntu 5.4.0-6ubuntu1~16.04.4) 5.4.0 20160609
GNU bash, version 4.3.48(1)-release (x86_64-pc-linux-gnu)
scipy v0.19.1
geopandas v0.3.0
MATLAB R2017b (Home Edition)
ADDA v.1.3b6

This article is related to ADDA (light scattering simulator based on the discrete dipole approximation).

関連: Superposition T-matrix codes > bisphere.f > 2つの球の光散乱計算コード (analytical orientation averaging)

ADDA側でbisphereの計算を実行しようとしているが、実行引数について不明な点が出ている。

@ manual.pdfのpage 15

bisphere - two identical spheres with outer diameter d and center-to-center distance R cc
(along the z-axis). It describes both separate and sintered spheres.

実行時ヘルプを見てみると

$ ./adda -shape bisphere
ERROR: Illegal number of arguments (0) to '-shape bisphere' option (1 expected)
Usage: -shape bisphere <R_cc/d> 
Type 'adda -h shape bisphere' for details

<R_cc/d>に指定する値により、２つの球の位置情報が変化する。
接触状態にするには0.5だろうか。

位置情報の可視化

Jupyter | Matplotlib > 凹凸の激しい形状の可視化 > v0.4, v0.5: ADDA InitField-Yファイル読込み > Sphere | Prism (N=3) | Chevyshev particles
で実装したコードを用いて接触状態かどうかの確認をする。

InitField-Yの生成

ADDAの実行時引数「-store_int_field」を指定することで、IntField-Yを作る。

$ ./adda -shape bisphere 0.5 -store_int_field
...
$ ls run22099_bisphere_g16_m1.5/
CrossSec-Y  IntField-Y  log  mueller

上記のリンク先のコード: showChebyshev_170910.ipynbを実行する。
描画までに30秒程度かかる。

結果

$./adda -shape bisphere 0.5 -store_int_field

$./adda -shape bisphere 1.0 -store_int_field

2つの球が接触状態になるbisphereパラメータ

上記の画像より、以下のパラメータにすることで２つの球が接触状態になることが分かった。

$./adda -shape bisphere 1.0

Possible correction to the manual.pdf page 15

bisphere - two identical spheres with outer diameter d and center-to-center distance Rcc
(along the z-axis). It describes both separate and sintered spheres.

shall better be:

"bisphere - two identical spheres with composing sphere diameter d and center-to-center distance Rcc
(along the z-axis). It describes both separate and sintered spheres."