• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.761-776
• /
• 2003

Numerical models were employed to investigate the hydrodynamics of water flow in the lake behind a dam and the spillway where supercritical flows and negative pressures are likely to occur. In this study, 2-D model, RMA2 was employed to examine the upstream flow pattern and 3-D CFD model, FLUENT was used to evaluate the three-dimensional flow in the approaching region and flow distributions in the spillways and discharge culverts. The bathymetry and the details of structures were carefully taken into consideration in building the models. The results from applying the 2-D model for the planned Hantan River Dam show that large eddies, the velocity of which reaches up to 1 m/s are occurring in several places upstream of the dam. That means that the 2-D numerical model could be utilized to investigate the two-dimensional flow patterns after the construction of a dam. Three-dimensional numerical results show that the approach flow varies depending on stages and discharge conditions, and velocities at spillways, discharge culverts, and sediment flushing tunnels are differently distributed. The velocity distributions obtained from the numerical model and a hydraulic model at the centerline of spillways 100 m upstream of the dam show reasonably similar results. It is expected that 2-D and 3-D numerical models ate useful tools to help optimize the dam design through investigating the flow patterns in the spillway and at the upstream of the dam, which is not always feasible in hydraulic modeling.

• Journal of the Korean earth science society
• /
• v.39 no.3
• /
• pp.211-227
• /
• 2018

This study evaluates the prediction skills of stratospheric polar vortex intensification events (VIEs) in Global Seasonal Forecasting System (GloSea5) model, an operational subseasonal-to-seasonal (S2S) prediction model of Korea Meteorological Administration (KMA). The results show that the prediction limits of VIEs, diagnosed with anomaly correlation coefficient (ACC) and mean squared skill score (MSSS), are 13.6 days and 18.5 days, respectively. These prediction limits are mainly determined by the eddy error, especially the large-scale eddy phase error from the eddies with the zonal wavenumber 1. This might imply that better prediction skills for VIEs can be obtained by improving the model performance in simulating the phase of planetary scale eddy. The stratospheric prediction skills, on the other hand, tend to not affect the tropospheric prediction skills in the analyzed cases. This result may indicate that stratosphere-troposphere dynamic coupling associated with VIEs might not be well predicted by GloSea5 model. However, it is possible that the coupling process, even if well predicted by the model, cannot be recognized by monotonic analyses, because intrinsic modes in the troposphere often have larger variability compared to the stratospheric impact.

• Fisheries and Aquatic Sciences
• /
• v.13 no.2
• /
• pp.167-181
• /
• 2010

To predict changes in the marine environment of the Beolgyo Stream Estuary in Jeonnam Province, South Korea, where cohesive tidal flats cover a broad area and a large bridge is under construction, this study conducted numerical simulations involving tidal flow and cohesive sediment transport. A wetting and drying (WAD) technique for tidal flats from the Princeton Ocean Model (POM) was applied to a large-scale-grid hydrodynamic module capable of evaluating the flow resistance of structures. Derivation of the eddy viscosity coefficient for wakes created by structures was accomplished through the explicit use of shear velocity and Chezy's average velocity. Furthermore, various field observations, including of tide, tidal flow, suspended sediment concentrations, bottom sediments, and water depth, were performed to verify the model and obtain input data for it. In particular, geologic parameters related to the evaluation of settling velocity and critical shear stresses for erosion and deposition were observed, and numerical tests for the representation of suspended sediment concentrations were performed to determine proper values for the empirical coefficients in the sediment transport module. According to the simulation results, the velocity variation was particularly prominent around the piers in the tidal channel. Erosion occurred mainly along the tidal channels near the piers, where bridge structures reduced the flow cross section, creating strong flow. In contrast, in the rear area of the structure, where the flow was relatively weak due to the formation of eddies, deposition and moderated erosion were predicted. In estuaries and coastal waters, changes in the flow environment caused by artificial structures can produce changes in the sedimentary environment, which in turn can affect the local marine ecosystem. The numerical model proposed in this study will enable systematic prediction of changes to flow and sedimentary environments caused by the construction of artificial structures.

