CFD中文网

@东岳

just for fun:

They used to say if man could fly, he'd have wings. But he did fly. He discovered he had to. Do you wish that the first Apollo mission hadn't reached the moon, or that we hadn't gone on to Mars and then to the nearest star? That's like saying you wish that you still operated with scalpels and sewed your patients up with catgut like your great-great-great-great-grandfather used to. I'm in command. I could order this. But I'm not because, Doctor McCoy is right in pointing out the enormous danger potential in any contact with life and intelligence as fantastically advanced as this. But I must point out that the possibilities, the potential for knowledge and advancement is equally great. Risk. Risk is our business. That's what the CFD is all about. That's why we're aboard her. You may dissent without prejudice. Do I hear a negative vote? (silence around the table) Engineer, solving that imcompressible N.S. equations using PISO algorithm.

求助各位。我的OF算例计算时间步和输出时间步设置的都比较小（0.002s），导致计算结果有4000个时程文件，请问导入paraview后处理的时候，可以间隔几个文件导入一个么，4000帧的话处理起来太费时费空间了。

替代文字

@bestucan 实际当然是有体积的，概化成无体积就是想简单点，网格逐列关闭也可以的。

实在不行就动网格

我先试试楼上李东岳老师的方法

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

:high: :high:

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

不知道当时的问题是否解决了，我去年做的一部分工作里面也需要相关的处理。当时采用的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

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

整个域中的摩擦速度用的是当地的还是壁面的呢？

壁面的，我刚刚验证过，确定

替代文字

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

@余正东 我都没听过 。。。

结构化网格非结构化网格按照网格是六面体还是体面体区分好像是不完全正确的。我看一个老师讲解说，结构化和非结构化区分是文件保存格式的区分，结构化是使用索引，非结构化是把每个网格都区分出来了。我看各位大佬上面回答可能是直接讨论的四面体网格，因为我最近也在学openfoam和icem，icem可以通过自动化分成四面体非结构网格，也可以分块化成6面体网格，但是我通常画完6面体网格后，icem有可以转化成非结构网格按钮，将画好的网格转化成非结构体网格（转化后依然是6面体）。
所以如果你要话非结构化网格也可以分块（画了6面体转化一下），如果你说的是画四面体非结构化网格分块的话，这方面知识我就不知道了

多谢老师解答，问题已经解决了

reFineMesh可以加密一个方向，这个需要指定区域，你试试topoSet+refineMesh呢》

@肖艺 只是看过相关的影片.. :zoule:

topo.png
我需要用ICEM画multiblock网格，几何类似于上图的模型，我知道三角形Y-block和圆形O-block，那请问这两个如何结合起来呢？在这种结构中它的edge并不共用相同的vertices
第二个问题是icem的multiblock的网格如何光顺？我尝试使用了pre-mesh smooth，但是这个功能很慢，容易卡死。我只知道可以用edit edge调整spline来改变edge形状，如果一个一个edge调是不是有点笨？

好像不能发动图，我试试截几个图吧（前两张是计算“正常”时，后面两张是出现了震荡和最后崩溃时的情况）
diver1.jpg diver2.jpg diver3.jpg diver4.jpg

参数波动跟监测点位置有关，放在分离区跟放在直管中间是不一样的，另外，感觉稳态计算是监测面上的通量吧，通量守恒就可以认为收敛了吧。

微信图片_20190306174704.png

微信图片_20190306174712.png

....