OF1.7中的p_rgh【翻译】
翻译自:CFD-online
帖子地址:http://www.cfd-online.com/Forums/openfoam-solving/80454-p_rgh-1-7-a.html
stawrogin:
亲爱的Foamer们
我对于OpenFOAM当中新的压力p_rgh有点困惑。它是怎样定义的?我有点担心的设置好的边界条件。比如如果我想要在出口就设置一个常数的"压力",这应该如何定义?如果定义p_rgh = 0,是否意味着在这个算例中pressure + (rho)*g*z是一个固定值?
我也很惊奇为什么在一些求解器中同时存在p—文件和p_rgh文件。P—文件仅仅用在后处理吗或者它也有可能在这里定义了边界条件?
先谢过
Stawrogin
juho:
在1.7.x的buoyantPimpleFoam求解器中,p和p_rgh文件都被读取。P被用在热力学模型中。
在createFields中,p_rgh被定义为:
p_rgh = p - rho*gh;
因此这个压力没有静压并且被从p文件中的压力场开始初始化。
在pEqn.H头文件中,压力方程被p_rgh来写和求解,因此p_rgh边界条件对于压力求解很重要。在解算完压力后,p通过下式计算:
p = p_rgh + rho*gh;
stawrogin:
亲爱的Juho
非常感谢你的帮助,你的回答对于我来说很清楚了。
Stawrogin
deniggo:
你好,Juho
p_rgh似乎是动压。我使用pimpleFoam求解器和pisoFoam求解器。是否可以调整这些代码获得p_rgh?
谢谢
Nico
Angela Wang:
你好,Nico
p_rgh代表动压,是否意味着p_rgh=1/2*rho*U^2?
我检查了案例中damBreak的算例,为什么p_rgh的值比p的值大?这是否意味着rho*g*h是负数?
我非常困惑。
Quote:
Originally Posted by deniggo Hi Juho, Thanks, Nico |
kjmaki:
你好,Angela
p_rgh是p-rho*g*h或动压。它在1.5版本中被叫做pd,它在版本1.6中是为了求解总压,在版本1.7中仅仅是为了求解动压。
Kevin
deniggo:
你好,Kevin
但是为什么p_rgh比p高?一个负的静压(rho*g*h)怎么可能?
不应该是1/2*rho*U^2(动压定义)导出p_rgh?贴出计算后U—场的分布。
nico
nbadano:
亲爱的朋友们,
够古怪了,gh在OpenFOAM中被定义:
gh = g & mesh.C()
换句话说,gh是g向量和单元中心位置向量的点积。g常常被定义为(0 0 -9.81),gh的值常常正z是负值!
结果的导致:
p_rgh = p-rho * gh
如果z>0,p_rgh比p大
希望这个有帮助
祝好
Nicolas
Angela Wang:
我的问题解决了,p_rgh是动压= p- rho*g*h
Logan Page:
关于这个一个问题:
使用buoyantBoussinesqSimpleFoam求解器,p和p_rgh是(m^2)/(s^2)[比如流动压力]
很明显的是这样定义太简单了:(p_rgh/rho)=(p/rho)-g*h
我的问题是当通过划分,密度如何使用。是rhok = 1 - beta (T - Tref)??
谢谢
kerim:
朋友们好
我将要继续Nico的想法。对于我来说,p_rgh不是动压。我将层流溃坝算例的结束时间改为500。那时初始流后将固定。因此在t=500时没有动压。但是interfoam求解器没有给定一个非0的p_rgh压力。这就是为什么我认为p_rgh不是动压。谁能告诉我发生了什么吗?或者我哪里错了吗?
Kerim
nbadano
你好,Kerim
p_rgh不是真正的动压,特别是在整个计算域内密度变化的多相流。是压力场和rho*g*z场。我认为这也是这个场在1.5版本没有被叫为pd的原因。
为了增加困惑,记住rho*g*z也不是真的静压!!甚至密度是常数(在z=0和大气压力p=0平面之间乘以rho*g)。
这里有一个p_rgh, rgh和p作为静压条件的一个草图(完全没有移动)。希望有帮助。
Cristiano:
这里的朋友你们好
我正在模拟一个在大气条件下测试过的燃烧器并且我有点困惑。
我应该如何为p和p_rgh设置合适的边界条件?
Inlet -->
Outlet->
Wall--->
实在谢谢你
Cristiano
linch:
Quote:
@nbadano I agree to your sketch: if the p-field is smooth, the h-field is also smooth, the g-field is uniform and there is a jump in the rho-field at h=h2, so the must be a jump in thep_rgh field. But if you set up your example in interFoam you'll notice, that the p_rghfield is uniform also. Though the p_rgh definition might be a littlte different. Well, theoretically you'll get an uniform field, if you divide p_rgh by rho, but in this case you wouldn't get the p-dimension [kg*m^-1*s^-2]. So somehow the hydrostatic pressure is being substracted and not the simple product of rho*g*z. |
对不起,我错了。P_rgh不是均匀的。
christianfrias:
尽管这个建议很简单,但还是值得说。
你有一个成功运行OpenFOAM1.6版本的算例并且想要运行在OpenFOAM1.7平台。记住你仅仅需要重命名你的p文件为p_rgh(在你的transportProperties 和fvSchemes文件里面改变一些其他参数)。然后将g的值从(0 0 -9.81)改为(0 0 0)。这个算例应该可以和OpenFOAM1.6一样运行,如果他们能顺利运行在OpenFOAM1.6下,你也无需担心你的边界条件是否适用于OpenFOAM1.7。
祝好
Christian F
为了理解不同压力,请看伯努利:
动压 --> 1/2*rho*v^2
液压--> rho*g*h
静压 --> p
1/2*rho*v^2 + rho*g*h + p = Constant
从OpenFOAM的站点,p_rgh=p-rho*g*h。
因此,p_rgh是静压减去液压,在任意高度
我不知道静压是从这个讨论的什么地方出来的。我认为p_rgh没有任何意思。动压是运动流体的压力并且在流体速度为0时,它将和静压转换。能量守恒回到伯努利。
很高兴知道:
为了让我的VOF模型运行,我将fvSolution文件放置在PISO求解器下。
pRefPoint (0.0 0.0 0.0);
pRefValue 1e5;
成功!
deniggo:
大家好
谢谢King清楚的解释。
为了解开压力疑惑:后处理工具"ptot"计算了每个时间步的总压(静压+动压):
http://www.cfd-online.com/Forums/openfoam/74978-where-define-fields-write.html
祝好
Nico
fasfcastro:
我认为p_rgh是扰动压强(从一个流体力学平衡状态)实际被用于模拟动量方程。
Oke'e:
为什么在bouyantBoussinesqPimpleFoam求解器和bouyantBoussinesqSimpleFoam求解器算例中p_rgh的单位是m^2s^-2?
但是在(多相流)interFoam算例中的单位是kgm^-1s^-2。
Andrea_85:
因为在前两个算例中除以了密度。为了获得"真实"压力必须乘以密度[kg/m^3]。
祝好
andrea
MOHAMMAD67:
亲爱的朋友们
你们好
这里已经有了许多关于p和p_rgh的讨论。但是自己没有从这些讨论中得出最后的结果。我想知道在多相流中为了和实验对比,我应该选择什么压力?如你所知,我们在实验室里只能得到液压。
p_rhg的值是否依赖于初始的计算域?
christianfrias:
我看既然p_rgh不是动压也就意味着p不是全压(如我们所想)实际它是静压。为了计算全压,你可以使用ptot(它将会计算1/2*rho*U^2 [动压] + p [静压)])。因此,为了将OpenFOAM的结果与实验数据对比,我将使用OpenFOAM中代表静压的p文件。
sharonyue:
Quote:
Originally Posted by christianfrias I see that since p_rgh is not the dynamic pressure it means that p is not the total pressure (as we could think) and is actually the static pressure. To calculate the total pressure you can use ptot (this will calculate 1/2*rho*U^2 [dynamic pressure] + p [static pressure)]). So, to compare experimental data with the results from OpenFOAM I will use the p file which is the static pressure in OpenFOAM. |
因此p(静压)加上1/2*rho*U^2(动压)常常在实验中被测量?为什么我记得在一个方盒或者其他什么东西,静压被一个压力表测量。
vonboett:
Quote:
Originally Posted by sharonyue so is the p(static pressure) plus 1/2*rho*U^2(dynamic pressure) usually measured in experiment? Why I remember usually in a tube or something ,static pressure is measured in a pressure gauge |
依据你测量的什么。水压力水平传感器(和压强计)眼下等于液压关于传感器,因此水平液压传感器并且你需要和OF中的p进行对比。如果你测量的是流动方向阻力的增加等。这个就可能很难直接对比测量和模拟了。
baedmaen:
大家好
我现在对于多相流的压力很困惑。有人能帮我把定性压力特征和附件中的文件关联起来吗?
第一幅图:
静止的两相系统(水和空气);我认为nbadano已经在2010年12月3日通过帖子给出了答案。
第二幅图:
气泡上升的表面,截图
第三附图:
一滴水滴入水的表面
非常感谢
附件:
Dream:
我们设置了分离器模型并且模拟了多相流。在我们运行了setFields,下面的信息出现了:
Setting field default values
Setting internal values of volScalarField alphaair
Setting internal values of volScalarField alphawater
--> FOAM FATAL IO ERROR:
Cannot find patchField entry for wall
file: /home/dell/OpenFOAM/damBreak4phaseFinelmxlmx/0/alphawater.boundaryField from line 34 to line 56.
From function GeometricField<Type, PatchField, GeoMesh>::GeometricBoundaryField::readField(const DimensionedField<Type, GeoMesh>&, const dictionary&)
in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 206.
FOAM exiting
-----------------------------------alphawater-------------------------------------------
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphawater;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField
{
Wall
{
type zeroGradient;
//value uniform 0;
}
new_new_out
{
type fixedValue;//;
value uniform 1;
}
new_out
{
type fixedValue;
value uniform 1;
}
out
{
type fixedValue;//outletInlet;
value uniform 1;
}
in
{
type outletOutlet;
outletValue uniform 0;
value uniform 0;
}
}
可能我们的边界条件设置不正确,但是我们不知道如何修改他们。
非常感谢帮助
Andrea_85:
你好
用
Wall
{
type zeroGradient;
//value uniform 0;
}
替换
wall
{
type zeroGradient;
//value uniform 0;
}
祝好
andrea
vonboett:
Quote:
Originally Posted by baedmaen Hi all, now i am confused about multiphase pressure. Could anybody put the corresponding qualitative pressure profiles into the file attached? Thank you very much! |
好的,p_rgh依赖于你的坐标系统,边界条件和最后的pRefValue和pRefPoint。然而,如果你有一些像大气pressure=0Pa在其下为空气通过指定大气边界或者pRefValue = 0 Pa,你的水密度和g被正确指定,你会得到一个的很好的从0到表面压力p三角液压。我在纳维斯托克斯方程方程中使用p_rgh但是我使用p作为一个比较结果。我的附件是两张设置大气环境为0Pa,3D计算域从底部z=2m到顶部z=3m,一个气泡在z=2.25m处,同时一个半径为0.01m的水滴在z=2.51m处。自由面在2.5m处。这个截图是在开始模拟T=0.01s抓取的。在各个方向上的单元尺寸为5mm。注意p_rgh的密度影响,20468Pa=2.51m*g*rho_water。
Dream:
Quote:
Originally Posted by Andrea_85 Hi, replace Wall with wall Best, andrea |
非常感谢!
我们解决了这个问题,原因在于我们大写了单词壁面的首字母W
John Handel Kennedy:
大家好,
我打算模拟直管中的换热
我使用buoyantBoussinesqSimpleFoam求解器
我设置g和beta为0
我的壁面温度是373K并且入口温度为293K
入口的速度为1m/s
管的直径为1m并且为了让雷诺数达到100设置nu的值为0.01
层流普朗特数为1.5
我用simpleFoam求解器得到了充分发展的流动。例如速度增加到2m/s
然而我用buoyantBoussinesqSimpleFoam求解器不能得到相同的速度分布,速度朝出口降低。
我应该如何解决这个问题?
我应该在alpha_t和p_rgh文件中指定什么吗?
祝好
John
转载于:https://www.cnblogs.com/liusuanyatong/p/11259732.html
OF1.7中的p_rgh【翻译】相关推荐
- 【视频联动】编译原理:写出布尔表达式A or (B and not(C or D)) 的四元式序列。说明:按照控制语句中的布尔表达式翻译
编译原理:写出布尔表达式A or (B and not(C or D)) 的四元式序列.说明:按照控制语句中的布尔表达式翻译 这里是总结的知识点.如果有问题可以下方留言提问,视频已经放到Bilibi ...
- 在word中如何设置翻译功能?原来word还能翻译多种语言!
在工作中我们经常要接触到word文件,当我们需要翻译word文件时,却不知道如何下手了!其实在word中是有翻译功能的,并且可以翻译多种语言.有需要的小伙伴可以跟着小编学习一下的啊! No.1在哪里设 ...
- PMP考试中常见的翻译问题
PMP考试的中文翻译,是近年来广受吐槽的一个问题.除了语句生涩.不通顺之外,还有一些词汇和PMBOK或者我们学习时常用的说法不一致甚至是错误翻译的问题.比如kick-off-meeting在 试题中永 ...
- 如何在listary中调用谷歌翻译
文章目录 如何在listary中调用谷歌翻译 1具体设置 2测试一下 3跳转到浏览器 如何在listary中调用谷歌翻译 listary是一个非常好用的电脑搜索软件,我们也可以用它来web定位搜索. ...
- 论文中英对照翻译--(Fusing Multiple Deep Features for Face Anti-spoofing)
[开始时间]2018.10.22 [完成时间]2018.10.22 [论文翻译]论文中英对照翻译--(Fusing Multiple Deep Features for Face Anti-spoof ...
- 【论文翻译】GoogleNet网络论文中英对照翻译--(Going deeper with convolutions)
[开始时间]2018.09.25 [完成时间]2018.09.26 [论文翻译]GoogleNet网络论文中英对照翻译--(Going deeper with convolutions) [中文译名] ...
- 手机上能翻译C语言的app,如何在手机上实现中英在线翻译?
现如今,手机已经成为我们身边不可缺少的一部分了,当然用手机也可以实现很多以前都不敢想象的东西.比如说学习,我们在遇到不认识的英文时,就可以借助手机来进行在线翻译.那么今天小编就来说一下在手机上如何实现 ...
- 如何在Excel/WPS表格中使用百度翻译公式?
之前小编讲了如何在表格中使用谷歌翻译公式,有的小伙伴便问有没有百度翻译公式呢?今天,它来了! 使用百度翻译公式需要五个参数,即要翻译的文本.翻译语言.目标语言.APP ID.百度翻译平台秘钥.前三个参 ...
- Mozilla 正在火狐中集成谷歌翻译;内存产能增长减缓,基本告别降价
(点击上方蓝字,快速关注我们) 转自:IThome.开源中国.solidot.cnBeta.腾讯科技等 0.翻译功能不好用,Mozilla 正在 Firefox 中集成谷歌翻译 你是否知道 Firef ...
最新文章
- 深入研究 System.out.println()
- 二叉树的中序遍历(递归和非递归版本)
- nn.Conv2d中padding详解【pytorch学习】
- 【中级软考】UML图(Unified Modeling Language统一建模语言、标准建模语言)
- 公共课计算机基础怎么样,公共课第一学期《计算机基础》
- Android 系统(21)---OPPO 刘海屏适配
- mysql msdtc 不支持_MSDTC”该伙伴事务管理器已经禁止了它对远程网络事务的支持”的错误(转) | 学步园...
- Python爬虫--爬取电影天堂网站电影下载地址
- matlab usb采集,MATLAB调用研华USB数据采集卡实现数据实时采集_董超.pdf
- Android一步步实现无痕埋点(3)-------虎躯一震
- SpringOAuth2-启动网关Factory method ‘jwtTokenEnhancer‘ threw exception;
- 美团后端2020.4.23笔试题目
- Latex论文中用到的花体字
- python创建目录
- Qt之流式布局(FlowLayout)
- 北邮计算机学硕读博,博士生考试_考博经验——说说我北邮北航考博经历_沪江英语...
- 你后悔进入IT行业吗?IT行业到底如何?
- 中小初创企业网站,该怎么做SEO优化
- MAC OS X 安装Scrapy爬虫工具
- Java实现邮箱验证