CFD中文网

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

@李东岳 好的，谢谢李老师耐心解答

http://www.cfd-china.com/topic/3005 看看这个 我不太熟悉超临界流体

缩小模型只能保证一个无量纲数相同，一般的流动考虑雷诺数就行，缩放会带来一些非真实的物理效应，但是计算量不够就不要讲究了

还有监测的一个点的密度变化
微信图片_20191214213757.png

微信图片_20190306174704.png

微信图片_20190306174712.png

....

@mark Hi mark,

错误提示，缺少twoPhaseMixtureThermo.H文件。可以看看这个文件有没有在你的Make/options中包含

@noodles 开个玩笑，我觉得你这个问题过于简洁，不知道是周期性边界条件还是什么意思
udf可以实现

18年7月底还是6月底投的一篇文章，审稿审了半年，
19年1月中旬给的意见，
19年4月9号修回去了，
现在还在Editor手里

关键是，这是E-QBMM的第一篇文章，part (1)
part (2)在2019年年前投的另外一个杂志

第二篇文章比第一篇文章投稿晚半年多
看起来第二篇文章见刊了
第一篇还没有

截图.JPG

不能的，目前只能二选一。

Hi,

我曾经调用过Fortran库来计算矩阵特征值，不过后来我采用Alberto的库而不用Fortran了。在OpenFOAM里面调用Fortran混编需要安装gfortran，可以在终端下通过：

sudo apt-get install gfortran

来安装。安装后下载你所需要的Fortran源代码（目前我还不会写fortran，我用的是现成的包），里面通过makeFile来制定规则，在终端通过make编译就可以。编译出来的是.so文件。OpenFOAM可以直接调用。

也可以在./Allwmake里面把make写为脚本。

Finite volume methods for hyperbolic problems

https://www.jianguoyun.com/p/DaAmei4Q9s3ZBhjBxKwB

Thermo-fluid dynamics of two-phase flow

https://www.jianguoyun.com/p/DaEadzUQ9s3ZBhjFxKwB

牛顿内摩擦定律与stokes定律间有什么关系？二者计算壁面应力有何区别？

Hi,

你的twoPhaseProperties没有定义，如果你是在interFoam中调用mu，在interFoam的createFields.H中已经定义了：

Info<< "Reading transportProperties\n" << endl; incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);

所以可以这样写：

volScalarField& muc(twoPhaseProperties.mu());

试一下。我没有测试，有问题再联系。

你用的什么类型的动网格？我之前用那个layering感觉挺好用的，一个本来很高的流场被我压到了高度是原来的五分之一，没出现什么大问题。我印象中那个remeshing好像挺难用，不知道你是不是用的这个。

另外我的网格变形都是提前指定的，如果是做流固耦合之类的，网格变形是根据流场计算结果更新的，这个就复杂了，没有研究过。

以后可以有积分制度了，回帖、在线等都有积分，积分可交易等，也和人民币挂钩，分越多，能换的人民币越多

:high: :high:

主要为了刺激大家积极性，也为了给这些一直帮助大家解决问题的人实质性奖励

有的时候icem就是会报错，只能没做一些操作之后多保存，然后再进行。。

@东岳 谢谢东岳老师，我明白如何处理了哈。东岳老师整理的文档和资料都很通俗易懂。

$\epsilon \rightarrow 0$ 并不是$\epsilon = 0$，否则不会有...if $u$ is smooth at $x_1$ and $x_2$...。在粘度趋向于很小的时候，不连续变成具备一定厚度的光滑解，同样承认有厚度的激波。所以
$$
\epsilon\int_\Omega\frac{\p}{\p u}\left(\frac{\p \eta}{\p u}\right)\left(\frac{\p u}{\p x}\right)^2\rd x\rd t \geq 0
$$
另外，
\begin{equation}
\int_{x_2}^{x_1}\left(\epsilon(\eta_q q_x)_x\rd x -\epsilon\eta(\eta_q)_q q_x^2\right)\rd x=\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)-\epsilon\eta(\eta_q)_q q_x^2\Delta x
\end{equation}
考虑一个非常小的$\epsilon=1e-10$，在控制体内$\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)\rightarrow 0$，$\epsilon\eta(\eta_q)_q q_x^2\Delta x$还是大于0.

我现在一打开并行运算模式就报错，请问怎么搞？百度经验那个没有资源包，看了教程也没用，求指导！