@fkingdom MPI有一个参数叫绑定核心,不同的mpi命令不一样,你可以查一下。绑定之后,可以同时跑2个20核心的算例,并且不互相争抢资源。
段志伟
帖子
-
-
@疏影横斜水清浅 在空气动力学领域,一般是大一个量级。
-
@shrine 没算错,这个算法就是会出现这种情况。可以尝试减小时间步长。
-
LES只能计算瞬态。在湍流计算中,RANS方法全部模化湍流项;LES方法则模化小尺度涡部分,解析大尺度涡部分, 所谓“解析”就是瞬态部分。LES的强制稳态没有意义,因为只计算了很少一部分的模化湍流,而大部分湍流(解析)并没有计算。因此LES的强制稳态结果,可能更接近于强制的层流稳态结果。
-
@cccrrryyy 求解器本身应该没啥问题,我也没弄清为啥foam的效率提高的不多,可能是给的算例里加了relaxation?线性求解器可以用开源库试试,比如petsc。
-
可能是foam本身的问题,也可能是设置问题。从我自己的经验来看,大概couple算法时间是segregate算法的一半左右吧。工业案例失败主要是线性求解器不太好用,需要修改。
-
@疏影横斜水清浅 没区别
-
@fu 实际上URANS计算的总的雷诺应力与真实雷诺应力是有差距的。图上仅仅是URANS能解析的那部分。
-
@fu 用湍流模型计算的时候,雷诺应力会分成两部分:模化+解析。 fieldAverage统计的只是解析的部分,也就是真实流体运动的应力。用Laminar的时候,模化的应力为零,所以显示的就是真实的应力;而用SST的时候,模化的应力占绝大多数,解析部分几乎为零。
-
这个翻译需要看场景吧,可以为“湍流间歇性”,也可以为“湍流间歇因子”。看你植入了一个和SA湍流方程耦合的转捩模型,增加的那个方程就是“间歇因子”输运方程。
-
@jyang
在decomposeParDict文件最后,加入:
constraints
{
constraint1
{
type preservePatches;
patches (region1_to_region2 region2_to_region1);
}
} -
在decompose的时候,可以指定region1和region2的交界面在一个核中,这样应该就没问题了
-
这个pRefCell和pRefValue只在封闭腔体时才会起作用,你的办法是对的,给压力边界值之后就不需要setReference了。我这个问题还比较特殊,不是封闭腔体,但在某几个cell的地方一直会有温度极值,这个值还会持续升高。
-
BoussinesqSimpleFoam计算,开启重力时有局部点温度一直升高,太奇怪了。有遇到过的大神吗?
-
@hotairballon1997 和steady Rans结果不同,因为每个时间步的结果不同。
-
@hotairballon1997 每时每刻都在变化。
-
@麦迪文 在固体区域认为粘性非常大,1E6量级大概,所以不能流动。
-
开源社区里类似的求解,好像是erfConvectiveMeltingFoam,简单的凝固求解器。
另外,fvOptions里有一个solidificationMeltingSource可以用。 -
你这是在耦合计算中没有找到transportProperties这个字典文件。可以仔细看一下,耦合计算和单场计算时,流体域是不是用的相同的求解器。
-
@程迪 大佬可以详细的解释一下cyclic的计算过程中吗?比如这个CL和CU到底对应哪一项,initMatrixInterfaces怎么计算了Cl*x0.
-
@Jacobian 我也是初用OpenFOAM算气动问题。
(1)对于bounded格式,实际上看代码的话,会减去phid对应的场量的散度再乘以p。如果散度为0,那么没有问题。但是phid的散度实际上不为0,所以我倾向于不用bounded。
(2)最小二乘其实是为了保证精度。前面的cellLimited是为了数值稳定性。
我用Simple-C算没问题,但Simple一直发散,目前还没找到原因。
(3)网格其实实用Pointwise做的,SHM不太容易做。 -
@东岳 并没有解决,还发现了新问题。simpleC和simple算的升力系数大概差20%
-
@kdfluxit rhoPimpleFoam里有 include "rhoEqn"这段代码,也就是显式求解rho方程。你再仔细看看。
-
@dzw05 抱歉网站上传附件有大小和格式限制。所以我把压缩文件分成三份。rhoSimple.z01和rhoSimple.z01文件,下载后请去掉zip后缀。rhoSimple.z02.zip
-
理论上rhoSimpleFoam可以计算可压缩流动问题,官方给的tutorial中也有计算超音速流动的算例。但是我在计算一个三段翼的亚音速流动(U=125m/s)时,却无论如何都不能收敛,但是我把速度降到60之后就很容易收敛。不知道各位大佬有何高见?附件是计算的网格和设置文件,各位有兴趣可以试一试。rhoSimple.zip
-
建议用排除法比较。(1)首先计算层流,对比fluent和OpenFOAM结果,这一步基本不会有问题,否则OpenFOAM不用混了。(2)再计算湍流,进行对比。个人觉得可能你的湍流模型可能有问题。
-
@林之流风 我计算的是两相流问题,把感兴趣区域中某一相的粘度设的特别大(比如1e5),该相就相当于凝固了。
-
@win 这不是错误,是连续性方程的误差。
-
@Samuel-Tu 试过,可行。
-
@upc_ngh 如@星星星星晴 所言,parcel就是n个particle的合集,具有相同的参数。
-
@东岳 哈哈,当时东岳老师还在国外。东岳现在哪高就?有机会去拜访您。
-
@张RP 你看看压力文件,即0/p文件,应该是压力的量纲写错了。
-
@张RP 建议用可压缩求解器。
-
不可压缩单相里面没有rho这个变量,所以你目前的质量源项无效。
-
@upc_ngh 可以,DPM和MPPIC就是高浓度粒子的求解器。
-
@upc_ngh MPPIC和DPM好像没有粘连?
-
@Calf-Z-DNS delta Y+的计算公式你已经写出来了,和delta X+是一样的。平板边界层的估算算法:https://www.cfd-online.com/Wiki/Y_plus_wall_distance_estimation
-
delta y就是指法向相邻网格间的距离。DNS计算一般要求法向第一层网格y+<1;delta y+取值大小取决于你想要在边界层内布置多少网格点。如是是计算平板,可以用经验公式估计一下摩擦系数,然后算出y+;如果模型比较复杂,那只能在实际计算后得到y+。
-
@东岳 类似算例我用自己写的代码算过,层流、湍流、网格、离散格式影响,基本都能算出分离,主要还是数值粘性的问题。
-
@QuietCheng 首先确定网格,y+的范围,流向网格分布;其次看一下离散格式,低阶格式数值粘性大,可能也算不出来分离。
-
@红豆沙 本质上这是由于插值引起的,如果自己写代码做计算和后处理展示,可以做到没有(或减少)瑕疵。减小网格尺寸,会让这些瑕疵显得小一些。
-
@小令 你可以在plot—>blanking—>value Blanking中,选择把右侧的数据隐藏掉。
-
@我是河滩 不需要显式设置密度,自己心里知道就行了。
-
@我是河滩 密度可以是任意值,不是1.225。计算雷诺数时,密度体现在粘性里,由miu,变为niu。
-
@zwl paraview,可以不需要重构,直接后处理。
-
@麦迪文 理论上来说,你用pimple算瞬态问题不应该加任何亚松弛(当然为了稳定可以加一点)。
-
@bishuiyouyou 是指允许的最大体积分数,数值上不会再大了。不影响流动。
-
代码中的描述的作用是为了使求解稳定,在原有格式的基础上减去了
fvm::Sp(fvc::surfaceIntegrate(faceFlux), vf)
,而计算收敛的时候该项为0,所以理论上应该不会影响精度。
关于并行计算效率问题
模型尺寸与计算域尺寸之间的关系
interFoam算出来α大于1是什么情况
LES只能算瞬态还是说稳态瞬态都能算?
有人玩儿过耦合求解器pUCoupledFoam么
有人玩儿过耦合求解器pUCoupledFoam么
通量表示的区别
使用SSTkomega和laminar算出来的雷诺应力问题
使用SSTkomega和laminar算出来的雷诺应力问题
intermittency怎么翻译成中文
多区域模型单核计算没问题,但无法并行计算的问题
多区域模型单核计算没问题,但无法并行计算的问题
BoussinesqSimpleFoam开启重力,温度计算发散
BoussinesqSimpleFoam开启重力,温度计算发散
OpenFOM里面的是RANS还是URANS
OpenFOM里面的是RANS还是URANS
请问OF里有没有融化和凝固相关的求解器
请问OF里有没有融化和凝固相关的求解器
流固耦合alphat使用alphatJayatilleWallFunction报错跪求大佬指点
openfoam中cyclic周期性边界的问题
rhoSimpleFoam计算翼型发散
rhoSimpleFoam计算翼型发散
OF可压流求解器
质量源项添加问题,参数rho无效
rhoSimpleFoam计算翼型发散
rhoSimpleFoam计算翼型发散
rhoSimpleFoam计算翼型发散
fluent总能比of得到更满意的结果?
OpenFOAM有方法能够使一部分网格不参与计算吗?
出现时间连续错误提醒
OpenFOAM有方法能够使一部分网格不参与计算吗?
请问有人用过sprayFoam中的下面coneNozzleInjection吗?
质量源项添加问题,参数rho无效
质量源项添加问题,参数rho无效
质量源项添加问题,参数rho无效
质量源项添加问题,参数rho无效
液固两相流
液固两相流
DNS中△y+的意义和计算方法
DNS中△y+的意义和计算方法
关于超声速双压缩拐角流动分离问题
关于超声速双压缩拐角流动分离问题
重叠(嵌套)网格讨论
各位大佬们,在tecplot中,我在使用切片工具的时候,都会把后面那一段切上,如何把红线勾画的这段去掉。谢了
不可压计算雷诺数的计算问题
不可压计算雷诺数的计算问题
并行计算的结果,不进行reconstructPar,用什么软件能整体处理?
耦合换热求解器chtMultiRegionBoussinesqPimpleFoam
颗粒流packing limit
openfoam中bounded和unbounded的区别