@stardust 用过,网格量上去的话,确实算不动。特别慢
好好科研 发布的最新帖子
-
RE: 在Openfoam中新加入weno格式
@shrine 时间太长了,我有些找不到当时调通的文件啦。印象比较深刻的是,天河超算上不好调通,我是在ubuntu上的openfoam里调通的,你可以考虑先在ubuntu上Allwmake一下。
-
rhoPimpleFoam计算飞机跑几步就发散
各位大佬,我参考文献 3rd High-Lift Workshop Summary Paper - OpenFOAM, STAR-CCM+ & LAVA simulations on Unstructured Grids链接https://arc.aiaa.org/doi/10.2514/6.2018-1253。用rhoPimpleFoam计算飞机,但是算了几步就发散了,希望各位大佬帮忙看看是fvSchemes和fvSolution的问题嘛,还是nut、nuTilda设置的问题。
目前用的是作者上传到会议网站的网格,按照文章里说的,采用的是SA模型,以及动量和湍流项采用的是二阶迎风。
Courant Number mean: 1.4881475e-08 max: 1.8237221 Time = 3.2e-09 * * PIMPLE: Iteration 1 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for Ux, Initial residual = 2.421018e-11, Final residual = 2.421018e-11, No Iterations 0 DILUPBiCGStab: Solving for Uy, Initial residual = 2.1538356e-10, Final residual = 2.1538356e-10, No Iterations 0 DILUPBiCGStab: Solving for Uz, Initial residual = 2.0979924e-10, Final residual = 2.0979924e-10, No Iterations 0 DILUPBiCGStab: Solving for e, Initial residual = 4.1760218e-08, Final residual = 2.1963671e-13, No Iterations 1 GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844072e-07, Final residual = 7.5844072e-07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 5.0250634e-11, global = 5.0210593e-11, cumulative = 3.7505015e-09 GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0 GAMG: Solving for p, Initial residual = 7.5844071e-07, Final residual = 7.5844071e-07, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 5.0250633e-11, global = 5.0210593e-11, cumulative = 3.8007121e-09 PIMPLE: Iteration 2 DILUPBiCGStab: Solving for Ux, Initial residual = 1.3321419e-08, Final residual = 1.002751e-14, No Iterations 1 DILUPBiCGStab: Solving for Uy, Initial residual = 6.2117311e-09, Final residual = 6.2117311e-09, No Iterations 0 DILUPBiCGStab: Solving for Uz, Initial residual = 1.4225061e-10, Final residual = 1.4225061e-10, No Iterations 0 DILUPBiCGStab: Solving for e, Initial residual = 1.5129023e-07, Final residual = 2.4037934e-13, No Iterations 1 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 3.900971e-09 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 GAMG: Solving for p, Initial residual = 1.5142431e-06, Final residual = 1.5142431e-06, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.0032647e-10, global = 1.0025895e-10, cumulative = 4.00123e-09 PIMPLE: Iteration 3 DILUPBiCGStab: Solving for Ux, Initial residual = 2.4133e-08, Final residual = 2.8587861e-14, No Iterations 1 DILUPBiCGStab: Solving for Uy, Initial residual = 1.1420948e-08, Final residual = 1.9949995e-14, No Iterations 1 DILUPBiCGStab: Solving for Uz, Initial residual = 1.4399722e-10, Final residual = 1.4399722e-10, No Iterations 0 DILUPBiCGStab: Solving for e, Initial residual = 1.3956424e-07, Final residual = 2.3199397e-13, No Iterations 1 [73] #0 Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found addr2line failed [73] #1 Foam::sigFpe::sigHandler(int)sh: addr2line: command not found addr2line failed [73] #2 ?sh: addr2line: command not found addr2line failed [73] #3 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate()sh: addr2line: command not found addr2line failed [73] #4 Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: addr2line: command not found addr2line failed [73] #5 ?sh: addr2line: command not found addr2line failed [73] #6 __libc_start_mainsh: addr2line: command not found addr2line failed [73] #7 ?sh: addr2line: command not found addr2line failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x3c75cd0, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=22, req_array=0x3966230, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=52, req_array=0x304c450, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=46, req_array=0x355b680, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=36, req_array=0x3dfabb0, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=40, req_array=0x349e540, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=26, req_array=0x3bd5ce0, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=12, req_array=0x3f18790, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=18, req_array=0x3f97b30, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=14, req_array=0x345ac20, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=32, req_array=0x3df8960, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x4316400, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=28, req_array=0x2fe0f30, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=50, req_array=0x3a36c60, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x3f6c5f0, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed Fatal error in PMPI_Waitall: Unknown error class, error stack: PMPI_Waitall(405).....................: MPI_Waitall(count=24, req_array=0x520edf0, status_array=0x1) failed MPIR_Waitall_impl(181)................: fail failed MPIDU_Complete_posted_with_error(1710): Process failed slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:15 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 *** slurmd[cn10696]: *** STEP 17962804.0 KILLED AT 2021-12-21T17:17:16 WITH SIGNAL 9 ***
fvSchmes
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; limited cellLimited Gauss linear 1; grad(U) $limited; grad(k) $limited; grad(omega) $limited; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U);//Gauss linearUpwind limited; turbulence Gauss linearUpwind default;//Gauss linearUpwind limited; energy Gauss linearUpwind default;//Gauss linearUpwind limited; div(phi,k) $turbulence; div(phi,omega) $turbulence; div(phi,e) $energy; div(phi,K) $energy; div(phi,Ekp) $energy; div(phi,nuTilda) $energy; div(phiv,p) Gauss upwind; div((phi|interpolate(rho)),p) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited 0.5;//Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default limited 0.5;//corrected; } wallDist { method meshWave; }
fvSolution
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; smoother GaussSeidel; tolerance 1e-5; relTol 0.01; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; ncellsInCoarsetLevel 1000; mergeLevels 1; } pFinal { $p; relTol 0; } "(rho|U|k|omega|e)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-8; relTol 0.0; } "(rho|U|k|omega|e|nuTilda)Final" { $U; relTol 0; } } SIMPLE { residualControl { p 1e-5; U 1e-5; "(k|omega|e)" 1e-5; } nNonOrthogonalCorrectors 3; pMinFactor 0.1; pMaxFactor 2; } PIMPLE { nCorrectors 2; nNonOrthogonalCorrectors 3; nOuterCorrectors 3; pMinFactor 0.1; pMaxFactor 2; } relaxationFactors { p 0.3; U 0.7; k 0.7; omega 0.7; } // ************************************************************************* //
nut
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { freestream { type calculated; value uniform 0; } wall { type nutkWallFunction; type nutkWallFunction; value uniform 0; } #includeEtc "caseDicts/setConstraintTypes" } // ************************************************************************* //
nuTilda
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { freestream { type inletOutlet; inletValue $internalField; value $internalField; } wall { type fixedValue; value uniform 0; } symmetry { type symmetry; } #includeEtc "caseDicts/setConstraintTypes" } // ************************************************************************* //
-
用fluent3DMeshToFoam转化网格的时候报错
请教一下,我在icem保存用acsii格式保存了msh文件,用fluent3DMeshToFoam转化网格时出现如下报错,不知道是什么原因呢,操作是在天河二号上进行的。
patch 15 from Fluent indices: 59793665 to: 59902377 type: wall #0 Foam::error::printStack(Foam::Ostream&)sh: addr2line: command not found addr2line failed #1 Foam::sigSegv::sigHandler(int)sh: addr2line: command not found addr2line failed #2 ?sh: addr2line: command not found addr2line failed #3 Foam::polyTopoChange::renumberCompact(Foam::List<int> const&, Foam::List<int>&)sh: addr2line: command not found addr2line failed #4 Foam::polyTopoChange::compact(bool, bool, int&, Foam::List<int>&, Foam::List<int>&)sh: addr2line: command not found addr2line failed #5 Foam::polyTopoChange::compactAndReorder(Foam::polyMesh const&, bool, bool, bool, int&, Foam::Field<Foam::Vector<double> >&, Foam::List<int>&, Foam::List<int>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::objectMap>&, Foam::List<Foam::Map<int> >&, Foam::List<int>&, Foam::List<int>&, Foam::List<Foam::Map<int> >&)sh: addr2line: command not found addr2line failed #6 Foam::polyTopoChange::changeMesh(Foam::polyMesh&, bool, bool, bool, bool)sh: addr2line: command not found addr2line failed #7 ?sh: addr2line: command not found addr2line failed #8 __libc_start_mainsh: addr2line: command not found addr2line failed #9 ?sh: addr2line: command not found addr2line failed
-
RE: 用rhoPimpleFoam计算飞机外流,为什么iteration为0
@number44 我这个p和U的tolerance分别设置的是1e-8和1e-9,感觉挺小的,还要调的更小嘛
-
用rhoPimpleFoam计算飞机外流,为什么iteration为0
各位大佬,想问一下我计算一个飞机模型的外流,设置的dt为2e-8,但是不知道为什么从最开始就是No Iterations 0,这是没在计算嘛?刚开始设置的dt为5e-8,计算一段时间就发散了。
Courant Number mean: 7.5311001e-06 max: 0.56432736 Time = 1.2e-07 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCGStab: Solving for Ux, Initial residual = 6.9178663e-11, Final residual = 6.9178663e-11, No Iterations 0 DILUPBiCGStab: Solving for Uy, Initial residual = 2.5448919e-10, Final residual = 2.5448919e-10, No Iterations 0 DILUPBiCGStab: Solving for Uz, Initial residual = 4.2785621e-10, Final residual = 4.2785621e-10, No Iterations 0 DILUPBiCGStab: Solving for e, Initial residual = 9.7719937e-09, Final residual = 1.0109782e-12, No Iterations 1 GAMG: Solving for p, Initial residual = 7.8947472e-09, Final residual = 7.8947472e-09, No Iterations 0 GAMG: Solving for p, Initial residual = 7.8947472e-09, Final residual = 7.8947472e-09, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.0995375e-11, global = -1.6176582e-13, cumulative = -8.4645739e-13 GAMG: Solving for p, Initial residual = 7.8947534e-09, Final residual = 7.8947534e-09, No Iterations 0 GAMG: Solving for p, Initial residual = 7.8947534e-09, Final residual = 7.8947534e-09, No Iterations 0 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 1.0995384e-11, global = -1.6176582e-13, cumulative = -1.0082232e-12 DILUPBiCGStab: Solving for nuTilda, Initial residual = 9.7105217e-10, Final residual = 9.7105217e-10, No Iterations 0
-
RE: OpenFOAM后处理问题:提取机翼截面位置的压力分布
@五好青年 用tecplot,先slice后extract,就可以了。具体参考这个https://www.bilibili.com/video/BV1Mr4y1w71c/
-
RE: snappyHexMesh生成三维外形添加边界层的相关问题?
@static 你好,我最近也在做三维飞机的相关数值模拟,可否给个联系方式,请教一些问题。谢谢,这是我的qq:3310792254
-
计算添加的方程想用更小的时间步,为什么不能在C文件下改动
希望各位大佬指点一下,就是我想用双时间步,应该怎么改呢。
如下图代码,Info出来的deltaT还是controldict里的deltaTrunTime.deltaT() =0.00001 * demensiont; Info<< "deltaTime = " << runTime.deltaT() << nl << endl; fvScalarMatrix FEqn ( fvm::ddt(H) + fvm::div(phi,H) - fvm::laplacian(DH,H) );