Figure 2. Scour cavity obtained in the laboratory test (Ramos, 2012).

이 기술 요약 자료는 Pedro Ramos, João Pedro Pêgo & Rodrigo Maia가 2014년 3rd IAHR Europe Congress에 발표한 논문 “NUMERICAL SIMULATION OF THE FLOW AROUND A PIER USING OPENFOAM”을 기반으로 하여 분석 및 요약되었습니다.

Keywords

  • Primary Keyword: 교각 세굴 시뮬레이션
  • Secondary Keywords: OpenFOAM, CFD, Large Eddy Simulation (LES), 와류(Vortex), 전단 응력(Shear Stress), 항력 계수(Drag Coefficient), Fluent

Executive Summary

  • The Challenge: 교량 기초 주변에서 발생하는 세굴(Scour) 현상은 교량 붕괴로 이어질 수 있는 심각한 문제이며, 이를 유발하고 유지하는 구동 메커니즘에 대한 이해는 여전히 부족합니다.
  • The Method: 오픈소스 CFD 툴박스인 OpenFOAM을 사용하여 고립된 원형 교각 주변의 유동을 수치 시뮬레이션했으며, 평탄한 하상(초기 상태)과 평형 세굴 상태의 두 가지 고정층 구성에 대해 LES(Large Eddy Simulation) 난류 모델을 적용했습니다.
  • The Key Breakthrough: 이전 연구(Fluent 사용) 대비 정규화된 구조 격자(structured mesh)를 적용한 OpenFOAM 시뮬레이션은 교각 후류에서 발생하는 와류(wake vortex)를 더 크고 현실에 가깝게 모사했으며, 이는 세굴 과정 해석의 정확도를 크게 향상시켰습니다.
  • The Bottom Line: 오픈소스 CFD 소프트웨어인 OpenFOAM이 적절한 격자 구성과 설정을 통해 상용 소프트웨어에 필적하는 고품질의 복잡한 수리 현상 시뮬레이션 결과를 제공할 수 있음을 입증했습니다.

The Challenge: Why This Research Matters for CFD Professionals

교각 세굴은 하상에 위치한 교량 기초에 심각한 손상을 입히는 주요 원인으로, 전체 교량 구조의 안정성 붕괴 위험을 증가시킵니다. 포르투갈에서는 도루강의 Entre-os-Rios 교량 붕괴 사고의 부분적인 원인이 교각 세굴이었습니다. 이처럼 교각 주변의 국부 세굴은 광범위하게 연구되어 왔지만, 대부분의 연구는 실험에 기반하여 최대 평형 세굴 깊이를 예측하는 데 중점을 두었습니다.

하지만 세굴 과정을 촉발하고 유지하는 핵심적인 구동 메커니즘에 대한 지식은 여전히 부족한 실정입니다. 기존의 수치 시뮬레이션 연구들은 존재했지만, 난류 모델의 한계로 인해 와류와 같은 미세한 유동 구조의 직접적인 영향을 연구하는 데 어려움이 있었습니다. 본 연구는 수치 시뮬레이션을 통해 이 분야의 지식을 확장하고, 특히 오픈소스 CFD 툴의 성능을 검증하는 것을 목표로 합니다.

The Approach: Unpacking the Methodology

본 연구는 오픈소스 소프트웨어인 OpenFOAM을 사용하여 교각 주변 유동의 3차원 수치 시뮬레이션을 수행했습니다. 세굴 과정에서 와류와 같은 유동 구조의 역할을 규명하기 위해 LES(Large Eddy Simulation) 난류 모델이 적용되었습니다.

시뮬레이션은 두 가지 주요 구성(Configuration)으로 진행되었습니다. 1. Configuration A (평탄한 하상): 세굴 과정이 시작되기 전의 초기 상태를 모사합니다. 2. Configuration B (세굴된 하상): 실험을 통해 얻은 평형 세굴 깊이에 도달한 후의 하상 형상을 나타냅니다.

이전 연구(Ramos, 2012)와의 가장 큰 차이점은 격자(mesh) 구성에 있습니다. 본 연구에서는 교각 주변과 채널 벽, 하상 근처에 더 미세하고 정규화된 분포를 갖는 구조 격자(structured mesh)를 사용했습니다(Figure 5, 6). 이는 데이터 처리 속도를 높이고 계산의 정확성을 향상시키는 핵심적인 개선 사항입니다. 유입 경계 조건으로는 평균 유속 0.32 m/s의 대수 법칙 프로파일을 적용했으며, 난류성 비압축성 유동 해석을 위해 interFoam 솔버를 사용했습니다.

