T字管路用のblockMeshDictについて

61 views
Skip to first unread message

梅澤巧

unread,
Nov 11, 2025, 9:03:25 PMNov 11
to OpenFOAM
お世話になっております。
初歩的な質問で申し訳ないのですが、ぜひアドバイスをいただきたいです。

現在大学の学部生でOpenFOAMを用いて、気水二相流の研究をしております。
私が使用しているOpenFOAMのバージョンは、OpenFOAM-12です。
研究室内でOpenFOAMを使用して研究をすることが初めてです。

私は今、気水二相流の解析を行うための管路を作成して、以下のようなblockMeshDictを自作してみたのですが、合流部分にdefaultが生成されてしまい、合流部分がふさがってしまっているのが現状です。

他に必要な情報がありましたら、添付致します。

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
    location    "system";
}
convertToMeters 1.0;


vertices
(
    (0   -0.02      -0.02)   //  0
    (0    0.02      -0.02)   //  1
    (0    0.02       0.02)   //  2
    (0   -0.02       0.02)   //  3
    (0   -0.0353553391  -0.0353553391) //  4
    (0    0.0353553391  -0.0353553391) //  5
    (0    0.0353553391   0.0353553391) //  6
    (0   -0.0353553391   0.0353553391) //  7
    (2.9646446609  -0.02      -0.02)   //  8
    (2.9646446609   0.02      -0.02)   //  9
    (2.9646446609   0.02       0.02)   // 10
    (2.9646446609  -0.02       0.02)   // 11
    (2.9646446609  -0.0353553391  -0.0353553391) // 12
    (2.9646446609   0.0353553391  -0.0353553391) // 13
    (2.9646446609   0.0353553391   0.0353553391) // 14
    (2.9646446609  -0.0353553391   0.0353553391) // 15
    (3.0353553391  -0.02      -0.02)   // 16
    (3.0353553391   0.02      -0.02)   // 17
    (3.0353553391   0.02       0.02)   // 18
    (3.0353553391  -0.02       0.02)   // 19
    (3.0353553391  -0.0353553391  -0.0353553391) // 20
    (3.0353553391   0.0353553391  -0.0353553391) // 21
    (3.0353553391   0.0353553391   0.0353553391) // 22
    (3.0353553391  -0.0353553391   0.0353553391) // 23
    (11.06  -0.02      -0.02)   // 24
    (11.06   0.02      -0.02)   // 25
    (11.06   0.02       0.02)   // 26
    (11.06  -0.02       0.02)   // 27
    (11.06  -0.0353553391  -0.0353553391) // 28
    (11.06  0.0353553391  -0.0353553391) // 29
    (11.06   0.0353553391   0.0353553391) // 30
    (11.06  -0.0353553391   0.0353553391) // 31
    (2.98    0.0353553391     -0.02)   // 32
    (3.02    0.0353553391     -0.02)   // 33
    (3.02    0.0353553391      0.02)   // 34
    (2.98    0.0353553391      0.02)   // 35
    (2.98   5.0353553391   -0.02)   // 36
    (3.02   5.0353553391   -0.02)   // 37
    (3.02   5.0353553391    0.02)   // 38
    (2.98   5.0353553391    0.02)   // 39
    (2.9646446609  5.0353553391   -0.0353553391) // 40
    (3.0353553391  5.0353553391   -0.0353553391) // 41
    (3.0353553391  5.0353553391    0.0353553391) // 42
    (2.9646446609  5.0353553391    0.0353553391) // 43
    (2.98  0.02  -0.02) //44
    (3.02  0.02  -0.02) //45
    (3.02  0.02   0.02) //46
    (2.98  0.02   0.02) //47
    (2.9646446609  0.02  -0.0353553391) //48
    (3.0353553391  0.02  -0.0353553391) //49
    (3.0353553391  0.02   0.0353553391) //50
    (2.9646446609  0.02   0.0353553391) //51
);

