• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.71-80
• /
• 1996

An iterative time marching procedure for solving incompressible turbulent flow has been applied to the flows around a high speed train including cross-wind effects. This procedure solves three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using first-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. Turbulent flows have been modeled by Baldwin-Lomax turbulent model. To validate present procedure, the flow around a high speed train at zero yaw angle was simulated and compared with experimental data. Generally good agreement with experiments was achieved. The flow fields around the high speed train at 9.2°, 16.7°, and 45° of yaw angle were also simulated.

• 한국전산유체공학회지
• /
• 제16권3호
• /
• pp.8-14
• /
• 2011

Although the geometry of circular cylinder is simple, the flow is complicate because of the flow separation and vortex shedding. In spite of many numerical and experimental researches, the flow around a circular cylinder has not been clarified even now. It has been known that the unsteady vortex shedding from a circular cylinder can vibrate and damage a structure. Lock-on phenomenon is very important in the flow around an oscillating circular cylinder. The lock-on phenomenon is that when the oscillation frequency of the circular cylinder is at or near the frequency of vortex shedding from a stationary cylinder, the vortex shedding synchronizes with the cylinder motion. This phenomenon can be recognized by the spectral analysis of the lift coefficient history. At the lock-on region the vortex is shedding by the modulated frequency to the body frequency. However, the vortex is shedding by the mixed frequencies of natural shedding and forced body frequency in the region of non-lock-on. In this paper, it was analyzed the relation between the frequency of rotary oscillating circular cylinder and the vortex shedding frequency.

• 대한조선학회논문집
• /
• 제29권3호
• /
• pp.33-44
• /
• 1992

2차원 및 3차원 sloshing 문제의 해석을 위해 수치기법들을 적용하여 보았다. 2차원 탱크 내의 유동을 source 분포법과 유한차분법을 이용해 해석하였고 실험과 선형해의 결과와 비교하였다. 이 방법들의 계산결과에는 큰 차이를 발견하지 못하였고, 다만 유한차분법을 이용하는 경우 탱크 내의 내부재에 대한 고려가 용이하나 계산시간이 많이 소요되었고 source 분포법을 이용하여 해석하는 경우에는 3차원으로의 확장이 용이하나 내부재의 고려와 대진폭 운동에 대해서는 적용하기 힘들 것으로 판단된다. 3차원 문제는 source 분포법을 이용하여 해석하였고 계산모델은 사각형 및 구형탱크이다. 그러나, 이러한 sloshing 해석기법을 보다 실용적으로 사용하기 위해서는 유체의 충격력에 대한 많은 연구가 있어야 할 것으로 사료된다.

• 대한기계학회논문집
• /
• 제16권11호
• /
• pp.2150-2158
• /
• 1992

본 연구에서는 복잡한 유동형태를 지닌 충돌분사류에 대한 유동특성을 연구하 기 위하여 단순화된 실험모델로써 형상이 동일한 두 원형분류의 충돌에 의한 충돌분류 의 혼합현상 및 유동구조 등을 질량유량비의 변화에 따라 유체역학적으로 구명하고자 하였으며, 본 연구 결과는 연소기관에서의 연소효율 증대 및 구조개선등의 공학적 응 용을 위한 기본자료로 활용하고, 이론적 연구에 의한 난류의 유동구조 및 유동특성 에 대한 타당성 입증과 이론적 모델의 보완을 위한 실험자료로 이용하고자 한다. 충돌유동에 영향을 미치는 주요인자는 노즐직경, 충돌각, 충돌질량유량비, 온도, 밀도 등이며, 이 인자들 중에서 충돌질량유량비와 출돌각이 충돌후 형성되는 난류혼합유동 에 지배적인 영향을 미치므로, 본 연구에서는 두 원형분류의 충돌질량유량비를 가변할 수 있는 장치를 고안하였으며, 두 분류의 충돌각을 45˚로 고정하고, 고속측과 저속측 노즐의 질량유량비를 1.0, 0.8, 0.6, 0.4로 설정하여 질량유량비에 따른 혼합 유동구 조의 구명을 위한 실험적인 연구를 수행하였다. 충돌후의 혼합유동의 특성을 연구하 기 위하여 유동중심궤적, 유동반폭, 유동단면, 2차원 및 3차원 유동장, 평균속도분포 등을 온라인 컴퓨터시스템을 이용하여 측정분석하였다.

