- 1、本文档共41页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
水动力学计算 薛雷平
SIP: 要求 M? ? A ? (A?)P = 且设 对比系数有 4大型稀疏方程组迭代求解 NW N W P E S SE 4大型稀疏方程组迭代求解 代入 4大型稀疏方程组迭代求解 可得 ?=0 时,为普通的非完全三角分裂法 4大型稀疏方程组迭代求解 SIP方法的迭代矩阵M,可结合共轭梯度类方法结合,作为CG法(对称系数矩阵)以及BiCGSTAB(非对称系数矩阵)的前置优化矩阵使用。尤其对于多块结构化网格,可以取得更快的迭代收敛速度。 * * 谢谢! * * BICGSTAB主要思想 * * * * BICGSTAB主要思想 BICGSTAB计算流程 * * BICGSTAB计算流程 * * CG 计算流程 * * 现代密码学理论与实践之五 水动力学计算若干问题2012年研究生暑期学校上海交通大学 薛雷平 提纲 数值计算要点 有限体积法 压力修正法 大型稀疏方程组迭代求解 * * 1 数值计算要点 基本模型的选择 计算方法的选择 边界条件和计算流程的选择 网格生成及其收敛性研究 结果分析和总结 * * 1 数值计算要点 基本模型的选择 可压缩非定场粘性流动方程 不可压缩NS方程 + 湍流模型(RANS) 势流方程(理想流体) 边界层方程 2D简化、几何简化、定解条件简化 … 适用范围 ? 简化模型的误差估计? 可通过典型算例不同层次的简化计算 得到合适的基本计算模型 * * 1 数值计算要点 计算方法的选择 离散方法:有限体积,有限差分,有限元… 计算格式: 伪可压缩模型的时间推进法(显式), 压力修正法(隐式)… 湍流模型、多相流及空化模型、数值造波\消波… * * 1 数值计算要点 边界条件和计算流程的选择 边界位置及给出的边界条件并不一定要与实际问题完全一致 通过物理意义的考虑和试算得到能够简单实现, 但又不使实际问题失真的边界条件 注意点:来流速度分布和湍流度 非远场出口边界条件 对称几何的非对称流动 空化流动时的压力条件的实现 计算流程:对于复杂问题,通常可以先计算物理、几何简化的初始解,逐步加入复杂性 * * 1 数值计算要点 网格生成 截断误差取决与于 网格间距 和 解的梯度,合理分布网格,通过初步计算重新分布网格通常是必要的 粘性流动计算时边界层网格(最好是结构化网格), 不同湍流模型对y+的要求 网格收敛性研究 网格收敛性研究指采用不断加密的网格使解收敛 不是指采用某一密度的网格时解与实验值一致 这是获得可靠数值计算结果必不可少的步骤 * * 1 数值计算要点 结果分析和总结 数值计算结果分析与实验结果分析类似 尽量采用无量刚表达式 CFD不是colour fluid dynamics,少用云图,多用曲线 分析总结的层次: 1:初步的定性规律(相关性和大小) 2:进一步的定性规律(线性、幂次、指数…,先用量纲分析和理论储备获得基本规律,多变量时考虑采用变量分离法) 3:定量规律(系数的拟合,加入特征点) 4:分析定量规律背后的物理机理 * * * * 2 有限体积法 流体力学基本方程的一般形式 对流扩散方程 随体变化=时间变化+对流 = 扩散项+体积源项 +面积源余项 * * 2 有限体积法 对流扩散方程的体积分形式 V为任意空间体积(控制体),F为体积V的表面 该输运方程可由 - 微分方程体积分后直接由GAUSS积分公式得到, - 也是控制体中物理量变化规律的描述 * * 2 有限体积法 将积分方程应用于一个微元有限体积 * * 2 有限体积法 对流项的离散 迎风格式 + 高阶修正 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 扩散项的离散——梯度格式 * * 2 有限体积法 有限体梯度项的数值计算 * * 2 有限体积法 时间变化项的离散 2 有限体积法 时间变化项的离散 与对流项合并(对任意f:
文档评论(0)