blocks
(
    hex (0 1 2 3   8 9 10 11)   (8 8 124) simpleGrading (1 1 1)
    hex (3 2 6 7  11 10 14 15) (8 6 124) simpleGrading (1 1 1)
    hex (1 5 6 2   9 13 14 10)  (6 8 124) simpleGrading (1 1 1)
    hex (4 5 1 0  12 13 9 8)    (8 6 124) simpleGrading (1 1 1)
    hex (4 0 3 7  12 8 11 15)   (6 8 124) simpleGrading (1 1 1)
    hex (8 9 10 11   16 17 18 19)   (8 8 3) simpleGrading (1 1 1)
    hex (11 10 14 15  19 18 22 23) (8 6 3) simpleGrading (1 1 1)
    hex (12 13 9 8   20 21 17 16)  (8 6 3) simpleGrading (1 1 1)
    hex (12 8 11 15  20 16 19 23)  (6 8 3) simpleGrading (1 1 1)
    hex (16 17 18 19  24 25 26 27) (8 8 335) simpleGrading (1 1 1)
    hex (19 18 22 23  27 26 30 31) (8 6 335) simpleGrading (1 1 1)
    hex (17 21 22 18  25 29 30 26) (6 8 335) simpleGrading (1 1 1)
    hex (20 21 17 16  28 29 25 24) (8 6 335) simpleGrading (1 1 1)
    hex (20 16 19 23  28 24 27 31) (6 8 335) simpleGrading (1 1 1)
    hex (36 37 38 39   32 33 34 35) (8 8 209) simpleGrading (1 1 1)
    hex (39 38 42 43   35 34 22 14) (8 6 209) simpleGrading (1 1 1)
    hex (37 41 42 38   33 21 22 34) (6 8 209) simpleGrading (1 1 1)
    hex (40 41 37 36   13 21 33 32) (8 6 209) simpleGrading (1 1 1)
    hex (40 36 39 43   13 32 35 14) (6 8 209) simpleGrading (1 1 1)
    hex (32 33 34 35  44 45 46 47) (8 8 1) simpleGrading (1 1 1)
    hex (35 34 22 14  47 46 50 51) (8 6 1) simpleGrading (1 1 1)
    hex (33 21 22 34  45 49 50 46) (6 8 1) simpleGrading (1 1 1)
    hex (13 21 33 32  48 49 45 44) (8 6 1) simpleGrading (1 1 1)
    hex (13 32 35 14  48 44 47 51) (6 8 1) simpleGrading (1 1 1)
);

edges
(
    arc 6 7  (0    0.0   0.05)
    arc 5 6  (0    0.05   0.0)
    arc 4 5  (0    0.0  -0.05)
    arc 7 4  (0   -0.05   0.0)
    arc 14 15 (2.9646446609   0.0   0.05)
    arc 12 13 (2.9646446609   0.0  -0.05)
    arc 15 12 (2.9646446609  -0.05   0.0)
    arc 22 23 (3.0353553391   0.0   0.05)
    arc 20 21 (3.0353553391   0.0  -0.05)
    arc 23 20 (3.0353553391  -0.05   0.0)
    arc 30 31 (11.06   0.0   0.05)
    arc 29 30 (11.06   0.05   0.0)
    arc 28 29 (11.06   0.0  -0.05)
    arc 31 28 (11.06  -0.05   0.0)
    arc 22 14 (3     0.0353553391    0.05)
    arc 21 22 (3.05  0.0353553391    0   )
    arc 13 21 (3     0.0353553391   -0.05)
    arc 14 13 (2.95  0.0353553391    0   )
    arc 42 43 (3     5.0353553391   0.05)
    arc 41 42 (3.05  5.0353553391    0.0)
    arc 40 41 (3     5.0353553391  -0.05)
    arc 43 40 (2.95  5.0353553391    0.0)
    arc 50 51 (3     0.02    0.05)
    arc 49 50 (3.05  0.02    0   )
    arc 48 49 (3     0.02   -0.05)
    arc 51 48 (2.95  0.02    0   )
);