The Breakthrough: Key Findings & Data

본 연구는 OpenFOAM을 활용한 시뮬레이션이 이전 Fluent 기반 연구 결과와 비교하여 몇 가지 중요한 개선점을 보여주었습니다.

Finding 1: 향상된 후류 와류(Wake Vortex) 모사

세굴 과정에서 교각 후류에서 발생하는 와류는 하상 입자를 들어 올려 하류로 이동시키는 중요한 역할을 합니다. Configuration A(평탄한 하상)에서 OpenFOAM 시뮬레이션 결과(Figure 10)는 Fluent 결과(Figure 9)보다 후류 와류 영역이 더 크고 난류성이 강하게 나타났습니다. 논문은 이 결과가 “실험실에서 관찰된 것에 더 가깝다”고 언급하며, 개선된 격자 구성이 더 현실적인 유동 구조를 모사하는 데 기여했음을 시사합니다. 이는 와류에 의한 하상 입자의 흡입 효과가 더 강하게 나타남을 의미하며, 세굴 현상 분석의 신뢰도를 높입니다.

좌: Figure 9 (Fluent 결과), 우: Figure 10 (OpenFOAM 결과). OpenFOAM 결과에서 더 넓고 강한 후류 와류가 관찰된다.

Finding 2: 항력 계수(Drag Coefficient)의 동적 거동 분석

교각에 작용하는 항력 계수는 구조물의 동적 안정성 분석에 필수적인 데이터입니다. OpenFOAM 시뮬레이션(Configuration A)을 통해 계산된 평균 항력 계수는 1.21이었습니다(Figure 15). 이는 이전 Fluent 시뮬레이션의 0.68이나 경험식에 따른 예상치인 약 0.73과는 차이가 있습니다.

하지만 더 중요한 것은 항력 계수가 시간에 따라 진동하는 패턴을 보인다는 점입니다. 이 진동은 교각 후류에서 발생하는 와류 방출(vortex shedding)로 인해 발생하는 힘의 변화를 나타내며, 이러한 진동의 주파수 정보는 교량 구조물의 동적 해석에 매우 중요한 정보를 제공합니다. OpenFOAM 결과는 이러한 동적 거동을 성공적으로 포착했습니다.

Figure 15. 시간에 따른 항력 계수의 변화. 와류 방출로 인한 주기적인 진동이 명확하게 나타난다.

Practical Implications for R&D and Operations

  • For Hydraulic Engineers: 이 연구는 OpenFOAM과 같은 오픈소스 툴이 적절한 격자 전략과 함께 사용될 때, 하상 전단 응력(bed shear stress)을 정확하게 예측할 수 있음을 보여줍니다. 시뮬레이션 결과(Figure 14)는 세굴이 교각 측면에서 시작될 가능성이 높다는 것을 나타내며, 이는 세굴 발생의 임계 지점을 파악하는 데 도움을 줍니다.
  • For Bridge Maintenance Engineers: 평형 세굴 상태(Configuration B)에 대한 시뮬레이션은 세굴이 최대로 진행되었을 때의 하상 형상과 유동 특성에 대한 통찰력을 제공합니다. 이는 교량의 정기 점검 시 중점적으로 확인해야 할 부분을 결정하고, 보수 및 보강 전략을 수립하는 데 활용될 수 있습니다.
  • For Bridge Design Engineers: 항력 계수의 시간에 따른 진동 데이터(Figure 15)는 와류 방출로 인한 유체-구조 상호작용(Fluid-Structure Interaction)을 고려한 교량의 동적 해석에 필수적입니다. 이 연구 결과는 설계 단계에서부터 교량의 안정성을 더욱 정밀하게 평가할 수 있는 기초 자료를 제공합니다.

Paper Details

NUMERICAL SIMULATION OF THE FLOW AROUND A PIER USING OPENFOAM

