• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.169-169
• /
• 2015

유사 입자의 크기는 유사의 특성 및 그에 따른 거동 변화에 중요한 영향을 미친다. Stokes 침강 속도 모형에서 유사의 침강 속도에 가장 많은 영향을 주는 인자는 유사의 크기인 것 또한 확인된다. 유사 입자의 크기가 약 $60{\mu}m$보다 작은 유사들은 알갱이 사이의 점착력을 무시할 수 없다. 이로 인해 유사들은 응집 현상을 겪으며 입자 본래의 크기보다 크기가 큰 플럭을 형성하는 점착성 유사로 분류된다. 응집 현상이란, 흐름 내 점착성을 띠는 일차입자(Primary Particle)가 응집과 파괴를 반복하며 플럭을 형성하는 현상을 뜻한다. 입자 간의 충돌을 통해 응집이 진행되며 난류 전단으로 인해 형성된 플럭의 파괴가 발생한다. 많은 연구에서 점착성 유사의 충돌을 야기하는 가장 지배적인 원리는 난류라 알려져 있다. 이러한 응집 현상으로 인하여 플럭의 크기와 밀도는 지속적으로 변화를 겪으며 비점착성 유사와 다른 특징들을 보인다. 흐름에 존재하는 유사의 이동은 이송-확산 방정식을 통해 표현된다. 이송-확산 방정식은 시간 변화에 따른 농도의 변화를 입자의 침강과 난류 및 유사 자체의 특징에 의한 확산으로 해석한다. 침강속도로 대변되는 이송과 달리, 확산은 난류흐름 내에서 유사가 확산되는 정도를 정량화하기 위한 인자가 요구된다. 난류에 의한 유사의 확산은 유사 자체 특성에 따른 물질 확산에 비하여 매우 큰 값을 가지며, 이를 확산 계수로 개념화 한다. 확산계수는 와점성계수와 Schmidt 수(${\sigma}_c$)의 비로 정의된다. ${\sigma}_c$는 난류의 점성과 난류로 인한 부유과정에 의해 유사가 확산되는 정도를 나타낸다. 이에 따라 ${\sigma}_c$의 변화가 유사의 부유 및 침강거동에 많은 영향을 미칠 것이라 판단되나, 국내외에서 수행된 연구 동향에서는 ${\sigma}_c$를 0.5부터 1.0 사이의 상수를 적용하여 수행되었다. 이에 본 연구에서는 ${\sigma}_c$의 크기에 따라 달라지는 유사의 부유 및 침강 변화에 의한 총 부유량을 살펴보고자 한다. 유사의 점착성을 고려할 수 있는 1DV 수치 모형을 이용하여 비점착성 유사와 점착성 유사를 대상으로 수치연구를 수행하며, 유사의 크기 및 ${\sigma}_c$의 변화에 따른 총 부유량 경향을 살펴본다. 그 결과, 점착성 유사는 ${\sigma}_c$의 증가에 따라서 유사의 총 부유량이 증가하는 현상이 나타난 반면 비점착성 유사는 ${\sigma}_c$의 증가에 따라 유사의 총 부유량이 감소하는 경향이 나타났다. 그러나 크기가 아주 작은 비점착성 유사를 대상으로 수치 연구를 수행한 결과, ??에 따른 총 부유량의 경향은 유사의 점착성에서 기인하는 것이 아닌 입자의 크기로부터 야기되는 특성이라는 결론이 도출되었다.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.100-106
• /
• 1997

Viscous flow around an actual ship is calculated by an use of RANS(Reynolds-averaged Navier-Stokes) solver. Reynolds stress is modelled by using k-$\varepsilon$ turbulence model and the law of wall is applied near the body. Body fitted coordinates are introduced for the treatment of the complex boundary of the ship hull form. The transformed equations in the computational domain are numerically solved by an employment of FVM(Finite Volume Method). SIMPLE(Semi-Implcit Pressure Linked Equation) method is adopted in the calculation of pressure and the solution of the disssssssscretized equation is obtained by the line-by-line method with the use of TDMA(Tri-Diagonal Matrix Algorithme). The subject ship model of actual calculation is 4,410 TEU class container carrier. For 4 geosim models the calculated viscous resistancce values are compared with the model test results and analyzed on their componentss. The resistance performance of an actual ship is predicted very resonably, so this mothod may be utilized as a design tool of hull form.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.1
• /
• pp.1-8
• /
• 2000