• Ocean and Polar Research
• /
• v.25 no.spc3
• /
• pp.361-371
• /
• 2003

The effect of bottom friction on the steady wind-driven circulation in the Yellow Sea and the East China Sea (YSECS) has been studied using a two-dimensional numerical model with and without tidal forcing. Upwind flow experiment in YSECS has also been carried out with a schematic time variation in the wind field. The surface water setup and circulation pattern due to steady wind forcing are found to be very sensitive to the bottom friction. When the effects of tidal currents are neglected, the overall current velocities are overestimated and eddies of various sizes appear, upwind flow is formed within the deep trough of the Yellow Sea, forming a part of the topographic gyre on the side of Korea. When tidal forcing is taken into account, the wind-induced surface elevations are smoothed out due to the strong tide-induced bottom friction, which is aligned almost normal to the wind stresses; weak upwind flow is farmed in the deep trough of the Yellow Sea, west and south of Jeju. Calculation with wind forcing only through a parameterized linear bottom friction produces almost same results from the calculation with $M_2$ tidal forcing and wind forcing using a quadratic bottom friction, supporting Hunter (1975)'s linearization of bottom friction which includes the effect of tidal current, can be applied to the simulation of wind-driven circulation in YSECS. The results show that steady wind forcing is not a dominant factor to the winter-time upwind flow in YSECS. Upwind flow experiment which considers the relaxation of pressure gradient (Huesh et al. 1986) shows that 1) a downwind flow is dominant over the whole YSECS when the northerly wind reaches a maximum speed; 2) a trend of upwind flow near the trough is found during relaxation when the wind abates; 3) a northward flow dominates over the YSECS after the wind stops. The results also show that the upwind flow in the trough of Yellow Sea is forced by a wind-induced longitudinal surface elevation gradient.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.2
• /
• pp.79-93
• /
• 2015

Three-dimensional structures of large ocean rings in the Gulf Stream region are investigated using the HYbrid Coordinate Ocean Model (HYCOM). Numerically simulated flow structures around four selected cyclonic and anticyclonic rings are compared with two different horizontal resolutions: $1/12^{\circ}$ and $1/48^{\circ}$. The vertical distributions of Lagrangian Coherent Structures (LCSs) are analyzed using Finite Size Lyapunov Exponent (FSLE) and Okubo-Weiss parameters (OW). Curtain-shaped FSLE ridges are found in all four rings with extensions of surface ridges throughout the water columns, indicating that horizontal stirring is dominant over vertical motions. Near the high-resolution rings, many small-scale flow structures with size O(1~10) km are observed while these features are rarely found near the low-resolution rings. These small-scale structures affect the flow pattern around the rings as flow particles move more randomly in the high-resolution models. The dispersion rates are also affected by these small-scale structures as the relative horizontal dispersion coefficients are larger for the high-resolution models. The absolute vertical dispersion rates are, however, lower for the high-resolution models, because the particles tend to move along inclined eddy orbits when the resolution is low and this increases the magnitude of absolute vertical dispersion. Since relative vertical dispersion can reduce this effect from the orbital trajectories of particles, it gives a more reasonable magnitude range than absolute dispersion, and so is recommended in estimating vertical dispersion rates.

• Journal of the Korean earth science society
• /
• v.29 no.3
• /
• pp.280-289
• /
• 2008

An extensive set of both in-situ and satellite data regarding oceanic sea surface temperatures in Northeast Asian seas, collected over a 10-year period, was collocated and surveyed to assess the accuracy of satellite-observed sea surface temperatures (SST) and investigate the characteristics of satellite measured SST errors. This was done by subtracting insitu SST measurements from multi-channel SST (MCSST) measurements. 845 pieces of collocated data revealed that MCSST measurements had a root-mean-square error of about 0.89$^{\circ}C$ and a bias error of about 0.18$^{\circ}C$. The SST errors revealed a large latitudinal dependency with a range of $\pm3^{\circ}C$ around 40$^{\circ}N$, which was related to high spatial and temporal variability from smaller eddies, oceanic currents, and thermal fronts at higher latitudes. The MCSST measurements tended to be underestimated in winter and overestimated in summer when compared to in-situ measurements. This seasonal dependency was discovered from shipboard and moored buoy measurements, not satellite-tracked surface drifters, and revealed the existence of a strong vertical temperature gradient within a few meters of the upper ocean. This study emphasizes the need for an effort to consider and correct the significant skin-bulk SST difference which arises when calculating SST from satellite data.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.233-246
• /
• 2009

