• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.85-88
• /
• 2010

This study has been mainly motivated to numerically model the supercritical mixing and combustion processes encountered in the liquid propellant rocket engines. In the present approach, turbulence is represented by the extended $k-{\varepsilon}$ turbulence model. To account for the real fluid effects, the propellant mixture properties are calculated by using SRK (Souve-Redlich-Kwong) equation of state. In order to realistically represent the turbulence-chemistry interaction in the turbulent nonpremixed flames, the flamelet approach based on the real fluid flamelet library has been adopted. Based on numerical results, the detailed discussions are made for the real fluid effects and the precise structure of gaseous methane/liquid oxygen coaxial jet flame.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.4
• /
• pp.1-8
• /
• 1999

Numerical study is carried out to investigate flow characteristics around an axisymmetric body with and without an integrated propulsor. The incompressible Reynolds-Averaged Navier-Stokes(RANS) equations are also solved using the finite volume method and the standard $k-\varepsilon$ turbulence model for turbulence closure. In order to investigate the propulsor-hull interaction, the induced velocity calculated by surface panel methods is utilized for the boundary condition at the propeller plane. The calculated results are compared to the experimental results. It is considered that the present numerical code can be used for design of an integrated propulsor.

• Composites Research
• /
• v.13 no.2
• /
• pp.30-39
• /
• 2000

It is difficult for civil engineer to apply composite materials of laminated type to structure for civil construction because of complicated theory although those have much merit. A simple method by which one can predict "exact" values of the natural frequency of vibration of laminated plates is presented. Many laminates with particular orientations have negligible values of $B_{16}$ and $B_{26}$ as the number of plies increases. Such laminates, with $D_{16}$ = $D_{26}$ ->0 can be solved by the same equation as for the special orthotropic laminates(1991,1995). If the quasi-isotrpic constants are used, It is possible to simplify analysis procedure since the equations for isotropic plates can be used. Use of some coefficients can produce "exact" values for laminates with such configuration. This coefficient, in fact, represents the effect of the anisotropy of the laminate. D. H. Kim proposed to use a correction factor, he developed, to produce "exact" solution out of the approximate solution obtained by using the quasi-isotropic constants(1995). In this paper, the fiber orienation studied is [$\alpha$/$\beta$/$\beta$/$\gamma$/$\alpha$/$\alpha$/$\beta$]r, with $\alpha$=-$\beta$, and $\gamma$${\gamma}$=$0^{\circ}$ or $90^{\circ}$ and the above-mentioned method is used to obtain the natural frequencies of plates with such orientations, and detailed result is presented for the use of engineers.nted for the use of engineers.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.18 no.5
• /
• pp.35-45
• /
• 1981

An asymptotic solution of electromagnetic waves scattered by a right-angled dielectric wedge for plane wave incidence is obtained. Scattered fields are constructed by waves reflected and refracted from dielectric interfaces (geometric-optical fields) and a cylindrical wave diffracted from the edge. The edge diffracted field is obtained by adding a correction to the edge diffraction of physical optics approximation, where the correction field is calculated by solving a dual series equation amenable to simple numerical calculation. Validity of this result is assured by two limits of relative dielectric constant $\varepsilon$ of the wedge. The total asymptotic field calculated results in a Rawlins' Neumann series solution for small $\varepsilon$, and the edge diffraction pattern is shown to approach that of a perfectly conducting wedge for large $\varepsilon$. Calculated field patterns are presented and the accuracy of physical optics approximation is discussed.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.333-340
• /
• 1999

A numerical study of turbulent condensing vapor jet submerged in subcooled liquids has been conducted. A physical model of the process is presented employing the locally homogeneous flow approximation of two phase flow in conjunction with a $\kappa$-$\varepsilon$-g model of turbulence properties. In this model the turbulence is represented by differential equations for its kinetic energy and dissipation. A differential equation for the concentration fluctuations is solved and a clipped normal probability distribution function is proposed for the mixture fraction. Effects of steam mass flux, pool temperature and nozzle internal diameter on the condensing vapor jet are also analyzed. The model is evaluated using existing data for turbulent condensing vapor jets. The agreement between the predictions and the available experimental data is good.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.541-550
• /
• 1997

Computational and experimental investigations on the three-dimensional flowfield through an automotive cooling fan are carried out in this work. Steady, incompressible, three-dimensional, turbulent flow through a rotating axial-flow fan is analyzed with Reynolds averaged Navier-Stokes equations and standard k-.epsilon. turbulence model. The governing equations are discretized with finite-volume approximations in non-orthogonal curvilinear coordinates. Computational static pressures on the casing wall agree well with the experimental data which are measured in this work. And, they are sensitive to the change of tip clearance. The flowfield is not significantly affected by the thickness of the blade. The k-.omega. model gives the static pressure rise on the casing wall which is similar to that with the k-.epsilon. model.

• Journal of Korea Water Resources Association
• /
• v.33 no.5
• /
• pp.581-592
• /
• 2000

