• 대한조선학회논문집
• /
• 제30권2호
• /
• pp.76-85
• /
• 1993

자유표면과 수면하 vortex 유동과의 상호작용으로 생기는 수면위의 특이 형상을 실험과 수치해석으로 얻어진 결과에 근거하여 논하였다. 먼저 2차원적인 실험이나 해석적인 연구로 이해될 수 있는 "striation"의 발생을 살펴 본다. "Scar"라 불리우는 striation과 수직한 형상은 2차원 vortex filament와 vortex pair 생성시 생겨나는 residual vortex의 상호작용으로 설명하였다. 간단하면서도 중요한 상호작용의 역학을 연구할 목적으로 vortex ring이 자유표면과 수직하게 충돌할 경우에 나타나는 현상들을 실험적으로 논하였고 이때 상호작용의 초기단계를 축대칭 문제로 하여 수치해석적으로 얻어진 결과와 비교하였다. 상호작용의 후기단계는 vortex reconnection을 포함하는 복잡한 과정인데 이때 수면에서의 짧은파 생성을 vortex ring의 수면과의 경사진 충돌로 설명하였다. 또한 위의 수면과 생성을 간단한 해석적 model로 재구성하여 비교해 보았다. 끝으로 수면위에 존재하는 오염물질이 수면하 vortex 유동에 미치는 영향을 수치해석적인 방법으로 free slip및 rigid wall 조건의 경우와 비교하여 살펴보았다.

• 한국산학기술학회논문지
• /
• 제16권11호
• /
• pp.7990-8000
• /
• 2015

식물공장은 환경적 조건을 조절하여 계절이나 장소에 관계없이 농산물을 일정하게 생산하는 식물재배 시스템이다. 식물공장에 있어서 내부의 열 유동은 식물공장의 중요한 변수이다. 본 논문에서는 다층 재배선반을 갖는 하이브리드 식물공장 내의 열유동 특성을 수치 시뮬레이션을 통해 연구하였다. 열 유동 특성의 수치해를 얻기 위해 유한체적법(Finite Volume Method)을 이용하였다. 수치해석 모델에 있어서는 저 레이놀드수(low Reynolds number) ${\kappa}-{\epsilon}$ 난류모델을 이용하였으며, 해석방법으로는 비압축성 점성유동 영역과 압력경계 조건을 사용하였다. 수치해석에 있어, 3종류의 유입 공기속도와 재배선반의 위치변화에 따른 식물공장 내의 열유동 특성을 상용 프로그램인 Solid Works Flow simulation을 사용하였다. 수치해석을 통한 결론은 다음과 같다. 첫째, 송풍기를 통한 식물공장 내의 유입 공기속도를 1.6 m/s로 하였을 때, Case 3의 재배선반 배치가 비교적 균일한 온도 분포를 나타내었고, LED광만을 사용한 Case 1의 배치는 시험베드(test bed) TB2와 TB3 영역의 낮은 온도분포로 인하여 서큘레이터 등을 통한 추가적인 온도 제어가 필요할 것으로 판단되었다. 둘째, 식물재배 LED 인공광(DYLED47)의 청색광과 적색광의 비율을 1:1로 100% 구동 시, 재배선반의 평균 온도가 약 $3^{\circ}C$ 증가하였다. 본 논문은 한국산학기술학회 논문지 심사용 투고요령입니다.

• 대한기계학회논문집B
• /
• 제24권1호
• /
• pp.21-28
• /
• 2000

Although commercial PIV systems have been widely used for the non-intrusive velocity field measurement of fluid flows, they are still under development and have considerable room for improvement. In this study, a single-frame double-exposure PIV system using a high-resolution CCD camera was developed. A pulsed Nd:Yag laser and high-resolution CCD camera were synchronized by a home-made control circuit. In order to resolve the directional ambiguity problem encountered in the single-frame PIV technique, the second particle image was genuinely shifted in the CCD sensor array during the time interval dt. The velocity vector field was determined by calculating the displacement vector at each interrogation window using cross-correlation with 50% overlapping. In order to check the effect of spatial resolution of CCD camera on the accuracy of PIV velocity field measurement, the developed PIV system with three different resolution modes of the CCD camera (512 ${\times}$ 512, lK ${\times}$ IK, 2K ${\times}$ 2K) was applied to a turbulent flow which simulate the Zn plating process of a steel strip. The experimental model consists of a snout and a moving belt. Aluminum flakes about $1{\mu}m$ diameter were used as scattering particles for the liquid flow in the zinc pot and the gas flow above the zinc surface was seeded with atomized olive oil with an average diameter of 1-$3{\mu}m$. Velocity field measurements were carried out at the strip speed $V_s$=1.0 m/s. The 2K ${\times}$ 2K high-resolution PIV technique was significantly superior compared to the smaller pixel resolution PIV system. For the cases of 512 ${\times}$ 512 and 1K ${\times}$ 1K pixel resolution PIV system, it was difficult to get accurate flow structure of viscous flow near the wall and small vortex structure in the region of large velocity gradient.

• 대한조선학회논문집
• /
• 제32권1호
• /
• pp.83-93
• /
• 1995

