计算力学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++的稳定版本的推出，故其并未用标准库，但其类库也随语言而进化 其核心代码不断的重写或提炼：清除冗余代码、使用新的语言知识和程序理念 其代表了在面向对象、模块化、代码重用方面的经典案例