CFD中文网

本人在使用openfoam自带的重叠网格进行计算。
按理来说在设置初场的时候两套网格重叠的区域是不参与计算的，但是本人在计算过程中，重叠网格区域好像也参与了计算，本人想就这个问题问一下各位大佬。
这个是被重叠的网格：
6384ab7c-e651-4dc9-9bd2-16298be04a95-432b47d32d1e3106b90bbfeb518ec39.png
这个是背景网格：
e2221afd-ff60-4bdc-b424-758b0087efb0-8d2e870096a5f44c86bbcf28b4318c3.png
![b6f02d62-3983-441a-9196-ca13443e58db-fb193eac98c7188e6dfc2a67c5ec6cf.png]
这个是重叠后的网格
4e960bd4-99a3-45d3-8c01-8f2c8888dabf-fb193eac98c7188e6dfc2a67c5ec6cf.png
这个是本人设置的初始条件：
b57ddbe1-1225-4cc4-847d-243257530ba2-913b4b3992b0d17eedda2894abcc303.png
由于初始条件只能给两个网格，并且重叠网格重叠区域不参与计算，所以就给了一个圆柱形的高温高压的初始条件。

计算过程中总是因为局部速度过大导致计算发散，这是计算结果。
这是没有进行Threshold的计算结果
777b8258-58a1-41ae-b84d-c6775eddfaca-4a3190545fd101c6ad980a5cbbf4455.png
这个是分离过后的结果
5fd95244-6427-43a8-a9d9-c8c1a0bb9b6f-bf547f2742f3a21d4d4fd8e035395ec.png
可以看出未进行后处理分离的结果出现了特别大的速度导致了计算的不稳定性，但是按理来说，重叠区域的网格是不参与计算的，上面的结果显示他好像参与计算，有没有大佬知道是为什么。
这个是这个算例的设置条件controlDict
dynamicMeshDict fvSchemes fvSolution
计算的求解器是自己仿写的，但是用他们自带的overcompressibleinterDyMFoam也是这样的结果，有没有大佬知道是为什么。

试试：mpirun --allow-run-as-root -np 24 compressibleInerIsoFoam

@wangfei9088 谢谢大佬能用了:xinxin:

@HIT_lly 好的，谢谢！我试试

@Prometheus10 谢谢您的回复，我找到原因了，原来我之前没有进入mnt文件夹，所以看不到自已本来的硬盘，进入mnt后就能正常使用了。

之前snappyHexMesh画的3D网格进行重叠网格计算没有问题。后来为了节省计算成本用的2D网格，用的snappyHexMesh画的网格，然后使用extrudeMesh将3D网格变成2D网格后就出现问题了

问了一下别人，好像是使用snappyHexMesh命令画2D网格就是会出现面贴合效果不好的情况

OpenFOAM重叠网格技术是不是不支持半模计算（其中有一个面是对称面）。

如图：
0eea5f23-83c2-458a-ba95-518789118336-95b7bc47d4925bb22540e673dda9fd11.JPG
本人想把这个六面体的一个侧面当成对称面，但是无论在blockMesh中选择任意侧面当成对称类型，当执行完blockMesh命令后，总是顶面被选择为对称面。
30dbbabe-3fa4-4563-aea3-80c2c3941a37-d1e92b73bb6e50611f2d95dc06abf563.JPG

这是blockMeshDict文件

/*--------------------------------*- 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 blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // scale 1; vertices ( (-0.2 -0.2 -0.2) (1.2 -0.2 -0.2) (1.2 0.3 -0.2) (-0.2 0.3 -0.2) (-0.2 -0.2 2) (1.2 -0.2 2) (1.2 0.3 2) (-0.2 0.3 2) ); blocks ( hex (0 1 2 3 4 5 6 7) (80 80 70) simpleGrading (1 1 1) ); edges ( ); boundary ( // Dummy patch to trigger overset interpolation before any other bcs oversetPatch { type overset; faces (); } stationaryWalls { type wall; faces ( (1 5 4 0) (3 7 6 2) (0 3 2 1) (2 6 5 1) ); } sys { type symmetry; faces ( (0 4 7 3); ) } atmosphere { type patch; faces ( (4 5 6 7) ); } ); mergePatchPairs ( ); // ************************************************************************* //

请问各位大神现在openfoam自带的重叠网格方法是不是不支持对称面(symmetry)呢

@学流体的小明 嗯嗯，谢谢大佬

@学流体的小明 恩呢，谢谢您的回答。
/--------------------------------- 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;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicRefineFvMesh;

// How often to refine
refineInterval 1;

// Field to be refinement on
field alpha.water;

// Refine field inbetween lower..upper
lowerRefineLevel 0.001;
upperRefineLevel 0.999;

// If value < unrefineLevel unrefine
unrefineLevel 1;

// Have slower than 2:1 refinement
nBufferLayers 6;

// Refine cells only up to maxRefinement levels
maxRefinement 1;

// Stop refinement if maxCells reached
maxCells 15000000;

// Flux field and corresponding velocity field. Fluxes on changed
// faces get recalculated by interpolating the velocity. Use 'none'
// on surfaceScalarFields that do not need to be reinterpolated.
correctFluxes
(
(phi none)
(nHatf none)
(rhoPhi none)
(alphaPhi_ none)
(ghf none)
(phi0 none)
(dVf_ none)
(alphaPhi0.water none)
(alphaPhiUn none)
);

// Write the refinement level as a volScalarField
dumpLevel true;

// ************************************************************************* //
这是我的dynamicMeshDict文件。
我按您说的方法去试一下。

最近在做多相流的仿真验证，但是一直找不到LS方法的多相流案例，想请教以下各位大佬openfoam中有没有关于LS方法的多相流案例

@yhdthu 好的 谢谢大佬