• Journal of Wetlands Research
• /
• v.16 no.3
• /
• pp.327-336
• /
• 2014

Submerged thin walls are extreme case of submerged rectangular blocks, and could be used for many purposes in rivers or coastal zones, e.g. to tsunami. To understand flow characteristics including flow and pressure fields around a specific submerged thin wall a numerical model was applied which includes computation of hydrodynamic pressure on ${\sigma}$-coordinate. ${\sigma}$-coordinate has strong merits for simulation of subcritical flow over mild-sloped beds. On the other hand ${\sigma}$-coordinate is quite poor to treat sharp structures on the bed. There have been a few trials to incorporate dynamic pressure in ${\sigma}$-coordinate by some researchers. One of the previous approaches includes process of sloving the Poisson equation. However, the above method includes many high-order terms, and requires long cpu for simulation. Another method SOLA was developed by Hirt et al. for computation of dynamic pressure, but it was valid for straight grid system only. Previous SOLA was modified for ${\sigma}$-coordinate for the present purpose and was adopted in a model system, CST3D. Computed flow field shows reasonable behaviour including vorticity is much stronger than the upstream and downstream of the structure. The model was verified to laboratory experiments at a 2DV flume. Time-average flow vectors were measured by using one-dimensional electro-magnetic velocimeter. Computed flow field agrees well with the measured flow field within 10 % error from the speed point of view at 5 profiles. It is thought that the modified SOLA scheme is useful for ${\sigma}$-coordinate system.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1263-1267
• /
• 2007

A more complete two-dimensional vertical numerical model has been developed to describe the saltwater intrusion in an estuary. The model is based on the previous studies in order to obtain a better accuracy. The non-linear terms of the governing equations are analyzed and the $\sigma$-coordinate system is employed in the vertical direction with full transformation which is recently issued in several studies because numerical errors can be generated during the coordinate transformation of the diffusion term. The advection terms of the governing equations are discretized by an upwind scheme in second-order of accuracy. By employing an explicit scheme for the longitudinal direction and an implicit scheme for the vertical direction, the numerical model is free from the restriction of temporal step size caused by a relatively small grid ratio. In previous researches, some terms induced from the transformation have been intentionally excluded since they are asked the complicate discretization of the numerical model. However, the lack of these terms introduces significant errors during the numerical simulation of scalar transport problems, such as saltwater intrusion and sediment transport in an estuary. The numerical accuracy attributable to the full transformation is verified by comparing results with a previous model in a simply sloped topography. The numerical model is applied to the Han River estuary. Very reasonable agreements for salinity intrusion are observed.

• Journal of Environmental Science International
• /
• v.20 no.6
• /
• pp.679-686
• /
• 2011

The characteristics of circulation in the coastal area of Jeju Harbor in Korea was examined using the Princeton Ocean Model(POM) with a sigma coordinate system. The result of numerical analysis well corresponded to the observed current data. The velocity at offshore was stronger compared to coastal area during the both period of in maximum flood and maximum ebb of spring tide. According to mean wind velocity, the tidal velocity at the shallow area of Jocheon was slightly increasing during maximum ebb. The effect of wind on the circulation was stronger in shallow area and showed rapid change with depth.

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

POM(Princeton Ocean Model) was utilized in this study because it took ${\sigma}-coordinate$ system which could predict the behavior of bottom water. The model has been increasingly applied to costal area although it was initially developed as the ocean flow model. The original POM did not correct computational errors in transformation of ${\sigma}-coordinate$ system. The trying to reduce conversional errors might improve accuracy of flow velocity in vicinities of bottom layer. Therefore, in this study it was proposed to modify the original POM by using error correction method suggested by $Sl{\Phi}rdal$(1997). The modified POM was applied to Young-rang Lake, one of the typical brackish lakes in Korea. It was found that the behavior of bottom water could be well predicted. Thus, it seems that the modified POM can be used as a useful tool to clarify the mechanism of formation and behavior of bottom water including oxygen-deficient water mass.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.14 no.3
• /
• pp.176-183
• /
• 2011

Applicabilities of two 3-D hydrodynamic models on Chungkye Bay in which tidal flats are well developed were examined. Both EFDC model and ESCORT model with the ${\sigma}$-coordinate showed fairly good results. However, their efficiencies were lowered especially on the tidal flats due to dense vertical grids. This inefficiency could be overcome by using the z-coordinate of the ESCORT model keeping similar accuracy.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1145-1158
• /
• 1994

A three-dimensional, finite difference, numerical model with free surface was developed on ${\sigma}$-coordinate. A semi-implicit numerical scheme in time has been adopted for computational efficiency. The scheme is essentially independent of the stringent stability criteria (CFL condition) for explicit schemes of external surface gravity wave. Implicit algorithm was applied for vertical shear stress, Coriolis force and pressure gradient terms. The reliability of the model with vertically variable grid system was checked by the comparison of simulation results with analytic solution of wind-driven currents in a one-dimensional channel. Sensitivity analysis of differencing parameters was carried out by applying the model to the calculation of wind-driven currents in a square lake.