1. Overview:

  • Title: NUMERICAL SIMULATION OF THE FLOW AROUND A PIER USING OPENFOAM
  • Author: PEDRO RAMOS, JOÃO PEDRO PÊGO & RODRIGO MAIA
  • Year of publication: 2014
  • Journal/academic society of publication: 3rd IAHR Europe Congress, Book of Proceedings
  • Keywords: Scour, OpenFOAM, Large Eddy Simulation, CFD.

2. Abstract:

세굴은 하상에 위치한 교량 기초 손상의 주요 원인으로, 전체 교량 구조의 안정성 붕괴 위험을 증가시킵니다. 포르투갈에서는 도루강의 Entre-os-Rios 교량 붕괴 사고가 부분적으로 교각 세굴 때문이었습니다. 교각 주변의 국부 세굴은 광범위하게 연구되었으며, 대부분 실험 연구를 기반으로 최대 평형 깊이를 예측하는 데 중점을 둡니다. 그럼에도 불구하고, 세굴 과정을 촉발하고 유지하는 구동 메커니즘에 대한 지식은 여전히 부족합니다. 본 논문은 수치 시뮬레이션을 통해 이 분야의 지식을 추가하고자 합니다. 이 연구의 주요 목표는 OpenFOAM 툴박스를 사용하여 고립된 원형 교각 주변의 유동을 두 가지 고정층 구성, 즉 세굴 과정의 시작에 해당하는 평탄한 고정층과 실험적으로 얻은 평형 세굴 깊이에 도달한 후의 하상 형상을 대표하는 구성에 대해 수치 시뮬레이션하는 것입니다. Large Eddy Simulation 난류 모델을 적용하여 두 구성에 대한 속도 및 와도장, 하상 전단 응력, 교각의 항력 계수를 특성화할 것입니다. OpenFOAM 결과는 Ramos(2012)가 Fluent를 사용하여 얻은 데이터와 비교 검증될 것입니다. Fluent와 OpenFOAM은 산업 및 유체 역학 연구에서 널리 사용되는 두 가지 CFD 도구로, 전자는 상용 소프트웨어이고 후자는 무료 오픈 소스 소프트웨어입니다. 사용된 두 3차원 수치 모델은 세굴 과정에서 중요한 역할을 하는 말굽 와류와 같은 유동 구조를 연구할 수 있게 합니다. 두 CFD 도구의 성능을 분석하고 비교할 것입니다. 수치 결과는 참고 문헌으로 인정된 서지에서 얻은 데이터와 비교 검증될 것입니다.

3. Introduction:

교량 기초 주변의 국부 세굴은 교각의 부분적 파손이나 붕괴로 이어질 수 있습니다. 수직 원형 교각이 일정한 흐름 속에서 하상에 놓이면, 하상과 상호 작용하는 유동 패턴에 변화가 생깁니다. 특징적인 유동 구조는 말굽 와류(교각 앞에서 형성됨)와 교각의 바람이 불어가는 쪽(lee-side)에서 형성되는 후류 와류 유동 패턴(보통 와류 방출 형태)입니다(Figure 1). 또한, 교각 상류에서 유동 감속의 결과로 하향류(downflow)가 존재합니다. 이러한 변화의 전반적인 효과는 일반적으로 퇴적물 이동을 증가시켜 교각 주변의 국부 세굴을 초래합니다.

4. Summary of the study:

Background of the research topic:

교각 세굴은 교량의 안정성을 위협하는 중요한 수리 현상입니다. 유동 속에 놓인 교각은 주변 유동장을 변화시키고, 이로 인해 발생하는 말굽 와류(horseshoe vortex), 하향류(downflow), 후류 와류(trailing vortex) 등이 하상의 퇴적물 이동을 촉진하여 국부적인 침식을 유발합니다.

Figure 1. Flow around pier (reproduced from Breusers and Raudkivi, 1991).
Figure 1. Flow around pier (reproduced from Breusers and Raudkivi, 1991).

Status of previous research:

세굴 현상에 대한 연구는 주로 실험적으로 이루어졌습니다. Twaites(1960), Bachelor(1967) 등의 공기역학 연구를 시작으로 Melville과 Raudkivi(1977) 등이 수중 세굴 현상으로 연구를 확장했습니다. 수치 시뮬레이션 분야에서는 Deng과 Piquet(1992)가 RANS 모델을 사용하여 3차원 시뮬레이션을 수행했으며, Olsen과 Melaaen(1993)이 처음으로 유동 시뮬레이션과 세굴 모델링을 결합했습니다. Roulund 등(2005)은 SST k-ω 모델을 사용했으나, 이는 와류의 직접적인 영향을 연구하기에는 한계가 있었습니다. Ramos(2012)는 이러한 한계를 극복하기 위해 LES 모델을 사용하여 Fluent로 시뮬레이션을 수행했으며, 본 연구는 이 결과를 비교 대상으로 삼습니다.

