CFD中文网

各位大佬，想请教下模拟搅拌喷气，气液用欧拉双流体模型，固体颗粒用dpm模型，最终颗粒受力是怎么计算的？是气固的力加上液固的力吗？欧拉欧拉双流体模型流体是相互渗透的，那么在计算颗粒受力是到底是哪一相的力还是说两相的力

块就是将你的整个计算区域分割成几个区域，每个区域有各自的网格分布规律；
模型应该是你计算的几何外形吧

请教各位老师和同学一个问题：我在使用sprayFoam求解器的时候，想把液滴初次雾化模型 atomizationModel设置为LISA模型，请问该模型得到的液滴直径分布是否会和injectionModel中设置的初始粒径相互冲突吗，foam中是怎么处理的？

@东岳 搜了下相关论文，以前叫HyperFlow，不知道是不是同一个软件

时间离散格式应该是计算时间导数的格式，模板没有限制；而推进格式应该要求模板是单侧偏心的。http://www.scholarpedia.org/article/General_linear_methods是一般的时间推进格式，偏心模板最多也就是隐式格式。

请教各位，我使用fluent计算二维LNG储罐发生密度引起的液层翻滚以及气象空间超压过程。气液之间使用VOF模型，液相内部使用组分输运模型，也就是多相多组分模型。其中气液相变之间使用Lee模型udf进行。

QQ20250311-113118.jpg

在计算的过程中我这么设置：假设流体是甲烷和乙烷的混合物，并且液相的密度是随着温度和组分发生变化，随着乙烷组分的增加密度增加，随着流体温度的增加密度减小。当流体温度小于甲烷饱和温度就二者都发生冷凝，当流体温度大于甲烷饱和温度小于乙烷饱和温度就甲烷发生蒸发乙烷发生冷凝，当流体温度大于乙烷饱和温度就甲烷和乙烷都发生蒸发。

在我模拟的过程中，开始未发生多组分液层混合过程，气相空间的压力是逐渐变大的，但是当混合完成后，虽然气相和液相空间的温度都在上升，液相空间的密度下降，但是气相空间开始压力下降了，我不明白这是为什么，所以想请教下各位大佬。

0DDFF043D7976A246C818A596DF9F6CD.png

同时，在我进行Lee模型使用过程中，我发现有人使用Hashemi-Wesson这个表面蒸发公式进行蒸发率的计算，Lee模型的量纲是kg/(m^{3} \cdot s)，而这个公式的量纲是kg/(m^{3} \cdot s)。那这个纯蒸发的模型是怎么应用在fluent中，我也想请教一下大佬们。

QQ20250311-094104.jpg

好吧！我先自己回复下自己。在稀疏矩阵存储中，通用有两种格式，如下矩阵

A = [ a00 a01 0 0 0 a11 a12 0 a20 0 a22 0 a30 0 0 a33];

引用一篇博士论文 CRS scheme stores the non-zero elements the matrix in left-to-right and top-to-bottom order in vector v (row-wise storge). 存储格式如下

v = [a00, a01, a11,a12,a20,a22,a30,a33] cI = [0,1,1,2,0,2,0,3] rI = [0,1,1,2,0,2,0,3]

但是在openfoam中，不是这么存储的
openfoam是一种 LDU 格式 对角线 lowerAddr() upperAddr()

A = [ d0 u0 u1 0 0 0 l0 d1 u2 0 0 0 0 l1 l2 d2 u4 u5 . . . ] l = [l0 l1 l2 ...] d = [d0 d1 d2 d3 d4 ...] u = [u0 u1 u2 u3 ...] L = [0 0 1 1 2 2 3 ...] U = [1 2 2 3 3 4 4 ..]

具体可参见lduAddressing.H 头文件
因为计算流体力学中，大部分矩阵都是对称（因为一个internal face对应两个体）除了特别的边界条件外，基本上都是对称矩阵。而矩阵对角元素就更好求了，大部分情况直接把该行相加就得到主对角元素。这种sparse matrix设计是符合cfd计算原理的，但是我的问题还是没解决。看到相关资料分享一下。
稀疏矩阵也分好多种，openfoam中的renumberMesh就是尽量保证在是多对角阵。
谢谢！

