• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 1997년도 제9회 학술강연회논문집
• /
• pp.17-17
• /
• 1997

극초음속 여객기와 군사용 항공기에 대한 수요가 증가함에 따라서 새로운 개념의 다양한 추진기관이 연구가 진행되고 개발되어 왔다. 초음속 항공기의 속도 영역은 마하 10-20 정도가 되는데 이 속도 한계를 극복하기 위하여 초음속 연소 램제트 엔진(SCRamjet; Supersonic Combustion Ramjet)이 제안되었다. 스크램 제트를 개발하기 위해서는 연료와 산화제의 혼합 효율 문제, 화염의 안정화 문제, 벽면의 냉각에 관한 문제 등 몇 가지 기본적인 문제들을 해결해야 한다. Univ of Michigan에서 실험한 연소기를 모델로 본 연구에서는 연료와 공기의 혼합에 관한 수치 연구를 수행하였다. 다원 혼합기체에 관한 축대칭 Navier-Stokes 방정식을 지배 방정식을 이용하였고 비평형 화학반응식을 고려하였다. 공간 차분에는 유한 체적법을 이용하였다. 대류 플럭스 항은 Roe의 Upwind FDS 기법을 사용하여 차분하였고 점성항에는 중심 차분법을 이용하였다. 시간 적분법으로는 근사 자코비안과 LU분할 기법을 이용한 완전 내재적 방법이 쓰였다. 난류 모델로는 Mentor에 의해 제안된 2 방정식 k-$\varepsilon$/k-$\omega$ 혼합모델을 사용하였다. 유동장이 실험에서의 찍은 사진과 유사한 모습의 충격파 간섭을 수치 모사하였고 수소가 확산되는 모습과 함께 노즐 lip 주위의 재순환 영역에 대해서 살펴볼 수 있었다.

• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

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

2개의 경사 충격파와 하나의 수직 충격파로 초음속 유동을 압축하는 초음속 공기 흡입구의 수치적 연구를 수행하였다. 지배방정식으로는 Navier-Stokes방정식을 사용하였고 난류모델로는 SST 모델을 사용하였다. 지배방정식의 점성항 계산에는 중심차분법을 사용하였고 대류항 계산에는 풍상차분법인 Roe의 FDS기법을 MUSCL기법과 결합하여 이용하였다. 유한 체적법을 이용하여 차분된 방정식은 LU분할 기법을 이용한 완전 내재적 방법으로 2차 정확도 시간 적분으로 비정상 과정의 연구를 수행하였다. 흡입구 배압을 정해주어야 하는 어려움을 해결하기 위해 흡입구 후면에 노즐을 달고 노즐의 면적을 조절하여 배압이 형성되도록 하였다.(중략)

• Journal of Biosystems Engineering
• /
• 제10권1호
• /
• pp.24-38
• /
• 1985

