• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

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

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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1209-1213
• /
• 2007

하천 수치모델링을 통한 흐름, 오염물질 거동, 지형변화 해석 등은 효율적인 하천 수질 관리를 위해서 상당히 중요한 부분은 차지한다. 수질이나 지형변화를 보다 정확하게 예측하기 위해서는 하천 흐름 예측의 정확도 향상이 중요한 역할을 하게 된다. 본 연구는 평면 이차원 하상변동 및 수질예측 수치모형인 KU-RLMS 모형을 이용하여 낙동강 상류의 반변천 합류부의 흐름 특성을 규명하고, 수질 모형을 수행하기 위한 흐름 계산 결과를 제공하기 위해 수행하였다. KU-RLMS 모형은 하천 및 저수지의 국부적인 수리, 수질, 유사이동 해석을 위해 개발된 평면 이차원 비정상 수치모형이다. 직사각형 격자를 사용하는 유한차분법의 단점을 보완하기 위해, 흐름 계산을 위한 지배방정식은 3차원 Reynolds 방정식으로부터 수심적분된 2차원 연속방정식과 운동량방정식을 불규칙한 경계를 현실적으로 모사할 수 있는 직교곡선 좌표계로 변환한 방정식을 사용한다. 수치모형 적용을 위한 현황분석으로 안동 및 임하 조정지댐의 방류량, 안동 수위관측소의 자료를 분석하였다. 흐름 모형을 보정하기 위해 안동대교 지점에서 횡유속 분포를 측정하였으며, 이 결과를 사용하여 흐름 모형의 매개변수인 Manning 계수와 공간가중계수를 추정 및 검증하였다. 안동다목적댐과 임하다목적댐의 방류량을 고려하여 수치모의조건을 결정하였으며, 각 조건에 대한 흐름 변화 특성을 분석하였다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.8
• /
• pp.71-78
• /
• 2006

When a tracking filter may be designed in the Cartesian coordinate, the covariance of the measurement errors varies according to the range and the bearing of an interested target. In this paper, interacting multiple model based tracking filter is formulated in the Cartesian coordinate utilizing the analytic solution of the steady state Kalman filter, which can be able to consider the variation of the measurement error covariance. 100 Monte Carlo runs performed to verify the proposed method. The performance of the proposed method is compared with the conventional fixed gain and Kalman filter based IMM tracking filter in terms of the root mean square error. The simulation results show that the proposed approach meaningfully reduces the computation time and provides a similar tracking performance in comparison with the conventional Kalman filter based IMM tracking filter.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.10 no.5
• /
• pp.535-549
• /
• 1998

A numerical study for performance analysis of a counterflow type forced draft tower and natural draft cooling tower has been performed based on the method using the finite volume method with non-orthogonal body fitted and non-staggered grid system. For solving the coupling problem between water and air, air enthalpy balance, moisture fraction balance, water enthalpy balance, and water mass balance equations are solved with Navier-Stoke’s equations simultaneously. For the effect of turbulence, the standard k-$\varepsilon$ turbulent model is implied in this analysis. The predicted result of the present analysis is compared with the experimental data and the commercial software result to validate the present study, The predicted results show good agreement with the experimental data and the commercial software result. To investigate the influence of the cooling tower design parameters such as approach, range and wet bulb temperature, parametric studies are also peformed.

• Composites Research
• /
• v.12 no.2
• /
• pp.82-90
• /
• 1999

The mechanical properties of generally orthotropic materials are conventionally measured by performing off-axis tensile and flexure tests. However, the inevitable coupling between tension and shear in case of tensile test or bending and twisting in flexure test case induces nonuniform displacement and stress fields. Consequential stress concentration along the boundary of specimens would result in inaccurate modulus and underestimated strength. This paper proposes the variation of specimen geometry in terms of appropriate obliquity of loaded boundary. For the purpose, classical lamination theory is transformed into skewed coordinate, and characteristic equations for both of unidirectional and laminated composite specimens are formulated. Finite element analysis is employed to show the validity of the skewedness in tensile and bending test specimens.

• The Bulletin of The Korean Astronomical Society
• /
• v.37 no.2
• /
• pp.160.2-160.2
• /
• 2012

Unified State Model(이후 USM)은 Altman(1972)에 의해 처음 제안된 이후 Chodas(1981), Raol & Sinha(1985), Vittaldev et al.(2012) 등을 거치며 연구 발전되어 왔다. 이 모델은 공간상 6개 성분의 위치, 속도 벡터를 이용해 위성의 운동을 기술하는 기존 계산 방법과 달리 4개의 Quaternion 변수를 도입하여 위성의 위치를, 3개의 Hodograph 변수를 도입하여 위성의 속도를 각각 기술한다. USM의 장점은 직교좌표계로 표현된 위성의 위치, 속도 변수에 비해 USM 변수의 변화량이 상대적으로 작기 때문에 수치 계산 시 계산의 안정도가 높다. 또한 원궤도(${\omega}$ : undefined)와 적도면 궤도(i = 0, ${\Omega}$ : undefined) 계산 시에 나타나는 특이성(singularity) 문제가 발생하지 않는다. 본 연구에서는 USM 계산방법과 기존 방법에 의한 위성궤도 계산결과의 차이를 비교 분석하였다. 지구궤도 위성의 정밀계산을 위해 이체항 이외에 지구타원체 섭동항과 대기 항력에 의한 섭동항을 추가 적용하였다. 비구형 지구 중력 포텐셜에 의한 섭동은 J4항까지 고려하였으며, 대기 항력은 간단한 exponential 모델을 적용하였다. 또한 수치계산 시 적분 간격과 정밀도 차수를 조절하여 각 모델의 계산 안정성을 테스트하였다. 본 연구의 궤도계산 결과 USM 모델을 이용한 계산방법은 그 정밀성과 계산효율성이 매우 우수한 것으로 검증되었다.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.58-65
• /
• 1994

It is very important to understand characteristics of solution due to the variation of computational grid sizes, especially when turbulence model not incorporating wall-function is used. The present paper performs numerical investigation on the grid dependency of numerical solution for three dimensional turbulent flow field around a ship. In the present study a finite volume method with a modified sub-grid scale turbulence model and a numerically constructed non-orthogonal curvilinear coordinate system capable of conforming complex ship geometries are used. Numerical studies are then performed for a mathematical Wigley hull and the Series 60, $C_B=0.8$ hull forms. The results for various grid sizes are compared with each other and with measured data to show grid dependencies of numerical solutions.

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.