Figure 2. Scour cavity obtained in the laboratory test (Ramos, 2012).

Purpose of the study:

본 연구의 주된 목적은 오픈소스 CFD 툴인 OpenFOAM을 사용하여 교각 주변 유동을 수치 시뮬레이션하고, 그 결과를 상용 소프트웨어인 Fluent를 사용한 이전 연구(Ramos, 2012) 및 실험 데이터와 비교 검증하는 것입니다. 이를 통해 오픈소스 CFD 툴의 성능을 평가하고, 개선된 격자 기법을 적용하여 세굴을 유발하는 와류 구조에 대한 더 깊은 이해를 얻고자 합니다.

Core study:

연구의 핵심은 두 가지 하상 조건(평탄한 하상, 세굴된 하상)에 대해 OpenFOAM을 이용한 LES 시뮬레이션을 수행하는 것입니다. 시뮬레이션을 통해 속도장, 와도장, 하상 전단 응력, 교각의 항력 계수를 분석합니다. 특히, 이전 연구 대비 개선된 구조 격자를 사용하여 결과의 질을 향상시키고, 이를 Fluent 결과와 비교하여 OpenFOAM의 신뢰성과 유용성을 검증합니다.

5. Research Methodology

Research Design:

본 연구는 실험 및 수치 해석 연구(Ramos, 2012)를 참조 사례로 사용했습니다. 실험은 폭 1m, 길이 32.2m의 수로에서 직경 0.05m의 PVC 실린더를 사용하여 수행되었습니다. 하상 재료는 중간 직경 0.86mm의 석영 모래였습니다. 유량 64 L/s, 수심 0.20m 조건에서 평균 유속은 0.32 m/s였으며, 이는 입자 이동이 시작되는 임계 속도에 가깝습니다. 이 조건을 OpenFOAM을 사용하여 3차원 수치 시뮬레이션으로 재현하고, LES 난류 모델을 적용했습니다.

Data Collection and Analysis Methods:

수치 시뮬레이션 데이터는 오픈소스 시각화 툴인 Paraview를 사용하여 분석되었습니다. 속도장, 하상 전단 응력, 항력 계수 등의 데이터를 추출하여 이전 Fluent 연구 결과와 정성적, 정량적으로 비교 분석했습니다. 항력 계수는 OpenFOAM의 forceCoeffs functionObject를 사용하여 계산되었습니다.

Research Topics and Scope:

연구 범위는 고정된 하상(fixed bed)을 가진 단일 원형 교각 주변의 난류 유동 시뮬레이션에 한정됩니다. 두 가지 시나리오, 즉 세굴 전 평탄한 하상(Configuration A)과 평형 세굴 후의 하상(Configuration B)을 다룹니다. 이동상(mobile bed)이나 세굴 과정 자체의 시뮬레이션은 포함되지 않습니다.

6. Key Results:

Key Results:

  • Configuration A (평탄한 하상): OpenFOAM 시뮬레이션은 Fluent 시뮬레이션보다 더 크고 난류성이 강한 후류 와류를 보여주었으며, 이는 실험 관찰 결과와 더 일치합니다. 하상 전단 응력은 양쪽 시뮬레이션에서 유사한 크기(최대 약 2 Pa)를 보였으나, OpenFOAM에서 와류의 영향이 더 명확하게 나타났습니다. OpenFOAM으로 계산된 평균 항력 계수는 1.21로, Fluent의 0.68보다 높았지만 와류 방출로 인한 주기적인 진동을 잘 포착했습니다.
  • Configuration B (세굴된 하상): 세굴된 하상 지형으로 인해 후류 와류의 강도는 Configuration A보다 약해졌습니다. 이는 평형 상태에 도달했음을 시사합니다. 하상 전단 응력 값은 임계 전단 응력(0.5 Pa) 이하로 낮아져 추가적인 세굴이 발생하지 않는 평형 상태와 일치하는 결과를 보였습니다. 평균 항력 계수는 0.053으로 계산되었습니다.
  • 전반적 평가: 개선된 구조 격자를 사용한 OpenFOAM 시뮬레이션은 이전 Fluent 연구에 비해 결과의 질이 향상되었으며, 특히 와류 구조 모사에서 뚜렷한 개선을 보였습니다. 두 소프트웨어의 시뮬레이션 결과는 정성적, 정량적으로 합리적인 일치를 보였습니다.

