@15735420621 把入口温度条件改成固定温度了，如果用zerogradien的话会在入口边界上产生巨大的温度梯度，导致算例崩溃，但是改成固定值后会稍微好点，但是也会产生负温度，但是将库郎数调整到0.75以下就可以算过去。
H
HIT_lly
@HIT_lly
帖子
-
多相流计算出现局部高压点 -
多相流计算出现局部高压点非常感谢东岳老师的回复，这个算例可能有点大，现在换成固定值入口条件了，到现在还没挂，还正在计算，还算稳定没有出现局部异常情况。
-
openfoam重叠网格出现问题本人在使用openfoam自带的重叠网格进行计算。
按理来说在设置初场的时候两套网格重叠的区域是不参与计算的，但是本人在计算过程中，重叠网格区域好像也参与了计算，本人想就这个问题问一下各位大佬。
这个是被重叠的网格：
这个是背景网格：
![b6f02d62-3983-441a-9196-ca13443e58db-fb193eac98c7188e6dfc2a67c5ec6cf.png]
这个是重叠后的网格
这个是本人设置的初始条件：
由于初始条件只能给两个网格，并且重叠网格重叠区域不参与计算，所以就给了一个圆柱形的高温高压的初始条件。计算过程中总是因为局部速度过大导致计算发散，这是计算结果。
这是没有进行Threshold的计算结果
这个是分离过后的结果
可以看出未进行后处理分离的结果出现了特别大的速度导致了计算的不稳定性，但是按理来说，重叠区域的网格是不参与计算的，上面的结果显示他好像参与计算，有没有大佬知道是为什么。
这个是这个算例的设置条件controlDict
dynamicMeshDict fvSchemes fvSolution
计算的求解器是自己仿写的，但是用他们自带的overcompressibleinterDyMFoam也是这样的结果，有没有大佬知道是为什么。 -
多相流计算出现局部高压点@15735420621 没有调试成功，我换了入口的边界调节入口由之前的零梯度换成了固定值，这样好像会好一点，但是也不知道会在什么时候发散，感觉很神奇。
-
root账户下不可以并行计算（使用不了mpirun）
如图，本人想在root账户下使用mpirun但是ubuntu不建议在root下使用mpirun，在网上找的方法说是要添加--allow-run-as-root进行强制执行，但是发现好像没有解决问题，有没有大神知道该怎么解决这个问题呢 -
OpenFOAM并行核数出现问题@wangfei9088 谢谢大佬能用了
-
OpenFOAM并行核数出现问题
我的电脑逻辑处理器是32个核数，但在使用并行时，如果使用30个核的时，就会提示核数不够的信息，如图所示。使用24个核数的时候会正常运行。之前在使用fluent的时候使用30个核数进行并行没有问题。有没有大佬知道是什么原因。
-
wsl安装完ubuntu但是文件资源管理器中不显示ubuntu@Prometheus10 谢谢您的回复，我找到原因了，原来我之前没有进入mnt文件夹，所以看不到自已本来的硬盘，进入mnt后就能正常使用了。
-
多相流计算出现局部高压点@李东岳 谢谢东岳老师的回复，我再重新调试调试
-
多相流计算出现局部高压点本文在使用compressibleInterIsoFoam进行过膨胀喷管的水下射流计算。开始计算较为稳定无异常情况发生。但是在换个电脑（从另一个电脑copy算例到其他电脑上）进行计算的时候出现了不正常的现象。在射流出口出现了不正常的局部高压，导致计算不稳定，报错内容是计算出现负温度。
这是报错内容
这是出现局部高压过程的压力云图。
时间分别是21.2ms，21.25ms，21.3ms
该算例的fvScheme文件和fvSolution文件如下。/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2212 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; limited cellLimited Gauss linear 1; grad(U) $limited; grad(k) $limited; grad(omega) $limited; } divSchemes { div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss upwind; div(phi,thermo:rho.water) Gauss upwind; div(phi,thermo:rho.air) Gauss upwind; div(Phirb,alpha) Gauss interfaceCompression; div(rhoPhi,U) Gauss limitedLinearV 1; div(rhoPhi,T) Gauss upwind; div(rhoPhi,K) Gauss upwind; div(phi,p) Gauss upwind; div(phi,k) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; default Gauss upwind; div(phi,U) Gauss linearUpwind limited; energy Gauss linearUpwind limited; div(phi,e) $energy; div(phi,K) $energy; div(phi,Ekp) $energy; turbulence Gauss linearUpwind limited; div(phi,k) $turbulence; div(phi,omega) $turbulence; div(rhoPhi,omega) $turbulence; div(rhoPhi,k) $turbulence; div(phiv,p) Gauss upwind; div((phi|interpolate(rho)),p) Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } // ************************************************************************* //这是fvSolution文件
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2212 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { "alpha.water.*" { nAlphaCorr 1; nAlphaSubCycles 5; cAlpha 1; reconstructionScheme isoAlpha; vof2IsoTol 1e-8; surfCellTol 1e-8; nAlphaBounds 3; snapTol 1e-12; clip true; minIter 1; } "pcorr.*" { solver GAMG; preconditioner { preconditioner GAMG; tolerance 1e-12; relTol 0.00001; nSweeps 4; smoother DICGaussSeidel; } tolerance 1e-12; relTol 0.00001; nSweeps 2; maxIter 100; smoother GaussSeidel; minIter 1; } ".*(rho|rhoFinal)" { solver diagonal; minIter 1; } p_rgh { solver GAMG; tolerance 1e-14; relTol 0.00000001; nSweeps 4; smoother GaussSeidel; minIter 2; } p_rghFinal { solver GAMG; preconditioner { preconditioner GAMG; tolerance 1e-14; relTol 0.0000001; nSweeps 4; nVcycles 2; smoother DICGaussSeidel; nPreSweeps 2; } tolerance 1e-14; relTol 0.0000001; nSweeps 3; smoother GaussSeidel; maxIter 100; minIter 2; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-10; relTol 0.00001; nSweeps 2; minIter 1; } UFinal { solver smoothSolver; smoother GaussSeidel; tolerance 1e-10; relTol 0.00001; nSweeps 2; minIter 1; } "(T|k|B|nuTilda).*" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-12; relTol 0.00001; nsweeps 2; limit 10.0; maxIter 100; minIter 1; } p { solver GAMG; smoother GaussSeidel; tolerance 1e-10; relTol 0.0000001; nsweeps 2; minIter 1; } "(U|e).*" { $p; tolerance 1e-9; } "(k|epsilon).*" { $p; tolerance 1e-10; } pFinal { $p; relTol 0.000001; nSweeps 2; minIter 1; } "(rho|U|k|omega|e)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-10; relTol 0.001; nSweeps 2; minIter 1; } "(rho|U|k|omega|e)Final" { $U; relTol 0.00001; nSweeps 2; minIter 1; } } SIMPLE { residualControl { p 1e-8; U 1e-8; "(k|omega|e)" 1e-8; } nNonOrthogonalCorrectors 1; pMinFactor 0.3; pMaxFactor 1; } PIMPLE { momentumPredictor yes; transonic yes; nOuterCorrectors 10; nCorrectors 6; nNonOrthogonalCorrectors 5; nAlphaCorr 1; nAlphaSubCycles 3; CAlpha 1; pMin 10000; pMax 1000000; } relaxationFactors { equations { ".*" 0.1; } } // ************************************************************************* // -
wsl安装完ubuntu但是文件资源管理器中不显示ubuntu
如图，ubuntu已经安装在windows上了，并且已经安装完wsl2，但是想要进行文件互相传递的时候，发现文件系统中找不到“linux”项，有没有懂哥知道这是怎么回事呢？ -
snappyHexMesh画网格之前snappyHexMesh画的3D网格进行重叠网格计算没有问题。后来为了节省计算成本用的2D网格，用的snappyHexMesh画的网格，然后使用extrudeMesh将3D网格变成2D网格后就出现问题了
-
snappyHexMesh画网格@李东岳 感谢东岳老师的回复，我的fvScheme是这样设置的，这用的是重叠网格技术。开始的isoadvector是我耦合的重叠网格算法，然后后面的compressibleinterFoam是openfoam 中自带的算法。
-
snappyHexMesh画网格/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2312 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default Gauss upwind; div(rhoPhi,U) Gauss upwind; div(U) Gauss linear; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss linear; div(phi,thermo:rho.water) Gauss upwind; div(phi,thermo:rho.air) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(rhoPhi,T) Gauss upwind; div(rhoPhi,K) Gauss upwind; div(phi,p) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,U) Gauss linearUpwind limited; energy Gauss linearUpwind limited; div(phi,e) $energy; div(phi,K) $energy; div(phi,Ekp) $energy; turbulence Gauss linearUpwind limited; div(phi,k) $turbulence; div(phi,omega) $turbulence; div(rhoPhi,omega) $turbulence; div(rhoPhi,k) $turbulence; div(phiv,p) Gauss upwind; div((phi|interpolate(rho)),p) Gauss upwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } oversetInterpolation { method inverseDistance; holeLayers 6; useLayer 3; } oversetInterpolationSuppressed { grad(p_rgh); surfaceIntegrate(phiHbyA); } fluxRequired { default no; p_rgh; pcorr; alpha.water; } // ************************************************************************* // -
snappyHexMesh画网格@李东岳 东岳老师，上面那个报错是因为使用了isoAdvector算法。
使用interFoam进行计算也会进行报错，报错是这样的FOAM Warning : From Foam::triFace Foam::tetIndices::faceTriIs(const Foam::polyMesh&, bool) const in file meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H at line 70 No base point for face 7847, 4(5577 25968 25969 5578), produces a valid tet decomposition. -
snappyHexMesh画网格@李东岳 谢谢东岳老师，这个我用interFoam也会报错，应该是snappyHexMesh切分网格的时候产生的问题。
-
snappyHexMesh画网格使用snappyHexMesh画网格时，画出网格时没有什么问题但是在计算的时候会报错。有没有大佬知道这种情况该怎么解决
FOAM Warning : From void Foam::cutFaceAdvect::cutPoints(Foam::label, Foam::scalar, Foam::DynamicList<Foam::Vector<double> >&) in file cellCuts/cutFace/cutFaceAdvect.C at line 927 cutPoints = 4((0 0.0449362922089 0.385062422146) (0 0.0436853188202 0.385062441697) (0.01 0.0436853188202 0.385062441697) (0.01 0.0449362922089 0.385062422146)) for pts = 4((0 0.0449362922089 0.385062422146) (0 0.0436853188202 0.385062441697) (0.01 0.0436853188202 0.385062441697) (0.01 0.0449362922089 0.385062422146)), f - f0 = 4(23405 23403 38935 38937) and f0 = 2.22618323362e-13 -
snappyHexMeshDict画二维网格出现锯齿状的区域问了一下别人，好像是使用snappyHexMesh命令画2D网格就是会出现面贴合效果不好的情况
-
snappyHexMeshDict画二维网格出现锯齿状的区域这个是所有的文件内容，请各位大佬指导一下
floatingBody.zip -
snappyHexMeshDict画二维网格出现锯齿状的区域@lalalala 我是用的snappyHexMesh -overwrite命令，好像是直接覆盖了