CFD中文网

@G8S7

究竟能不能用欧拉-欧拉模型来说明这个体系下固体颗粒的运动情况，

如果有实验图片的话更好判断。取决于固体颗粒粒径，相分数，粘度等。如果用你所说的VOF+欧拉-欧拉才能描述这一个过程，也是可以的。另一种方式是三相欧拉，在皮层熔体中对界面进行尖锐挤压处理。

另外，用什么模型取决你你关注的重点，如果关注内部的运动，皮层溶体不是重点用欧拉无妨，VOF也是有耗散。

At last, 几何很简单，是验证算法的利器，不用OpenFOAM可惜了。

是的，我觉得，本质上，还是应该走看书 - 手算 - 编程这个过程，如果看书直接去编程，可能理不清思路。

@东岳 多谢李老师了，是时间步长的原因，我开始按库伦数取得 听了您的建议往小了试，没有出现发散

:141: :142: :143: :144: :145: :baobao:

想设置一个壁面接触角随时间增加减小，用公式描述
QMF9EZM9NI7B{$AS17ASF{M.jpg
但是模拟过程中不到0.01s就已经完全铺开了
8d0dfcc4-1ada-4bec-8a31-d4000f59d818-0~FKZ}KBX0WEF3X7@2%@NTD.png
接触角为定值的时候没有问题，请问接触角随时间增加减小的设置哪里出了问题，该如何设置

@东岳 让孩子赢在起跑线上:high:

loop的是面网格上的所有网格。你的需求应该能实现，最简单的办法是，每个面网格中心去读一遍文档，把散点距离自己最近的那个点的速度赋给所在网格。但是会有其他问题，散点疏密是不是和面网格疏密比较匹配，也许一个面网格配多个散点或者反之。

不知道当时的问题是否解决了，我去年做的一部分工作里面也需要相关的处理。当时采用的Matlab对导出的数据进行处理，不过在数据处理之前需要将 x,y当中一些极小的值变成0，比如10的负十几次方的值，因为这些fluent输出的值可能会对处理产生影响，然后导入程序就可以变成矩阵形式了，程序大概如下。供以后有需要的朋友们参考。

data1 = xlsread('D:\User.....');

% 获取X、Y坐标和密度值的向量
X = data1(:, 1);
Y = data1(:, 2);
density = data1(:, 3);

% 确定X和Y的唯一值，并按升序排列
unique_x = unique(X);
unique_y = unique(Y);

% 创建行列编号映射表
x_mapping = containers.Map(unique_x, 1:length(unique_x));
y_mapping = containers.Map(unique_y, 1:length(unique_y));

% 初始化矩阵
num_rows = length(unique_y);
num_cols = length(unique_x);
densityMatrix1 = zeros(num_rows, num_cols);

% 将密度值填充到矩阵中
for i = 1:length(X)
row = y_mapping(Y(i));
col = x_mapping(X(i));
densityMatrix1(row, col) = density(i);
end

那位OpenFOAM经验丰富的前辈能讲讲 OpenFOAM的算例之间如何耦合？

\begin{equation}\label{MM3}
\rho\frac{\partial \mathbf{U}}{\partial t}+\rho\nabla \cdot (\mathbf{U} \otimes\mathbf{U})=-\nabla p+\nabla \cdot(\mu \nabla \mathbf{U})+\left(\rho_\rref\left(1-\beta\left(T-T_\rref\right)\right)\right)\bfg
\end{equation}

@bestucan 在别人电脑上是正常的 ，很奇怪，我最后重装电脑解决了

请参考无痛苦ns方程笔记的这一节：R、sigma、prime2Mean、以及雷诺应力相关量

@wwzhao 再次感谢:xinxin:

@wwzhao 谢谢两位老师，文献DOI: 10.1115/1.1471361中有提到离散求解的

ICEM画二维网格导出时报错：
WARNING: Mesh has uncovered edges. ANSYS Fluent needs a complete boundary (lines in 2D) or it will give a variety of errors and not read in the mesh! If this was 2D Hexa, perhaps your edges are not associated with perimeter curves

模型描述：二维管道流体的测试网格，使用ICEM导入模型后：
1、划分块
2、在preMesh 里 使用Edge Paras进行边界层的绘制，
3、preMesh
4、在preMesh里Convert to Unstruct Mesh
导出网格选二维，其他未改变，出现上述报错
新手从B站学了一点划分网格操作就尝试了这一个，不知道哪里出错了，还希望前辈们解答

好像是因为你发散了想看连续性误差？这个没什么用。不能表现发散的本质信息

如果你想看的话，openfoam里面这个叫continuity error，有个sum，有个global，有个cumulative，可以参考sum和global

完结帖子，后来找到了，我的openmpi不是装在系统文件夹的，换一下bash文件里的export WM_MPLIB=OPENMPI就好了（如果openmpi-x.xx.x装在系统文件夹就是SYSTEMOPENMPI)

x3d
用这个导出，导入到blender

我用子系统上的OpenFOAM，总感觉并行计算没用，和我不开并行计算时间差不多。。。而且还decomposePar不了，好像是安装的时候少了啥东西。。。

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