• 한국추진공학회지
• /
• 제23권1호
• /
• pp.24-32
• /
• 2019

유체 추력벡터 제어에서 이중목 노즐 개념의 이용 가능성을 조사하기 위하여, 초음속 노즐에서 수치해석을 수행하였다. 수치해석 검증에서 SST $k-{\omega}$ 난류모델을 사용하여 실험결과를 잘 구현하였다. 광범위한 노즐 압력비와 분사 압력비에서 편향각도, 시스템의 전체 추력비 및 추력 효율을 조사하였다. 본 연구에서 이중목 노즐의 추력벡터제어 시스템의 성능 변화는 2차원 계산영역에서 명확하게 설명되었다. 본 연구에서 얻어진 결과들은 유체추력벡터제어 분야에 중요한 기초자료를 제공할 것이다.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1997년도 춘계학술발표회논문집(1)
• /
• pp.298-303
• /
• 1997

실제의 격납건물의 구조는 하부 원통형의 구조를 가지는 영역과 상부 돔 형태와 굴뚝 형태의 구조를 가지는 영역으로 나눌 수 있다. 하부 원통형의 구조만을 고려한다면, 고온의 철제 벽면과 콘크리트 벽면 사이의 gap 크기에 비해서 원통의 반지름이 상대적으로 매우 큰 값을 가지기 때문에 2차원 무한평판으로 가정하는 것이 가능하다. 그러나 돔 및 굴뚝 영역에서는 높이가 높아질수록 돔 단면직경이 감소하고 굴뚝 영역도 유동단면적이 작은 원통의 구조를 가져 2차원 무한평판의 가정에 많은 무리가 따른다. 앞에서 명시한 세 가지의 격납건물 형태에 있어서 ASPWR의 경우는 굴뚝을 포함한 영역까지도 무한평판으로 가정하는 것이 가능하나(돔에서의 열전달 단면적이 하부의 열전달 단면적에 비해 매우 작다는 가정을 한다면) 나머지 AP600과 HWRF의 격납건물에 있어서는 상부까지도 무한평판 가정을 사용하는 것에는 무리가 있다. 본 연구에서는 일반적인 유체해석 코드인 FLUENT V4.3을 이용하여 실제 격납건물 구조에 대한 분석을 시도하여 무한평판 구조에 대한 가정이 과도한 열전달량을 예측하고 있음을 확인하였다.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2014년도 제46회 동계 정기학술대회 초록집
• /
• pp.231.1-231.1
• /
• 2014

열 플라즈마(thermal plasma) 는 저온 플라즈마(cold plasma)와 달리 이온과 전자와 중성입자들이 충분한 에너지 교환으로 인해 열평형 상태를 가진다. 열 플라즈마를 생성 시킬 때 전극 사이에서 아크방전을 시켜 제트 형태로 플라즈마를 발생시키는 것을 플라즈마 토치(plasma torch)라고 한다. 이러한 플라즈마 토치는 화학 원소 분해, 강판 절단, 유해 기체 분해 등으로 널리 사용되고 있다. 본 연구에서는 플라즈마 토치를 수치적으로 해석하여 플라즈마의 특성을 알아보았다. 수치해석적 접근방법으로 열 플라즈마는 LTE (local thermodynamic equilibrium)을 가정하였으며 one-fluid 이론을 적용하였다. 이때 사용된 코드는 DCPTUN으로서 $C^{+}^{+}$로 작성된 열플라즈마 유동의 특성해석 코드인 동시에 SIMPLE 알고리즘을 이용한 유체 코드이다. 시뮬레이션은 2차원 축대칭이며 정렬격자계 및 비정렬격자계 모두에서 사용이 가능하도록 되어있다. 또한 맥스웰 방정식을 통해 electromagnetic field를 풀도록 하여 RF 시뮬레이션이 가능하도록 하였다. 이와 같은 열 플라즈마 시뮬레이션을 통해서 플라즈마 토치의 특성을 알아보았다.

