• 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 the Korean earth science society
• /
• v.39 no.4
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.1 no.1
• /
• pp.20-28
• /
• 2009

In thus paper we validate a numerical model for wave-structure interaction by comparing numerical results with laboratory data. The numerical model is based on the Navier-Stokes (N-S) equations for an incompressible fluid. The N-S equations are solved by a two-step projection finite volume scheme and the free surface displacements are tracked by the volume of fluid (VOF) method The numerical model is used to simulate solitary waves and their interaction with a group of slender vertical piles. Numerical results are compared with the laboratory data and very good agreement is observed for the time history of free surface displacement, fluid particle velocity and wave force. The agreement for dynamic pressure on the cylinder is less satisfactory, which is primarily caused by instrument errors.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.28 no.7
• /
• pp.293-298
• /
• 2016

A three-dimensional (3D) numerical model of the vertical ground-coupled heat exchanger is useful for analyzing the modern ground source heat pump system. Furthermore, a detailed description of the inner side of the exchanger allows to account for the effects of the thermal capacity. Thus, both methods are included in the proposed numerical model. For the ground portion, a FDM (Finite Difference Method) scheme has been applied using the Cartesian coordinate system. Cylindrical grids are applied for the borehole portion, and the U-tube configuration is adjusted at the grid, keeping the area and distance unchanged. Two sub-models are numerically coupled at each time-step using an iterative method for convergence. The model is validated by a reference 3D model under a continuous heat injection case. The results from a periodic heat injection input show that the proposed thermal capacity model reacts more slowly to the changes, resulting in lower borehole wall temperatures, when compared with a thermal resistance model. This implies that thermal capacity effects may be important factors for system controls.

• 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.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.1
• /
• pp.29-41
• /
• 1996

The two-stage numerical model was used to study the relation between three-dimensional local wind model, advection/diffusion model of random walk method and second moment method on Pusan coastal area. The first stage is three dimensional time-dependent local wind model which gives the wind field and vertical dirrusion coefficient. The second stage is advection/diffusion model which uses the results of the first stage as input data. First, wind fields on Pusan coastal area for none synoptic scale wind showed typical land and sea breeze circulation, and convergence zone occured at 1200LST in northern of domain, in succession, moved northward of domain. Emissions from Sinpyeong industrial district were trasnported toward the inland by sea breeze during daytime, and reached the end part of domain about 1800LST. During nighttime, emissions return to sea by land breeze and vertical diffusion also contributes to upward transport. In order to use this model for forecast of air pollution concentration on the Pusan coastal area, it is necessary that computed value must be compared with measured value and wind fields model must also be dealt in detail.

• Water for future
• /
• v.26 no.4
• /
• pp.61-71
• /
• 1993

To investigate the structure of internal flow hydrodynamically, the complete vertical equation of motion should be assembled into the model. In the present study a numerical simulation model not hydrodynamically approximated is established. From the comparison of the predicted results with the computed results from k-$ two equation turbulence model by Huh et. al.(1991)and the experimental data by Nakatsuji(1984), the vertical acceleration and its effects on the development of buoyant surface jets are evaluated quantitatively.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.113-122
• /
• 1987

A numerical simulation technique of three-dimensional finite difference model is developed to study the groundwater flow system in Dcjima, an upland area which faces Kasumigaura Lake. For general perspectives of the groundwater flow system, a steady state three-dimentional model is simulated. For the sedimentary mud formations which are found in the representative formation, three situations of hydraulic conductivity are considered, representing an isotropic condition and situations where the horizontal permeability is equal to 10 times and 100times of the vertical one. The finite difference grid used in the simulation has 60x50x30=90,000 nodes. A converged solution with a tolerance of 0.001 meter of hydraulic head is set. Having determined the flow net by using a steady state three-dimensional model. the results for the three cases of hydraulic conductivity are compared with the results of tracer methods (Bae and Kayane 1987) With the aid of four representative vertical cross-sections, groundwater flow systems in the study area are assumed. Water balances for the three cases indicate very good agreement between total recharge and discharge in each case Analyses of groundwater flow system based on the tritium concentrations and water quality measurements (Bae and Kayane 1987) are confirmed by the numerical simulation and the results obtained by these two methods appeared to be in close agreement.

• Wind and Structures
• /
• v.10 no.4
• /
• pp.315-346
• /
• 2007

A preceding companion article introduced the slot jet approach for large-scale quasi-steady modelling of a downburst outflow. This article extends the approach to model the time-dependent features of the outflow. A two-dimensional slot jet with an actuated gate produces a gust with a dominant roll vortex. Two designs for the gate mechanism are investigated. Hot-wire anemometry velocity histories and profiles are presented. As well, a three-dimensional, subcloud numerical model is used to approximate the downdraft microphysics, and to compute stationary and translating outflows at high resolution. The evolution of the horizontal and vertical velocity components is examined. Comparison of the present experimental and numerical results with field observations is encouraging.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.573-574
• /
• 2006

This research deals with the wave transmission and dissipation problems for two dimensional regular waves and s vertical circular cylindersr. Using the unsteady mild slope equation, a numerical model has been developed to calculate the reflection and transmission of regular waves from a multiple-row vertical circular cylinders. In addition hydraulic model experiments have been conducted with different values of properties between the cylinders and opening ratio (distances) between the rows of the cylinders. It is found that the transmission coefficients decrease with decreasing the opening ratio and increasing the rows of vertical cylinders. Comparison between hydralic and numerical experiments results shows resonable agreement.