The partitioning of evapotranspiration (ET) into evaporation (E) and transpiration (T) is critical in understanding the water cycle and the couplings between the cycles of energy, water, and carbon. In forests, the total ET measured above the canopy consists of T from both overstory and understory vegetation, and E from soil and the intercepted precipitation. To quantify their relative contributions, we have measured ET from the floors of deciduous and coniferous forests in Gwangneung using eddy covariance technique from 1 June 2008 to 31 May 2009. Due to smaller eddies that contribute to turbulent transfer near the ground, we performed a spectrum analysis and found that the errors associated with sensor separation were <10%. The annual sum of the understory ET was 59 mm (16% of total ET) in the deciduous forest and 43 mm (~7%) in the coniferous forest. Overall, the understory ET was not negligible except during the summer season when the plant area index was near its maximum. In both forest canopies, the decoupling factor ($\Omega$) was about ~0.15, indicating that the understory ET was controlled mainly by vapor pressure deficit and soil moisture content. The differences in the understory ET between the two forest canopies were due to different environmental conditions within the canopies, particularly the contrasting air humidity and soil water content. The non-negligible understory ET in the Gwangneung forests suggests that the dual source or multi-level models are required for the interpretation and modeling of surface exchange of mass and energy in these forests.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.8 no.2
• /
• pp.94-110
• /
• 2003

This study utilizes the dataset of Topex/Poseidon(T/P) altimeter sea surface height (1992-2000 yr., 286 cycles)to investigate the tempore-spatial variability in the East (Japan) Sea. Optimal interpolation (Ol) technique was applied to the pre-processed T/P dataset (level 2) to produce sea surface height anomaly (SSHA) map on regular grids. Spectral analyses of the timeseries of the SSHA at chosen stations and empirical orthogonal function (EOF) analysis of the SSHA in the entire East Sea were made. Distribution of the SSHA can be divided by the southern and northern regions sharply by the polar front situated in the middle of the East Sea. The southern region under the direct influence of the Tsushima Current exhibits higher amplitude of the SSHA fluctuation, while the northern region does relatively smaller one. The spatio-temporal variability of the SSHA in the East Sea can be characterized by the five modes of the EOFs accounting for more than 85% of the total variance. The first mode dominates the SSHA variation in the entire domain with strong seasonal and inter-annual periods accounting for the 72.3% of the total variance. The other modes (up to 5th account for 14%) are responsible for the SSHA variation associated with the local current system, meandering of the polar frontal axis, and mesoscale eddies. Spectral peaks with significant confluence level show semi-annual, annual and interannual (2, 3-4 years) periods.

• Korean Journal of Remote Sensing
• /
• v.20 no.4
• /
• pp.245-252
• /
• 2004

The mechanism of cold core eddy formation was investigated using boundary conditions between the East China coastal cold water and the Kuroshio Warm Current, wind data related to the monsoon which was measured by QuikSCAT, and the bottom topography of the East China Sea. When winds blow from the southeast at an intensity comparable to that in the winter period in 1999 and 2003, the warm Kuroshio and Tsushima Current became stronger, and temperatures were considerably higher than those of the extended cold water of the coast of the East China. At that time, the cold water was captured by warm water from the Kuroshio and the Tsushima Current. This facilitated the formation of mesoscale cold core eddies with diameter of 150km in the East China Sea in May, 1999 and February, 2003. The cold core eddy which was detected by NOAA, SeaWiFS and QuikSCAT satellites. The East China Sea is considered to be important not only as a good fishing ground but also nursery and spawning area for many kinds of fishes. Therefore, it would be worth studying spatio-temporal variations of the cold core eddy in the environmental conditions of the northwestern East China Sea using systematic remote sensing techniques.