局部加密区域跨两个block，运行topoSet和refineMesh命令会报错

GenhouSUN

最近是用topoSet和refineMesh局部加密。如果加密区域在单个block内可以正常运行，但加密区域在两个内就会出错。
我的疑问：topoSet和refineMesh命令：topoSet 和 refineMesh 选定的区域不能包含两个block对么？如果想对横跨两个block的区域加密，怎么实现呢？是在各自的block内分别加密么？
topoSet和refineMesh 运行都是正确的

两个block信息如下
vertices
(
( 0.0 0.0 -10.0)
( 100.0 0.0 -10.0)
( 100.0 5 -10.0)
( 0.0 5 -10.0)

   ( 0.0 0.0 5.0)
   ( 100.0 0.0 5.0)
   ( 100.0 5 5.0)
   ( 0.0 5 5.0)

  ( 0.0 0.0 100.0)
  ( 100.0 0.0 100.0)
  ( 100.0 5 100.0)
 ( 0.0 5 100.0)

);

blocks
(
hex (0 1 2 3 4 5 6 7) (200 10 30) simpleGrading (1 1 0.2)
hex ( 4 5 6 7 8 9 10 11) (200 10 200) simpleGrading (1 1 5)

);

topoSet命令如下
actions
(
{
name c0;
type cellSet;
action new;
source boxToCell;
box (0 0 -3) (100 5 10);
}
);

refineMesh的命令
refineMesh -overwrite的命令运行也正常

checkMesh -allTopology -allGeometry结果如下
/---------------------------------------------------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2206 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
Build : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206
Arch : "LSB;label=32;scalar=64"
Exec : /BIGDATA1/app/OpenFOAM/2206/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/checkMesh -allTopology -allGeometry
Date : Jan 11 2024
Time : 11:19:36
Host : cn7359
PID : 10713
I/O : uncollated
Case : /BIGDATA1/nieer_cwfu_1/OpenFOAM/OpenFOAM22_RUN/tutorials/multiphase/interFoam/laminar/waves/stokesV_test/test_ABL03
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
points: 2168121
faces: 6267050
internal faces: 6092950
cells: 2050000
faces per cell: 6.02927
boundary patches: 7
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 2030000
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 20000
Breakdown of polyhedra by number of faces:
faces number of cells
9 20000

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Topological cell zip-up check OK.
Face-face connectivity OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology Bounding box
inlet_water 960 1031 ok (non-closed singly connected) (0 0 -10) (0 5 5)
inlet_air 3090 3331 ok (non-closed singly connected) (0 0 5) (0 5 100)
outlet_water 960 1031 ok (non-closed singly connected) (100 0 -10) (100 5 5)
outlet_air 3090 3331 ok (non-closed singly connected) (100 0 5) (100 5 100)
ground 2000 2211 ok (non-closed singly connected) (0 0 -10) (100 5 -10)
top 2000 2211 ok (non-closed singly connected) (0 0 100) (100 5 100)
sides 162000 164062 ok (non-closed singly connected) (0 0 -10) (100 5 100)

Checking faceZone topology for multiply connected surfaces...
No faceZones found.

Checking basic cellZone addressing...
No cellZones found.

Checking geometry...
Overall domain bounding box (0 0 -10) (100 5 100)
Mesh has 2 geometric (non-empty/wedge) directions (1 0 1)
Mesh has 2 solution (non-empty) directions (1 0 1)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-1.6938e-16 1.59973e-16 0) OK.
Max cell openness = 1.08181e-16 OK.
Max aspect ratio = 4.84872 OK.
Minimum face area = 0.015625. Maximum face area = 0.49938. Face area magnitudes OK.
Min volume = 0.00298312. Max volume = 0.24969. Total volume = 55000. Cell volumes OK.
Mesh non-orthogonality Max: 29.5789 average: 1.99555
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.497993 OK.
Coupled point location match (average 0) OK.
Face tets OK.
Min/max edge length = 0.125 0.99876 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : min = 1 average = 1
All face flatness OK.
Cell determinant (wellposedness) : minimum: 0.0490101 average: 0.322733
Cell determinant check OK.
***Concave cells (using face planes) found, number of cells: 20000
<<Writing 20000 concave cells to set concaveCells
Face interpolation weight : minimum: 0.486115 average: 0.497888
Face interpolation weight check OK.
Face volume ratio : minimum: 0.236503 average: 0.982199
Face volume ratio check OK.

Failed 1 mesh checks.

End

csj1246957849

您好，我也在研究topset和refinemesh的问题，您知道在blockmesh生成一个梯形（具有斜面）的网格时，使用refinemesh会出现报错是什么原因吗？

CFD中文网

局部加密区域跨两个block，运行topoSet和refineMesh命令会报错