python进行数值模拟代码-数值模拟方法.pdfVIP

python进⾏数值模拟代码_数值模拟⽅法 数值模拟⽅法 科学研究与解决⼯程问题的基础在于物理实验与实物观测， 例如对⾦属材料的 凝固过程 进⾏物理实验、对天体运⾏进⾏观测。现代科学研究⽅法的核⼼是通 过实验或观测建⽴研究 对象的数学模型，基于数学模型进⾏研究与分析。这种 研究⽅法可以追溯到伽利略的⼯作， 成熟于⽜顿的三⼤定律与微积分。 采⽤实物 模型进⾏物理实验的研究周期 长、投⼊⼤，有时甚⾄⽆法在实物上进⾏，如天 体物理的研究。在数学模型上进⾏的数值模拟研究具有研究周 期短、安全、投 ⼊少，已经成为不可或缺的⼯具。 数值模拟⽅法的应⽤对象分为三个层次 1) 宏观层次：常见的⼯程建筑、制造设备、零件等 ; 2) 界观层次： 材料的微观组织与性能， 如⾦属材料的晶粒度影响其屈服强度 ; 3) 微观层次：基本物理现象与机理，如⾦属材料凝固时的结晶与晶粒⽣长 过程。宏观与界观层次的数值模拟⽅法包括： 1) 有限差分⽅法 (Finite Difference Method, FDM) ：微分⽅程的直接离散 ⽅法 ; 2) 有限元单法 (Finite Element Method, FEM) ：⽤有限尺度的单元的集合来 代替连续体， 分 Lagrange ⽅法， Euler ⽅法， ALE ⽅法 ; 3) 边界单元⽅法 (Boundary Element Method, BEM) ：⼀种半解析⽅法 ; 4) 有限体积⽅法 (Finite Volume Method, FVM) ：把空间划分成有限尺度的 体积单元，连 续体通过这些在空间上固定的体积单元，单元的空间位置不变 ; 5) ⽆⽹格⽅法 (Meshless Method) ：只布置结点，不需要划分单元⽹格，有 权函数。 微观层次的数值模拟⽅法包括： 1) 第⼀原理法 (First Principle Simulation) ：量⼦⼒学⽅法，直接计算原 ⼦的电⼦结构 ; 2) 元胞⾃动机⽅法 (Cellular Automata) ：把空间⽤元胞演化、元胞的局部 相互作⽤来描述复杂的、全局的系统。 3) 蒙特卡洛⽅法 (Monte Carlo Method ) ：把颗粒运动定义为随机过程，⽤ 势能的变化来 判断颗粒运动能否被接受。 4) 分⼦动⼒学⽅法 (Molecular Dynamics) ，分为经典⽅法、嵌⼊原⼦模型 (Embedded Atom