CFD中文网

@Wayne 对，电磁场里面会用到

看一些边界条件，都会定义多个构造函数，每个构造函数的输入量不同，为什么需要多个构造函数呢？

我们都知道，速度可以做Helmholtz分解，用速度势的梯度和流函数旋度叠加的形式，求解二维粘性可压缩翼型绕流中的速度势和流函数相当于求两个泊松方程，但是需要给出边界条件，一般情况下，在速度场已知的情况下，即涡量场和胀量场已知，该如何分别给速度势和流函数的物面和远场边界条件？

半新手，看其他论坛上可以实现这样的操作，即把一个复杂的几何体分成多部分分别划分网格之后再进行拟合，但是我自己试了下建块merge后导入到fluent中出错，请问该如何操作呢？不胜感激

这个问题我解决了，就是在周期性关联的时候，在edit periodic/periodic veteric 时要在中心轴上把所有的veterivic都要关联。即一个点要点两次

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

@cfd-china 我觉得不是，并不是大于$\Delta$的窦唯大涡。

ffmpeg -i input.MP4 -c:v libx264 -preset slow -crf 40 -c:a copy output.mp4

Jacobi是并行性最好的算法，LU-SGS并行性较差。

首先在稳态条件下模拟一个风流场，然后利用dpm模型模拟压力旋流喷雾过程，试跑的过程有的方向能喷， 有的方向喷不出来是为什么，新手求帮助ebb4c6e2-270a-497e-a003-fb08d398d4bd-image.png
z两个方向可以喷，但是换到-y方向时候就喷不出来

@东岳 李老师，我试了一下，把出口和入口速度设置的速度一样，在开始的时候水池的水就不稳定，在入口处水面上升，在出口处水面下降。得好久才能接近一种平衡，而且导致液面晃荡。
1.png

2.png

大家好，请问各位老哥有使用过salome前处理的吗，最近打算研究一下，可以交流交流:chouchou:

@搬运工

高斯迭代法处理20阶矩阵就已经非常慢了。通常和多重网格+高斯赛德尔求解器一起处理。并且OpenFOAM这种成熟的商软里面矩阵求解器算法非常精炼。通常由主求解器+预处理+光顺+重组构成。多重网格更是复杂。另外，在编程的时候使用动态、寻址的方法寄存系数。

所以如果你要写一个快速的求解器：

使用外挂包； 潜心研究稀疏线性系统求解理论；我建议花大量时间研读：YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:

替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程，现在模型验证一直没弄好。首先是编写了一个简单的UDF，采用热平衡模型，导入fluent后提醒我定义比热的那一块有错误，求各位前辈帮我看一下；其次，我添加了Boussinesq项后，优化了一点点液相云图中熔融线一直平行的问题，但是与论文中的还是有差距，且我的熔融线上下厚度保持一致（实际上应该下面厚很多），UDF中还需要添加什么项来解决这个问题呢？以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图，求各位帮我解惑，不甚感激。2.png 3.pn替代文字g 4.png

谢谢，那个我也看了（但没仔细看，只是特意找了扩散项）。应该也是和Menon一样的结论。因为Menon的文章截图比较方便，才引用他的

2422ba2b-f9b8-402e-936e-8075325aae10-image.png 例如我生成一个part名字叫多“cavity1”，但是导入到fluent里面会多出“cavity1：007，cavity1：008，cavity1：008shadow，”这几个边界fb6ad2d6-0c7d-4266-a7b4-fa752482f0a6-image.png 。相应的网格也被分裂出好多细块。微信图片_20190126110946.jpg 微信图片_20190126110950.jpg

a80fed28-45ef-4507-a08d-13abf19c7ccc-image.png 这个udf具体该怎么编写和使用？

@李东岳 十分感谢！

终于注册成功了，感谢李老师！