Turbulence structure and suspended sediment transport capacity in vegetated open-channel flows are investigated numerically in the present paper. The $\textsc{k}-\;\varepsilon$ model is employed for the turbulence closure. Mean velocity and turbulence characteristics including turbulence intensity, Reynolds stress, and production and dissipation of turbulence kinetic energy are evaluated and compared with measurement data available in the literature. The numerical results show that mean velocity is diminished due to the drag provided by vegetation, which results in the reduction of turbulence intensity and Reynolds stress. For submerged vegetation, the shear at the top of vegetation dominates turbulence production, and the turbulence production within vegetation is characterized by wakes. For emergent condition, it is observed that the turbulence generation is dominated by wakes within vegetation. In general, simulated profiles compares favorably to measured data. Computed values of eddy viscosity are used to solve the conservation equation for suspended sediment, yielding sediment concentration more uniform over the depth compared with the one in the plain channel. The simulation reveals that the suspended load decreases as the vegetation density increases and the suspended load increases as the particle diameter decreases for the same vegetation density.

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

수제는 유수에 의한 하안 침식을 방지하기 위한 목적의 수공구조물이다. 최근에는 하천복원 및 생태계복원에 대한 관심이 증가되면서, 하안 침식 방지뿐만 아니라 수제에 의한 흐름분리 및 재순환영역의 생태 서식처 제공에 대한 역할이 부각되고 있다. 따라서 수제 설계 시, 수제 설치에 따른 수리특성 변화를 미리 파악하는 것이 필요하지만, 수제 간격, 설치 방향, 수제 높이, 수제 길이, 하도 특성 및 접근 유수 특성 등 다양한 인자에 의해 수리특성 변화를 예측하기란 쉽지 않으며, 비용과 시 공간 측면에서 물리 모형을 이용한 접근법도 용이하지 못하다. 따라서 본 연구에서는 2차원 수치모형을 이용하여 수제 설치에 따른 흐름특성 변화를 모의하고, 모형의 적용성을 판별하고자 하였다. 지배 방정식은 기본적으로 천수방정식을 적용하게 되는데, 수제 주변 흐름 특성은 난류에 의한 횡방향 확산이 중요한 인자로작용하게 된다. 상수 와점성(constant eddy viscosity) 모형, 포물선형 와점성 (parabolic eddy viscosity) 모형을통해 개발 모형의예측 결과와 실험수로 결과를 비교하였다. 수제하류의 재순환영역의 길이가 실험수로의 결과에 비해서 5~6 배까지 확대된 결과를 나타내고 있어, 만족할 만한 결과를 현재까지 얻지 못하였다. 횡방향 확산이 중요한 만큼 k-$\varepsilon$ 모형 등 다양한 난류 모형의 적용이 고려되어질 필요가 있으며, 또한 1차원 요소가 아닌 2차원 요소의 적용을 고려할 필요도 있다.

• Proceedings of KOSOMES biannual meeting
• /
• 2017.11a
• /
• pp.162-163
• /
• 2017

해양 HNS(Hazardous and Noxious Substances)의 유출 사고 시, 막대한 인명 피해와 환경 훼손을 피하기 위해 유출 사고 조기 예측과 정확한 확산 경로를 예측하는 것이 필수적이다. 본 연구의 최종목적은 전산유체역학을 이용하여 HNS사고가 발생하였을 때 위험구역을 적절히 예측할 수 있는 수치해석기법을 개발하고, 다양한 해양사고조건과 환경영향을 고려하여 근접역에서의 2차원 확산 특성을 고찰하고 확산 현상을 예측하기 위한 모델을 개발하는 것이다. 본 연구에서는 상용코드인 ANSYS FLUENT(V. 17.2)을 사용하여 근접역에서의 2차원 확산특성을 모사하고 분석하였다. 특히, 누출된 HNS의 위치별 농도를 예측하기 위해 종수송방정식(Species Transport Equation)을 이용하였으며 RANS(Reynolds-Averaged Navier-Stokes) 방정식과 표준 $k-{\varepsilon}$ 모델을 이용하여 난류유동을 모사하였다. 해석된 결과는 문헌에서 얻어진 실험데이터와 상호비교하였으며 해수의 유속, HNS의 밀도에 따른 유층 두께, 해수면 HNS 평균 농도 그리고 HNS 전파 속도를 분석하였다. 유층 두께는 해류 유속에 따라 변화하며 변화 경향에 따라 두 구간으로 나눌 수 있다. 해류 전파 속도는 대체로 해류 유속과 선형적 비례관계를 갖는 것으로 나타났다. 해수면 평균 HNS 농도는 해류 유속에 선형적으로 비례하여 감소하며, HNS 밀도가 큰 경우 해수면 평균 HNS 체적 농도는 더 빠르게 감소하게 된다. 이러한 결과는 HNS 확산 특성을 분석하고 관련된 예측모델을 개발하는 데에 기여할 수 있다.