boundary
(
    inlet1
    {
        type patch;
        faces
        (
            (0 1 2 3)
            (3 2 6 7)
            (1 5 6 2)
            (4 5 1 0)
            (4 0 3 7)
        );
    }

    outlet
    {
        type patch;
        faces
        (
            (24 25 26 27)
            (27 26 30 31)
            (25 29 30 26)
            (28 29 25 24)
            (28 24 27 31)
        );
    }

    inlet2
    {
        type patch;
        faces
        (
            (36 37 38 39)
            (39 38 42 43)
            (37 41 42 38)
            (40 41 37 36)
            (40 36 39 43)
        );
    }
    wall
    {
        type wall;
        faces
        (
            (7 6 14 15)
            (6 5 13 14)
            (4 5 13 12)
            (7 4 12 15)
            (15 14 22 23)
            (12 13 21 20)
            (23 22 30 31)
            (22 21 29 30)
            (20 21 29 28)
            (23 20 28 31)
            (43 42 22 14)
            (42 41 21 22)
            (40 41 21 13)
            (43 40 13 14)
            (14 22 50 51)
            (13 21 49 48)
        );
    }
);

defaultPatch
{
    name default;
    type patch;
}

mergePatchPairs();


E.Mogura

unread,
Nov 11, 2025, 10:11:05 PMNov 11
to open...@googlegroups.com
E.Mogura です

OpenFOAM-12を使っているのであれば、

paraFoam -blockMesh

にて、ParaViewでblockMeshを図化表示できます。

そうすると、たとえば下図のように、怪しい箇所を比較的簡単に見つけることができるので、お試しください。

image.png

2025年11月12日(水) 11:03 梅澤巧 <ume...@wem-uu.org>:
--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
このディスカッションを表示するには、https://groups.google.com/d/msgid/openfoam/a680419e-d180-41fd-bc57-991898342e28n%40googlegroups.com にアクセスしてください。

梅澤巧

unread,
Nov 11, 2025, 10:45:43 PMNov 11
to OpenFOAM
E.Mogura 様 

  ご返信ありがとうございます。

図化表示の仕方を教えていただきありがとうございます。
自分の管路がどのようになっているのか理解しやすくなりました。

T字管路の現状として、x方向に延びる管路を主管、-y方向に延びる管路を支管としたとき、
主管に支管が同じpointを使っているからくっついているように見えるけど、メッシュが連続していないから、実際には合流部分が繋がっていないという理解で正しいでしょうか?

教えていただけると幸いです。

2025年11月12日水曜日 12:11:05 UTC+9 E.Mogura:

E.Mogura

unread,
Nov 12, 2025, 1:09:43 AMNov 12
to open...@googlegroups.com
E.Mogura です

言葉で説明するのは何とも難しいのですが、
隣接するブロックで、面を共有するように構成する必要があるということですかね。

赤のブロックと水色のブロックの接合面の形状が違うので、つながりません。
水色のブロックで、「51」でなく、「10」を使うとかすれば良いです。
(下の方も同様に変更が必要ですが、ちょっとこの図からは詳細不明です)
あとhex定義ブロックでは、8つの節点の並び順にも注意する必要があります。



2025年11月12日(水) 12:45 梅澤巧 <ume...@wem-uu.org>:
このディスカッションを表示するには、https://groups.google.com/d/msgid/openfoam/4f2ddba2-5fc3-42c5-8c6c-7364a1d80bafn%40googlegroups.com にアクセスしてください。

梅澤巧

unread,
Nov 12, 2025, 2:53:20 AMNov 12
to OpenFOAM
E.Mogura 様

お返事いただきありがとうございます。

いただいたアドバイスを参考に、もう一度考えてみます。

2025年11月12日水曜日 15:09:43 UTC+9 E.Mogura:

梅澤巧

unread,
Nov 13, 2025, 2:08:34 AMNov 13
to OpenFOAM

いただいたアドバイスを元に、考えてみたのですが、自分の中で理解できませんでした。

T字管路を作成する上で、おすすめのやり方や参考になる資料を教えていただきたいです。

何卒よろしくお願いいたします。


2025年11月12日水曜日 16:53:20 UTC+9 梅澤巧:

ttsy shmz

unread,
Nov 13, 2025, 3:45:09 AMNov 13
to open...@googlegroups.com, OpenFOAM

梅澤巧

unread,
Nov 14, 2025, 1:36:41 AMNov 14
to OpenFOAM

ttsy shmz 様

お返事いただきありがとうございます。

試してみます。



2025年11月13日木曜日 17:45:09 UTC+9 ttsy shmz:
Reply all
Reply to author
Forward
0 new messages