• Journal of Environmental Science International
• /
• v.14 no.1
• /
• pp.41-51
• /
• 2005

Using three-dimensional non-hydrostatical numerical model with one way double nesting technique, atmo­spheric circulation in the mountainous coastal region in summer was investigated from August 13 through 15, 1995. During the day, synoptic westerly wind blows over Mt. Mishrung in the west of a coastal city, Sokcho toward the East Sea, while simultaneously, easterly upslope wind combined with both valley wind from plain (coast) toward mountain and sea-breeze from sea toward inland coast blows toward the top of the mountain. Two different directional wind systems confront each other in the mid of eastern slope of the mountain and the upslope wind goes up to the height over 2 km, becoming an easterly return flow in the upper level over the sea and making sea-breeze front with two kinds of sea-breeze circulations of a small one in the coast and a large one in the open sea. Convective boundary layer is developed with a thickness of about 1km over the ground in the upwind side of the mountain in the west and a thickness of thermal internal boundary layer from the coast along the eastern slope of the mountain is only confined to less than 200 m. On the other hand, after sunset, no prohibition of upslope wind generated during the day and downward wind combined with mountain wind from mountain towardplain and land-breeze from land toward under nocturnal radiative cooling of the ground surfaces should intensify westerly downslope wind, resulting in the formation of wind storm. As the wind storm moving down along the eastern slop causes the development of internal gravity waves with hydraulic jump motion in the coast, bounding up toward the upper level of the coastal sea, atmospheric circulation with both onshore and offshore winds like sea-breeze circulation forms in the coastal sea within 70 km until midnight and after that, westerly wind prevails in the coast and open seas.

• Journal of Environmental Science International
• /
• v.11 no.9
• /
• pp.865-880
• /
• 2002

As prevailing synoptic scale westerly wind blowing over high steep Mt. Taegulyang in the west of Kangnung coastal city toward the Sea of Japan became downslope wind and easterly upslope wind combined with both valley wind and sea breeze(valley-sea breeze) also blew from the sea toward the top of the mountain, two different kinds of wind regimes confronted each other in the mid of eastern slope of the mountain and further downward motion of downlsope wind along the eastern slope of the mountain should be prohibited by the upslope wind. Then, the upslope wind away from the eastern slope of the mountain went up to 1700m height over the ground, becoming an easterly return flow in the upper level of the sea. Two kinds of circulations were detected with a small one in the coastal sea and a large one from the coast toward the open sea. Convective boundary layer was developed with a thickness of about 1km over the ground in the upwind side of the mountain in the west, while a thickness of thermal internal boundary layer(TIBL) form the coast along the eastern slope of the mountain was only confined to less than 200m. After sunset, under no prohibition of upslope wind, westerly downslope wind blew from the top of the mountain toward the coastal basin and the downslope wind should be intensified by both mountain wind and land breeze(mountain-land breeze) induced by nighttime radiative cooling of the ground surfaces, resulting in the formation of downslope wind storm. The wind storm caused the development of internal gravity waves with hydraulic jump motion bounding up toward the upper level of the sea in the coastal plain and relatively moderate wind on the sea.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• Journal of Environmental Science International
• /
• v.13 no.8
• /
• pp.721-731
• /
• 2004

