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,

p_rgh seems to be the dynamic pressure. I work with pimpleFoam and
pisoFoam. Is it possible to adjust these codes to obtain
p_rgh?



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?

First figure:

stationary two phase system (water and air); i think nbadano
already posted the answer at December 3, 2010, via this
thread

Second figure:

bubble rising to surface; snap shot

Third picture; first contact of a water drop with water
surface



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

{

type zeroGradient;

//value uniform 0;

}





with



wall

{

type zeroGradient;

//value uniform 0;

}



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