@搬运工
我们可以用涡和波的方法，描述湍流。说湍流是涡，比如，Kolmogorov理论：湍流是由大涡，中涡，小涡，和最小涡组成。涡的能量的转递维持这些涡的生存和运动。每种涡都有时间尺度和空间尺度。时间尺度指的是涡旋一周用的时间；空间尺度是指涡的空间大小。可见，湍流有很多时间尺度和空间尺度。

湍流也可以用波来描绘。波理论来描绘湍流，湍流是由各种波组成。分解湍流到频谱空间，你可以得到各种波福和波频。这些波福是湍流的空间尺度，而波频的倒数是湍流的时间尺度。

@搬运工 后来我发现我这个问题不是周期性问题，虽然结构具有周期性，谢谢你的解答:happy:

有的时候，某些内部的情况，没有几何，但是需要对网格进行一些处理，就需要手动的在不能关联的情况下，切一下。 某些需要投影的情况，比如地面上悬空的一个圆筒。在地面上只有圆筒的投影，也不存在几何关系。需要把edge投影过去，就需要spline方法。

@深蓝 请问您找到输入粒径分组的方式了吗？我最近也要做这方面，如果您看到了请您教我一下

按照文献中的案例模拟了一个三维方腔流
我的垂直速度曲线图和水平速度曲线图
BOT8L352FD5IWM~S5IFVX42.png MS0N{_5W.png
和实验结果不一样
QQ图片20190306104704.png
是什么原因造成的？

udf有问题，你试试这个

#include "udf.h" DEFINE_PROFILE(inlet_velocity, t, i) { real x[ND_ND],y,vel; face_t f; begin_f_loop(f,t) { F_CENTROID(x, f, t); y = x[1]; if(y>0.0393) vel = 29.5; else vel = 29.5 * pow((y / 0.0393), 1/7); F_PROFILE(f, t, i) = vel; } end_f_loop(f,t) }

new cannot satisfy memory request. This does not necessarily mean you have run out of virtual memory. It could be due to a stack violation caused by e.g. bad use of pointers or an out of date shared library

如果你遇到这个问题，一种解决方式是把你的求解器重新编译一下，可能是bad use of out of date shared library。原因是你的库重新写了一些代码并重新编译了，但求解器没跟进

痛苦，结题更痛苦

fluent LES模拟螺旋列板抑制涡激振动时高度0.1D和0.2D的螺旋列板升力系数幅值在0.1和0.15左右，但0.15D的螺旋列板的升力系数幅值却到了0.5。三个都是同样的网格划分与时间步，为什么差异会这么大？

相变，多相气液固 、传质、传热 、
要是能把天然气水合物这玩意 模拟了，，那也差不多千古流芳了 哈哈哈

一个小白 小白 不能再白的人 在一群大佬环绕中 悄摸摸的说到

这更像是用代码格式化软件统一生成的。高级编辑器自带的代码格式化一般只改行前缩进，代码格式化软件是全部重新排版。

一般大型代码库都有非常细致的格式要求，代码格式化软件可以非常细致的定义每一个格式化要求：每行最多字符；算符前后加不加空格、换不换行，前换 后换还是都换……

然后用被处理过的代码非常满足每一个细致要求，所以某些地方就看着有一点点别扭了。 indent是GNU出的格式化软件。里面有一些常用的风格模板：GNU、K&R、linux。也可以跟“像一篇文章一样长的”参数来定义自己的代码排版风格。

在多区域（液体和固体换热）计算中，某一个物理量（温度T）在两个区域间的interface是一个连续的量，即没有热阻，这个interface边界条件该怎么设置？
在OpenFOAM自带的多区域算例中，这种耦合面用的边界条件是turbulentTemperatureCoupledBaffleMixed，这个其实是根据边界面的值是根据边界上的值来计算传热系数的，里面写的很麻烦，有点看不懂。有没有哪位老师指点一下？多谢多谢~

很普适性的问题。如果你说的是如何判定非稳态的收敛性的话，可以监控一些参数，如果稳定波动一定周期，就可以认为收敛了。