@波流力 网格差不多长这样

Latest posts made by wuyukai
-
RE: 关于nOuterCorrectors影响计算速度和计算精度的问题
@cccrrryyy 感谢感谢,我在正式计算之前,算了单个圆柱的,和Oseen公式推出来的理论解能够比上。
您说的方法我试试,由于OpenFoam我也只是初学,没有特别系统的学过,您说的relaxationFactors 我还不太清楚,我去好好了解一下。 -
关于nOuterCorrectors影响计算速度和计算精度的问题
计算一个简单二维多圆柱绕流的动网格算例,网格数量大概40-50万,fvsolution中的nOuterCorrectors和nCorrecters我尝试了不同的值,发现随着该值设置的越来越大,结果才会收敛,最终分别等于20和10。但是这个条件下,计算时间特别久,一个二维40-50万网格的算例要计算将近90个小时,感觉很不合理。想请教一下各位大佬,我是不是哪里设置出错了?
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2006 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { "(p|pcorr)" { solver PCG; preconditioner DIC; tolerance 1e-04; relTol 0.1; } "(p|pcorr)Final" { $p; relTol 0; } U { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-04; relTol 0.1; } UFinal { $U; tolerance 1e-04; relTol 0.0; } } PIMPLE { momentumPredictor no; correctPhi yes; nOuterCorrectors 20; nCorrectors 10; nNonOrthogonalCorrectors 1; ddtCorr true; pRefPoint (20.0 0.0 0.5); pRefValue 0.0; } relaxationFactors { fields { } equations { "U.*" 1; } } // ************************************************************************* //
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2006 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default none; div(U) Gauss linear; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss linearUpwind grad(U); div(phi,K) Gauss linearUpwind grad(U); div(phi,omega) Gauss linearUpwind grad(U); div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } // ************************************************************************* //
-
RE: 请问pointwise 生成的网格,保存成什么格式,可以使用openfoam的网格转换工具?
@队长别开枪 您好,我想请问一下,pointwise导入openfoam时,如何设置网格各个面的boundary类型呢
-
网格整个运动的边界条件问题
运用pimpleFoam,求解一个平板的运动问题,想让整个网格一起运动,程序顺利运行了,但是结果明显是错误。猜测是边界条件设置的问题。想请教一下各位大佬,这种情况下,边界条件如何设置?或者是否有其他方面的问题。
-
RE: reconstructPar合并过程出现Segmentation fault是什么原因?
补充记录一下今天的发现,只有在合并压力p的时候会出现段错误,合并其他物理量的时候则不会。
-
RE: foam-extend 如何对特定区域采用特定的分区方法?
我在尝试使用foam-extend4.0的浸入边界法时,并行也总是出错。后来改用foam-extend4.1的浸入边界法pimpleDyMIbFoam,分块并行的时候没问题,reconstructPar的时候出现段错误