OpenFOAM中动量方程 UEqn 浅析

OpenFOAM中大多数求解器（不可压缩流）采用了以下的动量方程：

fvVectorMatrix UEqn
(fvm::ddt(U) + fvm::div(phi, U)+ MRF.DDt(U)    //MRF是多参考坐标系，在动网格中会用到（像旋转机械等），一般情况下可忽略+ turbulence->divDevReff(U)==fvOptions(U)//源项或约束
);

fvm::ddt(U) ：时间导数项；
fvm::div(phi, U)：对流项，涉及到非线性问题的线性化处理；
MRF.DDt(U)：涉及动网格技术（旋转机械等）需要引用MRF，否则可忽略。具体可参考MRF.DDt(U)—MRF模型简介；
turbulence->divDevReff(U)：剪切力张量项，可查看turbulence->divDevReff(U)代码解读；
fvOptions(U)：源项或约束，如：多孔介质、磁流变液和体积力等。具体可参考fvOptions(U)；
在OF中，压力梯度项一般不直接放入UEqn中，而是根据算法字典中的关键字momentumPredictor的状态，对动量方程进行修正。

if (pimple.momentumPredictor())//动量预测开关，默认为关闭，这里的p由上一时间步求得
{solve(UEqn == -fvc::grad(p));fvOptions.correct(U);
}

OpenFOAM中动量方程 UEqn 浅析相关推荐

OpenFOAM中参考压力p_rgh的由来
在OpenFOAM的动量方程UEqn.H中经常能看到以下代码: solve (UEqn==fvc::reconstruct((- ghf*fvc::snGrad(rho)- fvc::snGrad(p ...
动量方程基本表达形式
不同文献中动量方程的表达形式也多种多样,现总结如下: 守恒形式: ∂(ρU)∂t=−∇⋅(ρUU)−∇⋅τ−∇p+ρg\frac{\partial( \rho \mathbf U)}{\partial ...
openFOAM动量方程中张量的处理
动量方程有许多种形式.本文以rhoPimpleFoam为例子进行说明. rhoPimpleFoam/UEqn.h内容如下: 文件内容 tmp<fvVectorMatrix> UEqn (f ...
有限体积法求解二维方腔流（一）——动量方程和连续性方程的离散
1. 有限体积法求解二维方腔流–理论手册 1.1. 不可压缩流体控制方程连续性方程 ∇⋅U=0(1)\nabla \cdot U=0 \tag{1} ∇⋅U=0(1) 动量方程 ∂U∂t+∇⋅(UU ...
OpenFOAM 中边界条件的设定【转载】
转载自:http://blog.sina.com.cn/s/blog_a0b4201d0102v7jt.html 用习惯了FLUENT的操作界面,再使用OpenFOAM就会觉得非常繁琐.遇到的第一个问 ...
流体力学-----动量方程的推导和理解---NS方程
文章目录前言推导动量方程的流动模型推导过程书中给的剪切力分析前提条件:速度的三个分量u.v.w的正增量和坐标轴一致前言可以参考之前的博客计算流体力学1-流体力学的控制方程推导动量方程的流 ...
工程流体力学笔记暂记7(动量方程)
积分形式的动量方程其为动量定理的表达形式对于一般刚体有高中所学的牛顿第二定律但对于流体还有流入流出流体所携带的动量定常流动的性质不随时间变化 (注意:牛顿第二定律仅适用于惯性坐标系) 对于定常 ...
连续方程能量方程动量方程_简单的可伸缩性方程
连续方程能量方程动量方程排队论排队理论使我们能够预测队列长度和等待时间,这对于容量规划至关重要. 对于架构师来说,这是一个非常方便的工具,因为队列不仅仅是消息传递系统的工具. 为了避免系统过载,我 ...
计算流体力学三大基本方程：连续性方程、动量方程（纳维-斯托克斯方程）和能量方程推导
方法取自安德森<计算流体力学基础> 连续性方程固定流体微元内质量变化率=流体从笛卡尔坐标三个方向流出量因此可得: 质量变化率: 则: 连续性方程: 用散度表示则可得到: 对于不可压缩流 ...

OpenFOAM中动量方程 UEqn 浅析

OpenFOAM中动量方程 UEqn 浅析相关推荐

最新文章

热门文章