二十五对角系数矩阵的求解过程：三维线性福克-普朗克方程的直接数值解
计算物理
2009-11-11 v2 天体物理学
摘要
我们描述了一种用于求解三维（七个变量）线性福克-普朗克方程离散化所得线性方程组的隐式过程。福克-普朗克方程的离散化采用二十五点分子，导致系数矩阵具有相同数量的对角线。该方法是斯通隐式过程的扩展，涵盖了大量碰撞项，可应用于具有不同时间离散化的稳态或非稳态问题。我们在两个稳态示例中展示了测试计算及与其他方法的比较，包括一个针对典型星系的宫本-永井盘势的天体物理应用。
引用
@article{arxiv.physics/0505093,
title = {Solving procedure for a twenty-five diagonal coefficient matrix: direct numerical solutions of the three dimensional linear Fokker-Planck equation},
author = {Maximiliano Ujevic and Patricio S. Letelier},
journal= {arXiv preprint arXiv:physics/0505093},
year = {2009}
}
备注
20 pages, RevTex, no proofreading, accepted in Journal of Computational Physics