• Proceedings of the KIEE Conference
• /
• 2001.11a
• /
• pp.284-287
• /
• 2001

유한체적법(Finite Volume Method)은 기계공학분야에서 열, 유체현상을 수치해석하는 방법으로 널리 쓰이고 있다. 열전달(Heat Transfer)문제의 지배방정식은 에너지방정식(Energy equation)으로써 그 중 순수 전도문제의 경우 지배방정식은 라플라스(Laplace)방정식의 형태를 띄고 있다. 전계해석의 지배방정식 역시 라플라스 방정식의 형태이다. 수학적으로 동일한 지배방정식을 갖는 다는 것은 물리적으로 같은 현상을 나타내고 있다는 것을 의미하며 이러한 점에 착안을 하여 본 논문에서는 유한체적법을 사용하여 2차원 모델에 대한 전계해석을 수행하였다.

• Journal of the KSME
• /
• v.34 no.9
• /
• pp.688-697
• /
• 1994

유체유동이나 열전달 그리고 물질전달 (물질의 혼합 및 확산) 또는 이들 현상이 복합적으로 나 타나는 각종 기계의 설계와 성능 해석을 하기 위해서는 그 현상을 지배하는 편미분 방정식들의 해를 수치적으로 구해야 한다. 유동 상태가 충류 유동인 경우는 지배 방정식의 수가 알고자 하는 미지변수 즉 속도, 압력, 온도, 농도 등의 개수와 같고 또한 이들 변수들의 변동이 그리 심하지 않기 때문에 적절한 수치 해법을 사용하면 그 해를 구할 수 있다. 그러나 난류유동의 경우에는 변수들이 시간상으로 또한 공간적으로 대단히 심하게 변동(fluctuation)하기 때문에 공 학적으로 우리가 원하는 정보들, 즉, 표면 마찰저항이나 양력, 얼전달 계수, 물질 확산계수 등을 현재 수준의 전자계산기로 계산하는 데는 계산시간이 엄청나게 소요될 뿐만 아니라 변수 저장 메모리도 과도하게 차지하기 때문에 실제적인 계산 방법이 되지 못하고 있다. 이러한 이유로 변수들의 순간 변화 상태를 나타내는 지배 방정식들을 해석하는 대신에 이들 지배 방정식의 시 간평균을 취하여 유도한 난류 방정식들을 사용하게 된다. 그러나 이 시간 평균 과정에서 파생 되는 또 다른 미지의 난류 변수들 때문에 난류 지배 방정식에 있어서는 그 지배 방정식의 개수 보다 미지 변수의 개수가 많아져서 난류 지배 방정식을 풀기 위해서는 시간평균 과정에서 나타난 난류 변수들을 원래 있던 미지 변수들의 함수나 방정식의 형태로 가정할 필요가 있게 되는데 이 가정되는 함수 관계들을 난류 계산 모형이라고 한다. 난류 계산 모형은 물리적인 통찰과 직관에 의해서 실용적인 형태로 가정되기도 하지만 최근에는 논리적으로 엄격한 모형 원칙에 따른 수 학적인 방법으로 유도되고 있는데 이 글에서는 일반 독자들이 쉽게 이해할 수 있도록 마하수가 낮은 2차원 비압축성 난류 유동을 예로 들어 x-y 직교 좌표계에서 표현되는 난류 계산 모형들을 소개하고 앞으로듸 발전 방향을 개관하며 현재의 응용 사례들을 예로 들어 모형의 성능을 비교 하여 보기로 한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.660-665
• /
• 2004

