船舶在波浪中的运动-ch1_引论.ppt

上述三种情形，总结为： 称之为格林第三公式。其中， 代表偶极子， 代表点源，而 和 则分别相当于偶强和源强的分布密度。上式表明，场内任一点处 可 以由边界面上的分布奇点来表示。 * （4）赫斯-斯密斯方法 ——格林函数 只是 的一个特解，故 的具体形式并不唯 一。如在域内 存在处处调和的函数 , 则在域内处处满足 , 则 都是泊松方程 的解。因此，更一般地，将格林第三公式写成 为了求解场内 ，现在的问题就是要寻找既满足Laplace方程又满足给定边界条件的Green函数。上述是混合分布模型，解算存在困难，设法建立单一分布模型（如只有源分布）。 Green Function * （4）赫斯-斯密斯方法 ——面源分布方法 设 和 为外部域和内部域（虚构），分界面为S，法线均指向域外。按格林第三公式，对 内 和 内 ，分别有 两式相减，有 * （4）赫斯-斯密斯方法 ——面源分布方法 因在S上有 , 并记 。于是，上式成为 上式左端为含参数积分，相当于在S上有密度为 的源分布在场内某点 p 引起的速度势（单层势），这样的速度势在S上连续 ，其法向导数则不连续 ，即法向导数的阶跃构成了源强的分布密度 。 * （4）赫斯-斯密斯方法 ——面源分布方法 上式对内、外域及S 均成立，式中的 由边界条件确定，而 而 上的积分即是由其上的分布源诱导的法向速度，由 ， 即当 p 趋于 时，有 于是 * （4）赫斯-斯密斯方法 ——Hess-Smith方法 现讨论无界绕流问题。设远处来流速度为 ，在固定于物体的直角坐标系下，流场扰动速度势 满足 （物体外） （物面S上） （无穷远处） 这里的单位法线指向物体内（即域外）。在物面上布置源汇，由物面条件，有 上式的右端项为已知项，构成了关于分布源密度 的线性积分方程。解此方程可获得物面上的源密度 的分布，继而由下式确定流场内的速度势和诱导速度 这就是Hess-Smith方法的基本思想。 * （4）赫斯-斯密斯方法 ——Hess-Smith方法 将物面S 离散成N个单元小块 ，只要单元足够小，可认为每个单元上的源密度 为一常值，则有 并让物面条件在每个单元的某一点处（控制点：控制物面条件使之成立的点）成立，于是，积分方程离散成线性代数方程组： 这里， 称 为影响系数，即第j 个单元上的分布源在第i 个控制点上的影响。上述数值离散求解方法称为基于分布源模型的面元法，又称Hess-Smith法。 i :场点 j :源点 * （4）赫斯-斯密斯方法 —Hess-Smith方法 对于无界绕流问题，物面S外的流场内一点p 处的速度势 为 同理，令点 p 在流场中沿物面上点 处法线方向趋于 ，得到 上式既有分布源，又有分布偶极，称为混合分布模型。离散后，有 其中， * （4）赫斯-斯密斯方法 ——诱导速度公式 总结上述涉及到的单元曲面 上的积分有以下三类： 具体计算时，可用平面四边形或三角形（三角形也可以看成四边形）来近似代替小曲面 。具体处理如下： ①取其四个顶点坐标之算术平均值为中心 点 的坐标，该点视为控制点； ②计算对角线向量的向量积并取其单位向量 作为该单元的单位法线 ； ③ 过控制点 作垂直于 的平面，并以四个 顶点向该平面投影，得四个投影点构成平面四边形（或三角形） 。 pj * （4）赫斯-斯密斯方法 ——诱导速度公式 直接给出 上常值分布源密度为1时在 p x,y,z 点的诱导速度计算公式： 其中 * （4）赫斯-斯密斯方法 ——诱导速度公式 顺便，直接给出以下两组的计算公式： ☆ ☆ 单层位势 分布偶极诱导速度 * （4）赫斯-斯密斯方法 ——模型比较 分布源模型 源-偶混合分布模型 单元几何：三角形代替四边形。 分布源：高阶分布源代替均匀分布源。 边界面元上的 ，精度高 流场速度 改善途径： * 在t 时刻，流场中某一标量函数 在空间某点P 处沿 的方向导数： 表示沿该方向的变化率。而 在P点的最大方向导数则称为 的梯度（gradient），以 或grad 表示之。符号“ ”为Hamilton算符： 梯度是矢量 梯度 在直角系下的运算关系为 在流体微团界面法线上的投影即为法向导数，即 。如果 是流场的速度势，则微团界面的运动速度 满足： （1）场论基础 称之为物面条件 * 如流体以速度 流出或流进微团，则单位时间内流体进出的通量为 浮体水动力学中经常用到下面的关系式： 当微团趋于一点时，通量的变化率 称为散度（divergence ），表示流体微团或流体质点处的流体扩散率。显然，如微团内无源汇或流体进出平衡，则散度为零。 散度 也可以 表