경사진 밀폐 공간에서 마주 보는 두 벽면의 온도 차로 인하여 발생되는 자연 대류 현상은 여러 공학 분야에서 볼 수 있는 중요한 열전달 현상으로서, 최근 들어 평판형 태양열 집열기를 설계하려는 사람들에게 많은 관심의 대상이 되고 있다. 평판형 태양열 집열기의 경우 덮개판으로 부터의 대류 열손실을 감소시킴으로서 집열 효율을 높일 수 있을 뿐만 아니라 사용목적에 따라 소형 집열기를 제작할 수 있어 경제적으로 유리하게 될 것이다. 밀폐된 공간에서 최초에 정지 상태에 있는 얇은 유체층을 하부에서 가열시켜 주면 열팽창 현상이 일어나고, 이것에 의한 부력이 점도나 열전도도 등의 안정화 요인을 극복할 수 있을 정도로 커지면 System이 불안정하게 되어 자연 대류 현상이 수반되며 이 때문에 열전달율이 급격히 증가하게 된다. 이러한 현상의 지배 방정식은 연립 비선형 편미분 방정식으로 특수한 경계 조건외에는 일반적으로 해석적 해를 구하기가 어렵기 때문에 실험적 연구가 많이 실시되어 왔고 이들 결과의 대부분은 전반적인 열전달 특성치만을 구하는데 집중되어 왔다. 본 연구에서는 수치 해석법의 하나인 유한 차분법을 도입하여 이차원으로 가정한 경사진 평판형 밀폐 공간에서의 자연 대류 현상의 지배 방정식을 유한 차분화시켜, $$2.74{\times}10^3\leq_-Gr\leq_-2.0{\times}10^6$$, Pr=0.73, $$15^{\circ}\leq_-a\leq_-150^{\circ}$$, 종횡비는 1, 2, 3, 5, 9에 대하여 정상 상태에서의 해를 구하면서 수치적으로 실험하였다. 본 연구에서 얻어진 결론을 요약하면 다음과 같다. (1) 해석적으로 구하기 어려운 경사진 밀폐 공간에서 자연대류현상의 지배 방정식을 유한 차분법으로 해결할 수 있으며, 대류항과 확산항을 따로 고려한 유한차분법이 효과적임을 확인하였다. (2) 저온과 고온 벽면에서의 온도를 각각 균일하게 놓고 단변으로 이루어진 벽면은 완전히 절연되어 있는 경우에 대하여 수치해를 구한결과, 이전의 해석적 및 실험적 결과와 일치하였으며, 시간의 경과에 따른 온도 및 유선의 변화를 현상학적으로 관찰할 수 있었다. (3) 평균 열전달 계수에 미치는 경사각의 효과를 살펴본 결과 종횡비가 1 인 경우 경사각이 $45^{\circ}$에서, 종횡비가 2, 3, 5, 9인 경우 경사각이 $60^{\circ}$에서 각각 평균 열전달 계수 최대치가 나타났다. (4) Ra수(Rayleigh number) 가 증가될수록, 경사각에 상관없이 평균 열전달 계수도 증가되었다. 한편 Ra수 및 경사각의 변화에 따라 종횡비가 증가될수록 평균 열전달 계수는 경사각이 $90^{\circ}$인 경우를 제외하고는 감소됨을 볼 수 있었다. 경사각이 $90^{\circ}$인 경우, 평균 열전달 계수는 종횡비가 2 인 곳에서 최대치를 얻을 수 있었으며, 종횡비가 계속 증가될수록 평균 열전달 계수는 점차 감소되어짐을 볼 수 있었다.

• 한국컴퓨터그래픽스학회논문지
• /
• 제15권1호
• /
• pp.1-6
• /
• 2009

이 논문은 비격자메쉬에서의 비물리적 감쇄현상을 감소시키는 새로운 CIP 대류방법을 제안한다. 제시되는 방법은 물리량과 물리량의 미분정보를 함께 사용하여 기존의 선형방법보다 더 정확한 물리현상을 반영할 수 있다. 비격자메쉬에서 해결하기 힘든 CIP의 비대류항을 풀기 위해 기존의 격자 메쉬에서 사용되었던 유한차분법 방법과는 다른 유한체적법을 사용하였다. 그리고 대류항의 보다 나은 정확도를 위해 여러가지 보간법에 대해 정확도 실험을 수행하였으며, 선형방법과의 비교 실험을 통해 보다 적은 비물리적 감쇄현상을 보임을 확인할 수 있었다.

• 조명전기설비학회논문지
• /
• 제19권4호
• /
• pp.108-116
• /
• 2005

1계 미분항이 포함되는 미분방정식의 수치해를 구하고자 할 때 중앙차분을 사용한 유한차분법이나 Galerkin법을 사용한 유한요소법은 그 해가 매우 불안하여 요소분할을 세밀하게 하여야만 해를 얻을 수 있다. 이러한 해의 불안 정성이 일어나는 이유는 대류항의 크기가 커질수록 후류에서의 경계조건이 해의 급격한 변화를 요구하는데 수치해가 급격한 변화에 적응하지 못하기 때문이다. 이러한 문제를 해결하기 위해 1970년대부터 upwind법이 개발되어 왔다. 본 논문은 1계 미분항이 표현되는 속도기전력이 발생하는 전자계 문제를 유한요소법을 이용하여 해석할 때 발생하는 해의 진동 문제를 해결하기 위해 Heinrich에 의해 제안된 upwind법을 적용하였다.

• 대한조선학회논문집
• /
• 제31권3호
• /
• pp.80-86
• /
• 1994