본 연구에서는 비선형적인 메커니즘을 갖는 수문현상의 불확실성을 해석하고자 하는 목적으로 새로운 개념의 지배방정식이 유도된다. 제안된 모형의 불확실성은 토양 특성치의 공간적 변동성에 기인하고 있는 것으로 가정하여, 유도된 방정식은 Fokker-Planckl 방정식의 형태를 가지고 있다. 실제 유역단위에서 토양 내 수분 흐름의 연직방향 흐름을 모의하기 위해 미소단위에서 유도된 Richards 방정식은 토양의 공간적 변동성으로 말미암아 불확실한 매개변수를 갖는 비선형 추계학적 편미분방정식의 형태를 갖게 된다. 이는 먼지 수직 방향적분을 통하여 단순화된 비선형 추계학적 상미분방정식으로 전환되고, 이렇게 전환된 비선형 추계학적 상미분방정식은 다시 추계학적 Liouville 방정식을 이용하여 선형 추계학적 편미분방정식으로 전환되어진다. 최종적으로 cumulant 급수방법을 이용하여 상기 방정식을 선형 결정론적 편미분방정식으로 전환시킴으로써, 강우 시 토양 내 수분 침투현상을 모형화할 경우 유역단위에서 토양의 공간적 변동성을 설명할 수 있는 지배방정식을 유도할 수 있다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1989.07a
• /
• pp.163-172
• /
• 1989

강우로 인한 급경사면의 흐름 해석은 지배방정식인 연속방정식과 운동량 방정식에 의해 수치적으로 해석된다. 이에대한 해석 방법에는 수문학적인 방법과 수리학적인 방법이 있는데, 본 연구에서는 세인트 베난트 방정식을 초기 및 경계조건에 맞게 푸는 수리학적인 방법을 이용하였고, 일반적인 완경사면에서 그 영향이 작아 무시되었던 경사의 영향을 고려한다. 수치해석 과정은 지배방정식에 무차원 매개변수를 도입하여 6점 음해법을 이용하여 해석하였다. 수치실험은 마닝(Manning)의 조도계수 n을 변화시키면서, 각각 검사에서 강우에 의한 평형상태에 도달한 이후의 수심과 유량, kinematic flow 수, 후르드 수, 평형상태 도달시간등 변수 상호간의 관계를 고찰하였다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.49 no.6
• /
• pp.449-456
• /
• 2021

In this paper, we present a strategy to extend solution capability of an existing low Mach number preconditioned compressible solver to incompressible flows with a little modification. To this end, the energy equation that is of the same form of the total energy equation of compressible flows is used. The energy equation is obtained by a linear combination of the thermal energy equation, the continuity equation and the mechanical energy equation. Subsequently, a modified artificial compressibility method in conjunction with a time marching technique is applied to these incompressible governing equations for steady flow solutions. It is found that the Roe average of the common governing equations is equally valid for both the compressible and incompressible flow conditions. The extension of an existing compressible solver to incompressible flows does not affect the original compressible flow analysis. Validity for incompressible flow analysis of the extended solver is examined for various inviscid, laminar and turbulent flows.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.57-57
• /
• 2011

지진해일은 진행속도가 빠르고 파장이 길며 파형의 변화 없이 먼 거리를 진행 할 수 있어 주변지역은 물론 멀리 떨어진 지역에도 심한 범람피해를 야기시킨다. 지진해일의 일반적인 특징으로 장파와 단파가 합성되어 있고 먼 거리를 전파할 경우 분산효과의 역할이 중요하게 된다. 특히 우리나라의 동해안에 영향을 주는 지진해일은 단주기파 성분이 강하고 파장에 비해 먼 거리를 전파하기에 분산을 고려하는 선형 Boussinesq 방정식을 지배방정식으로 사용하는 것이 바람직하다. 하지만 지금까지의 지진해일 전파모의를 위한 모형은 선형 Boussinesq 방정식의 복잡한 계산과 계산시간이 길다는 단점 때문에 선형 천수방정식을 지배방정식으로 사용하고 분산효과는 수치분산을 이용하여 고려해왔다. 지진해일 해석 시 일반적으로 사용되어 오던 기존의 leap-frog 유한차분 모형(Imamura et al., 1988; 조용식, 1996)은 지배방정식으로 선형 천수방정식을 사용하고 파의 분산효과는 수치분산을 이용하여 고려하므로 정해진 시간 간격에 대해 수심에 따라 격자 간격을 적절히 선택해야 하는데 수심이 복잡하게 변하는 경우 격자간격 조정이 불가능하여 분산효과를 정도 높게 고려할 수 없다. 이 문제점을 해결하기 위하여 윤성범 등(2004)은 파동방정식의 인위적인 분산항을 이용하여 Boussinesq 방정식의 분산효과를 고려할 수 있는 능동적인 분산보정기법을 제안하였고 Cho et al.(2007)는 일정한 수심에서 수치적인 분산오차가 Boussinesq 방정식의 물리적인 분산항을 대체하도록 수심, 격자 간격 및 계산 시간 간격 사이의 관계식을 유도하고 Boussinesq 방정식의 분산항과 일치하는 수치분산을 이용하여 실용적인 분산보정기법을 개발하였다. 이에 Ahn(2010)은 현재 컴퓨터의 계산 능력이 향상되어 선형 Boussinesq 방정식을 직접 차분하여 계산하는데 무리가 없다고 판단하여 선형 Boussinesq 방정식을 직접 차분한 모형을 개발하였다. 본 연구에서는 기존의 원해 지진해일 전파모의에 이용되어왔던 선형 천수방정식에 수치분산을 고려한 모형 대신 선형 Boussinesq 방정식의 유한차분 모형을 제안하였으며 기존의 선형 Boussinesq 방정식 모형의 격자와 수심간의 제약을 없애기 위해 차분 기법을 달리 한 2차 정확도의 유한차분 모형을 제안하였다. 검증을 위하여 선형 Boussinesq 방정식의 해석해(Carrier, 1991)와 비교하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.177-177
• /
• 2022

