CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

No interfaces between the solid and fluid after splitMeshRegions

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 17, 2021, 07:39
Default No interfaces between the solid and fluid after splitMeshRegions
  #1
New Member
 
a
Join Date: Apr 2019
Posts: 2
Rep Power: 0
a2019 is on a distinguished road
Hi All,

I am very new to OpenFoam and trying to run the case with chtMultiRegionFoam. The case is very simpele geometrically (fluid layer "r4" between the solid layers). The mesh is created by blockMesh and then command splitMeshRegions -cellZones -overwrite is used. The problem is that the interface is not created between the fluid and solid regions. How to deal with this problem?

blockMeshDisct:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

r1_depth0 10;
r2_depth0 180;
r3_depth0 15;
r4_depth0 28;
r5_depth0 100;

r4_width0 45;

case_width #calc " 2 * $r4_width0 ";
case_height 21; //75;

r2_depth #calc " $r2_depth0 + $r1_depth0 ";
r3_depth #calc " $r3_depth0 + $r2_depth ";
r4_depth #calc " $r4_depth0 + $r3_depth ";
r5_depth #calc " $r5_depth0 + $r4_depth ";

bcs 3;
bcf 1;
ratio 1.5;
nnw #calc " ceil( $case_width / $bcs ) ";
nnh #calc " ceil( $case_height / $bcs / $ratio) ";
nn_r1 #calc " ceil( $r1_depth0 / $bcs ) ";
nn_r2 #calc " ceil( $r2_depth0 / $bcs ) ";
nn_r3 #calc " ceil( $r3_depth0 / $bcs ) ";
nn_r5 #calc " ceil( $r5_depth0 / $bcs ) ";
nn_r4 #calc " ceil( $r4_depth0 / $bcf ) ";

convertToMeters    0.001;

vertices
(
// r1
  (0            0           0)  // vert 0
  ($r1_depth0 0           0)  // vert 1
  ($r1_depth0 $case_width 0)  // vert 2
  (0            $case_width 0)  // vert 3
  (0            0           $case_height)  // vert 4
  ($r1_depth0 0           $case_height)  // vert 5
  ($r1_depth0 $case_width $case_height)  // vert 6
  (0            $case_width $case_height)  // vert 7

// r2
  ($r1_depth0 0           0)  // vert 8
  ($r2_depth  0           0)  // vert 9
  ($r2_depth  $case_width 0)  // vert 10
  ($r1_depth0 $case_width 0)  // vert 11
  ($r1_depth0 0           $case_height)  // vert 12
  ($r2_depth  0           $case_height)  // vert 13
  ($r2_depth  $case_width $case_height)  // vert 14
  ($r1_depth0 $case_width $case_height)  // vert 15

// r3
  ($r2_depth 0           0)  // vert 16
  ($r3_depth 0           0)  // vert 17
  ($r3_depth $case_width 0)  // vert 18
  ($r2_depth $case_width 0)  // vert 19
  ($r2_depth 0           $case_height)  // vert 20
  ($r3_depth 0           $case_height)  // vert 21
  ($r3_depth $case_width $case_height)  // vert 22
  ($r2_depth $case_width $case_height)  // vert 23

// r4
  ($r3_depth   0           0)  // vert 24
  ($r4_depth 0           0)  // vert 25
  ($r4_depth $case_width 0)  // vert 26
  ($r3_depth   $case_width 0)  // vert 27
  ($r3_depth   0           $case_height)  // vert 28
  ($r4_depth 0           $case_height)  // vert 29
  ($r4_depth $case_width $case_height)  // vert 30
  ($r3_depth   $case_width $case_height)  // vert 31

// r5
  ($r4_depth   0           0)  // vert 32
  ($r5_depth      0           0)  // vert 33
  ($r5_depth      $case_width 0)  // vert 34
  ($r4_depth   $case_width 0)  // vert 35
  ($r4_depth   0           $case_height)  // vert 36
  ($r5_depth      0           $case_height)  // vert 37
  ($r5_depth      $case_width $case_height)  // vert 38
  ($r4_depth   $case_width $case_height)  // vert 39

// plaster


);

blocks
(

    // Fluid regions
    hex (24 25 26 27 28 29 30 31) r4 ($nn_r4 $nnw $nnh) simpleGrading (1 1 1)

    // Solid regions
    hex ( 0  1  2  3  4  5  6  7) r1 ($nn_r1 $nnw $nnh) simpleGrading (1 1 1)
    hex ( 8  9 10 11 12 13 14 15) r2 ($nn_r2 $nnw $nnh) simpleGrading (1 1 1)
    hex (16 17 18 19 20 21 22 23) r3 ($nn_r3 $nnw $nnh) simpleGrading (1 1 1)
    hex (32 33 34 35 36 37 38 39) r5  ($nn_r5  $nnw $nnh) simpleGrading (1 1 1)
);

boundary
(
    inlet
    {
        type patch;
        faces
        (
          (24 25 27 26)
        );
    }

    external
    {
        type patch;
        faces
        (
            (0 3 7 4)

            (33 34 38 37)
        );
    }

    walls
    {
        type wall;
//        inGroups (externalWall);
        faces
        (
           (17 18 22 21)
        );
    }

    outlet
    {
        type patch;
        faces
        (
          (28 29 30 31)
        );
    }
);



mergePatchPairs
(
);


// ************************************************************************* //
splitMeshRegions outprint:
Code:
Region    Zone    Name
------    ----    ----
0    0    r1
1    1    r2
2    2    r3
3    3    r4
4    4    r5

Sizes of interfaces between regions:

Interface    Region    Region    Faces
---------    ------    ------    -----

Reading geometric fields

Reading volScalarField cellToRegion
a2019 is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[solids4Foam] How to calculate drag coeff when using solids4Foam amuzeshi OpenFOAM CC Toolkits for Fluid-Structure Interaction 15 November 7, 2019 12:50
Question about adaptive timestepping Guille1811 CFX 25 November 12, 2017 17:38
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 05:15
Problem in setting Boundary Condition Madhatter92 CFX 12 January 12, 2016 04:39
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20


All times are GMT -4. The time now is 03:01.