浮力流的初始化问题

tyhero123

有人知道浮力流里面的prgh该怎么设置边界条件和初始化吗，我是一个这样的模型，用的是of10的buoyantFoam：

高度有两米，右边两个口底下那个是进口，上面那个是出口。
流体的密度有10000。
现在就是刚开始算流体就会一直被加速，感觉就是跟压力的设置有关系，初场或者边界，很多方法试了都没法解决，想问问有没有大佬懂这一块的。

tyhero123

入口是固定的速度

李东岳

fixedFluxPressure不对么

tyhero123

@李东岳在浮力流的初始化问题中说：

fixedFluxPressure不对么

我给壁面和入口设置的fixedFluxPressure，或者只给壁面设置，都不行，会发散掉。

李东岳

算例网格少的发上来我给你看下

tyhero123

test.zip
其他一些物理场都暂时无关，只要buoyantFoam能跑通就行，辛苦老师。

李东岳


blocks
(
    //- block0
    hex (0 1 2 3 4 5 6 7) (33 33 56) 
    simpleGrading 
    (
		(
			(0.39716 0.45455 3.7975) (0.20569 0.09091 1) (0.39716 0.45455 0.26333)
		)
		(
			(0.39716 0.45455 3.7975) (0.20569 0.09091 1) (0.39716 0.45455 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block1
    hex (3 2 9 8 7 6 13 12) (33 365 56) 
    simpleGrading 
    (
		(
			(0.39716 0.45455 3.7975) (0.20569 0.09091 1) (0.39716 0.45455 0.26333)
		)
		(
			(0.01826 0.0411 3.7975) (0.96348 0.91781 1) (0.01826 0.0411 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block2
    hex (8 9 10 11 12 13 14 15) (33 41 56)
    simpleGrading 
    (
		(
			(0.39716 0.45455 3.7975) (0.20569 0.09091 1) (0.39716 0.45455 0.26333)
		)
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block3
    hex (1 36 37 2 5 39 38 6) (129 33 56)
    simpleGrading
    (
		(
			(0.05252 0.11628 3.7975) (0.94748 0.88372 1) (0 0 1)
		)
		(
			(0.39716 0.45455 3.7975) (0.20569 0.09091 1) (0.39716 0.45455 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block4
    hex (9 16 19 10 13 20 23 14) (15 41 56)
    simpleGrading 
    (
		1
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block5
    hex (16 17 18 19 20 21 22 23) (41 41 56)
    simpleGrading 
    (
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block6
    hex (27 26 17 16 31 30 21 20) (41 338 56)
    simpleGrading 
    (
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.0198 0.04438 3.7975) (0.96041 0.91124 1) (0.0198 0.04438 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block7
    hex (24 25 26 27 28 29 30 31) (41 41 56)
    simpleGrading 
    (
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
    //- block8
    hex (25 32 33 26 29 35 34 30) (102 41 56)
    simpleGrading 
    (
		(
			(0.0676 0.14706 3.7975) (0.9324 0.85294 1) (0 0 1)
		)
		(
			(0.26477 0.36585 3.7975) (0.47046 0.26829 1) (0.26477 0.36585 0.26333)
		)
		(
			(0.17054 0.26786 3.7975) (0.65891 0.46429 1) (0.17054 0.26786 0.26333)
		)
    )
);

你把上面的blockMesh的block网格减少，处理成1万个网格一下，然后用跟我一样的形式发上来

tyhero123

blocks
(
    //- block0
    hex (0 1 2 3 4 5 6 7) (10 10 10) simpleGrading (1 1 1)
    //- block1
    hex (3 2 9 8 7 6 13 12) (10 20 10) simpleGrading (1 1 1)
    //- block2
    hex (8 9 10 11 12 13 14 15) (10 10 10) simpleGrading (1 1 1)
    //- block3
    hex (1 36 37 2 5 39 38 6) (10 10 10) simpleGrading (1 1 1)
    //- block4
    hex (9 16 19 10 13 20 23 14) (10 10 10) simpleGrading (1 1 1)
    //- block5
    hex (16 17 18 19 20 21 22 23) (10 10 10) simpleGrading (1 1 1)
    //- block6
    hex (27 26 17 16 31 30 21 20) (10 20 10) simpleGrading (1 1 1)
    //- block7
    hex (24 25 26 27 28 29 30 31) (10 10 10) simpleGrading (1 1 1)
    //- block8
    hex (25 32 33 26 29 35 34 30) (10 10 10) simpleGrading (1 1 1)
);

这样网格不会太粗吗

李东岳

test22.tar.xz

直接blockMesh，然后buoyantFoam跑就可以，我跑到了20多秒没啥问题，流场看起来也正确。

tyhero123

@李东岳在浮力流的初始化问题中说：

test22.tar.xz

直接blockMesh，然后buoyantFoam跑就可以，我跑到了20多秒没啥问题，流场看起来也正确。

谢谢，我试试看，我一直是发散的。

CFD中文网

浮力流的初始化问题