好哥哥,问题解决了没有。用硫酸老师的方法,总是报错误,可能是我电脑有什么插件没有安装上。
希望先生
帖子
-
-
@李东岳 好的李老师,受教了
-
@bestucan 嗯嗯,感谢,非常仔细的教程。只可惜手头上的系统是win10系统,后面有机会尝试一下在换成Linux系统。
-
@bestucan 对的,fluent是可以直接读一个journal,这样可以进行大量相似case的批量计算,很方便。但是,经常会有大量完全不一样的case的情况,而且有时候会很多人用一台服务器,这样使用journal的话,就有点力不从心的感觉。
我之前了解到,出现前面同时打开两个case计算速度变慢的原因,好像是和没有绑定CPU有关,就是没有把一个进程绑定到一个CPU核心上,导致在计算的时候,fluent的很多个进程会在不同的CPU核心上来回切换,导致即使没有用满全部的CPU核心,计算两个case,速度也会变慢。但是就是不知道怎样能够实现批量绑定fluent进程到CPU核心上,也不知道BIOS里面有没有相关的设置。 -
手头上有一台56物理核的windows10系统的服务器,8273cl2,内存是16G8,算了两个case,一个case用了40核,另外一个只用了10个核。在打开10核case之前,40核case的计算速度非常非常猛 ,但是打开10核case计算之后,40核case计算速度瞬间就降下来了,而且降了非常多 各位大佬,请问这是怎么回事,应该怎么解决这种问题
-
这里给出了一种方法:fluent导入形成流场的外部结构
-
大家好,如果单纯给第一个图的话,大部分人可能不知道第一个图是什么流场;如果是用第二个图的话,可能就会有一部分人知道这是一个柱塞阀的内部流场,主要是因为加上了形成流场的外部结构。对于这一部分,即便是用fluent计算的流场,也可以在paraview中实现。但是,如果考虑到其他的一些需要考虑壁面的情况,比如复杂结构的壁面上的冲蚀磨损(气力输送阀门中含有固体颗粒的情况),最好是将每一个壁面都在划分网格的时候列为一个part,使用fluent计算完成之后,便可以直接查看每个part壁面的上的冲蚀磨损。但是想要像上面第二个图一样,加上外部结构的话,应该怎么办呢 ?像上面一样,通过paraview进行后处理,但是fluent计算结果导入到paraview之后,似乎是无法对fluent中的每一个part壁面单独显示。这个应该怎么解决呢 ?(OpenFOAM+LIGGGHTS计算冲蚀并通过paraview后处理是完全可以的,并不存在上面所说的问题,但是考虑每一个真实颗粒和碰撞需要非常长的计算时间,时间成本有点高)或者还有什么后处理软件可以对fluent的计算结果进行后处理,并且经过一系列的操作,能够实现显示形成流场的外部结构,让得到的云图不会太抽象,并能显示每个part壁面上的冲蚀磨损 ? -
各位老师,大家好,想请教一下,在拥有管理员账号的前提下,怎样在联想集群的LSF作业调度系统中添加新的普通账户?
之前是按照这样的方式,也忘记了是在哪里找到的了://在管理节点上添加新用户(root): useradd -d /gpfsdata/home/cfdem -m cfdem //在管理节点上给新用户添加密码(root) passwd cfdem //添加到用户组 vim /gpfsdata/apps/lsf/conf/lsbatch/cluster1/configdir/lsb.users //添加完后运行 badmin ckconfig //在计算节点上添加新用户 ssh c01n01 useradd -m cfdem passwd cfdem并输入密码 //同理在c01n02~c01n03同样上述操作 //ssh到管理节点mgt01并切换到cfdem ssh mgt01 su cfdem ssh-keygen//然后一直回车
但是这种方式感觉有点繁琐。
后来了解到PBS管理系统里面可以直接通过命令clusconf -au 账号名
就可以完成普通账户的添加。
现在想问一下,LSF管理系统添加账户的时候,有没有类似于PBS管理系统的这种一气呵成的简单命令? -
@李东岳 东岳老师,这个就已经够用了
-
东岳老师要不给大家整一个收藏功能 让大家把这些牛逼的小窍门给收藏下来,方便查找
-
@李东岳 感谢李老师
-
@李东岳 谢谢李老师提醒,我去尝试一下
-
@李东岳 感谢李老师 ,但是没有T型管的试验数据。原本是要计算一个掺合阀的,原理类似于T型管,但是那个掺合阀模型比这个T型管要复杂很多,尤其在小口径管和大口径管结合的地方,那个模型一直计算不下去就画了个简单的T型管。这个T型管的边界条件,是我按照厂子里面测试掺合阀的实验条件放上的,如果按照他们的试验的话,热流进口是应该大概是在930-1600摄氏度范围内,冷流大概是160-170摄氏度,掺合阀调控的混合流出口温度大概在220-260摄氏度。李老师,厂子里面实验条件非常简陋,仅作小小参考 只要结果不太一样,那就勇敢地相信是他们试验做得不对
-
@李东岳 李老师,会不会还是网格的原因?有可能是我网格划分的有问题
-
@李东岳 李老师,不好意思,没有及时关注回复。我这个算例使用OpenFOAM-5.x,因为需要使用CFDEM,所以一直没有编译安装新版本的OpenFOAM。
-
@wangfei9088 多谢老师提醒,这个地方确实是把INLET_COLD和OUTLET的边界条件写反了,我尝试改了过来,但是仍旧没有解决问题。
-
@hurricane007 多谢老师提醒,我这里并没有特别地去做一些明显的过渡,我去试一下看能不能解决问题。
-
@李东岳 感谢李老师,我不着急的
-
@李东岳 李老师,你好。我已经将网格数量调到了77万,但是上传不到这上面来,我发送到您的邮箱里面去了,不知道符合不符合您的要求。
-
@李东岳 李老师,你好。我上面画的T型管的网格数量大概为360万。现在没有sci数据,原本是需要计算一个类似于T型管的掺合阀,那个掺合阀的模型比这个T型管要复杂的多。 因为一直计算不下去,开始以为是掺合阀用的是非结构网格的原因,于是就先计算一个这样比较简单的T型管,划分了结构网格,试了一下还是计算不下去。现在感觉可能是边界条件的原因,或者是什么其他的原因。
-
各位老师,大家好,我在使用rhoPimpleFoam求解器计算一个T型管冷热空气掺混算例的时候,总是计算几步就停止掉了。之前从来没有计算过这种和热有关的算例,我感觉有可能是边界条件设置的问题,烦请有经验的老师给指导一下。
这是T型管算例的大概示意图
这是0/p文件:dimensions [1 -1 -2 0 0 0 0]; internalField uniform 101325; boundaryField { INLET_HOT { type fixedValue; value uniform 101550; } INLET_COLD { type zeroGradient; } OUTLET { type fixedValue; value uniform 101350; } WALL_HOT { type zeroGradient; } WALL_MID { type zeroGradient; } }
这是0/T文件:
dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { INLET_HOT { type fixedValue; value uniform 1673.15; } INLET_COLD { type fixedValue; value uniform 438.15; } OUTLET { type inletOutlet; value $internalField; inletValue $internalField; } WALL_HOT { type zeroGradient; } WALL_MID { type zeroGradient; } }
这是0/U文件:
dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { INLET_HOT { type pressureInletVelocity; value uniform (0 0 0); } INLET_COLD { type pressureInletVelocity; value uniform (0 0 0); } OUTLET { type pressureInletOutletVelocity; value uniform (0 0 0); } WALL_HOT { type noSlip; } WALL_MID { type noSlip; } }
这是0/alphat文件:
dimensions [1 -1 -1 0 0 0 0]; internalField uniform 1e-3; boundaryField { INLET_HOT { type calculated; value $internalField; } INLET_COLD { type calculated; value $internalField; } OUTLET { type calculated; value $internalField; } WALL_HOT { type compressible::alphatWallFunction; value $internalField; } WALL_MID { type compressible::alphatWallFunction; value $internalField; } }
这是0/epsilon文件:
dimensions [0 2 -3 0 0 0 0]; internalField uniform 2.2667; boundaryField { INLET_HOT { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.005; value uniform 2.2667; } INLET_COLD { type turbulentMixingLengthDissipationRateInlet; mixingLength 0.005; value uniform 2.2667; } OUTLET { type inletOutlet; inletValue uniform 2.2667; value uniform 2.2667; } WALL_HOT { type epsilonWallFunction; value uniform 2.2667; } WALL_MID { type epsilonWallFunction; value uniform 2.2667; } }
这是0/k文件:
dimensions [0 2 -2 0 0 0 0]; internalField uniform 0.054; boundaryField { INLET_HOT { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.054; } INLET_COLD { type turbulentIntensityKineticEnergyInlet; intensity 0.05; value uniform 0.054; } OUTLET { type inletOutlet; inletValue uniform 0.054; value uniform 0.054; } WALL_HOT { type kqRWallFunction; value uniform 0.054; } WALL_MID { type kqRWallFunction; value uniform 0.054; } }
这是0/nut文件:
dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { INLET_HOT { type calculated; value uniform 0; } INLET_COLD { type calculated; value uniform 0; } OUTLET { type calculated; value uniform 0; } WALL_HOT { type nutkWallFunction; value uniform 0; } WALL_MID { type nutkWallFunction; value uniform 0; } }
这是constant/thermophysicalProperties文件:
thermoType { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } mixture { specie { molWeight 28.9; } thermodynamics { Cp 1007; Hf 0; } transport { As 1.4792e-06; Ts 116; } }
这是system/fvScheme文件:
ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss upwind; div(phid,p) Gauss upwind; div(phiv,p) Gauss linear; div(phi,K) Gauss linear; div(phi,e) Gauss upwind; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,R) Gauss upwind; div(phi,omega) Gauss upwind; div((rho*R)) Gauss linear; div(R) Gauss linear; div(U) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; }
这是system/fvSolution文件:
solvers { p { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-05; relTol 0; smoother GaussSeidel; cacheAgglomeration true; } tolerance 1e-05; relTol 0.1; } pFinal { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-06; relTol 0; smoother GaussSeidel; cacheAgglomeration true; } tolerance 1e-06; relTol 0.1; } "(rho|U|e|k|epsilon|omega)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-06; relTol 0; } "(rho|U|e|k|epsilon|omega)Final" { $U; relTol 0; } } PIMPLE { momentumPredictor yes; transonic no; nOuterCorrectors 50; nCorrectors 1; nNonOrthogonalCorrectors 0; consistent yes; SIMPLErho yes; pMaxFactor 1.5; pMinFactor 0.9; residualControl { "(U|k|epsilon)" { relTol 0; tolerance 0.00001; } } turbOnFinalIterOnly no; } relaxationFactors { fields { "p.*" 0.5; "rho.*" 1; } equations { "U.*" 0.9; "e.*" 0.7; "(k|epsilon|omega).*" 0.8; } }
然后,这是运行rhoPimpleFoam时保存的log文件,看不出非常明显的报错信息:
PIMPLE: iteration 25 DILUPBiCGStab: Solving for Ux, Initial residual = 0.188633, Final residual = 6.22293e-07, No Iterations 4 DILUPBiCGStab: Solving for Uy, Initial residual = 0.108362, Final residual = 3.91435e-07, No Iterations 4 DILUPBiCGStab: Solving for Uz, Initial residual = 0.36688, Final residual = 5.86313e-07, No Iterations 4 DILUPBiCGStab: Solving for e, Initial residual = 0.170771, Final residual = 1.74252e-07, No Iterations 3 ------------------------------------------------------- Primary job terminated normally, but 1 process returned a non-zero exit code. Per user-direction, the job has been aborted. -------------------------------------------------------
这是在集群上计算停止后保存的errput文件,不知道是否有用,先放上 :
-------------------------------------------------------------------------- [[18352,1],4]: A high-performance Open MPI point-to-point messaging module was unable to find any relevant network interfaces: Module: OpenFabrics (openib) Host: c01n02 Another transport will be used instead, although this may result in lower performance. NOTE: You can disable this warning by setting the MCA parameter btl_base_warn_component_unused to 0. -------------------------------------------------------------------------- [c01n02:275299] 31 more processes have sent help message help-mpi-btl-base.txt / btl:no-nics [c01n02:275299] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages [16] #0 Foam::error::printStack(Foam::Ostream&)[15] #0 Foam::error::printStack(Foam::Ostream&) at ??:? at ??:? [15] #1 Foam::sigFpe::sigHandler(int)[16] #1 Foam::sigFpe::sigHandler(int) at ??:? [15] #2 ? at ??:? [16] #2 ? in /usr/lib64/libpthread.so.0 [15] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool) in /usr/lib64/libpthread.so.0 [16] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool) at ??:? [15] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? [16] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:? [15] #5 at ??:? [16] #5 ?? at ??:? [15] #6 __libc_start_main at ??:? [16] #6 __libc_start_main in /usr/lib64/libc.so.6 [15] #7 in /usr/lib64/libc.so.6 [16] #7 ?? at ??:? [c01n02:275348:0:275348] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2)) at ??:? [c01n02:275349:0:275349] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2)) ==== backtrace ==== 0 0x000000000000f4ab raise() ???:0 1 0x000000000000f5e0 _L_unlock_13() funlockfile.c:0 2 0x00000000002addcf Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() ???:0 3 0x00000000002aea19 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() ???:0 4 0x00000000004255e6 main() ???:0 5 0x0000000000021c05 __libc_start_main() ???:0 6 0x000000000042d711 _start() ???:0 =================== ==== backtrace ==== 0 0x000000000000f4ab raise() ???:0 1 0x000000000000f5e0 _L_unlock_13() funlockfile.c:0 2 0x00000000002addcf Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() ???:0 3 0x00000000002aea19 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() ???:0 4 0x00000000004255e6 main() ???:0 5 0x0000000000021c05 __libc_start_main() ???:0 6 0x000000000042d711 _start() ???:0 =================== -------------------------------------------------------------------------- mpirun noticed that process rank 16 with PID 0 on node c01n02 exited on signal 8 (Floating point exception). --------------------------------------------------------------------------
烦请有时间、有经验的老师给提一下建议
-
@冠竹 可以在ICEM导出.msh的时候,scale一下;也可以在转化网格的时候,用fluent3DMeshToFoam XXX.msh -scale 0.001缩放一下。
cgns格式下的多面体网格,在paraview中查看
fluent同时打卡两个case计算速度会变的非常慢
fluent同时打卡两个case计算速度会变的非常慢
fluent同时打卡两个case计算速度会变的非常慢
fluent同时打卡两个case计算速度会变的非常慢
后处理问题:添加形成流场的外部结构
后处理问题:添加形成流场的外部结构
联想集群的LSF管理系统添加普通账户
reconstructPar 一边组一边删除数据
reconstructPar 一边组一边删除数据
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
用rhoPimpleFoam计算T型管冷热流混合,计算几步就停掉
OpenFOAM 尺寸疑问