[论文解读] User Guide for the Discrete Dipole Approximation Code DDSCAT.6.0

B. T. Draine, Piotr J. Flatau|arXiv (Cornell University)|Sep 2, 2003

Particle accelerators and beam dynamics参考文献 1被引用 93

一句话总结

本文介绍了DDSCAT.6.0，这是一套免费的Fortran软件包，采用离散偶极近似（DDA）计算任意形状和复折射率的粒子的电磁散射与吸收。该软件利用快速傅里叶变换（FFT）和MPI实现高性能计算，可对尺寸参数小于15且折射率接近1的靶标进行精确模拟。

ABSTRACT

DDSCAT.6.0 is a freely available software package (http://www.astro.princeton.edu/~draine/DDSCAT.6.0.html) which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. DDSCAT.6.0 allows accurate calculations of electromagnetic scattering from targets with ``size parameters'' 2*pi*a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). DDSCAT.6.0 includes the option of using the FFTW (Fastest Fourier Transform in the West) package. DDSCAT.6.0 also includes MPI support, permitting parallel calculations on multiprocessor systems. DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to ``import'' arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.6.0 to carry out electromagnetic scattering calculations. CPU and memory requirements are described.

研究动机与目标

• 为DDSCAT.6.0提供全面的用户指南，该软件包是模拟复杂粒子电磁散射的广泛使用工具。
• 使研究人员能够对具有任意几何形状以及非均质或各向异性材料的非球形粒子执行精确的散射与吸收计算。
• 通过与FFTW和MPI集成，支持高性能计算，提升大规模模拟的计算效率。
• 通过提供模块化、可移植的Fortran源代码和清晰的文档，确保广泛可访问性与可扩展性，便于用户自定义和导入目标。
• 通过公共代码仓库和用户注册机制，鼓励引用、反馈与贡献，促进可重现性与社区使用。

提出的方法

• 离散偶极近似（DDA）将目标建模为极化点偶极子的晶格，通过迭代求解器求解由此产生的线性方程组。
• 该代码使用West地区最快傅里叶变换（FFTW）库，加速迭代求解过程中的矩阵-向量乘法。
• MPI（消息传递接口）支持使得程序可在多处理器系统上并行执行，显著缩短大规模偶极子阵列的计算时间。
• 该软件支持多种目标几何形状（如椭球体、圆柱体、四面体），并通过自定义输入例程支持用户定义的目标导入。
• 输出包括总吸收与散射截面、指定散射方向和入射偏振态下的穆勒散射矩阵元素和斯托克斯参数。
• 代码生成与机器无关的netCDF输出文件，便于数据交换，并通过readnet.pro和mie.pro等IDL实用程序支持基于IDL的后处理。

实验结果

研究问题

• RQ1 如何在可移植的高性能Fortran代码中高效实现用于任意粒子几何形状的离散偶极近似？
• RQ2 哪些计算策略（如FFT和MPI）是加速电磁散射DDA模拟的最佳选择？
• RQ3 DDSCAT.6.0对尺寸参数小于15且折射率与1的差异小于1的靶标，其散射与吸收建模的精度如何？
• RQ4 在不同计算平台上配置、编译和运行DDSCAT.6.0的实际考虑因素有哪些？
• RQ5 用户如何在不修改核心算法的前提下，扩展代码以支持新的目标形状、材料或散射配置？

主要发现

• DDSCAT.6.0可对尺寸参数小于15且折射率与1的差异小于1的靶标实现精确的电磁散射计算。
• FFTW的集成显著加速了迭代求解过程，使得在标准系统和高性能计算系统上均可实现大规模模拟。
• MPI支持实现了可扩展的并行执行，在多核和集群系统上显著缩短运行时间，且在多处理器环境下观察到良好的性能扩展性。
• 该代码支持ASCII和netCDF两种输出格式，推荐使用netCDF格式以实现可移植性和随机访问的数据处理。
• 包含readnet.pro和mie.pro等IDL实用程序，可实现有效的后处理，并能与球形单分散粒子的Mie理论结果进行对比。
• 该软件高度模块化且可扩展，用户可仅通过少量代码修改即可导入自定义目标几何形状并扩展功能。