东岳老师，请问您的SCI文章《Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model》里模拟鼓泡塔时的入口速度是如何设置的呢？直接速度入口，表观气体速度和论文里的相同吗？文中的这句话“At the inlet, the velocity of the gas phase was calculated from the superficial velocity in different test cases;”如何理解呢？谢谢！
您好，请问我用fluent计算这个算例，为什么整体gas hold up相差一个数量级。入口气体速度设置的0.0213m/s。设置基本和论文里相同。
https://pubs.cstam.org.cn/data/article/mie/preview/pdf/lxysj2011-045.pdf
感觉就是附近涡造成的。我在哈尔滨上学的时候也见过树周围的这种现象，确实都是出现在风比较大的地方。
不知道当时的问题是否解决了，我去年做的一部分工作里面也需要相关的处理。当时采用的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
