• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.575-581
• /
• 1988

To integrate the two-dimensional Navier-Stokes equations numerically in multiply-connected flow regions, the vorticity-stream function formulation is used. The steady stream function value at the surface of the multibody, initially unknown, has been determined interactively by introducing a line integral which requires the single-valuedness of pressure at each interaction step. This procedure is relatively simpler and more efficient than the primitive variable formulation which requires much more computing time and shows poor convergence. Three doubly-connected flow problems are defined and numerically analyzed by the present method. The results have been compared either with earlier existing ones or with the experimental interferograms to demon-strate the validity of the presented method.

• Journal of the Korean earth science society
• /
• v.35 no.5
• /
• pp.333-341
• /
• 2014

A new set of basis functions was constructed using the Rossby-Haurwitz waves, which are the eigenfunctions of nondivergent barotropic vorticity equations on the sphere. The basis functions were designed to be non-separable, that is, not factored into functions of either the longitude or the latitude. Due to this property, the nodal lines of the functions are aligned neither along with the meridian nor the parallel. The basis functions can be categorized into groups of which members have the same degree or the total wavenumber-like index on the sphere. The orthonormality of the basis functions were found to be close to the machine roundoffs, giving the error of $O(10^{-15})$ or $O(10^{-16})$ for double-precision computation (64 bit arithmetic). It was demonstrated through time-stepping procedure that the basis functions were also the eigenfunctions of the non-divergent barotropic vorticity equations. The projection of the basis functions was carried out onto the low-resolution geopotential field of Gaussian bell, and compared with the theory. The same projections were performed for the observed atmospheric-geopotential height field of 500 hPa surface to demonstrate decomposition into the fields that contain disturbance of certain range of horizontal scales. The usefulness of the new basis functions was thus addressed for application to the eigenmode analysis of the atmospheric motions on the global domain.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.640-645
• /
• 2003

In this paper, we report the numerical and experimental solutions of the vortical flows driven by an impingement of fluid from the bottom wall of a circular cylinder. We managed to visualize successfully the flow pattern shown on the vertical plane through the container axis. The numerical results are shown to compare well with the experimental results for the case of infinity Rossby number. The satisfactory agreement between the two results was possible when in the numerics the free surface was treated as a solid wall so that a no-slip condition was applied on the surface. The numerical solutions reveal that inertial oscillation plays an important role at small Rossby numbers, or at a large background rotation.

• Journal of the Korean earth science society
• /
• v.30 no.5
• /
• pp.630-645
• /
• 2009

The tropopause fold event that took place on January 1, 1997 over mid-region on the Korean Peninsula is examined by means of a numerical simulation based on a Mesoscale Model (MM5). The purpose of this paper is to investigate the effects of moisture in reinforcing a tropopause fold linked to an explosive cyclone. Two types of simulations were carried out; 1) simulations for moist conditions in which full physical and dynamic processes are considered and 2) simulations for dry conditions in which cumulus parameterization and cloud microphysics process are excluded. The results of the moist condition simulations demonstrate that the intensity of the central pressure of the cyclone was overestimated compared with the observed values and that the location of the center and the pressure deepening rates (-17 hPa/12 hr) complied with the observed values. The potential vorticity (PV) anomaly on the isentropic surface at 305 K continued to move in a southeast direction on January 1, 1997 and thus created a single tube of tropopause fold covering the northern and the middle area of the Korean Peninsula and reaching the ground surface at 0300 UTC and 0600 UTC. The results of the dry condition simulations show that the tropopause descended to 500 and 670 hPa in 0300 and 0600 UTC, respectively at the same location for the moist condition simulation; however, there was no deep tropopause fold observed. A comparison of the simulated data between the moist and the dry conditions suggests that a deep tropopause fold should happen when there is sufficient moist in the atmosphere and significantly large PV in the lower atmosphere pulls down the upper atmosphere rather than when the tropopause descends itself due to dynamic causes. Thus, it is estimated that moisture in the atmosphere should have played a crucial role in a deep tropopause fold process.

• Atmosphere
• /
• v.32 no.3
• /
• pp.207-221
• /
• 2022