기후변화 및 해안 구조물의 증가 등 여러 원인이 연안침식 및 해안선 변화와 같은 연안의 지형변화를 가속하고 있다. 빠르게 변화하는 연안의 지형변화예측 및 대응책 강구를 위해서는 연안의 유사이송 현상에 대한 신속한 예측이 필요하다. 본 연구에서는 GPU 엔진 기반 파랑해석모형인 Celeris Advent를 활용하여 실시간으로 연안의 유사이송 모의가 가능한 수치모형을 개발하였다. Celeris Advent는 GPU의 병렬코어를 활용해 실시간 연산과 GUI를 통한 사용자와의 실시간 상호작용이 가능한 모형이다. 지배방정식은 확장형 Boussinesq 방정식에 유사이송방정식을 양방향 결합하여 구성하였고, 지배방정식에는 하이브리드 유한체적-유한차분 수치기법을 적용하여 이송항은 유한체적법(Kurganov & Petrova, 2007), 소스항은 유한차분법을 통해 이산화하여 해석한다. 유사이송방정식은 수심적분형 이송확산방정식에 침식 및 퇴적 플럭스를 반영하는 소스항을 결합하여, 이송항 및 확산항을 통해 유사의 이송/확산을 고려함과 동시에 소스항을 통해 하상과의 상호작용을 고려하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.40-41
• /
• 2011

동해를 전파하는 지진해일은 세계적으로 다른 지역에서 발생하는 지진해일과 비교하였을 때 상대적으로 파장이 짧고, 이에 비해 먼 거리를 전파한다. 그러므로 동해에서 발생한 지진해일의 전파에 대한 해석을 수행할 때 물리적인 분산효과가 매우 중요하다. 따라서 지배방정식으로 분산 효과가 충분히 고려된 선형 Boussinesq 방정식을 사용한다. 기존의 연구에서는 leap-frog 기법을 사용하여 선형 천수방정식을 차분할 때 발생하는 수치분산항에 분산 보정계수를 이용하여 선형 Boussinesq 방정식의 물리적 분산항과 같은 형태로 나타나도록 유도하여 수치모의를 수행하였다. 그러나 기존에 사용한 지배방정식은 수심이 일정하다는 가정을 통하여 유도된 것으로, 수심에 변화가 있는 실제 지형을 통과하는 지진해일에 대한 수치모의를 수행한 결과의 정확도에 문제가 생길 수 있다. 본 연구에서는 기존의 연구에서 발생할 수 있는 수심 변화에 따른 오류를 개선하기 위하여 바닥 지형이 1차원으로 변한다는 가정을 이용하여 지배방정식을 유도하였으며, 이로 인해 발생하는 수심 변화가 고려된 항을 기존의 분산보정기법에 추가하였다. 그리고 적용성을 높이기 위하여 수치모의 기법의 제한을 최소화하는 연구를 진행하였다. 본 연구에서 제안한 수정 기법이 수심이 변화하는 지형을 전파하는 지진해일 수치모의 과정에서 경사에 대한 분산효과가 충분히 고려되는지 확인하기 위하여 Gaussian hump를 이용한 가상 지진해일을 원형 천퇴 지형에 통과시켰다. 본 연구에서 사용한 지형을 통과하는 Gaussian hump에 대한 해석해를 구하는 방법이 존재하지 않으므로, Boussinesq 방정식을 직접 차분하여 푸는 FUNWAVE를 사용하여 동일한 조건 하에서 수치모의를 수행하였다. 비교 결과를 통하여 본 연구에서 제안한 기법의 정확도 향상을 확인하게 되면, 실제 지형을 통과하는 지진해일의 수치모의에 대한 활용성을 높일 수 있을 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.98-98
• /
• 2019

