特征线法及其在流固耦合系统中的应用.pdfVIP

特征线法及其在流固耦合系统中的应用+ 陈责清董保珠 唐山学院基础教学部，河北唐山063000 摘要：特征线法以偏微分方程的特征理论为基础，是求解偏微分方程组的 一种近似计算方法。本文对该方法进行了简介，并介绍了其应用在流固耦合 系统中而推导出的几种典型特征方程，晟后按自然进程总结并分析了其在工 程中的应用。 关键词：特征线法，水击，管道 特征线法的基本原理是将偏微分方程组化为特殊的全微分方 程组，即特征方程，然后再转变为一阶有限差分方程，求其近似 解。 特征线法具有独特的优点，可处理：(1)流固耦合系统和边 界条件： (2)层流和紊流摩擦损失以及局部损失：(3)汽穴和 液柱分离等复杂物理过程。它的理论严密，物理意义明确，容易 编程，且具有足够的计算精度，从而成为进行数字仿真的一种有 效工具。 1 特征线法方法简介 对于含有两个自变量z，Y和两个因变量甜，1，的拟线性偏微 分方程组，其一般形式为… ★ 联合资助 陟贵清(19“一)，男，教授。博士，研究方向：结构流固耦合振动研究 扣iqmgch∞@sma．com 130 厶=4ux+蜀％+Cl叱+Dl匕+巨 (1) L2=／hu，+B2uy+C2V，+D2W+易(2) 且条件 其中4，4，…，E，易是x，y，U，V的函数， 线性组合 三=五厶+鸽=以4+见4n+似墨+丑垦心+研C：+是G也 (3) 、7 +(a日+友皿)K+(a乓+五岛) 选择常数A，五，使得 查：刍兰z±垒生：刍刍±垒刍 一=——L—二————=—二=——L—_=————_=-_： 砂 AE+五色丑q+五B 则微分表达式三可以写为 (5) 这时在微分表达式三中，U，1，的导数是如此组合的以至他 们的导数取在同一方向上，即砂／dx上，这个方向称为特征方向。 由(4)式可以得到 (6) a(dy)2—2bdxdy+c(dx)2=0 Cl，2b=4D2—4Dl，C=Bi砬一B2Dl，为了 其中口=4C2—4 131 方便起见，假定口≠o，可引进斜率孝：ey．’，并且方程 dx 口f2—2b孝+c=0有两个不同的实数解六和乒。那么，在(x，Y) 点上，由掣=六，孕=乒给出两个不同的特征方向，且六和乒 aX反X 是关于x，Y，z，，V的函数。一旦求得(1)和(2)的解u(x，y)， v(x，Y)，则六和乒是关于x，y的函数，这确定x，Y平面上 属于解u(x，少)和v(x，y)的两个单参数变量的特征线组C+和 Q。这两族特征线形成一个曲线坐标网络。可得 Fdu+(口六一G)幽+(K鼻-H)dx=o(7) (8) Fdu+(口乒一G)咖+(K乒-H)dx=0 此处 ∽ ，K兹二箍2GH嚣篙=4易一4日，=骂岛一垦巨 一 于