• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.13 no.3
• /
• pp.145-151
• /
• 2003

Using computational fluid dynamics (CFD) method, we investigated three-dimensional fluid flow field and particle movement with respect to the injected gas flow rate variation in typical cyclone separator system. The results of numerical investigation were deduced by coupling the analysis of fluid flow field with Wavier-stokes equation and the tracking of the particle trajectory with Langrangian approach. It was shown that the increasing of injected gas flow rate resulted in the increasing of pressure loss in the separator. This change of inner pressure had an effect on an aspect of the fluid flow in the separator. Particle movement was determined by fluid flow in the separator and was fully depended on a diameter of particles under the fixed flow rate. Increasing of injected gas flow rate was led to an increasing of the trace of particle, so the particles moved to the lower part of the separator. For this reason, the minimum diameters of the particles were decreased and increased the separation rate under the fixed particle diameter. In conclusion, the changes of injected gas flow rate have an important factor to the fluctuation of the fluid flow field and particle trajectory in the separator.

• 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.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.6
• /
• pp.524-532
• /
• 2020

In this study, noise radiated from a high-speed fan-motor unit for a cordless vacuum cleaner is reduced by designing splitter blades on the existing impeller. First of all, in order to investigate the flow field through a fan-motor unit, especially impeller, the unsteady incompressible Reynolds-Averaged Navier-Stokes (RANS) equations are numerically solved by using computational fluid dynamic technique. With predicted flow field results as input, the Ffowcs Williams-Hawkings (FW-H) integral equation is solved to predict aerodynamic noise radiated from the impeller. The validity of the numerical methods is confirmed by comparing the predicted sound pressure spectrum with the measured one. Further analysis of the predicted flow field shows that the strong vortex is formed between the impeller blades. As the vortex induces the loss of the flow field and acts as an aerodynamic noise source, supplementary splitter blades are designed to the existing impeller to suppress the identified vortex. The length and position of splitter are selected as design factors and the effect of each design factor on aerodynamic noise is numerically analyzed by using the Taguchi method. From this results, the optimum location and length of splitter for minimum radiated noise is determined. The finally selected design shows lower noise than the existing one.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.11
• /
• pp.897-903
• /
• 2015

In this paper, the contrast therapy of skin was numerically investigated by solving the conjugate heat transfer problem. A finite volume method based on the SIMPLE algorithm was adopted to solve the axisymmetric incompressible Navier-Stokes equations, coupled with an energy equation. These equations are strongly coupled with the Pennes bio-heat equation in order to consider the effect of blood perfusion rate. We investigated the thermal response of skin at some selected depths for various input temperature profiles of a stimulator for contrast therapy. From the numerical simulations, the regions with cold/hot threshold temperatures were found for five input temperature profiles. It was shown that the temperature varies mildly for different input profiles as the depth increases, owing to the Pennes effect. The input temperatures for effective hot/cold stimulation of dermis layer were found to be $47^{\circ}C$ and $7^{\circ}C$, respectively. The present numerical results will be used for finding an optimal temperature profile of a stimulator for contrast therapy.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.224-232
• /
• 1999

Latent heat storage system using micro-encapsuled phase change material is effective method for floor heating of house and building. The temperature profile in capsule block and flow rate of hot water are important parameters for the development of heat storage system. In the present study, a mathematical model based on 3-D, non-steady state, Navier-Stokes equations, scalar conservation equations and turbulence model ($\kappa$-$\varepsilon$), is used to predict the temperature profiles in capsule and the velocity vectors in hot water pipe. The multi-block grids and fine grids embedding are used to join the circle in hot water pipe and square in capsule block. The phase change process of the capsule is quite complex not only because the size of phase change material is very small, but also because phase change material is mixed with the cement to form thermal storage block. In calculation, it's assumed that the phenomena of phase change is limited only the thermal properties of phase change material and the change of boundary is not happened in capsule. The purpose of this study is to calculate the temperature profiles in capsule block and velocity vectors in hot water pipe using the numerical calculation. Two kinds of thermal boundary condition were considered, the first (case 1) is the adiabatic condition for the both outside surfaces of the wall, the second (case 2) is the case in which one surface is natural convection with atmosphere and another surface is adaibatic. Calculation results are shown that the temperature profile in capsule block for case 1 is higher than that for case 2 due to less heat loss in adaibatic surface. Specially, in the domain of near Y=0, the difference of temperature is greater in case 1 than in case 2. The detailed experimental data of capsule block on the temperature profile and the thermal properties such as specific heat and coefficient of heat transfer with the various temperature are required to predict more exact phenomena of heat transfer.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.9
• /
• pp.597-604
• /
• 2016

A parallel algorithm of bi-conjugate gradient method was developed based on CUDA for parallel computation of the incompressible Navier-Stokes equations. The governing equations were discretized using splitting P2P1 finite element method. Asymmetric stenotic flow problem was solved to validate the proposed algorithm, and then the parallel performance of the GPU was examined by measuring the elapsed times. Further, the GPU performance for sparse matrix-vector multiplication was also investigated with a matrix of fluid-structure interaction problem. A kernel was generated to simultaneously compute the inner product of each row of sparse matrix and a vector. In addition, the kernel was optimized to improve the performance by using both parallel reduction and memory coalescing. In the kernel construction, the effect of warp on the parallel performance of the present CUDA was also examined. The present GPU computation was more than 7 times faster than the single CPU by double precision.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.3
• /
• pp.189-197
• /
• 2006

The influence of leodo Ocean Research Station structure to surrounding atmospheric flow is carefully investigated using CFD techniques. Moreover, the validation works of computational results are performed by the comparison with the observed data of leodo Ocean Research station. In this paper, we performed 3-dimensional CAD modelling of the station, generated the grid system for numerical analysis and carried out flow analyses using Navier-Stokes equations coupled with two-equation turbulence model. For suitable free stream conditions of wind speed and direction, the interference of the research station structure on the flow field is predicted. Beside, the computational results are benchmarked by observed data to confirm the accuracy of measured date and reliable data range of each measuring position according to the wind direction. Through the results of this research, now the quantitative evaluation of the error range of interfered gauge data is possible, which is expected to be applied to provide base data of accurate sea surface wind around research stations.

• Journal of Navigation and Port Research
• /
• v.40 no.4
• /
• pp.173-181
• /
• 2016

From 2016, controls on reduction of NOx and SOx emissions from the vessels that are operated in the emission control area were tightened. The selectivity catalytic reduction system of the denitrification equipment which NOx among the above controlled materials is very effective and used commercially very much. But it has the disadvantage that CSR is activated at high temperatures. Therefore, the SCR and SCR activation instrument that can react even at low temperatures by using micro-nano bubbles so that the above problems can be minimized were developed. And the computational fluid dynamics technique was used by ANSYS-CFX package to prepare the plan that improves the SCR system's efficiency. Simulation for the viscous flow analysis of the SCR system was executed by applying the Navier-Stokes equation to it as a governing equation. For the SCR system's shape, 3D modeling was done by using CATIA V5. SCR jet nozzle's position was checked by changing it to the intervals of 1/3, 1/2, and 2/3 from the inlet of the vent pipe to compare the SCR system's efficiency. And the number of nozzles was compared and analyzed by simulating 4, 6, and 8 holes to check an effect of the number on the SCR system's efficiency. The simulation result has found that the closer nozzles are to the inlet of the vent pipe and the more nozzles are, the more efficiency is improved.