In Jeju, on January 23, 2016, a cold surge accompanied by heavy snowfall with the most significant amount of 12 cm was the highest record in 32 years. During this period, the temperature of 850 hPa in January was the lowest in 2016. Notably, in 2016, the average surface temperature of January on the Polar cap was the highest since 1991, and 500 hPa geopotential height also showed the highest value. With this condition, the polar vortex in the northern hemisphere meandered and expanded into the subtropics regionally, covering the Korean Peninsula with very high potential vorticity up to 7 Potential Vorticity Unit. As a result, the strong cold advection, mostly driven by a northerly wind, around the Korean Peninsula occurred at over 2𝜎. Previous studies have not addressed this extreme synoptic condition linked to polar vortex expansion due to the unprecedented Arctic warming. We suggest that the occurrence of a strong Ural blocking event after the abrupt warming of the Barents/Karas seas is a major cause of unusually strong cold advection. With a specified mesoscale model simulation with SST (Sea Surface Temperature), we also show that the warmer SST condition near the Korean Peninsula contributed to the heavy snowfall event on Jeju Island.

• Journal of the Korean Institute of Gas
• /
• v.25 no.1
• /
• pp.14-19
• /
• 2021

The change of the vorticity and the temperature distribution in heat exchanger tube bank were analyzed for the flows with the constant inlet velocity and the sinusoidal inlet velocity. The flow frequency characteristics were examined by analyzing power spectral density of lift and drag at a typical circular tube in the tube bank. Karman vortex street could be seen at the upstream region of tube bank for the case of constant inlet velocity. It could be seen that the Karman vortex street was affected by the change of inlet velocity near the circular tubes for the case with the sinusoidal inlet velocity. It was observed that the unsteady temperature distributions for both inlet velocity conditions had almost the same motion as the flow vorticity behavior. The flow frequency for the case with the constant inlet velocity is 37.25Hz, and that with the sinusoidal inlet velocity, the flow frequency is 18.63Hz, which is equal to the sinusoidal inlet velocity. The mean surface Nusselt number(Nu) for overall heat exchanger tube bank was 1051 for the case with the constant inlet velocity and 1117 for the case with the sinusoidal inlet velocity. From the result of heat transfer analysis, it could be seen that Nu with the sinusoidal inlet velocity showed 6.3% increase than that with the constant inlet velocity.

• Journal of Ship and Ocean Technology
• /
• v.7 no.3
• /
• pp.34-55
• /
• 2003

Explicit Kutta condition approximation, proved useful in existing time-domain solver of the unsteady propeller problem, requires a specified functional behavior of the vorticity in space near the trailing edge. In this paper, the strength of the discrete vortices is controlled to have a specified behavior in space in the frequency domain approach. A new formulation is introduced and is implemented for analysis of a lifting surface of a rectangular planform. Sample computations carried out according to the new formulation compares well with that of existing unsteady lifting problem in the time domain.

• Korean Journal of Computational Design and Engineering
• /
• v.1 no.1
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1377-1382
• /
• 2004

Flow fields in a half ducted propeller fan have been investigated by three-dimensional Reynolds-averaged Navier-Stokes (RANS) simulations and a vortex core identification technique. The simulation at the design operating condition shows that the tip vortex onset point is located at 30 percent tip chord of the suction surface on the blade tip. There is no interaction between the tip vortex and the adjacent blade, so that the tip vortex smoothly convects to the rotor exit. However, the high vorticity in the tip vortex causes the wake and the tip leakage flow to be twined around the tip vortex and to interact with the pressure surface of the adjacent blade. This flow behavior corresponds well with experimental results by Laser Doppler Velocimetry. On the contrary, the simulation at the low-flowrate operating condition shows that the tip vortex onset point is located at the 60 percent tip chord of the suction surface. In contrast to the design operating condition, the tip vortex grows almost tangential direction, and impinges directly on the pressure surface of the adjacent blade.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2009.10a
• /
• pp.243-244
• /
• 2009

From large to small vessels of the surface is not smooth. and The surface ship has a surface roughness. Because surface roughness increases the surface resistance and heat transfer, be considered when designing a ship that is an important design factor. Due to surface roughness study on flow around and due to changes in flow and turbulence intensity for the ongoing research is conducted. Roughness of the surface ships from the ship by air as well as machines can be widely applied. In this study, the surface roughness of the leisure marine interval, any change will affect the surface flow, area due to surface roughness for boundary-experimental study was performed.