The dispersion of suspended particulates in the coastal complex terrain of mountain-inland basin (city)-sea, considering their recycling was investigated using three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). Convective boundary layer under synoptic scale westerly wind is developed with a thickness of about I km over the ground in the west of the mountain, while a thickness of thermal internal boundary layer (TIBL) is only confined to less than 200m along the eastern slope of the mountain, below an easterly sea breeze circulation. At the mid of the eastern slop of the mountain, westerly wind confronts easterly sea breeze, which goes to the height of 1700 m above sea level and is finally eastward return flow toward the sea. At this time, particulates floated from the ground surface of the city to the top of TIBL go along the eastern slope of the mountain in the passage of sea breeze, being away the TIBL and reach near the top of the mountain. Then those particulates disperse eastward below the height of sea-breeze circulation and widely spread out over the coastal sea. Total suspended particulate concentration near the ground surface of the city is very low. On the other hand, nighttime radiative cooling produces a shallow nocturnal surface inversion layer (NSIL) of 200 m thickness over the inland surface, but relatively thin thickness less than 100m is found near the mountain surface. As synoptic scale westerly wind should be intensified under the association of mountain wind along the eastern slope of mountain to inland plain and further combine with land-breeze from inland plain toward sea, resulting in strong wind as internal gravity waves with a hydraulic jump motion bounding up to about 1km upper level in the atmosphere in the west of the city and becoming a eastward return flow. Simultaneously, wind near the eastern coastal side of the city was moderate. Since the downward strong wind penetrated into the city, the particulate matters floated near the top of the mountain in the day also moved down along the eastern slope of the mountain, reaching the. downtown and merging in the ground surface inside the NSIL with a maximum ground level concentration of total suspended particulates (TSP) at 0300 LST. Some of them were bounded up from the ground surface to the 1km upper level and the others were forward to the coastal sea surface, showing their dispersions from the coastal NSIL toward the propagation area of internal gravity waves. On the next day at 0600 LST and 0900 LST, the dispersed particulates into the coastal sea could return to the coastal inland area under the influence of sea breeze and the recycled particulates combine with emitted ones from the ground surface, resulting in relatively high TSP concentration. Later, they float again up to the thermal internal boundary layer, following sea breeze circulation.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.105-114
• /
• 1993

The ciculation due to bottom density current produced by a dense inflow into a small reservoir is analysed by numerical scheme. Before the front of the density current arrives at the downstream end, the mixing in the reservoir is mainly caused by the anticlockwise vortex formed at the downstream of plunging point along the movement of bottom density current. Upon the arrival of the front of the density current at the downstream end an internal surge is created through an internal hydraulic jump. With repeated propagation of the internal surge back and forth the mixing in the reservoir is progressed and the thickness of dense layer is increased upward. The dilution of the overflow at downstream end is found to depend on inflow densimetric Froude number, reservoir length and elapsed time. The time required for the overflow to attain a specified dilution increases as reservoir length increases and Fre decreases.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.19-26
• /
• 2003

The dispersion of recycled particulates in the complex coastal terrain containing Kangnung city, Korea was investigated using a three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). The results show that particulates at the surface of the city that float to the top of thermal internal boundary layer (TIBL) are then transported along the eastern slope of the mountains with the passage of sea breeze and nearly reach the top of the mountains. Those particulates then disperse eastward at this upper level over the coastal sea and finally spread out over the open sea. Total suspended particulate (TSP) concentration near the surface of Kangnung city is very low. At night, synoptic scale westerly winds intensify due to the combined effect of the synoptic scale wind and land breeze descending the eastern slope of the mountains toward the coast and further seaward. This increase in speed causes development of internal gravity waves and a hydraulic jump up to a height of about 1km above the surface over the city. Particulate matter near the top of the mountains also descends the eastern slope of the mountains during the day, reaching the central city area and merges near the surface inside the nocturnal surface inversion layer (NSIL) with a maximum ground level concentration of TSP occurring at 0300 LST. Some particulates were dispersed following the propagation area of internal gravity waves and others in the NSIL are transported eastward to the coastal sea surface, aided by the land breeze. The following morning, particulates dispersed over the coastal sea from the previous night, tend to return to the coastal city of Kangnung with the sea breeze, developing a recycling process and combine with emitted surface particulates during the morning. These processes result in much higher TSP concentration. In the late morning, those particulates float to the top of the TIBL by the intrusion of the sea breeze and the ground level TSP concentration in the city subsequently decreases.

• Journal of Korea Water Resources Association
• /
• v.50 no.1
• /
• pp.55-64
• /
• 2017

This study presents a numerical model for simulating turbidity currents using the ULTIMATE scheme. For this, the layer-averaged model is used. The model is applied to laboratory experiments, where the flume is composed of sloping and flat parts, and the characteristics of propagating turbidity currents are investigated. Due to the universal limiter of the ULTIMATE scheme, the frontal part of the turbidity currents at a sharp gradient without numerical oscillations is computed. Simulated turbidity currents propagate super-critically to the end of the flume, and internal hydraulic jumps occur at the break-in-slope after being affected by the downstream boundary. It is found that the hydraulic jumps are computed without numerical oscillations if Courant number is less than 1. In addition, factors that affect propagation velocity of turbidity currents is studied. The particle size less than $9{\mu}m$ does not affect propagation velocity but the buoyancy flux affects clearly. Finally, it is found that the numerical model computes the bed elevation change due to turbidity currents properly. Specifically, a discontinuity in the bed elevation, arisen from the hydraulic jumps and resulting difference in sediment entrainment, is observed.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.24 no.1
• /
• pp.29-45
• /
• 1996

