OpenFOAM里fixedJump命令

文章介绍了在二维泊肃叶流动中,如何通过OpenFOAM的fixedJump命令在周期性边界条件下创建压力差来驱动流体流动。cyclic和cyclicAMI两种边界条件被提及,其中cyclic要求网格匹配,而cyclicAMI则不需要。通过设置Jump量,使得入口压力大于出口压力,从而模拟压力差驱动的流动。paraview的分析结果显示压力差与设定值相符。

问题描述

对于二维泊肃叶流动来说,驱动流动一般有两种方式,一种是体积力驱动,另一种是进出口压力差驱动。这里我们主要对对第二种方式进行说明。

周期性边界条件

对于channel flow等流动来说,在进口和出口出我们一般采取周期性边界条件,在OpenFOAM种有cyclic和cyclicAMI两个命令,其中cyclic要求进口和出口网格完全一致,而cyclicAMI则没有这个要求。
在这里,周期性边界条件的意义是进口和出口面上物理量完全一致,简单的来说,即出口流出的继续从入口流入,采用这种边界条件的意义主要是用来模拟无限长的通道。这样的话会产生一个问题,对于压力p来说,进口和出口完全一致的话,就没有压力差了,那我们如何用压力差进行驱动流体流动呢?

利用fixedJump命令设置压力差

事实上,fixedJump命令就是解决这个问题的,我们可以设置一个Jump量,使得入口处的压力等于出口处的压力加上一个常数,这样就可以产生压力差。其具体设置在0文件夹下的p文件中。

    inlet
    {
   
   
    type            fixedJump;
	patchType       cyclic;
	jump            uniform -0.001220377373900;
	value           $internalField;
    }
    outlet
    {
   
   
    type            fixedJump;
	patchType       cyclic;
	jump            uniform 0;
	value           $internalField;
    }

上面是一种设置方法,这里我们需要注意两点:
1、我们想要的结果是inlet处压力更大,因此需要设置一个Jump量,特别需要注意这里是负号
2、对于出口处,我们的Jump量可以设置为0

验证

在paraview中打开已经算好的算例,只打开inlet面,并利用工具探测处inlet面上的p的上下界;如下图
在这里插入图片描述

pinletmax=0.0245446;pinletmin=−0.00250727 p_{inletmax} = 0.0245446;\quad p_{inletmin} = -0.00250727 p

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值