A modified low-Reynolds-number Reynolds stress model is developed for the calculation of drag-reducing turbulent flows induced by polymer injection. The results without polymer injection are compared with the results of direct numerical simulation to ensure the validity of the basic model. In case of drag reduction, profiles of mean velocity and Reynolds stress components, in two-dimensional channel flow, obtained with a proper value of viscosity ratio are presented and discussed. Computed mean velocity profile is in very good agreement with experimental data. And, the qualitative behavior of Reynolds stress components with the viscosity ratio is also reasonable.

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

Injection characteristics of a kerosene swirl injector of liquid rocket engine operating at supercritical environment have been investigated. Kerosene surrogate models are proposed to model the kerosene properties. Turbulent numerical model is based on large eddy simulation and contains Soave modification of Redlich-Kwong equation of state and Chung's model. Numerical analysis results at supercritical environment are compared with the one at transcritical condition. Differences of density and viscosity are analyzed at both liquid film and core gas in the swirl injector.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.733-739
• /
• 1999

A modified low-Re $k-\varepsilon$ model is used for the calculation of drag-reducing turbulent flow by polymer injection in a pipe. With the viscoelastic model, molecular viscosity in the definition of turbulent viscosity is related to elongations viscosity of the solution to account for the effects of drag reduction. Finite volume method is used for the discretization, and power-law scheme is used as a numerical scheme. Computed dimensionless velocity profiles are in good agreements with the experimental data in case of low drag reductions. However, in case of high drag reductions, they deviate largely from the measurements in the central zone of the flow field.

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

The multiphase simulation of a liquid jet in a lab-scale ramjet combustor with a plain orifice type injector was studied with a commercial CFD tool, a FLUENT program. The objectives of the current study are to analysis the breakup characteristics of a hexane liquid jet in a cross flow and to derive the correlation between flow conditions and drag force coefficients in a test section. From the result of a numerical simulation, we concluded that a DPM and Realizable $k-{\varepsilon}$ model with an enhanced wall treatment were available to simulate the multiphase flow simulation. And the calculated distribution of a hexane vapor concentration was well-matched with experimental results.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.103-108
• /
• 2011

The turbulent mixing of a kerosene/liquid oxygen coaxial swirl injector under supercritical pressures have been numerically investigated. Kerosene surrogate models are proposed for the kerosene thermodynamic properties. Turbulent numerical model is based on LES(Large Eddy Simulation) with real-fluid transport and thermodynamics over the entire pressure range; Soave modification of Redlich-Kwong equation of state, Chung's model for viscosity/conductivity, and Fuller's theorem for diffusivity to take account Takahashi's compressible effect. The effect of operating pressure on thermodynamic properties and mixing dynamics inside an injector and a combustion chamber are investigated. Power spectral densities of pressure fluctuations in the injector under various chamber pressure are analyzed.

• Journal of Energy Engineering
• /
• v.2 no.3
• /
• pp.268-275
• /
• 1993

A series of parametric investigations are performed in order to resolve the flow characteristic of a Quebec city stoker incinerator. The parameters considered in this study are five internal configurations of the Quebec city stoker itself and its modified ones, primary air velocity, the injection velocity and angle of the secondary air, and the reduction of the stoker exit area. A control-volume based finite-difference method by Patankar together with the power-law scheme is employed for discretization. The resolution of the pressure-velocity coupling is made by the use of SIMPLEC algorithm. The standard, two equation, k-$\varepsilon$ model is incorporated for the closure of turbulence. The size of recirculation region, turbulent viscosity, the mass fraction of the secondary air and pressure drop are calculated in order to analyze the characteristics of flow field. The results are physically acceptable and discussed in detail. The flow field of the Quebec city stoker shows the strong recirculation zone together with the high turbulence intensity over the upper part of the incinerator.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.5
• /
• pp.1-9
• /
• 2017

The cryogenic nitrogen spray of a swirl injector has been numerically investigated using three dimensional LES turbulence model to analyze the dynamic characteristics under supercritical condition. To predict the precise nitrogen properties under supercritical condition, SRK equation of state, Chung's method for viscosity and thermal conductivity and Takahashi's correlation based on Fuller's theory for diffusion coefficient are implemented. The complex flow structures due to interaction between flow field and acoustic field are observed inside and outside the injector under supercritical condition. FFT, POD, and DMD techniques are employed to understand the coherent structures. By implementing the FFT, the dominant frequencies are identified inside and outside the injector. The coherent flow structures related to the dominant frequencies are visualized using the POD and DMD techniques. In addition, the DMD provides the damping coefficient which is related with the instability prediction.