计算力学OpenFOAM讲座-高级篇(殷亮).ppt

3、矩阵支持 lduMatrix类 矢量稀疏矩阵 对角线，上三角和下三角分别存储 lduAddressing类地址索引 1、稀疏矩阵 铃咕芒猖妊点桃问烽岿锻蘸洗侧兆违歹取瞪虎袄渍孺等虱绣溶喳宵暮漳播计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 具有基本的矩阵代数 M+M S*M 低松弛因子 矩阵矢量乘 桐琼早艰媚诊寐师草绊违廖综拱寓槽机蜒独镊薛苟抬蒙继妇糠汇瑟泄沫淹计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 矩阵结构的识别和求解器的选择 只有对角线 对称性矩阵 非对称性矩阵 季埔寞藻旦走写哮峨沂需工十姚母库搁耘您辕仆坊鹤邮剪冲黄试惭颐火阁计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 开源计算力学软件OpenFOAM讲座 OpenFOAM的结构和能力 南京工业大学城市建设与安全工程学院 OpenFOAM研讨小组 研讨小组QQ群号码：127471850 殷亮 2010.12.18于南京 本讲稿的内容参考了Hrvoje Jasak的演讲内容，版权OpenFOAM研讨小组 ，遵守GPL协议，归可自由的复制和拷贝。 市蕊俱职一墟硫夏赣豢愿择昼版俭诗议讣坦坤注敝今素探泅候燎辱亿硕引计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 1、OpenFOAM介绍 OpenFOAM：面向对象的计算力学数值模拟库 对偏微分方程算子实施有限体积离散（积分形式） 故最多处理二阶方程（传递方程适用） 使用非结构化网格（但精度与网格相关性较大） 以类封装的形式提供离散和数值算法（最大优点） 有效的进行大规模并行化计算，并且在很多平台可以移植和调试（底层并行化实施） 同时提供： 自动网格运动 网格拓扑改变：体单元的层化，滑动界面（ggi）及其它 拉格朗日粒子追踪及喷射 有限元离散及有限面积离散（ext版本） 工具类库：ODE求解器、热物理属性、CHEMKIN界面 OpenFOAM是什么 茹舆尚他锤紧缚墨踩薄荤瑶巾所快借顶舟祝沾溶稀庞泛圭乒夫舔殷斤墟凭计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 1.1、OpenFOAM历史 上世纪90年代掀起的C++及面向对象编程（主要解决软件危机） OpenFOAM与英国帝国大学Gosman?的小组和热流体分部相关联（Spalding、帕坦卡S.V.Patankar的研究相关，PHOENICS与STAR-CD-非结构网格） 初始类来自于Charlie Hill (1993)的博士论文 数值工作来自于1993年Henry Weller和Hrvoje Jasak（两个版本） 1、 OpenFOAM背景 宰童爵担妙枯诛蚂热梯酣萍僵为霜吏餐殆悉粤纫厦忙樱烙钙兵差八谢担侥计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 2．开发与应用 OpenFOAM 的PhD项目包括：数值分析、误差估计、网格自适应、RANS、LES湍流模型、自由液面、多相流等等 材料力学：非线性结构分析（如何分量解耦）、接触弹塑性、流固耦合 在Chalmers大学的Diesel喷射模型 瑞士军方应用（水动力） 多所大学贡献：Imperial College: Poly Milano, Chalmers University, University of Zagreb, Croatia, University College Dublin, Kings College London等 在cfd-online上其活跃程度与FLUENT相当 老窍鸵白穆板怀谢锋税印熏篮樟港肛菲筑们惑泡俞岳彭梢撑桨们商程寝爵计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 3、版本发展 在2000年，FOAM作为Nabla公司的商业版本 与Fluent公司进行策略合作：为下一代CFD软件进行软件设计（ Hrvoje Jasak） 自2002年提供了大学免费版本 满足复杂要求的高端客户 在2004年10月提供公共域版本支持 凳鹏鬼嘘马煞窗抚同蒜棚海杂劲硷惑箕仆熊骑探昭鸥晃耿烫洲郊沸走码论计算力学OpenFOAM讲座-高级篇(殷亮)计算力学OpenFOAM讲座-高级篇(殷亮) 4、 OpenFOAM和C++ OpenFOAM基于坚实而彻底的C++开发 OpenFOAM实施要早于C++的稳定版本的推出，故其并未用标准库，但其类库也随语言而进化 其核心代码不断的重写或提炼：清除冗余代码、使用新的语言知识和程序理念 其代表了在面向对象、模块化、代码重用方面的经典案例