• 대한기계학회논문집B
• /
• 제41권7호
• /
• pp.455-462
• /
• 2017

외부자기장에 의한 타원형 야누스 자성입자 사이의 자성 상호작용을 직접수치해석을 사용하여 분석하였다. 유한요소법에 기초한 가상영역법을 사용하여 입자계 유동해석을 수행하였고, 자기장 문제에서는 자성 포텐셜에 대한 지배방정식을 입자와 유체를 포함하는 전체영역에 대하여 풀어 자기장을 구하였다. 이 때 구해진 자기장으로부터 구한 맥스웰 응력을 사용하여 개별 입자에 작용하는 자기력이 계산된다. 입자의 운동과 최종적인 조립구조는 입자의 형상비, 개별 입자의 배향, 입자의 초기 분포에 크게 영향을 받는 것이 확인되었다. 또한 입자의 배향은 입자 주위의 유체 유동에도 영향을 주었다. 외부자기장에 의한 타원형 야누스 입자의 최종 조립구조는 앞서 언급한 인자들에 의해서 영향을 받은 것을 수치해석을 통해 확인할 수 있었다.

• 한국음향학회지
• /
• 제37권4호
• /
• pp.163-173
• /
• 2018

본 논문에서는 비행 중 비행체 표면에 작용하는 음향하중 예측을 수행하였다. 비행 중 음향하중은 비행체 표면의 압력 변동에 의해 발생한다. 기존의 비행 중 음향하중 예측방법은 반경험적 방법으로 이론과 실험 결과를 기반으로 도출한 경험식을 활용한다. 하지만 경험식의 입력 값으로 사용되는 비행체 주변 유동특성 및 경계층 파라미터를 매번 실험을 통해 얻는 것에는 한계가 있다. 따라서 본 논문에서는 전산유체해석(Computational Fluid Dynamics, CFD) 결과를 반경험적 방법과 혼합하는 하이브리드 방법을 이용하여 비행 중 비행체에 작용하는 음향하중을 예측하였다. Cone-cylinder-flare 형상 비행체에 대해 아음속, 천음속, 초음속, 최대동압도달(Maximum dynamic pressure, Max-q) 시점의 비행 환경에 대한 음향하중 예측을 수행하였다. 하이브리드 방법 적용 시 전산유체해석결과를 기반으로 한 경계층 끝단 영역 판단 방법에 대해 비교하였고 여러 연구자에 의해 제시된 경험식에 따른 음향하중 예측결과를 비교하였다.

• 대한기계학회논문집B
• /
• 제39권5호
• /
• pp.425-433
• /
• 2015

본 논문에서는 외부에서 균일한 직류전기장이 인가될 때 점성유체에 자유롭게 잠겨있는 단일 입자가 근처의 비전도성 평면 벽과의 상호작용 때문에 유발되는 2차원 유전영동 운동에 대하여 수치연구를 수행하였다. 특히 입자-유체 경계면에서 불연속적으로 급격히 변화하는 전기전도도를 가진 Maxwell 방정식을 해석하고 전기장을 구한 후 Maxwell 응력텐서를 적분하여 입자에 작용하는 유전영동 힘을 계산하였다. 해석 결과 전기장이 벽과 평행하게 인가될 때 입자는 항상 반발력이 유도되어 벽으로부터 멀어지는 방향으로 유전영동 운동이 발생하였으며, 그 운동특성은 입자와 벽 사이 간격과 입자의 전도도에 따라 크게 달라졌다. 운동 강도는 입자와 유체의 전도도가 서로 같으면 사라지나, 전도도가 서로 다르면 그 차이가 클수록 강도는 증가하였다.

• 대한조선학회지
• /
• 제22권3호
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• 공업화학
• /
• 제10권2호
• /
• pp.324-329
• /
• 1999

