计算流体力学第四部分天津大学.pdfVIP

有限元与有限差分的主要区别 （1）有限差分以结点近似为基础，构造差分、导数之间的关系， 用网格结点值近似连续函数解；有限元以单元近似为基础，通过单 元积分构造函数、导数之间的关系，用插值函数近似连续函数解。 （2 ）有限差分只保证解在结点处的连续光滑，而有限元保证了解 在单元上的连续和光滑。 （3）有限差分解的收敛性取决与结点处的最大误差，而有限元解 的收敛性由单元上的误差权积分来决定，具有单元平均误差的含义。 （4 ）有限差分所适应的边界规整，边界处易形成锯齿状格点分， 局部加密较困难；有限单元能适应复杂边界要求，边界拟合程度较 好，局部加密方便灵活。  差分方法与有限元法的网格信息比较 信息量较小 信息量较大 25 20 15 10 5 0 0 5 10 15 20 25  1试函数与基函数  Au f ∈Ω  ∂u  (au +b ∂n ) au0 +bg  Γ 0 u u  t  0  第一类边界条件（狄里克雷即本质性边界条件） a 1,b 0  第二类边界条件（诺伊曼即自然边界条件 ） a 0,b 1  第三类边界条件（混合边界条件 ） a 1,b 1  时间上的初始条件  2余量法  （1）配置法  设方程为 Au −f 0  数值解表达的真解为 u(x ,y ,z ,t) uM (x ,y ,z ,t) +R (x ,y ,z ,t)  代入方程后 Au −f AuM +AR −f 0  剩余量为 Re −AR  方程为 Au −f R M e  数值解是单元解的线性组合，空间与时间变量可分离 M ( ) ( , , ) uM ∑Ci t ϕi x y z i 1 基函数  选择M个单元，确定M个uM ，最简单方法是在单 元中确定一点i配置给uMi ，则  M Au −f A∑C (t )ϕ (x , y , z ) −f R  M 1 1 i 0 i 1 1 1 1 e1 包含时间和空间变量属于初边值问题  i 1  M 只包含空间变量属于边值问题 Au −f A∑C (t )ϕ (x , y , z