• 동력기계공학회지
• /
• 제2권3호
• /
• pp.21-26
• /
• 1998

Numerical simulations for various fluid flows require enormous computing time during iterations. In order to solve this problem, several techniques have been proposed. A SOR method is one of the effective methods for solving elliptic equations. However, it is very difficult to find the optimum relaxation factor, the value of this factor for practical problems used to be estimated on the basis of expertise. In this paper, the implication of the relaxation factor are translated into fuzzy control rules on the basis of the expertise of numerical analysers, and fuzzy controller incorporated into a numerical algorithm. From two cases of study, Poisson equation and cavity flow problem, we confirmed the possibility of computational acceleration with fuzzy logic and qualitative reasoning in numerical simulations. Numerical experiments with the fuzzy controller resulted in generating a good performance.

• Journal of Mechanical Science and Technology
• /
• 제16권11호
• /
• pp.1522-1539
• /
• 2002

In the present study, an explicit algebraic stress model is shown to be the exact tensor representation of algebraic stress model by directly solving a set of algebraic equations without resort to tensor representation theory. This repeals the constraints on the Reynolds stress, which are based on the principle of material frame indifference and positive semi-definiteness. An a priori test of the explicit algebraic stress model is carried out by using the DNS database for a fully developed channel flow at Rer = 135. It is confirmed that two-point correlation function between the velocity fluctuation and the Laplacians of the pressure-gradient i s anisotropic and asymmetric in the wall-normal direction. Thus, a novel composite algebraic Reynolds stress model is proposed and applied to the channel flow calculation, which incorporates non-local effect in the algebraic framework to predict near-wall behavior correctly.

• Journal of Mechanical Science and Technology
• /
• 제18권7호
• /
• pp.1140-1149
• /
• 2004

The absorption characteristics of water vapor into a LiBr-H$_2$O solution flowing down on finned inclined surfaces are numerically investigated in order to study the absorbing performances of different surface shapes of finned tubes as an absorber element. A three-dimensional numerical model is developed. The momentum, energy, and diffusion equations are solved simultaneously using a finite difference method. In order to obtain the temperature and concentration distributions, the Runge-Kutta and the Successive over relaxation methods are used. The flat, circular, elliptic, and parabolic shapes of the tube surfaces are considered in order to find the optimal surface shapes for absorption. In addition, the effects of the fin intervals and Reynolds numbers are studied. The results show that the absorption mainly happens near the fin tip due to the temperature and concentration gradient, and the absorbing performance of the parabolic surface is better than those of the other surfaces.

• 대한설비공학회지:설비저널
• /
• 제10권3호
• /
• pp.185-219
• /
• 1981

To analyze two dimensional incompressible laminar natural convection in a rectangular enclosure heated from below and cooled by a horizontal ceiling and two vertical walls, he primitive Navier-Stokes equations and the energy equation were solved numerically in time dependent form by a marker and cell method. A successive over-relaxation method for the elliptic portion of the problem and an explicit method for the parabolic portion were applied for the range of Grashoff number of $5{\times}10^3\;to\;5{\times}10^4$ to get the transient and steady state dimensionless temperature and velocity profiles. For the range of aspect ratio $L/H{\leqq}2.4$ in which only a pair of convection rolls exists mean Nusselt number calculated are as follows : $$N_{NU}0.89\;N_{Gr}^{0.2}(H/L)^{0.45}$$ By path lines drawn by marker particle trajectories roll number of cellular motion were observed for various aspect ratio of the enclosure.

• 한국에너지공학회:학술대회논문집
• /
• 한국에너지공학회 1995년도 추계학술발표회 초록집
• /
• pp.50-55
• /
• 1995

향상되어진 93-PCGC-2는 기존의 PCGC-2와 같이 미분탄 연소를 포함하는 다양한 반응성흐름과 비반응성 흐름을 설명하기 위해 2차원 정상상태 모델로 제시되어 졌다. 93-PCGC-2는 실린더형의 축 대칭계에 응용되어질 수 있고, 난류(Turbulence)는 유체역학식과 연소기구 양쪽을 위해 고려되어졌으며, 불연속 세로좌표 방법(Discrete Ordinates Method)을 이용하여 기체, 벽 및 입자들로부터의 복사열(Radiation)을 모사하였다. 입자상은 입자 무리들의 평균 경로들을 따라 해석하는 Lagrangian계의 해석법으로 모델화되어졌다. 석탄의 팽윤(Swelling)과 촤의 반응성에 관한 부모델과 더불어 새롭게 일반화된 석탄 탈휘발화 부모델 (FG-DVC)도 첨가되어졌다. 비균일 반응기구는 확산과 화학반응 둘 모두를 고려하였다. 주요 기상반응은 국부 순간 평형을 가정하여 모델화하였다. 그래서 반응속도는 혼합의 난류속도에 의해 제한되어진다. Thermal NOx과 Fuel NOx의 유한속도 화학론(Finite Rate Chemstry)에 대한 부모델은 화학반응속도론와 난류성의 통계치를 통합하여 만들어져 있다. 기상은 반복적인 line-by-line기교에 의해 풀려지는 elliptic partial differential equation으로 묘사되어진다. 수치적인 안정을 고려하기 위해 under-relaxation이 이용되어졌다. 이렇게 코드화된 93-PCGC-2는 연소를 위해 모사되어졌다. 또한 더 나아가 이 수치모델의 활용범위는 미분탄의 가스화에도 활용되어질 것으로 기대되어진다.

