@merlinizlf 这就是发散了 但是我遇到这种问题时候 经常跟网格有关 我网格加密 然后尽量使他每个方向长度均匀 就好一些

大佬，我用的matlab，不知道大佬用不用，写了个简短的字符串处理，读取log文件，文件名aaa.log。查找Time=后的数据和Inconsistency value, average, min, max后面的数据，然后输出到data.txt，第一列时间，第二列是那个数据。用大佬提供的这三个时间步的log信息试了下，应该没啥问题。只是每个时间步Inconsistency value这个有两组，我取的第二组。

clear;
clc;

filein=['aaa.log'];
fidin=fopen(filein,'r');
ansn=zeros(500000,1);
anst=zeros(500000,1);
totn=0;
tott=0;

line=0;
flagn=0;
flagt=0;


while ~feof(fidin) % 判断是否为文件末尾  
    line=line+1;
    temp=fgetl(fidin); % 从文件读行
    if flagn ==2 
        flagn=0;
    end
    flagt=0;
    str = deblank(temp);
    s = regexp(str,  '\s+', 'split');
    if length(s) <=2 
        continue;
    end
    if (strcmp(s{2},'Inconsistency')) && (strcmp(s{3},'value,'))
        flagn=flagn+1;
    end 
    if (strcmp(s{1},'Time')) && (strcmp(s{2},'='))
        flagt=1;
    end 
    if (flagn==2)
        ss = char(s{10});
        sss=str2double(ss);
        totn=totn+1;
        ansn(totn) = sss;
    end
    if (flagt>0)
        ss = char(s{3});
        sss=str2double(ss);
        tott=tott+1;
        anst(tott) = sss;
    end 
end

fclose(fidin);
fid=fopen('data.txt','wt');
for i=1:tott
fprintf(fid, '%.6f %.8f\n',anst(i),ansn(i));
end
fclose(fid);

输出的数据：
0.502963 0.00282694
0.503704 0.00283924
0.504444 0.00285154

您好，想问下基于OpenFOAM写新模块 mpi并行那些需要考虑吗？因为课题组的代码面向过程的，写每个部分时候都要考虑mpi的消息传递。OpenFOAM这么高度面向对象，是不是按照串行写就行？

难道勾选这个不是显示圆柱实体吗？

@zp气固两相流 在 IBM浸没边界法生成拉格朗日标记点 中说：

A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries

好的谢谢，我去看看这两篇文章

@zp气固两相流 好的谢谢，其实我对这几种IBM的类型不是了解，只是看了下课题组一个IBM的代码，固体对流体反作用是加在NS的体积力里的，就是你原来想植入的那种吗？ 但是看extend那个版本IBM的描述 Discrete forcing approach with direct imposition of boundary conditions， 感觉也像这种方法哎。 里面ppt讲的Continuous forcing approach 和 Discrete forcing approach的说法都把我搞懵逼了，和很多IBM的文章叫法又不一样。

@zp气固两相流 同学你说那个extend版本的IBM不是你想要的IBM是什么意思呢？

@东岳 感谢老师解答！

查了一些资料，表示还是很疑惑，希望有老师能给我解答一下，谢谢！