Figure List:

  • Figure 1. Flow around pier (reproduced from Breusers and Raudkivi, 1991).
  • Figure 2. Scour cavity obtained in the laboratory test (Ramos, 2012).
  • Figure 3. Final bed geometry (Ramos, 2012).
  • Figure 4. a) Velocity vs. time on a random point in a turbulent flow b) Different eddies sizes in turbulence (reproduced from Fergizer, 2002).
  • Figure 5. 3D representation of the meshes (left – Ramos (2012); right – OpenFOAM).
  • Figure 6. Details of the mesh used in OpenFOAM simulations.
  • Figure 7. Schematic description of the geometry and its boundary-conditions.
  • Figure 8. Schematic description of the configurations and the planes where are presented the velocity fields (left – side view; right – top view).
  • Figure 9. Fluent results to velocity field (m/s) (t=20s; Configuration A) in the parallel plane to the bed channel, on half depth of the flow (Ramos, 2012).
  • Figure 10. OpenFOAM results to velocity magnitude (m/s) (t=20s; Configuration A) in the parallel plane to the bed channel, on half depth of the flow.
  • Figure 11. Fluent results to velocity magnitude (m/s) (t=20s; Configuration A) in the symmetry plane of the channel (Ramos, 2012).
  • Figure 12. OpenFOAM results to velocity magnitude (m/s) (t=20s; Configuration A) in the symmetry plane of the channel.
  • Figure 13. Fluent results to bed shear stress (Pa) (t=20s; Configuration A).
  • Figure 14. OpenFOAM results to bed shear stress (Pa) (t=20s; Configuration A).
  • Figure 15. Drag coefficient vs. flow time (Configuration A).
  • Figure 16. Fluent results to velocity magnitude (m/s) (t=20s; Configuration B) in the parallel plane to the bed channel, on half depth of the flow (Ramos, 2012).
  • Figure 17 – OpenFOAM results to velocity magnitude (m/s) (t=20s; Configuration B) in the parallel plane to the bed channel, on half depth of the flow.
  • Figure 18. Fluent results to velocity magnitude (m/s) (t=20s; Configuration B) in the symmetry plane of the channel (Ramos, 2012).
  • Figure 19. Fluent results to longitudinal velocity (m/s) (t=20s; Configuration B) in the symmetry plane of the channel (Ramos, 2012).
  • Figure 20. Fluent results to bed shear stress (Pa) (t=20s; Configuration B).
  • Figure 21. OpenFOAM results to bed shear stress (Pa) (t=20s; Configuration B).
  • Figure 22. Drag coefficient vs. flow time (Configuration B).

7. Conclusion:

OpenFOAM 툴박스를 사용하여 교각 주변의 유동에 대한 3차원 수치 시뮬레이션이 개발되었습니다. 평탄한 하상과 세굴된 하상 구성에 대해 LES 난류 모델을 사용하여 난류 유동을 시뮬레이션했습니다. 결과 중 일부는 상용 소프트웨어 Fluent를 사용하여 Ramos(2012)가 얻은 동일한 유동의 시뮬레이션과 비교되었습니다. Configuration A(평탄한 하상)에서는 Ramos(2012)의 연구와 달리 정규 격자를 사용하여 결과의 질을 명확하게 향상시켰습니다. 두 번째 구성(평형 세굴 깊이)에서도 정규 격자를 다시 적용했으며, 이를 위해 복잡한 모래 하상 형상으로 인한 문제를 방지하기 위해 불규칙한 하상을 단순화해야 했습니다. 와류, 특히 후류 와류 방출의 계산이 개선되었습니다. 속도장에서 말굽 와류의 출현은 명확하지 않지만, 세굴 구덩이 바로 상류에서 관찰된 교란은 긍정적인 신호입니다. 저자들의 관점에서, 각 소프트웨어(Fluent 및 OpenFOAM)로 얻은 시뮬레이션은 Ramos(2012)가 수행한 실험실 테스트에서 관찰된 예상 결과와 정성적, 정량적으로 합리적인 일치를 보입니다. 교각의 항력 계수 결과는 경험식(White, 2006)에서 주어진 값과 일치합니다.

