




已阅读5页,还剩50页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1,(三)偏微分方程的数值离散方法,3.1有限差分法3.2有限体积法(有限元,谱方法,谱元,无网格,有限解析,边界元,特征线),2,3.1有限差分法,3.1.1模型方程的差分逼近3.1.2差分格式的构造3.1.3差分方程的修正方程3.1.4差分方法的理论基础3.1.5守恒型差分格式3.1.6偏微分方程的全离散方法,3,3.1.1模型方程的差分逼近,4,3.1.2差分格式的构造,5,3.1.3差分方程的修正方程,差分方程所精确逼近的微分方程称为修正方程对于时间发展方程,利用展开的方程逐步消去带时间的高阶导数,只留空间导数。Warming-Hyett方法:差分方程(2)写成算子的形式:,6,3.1.3差分方程的修正方程(续),7,3.1.3差分方程的修正方程(续),8,3.1.4差分方法的理论基础,相容性,稳定性,收敛性等价性定理Fourier稳定性分析,9,3.1.4差分方法的理论基础(续),Fourier(VonNeumann)稳定性分析,10,3.1.4差分方法的理论基础(续),Fourier(VonNeumann)稳定性分(续)称为CFL条件(Courant,Friedrichs,Levy),11,3.1.5守恒型差分格式,流体力学方程组描述物理量的守恒性;守恒律组:定义,12,3.1.5守恒型差分格式(续),守恒性质:非守恒的差分格式一般没有对应于原始守恒律的“离散守恒律”。,13,3.1.5守恒型差分格式(续),守恒型差分格式的Lax-Wendroff定理:如果守恒型差分格式是和守恒律相容的,且当时间和空间步长趋于零时,差分解一致有界,几乎处处收敛于分片连续可微的函数,则这个收敛的函数就是守恒律的一个弱解。推论:守恒型差分各式的收敛解能自动满足间断关系。用途:(加上熵条件)可以得到正确的激波,研究中大量使用例如:Lax-Friedrichs格式,Lax-Wendroff格式,MacCormack格式,14,3.1.6偏微分方程的全离散方法,对差分格式的一般要求:有精度、格式稳定、求解效率高特殊要求物理定律(守恒性)、物理特征(激波、湍流、旋涡、多介质、化学反应等)、有界性(正密度、正温度、正湍动能、正组分浓度等)主要指非定常方程的时间离散,15,3.1.6偏微分方程的全离散方法(续),两层格式Crank-Nicolson格式、P-C格式、Lax-Wendroff格式、MacCormack格式Runge-Kutta方法时空全守恒:如Godunov格式、central-upwind格式、CESE方法多层格式Leap-Frog格式、Adams-Bashforth格式、后三点隐格式,16,3.1.6.1两层格式,Crank-Nicolson格式Predictor-Corrector格式Lax-Wendroff格式MacCormack格式Runge-Kutta方法,17,3.1.6.1两层格式(cont.),Lax-Wendroff格式一步LW格式,18,3.1.6.1两层格式(cont.),Lax-Wendroff格式两步LW格式常系数Jacobian时与单步LW等价。但计算更简单,不涉及矩阵相乘。,19,3.1.6.1两层格式(cont.),MacCormack格式(1969)两步格式比LW更简单,不需要计算函数在半点上的值。LW两步格式和MC各式的缺点:定常解的误差依赖于时间步长。,20,MacCormack格式的构造,21,3.1.6.2三层格式,Leap-Frog格式Adams-Bashforth格式,22,第二课后阅读提示,傅德薰计算流体力学,3.13.3水鸿寿一维流体力学数值方法3.1ComputationalMethodsforFluidDynamics,FerzigerandPeric,SpringerChap.6,23,作业2,1.用Fourier法分析3.1.6.1节中Crank-Nicolson格式的稳定性。2.分析前面3.1.6节中MacCormack格式是几阶精度。,24,3.2有限体积法,出发方程为积分型守恒方程(直角坐标、柱坐标、球坐标)以控制体为离散量计算体积分和面积分需要适当的插值公式和积分公式(quadratureformula)适用于任意形状的网格,复杂几何形状缺点:难以构造大于二阶以上的格式,25,3.2.1定常守恒型方程和控制体,26,3.2.2面积分的逼近,面积分用积分点的值表示(quadrature)积分点的值用CV的值表示(interpolation)对于Simpson公式,对积分点的插值需要四阶精度,27,3.2.4体积分的逼近,当被积函数为某种型函数时,可以得到精确的积分,逼近精度取决于型函数的精度。,28,3.2.4体积分的逼近,四阶精度:2D直角坐标网格最后一式可以四阶精度逼近3D的面积分,29,3.2.5插值和微分,积分点的函数值和其法向梯度1stUDS:取上风点的值,30,插值,2ndorder:向积分点线性插值等价于中心差分(CDS),31,插值,当积分点的函数是线性插值时Secondorder,32,插值,QUICK(quadraticupwindinterpolationforconvectivekinematics)插值三阶精度,但积分(差分)往往只有二阶精度。,33,插值,高精度:N阶精度的quadrture需要N-1阶多项式插值公式。界面上导数可以用插值公式的微分求出。,34,3.2.5有限体积法的边界条件,用边界条件替代面积分入口:通常给定对流通量(mass,momentum,energy,etc.)壁面和对称面:通量为零边界上函数值给定:和内部CV的值共同构建边界上的导数,35,FV例子,36,3.2.6守恒律的有限体积方法Godunov格式,37,38,3.2.6.1Godunov方法的思想,39,一阶迎风格式(CIR格式),40,用Godunov思想说明CIR格式=Godunov格式,41,42,Riemann解图示,43,44,3.2.6.11DEuler方程组的Godunov格式,Godunov格式是基于积分形式的方程组,间断关系自动满足,不需要另外考虑间断线上的间断关系,45,移动网格上的积分回路,46,移动网格上的Godunov格式,47,固定网格上的Godunov格式,48,Lagrange网格上的Godunov格式,49,Euler方程组的Riemann问题的解理想气体的5种解,50,51,二维Euler方程组的Riemann问题,52,53,仅是局部化的1DRP,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025至2030中国电焊头盔行业产业运行态势及投资规划深度研究报告
- 2025至2030中国电子设备维修服务行业产业运行态势及投资规划深度研究报告
- 2025至2030中国甜茶叶子市场投资风险及运行状况预测研究报告版
- 2025至2030中国环境质量检测行业市场发展分析及竞争格局与投资前景报告
- 培训需求调查课件
- 餐饮服务培训课件
- 儿童健康成长之路从骨关节健康知识普及开始
- 智慧教育新篇章技术如何重塑学习成效
- 学习者的创新思维培养与实践
- 那智智能技术助力商业高效运营与决策
- 初一生活学习指导
- 下肢静脉曲张
- 2024年露营帐篷项目可行性研究报告
- 《公务员录用体检操作手册(试行)》
- 2024粤东西粤北地区教师全员轮训培训心得总结
- 2024-2025学年华东师大版数学七年级上册计算题专项训练
- 福建省机关工作人员年度考核登记表
- JBT 7808-2010 无损检测仪器 工业X射线探伤机主参数系列
- DB44-T 2474-2024 自然教育标识设置指引
- 研学基地合作协议
- 驾驶员行为规范管理制度
评论
0/150
提交评论