본 연구에서는 고압하에서 이산화탄소와 이소프로필 알코올과의 이성분계에 대한 상평형데이터를 얻기 위해 실험을 수행하였다. 실험장치는 정지형 (static type)을 사용하였으며, 정확도를 실험하기 위해 $80^{\circ}C$에서 이산화탄소-이소프로필 알코올계의 실험을 수행하여 Radosz의 실험결과와 비교하였다. 이산화탄소-이소프로필 알코올과의 이성분계 상거동 실험은 온도 40, 60, 80, 100, 그리고 $120^{\circ}C$에서 실험하였으며, 이때 압력은 41~133 bar 범위였다. 이산화탄소-이소프로필 알코올계에 대해 동일한 압력에서 용해도는 온도가 증가함에 따라 증가함을 알 수 있다. 또한 순성분 이산화탄소와 이소프로필 알코올의 증기압을 서로 연결하는 혼합물 임계곡선을 나타내었다. 본 연구에서 실험한 결과를 statistical associating fluid theory (SAFT)상태방정식에 의해 계산하였으며, 그 결과 온도에 독립적인 두 파라미터에 의해 곡선을 결정하였다.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.

• 대한조선학회논문집
• /
• 제47권3호
• /
• pp.377-387
• /
• 2010

In this study, the turbulent free surface around KVLCC1 employed in the circular motion test simulation is numerically calculated using a commercial CFD(Computational Fluid Dynamics) code, FLUENT. Also, hydrodynamic forces and yaw moments around a ship model are calculated during the steady turning. Numerical simulations of the turbulent flows with free surface around KVLCC1 have been carried out by use of RANS equation based on calculation of hydrodynamic forces and yaw moments exerted upon the ship hull. Wave elevation is simulated by using the VOF method. VOF method is known as one of the most effective numerical techniques handling two-fluid domains of different density simultaneously. Boundary layer thickness and wake field are changed various yaw velocities of ship model during the steady turning. The calculated hydrodynamic forces are compared with those obtained by model tests.

• 한국유체기계학회 논문집
• /
• 제4권1호
• /
• pp.22-29
• /
• 2001

Centrifugal fans are widely used in industrial practices but the noise generated by these machines causes one of the most serious problems. In general, the centrifugal fan noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cutoff in the easing. However, only a few researches have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan, and to calculate the effects of small vanes that are attached in original impeller - Splitter impeller. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The splitter impeller changes the acoustic characteristics as well as performance. Two-splitter type impeller and splitter impeller which splitter locates in jet region are good for acoustic characteristics.

• International Journal of Fluid Machinery and Systems
• /
• 제8권3호
• /
• pp.169-182
• /
• 2015

When simulating the dynamic behaviour of a hydro power plant, it is essential to have a good representation of the turbine behaviour. The pressure transients in the system occurs because the flow changes, which the turbine defines. The flow through the turbine is a function of the pressure, the speed of rotation and the wicket gate opening and is, most often described in a performance diagram or Hill diagram. In the Hill diagram, the efficiency is drawn like contour lines, hence the name. A turbines Hill diagram is obtained by performance tests on scaled model in a laboratory. However, system dynamic simulations have to be performed in the early stage of a project, before the turbine manufacturer has been chosen and the Hill diagram is known. Therefore one have to rely on diagrams for a turbine with similar speed number. The Hill diagram is drawn through measured points, so for using the diagram in a simulation program, one have to iterate in the diagram based on curve fitting of the measured points. This paper describes an alternative method. By means of the Euler turbine equation, it is possible to set up two differential equations which represents the turbine performance with good enough accuracy for the dynamic simulations. The only input is the turbine's main geometry, the runner blade in- and outlet angle and the guide vane angle at best efficiency point of operation (BEP). In the paper, simulated turbine characteristics for a high head Francis turbine, and for a reversible pump turbine are compared with laboratory measured characteristics.

• Journal of applied mathematics & informatics
• /
• 제29권5_6호
• /
• pp.1583-1602
• /
• 2011

A special system of partial differential equations (PDEs) occur in a natural way while studying a class of irrotational inviscid fluid flow problems involving infinite channels. Certain aspects of solutions of such PDEs are analyzed in the context of flow problems involving multiple layers of fluids of different constant densities in a channel associated with arbitrary bottom topography. The whole analysis is divided into two parts-part A and part B. In part A the linearized theory is employed along with the standard Fourier analysis to understand such flow problems and physical quantities of interest are derived analytically. In part B, the same set of problems handled in part A are examined in the light of a weakly non-linear theory involving perturbation in terms of a small parameter and it is shown that the original problems can be cast into KdV type of nonlinear PDEs involving the bottom topography occurring in one of the coefficients of these equations. Special cases of bottom topography are worked out in detail and expressions for quantities of physical importance are derived.