Before (March 26, 1994) or after the occurrence of a downslope windstorm (March 29), the NO, $NO_2$, and $SO_2$ at the ground level of Kangnung city were monitored with high concentrations in the afternoon, due to a large amount of gases emitted from combustion of motor vehicle and heating apparatus, especially near 1600-1800 LST and 2000-2100 LST, but at night, they had low concentrations, resulting from small consumptions of vehicle and heating fuels. When both moderate westerly synoptic-scale winds flow over Mt. Taegwallyang and easterly meso-scale sea breeze during the day, atmospheric pollutants should be trapped by two different wind systems, resulting in higher concentration at Kangnung city in the afternoon. At night, the association of westerly synoptic wind and land breeze can produce relatively strong winds and the dissipation by the winds cause these low concentrations to lower and lower, as nightime goes on. From March 27 through 28, an enforced localized windstorm could be produced along the lee side of the mountain near Kangnung, generating westerly internal gravity waves with hydraulic jump motions. Sea breeze toward inland appartantly confines to the bottom of the eastern side of the mountain, due to the interruption of eastward violent internal gravity waves. As the windstorm moves down toward the ground, an encountering point of two opposite winds approaches Kangnung, and a great amount of NO and $NO_2$ were removed by the strong surface winds. Thus, their maximum concentrations are found to be near 18 and 20 LST, 17 and 21 LST. In the nighttime, the more developed storm should produce very strong surface winds and the NO and $NO_2$ could be easily dissipated into other place. The $SO_2$ concentration had no maximum value, that is, almost constant one all day long, due to its removal by the strong surface winds. Especially, the CO concentrations were slightly lower during the strom period than both before or after the strom, but they were nearly constant without much changes during the during the daytime and nighttime.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.177-184
• /
• 2003

Numerical simulation on marine wind and sea surface elevation was carried out using both three-dimensional hydrostatic and non-hydrostatic models and a simple oceanic model from 0900 LST, August 13 to 0900 LST, August 15, 1995. As daytime easterly meso-scale sea-breeze from the eastern sea penetrates Kangnung city in the center part as basin and goes up along the slope of Mt. Taegullyang in the west, it confronts synoptic-scale westerly wind blowing over the top of the mountain at the mid of the eastern slope and then the resultant wind produces an upper level westerly return flow toward the East Sea. In a narrow band of weak surface wind within 10km of the coastal sea, wind stress is generally small, less than l${\times}$10E-2 Pa and it reaches 2 ${\times}$ 10E-2 Pa to the 35 km. Positive wind stress curl of 15 $\times$ 10E-5Pa $m^{-1}$ still exists in the same band and corresponds to the ascent of 70 em from the sea level. This is due to the generation of northerly wind driven current with a speed of 11 m $S^{-1}$ along the coast under the influence of south-easterly wind and makes an intrusion of warm waters from the southern sea into the northern coast, such as the East Korea Warm Current. On the other hand, even if nighttime downslope windstorm of 14m/s associated with both mountain wind and land-breeze produces the development of internal gravity waves with a hydraulic jump motion of air near the coastal inland surface, the surface wind in the coastal sea is relatively moderate south-westerly wind, resulting in moderate wind stress. Negative wind stress curl in the coast causes the subsidence of the sea surface of 15 em along the coast and south-westerly coastal surface wind drives alongshore south-easterly wind driven current, opposite to the daytime one. Then, it causes the intrusion of cold waters like the North Korea Cold Current in the northern coastal sea into the narrow band of the southern coastal sea. However, the band of positive wind stress curl at the distance of 30km away from the coast toward further offshore area can also cause the uprising of sea waters and the intrusion of warm waters from the southern sea toward the northern sea (northerly wind driven current), resulting in a counter-clockwise wind driven current. These clockwise and counter-clockwise currents much induce the formation of low clouds containing fog and drizzle in the coastal region.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.