CFD中文网

TaXLee

首先很感谢东岳大大对新CFDer的帮助, 您的文章与提供的资料真的对我、对新的CFDer有很大的帮助。
今天在研究东岳流体上<SIMPLE PISO算法的评价>时，感觉一句话不太通顺，如下：
“继续把方程(5)还是一个关于速度的方程，同时把方程(4)转化为限制性条件，即转化为压力泊松方程；”
这句话在公式(6)上方。
也可能是个人没有读懂，希望能够补充一下。
非常感谢！
以上

TaXLee

@李东岳好的，谢谢东岳大大的建议！我再琢磨琢磨！

TaXLee

大家好，刚开始学OpenFoam，在用porousSimpleFoam。想用blockMesh生成一个这样的几何图形：

图片中的数字是点的编号，中心的正方体是一个多孔介质区域，外围就是正常的水域。想研究流体从左侧流入后，遇到多孔介质后的流体变化。
本来是建立了一个大正方体，里面有一个小正方体，但是发现流体直接传过去了，没有发生干扰的现象，不知道该怎么修改，所以换了这个麻烦的方法。
blockMeshDict中主体代码是这样的：

scale   0.01;

vertices
(
	(0 0 0) //0//
	(0 0 100) //1//
	(36 0 0) //2//
	(36 0 36) //3//
	(36 0 64) //4//
	(36 0 100) //5//
	(64 0 0) //6//
	(64 0 36) //7//
	(64 0 64) //8//
	(64 0 100) //9//
	(100 0 0) //10//
	(100 0 100) //11//
	(36 36 36) //12//
	(36 36 64) //13//
	(64 36 36) //14//
	(64 36 64) //15//
	(36 64 36) //16//
	(36 64 64) //17//
	(64 64 36) //18//
	(64 64 64) //19//
	(0 100 0) //20//
	(0 100 100) //21//
	(36 100 0) //22//
	(36 100 36) //23//
	(36 100 64) //24//
	(36 100 100) //25//
	(64 100 0) //26//
	(64 100 36) //27//
	(64 100 64) //28//
	(64 100 100) //29//
	(100 100 0) //30//
	(100 100 100) //31//
);

blocks
(
    // wall block
    hex (0 1 5 2 20 21 25 22) 
	inlet ( 15 20 20 ) simpleGrading (1 1 1)
	
	hex (2 3 7 6 22 23 27 26)
	middle1 ( 20 20 20 ) simpleGrading (1 1 1)
	
    hex (16 17 19 18 23 24 28 27)
	middle2	( 20 20 20 ) simpleGrading (1 1 1)
	
	hex (3 4 8 7 12 13 15 14)
	middle3 ( 20 20 20 ) simpleGrading (1 1 1)
	
	hex (4 5 9 8 24 25 29 28)
	middle4 ( 20 20 20 ) simpleGrading (1 1 1)
	
	hex (6 9 11 10 26 29 31 30) 
	outlet ( 20 20 20 )  simpleGrading (1 1 1)
	
    // porosity block
    hex (12 13 15 14 16 17 19 18) 
	porosity ( 20 20 20 ) simpleGrading (1 1 1)
);

edges
(
);

patches
(

    wall front
    (
		(1 5 25 21)
		(5 9 29 25)
		(9 11 31 29)
    )

    wall back
    (
		(0 2 22 20)
		(2 6 26 22)
		(6 10 30 26)
    )

    wall walls
    (
		(0 1 5 2)
		(4 5 9 8)
		(3 4 8 7)
		(2 3 7 6)
		(6 9 11 10)
		(20 21 25 22)
		(22 23 27 26)
		(23 24 28 27)
		(24 25 29 28)
		(26 29 31 30)
    )

    wall porosityWall
    (
	(16 17 19 18)
	(16 17 13 12)
	(13 12 14 15)
	(14 15 19 18)
	(17 13 15 19)
	(16 12 14 18)
    )

    patch inlet
    (
	(0 1 21 20)
    )

    patch outlet
    (
	(10 11 31 30)
    )
);

mergePatchPairs
(
);

运行blockMesh后报错是这样的：

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1712                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1712
Arch   : "LSB;label=32;scalar=64"
Exec   : blockMesh
Date   : Feb 25 2018
Time   : 19:32:16
Host   : "DESKTOP-PB8VIPU"
PID    : 839
I/O    : uncollated
Case   : /mnt/e/tutorials/incompressible/porousSimpleFoam/Try1
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Deleting polyMesh directory
    "/mnt/e/tutorials/incompressible/porousSimpleFoam/Try1/constant/polyMesh"
Creating block mesh from
    "/mnt/e/tutorials/incompressible/porousSimpleFoam/Try1/system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Reading patches section

Creating block mesh topology

Reading physicalType from existing boundary file

Default patch type set to empty


--> FOAM FATAL ERROR:
Trying to specify a boundary face 4(16 17 19 18) on the face on cell 2 which is either an internal face or already belongs to some other patch.  This is face 0 of patch 3 named porosityWall.

    From function void Foam::polyMesh::setTopology(const cellShapeList&, const faceListList&, const wordList&, Foam::labelList&, Foam::labelList&, Foam::label&, Foam::label&, Foam::cellList&)
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 322.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&)sh: 1: addr2line: not found
 addr2line failed
#1  Foam::error::abort()sh: 1: addr2line: not found
 addr2line failed
#2  Foam::polyMesh::setTopology(Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&)sh: 1: addr2line: not found
 addr2line failed
#3  Foam::polyMesh::polyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool)sh: 1: addr2line: not found
 addr2line failed
#4  Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&)sh: 1: addr2line: not found
 addr2line failed
#5  Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&)sh: 1: addr2line: not found
 addr2line failed
#6  ?sh: 1: addr2line: not found
 addr2line failed
#7  __libc_start_mainsh: 1: addr2line: not found
 addr2line failed
#8  ?sh: 1: addr2line: not found
 addr2line failed
Aborted (core dumped)

刚开始学OpenFoam，还比较小白，问题应该是出在(16 17 19 18)，也就是多孔介质的上层，和流域的block重叠了，也就是这一块：

求助，应该怎么修改？万分感谢！

TaXLee

现在Word的公式可以直接复制转换成LaTeX代码（Word 2016或者Office 365版本），在公式页面勾选“{}LaTeX”选项，直接copy，然后paste到LaTeX编译器里就是代码了。另外"Alt"+"="是插入公式的快捷键，直接输入\Alpha_l，\Delta这样的转义字符可以直接变成相应的符号，很多公式里的代码和LaTeX一样。公式结尾输入"#"+“（数字）”，再按Enter，会直接排版编号，但是仍然不能自动编号，有些麻烦。

CFD中文网

TaXLee

帖子