사각형의 MHD 추진 덕트 내부유동을 수치해석과 실험적인 방법으로 연구하였다. 수치해석연구에서는 전기장과 자기장의 영향하에 있는 비압축성 3차원 통전유체에 대하여 유한 차분법으로 계산하였다. 수치계산의 결과 전자기력이 약할때 층류유동의 전형적인 포물선 유동 형태가 전극 부근에서 M자 형상으로 변하였고, 균일분포의 일정 전자기력하에서 MHD 덕트 내부 압력은 입구에서부터 하류로 나아감에 따라 선형적으로 증가하였다. 실험에서는 MHD 추진 덕트 내부 유동을 해석할 수 있는 실험장치를 제작 하였으며, MHD 추진은 전류를 변화 시킴으로써 추력을 쉽게 조정할 수 있음을 알았다. 또한 MHD 덕트내 유동방향의 압력구배는 전자기력에 비례하여 증가하며 수치계산 결과와 잘 일치하였다.

• Korea-Australia Rheology Journal
• /
• 제15권3호
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• 한국산학기술학회논문지
• /
• 제18권3호
• /
• pp.545-550
• /
• 2017

대동맥이나 협착된 경동맥에서는 심장수축기에 간헐적으로 난류현상이 발생하고 있으며, 혈액의 점성특성으로 인해 기존 난류모델로는 정확한 해석이 어려운 실정이다. 혈류는 점탄성 유체의 성질을 가지고 있어 유체의 전단 변형률 증가에 따라 점도가 감소하는 점탄성 유체이며, 이러한 점탄성 유체는 난류 유동시 저항 감소 현상이 발생한다. 기존의 난류해석 모델들은 점성변화가 없는 뉴턴 유체에 적합한 모델들이 대부분이기 때문에, 점탄성 유체의 저항 감소 현상을 고려한 비뉴턴 유체 해석에 적합한 난류 모델개발이 필요하다. 본 논문은 난류 모델 가운데 수렴성이 좋고 해석시간이 짧은 표준 $k-{\varepsilon}$ 모델을 기반으로 저항 완충 함수를 이용하여 비뉴턴 유체의 저항감소 현상을 해석할 수 있는 수정된 난류모델을 제시하였으며, 이를 기존 난류모델들과 비교하여 제시된 난류 모델을 검증하였다. 새로 제시된 수정된 난류모델은 벽함수 및 점성저층을 고려하지 않았기 때문에 해석시간이 대폭적으로 감소하였으며, 적은 격자수를 이용하여 효율적으로 비뉴턴 유체의 난류 현상을 해석할 수 있기 때문에 향후 혈류해석 및 점탄성유체 해석에 적용할 예정이다.

• Journal of Nutrition and Health
• /
• 제26권6호
• /
• pp.793-803
• /
• 1993

This study was designed to measure viscosity, osmolality and in vitro flow rates via nasogastric tubes for 6 types of commercially available and 9 hospital-blenderized enteral solutions and to examine the effect of viscosity and osmolaility of enteral formula on the flow rates in gravity drip administration. Each solution was infused through 18, 16, 14, 12 French sizes of silicone rubber tube. Flow rates were measured six times at $25^{\circ}C$ using formula bags and drip sets hung at a uniform height on a intravenous drip stand with tube uniformly positioned in collecting container. Viscosity ranged widely from 16.0 to 195.5 cps with mean, 64.61$\pm$64.42 for hospital-blenderized formula while mean viscosity of commercial formula was 7.60$\pm$4.84 cps. Mean osmolality of commercial formula and hospital-blenderized formula were 370$\pm$100.80, 540.33$\pm$89.37 mOsm/kg respectively. There was negative relationship between viscosity of formula and flow rates through tubes but no significant relationship between flow rates and osmolalty. Some of hospital-blenderized formula was too viscous to be infused througth tube with gravity drip administration and the recipe of formula requires to be modiifed. On the other hand, commercial formula with the low viscosity flows too rapidly with large bore size tubes. Smaller size of tube must be selected for hyperosmolar solution to decrease possible side effects associated with tube feeding. Two kinds of regression equations for flow rates obtained according to viscosity and tube sizes were also presented for the purpose of practical uses. In conclusion, this study emphasizes that viscosity of fomula, osmolality, patient's tolerance and comfort, caloric density should be considered in the selection of tubes for gravify drip administration.

• 한국전산유체공학회지
• /
• 제18권2호
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• 한국전산유체공학회지
• /
• 제18권1호
• /
• pp.83-90
• /
• 2013

OpenMP is becoming more and more useful as a simple parallel processing paradigm on SMP (Shared Memory Multi-Processors) computing environment with the development of multi-core processors. However, very few data is available publically regarding the OpenMP performance in CFD (Computational Fluid Dynamics). In the present study a CFD test suite is prepared for the performance evaluation of OpenMP on various multi-core systems. The test suite is composed of two-dimensional numerical simulations for inviscid/viscous and reacting/non-reacting flows using three different levels of grid systems. One to five test runs were carried out on various systems from dual-core dual threads to 16-core 32-threads systems by changing the number of threads engaged for each test up to 80. The results exhibit some interesting results and the lessons learned from the tests would be quite helpful for the further use of OpenMP for CFD studies using multi-core processor systems.