• 한국항공우주학회지
• /
• 제33권9호
• /
• pp.34-40
• /
• 2005

본 연구에서는 비정상 유동해석을 위한 CFD 코드의 개발을 위해 대각화 ADI 기법을 적용한 정상 해석기법과 내재적 이중시간 전진기법을 도입하였다. 정상상태 Navier-Stokes 방정식의 Jacobian 행렬은 비점성항에 대해서만 적용하였고 여기에 내재적 인공점성 연산자를 첨가하여 블록 5대각 행렬을 유도하였다. 시간단축을 위해 스칼라 5대각 행렬로 대체하였다. 가상시간에 대한 정상상태기법에 실시간에 대한 미분항이 포함된 새로운 잔류항을 정의하였다. 가상시간에 대해 수렴된 해로부터 실시간 해를 구하고 시간에 대해 적분을 수행하는 내재적 이중시간 전진기법을 이용한 비정상 Navier-Stokes 코드를 개발하였다. 이에 대한 검증으로 정지한 유체속에 진동하는 평판문제, 원기둥 후방의 주기적인 Karman 와류생성, 이중원호 익형주위의 충격파 진동문제등을 수치해석하여 이론치, 실험치, 타연구자의 계산결과와 비교, 분석하였다.

• 한국농림기상학회지
• /
• 제12권1호
• /
• pp.45-62
• /
• 2010

영면변위(d)는 거친 식생군락에 의해 운동량이 모두 흡수되어 군락 내의 대수적 풍속 프로파일이 0이 되는 높이를 말한다. 군락의 표면 거칠기의 구조를 나타내는 영면변위는 군락난류의 분석과 지표 스칼라 플럭스의 계산에 매우 중요하다. 본 단보에서는 Monin-Obukhov 상사이론에 기반을 두고 단일층에서 관측된 평균수평풍속 자료를 사용하는 두 가지 다른 방법을 사용하여 광릉침엽수림에서 d 값을 추정하였다. 관측지의 비균질성과 복잡성을 고려해서, 표면거칠기와 바람체계가 d에 미칠수 있는 영향을 살펴보기 위해, 자료를 매 $30^{\circ}$ 간격의 풍향별로 나누었다. 전반적으로 두 방법을 사용한 결과는 서로 비슷했는데, $d/h_c$ (여기서 $h_c$는 군락의 높이로서 약 ~23m)는 풍향에 따라 0.51~0.97의 범위를 보였다. 이러한 $d/h_c$의 값의 범위는 문헌에 보고되어 있는 범위(0.64~0.94)와 크게 다르지 않았으나, 다소 높은 쪽에 분포되어 있었다. 이러한 원인의 하나로는 관측이 두 방법의 전제인 Monin-Obukhov 상사이론이 성립하지 않는 거칠기아층에서 이루어졌기 때문인 것으로 사료된다. 따라서 관측높이가 거칠기 아층에 존재할 경우에는 단일층 풍속으로부터 영면변위를 추정하는 방법을 적용하고 그 결과를 해석하는 데에 세심한 주의가 필요하다.