• Journal of Climate Change Research
• /
• v.2 no.2
• /
• pp.79-91
• /
• 2011

In order to study interactions between climate change and air quality, a modeling system including the downscaling scheme has been developed in the integrated manner. This research focuses on the development of a downscaling method to utilize CCSM3 outputs as the initial and boundary conditions for the regional climate model, MM5. Horizontal/vertical interpolation was performed to convert from the latitude/longitude and hybrid-vertical coordinate for the CCSM3 model to the Lambert-Conformal Arakawa-B and sigma-vertical coordinate for the MM5 model. A variable diagnosis was made to link between different variables and their units of CCSM and MM5. To evaluate the dynamic downscaling performance of this study, spatial distributions were compared between outputs of CCSM/MM5 and NRA/MM5 and statistic analysis was conducted. Temperature and precipitation patterns of CCSM/MM5 in summer and winter showed a similar pattern with those of observation data in East Asia and the Korean Peninsula. In addition, statistical analysis presented that the agreement index (AI) is more than 0.9 and correlation coefficient about 0.9. Those results indicate that the dynamic downscaling system built in this study can be used for the research of interaction between climate change and air quality.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.3
• /
• pp.253-265
• /
• 2007

This study applied the previous results of the NAO model, a tidal correction model, to the open boundary condition for the behavior of tidal elevation and tidal currents in the South Sea. This study used the EFDC model considering the wetting and drying problem and using the $\sigma-coordinate$ as a vertical coordinate and generated two mesh cases of the constant grid size of 2.0 km and the variable grid size of $0.5\sim2.0km$. The numerical results for the tides showed that the predicted results were in quite good agreements with the observational data acquired from the tidal stations of the NORI. The predicted tides were observed to propagate from the east area to the west area in the South Sea. The verification results reveal that the numerical results are more correlated with the measured tidal data as the grid size decreases. The grid size of 2 km results in proper simulation of tidal currents in wide waterway and offshore area whereas the numerical results from the grid size of 0.5 km tend to somewhat underestimate the tidal currents affected by narrow waterway and topography in inner-bay.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.5
• /
• pp.217-227
• /
• 2017

The object of this study is to estimate the net volume transport and the residual flow that changed by space and time at southern part of Yeomha channel, Gyeonggi Bay. The cross-section observation was conducted at the mid-part (Line2) and the southern end (Line1) of Yeomha channel for 13 hours during neap and spring-tides, respectively. The Lagrange flux is calculated as the sum of Eulerian flux and Stokes drift, and the residual flow is calculated by using least square method. It is necessary to unify the spatial area of the observed cross-section and average time during the tidal cycle. In order to unify the cross-sectional area containing such a large vertical tidal variation, it was necessary to convert into sigma coordinate system by horizontally and vertically for every hour. The converted sigma coordinate system is estimated to be 3~5% error when compared with the z-level coordinate system which shows that there is no problem for analyzing the data. As a result, the cross-sectional residual flow shows a southward flow pattern in both spring and neap tides at Line2, and also have characteristic of the spatial residual flow fluctuation: it northwards in the main line direction and southwards at the end of both side of the waterway. It was confirmed that the residual flow characteristics at Line2 were changed by the net pressure due to the sea level difference. The analysis of the net volume transport showed that it tends to southwards at $576m^3s^{-1}$, $67m^3s^{-1}$ in each spring tide and neap tide at Line2. On the other hand, in the control Line1, it has tendency to northwards at $359m^3s^{-1}$ and $248m^3s^{-1}$. Based on the difference between the two observation lines, it is estimated that net volume transport will be out flow about $935m^3s^{-1}$ at spring tide stage and about $315m^3s^{-1}$ at neap tide stage as the intertidal zone between Yeongjong Island and Ganghwa Island. In other words, the difference of pressure gradient and Stokes drift during spring and neap tide is main causes of variation for residual current and net volume transport.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.3
• /
• pp.116-121
• /
• 2005

In the present paper, a new calculation algorithm far solving large scale environmental or geophysical flows with free surface is proposed where the non-hydrostatic pressure component is taken into consideration. Predictor-corrector fractional step approach with explicit, forward time marching scheme in the sigma coordinate system is employed. In order to validate the present calculation algorithm and to estimate the effects of non-hydrostatic pressure on resultant flow and free surface movements, example calculations are carried out for typical steady and unsteady flow problems. Present method can be applied to the meso-scale free surface flows with complex bottom topography where MAC-like 3-d hydrodynamic calculations are quite ineffective and uneconomic.