1) Porous Media 개요

* 형상 파일 링크


형상 및 유동장

본 에제는 porous media 조건을 사용한 유동해석 예제이다. 상부 왼쪽에서 유동이 유입되고 porous 영역을 지나 아래로 유동이 흐르는 문제이다.(위 그림 왼쪽에서 파란색 부분이 porous 영역)

OpenFOAM의 porous media 모델은 porous 영역에서 불연속적인 속도 분포가 나타나고 압력손실도 입력 조건과 조금 다른 결과를 보이는 문제가 있다. 아래 그림의 왼쪽이 Baram-v23의 결과이며, 오른쪽이 OpenFOAM 2306의 standard 솔버를 사용한 결과이다.

Baram이 사용하는 NextFOAM에서는 porous 영역에서 압력의 interpolation 방법을 개선하여 이 문제를 해결하였다(이에 대한 자세한 내용은 아래 링크의 문서를 참고). 결과의 정확성과 함께 수렴성도 많이 좋아진 것을 확인할 수 있다.

*Porous Media 참고 문헌


결과 (좌)Baram v23, (우) openfoam 2306 standard solver


Residual (좌)Baram v23, (우) openfoam 2306 standard solver


2) 형상(Geometry)

덕트의 형상은 stl 파일을 사용한다.

Import 버튼을 눌러 duct.stl 파일을 선택하고 ‘Split Surface’ 옵션을 선택하여 경계면을 구분한다. 하나의 면으로 되어있는 형상 파일이 feature angle에 따라 20개의 면으로 나누어진다.


덕트 파일 불러오기

20개의 면들 중에 입구와 출구를 구분하기 쉽게 이름을 바꿔준다. 그래픽창에서 마우스로 입구를 클릭하면 Geometry 리스트에서 해당 면이 활성화된다. 마우스 오른쪽 버튼으로 Edit를 누르고 이름을 inlet, outlet으로 바꾸어 준다.

Porous 영역은 Add 버튼을 눌러 Hex를 이용해서 만든다. Hex의 설정은 다음과 같다.

  • Name : porousZone

  • Type : CellZone

  • Min. : (-0.3 -0.2 0.8)

  • Max. : (0.3 0.2 0.9)


Porous 영역 생성

Porous 영역을 생성하면 그 하위에 Hex_1_surface라는 것이 생성되는데 이것의 Type을 None으로 설정한다.

최종적인 Geometry 리스트는 다음 그림과 같이 된다.


Geomerty 설정

Next 버튼을 눌러 다음 단계로 넘어간다.

3) Region

상단 (+) 버튼을 눌러 영역을 생성하고 그래픽 창에 연두색으로 나타나는 선의 교차점을 마우스로 이동하여 유체 영역에 위치시킨 후 Add 버튼을 클릭하면 설정이 완료된다.


Region 설정

Next 버튼을 눌러 다음 단계로 넘어간다.

4) Base Grid

격자수를 60, 40, 120으로 설정한다. Generate 버튼을 누르면 배경 격자가 생성된다.


Base Grid 설정

Next 버튼을 눌러 다음 단계로 넘어간다.

5) Castellation

전체 영역의 격자를 균일하게 만들 것이기 때문에 별도의 refinement 설정 없이 Refine 버튼을 누르면 castellation이 진행된다.

작업이 끝나면 Next 버튼을 눌러 다음 단계로 넘어간다.

6) Snap

설정은 다음과 같이 변경한다.

  • Smoothing for Surface : 1

  • Smoothing for Internal : 3

  • Mesh Displacement Relaxation : 30

  • Snapping Relaxation : 15

  • Tolerance : 1

  • Concave Angle (degree) : 45

  • Min.Area Ratio : 0.3


snap 설정

7) Boundary Layer

덕트 벽면에 경계층 격자를 생성한다.

Configuration에서 (+)를 눌러 다음과 같이 설정한다.


Boundary Layer Configuration

  • Number of Layers : 3

  • Thickness Model Specification : Final and Expansion

  • Size Specification : Relative

  • First Layer Thickness : 0.3

  • Expansion Ratio : 1.2

  • Min. Total Thickness : 0.1

  • Boundary : 모든 덕트 벽면 선택

나머지는 Dafault 설정 그대로 적용하고 apply 버튼을 누른다.

작업이 끝나면 Next 버튼을 눌러 다음 단계로 넘어간다.

8) Export

마지막으로 porous라는 이름으로 Export 하면 baramFlow v23에서 열 수 있는 Project 폴더가 생성된다.


inlet 근처 격자


outlet 근처 격자


전체 격자