• 한국전산유체공학회지
• /
• 제21권2호
• /
• pp.25-31
• /
• 2016

In this paper, a study on the development of OpenFOAM grid generation program for two-dimensional flow analysis is described. By using the pre-processor(eMEGA) of EDISON_CFD system, grids for OpenFOAM flow calculation were obtained. Resultant two-dimensional grids were used to calculate flow fields by applying simpleFoam, one of the OpenFOAM's popular solvers, and the obtained flow results were compared with theoretical and experimental data available. Also grids generated by present program were compared with grids by a commercial pre-processor Pointwise for the purpose of verification. Verification work includes three cases(single block, O-type single block, and multi block grid), and all results show reasonable matches. According to the current achievement, it can be concluded that OpenFOAM grid can be constructed conveniently by using eMEGA with GUI.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2000년도 제14회 학술강연논문집
• /
• pp.5-5
• /
• 2000

액체 추진 로켓 엔진의 고주파 연소 불안정 관련 이론은 대체로 연소기 내부의 음향 공명 모드와 분무 연소 과정의 상호 작용을 구동 메커니즘으로 전제하며 Rayleigh Criterion의 재해석에 기초하여 불안정성 평가를 위한 매개변수를 도입하고 연소 불안정성을 예측한다. 여기에는 음향장 분석 이론, 음향 불안정 이론, 연소응답 및 기화반응 이론 등이 포함된다. 본 연구에서는 LOX/RPl 추진제 조합의 액체 추진 로켓 엔진 연소기를 대상으로 다차원 순수 음향장 해석과 연소-음향장 분석을 통해 대상 엔진의 고주파 연소 불안정 특성을 예측하였다. 수동 제어 기기인 음향공 설치에 따른 연소기의 음향장 및 연소-음향장의 특성 변화를 고찰하고 위 결과를 종합하여 음향공의 연소 불안정 억제 성능 및 대상 엔진의 연소 불안정성을 평가하였다. 연소기 형상 및 음향공 설치에 따른 다차원 순수 음향장 해석은 상용코드인 ANSYS를 사용하여 수행하였다. 내부 유체는 압축성, 비점성 유체로 유체의 평균 유동은 무시하며 위치에 관계없이 균일한 물성치를 부여하였다. 정상상태 연소과정을 가정하고 평형 화학을 이용한 분석 결과로부터 연소 기체의 관련 물성치를 결정하였다. 연소기 길이 방향, 반경 방향, 원주 방향 격자점들의 음향 특성을 주파수 영역에 대해 해석하고 3차원 음향 모드 형상을 토대로 음향장을 분석하였다. 연소-음향장 해석은 음향 불안정 이론 중 n- $\tau$ 2 매개변수 기법을 사용하였다. 연료 액적의 분무 연소 과정을 1차원적으로 가정하고 정상상태의 평형 화학 계산 결과를 이용하여 엔진의 연소면을 1차원적으로 설정하였다. 상류 연소응답과 중립 안정 곡선을 토대로 대상 엔진의 연소 불안정 특성을 분석하였다.구 분석 결과 기술적 문제점으로는 배기 가스온도가 낮은데 따른 출구 부분의 Bearing, Sealing이 문제가 될 수 있다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14)

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.27-33
• /
• 2009

In this paper, the flow past a rotary oscillating circular cylinder is simulated. The high-order and high-resolution numerical schemes with the characteristic boundary conditions are used for the compressible Navier-Stokes equation. The frequencies of rotating oscillation are $0.19\;{\leq}\;S_f\;{\leq}\;0.25$ for the maximum angular $\theta_{max}=10^{\circ}$ and $17^{\circ}$. The flow conditions are Mach number of 0.3 and Reynolds number of 1000. At Lock-on and Non-lock-on region which are defined by the relation between the vortex shedding frequency and the oscillating frequency, the drag and lift coefficient are analyzed.