一种基于修正流体运动方程线性迭代的稳态求解方法技术

技术编号:10813821 阅读:95 留言:0更新日期:2014-12-24 18:29
一种基于修正流体运动方程线性迭代的稳态求解方法,属于优化方法技术领域。首先采用有限体积法,以经典柯西运动方程为基础,在速度方程中引入对流速度的概念解构方程中原有的非线性对流项,推导出线性的修正流体运动方程,进而在有限体积离散方法环境下,建立一个专属于修正流体运动方程的稳态数值迭代求解器。在对修正流体运动方程对流项进行数值离散后,该求解器的线性特征充分显示出了其在处理流场问题中的优势,即所得到的方程组为一线性方程组,其中速度项系数a不再参与迭代过程,而且可以直接通过对流速度w解出,因此数值计算过程变成了简单的速度项U与压强项p的迭代求解过程。本发明专利技术具有运算效率高,计算时间短,运行成本低的特点。

【技术实现步骤摘要】
【专利摘要】,属于优化方法
。首先采用有限体积法,以经典柯西运动方程为基础,在速度方程中引入对流速度的概念解构方程中原有的非线性对流项,推导出线性的修正流体运动方程,进而在有限体积离散方法环境下,建立一个专属于修正流体运动方程的稳态数值迭代求解器。在对修正流体运动方程对流项进行数值离散后,该求解器的线性特征充分显示出了其在处理流场问题中的优势,即所得到的方程组为一线性方程组,其中速度项系数a不再参与迭代过程,而且可以直接通过对流速度w解出,因此数值计算过程变成了简单的速度项U与压强项p的迭代求解过程。本专利技术具有运算效率高,计算时间短,运行成本低的特点。【专利说明】
本专利技术涉及,属于优化方法

技术介绍
由于流体机械等工业研究领域的迫切需求,尤其是近30多年来计算机技术的迅 速发展,使计算流体力学(CFD)取得了很大的成就。今天,以数值求解N-S方程为代表 的CFD技术已经被广泛应用到流体机械的研发与设计过程中,取得了令人瞩目的成就。但 是,随着技术发展的逐渐深入,CFD也面临着越来越多的困惑,尤其是在N-S方程的求解问 题上。这一偏微分方程在经有限体积法离散后,往往会生成一个庞杂的非线性代数方程组。 从数学角度上来说,这一方程组的处理,不仅增加了求解CFD问题的难度,同时也严重阻碍 了相关的数值模拟计算方法的研究与发展。另外,N-S方程也对网格的质量有着苛刻的要 求,过高质量的网格无疑会严重增加计算机的负荷。这些都导致了 CFD模拟计算的超高耗 时,从而影响了流体机械的设计效率,增加了其研发成本。 经典N-S流体运动方程可以描述为:U 【权利要求】1. ,其特征在于,包括以下步 骤: (1) 通过计算机构建流体机械设计的物理结构模型,设定对流速度We、运动粘度V0、 流体密度Pe的初始化参数,计算出求解点速度项系数aP和求解邻近点速度项系数aN; (2) 根据流体机械设计要求设定时间步长,按设定时间步长进行循环迭代计算: d , Vpe 根据修正流体运动方程:Ιυ〃+,·νυ〃-ν/?νU〃=-I 式中:Ue为求解速度,pe为求解的压力; 进行时间步长的迭代计算过程,首先计算矩阵H(U) = '其中Un代表P 点邻近处各点的速度; 然后计算压力P:相对应的P点处的动量方程为:心Up =If(U)-V/;、速度方程为: 及、压力方程为:▽I丄=P ; Clp Ctp \apJ f\ aPJf 对压力P和速度Up依次进行求解,其中:压力P即为方程中的求解压力Pp速度Up即 为方程中的求解速度Ue ; (3) 最后计算速度UP,判断速度Up是否收敛, 当循环迭代计算数值的迭代时间步长-待求解项坐标曲线图的待求解项数值趋于稳 定时,则判定速度Up满足收敛条件,运算结束,提取运算结果为流体机械设计提供具体设计 参数; 当循环迭代计算数值的迭代时间步长-待求解项坐标曲线图仍有变化不稳定时,则判 定速度Up不满足收敛条件,则重复步骤(2),再进行时间步长的迭代计算,直至速度Up满足 收敛条件为止。2. 根据权利要求1所述的,其特 征在于:所述的设定时间步长根据不同流体机械的实际需求确定。【文档编号】G06F17/10GK104239625SQ201410454014【公开日】2014年12月24日 申请日期:2014年9月5日 优先权日:2014年9月5日 【专利技术者】万波, 张昕喆, 刘聪, 卞庆飞, 何炆峰, 彭懿延 申请人:中国矿业大学本文档来自技高网...

【技术保护点】
一种基于修正流体运动方程线性迭代的稳态求解方法,其特征在于,包括以下步骤:(1)通过计算机构建流体机械设计的物理结构模型,设定对流速度wβ、运动粘度νβ、流体密度ρβ的初始化参数,计算出求解点速度项系数aP和求解邻近点速度项系数aN;(2)根据流体机械设计要求设定时间步长,按设定时间步长进行循环迭代计算:根据修正流体运动方程:∂∂tUβ+wβ·▿Uβ-vβ▿2Uβ=-▿pβρβ,]]>式中:Uβ为求解速度,pβ为求解的压力;进行时间步长的迭代计算过程,首先计算矩阵H(U):其中UN代表P点邻近处各点的速度;然后计算压力p:相对应的P点处的动量方程为:速度方程为:UP=H(U)aP-▿paP,]]>压力方程为:▿·(1aP▿p)=ΣfA·(H(U)aP)f;]]>对压力p和速度UP依次进行求解,其中:压力p即为方程中的求解压力pβ,速度UP即为方程中的求解速度Uβ;(3)最后计算速度UP,判断速度UP是否收敛,当循环迭代计算数值的迭代时间步长‑待求解项坐标曲线图的待求解项数值趋于稳定时,则判定速度UP满足收敛条件,运算结束,提取运算结果为流体机械设计提供具体设计参数;当循环迭代计算数值的迭代时间步长‑待求解项坐标曲线图仍有变化不稳定时,则判定速度UP不满足收敛条件,则重复步骤(2),再进行时间步长的迭代计算,直至速度UP满足收敛条件为止。...

【技术特征摘要】

【专利技术属性】
技术研发人员:万波张昕喆刘聪卞庆飞何炆峰彭懿延
申请(专利权)人:中国矿业大学
类型:发明
国别省市:江苏;32

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1