• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.1
• /
• pp.131-141
• /
• 1994

Three characteristics-based split-operator methods were applied to a longitudinal pollutant dispersion problem, and the results were compared with those of several Eulerian schemes. The split-operator methods consisted of generalized upwind, two-point fourth-order and sixth-order Holly-Preissmann schemes, respectively, for the advection calculation, and the Crank-Nicholson scheme for the diffusion calculation. Compared with the Eulerian schemes tested, split-operator methods using the Holly-Preissmann schemes gave much more accurate computational results. Eulerian schemes using centered difference approximations for the advection term resulted in numerical oscillations, and those using backward difference resulted in numerical diffusion, both of which were more severe for smaller value of the longitudinal dispersion coefficient.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.240-240
• /
• 2015

As Simpson et al. (1990) emphasized the importance of the straining process in the stratification and mixing in the estuarine circulation process, various researches have investigated on the relative contribution of each process to the overall potential energy anomaly dynamics. However, many numerical works have done only for two dimensional modeling along channel or the short distance cross sectional three dimensional simulations as Burchard et al. (2008) and the estuarine channel was not simulated so far. But, in the study on the physics of shallow coastal seas, spatial dimension in the three dimensional way affects significantly on results of a particular numerical model. Therefore, the comparison of two and three dimensional models is important to understand the real physics of mixing and stratification in an estuary. Also, as Geyer and MacCready (2013) pointed out that the lateral process seems to be important in determining the periodic stratifications, to study such process the three dimensional modeling must be required. The present study uses a numerical model to show the signification roles of each term of the time-dependent dynamic equation for the potential energy anomaly (PEA) in controlling along and lateral channel flows and different stratification structures. Moreover, we present the relationships between the ${\Phi}$-advection, the depth mean straining, vertical mixing and vertical advection can explain well how water level, salinity distribution and across velocity 2D model are slightly different from 3D.

• Fisheries and Aquatic Sciences
• /
• v.26 no.11
• /
• pp.639-648
• /
• 2023

The sea surface temperature (SST) and ocean heat content in the Korea Waters are gradually increased. Especially the increasing trend of annual mean SST in the Korea Water is higher about 2.6 times than the global mean during past 55 years (1968-2022). Before 2010s, the increasing trend of SST was led by winter season in the Korea Waters. However, this pattern was clearly changed after 2010s. The increasing trend of SST during summer is higher about 3.9 times than during winter after 2010s. We examine the long-term variations of several ocean and climate factors to understand the reasons for the long-term pattern changes of SST between summer and winter in recent. Tsushima warm current was significantly strengthened in summer compare to winter during past 33 years (1986-2018). The long-term patterns of Siberian High and East Asian Winter Monsoon were definitely changed before and after early- or mid-2000s. The intensities of those two climate factors was changed to the increasing trend or weakened decreasing trend from the distinctive decreasing trend. In addition, the extreme weather condition like the heatwave days and cold spell days in the Korea significantly increased since mid- or late-2000s. From these results, we can consider that the occurrences of frequent and intensified marine heatwaves during summer and marine cold spells during winter in the Korea Waters might be related with the long-term pattern change of SST, which should be caused by the long-term change of climate factors and advection heat, in a few decade.

• Journal of the Korean earth science society
• /
• v.27 no.3
• /
• pp.302-312
• /
• 2006

This study has investigated the effect of Arctic Oscillation (AO) on cold surge through atmospheric circulation and heat source analysis for the past winters from 1979 to 2004. The number of occurrence of cold surge in the negative AO phase is about 14.3% larger than that in the positive AO phase. The number of occurrence of cold surge per a month in the negative (positive) AO phase is about 1.33 (1.05), respectively, indicating that the negative AO phase has about 26.6% larger occurrence than the positive AO phase. It means that the cold surge has occurred frequently in particular months with the negative AO phase. And it also shows that surface temperature in the negative AO phase is about $0.6^{\circ}C$ lower than in positive AO phase. As a result of the analysis for the difference of heat source according to the intensity of AO, it shows that surface air temperature around the Korean peninsula in the negative AO phase is more lower than in positive AO phase by the intensification of cold advection term. However, heat source term cancels out the cooling effect by cold advection term, indicating that it suppresses the decrease in surface air temperature. It results in a small difference of $0.6^{\circ}C$ in surface air temperature between the positive and negative AO phase in spite of the significance of atmospheric circulation change.

• Journal of Korean Port Research
• /
• v.9 no.1
• /
• pp.45-56
• /
• 1995

Investigation of the physical behavior of dredged material disposal in coastal water includes estimations of water column concentration in the receiving water, exposure time, the initial deposition pattern as well as thickness of material at the dumping fields near the estuary area. Calculation based on vertical setting and horizontal advection of single particles ignore the effects of bulk properties of the disposed material, vertical and horizontal diffusion, and material dilution due to the entrainment of ambient water during descent. This paper focuses on the spatial and temporal changes in the dumping fields for the water column and bottom at a hypothetically confined coastal water, where the ambient time-invariant velocity and density profiles are applied, within the initial time period following the instantaneous release of the dredged material. This model accounts the behavior of material after release divided into three phases: convective descent, dynamic collapse and long-term passive dispersion

