找不到资源
您好，欢迎来到欧浦钢网 [] []
综合频道：
意向订单(0)
您请求的页面无法找到
，10秒后自动返回首页
全国服务热线：400-887-8388 对本站有任何建议、意见或投诉，请
Copyright &
Europol Corporation, All Rights Reserved.
佛山公安网站备案编号:09112____求解N―S方程的压力修正
上亿文档资料，等你来发现
9112____求解N―S方程的压力修正
§12；求解N―S方程的压力修正方法；1、压力修正的意义：对动量方程的离散形式进行求；们迭代的压力场P*计算出来的速度场，未必可以满足；方程；是给定的或是上一层次计算出来的；度场未必满足连续性方程，所以要对压力场进行修正；4、压力修正法步骤：；1）假设一个压力场p*；aeu??au?b??p?p?Ae；*e*nbe*P*E；3）改进压力场，要求与改进后的
§12求解N―S方程的压力 修正方法1、压力修正的意义：对动量方程的离散形式进行求解时，我们迭代的压力场P*计算出来的速度场，未必可以满足连续方程。所以我们要寻找一个改进压力估计值P*的方法，以 使所得的带*的速度场逐渐的满足连续方程。 2、压力修正方法的基本思想：在对N―S方程的离散形式进 行迭代求解的任一层次上，可以给定一个压力场，它可以是给定的或是上一层次计算出来的。但由它所计算出的速度场未必满足连续性方程，所以要对压力场进行修正。 3、修正原则：与改进后的压力场相对应的速度场能满足这一 迭代层次上的连续性方程，据此来导出压力与速度修正值， 并以修正后的压力与速度开始下一层次的迭代计算。4、压力修正法步骤：1）假设一个压力场p*。 2）利用p*，求解动量离散方程，得出相应的u*,v*(不完善的速度场)。ae u ? ? a u ? b ? ? p ? p ?Ae* e * nb e * P * E3）改进压力场，要求与改进后的压力场相对应的速度场能满足 连续性方程。即用p’,u’,v’分别表示压力与速度的修正量，称之为压力修正值（pressure correction）与速度修正量 (velocity correction)。4）以正确的压力（p*+ p’），及相应的（u*+ u’），（v*+ v’）作为本层次的解，并据此开始下一层次的迭代计算。引出两个关键问题：1、如何获得压力修正值p’ 以使与p相应的u、v满足连续性方程?2、 获得压力修正值p’后，如何确定u*,v*?12.1 速度修正值的计算公式改进后的压力场和速度场满足这一迭代层次上的动量离散方程：* ' * ' * ' * ' ? ae ? ue ? ue ? a u ? u ? b ? p ? p ? p ? p ? ? nb ? nb nb ? ?? p p ? ? E E ?? ? Ae因为u*e是据p*从这一动量离散方程中导出的，因此也满足 动量离散方程，虽可能不满足连续性方程。ae u ? ? a u ? b ? ? p ? p ?Ae* e * nb e * P * E两式相减得：a u ? ?a u ?? p ? p' e e ' nb nb ' p' E?Ae说明：任意点速度的改进值u’e由两部分组成：1) 该速度在同一方向上的相邻两节点间压力修正值之差――产 生速度修正值的直接动力。 2) 由邻点速度修正值所引起的，可视为是四周压力的修正值对 所讨论位置上的速度改进的间接影响。为了简化，忽略间接影响，相当于假设∑ anbu’e 中anb=0的。 得到速度修正方程：' ' aeue ? ? p 'p ? pE ? Ae或? Ae ? ' ' ' ' u ? ? ? ? p p ? pE ? d p ? p ? e? p E? ? ae ?' e? An ? ' ' ' ' 类似可得 v ? ? ? ? p p ? pE ? ? d n ? p p ? p N ? ? an ?' n于是改进后的速度为：ue ? u ? de ? p ? p* e ' p' Evn ? v ? d n ? p ? p* n ' p' N? ?(1)12.2 求解压力修正值的代数方程压力修正值p’应满足的条件：由p’而改进的速度场u、v能满足连续性方程。将式1代入连续性方程的离散格式即可获得p’的代数方程。?? ? ? ? u ? ? ? ? v ? ? ? ?0 ?t ?x ?y在时间间隔Δt内对主控制体做积分，得：0 ?p ? ?p?t? x? y ? ? ?? ? u ? e ? ? ? u ? w ? ?? y ? ? ?? ? v ? n ? ? ? v ? s ? ?? x ? 0将（1）代入并整理可得：' ' ' ' a p p 'p ? aE pE ? aW pW ? aN pN ? aS pS ?b（2）其中：aE ? ?e d e ?y，aW ? ? w d w ?y aN ? ? n d n ?x，aS ? ? s d s ?xaP ? aE ? aW ? aN ? aS0 （? P -? P）?x?y * * ? ? b? ?（ ? u ） ? （ ? u ） ?y w e ? ? ?t * * ? ? ?（ ? v ） ? （ ? v ） s n ? ? ?x--------压力修正值的代数方程。说明：1）b项是按带星号的速度取值的离散化连续性方程左侧的负值。如b=0，说明该速度场已满足连续性方程，不必对压力进一步进行 修正，p’方程迭代已收敛，即由该压力场求解的速度场定满足二 守恒方程； 2）可用各控制容积b值（不满足连续性方程的剩余质量的大小）作为判断速度场收敛的判据。收敛速度与p’的特定方程有关，如忽略项过多，会影响收敛速度。如b≠0，重新校正压力场p’后，再求b值，直至 b≈0时， p’ 压力场可使速度分布满足连续性方程； 3）由p’计算而得的u’，v’可使u、v能满足连续性方程，以u，v作 为本层次上的速度场的解，并用它去改进动量离散化方程系数，开始进行下一层次的迭代计算；4）因求解的问题是非线性的，没有必要在每一层次迭代上去获得p’真解，采用ADI方式对p’场作3-5次扫描后即可 停止迭代（初始令各点p’=0）。12.3 压力修正值方程的边界条件工程流场计算中边界条件一般有两种。 1) 边界上的压力已知：则p’=0 2) 边界上的法向速度已知，设ue已知，则u’e=0，则不必引入p’E。边界法向速度ue已知:* ' ' / u e ? ue ? ue u ? de ? p aE ? ?e d e Ae ' / e e ? ue ? 0 ? de?pe ? 0 ? d e ? 0 ? aE ? 0 * ue 用已知值边界压力已知:* / pE ? pE ? pE 在压力修正值方程中 相当于 / / ?? pE ? 0 ? aE p E ? 0 ? aE ? 0 * / pE 取已知值 以aE p E 形式出现说明：不论是边界压力已知，还是法向速度已知，都没有必要 引入关于边界上压力修正值的信息。在计算中，可令与边界相 邻的主控制容积的p’方程相应的影响系数为零。12.4 SIMPLE算法（Semi-Implicit Method forPressuer Linked ）（解压力耦合方程的半隐方法） SIMPLE法就是一种计算压力修正的一种方法。1、半隐含义：所谓半隐式是在公式中略去了某一项的处理方法。略去间接影响保留直接影响的称为半隐式,全保留是全隐式/ / aeue/ ? ? anbunb ? ( p /p ? pE ) Ae1）后部分是压力修正对u’e的直接影响；前一部分是压力修正对u’e的间接影响，忽略该项称为半隐；2）计算中若保留∑ anbu’nb项，则u’e方程就是全隐的代数方 程。势必存在u’e的相邻点的压力修正和速度修正，相邻 点又引入新的相邻点，这样使计算变复杂，无法进行逐 点计算。但此项的忽略并不影响最终收敛解。 2、simple算法的计算步骤 1）假定一个速度分布，记为u0，v0，以此计算动量方程中的系 数及常数项。 （源项中压力梯度无法求）2）假定一个压力场p*；3）依次求解两个动量方程，得u*,v*（初始场）； 4）求解压力修正值方程，得p’；（ u*满足动量守恒但不满足连续性方程） 5）据p’改进速度值；（以使新速度u满足连续性方程但可能不满足动量守恒） 6）利用改进后的速度场（u→u0，v→v0）重新计算动量离散方程 的系数及源项中有影响速度场耦合的Φ变量，如果Φ并不影响 流场，应在速度场收敛后再求解。7）利用改进后的压力场（p→p*）作为下一次迭代计算的初值。重复上述步骤，直到得到收敛解（即满足质量方程b→0，又满足动量方程ε→0）。分析讨论：1）为什么在获得压力的改进值p（p*+ p’）后，不直接利用这一改进值及u*,v*去开始下一层次的迭代，而还要计算这一层次的 u’，v’？ 因u*,v*还不满足这一层次的连续性方程，如果用它们去确 定新的系数，开始下一层次迭代，会影响迭代的收敛速度，也 会使ap = ∑ anb 的关系得不到保证。（注（Fe- Fw）项，会使代数方程组系数矩阵对角占优的条件遭到破坏。）2）如何获得p’的离散化方程？盲目试探压力场p*是不明智的，如能找到一种调节压力方法，即压力修正p’方程，压力校正方程仅是起引向正确压力场的一个中间算法，作用相当于一 个负反馈器，对最终结果无直接影响，如推导p’方程略去项太多，可能导致 发散，起正反馈作用。12.5 simple算法的讨论1、采取的简化假定 1）速度场u0,v0 及p*的假定是各自独立进行的，两者间没有任 何联系。但实际上如果给定一个速度场，压力也随之而定。 独立设定可能会使两场不匹配。2）略去∑ anbu’nb项意味着在对速度修正值计算时没有把邻点速度修正值的影响考虑在内，如果对等号前的项作相应的变 化，可以减少因为略去邻点修正值的影响。 3）采用线性化的动量离散方程，也就是在一个层次的计算中， 动量离散方程中的系数及源项假定为定值。实际上它们与速度有关。2、上述简化只会影响收敛的快慢，不会影响收敛的解simple算法是以迭代方式进行求解的过程，它必须满足一个基本要求，即迭代模式的组织不能影响最终的收敛解。Simple算 法中采用的简化处理都满足这一基本要求。 1）u0,v0与p*之间独立的假设，可能造成速度场与压力场之间的 不协调，但只要迭代过程是收敛的，随着迭代的进行，这种不协调会逐渐减小，以致消失。2）略去∑ anbu’nb，∑ anbv’nb对最终的收敛解也没有影响。因为如 果迭代收敛，则u’→0,因而∑ anbu’nb，∑ anbv’nb均→0， 3）当迭代趋近于收敛时，两层间u，v也趋于相等。因而离散方 程的系数及源项保持不变。12.6 速度与压力修正值的亚松弛速度修正时略去邻界点的影响，使压力修正本身被夸大了，从而使叠代过程中速度场与压力场不能同步，需采用亚松弛。 压力修正：p ? p ??p p*'ap≈0.8为压力亚松弛因子为限制两层之间的变化，以有利于非线性问题的迭代收敛，将本 层次的计算结果与上一层次结果的差值做适当缩减，以避免由于差值过大而引起非线性迭代过程发散，需采用亚松弛。速度修正：? ? anb?nb ? b ? 0 0? ?p ? ?p ? ? ? nb ? ?p ? ? ap ? ? ? ?Φp0为上一层次之解，a为松弛因子式子改为：（ap?）?p ? ? anb?nb ? b ? (1 ? ? )nbap??0 pap主对角线的系数是?而不是ap，源项改为 b ? (1 ? ? )ap?0 ?p此时代数方程的解就是亚松弛后的解。12.7 流场迭代收敛的判据作用:在算法中,对αE、αW、αN、αS和b假定为定值。但实际受流速的变化而变化。因而没有必要把一组临时的系数与源项的代数方程准确地求出来，采用迭代法可以适时地停止迭 代，及时地用所得到的解去更新系数与源项，以进入下一层 次的计算。 1、两种迭代收敛问题1）同一层次上代数方程迭代求解――内迭代2）非线性问题，从一层次向另一层次推进的问题――外迭代2、终止内迭代的判据一般有三种方法1) 简单的规定迭代的轮数. 2) p’方程余量的范数(k ) Rp小于某一数值.1 2(k ) Rp ? {?[(? E p E '?? w pW '?? N p N '?? S p S '?b ? ? P p P ' ) ( k ) ]2 }R(k ) p? ?P(k ) RP 3) 终止迭代时的范数与初始范数之比小于允许值. ??P (0) RP注:①γp成为余量下降率.其值一般取0.25~0.05. ②迭代次数从始至终近似相同. ③缺点:增加了计算范数的工作量.说明：在内迭代中p’方程求解是关键，在每一层次的迭代上：1）若p’迭代停止的太早，则速度修正值u,v不能较好地满足连续性方程，由于u,v要用于下一层次迭代方程系数的确定，于是，误差就会传播出去，以至于使迭代发散， 2）若p’终止迭代的太晚，将不经济。即占计算机内存和求解 时间，因为p’方程的求解时间占到整个求解时间的80%以上。 3） p’方程组相当于一个扩散方程的离散形式，边界条件处理 相当于绝热，一个绝热体系只有当源项总合为零时，才能维 持稳定的场，这相当于要求计算区域满足总体质量守恒。3、终止非线性问题迭代的判据1）特征量在连续若干个层次迭代中的相对偏差小于允许值。(如Num平均努谢尔特数)(k ?n) (k ) Nu m ? Nu m ?? (k ?n) Nu m2）内节点上连续方程余量的代数和(Rsum)及余量的最大绝对值(Rmax)小于一定数值。(考虑有流量作参考)Rsum ?? qm且Rmax ?? qm3）连续性方程余量范数的相对值小于允许值。? (b ) ?? qmp 24）要求在整个求解区域内动量方程余量之和的范数与参考 动量之比小于一定值(?{a e u e ? [? a nb u nb ? b ? Ae ( p P ? p E )]} )nb1 2 2?u2 in??注:为确保可靠性,可以几个条件同时使用， 相对偏差ε一般取10-3~10-5思考题1、压力修正的基本思想，步骤，如何用p*求 解u*,v* ？ 2、压力修正速度修正如何得到？ 3、为什么求压力修正量，b为什么可以判断 收敛，b大说明什么？ 4、SIMPLE步骤中的Φ指的是什么？改进后 的值有哪些？ 5、 ∑ anbu’nb为什么与收敛解无关？怎么理 解收敛的时候u’，v’在收敛时相等？包含各类专业文献、生活休闲娱乐、中学教育、外语学习资料、文学作品欣赏、应用写作文书、高等教育、行业资料、9112____求解N―S方程的压力修正等内容。
　流体计算理论基础_物理_自然科学_专业资料。主要介绍...0 ,动量守恒,简称动量方程,也称 N-S 方程。 关于...目前工程上最为广泛的流场数值解法是压力修正法, ... 　s ? 0.44kg / m3 , g ? 9.8N / kg , ...9.8 ? 0.12 ? 16005.24 压力表的读数 pg ? ...p水 g (h5 ? h4 ) 联立方程,代入数据解得 p ... 　若导线以速率 50m/s 拉过模具, 试求所需牵拉力。 (1.O1N) [解] ? A...方程中两边的压强性质要相同,由于 Ps 为绝对压 强,因此出口处也为绝对压力。 ... 　3、湍流控制方程 三维的 N-S 方程无论对于层流还是湍流都是是使用的, 但由于...常有的原始变量法包含的常见解法有压力泊松方程法, 人为压缩法和压力修正 法。... 　压力损失计算 采用恒定总流能量方程式 P1/ν+Z1+...常温下:空气 ν:12.02N/m3,9.8 kN/m3。 P1/...m/s a1:介质的速度修正系数。 h’:压力损失,m。... 　必须使用严格计算法,以便确定各级上的温度、压力、流...(n+2)个方程,c 为组分数,简称 M 方程 组。 ...分率加和式, 共 3n+4 个方程, 简称 S 方程组... 　(u ) 、压力 ( p ) 和单位质量内能 (e) 的...空间步长 -B.12- xr=1.0 rou0=1.0 u0=1.0...N-S方程 19页 1下载券 一种求解N-S方程的自适应... 　6、直接求解每个网格压力波动的泊松方程: ? 需要注意的是矢量 F = F n ??...7、应用 Runge-Kutta 的第二部分更新速度: 再次注意,修正的速度场的散度和前面... 　由数值解出所有的 N-S 方程来模拟; b)解出尺寸...12%的情况 VOF 法用来处理没穿插的多相流问题 ...求解压力修正方程,得 p’ ; 5、据 p’改进速度...为了没有返程票的人生旅途——镇江市京口区政府实施“出生缺陷干预工程”和“特殊..
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
为了没有返程票的人生旅途——镇江市京口区政府实施“出生缺陷干预工程”和“特殊教育”纪实 [可编辑]
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口