STAR-CCM+ CCM mesh를 이용한 CHT + Porous 해석 시 solver가 에러 메시지 없이 종료되는 문제

안녕하세요.

STAR-CCM+에서 생성한 mesh를 OpenFOAM의 CCMtoFOAM utility를 이용하여 OpenFOAM mesh로 변환한 뒤, 해당 mesh를 BARAM으로 불러와 porous 영역이 포함된 복합열전달 해석​을 수행해보고 있습니다.

1. Mesh 구성 및 변환 과정

해석하고자 하는 케이스의 STAR-CCM+의 유체 영역은 다음과 같은 구조입니다.(고체 영역은 단순 파이프 및 부속입니다)

FLUID - Porous - FLUID

STAR-CCM+에서는 FLUID와 POROUS가 각각 별도의 region으로 구분되어 있습니다.

FLUID와 POROUS region을 하나의 CCM 파일로 export한 뒤 CCMtoFOAM으로 변환하였으며, 변환 결과 OpenFOAM에서는 전체 영역이 하나의 Fluid mesh로 구성되고, 기존 POROUS 영역은 별도의 cellZone으로 정상적으로 생성된 것을 확인했습니다.

이후 변환된 polyMesh들을 BARAM에서 Multiple polyMesh 방식으로 불러와 Fluid/Solid region 및 porous 조건을 설정한 뒤 해석을 실행했습니다.

2. Solver 종료 현상

계산을 시작하면 첫 번째 iteration에서 Fluid region의 momentum equation과 pressure equation까지는 계산되는 것으로 보이지만, 그 직후 solver가 terminated 됩니다. Fluid region의 energy equation 관련 출력이 나오기 전에 종료됩니다.

stderr.log, log.txt 파일은 비어있고 stdout.log는 console과 동일 내용들이 저장되는 것 같았습니다.

[Baram Console]

Time = 1

Solving for fluid region FLUID

DILUPBiCGStab: Solving for Ux, Initial residual = 0.00283527498894, Final residual = 1.27478738423e-05, No Iterations 1

DILUPBiCGStab: Solving for Uy, Initial residual = 0.00152214992613, Final residual = 2.36942452232e-06, No Iterations 1

DILUPBiCGStab: Solving for Uz, Initial residual = 0.00320106280407, Final residual = 1.16962783583e-05, No Iterations 1

GAMGPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0625643859075, No Iterations 2

GAMGPCG: Solving for p_rgh, Initial residual = 0.0780646675445, Final residual = 0.00296566495384, No Iterations 1

3. 문의 사항

 과정자체는 문제없이 정상적인지와 기능상 아직 미구현인 부분이라거나 등 혹여나 짐작가시는 부분 있는지 여쭤보고자 문의드립니다.