8. References:

  • Batchelor, G. K. (1967) An introduction to fluid dynamics. University Press, Cambridge.
  • Breusers, H.; Raudkivi, A. J. (1991) Scouring. IAHR Hydraulic Structures Design Manual. A. A. Balkema, Rotterdam, Netherlands.
  • Dargahi, B. (1990) Controlling Mechanism of Local Scouring, J. Hydr. Engrg, ASCE, Vol 116, No.10
  • Deng, G.; Piquet, J. (1992) Navier-Stokes computations of horseshoe vortex flows. Int. J. Numer. Meth. Fluids, 15: 99-124.
  • Ferziger, J.H.; Peric, M. (2002) Computational Methods for Fluid Dynamics. Springer, Berlim.
  • Jongho L., Jungwoo K., Haecheon C. (2011) Sources of spurious force oscillations from an immersed boundary method for moving-body problems, Journal of Computational Physics, Vol. 230, Issue 7.
  • Melville, B.; Coleman S. (2000) Bridge Scour; Water Resources Publications LLC, (2000) Colorado, U.S.A.
  • Melville, B.; Raudkivi, A.J. (1977) Flow characteristics in Local Scour at Bridge Piers, J. Hydr. Research, ASCE, Vol 15, pp. 373-380.
  • Olsen, N.R.B.; Melaaen, M.C. (1993) Three-dimensional calculation of scour around cylinders. Journal of Hydraulic Engineering (ASCE), 119, 1048–1054.
  • Ramos, P.X. (2012, Modelação numérica do escoamento em torno de um pilar, MSc Thesis at Faculty of Engineering of University of Porto, Portugal.
  • Richardson, J.; Panchang, V. (1998) Three-dimensional simulation of scour-inducing flow at bridge piers. Journal of Hydraulic Engineering.
  • Roulund, A.; Sumer, M. Fredsøe, J. Michelsen, J. (2005) Numerical and experimental investigation of flow and scour around a circular pile. Journal of Fluid Mechanics, 534, pp 351-401.
  • Thwaites, B. (1960) Incompressible aerodynamics. University Press, Oxford.
  • White, F. (2006) Fluid Mechanics. McGraw-Hill Companies; 6th edition.

Expert Q&A: Your Top Questions Answered

Q1: OpenFOAM 시뮬레이션에서 구조 격자(structured mesh)를 사용한 이유는 무엇이며, 이전 연구와 비교하여 어떤 개선을 가져왔나요?

A1: 구조 격자는 셀이 정규적인 방식으로 정렬되고 분포되어 있어 데이터 처리 속도를 높입니다. 특히 본 연구에서는 교각 주변과 경계층 근처에 셀을 미세하게 배열하여 유동 구조를 더 정확하게 포착할 수 있었습니다. 논문에 따르면, 이 개선된 격자 구성은 후류 와류(wake vortex)를 더 현실에 가깝게 모사하는 데 결정적인 역할을 했으며, 이는 이전 연구(Ramos, 2012) 대비 결과의 질을 명확하게 향상시킨 핵심 요인입니다.

Q2: OpenFOAM으로 계산된 항력 계수(평균 1.21)가 경험식 값(약 0.73)과 상당한 차이를 보입니다. 이 결과를 어떻게 해석해야 하나요?

A2: 수치 자체의 차이도 중요하지만, 더 주목할 점은 Figure 15에서 보이는 항력 계수의 주기적인 진동입니다. 이 진동은 교각 후류에서 와류가 주기적으로 방출되면서 발생하는 힘의 변화를 물리적으로 올바르게 표현한 것입니다. 이 진동의 주파수와 진폭 정보는 교량 구조물의 피로 파괴나 공진 현상 등을 분석하는 동적 해석에 매우 중요한 데이터입니다. 따라서 평균값의 차이보다는 이러한 동적 거동을 성공적으로 포착했다는 점에 더 큰 의미를 둘 수 있습니다.

Q3: OpenFOAM 시뮬레이션에서 VOF(Volume of Fluid) 기법을 적용한 이유는 무엇인가요?

A3: 논문에 따르면, OpenFOAM 시뮬레이션은 20cm의 물층과 5cm의 공기층으로 구성된 2상 유동(two-phase flow)을 모사했습니다. VOF 기법은 이처럼 서로 섞이지 않는 두 유체(물과 공기)의 경계면을 추적하는 데 사용됩니다. 이는 이전 Fluent 연구에서 물만 시뮬레이션했던 것과 차이가 있으며, 자유 수면의 효과를 더 현실적으로 고려하기 위한 접근 방식입니다.

Q4: 복잡한 형상의 세굴된 하상(Configuration B)에 대해 어떻게 정규 격자를 적용했나요?

A4: 논문에서는 실제 세굴된 하상의 복잡하고 불규칙한 형상을 그대로 격자로 만드는 대신, 형상을 단순화(simplified)했다고 명시하고 있습니다. 이는 정규적인 셀 분포를 가진 구조 격자를 적용하기 위한 일종의 타협이었습니다. LES 모델이 많은 계산 시간을 요구하기 때문에, 격자 생성의 복잡성을 줄여 수치 해석 문제를 덜 복잡하게 만들기 위한 중요한 단계였습니다.

Q5: 세굴 과정에서 중요한 말굽 와류(horseshoe vortex)는 시뮬레이션에서 성공적으로 관찰되었나요?

A5: 논문에 따르면, 속도장(Figure 18, 19)에서 말굽 와류의 존재가 명확하게 나타나지는 않았습니다. 말굽 와류는 작고 복잡한 와류 시스템이라 계산하기 어렵기 때문입니다. 하지만 해당 영역에서 상당한 난류(turbulence)가 관찰되었으며, 저자들은 이를 “좋은 신호(a good signal)”라고 평가했습니다. 이는 비록 와류의 형태가 뚜렷하지 않더라도, 해당 물리 현상이 어느 정도 모사되고 있음을 시사합니다.

Conclusion: Paving the Way for Higher Quality and Productivity

이 연구는 교량 붕괴의 주된 원인인 교각 세굴 시뮬레이션에서 오픈소스 CFD 소프트웨어인 OpenFOAM의 높은 잠재력을 명확히 보여주었습니다. 특히, 체계적인 구조 격자 설계를 통해 상용 소프트웨어인 Fluent에 필적하거나 일부 측면(예: 와류 모사)에서는 더 현실적인 결과를 도출할 수 있음을 입증했습니다. 이는 고가의 상용 소프트웨어에 대한 효과적인 대안을 제시하며, 더 많은 엔지니어와 연구자들이 복잡한 수리 현상 해석에 접근할 수 있는 길을 열어줍니다.

“STI C&D는 최신 산업 연구 결과를 적용하여 고객이 더 높은 생산성과 품질을 달성할 수 있도록 최선을 다하고 있습니다. 이 논문에서 논의된 과제가 귀사의 운영 목표와 일치한다면, 저희 엔지니어링 팀에 연락하여 이러한 원칙을 귀사의 구성 요소에 어떻게 구현할 수 있는지 알아보십시오.”

(주)에스티아이씨앤디에서는 고객이 수치해석을 직접 수행하고 싶지만 경험이 없거나, 시간이 없어서 용역을 통해 수치해석 결과를 얻고자 하는 경우 전문 엔지니어를 통해 CFD consulting services를 제공합니다. 귀하께서 당면하고 있는 연구프로젝트를 최소의 비용으로, 최적의 해결방안을 찾을 수 있도록 지원합니다.

  • 연락처 : 02-2026-0442
  • 이메일 : flow3d@stikorea.co.kr

Copyright Information

  • This content is a summary and analysis based on the paper “NUMERICAL SIMULATION OF THE FLOW AROUND A PIER USING OPENFOAM” by “PEDRO RAMOS, JOÃO PEDRO PÊGO & RODRIGO MAIA”.
  • Source: ISBN 978-989-96479-2-3, 3rd IAHR Europe Congress, Book of Proceedings, 2014, Porto – Portugal.

This material is for informational purposes only. Unauthorized commercial use is prohibited. Copyright © 2025 STI C&D. All rights reserved.

PDF View