자유수면을 항주하는 선박에 의하여 발생되는 비선형 조파현상을 해석하기 위한 수치해석법을 개발하였다. 유동은 비점성, 비압축성으로 가정하고 선체 및 자유수면 형상과 일치하는 좌표계의 생성을 위하여 타원형 편미분방정식을 수치해석하여 물체적합 좌표계를 생성하였으며 변환된 정규격자 물체적합 좌표계에 대한 Euler방정식을 유한차분법(Finite Difference Method)을 이용하여 계산하였다. 수치해석을 위하여 시간에 대한 미분항은 전진차분, 공간에 대한 미분항은 중심차분법으로 이산화하였고 대류항에는 수치계산의 안정을 위해 인위적인 소산(dissipation)항을 첨가하였다. 자유수면의 형상은 매 시간 단계마다 자유수면 경계조건들을 만족하도록 다시 계산되었고 격자점들은 자유수면형상의 변화에 적합하게 다시 생성되도록 하였으며 압력에 대한 Poisson방정식은 반복연산법에 의하여 풀고 그 결과를 이용하여 속도를 외삽하였다. 개발된 수치해석법의 검증을 위해 수식선형인 Wigley 모형에 대한 계산을 Fn=0.250-0.408에 대하여 수행하고, 그 결과를 실험 결과와 비교하여 잘 일치함을 보였다.

• 대한조선학회논문집
• /
• 제35권1호
• /
• pp.1-14
• /
• 1998

2차원 천수방정식에 기초한 조류계산 수치모델을 개발하였다. 복잡한 해안선의 입력을 단순화할 수 있는 직각격자를 사용하였으며, 방정식은 효율이 뛰어난 것으로 알려져 있는 인수분해음해법(Implicit Factorization Method)에 의하여 수치적으로 풀었다. 대류항은 상류차분법으로 처리하여 수치소산이 자동적으로 들어가도록 하였으며, 따라서 수치불안정에 의한 발산을 억제할 수 있었다. 모델의 검증을 위하여 선형화된 천수방정식을 풀고 계산결과를 해석해와 비교하였으며 잘 일치하고 있음을 확인하였다. 수치모델의 응용 예로서 인천항부근의 경인운하 개발예정인 수로에서의 조류 계산결과를 보였다.

• 대한기계학회논문집B
• /
• 제30권7호
• /
• pp.671-681
• /
• 2006

This study focuses on the development of numerical procedure to analyze the nonlinear combustion instabilities in liquid rocket engine. Nonlinear behaviors of acoustic instabilities are characterized by the existence of limit cycle in linearly unstable engines and nonlinear or triggering instability in linearly stable engines. To discretize convective fluxes with high accuracy and robustness, approximated Riemann solver based on characteristics and Euler-characteristic boundary conditions are employed. The present procedure predicts well the transition processes from initial harmonic pressure disturbance to N-like steep-fronted shock wave in a resonant pipe. Longitudinal pressure oscillations within the SSME(Space Shuttle Main Engine) engine have been analyzed using the pressure-sensitive time lag model to account for unsteady combustion response. It is observed that the pressure oscillations reach a limit cycle which is independent of the characteristics of the initial disturbances and depends only on combustion parameters and operating conditions.

• 대한조선학회논문집
• /
• 제37권4호
• /
• pp.19-30
• /
• 2000

선수와 선미에 벌브를 가진 현대적인 상선 주위의 난류유동을 해석하기 위해 유한체적법을 이용한 RANS 방정식의 해법이 개발되었다. 복잡한 선수미의 격자계 생성을 위해 다중블록 선체표면 격자계 생성 기법이 도입되었고, 타원형 미분 방정식의 해를 이용하여 O-H 형태의 공간 격자계생성 기법이 이용되었다. 지배방정식의 대류항과 확산항은 각각 QUICK과 중앙차분법을 이용하여 근사되었고, 속도-압력 연성을 위해 SIMPLEC법이 채택되었다. 개발된 방법은 KRISO 3600TEU 컨테이너선에 적용하여 유선가시화 실험 및 반류계측 결과와 비교하였다. 계산 결과는 선수와 선미에 벌브를 가진 현대적인 선형 주위의 유동 해석을 위해서도 이러한 수치계산 방법이 적용될 수 있음을 보여주고 있다.