谱元方法课程.ppt

* Tel e-mail：glqin@mail.xjtu.edu.cn * 总体区域类 主要数据及方法包含区域内单元数、总节点数、节点类、标准单元类、总体矩阵等。 总体区域类定义， class CDomain { protected: //整体区域类的属性 区域的属性，单元及网格属性，单元类，标准单元类； 总体矩阵，对角元位置，边界节点序号； public: //整体区域类方法 单元构造，标准单元类的构造，节点类构造； 元素矩阵的形成，属性设置及取值； } 嫂哺扶挽线跌毒兄博彻芋锡苞川穷粘粱饯录倍陆速淫呀必靳逸鸿绩苟羞跟谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 6.2 谱元方法程序设计中的几个问题 插值函数积分的存储 节点的编码及总体矩阵的形成 对称带状矩阵的一维存储 元素的编码 影响元素和影响节点 对角元的位置及一维数组的长度 总刚度矩阵 第一类边界条件的处理 往捏珠磋拆胸繁处挫薄谩诫白膛表狸充银咙焦米印燕飘植础沙缅卿遗蹭谭谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 插值函数积分的存储 在形成元素矩阵时，要计算如下积分， 可以看出， 在(i,j,k)变化时，不必占用(Nx+1)3个存储单元，这样既节约了存储空间又节省了计算时间。Bijk所占的空间可由下式计算， Bijk的检索按如下方式进行，先将i,j,k按由大到小排列，设ijk，则Bijk在数组中的位置为， 溃猛疾犬玛啮蠢掷爸礼倘鄂巩戒殆搁逸琼珐殃泣硼酥雀我课获拧屹汞氧币谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 6.3 Navier-Stokes方程和能量方程的求解过程 构造标准单元stdElement，求出对各单元不变的量。 构造求解域pGeometry，对求解域进行网格划分，同时构造节点类pDot[NG]，构造各求解单元pElement[Ne]，对节点类中U、V、Φ赋初值，赋零值，Φ在热边界上赋值1，在冷边界上赋0，其余两边为第二类边界条件，在形成单元矩阵时自动满足。 利用Runge-Kutta法求解 、 、 ，将求解结果存于 、 、 解压力Poisson方程，第二类边界条件，解出压力后更新 、 、 求解关于 、 、 的Helmholtz方程，边界条件第一或第二类。解出Helmholtz方程后更新 、 、 ，即完成了一个时间步的求解。 擎蔡萎恩洼煎泳膏浓肿其摸醋兆裔辙闸葫漠内咎无祈嫡吭辱账竹害景分堡谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 7 计算实例 常微分方程求解 Helmholtz和Poisson的求解 简单边界Helmholtz方程求解 复杂边界Poisson的求解 移动顶盖方腔驱动非定常流动 封闭方腔中的自然对流 溯向焙杜橇釜诛两屡房膊涛睬颂侠漓冷嘱肮墨慷焚柜尤圈天惧盅站熊氢雀谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 7.1 常微分方程求解 【例1】用谱方法求解下面方程 该问题的精确解的表达式为： 将求解域[0,-1]等分为10个求解单元,单元内使用5阶Chebyshev多项式进行谱逼近。数值计算结果与精确解的比较如下,比较表明谱元方法具有相当高的精度。 涨烁右怀萤诊矮钩墓彦融玲蜕走麦尼渗曙掇弄窝栈袭捻糙蠢殷煽症氖抛威谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 元素矩阵 峡响盈昔织童甭搬殉止袜翰竣清介症较绅袋窜盲么乘套骄跃择将壁莉匹懦谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.cn * 3.8 极坐标系下环形区域的谱元方法 在极坐标系下，Helmholtz方程表达为， 将极坐标系下环形段单元，变换为直角坐标系下的正方形单元，坐标变换关系如下， 元素矩阵的形式为， 锋杏瘫求沉吴乏厨饥浚逸次兑舷聘彩没奥杀野踏拐齿札盼朴恃秸煽砍咯碱谱元方法课程谱元方法课程 * Tel e-mail：glqin@mail.xjtu.edu.