• Honam Mathematical Journal
• /
• v.38 no.1
• /
• pp.47-57
• /
• 2016

In the discrete formulation of the bubble stabilized Legendre Galerkin methods, the system of equations includes the artificial viscosity term as the parameter. We investigate the estimation of this parameter to get the least-squares solution which minimizes the sum of the squares of errors at each node points. Some numerical results are reported.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.49-54
• /
• 1997

Rice's monotone streamline upwind finite element method, which was proposed to treat convection-dominated flows, is applied to the linear triangular element. An alignment technique of unstructured grids with given velocity fields is used to prevent the interpolation error produced in evaluating the convection term in the upwind method. The alignment of grids is accomplished by optimizing a target function defined with the inner-product of a properly chosen side vector in the element with the velocity field. Two pure advection problems are considered to demonstrate the superiorities of the present approach in solving the convection-dominated flow on the unstructured grid. Solutions obtained with aligned grids are much closer to the exact solutions than those with initial regular grids. The capability of the present approach in predicting the appearance of the secondary vortex in the laminar confined jet impingement is shown by comparing streamlines to those produced by SIMPLE on a highly stretched grid toward the impingement plate.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Ocean and Polar Research
• /
• v.38 no.2
• /
• pp.89-102
• /
• 2016

In this study, we evaluated the model performance with respect to Sea Surface Temperature (SST) and Net Heat Flux (NHF) by considering the characteristics of seasonal temperature variation and contributing factors and by analyzing heat budget terms in the Northwestern Pacific and East Asian Marginal Seas ($110^{\circ}E-160^{\circ}E$, $15^{\circ}N-60^{\circ}N$) using the HadGEM2-AO historical run. Annual mean SST of the HadGEM2-AO is about $0.065^{\circ}C$ higher than observations (EN3_v2a) from 1950 to 2000. Since 1960, the model has simulated well the long-term variation of SST and the increasing rate of SST in the model ($0.014^{\circ}C/year$) is comparable with observations ($0.013^{\circ}C/year$). Heat loss from the ocean to the atmosphere was simulated slightly higher in the HadGEM2-AO than that in the reanalysis data on the East Asian Marginal Seas and the Kuroshio region. We investigated the causes of temperature variation by calculating the heat budget equation in the two representative regions. In the central part of the Kuroshio axis ($125^{\circ}E-130^{\circ}E$, $25^{\circ}N-30^{\circ}N$: Region A), both heat loss in the upper mixed layer by surface heat flux and vertical heat advection mainly cause the decrease of heat storage in autumn and winter. Release of latent heat flux through the heat convergence brought about by the Kuroshio contributes to the large surface net heat flux. Positive heat storage rate is mainly determined by horizontal heat advection from March to April and surface net heat flux from May to July. In the central part of the subtropical gyre ($155^{\circ}E-160^{\circ}E$, $22^{\circ}N-27^{\circ}N$: Region B), unlike Region A, vertical heat advection predominantly causes the decrease of heat storage in autumn and winter. In spring and summer, surface heat flux contributes to the increase of heat storage in Region B and the period is two times longer than the period for Region A. In this season, shoaling of the mixed layer depth plays an important role in the increase of SST.

• 한국해양학회지
• /
• v.30 no.6
• /
• pp.565-583
• /
• 1995

We have used a nonlinear quasi-geostrophic model to study effects of lateral friction and bottom topography on the motion of warm eddies. The two empirical orthogonal functions of the stream function, accounting for the vertical structure, represent the barotropic and first baroclinic dynamic modes. This model is integrated 360 days on a 1000 km ${\times}$ 1000 km domain with a resolution of 10 km ${\times}$ 10 km including both the thermocline and idealized topography of the East Sea. Prescribed inflow through the Korea Strait is compensated by outflow through the Tsugaru Strait. The balance between the nonlinear advection term and the planetary ${\beta}$-effect tends to make northward movement of warm eddy over a flat bottom. The motion of a warm eddy over a sloping topography can be dominated by the nonlinear advection, while nonlinearity plays a secondary role over a flat topography. For eddies dispersing over topography, the nonlinear tendency is a function of time. For a strong warm eddy, northward propagation can occur. For intermediate strength of eddies one might expect a balance between the nonlinear term and the topographic ${\beta}$-effect. As nonlinearity decreases with eddy dispersion, southward motion along the slope may occur by such as a topographic Rossby wave. Our numerical simulations have confirmed the importance of lateral friction on eddy motions, in such a way that the northward penetration of the warm eddy increases drastically by the decrease of the lateral friction. The northward motion of warm eddy can be prevented by reducing the Reynolds number sufficiently. We have also demonstrated the crucial role of topographic effects in the eddy motion process.