중력류 또는 밀도류는 주변 유체에 비해 상대적으로 밀도가 큰 유체가 밀도차에 의한 추진력으로 흐르는 것이다. 중력류의 수치모델링에는 두 가지 어려움이 있다. 즉, 적합한 지배방정식을 구성하여 적용하는 것 그리고 난류의 영향을 합리적으로 반영하는 것이다. 기존 중력류 해석을 위한 지배방정식들은 유체의 연속방정식과 운동량 방정식 그리고 밀도 또는 농도의 이송방정식을 조합하여 구성된다. 이들 지배방정식을 이용한 연구들은 대부분 두 유체 사이의 밀도차가 충분히 작아서 밀도 변동(variations)의 영향은 오로지 부력항에서만 유지된다는 Boussinesq 근사에 근거를 둔다. 그리고 이송방정식에서 밀도 또는 농도의 확산계수을 점성계수의 함수로 표현하기 위해서 Schmidt 수를 이용한다. 수치모델링에서 Schimdt 수는 상수값을 적용하지만, 이 값은 밀도의 연직방향 경사에 근거한 부력빈도(buoyancy frequency)와 난류량의 따라 큰 차이를 보이는 것으로 알려져있다. 한편, 표준 통계학적 난류모델과 벽함수를 적용한 수치모델링은 초기 중력에 의해서 무너지는(slumping) 단계를 넘어 관성력으로 추진되는 단계와 점성 효과가 지배적인 단계에서는 정확도에 현저히 낮아지기 때문에 대부분 큰와모의(large-eddy simulation, LES) 또는 DNS(direct numerical simulation)수준의 고해상도(high-resolution) 해석기법을 적용하여 공학적인 문제에 적용하는 데는 한계가 있다. 이 연구에서는 Boussinesq 근사와 Schmidt 수를 사용하지 않으며, LES 보다 적용이 용이한 DES (detached-eddy simulation)기법을 조합한 다상흐름 수치모델을 적용하여 중력류를 해석을 시도하였다. 수치해석결과를 실험값과 함께 기존 수치모델링 기법으로 구한 수치해와 비교분석하여 이 연구에서 개발 및 적용된 수치모델링 기법의 적용성을 평가한다.

• Journal of Korean Society of Steel Construction
• /
• v.22 no.3
• /
• pp.241-251
• /
• 2010

This paper presents the development of a new analytical solution to the governing differential equation for isotropic annular sector plates subjected to uniform loading in a three-dimensional polar coordinate system. The 4th order governing partial differential equation (PDE) was converted to an ordinary differential equation (ODE) by assuming the Levy-type series solution form and the subsequent mathematical operations. Finally, a series-type solution was assembled with homogeneous and nonhomogeneous solution parts after operating real values and complex conjugates derived from the characteristic equation. To demonstrate the convergence rate and the accuracy of the featured method, several examples with various sector angles were selected and solved. The deflections and internal moments in the example annular sector plates that were obtained from the proposed solution were compared with those obtained from other analytical studies and numerical analyses using the finite element analysis package program, ABAQUS. Very good agreement with the results of other analytical and numerical methodologies was shown.