xjwang
帖子
-
@bestucan 感谢您的回复,抱歉之前不知道为什么没看到。好几年过去了,今天刚看到这个之前问的问题,我那个卡住的原因是:(拿好笔记本记好了,这个问题值100万)我刚开始不会用openfoam,然后硬盘满了,所以算不下去了。
-
-
@dxl 注意你用的是reconstructParMesh 还是reconstructPar.
比如,你可能用60核来用snappyHexMesh画网格:decomposePar mpirun -np 60 snappyHexMesh -overwrite -parallel
下一步可以用
reconstructParMesh -constant
来合并成最终的网格,下一步需要做模拟的时候,你可能用300核做并行计算
decomposePar mpirun -np 300 pisoFoam -parallel
最终你想把某几个时间步的结果画出来,所以你需要合并一下
reconstructPar -time xxx
-
@coolhhh 也许DSRFG方法是说用的能谱,但最终还是用的von karman的频谱,就像你说的:
@coolhhh 在 OpenFOAM怎么保存指定不等间距时刻的结果 中说:
用频谱叠加来表示三维能谱
关于CDRFG方法,他们的code里是对每一个点都生成一组PQK,所以我才说有点鸡肋。
@coolhhh 在 OpenFOAM怎么保存指定不等间距时刻的结果 中说:CDRFG考虑的是频谱,确实没必要类似kraichnan生成球面分布的。非均匀问题感觉是很难解决,不仅RFG系列,其他合成方法也都有类似问题。就如RFG文章所说,如果非均匀变化足够缓慢,近似满足零散度。但从模拟过的工况看,影响更大的是入口质量通量平衡,这决定了流场中是否会产生非物理压力脉动。
关于生成inhomogeneous的inflow,后来我想着试一下vector potential这种思路,貌似也行不通。
-
@coolhhh 我几天内发给你啊,这两天有点事,还需要整理一下!
-
@coolhhh 论文发表出了点问题,审了两年了,貌似杂志那边有点变动,然后还有些其他因素,你如果需要的话,我可以私发你一个pdf看看,但我只是用latex写道了我的dissertation里边了,也许需要一两天的时间整理成一个单独的pdf。
-
我这个模拟就是用的smagorinsky SGS model, 和前边某大佬说的一样,压力松弛是1.
-
TPU的实验数据是一个小尺度的模型(缩尺比比较大),所以应该是能得到不错的模拟结果的,你可以参考一下Patruno之前发表的文章,有两篇是关于TPU数据的。我做的是NIST数据库的,长度比1:100,速度比1:4,也是个小尺度的模拟,结果还算可以,peak suction pressure 有的地方不太一样。
-
@coolhhh 是,对于谱方法来说,类似于DSRFG方法,这个确实是比较难满足。还有个问题是inhomogeneous的问题,虽然那个CDRFG方法声称能产生inhomogeneous wind filed,但是感觉那个方法很鸡肋,完全没考虑divergence-free,所以他们根本就不需要用kraichnan那个方法。后来看到过Xuesong Bai一篇文章,考虑了inhomogeneous的问题,但是是一个沿三个方向同时变化的,对于做风工程,只有在竖向是变化的。再就是空间相关性的问题,后来我还试着对DSRFG做一些改进,但发现都是徒劳的。Kraichnan方法演变出来的一系列方法,最终看来因为方法过于简单,很难做出一些新的变化。我感觉DSRFG最适合用来做桥梁,因为桥梁高度不高,那一段风场可以看做homogeneous,类似于BARC的那些研究也很适合。
-
@coolhhh 第二条你说的满足NS方程,感觉现在大多数都是满足divergence-free, 不知道有没有同时满足momentum equation的?一个简化的满足momentum equation的方式是满足taylor frozen hypothesis. 之前看patruno有一篇文章说的是这个,大佬有没有其他的推荐?
-
@turbulen 没有见过你说的AR线性滤波, RFG 算是一种谐波合成法吧, RFG DSRFG 以及G.T. Bitsuamlak "A divergence-free inflow turbulence generator using spectral representation method for large-eddy simulation of ABL flows"这个方法,都应该算是一类的。做LES还有用Synthetic Eddy Method 和Digital Filter Method的,这两种方法在openfoam.com中有,如果你在用openfoam.org,可以去找一下这个网址:https://github.com/NHERI-SimCenter/TurbulenceInflowTool/tree/master/openFOAM_code/OpenFOAM6/turbulentInflow
-
个人感觉,如果犯了一个不太常见的错误,在cfdonlie这种论坛也搜不到太多类似的帖子,很可能是犯了一些很基本的错误。就好像我一直去搜autoMap 和rmap这两个member function怎么用一样,发现讨论的人特别少,那很可能别人在这两个function都没犯过错误,也意味着,这两个function可能不会出错,虽然最终搞明白了这两个function的用法,但对解决问题完全没有任何帮助。
-
问题解决了,原来是有个constructor写错了,本来这个错误可以避免的,用foamNewBC生成的模版就可以了,但是我的code是根据一个of v2112版本的code改的,我用的是openfoam v6,所以没注意,把一个constructor改了一下,导致出现了这个错误。
-
最近写了一个inlet boundary condition的code, 并行计算的时候一直报错。具体问题是:
我在0文件下定义了:
inlet
{
type DSRFGVelocity;
M 600;
N 100;
uref 22.5091;
alpha 0.1538;
zref 1.3860;
inflowPar pkqomega;
value uniform (0 0 0);
}
如果单核计算没有问题,但是等到并行计算renumberMesh的时候 会报错。我是这么做的,当mesh都处理好以后:
decomposepar
mpirun -np xx renumberMesh -overwrite -parallel
mpirun -np xx pisoFoam -parallel.
报错信息是:
--> FOAM FATAL IO ERROR:[0] compound has already been transferred from token
on line 33 the empty compound of type List<vector>
[0]
[0] file: /ddnA/work/xuwang/DSRFGinflow/processor0/0/U.boundaryField.inlet.value at line 33.
[0]
[0] From function Foam::token::compound& Foam::token::transferCompoundToken(const Foam::Istream&)
[0] in file lnInclude/token.C at line 99.
[0]
FOAM parallel run aborting
然后打开processor0/0/U 文件后,那个inlet是这样的:
inlet
{ type DSRFGVelocity; M 600; N 100; zref 1.386; alpha 0.1538; uref 22.5091; inflowPar pkqomega; value nonuniform List<vector>
1372
(
(27.414161 3.7652743 -1.313935)
我能感觉到的就是,在0文件下U的value被定义成了uniform (0 0 0),但是程序一旦运算,就把inlet 上的速度值给重新计算了一次,导致变成了nonuniform。但是of6中自带的turbulentInlet算例,也是这样做的,不知道为什么我这里出了问题。我自己写的inflow的updateCoeffs()函数如下:
void Foam::DSRFGVelocityFvPatchVectorField::updateCoeffs(){
if (this->updated()) { return; } if (curTimeIndex_ != this->db().time().timeIndex())
{
scalar time = this->db().time().value(); vectorField Uavg(patch().size(),vector::zero); vectorField Uprim(patch().size(),vector::zero); vectorField Utotal(patch().size(),vector::zero);
xxxxxxx 省略中间一些繁琐的步骤
this->operator==(Utotal);curTimeIndex_ = this->db().time().timeIndex();
}
fixedValueFvPatchVectorField::updateCoeffs();
}
请各位高人指点!
-
@coolhhh 感谢回复!
-
@coolhhh 请问你这个vandriestCoeffs中的delta是怎么设置的?用的cubeRootVol吗?你在最开始贴的那个代码,应该是复制的Openfoam basic trainning那个文档中的设置吧:
@coolhhh 在 Smagorinsky模型及边界条件设置问题 中说:
vanDriestCoeffs
{
delta cubeRootVol;官方文档中这么是这样的:
@coolhhh 在 Smagorinsky模型及边界条件设置问题 中说:
vanDriestCoeffs
{
delta <geometricDelta>;搜了一下也没找到具体的例子,希望不吝赐教!
如果用van driest的话,是不是nut在你的building表面是不是可以写成fixedvalue 0?还是只能在网格比较精细的时候才能使用fixedvalue 0? -
@vortex
感谢分享! -
@sunss
感谢分享,前年看到他们那个工具包的时候还很简单,最近看了看已经变得很完整了,还有一些算例。你也是做风场模拟的吗? -
@霜染丹枫
我做的风场模拟,所以我检察风洞模拟是否正确是通过检察计算域内的风速时程的功率谱是否和目标谱一致,也就是是否和在入口处生成的风速满足同样的能量分布。
不知道你们槽流最关心的是哪一个数值特征。比如,你可以验证turbulence intensity也可以验证power spectra,但不能通过一个瞬态的流场来判断模拟的好坏。一直不清楚你通过离散能量谱得到的脉动速度为什么要通过观察瞬时速度场来判断模拟的好坏,就算有很强的紊流特性也不能保证这个速度场的能量分布就和你目标谱一样。
上文中@sunss 同学说的那个tinf,你可以看一下,其中有dfsem, dfm等方法,特别是对于槽流,openfoam的tutorial中有一个RE395的算例(也许把名字记错了),应用的是DFSEM方法。
如@sunss 同学所说,也许高斯谱更适合你的功率谱,这样的话你可以直接应用fluent自带的那个对应RFG(Smirnov et al,2001)的生成方法,具体名字不记得了,你可以自己去查一下。
在其他回复中我提到过,DSRFG方法虽然理论上来说是无散度的,但如果应用在数值模拟中,就不是严格的无散度了,需要考虑网格尺寸,所以我一开始就问你网格尺寸和紊流积分尺度的数值。我们完成了一篇这样的论文,希望早点可以和你们分享。 -
@霜染丹枫
不知道你想模拟什么,想和你确认一下是不是模拟ABL,没有得到回复,单凭一张瞬态的速度场图,没法判断你的模拟结果好与坏。 -
@sunss
文件格式应该是:
points:
.......................................................
(
(1 1 0)
(0 1 0)
(1 0 0)
(0 0 0)
)
...........................................................
U:
..........................................................
(
(三维速度向量)
(三维速度向量)
(三维速度向量)
(三维速度向量)
)
.....................................................
所以你的设置有两个问题,一个是标出了4或者8这种数量,这个需要删掉,一个是需要把速度分别存储在对应的文件夹下边,不能在一个文件中。还有那些openfoam的星号可以删掉哈哈 -
@sunss
应该是在constant/boundaryData/inlet文件夹下边有一个points的文件,然后每一个时间步都是一个文件夹,其中包含一个U文件。比如说,我在inlet文件夹下边在terminal中输入ls,会显示:0 0.001 0.002 ...... points
而在每一个时间步文件夹下输入ls,比如在0.001文件夹下,会显示:
U
-
@李东岳
哦哦,因为知道您不用fluent 看到您贴出这个图 我还以为paraview的后处理也能这么画! -
@李东岳
李老师 突然想起来,您这个图是fluent的图还是openfoam的图?如果是openfoam的图能不能提示一下怎么画出来的?用的什么filter? -
@sunss
这个设置不是很难 你可以看一下这个视频 以及Eddylicious这个软件的document,他们里边也是用的这个边界条件
https://www.bilibili.com/video/BV1CV411C7Jr -
- 个人感觉是能画等号的,当f=kU. CDRFG方法中,对于定义p和q的时候,直接代入了Su,Sv,Sw, 而且还说有个系数halved,但后来检察代码发现没有减半。还有一个就是CDRFG方法是用来生成inhomogeneous inflow的,虽然从原始的RFG(Smirnov)中说,这个方式在生成inhomogeneous inflow的时候也可以近似看做零散度,但从CDRFG方法那个计算k,p,q的过程来看,应该能造成divergence.个人感觉这个方法之所以很多人用,也许是因为他们把程序贴在了文章后边。不过这也只是我个人感觉,如果我有错误的地方,还请指出来!!!
- Lu,Lv,Lw虽然理论上来说不影响turbulence intensity,只是按你所说的那样,造成峰值的不同,但具体用到这个方法中,会有所影响。我们最近正完成了一篇论文讨论了这个情况,希望能快点见刊与大家分享。
- 这个问题比较常见,虽然理论上来说,LES是对“大”涡进行解析模拟,而“小涡”进行SGS模拟,但还是会对高频处那些本应该解析模拟的部分不能精确模拟。
4.严格说来,你看到的DSRFG方法并不是严格零散度,因为要考虑网格离散的影响。
-
@李东岳 对,这个可以理解,但是衰减太快就是有问题了,一般处理这种瞬态问题用的piso算法,会根据质量守恒去修正求解的速度,所以如果入口处不满足质量守恒,速度场会有很大的改变。
-
@霜染丹枫 还有你说的那三种Ls的计算方式很重要,但这个是影响空间相关性的。这个参数不是自己凭空造出的,我是根据试验来确定的
-
@李东岳 internalField应该没啥影响,毕竟模拟时间比较长,初始条件可以忽略不计。我试过internalField是(0 0 0)和(U 0 0) 都没啥影响
-
@霜染丹枫 Ls太小吧,你这个算例感觉设置的不是太正确。
题外话:
1.你的首层网格高度0.0001m,不知道你具体咋划分的网格,沿高度方向如果增长太快的话,应该会对计算结果有影像,而且你网格的三个方向的尺寸比例如果出现这种1:1:1000这种情况的话,也会影响你的计算结果。
2.不知道你是不是想摸你ABL, 如果是的话,上部设置成无滑移边界条件不太好吧。 -
@李东岳 李老师,请问那个公众号下边回复说用CDRFG方法的同学或老师注册了吗?能否邀请他一起讨论一下那个方法?我一直感觉那篇论文中有一个错误,给作者发邮件也没回复。
-
@霜染丹枫 “根据自己的工况” 这句话,如果你的参数是从风洞试验中或者其他数据库中得到的,那应该没啥问题。不过我觉得,你貌似没有把紊流积分尺度等参数设置正确,一个比较快速的检查:你的网格尺寸是多少,你的紊流积分尺度是多少?看你发的那个图片,貌似是遇到了我之前那样的问题,脉动部分急剧减少。我现在手边没有我计算的图片,不过我的计算结果的瞬时速度场在整个计算域都能表现出明显的湍流特征。
-
@李东岳 这个不仅仅是流槽了,算是风洞。虽然边界条件会有影响,但是不太影响计算域中的一些地方。
-
@霜染丹枫 换了台电脑回复你^-^
能不能贴一下你不同位置的速度的紊流度以及功率谱的图片?
还有一个问题是,比如你的计算域顺风向长度是15,你目标结构应该放在5-8左右,这样紊流度等衰减都不是很大。
我之前犯的错是,入口的参数搞错了,紊流积分尺度这些设置错了。
以下是在公众号下边的留言:
抱歉在论坛上提问问题后没有及时更新后续。
很惭愧,这个问题后来我发现是参数设置错了,生成风场的参数是按照缩尺的参数选取的,而DSRFG这个方法,本来是对网格尺寸不敏感的,但因为参数选错了,所以导致了在入口处引起了divergence-free问题。至于后续问的那个沿顺风向紊流度减小的问题,确实是会减小,不过不是太过于明显,因为我错误选取了参数,导致顺风向的衰减程度特别大。之所以会问在两个面定义成inlet,是因为看过一篇文章,Kim, Y., Castro, I. P., & Xie, Z. T. (2013). Divergence-free turbulence inflow conditions for large-eddy simulations with incompressible flow solvers. Computers & Fluids, 84, 56-68.
看楼上有位同学或老师回复过这个问题,他用的CDRFG方法,这个方法是产生inhomogeneous的inflow的,所以在入口处不满足divergence-free,而且文章中感觉有一处错误(也许是我理解不够吧),所以就没有用这个方法。
在入口处加尖劈或者粗糙元,我个人觉得不太好控制想要的风场参数(如紊流积分尺度等),而且由此导致的网格加密会增加计算负担,虽然DSRFG方法也是很耗费计算资源。
-
@霜染丹枫 yes i have solved this problem. Firstly, which method are you using to generate the inflow? Is this method sensitive to the grid sizes?
Sorry for replying you in English because I don't have Chinese input in my desktop.
-
@ibelief 你这个是pt站吗?我现在用的是hdsky,然后让我国内上班的同学用公司的服务器给我搭了个梯子,utorrent挂代理下载,哈哈
-
@liujunCFD 在 一套鼠标控制多个电脑的软件,MARK 中说:
@xjwang 是只要将两台电脑的ip输入进去就可以吗?
我笔记本用的无线,工作站插的网线,能用不
yes you can! 就是输入一个电脑的名字和ip 至于ip是多少,上边会有显示的
-
@李东岳 种子网站 www.rarbg.com,外语资源比较全,电影电视剧游戏等都有,就是刚进去的时候,点一下,会弹出一两个广告界面,但没啥病毒之类的。种子下载好了,可以传到百度云中然后离线下载,所以下载速度快不快就看有没有充会员了,很长时间没用过迅雷了,不知道因为版权问题会不会迅雷禁止下载,但可以用utorrent这种下载软件下载。这个网址挺神奇的,就是成人内容有点多。。。。。
-
@Vortex 非常感谢!!!
-
@Vortex 上边说的损失到一定程度后不变不够准确,因为画在对数坐标中,看起来不那么明显。
-
@Vortex 感谢大佬回复!我是在模拟大气边界层流!我现在是平均风速能拟合地很好,但是脉动成分有所降低,尤其是垂直风向的spanwise和vertical方向,损失最为严重,但发现它们损失到一定程度后就不变了。这个压力项应该是顺风向的压力吧,我在两篇论文中看到过,是修改了求解算法,就是pisoFoam或者pimpleFoam这种,加了一个恒定的力。不知道这个压力项的大小怎么取值?能不能给一些参考?还有类似修改边界条件代码的,能不能也给一些参考?非常感谢您的回复!!!
-
@李东岳 只要网络比较好,用起来没有任何卡顿。最近在家办公,把windows的电脑扔办公室了,所以不记得能不能两台电脑之间复制文件了,我记得应该是可以的。不过东哥都推荐坚果云了,想复制文件,最起码可以拖到坚果云里边,我在国外一直用dropbox。
-
@aiweimo setfields 应该是定义cell centers 的值吧(这个我不是太确定,如果错误希望各位大神指正),这个应该是定义的入口边界处的,还有那个writeCellCentres现在应该是在postProcessing中了,不能直接用。如果你想定义随时间改变的inlet,可以参考一下timeVaringMappedFixedValue(也许拼错了哈哈)。具体实现的例子,可以参考一个叫eddylicious的工具包,或者查看openfoam的代码,或者去cfdonline上查看对应的例子。
-
大家好,最近遇到一个问题,我用了某种方法(DSRFG,Huang et al,2010)产生了入口风速,然后用openfoam中的timeVaryingMappedFixedValue作为入口的边界类型,把数据写在constant/boundaryData/inlet文件中,这样就可以读取数据了。但是这个边界类型,只能把数据定义在入口边界的各控制体的面中心上。其中有个offset选项,可以定义在平移后的一个界面上,但是这样的话,在入口处就没有风速了(不知道我这样理解对不对)。我现在想,能不能有什么方法,对入口边界面和入口边界的控制体中心同时定义速度。这样做的原因是,我发现,如果只是把速度定义在面上的话,那横风向和竖风向的脉动成分不能很好地传递进去。
还有一个常识问题,通过分析顺风向不同位置的脉动风功率谱,发现沿着顺风向,功率谱密度不断减小,这个应该是正常的吧(当然能量损失和采用的离散方式有关,但存在一定的损失应该符合物理现象吧)?具体就像这样,定义两个入口边界。不知道能否实现?
-
推荐一个软件 barrier,可以同时控制linux和windows。之前用mousewithoutborders 但是只能在windows系统之间的电脑用。
-
补充一点:在这个问题中,我删掉了0.17s的数据,然后让pisoFoam 从0.168s开始运行,可以再运行一会,但是又在0.196s卡住了,卡住的时候,courant number平均位0.2255,最大为0.348.
十分感谢!
-
在HPC上运行,一开始用一个节点16个处理器运行的时候,花了12个消失得到了3.6s的数据,后来为了提高速度,就换了16个节点一共256个处理器。
运行了两次,一次在0.163秒卡住了,一次在0.17秒卡住了。没有报错,只是卡在那里不动,查看了一下各个处理器还是在满载运行。
一共有560万个单元,入口风速的边界条件是自己写的,能产生随机风俗时程,之前没有加密沿高度向的底部单元的时候,没有出现过这样的问题。
我现在觉得问题可能出现在:1.在那一个具体时间步进入死循环出不来了,可能引起的原因是不收敛了,但是那一步的courant number还是很小。2.我写的入口风速有问题。
感谢各位能提供一些信息!
如下是运行到最后时候的
Time = 0.168 Courant Number mean: 0.2252957 max: 0.3151702 Patch: inlet recycled 3 eddies mass flow correction coefficient: 1.005238 smoothSolver: Solving for Ux, Initial residual = 0.005431484, Final residual = 3.443522e-07, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.00986471, Final residual = 6.099746e-07, No Iterations 4 GAMG: Solving for p, Initial residual = 0.1917015, Final residual = 0.01027828, No Iterations 2 time step continuity errors : sum local = 8.542895e-07, global = 6.909716e-10, cumulative = -8.085466e-07 GAMG: Solving for p, Initial residual = 0.01235708, Final residual = 8.729178e-07, No Iterations 21 time step continuity errors : sum local = 7.655085e-11, global = -5.037762e-13, cumulative = -8.085471e-07 ExecutionTime = 440.13 s ClockTime = 426 s Time = 0.169 Courant Number mean: 0.2252872 max: 0.3151045 Patch: inlet recycled 4 eddies mass flow correction coefficient: 1.005221 smoothSolver: Solving for Ux, Initial residual = 0.004916222, Final residual = 3.063116e-07, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.009755776, Final residual = 6.022459e-07, No Iterations 4 GAMG: Solving for p, Initial residual = 0.1274995, Final residual = 0.008578488, No Iterations 2 time step continuity errors : sum local = 7.616146e-07, global = -1.088968e-10, cumulative = -8.08656e-07 GAMG: Solving for p, Initial residual = 0.009608369, Final residual = 8.486479e-07, No Iterations 21 time step continuity errors : sum local = 7.726608e-11, global = 7.693458e-13, cumulative = -8.086552e-07 ExecutionTime = 442.66 s ClockTime = 429 s Time = 0.17 Courant Number mean: 0.2252813 max: 0.3153122 Patch: inlet recycled 7 eddies mass flow correction coefficient: 1.00542
-
没有报错,courant number 0.225.
OpenFOAM外挂libtorch
运行pisoFoam的时候,卡住不动了
reconstructPar合并问题
reconstructPar合并问题
OpenFOAM怎么保存指定不等间距时刻的结果
钝体建筑扰流的大涡模拟
钝体建筑扰流的大涡模拟
钝体建筑扰流的大涡模拟
钝体建筑扰流的大涡模拟
OpenFOAM怎么保存指定不等间距时刻的结果
OpenFOAM怎么保存指定不等间距时刻的结果
LES定义入口速度的问题(DSRFG方法)
OpenFOAM 并行计算renumberMesh 报错
OpenFOAM 并行计算renumberMesh 报错
OpenFOAM 并行计算renumberMesh 报错
vanDirest模型及边界条件设置问题
vanDirest模型及边界条件设置问题
介绍一种新的大气边界层入口生成方法DFSR
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
不在学校的话你们都在哪里看电影
一套鼠标控制多个电脑的软件,MARK
不在学校的话你们都在哪里看电影
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
一套鼠标控制多个电脑的软件,MARK
【分享+搬运】自定义非均匀inlet U
LES定义入口速度的问题(DSRFG方法)
LES定义入口速度的问题(DSRFG方法)
一套鼠标控制多个电脑的软件,MARK
运行pisoFoam的时候,卡住不动了
运行pisoFoam的时候,卡住不动了
运行pisoFoam的时候,卡住不动了