[論文レビュー] URDME v. 1.1: User's manual

Brian Drawert, Stefan Engblom+1 more|arXiv (Cornell University)|Feb 17, 2009

Gene Regulatory Network Analysis参考文献 28被引用数 1

ひとこと要約

URDME v.1.1 は、有限要素法に基づくジャンプ係数を用いて拡張された次世代サブボリューム法（NSM）を用いて、非構造格子上での確率的反応拡散シミュレーションを向上させた。主な改善点には、Comsol Multiphysics 3.5a へのネイティブ対応および MATLAB インターフェースから分離されたコアソルバーのスタンドアロン実行可能ファイル化が含まれ、複雑な幾何形状における大規模なアンサンブルシミュレーションのための効率的でスケーラブルなワークフローを実現した。

ABSTRACT

We have developed URDME, a general software for simulation of stochastic reaction-diffusion processes on unstructured meshes. This allows for a more flexible handling of complicated geometries and curved boundaries compared to simulations on structured, cartesian meshes. The underlying algorithm is the next subvolume method (NSM), extended to unstructured meshes by obtaining jump coefficients from the finite element formulation of the corresponding macroscopic equation. This manual describes the changes in the software compared to the previous version, URDME 1.0. Notable changes include support for Comsol Multiphysics 3.5a and a stronger decoupling between core simulation routines and the Matlab interface by making core solvers stand-alone executables. This opens up for more flexible simulation workflows when many realizations of the stochastic process need to be generated.

研究の動機と目的

• 複雑で直交座標系でない幾何形状における確率的反応拡散プロセスのシミュレーションの柔軟性を向上させること。
• 有限要素法による定式化を用いて、非構造格子上への次世代サブボリューム法（NSM）の拡張を実現し、正確なジャンプ係数計算を可能とすること。
• MATLAB インターフェースからコアシミュレーションルーチンを分離することで、ソフトウェアの拡張性を高めること。
• 複数の実現を効率的に生成できるようにすることで、高スルーレットのシミュレーションワークロードをサポートすること。
• メッシュ生成および事前処理のための Comsol Multiphysics 3.5a との互換性を確保すること。

提案手法

• マクロな反応拡散方程式の有限要素法による定式化から、非構造格子上におけるジャンプ係数を導出することで、次世代サブボリューム法（NSM）を非構造格子に適応させた。
• 空間的に変化する拡散ジャンプレートを計算するために有限要素法を用い、これにより元となる確定的 PDE と整合性を保った。
• コアシミュレーションソルバーをスタンドアロン実行可能ファイルにコンパイルすることで、実行時の MATLAB 依存を低減した。
• ソフトウェアアーキテクチャにおいて、シミュレーションロジックと MATLAB インターフェースを分離し、モジュラリティとパフォーマンスを向上させた。
• 複雑なドメイン表現を可能にするために、メッシュのインポートおよび幾何形状処理に Comsol Multiphysics 3.5a を統合した。
• スタンドアロンソルバーの効率的かつ並列実行を可能にする仕組みにより、複数の実現に対する確率的プロセスのシミュレーションをサポートした。

実験結果

リサーチクエスチョン

• RQ1 次世代サブボリューム法は、非構造格子上に拡張可能であり、かつ正確性を保持できるか？
• RQ2 コアソルバーを MATLAB から分離することで、シミュレーションのパフォーマンスおよびワークフローの柔軟性にどのような影響を与えるか？
• RQ3 Comsol Multiphysics 3.5a の統合によって、生物学的および物理的幾何形状のメッシュ処理がどの程度向上するか？
• RQ4 スタンドアロンコアソルバーは、スケーラブルで高スルーレットの確率的プロセスシミュレーションを可能にするか？
• RQ5 有限要素法に基づくジャンプ係数計算は、不規則なグリッド上での数値的安定性および正確性をどのように向上させるか？

主な発見

• 次世代サブボリューム法の非構造格子上への拡張は、有限要素法に基づくジャンプ係数の導出により実現され、マクロな拡散方程式と整合性を保った。
• コアソルバーをスタンドアロン実行可能ファイルに分離することで、MATLAB 依存のない効率的な実行が可能となり、大規模シミュレーションにおけるパフォーマンスとスケーラビリティが向上した。
• Comsol Multiphysics 3.5a への対応により、複雑な非構造メッシュのシームレスなインポートが可能となり、生物学的および物理的ドメインへの適用性が向上した。
• 新しいアーキテクチャにより、MATLAB のオーバーヘッドなしに複数の実現を並列実行できるため、高スルーレットのシミュレーションワークフローが実現された。
• 有限要素法による定式化により、不規則なメッシュ上でも正確で安定したジャンプ係数が得られ、シミュレーションの忠実性が向上した。
• アップデートされたソフトウェアにより、以前のバージョンと比較して、複雑な幾何形状におけるより柔軟でスケーラブルな確率的反応拡散シミュレーションが可能になった。