• 한국정보처리학회논문지
• /
• 제2권4호
• /
• pp.535-542
• /
• 1995

이 논문에서 우리는 CM-5와 같은 초대형병렬컴퓨터에서 대형 이산선형체제를 풀기 위한 준비행렬로써 두 가지를 소개한다. 대다수의 초대형병렬컴퓨터들은 프로세서간의 통신을 메세지패씽(messagepassing)에 의존하는데 현재의 기술수준하에서는 이 통신속 도가 실수계산속도에 비해 매우 느리므로 종래의 메모리공유컴퓨터에서와는 달리 데이 터통신량을 최소화하는 알고리듬이 요구된다. 블록 SOR에 다중색채기법을 가미한 알고 리듬이 그 한 예로써 우리는 이를 CM-5에서 구현한 결과 N=512x512 행렬에서 프로세서 의 수가 16에서 512의 범위 하에서 50%의 효율을 실현하였다. 반면 종래의 효율적인 병렬 준비행렬로 알려진 AKI알고리듬은 방대한 량의 데이터통신 때문에 매우 열등한 결과를 보여준다.

• 대한기계학회논문집B
• /
• 제24권2호
• /
• pp.204-213
• /
• 2000

This study deals with jet impingement, which is extensively used in the process industries to achieve intense heating, cooling or drying rates and also widely employed as a test flow for turbulent models due to its complex flow configuration, on a flat plate by numerical methods. In this calculation, the finite volume method was employed to solve the Navier-stokes equation based on the non-orthogonal coordinate with non-staggered variable arrangement. To get a better understanding for the fluid flow and heat transfer characteristics of the turbulent jet impingements, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model was adapted and compared with the experimental data and the result of standard $k-{\varepsilon}$ turbulent model. Numerical calculations were carried out with various flow rates, nozzle to plate distances. In the case of the axisymmetric jet impingement on a flat plate, $k-{\varepsilon}-{\overline{v^{'2}}}$ turbulent model showed better agreement with the experimental data than the standard $k-{\varepsilon}$ turbulent model in the prediction of the mean velocity profiles, the turbulent velocity profiles. the turbulent shear stress and the heat transfer rate. The highest heat transfer rate can be obtained when the impingement occurs within the potential core..

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• 한국해안해양공학회지
• /
• 제1권1호
• /
• pp.22-30
• /
• 1989

파랑은 천해역을 진행하면서 천수효과에 의한 굴절 및 구조물에 의한 회절, 흐름의 영향, 반사등에 의한 변형을 일으킨다. 굴절 및 회절에 의한 파랑변형을 예측하기 위하여 흐름이 고려된 비정상 상태의 Kirby(1986a)의 식을 이용하여 회절이 고려된 기초방정식을 유도하고, 좌표계를 파향선식에 의한 직교곡선 좌표계로 변환한 후 양해법으로 반복계산을 행하여 해를 구하는 수치 모형을 개발하였다. 수치 모형의 적용성을 검토하기 위하여 다음과 같은 조건에서 계산하였으며, 기존의 계산결과 및 실험결과와 비교 검토하였다. 1) 일정 수심상에서 원형 천퇴 주변의 파랑변형 2) 사면상에서 타원형 천퇴 주변의 파랑변형 3) 사면상에서 파랑이 입사각을 갖고 입사할 경우 구조물 주변의 파랑변형 본 모형의 수치해를 검토한 결과 기존의 연구결과와 잘 일치하였다.

• 정보처리학회논문지A
• /
• 제10A권2호
• /
• pp.137-140
• /
• 2003

최근에 CG 반복법을 이용하여 레이레이 계수를 최소화함으로써 대칭행렬의 내부고유치를 구하는 방법이 개발되었다 그리고 이 방법은 병렬계산에 매우 적합하다. 적절한 준비행렬의 선택은 수렴속도를 향상시킨다. 우리는 본 연구에서 이를 위한 병렬준비행렬들을 비교한다. 고려된 준비행렬들은 Point-SSOR, 다중색채하의 ILU(0)와 Block SSOR이다. 우리는 128개의 노드를 가진 CRAY-T3E에서 구현하였다. 프로세서간의 통신은 MPI 리이브러리를 사용하였다. 최고 512$\times$512 행렬까지 시험하였는데 이 행렬들은 타원형 편미분방정식의 근사화에서 얻어졌다. 그 결과 다중색채 Block SSOR이 가장 성능이 우수한 것으로 판명되었다.