• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.52-57
• /
• 2001

Laminar two-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for the low Reynolds number (Re=164∼280). The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. The convection term is applied by the 7th order up wind scheme and the pressure and viscosity terms are applied by the 4th order central difference. The grid system makes use of the regular grid system and it is generated by an equation. The calculated results of drag coefficients, lift coefficients, pressure distributions, and vorticity contours and other information are compared with experimental and numerical ones. These results obtained by the present DNS show good agreement with the previous studies.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Journal of Computational Design and Engineering
• /
• v.3 no.4
• /
• pp.330-336
• /
• 2016

The present work is concerned with heat and mass transfer of an electrically conducting second grade MHD fluid past a semi-infinite stretching sheet with convective surface heat flux. The analysis accounts for thermophoresis and thermal radiation. A similarity transformations is used to reduce the governing equations into a dimensionless form. The local similarity equations are derived and solved using Nachtsheim-Swigert shooting iteration technique together with Runge-Kutta sixth order integration scheme. Results for various flow characteristics are presented through graphs and tables delineating the effect of various parameters characterizing the flow. Our analysis explores that the rate of heat transfer enhances with increasing the values of the surface convection parameter. Also the fluid velocity and temperature in the boundary layer region rise significantly for increasing the values of thermal radiation parameter.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.522-522
• /
• 1996

An investigation of the effects of transverse magnetic field and Peltier effect on melt convection and macrosegregation in vertical Bridgman crystal grosth of Te doped InSb was been carried out by means of microstructure observation, Hall measurement, electrical resistivity measurement and X-ray analysis. Before the experiments, Interface stability, convective instability and suppression of convection by magnetic field were calculated theoretically. After doping 1018, 1019 cm-3 Te in InSb, the temperature of Bridgman furnace was set up at $650^{\circ}C$. The samples were grown in I.D. 11mm, 100mm high quartz tube. The velocity of growth was about 2${\mu}{\textrm}{m}$/sec. In order to obtain the suppression of convection by magnetic field in the middle of growth, 2-4KG magnetic field was set on the melt. For searching of the shape of solid-liquid interface and the actual velocity of crystal growth, let 2A current flow from solid to liquid for 1second every 50seconds repeatedly (Peltier effect). The grown InSb was polycrystal, and each grain was very sharp. There was no much difference between the sample with and without magnetic field at a point of view of microstructure. For the sample with Peltier effect, the Peltier marks(striation) were observed regularly as expected. Through these marks, it was found that the solid-liquid interface was flat and the actual growth velocity was about 1-2${\mu}{\textrm}{m}$/sec. On the ground of theoretical calculation, there is thermosolutal convection in the Te doped InSb melt without magnetic field in this growth condition. and if there is more than 1KG magnetic field, the convection is suppressed. Through this experiments, the effective distribution coefficients, koff, were 0.35 in the case of no magnetic field, and 0.45 when the magnetic field is 2KG, 0.7 at 4KG. It was found that the more magnetic field was applied, the more convection was suppressed. But there was some difference between the theoretical calculation and the experiment, the cause of the difference was thought due to the use of some approximated values in theoretical calculation. In addition to these results, the sample with Peltier effect showed unexpected result about the Te distribution in InSb. It looked like no convection and no macrosegregation. It was thought that the unexpected behavior was due to Peltier mark. that is, when the strong current flew the growing sample, the mark was formed by catching Te. As a result of the phenomena, the more Te containing thin layer was made. The layer ruled the Hall measurement. The values of resistivity and mobility of these samples were just a little than those of other reference. It was thought that the reason of this result was that these samples were due to polycrystal, that is, grain boundaries had an influence on this result.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.179-182
• /
• 2008

This study simulated the flood inundations of the Nakdong River catchment running through Yangsan, a small city located in the south eastern area of Korea by using the depth averaged two-dimensional hydrodynamic numerical model. The numerical model employs the staggered grid system including moving boundary and a finite different method to solve the Saint-Venant equations. A second order upwind scheme is used to discretize the nonlinear convection terms of the momentum equations, whereas linear terms are discretized by a first order leap-frog scheme(Cho and Yoon, 1998). The numerical model was applied to a real topography to simulate the flood inundation of the Yangsan basin. The numerical results for urban district are visualized in three dimension. These results can be essentially utilized to construct the three dimensional inundation map after building the GIS-based database in local public organizations in order to protect the life and property safely.

• Wind and Structures
• /
• v.17 no.3
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1035-1054
• /
• 2010

We consider a mixed type singularly perturbed one dimensional elliptic problem with discontinuous source term. The domain under consideration is partitioned into two subdomains. A convection-diffusion and a reaction-diffusion type equations are posed on the first and second subdomains respectively. Two hybrid difference schemes on Shishkin mesh are constructed and we prove that the schemes are almost second order convergence in the maximum norm independent of the diffusion parameter. Error bounds for the numerical solution and its numerical derivative are established. Numerical results are presented which support the theoretical results.

• Atmosphere
• /
• v.30 no.4
• /
• pp.335-346
• /
• 2020

Three free parameters included in a cumulus parameterization are optimized by using micro-genetic algorithm for three precipitation cases occurred in the Korea Peninsula during the summer season in order to reduce biases in a regional model associated with the uncertainties of the parameters and thus to improve the predictability of precipitation. The first parameter is the one that determines the threshold in convective trigger condition. The second parameter is the one that determines boundary layer forcing in convective closure. Finally, the third parameter is the one used in calculating conversion parameter determining the fraction of condensate converted to convective precipitation. Optimized parameters reduce the occurrence of convections by suppressing the trigger of convection. The reduced convection occurrence decreases light precipitation but increases heavy precipitation. The sensitivity experiments are conducted to examine the effects of the optimized parameters on the predictability of precipitation. The predictability of precipitation is the best when the three optimized parameters are applied to the parameterization at the same time. The first parameter most dominantly affects the predictability of precipitation. Short-range forecasts for July 2018 are also conducted to statistically assess the precipitation predictability. It is found that the predictability of precipitation is consistently improved with the optimized parameters.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.165-173
• /
• 1998

Analytical solutions to predict the movement rate of submerged mound are derived using the convection coefficient and the joint distribution function of wave heights and periods. Assuming that the sediment is moved onshore due to the velocity asymmetry of Stokes' second order nonlinear wave theory, the micro-scale bedload transport equation is applied to the sediment conservation. The nonlinear convection-diffusion equation can then be obtained which governs the migration of submerged mound. The movement rate decreases exponentially with increasing the water depth, but the movement rate tends to increase as the spectral width parameter, $ u$ increases. In comparison of the analytical solution with the measured data, it is found that the analytical solution overestimates the movement rate. However, the agreement between the analytical solution and the measured data is encouraging since this over-estimation may be due to the inaccuracy of input data and the limitation of sediment transport model. In particular, the movement rates with respect to the water depth predicted by the analytical solution are in very good agreement with the estimated result using the discritization technique with the hindcast wave data.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.