• Journal of Digital Convergence
• /
• v.13 no.9
• /
• pp.209-218
• /
• 2015

Computer simulation plays an important role for a theoretical foundation in convergence technology and the interpolation is to know the unknown values from known values on grid points. Therefore it is an important problem to select an interpolation method for digital simulation. The aim of this paper is to compare analysis of interpolation methods for digital simulation. we test six different interpolation methods namely: Quartic-Lagrangian, Cubic Spline, Fourier, Hermit, PWENO and SL-WENO. Through digital simulation of a linear advection equation, we analyse pros and cons for each method. In order to compare performance, we introduce accuracy computing and Error functions. The accuracy computing is used well-known $L^1-norm$ and the Error functions are dispersion function, dissipation function and total error function. High-order methods well apply to computer simulation, unfortunately, side-effects (Oscillation) happen.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.38-47
• /
• 1995

The three-dimensional multi-layer computer simulation technology for tidal current developed in the previous study is updated to a new version. many improvements are achieved by following changes : (1) No-reflection condition is adopted instead of no-gradient condition as an open boundary condition. (2) Time marching algorithm is changed so that velocity and pressure(surface movement) might be salved in turn at different time step (3) Convection term in equation of motion is estimated by upwind differencing scheme instead of central differencing. The stability is improved considerably and the steady state is achieved within 2 tidal periods which is about 3 times shorter than that of the old version. Moreover, fluctuations in time disappeared by introducing the new time marching technique. An application to the real near shore area(near Inchon harbor) is performed by the new version. Simulated results are compared with those by the simulation total developed in the University of Tokyo. Validity and effectiveness of the two simulation technologies are chocked through the comparative research works.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.22 no.1
• /
• pp.1.1-15
• /
• 1994

The transports of the seasonal freshwater and salt from surface to 500 m depth in the tropical Atlantic Ocean are derived from the equations of the continuity and saltconservation, respectively. The freshwater transport is obtained by southward integration of the divergence of surface freshwater flux, using climatological freshwater(i. e. precipitation, evaporation, and river discharge) data. The annual freshwater transport is northward, ranging from 0 Sv near the equator to 0.3 Sv at $12^{\circ}{\;}N{\;}and{\;}20^{\circ}{\;}S$. The seasonal meridional transport amounts of freshwater range from 1.35 Sv to-0.45 Sv. The strong northward freshwater transports prevail for the intraseasonal period summer to fall. This seasonal cycle is caused by the shifts of the ITCZ as well as the changes in the local freshwater storage. Annual and seasonal salt transports are calculated from objectively analyzed historical (1900-86) salinity observations. The annual salt flux in the ocean zero, showing that the salt flux by horizontal advection balances the flux by horizontal diffusion. The salt flux due to the diffusion is northward, and has a maximum of $5{\;}{\times}{\;}10^6kg/s$ at 15oN. Seasonal transport amounts of salt range from $30{\;}{\times}{\;}10^6kg/s{\;}to{\;}-35{\;}{\times}10^6kg/s$. The direction of the seasonal salt transports is northward except for the intraseasonal period summer to fall.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.289-301
• /
• 2009

In this work, wave run-up heights and resultant wave forces on a vertical revetment due to tsunami (solitary wave) are investigated numerically using a numerical wave tank model called CADMAS-SURF (CDIT, 2001. Research and Development of Numerical Wave Channel (CADMAS-SURF). CDIT library, No. 12, Japan.), which is based on a 2-D Navier-Stokes solver, coupled to a volume of fluid (VOF) method. The third order approximate solution (Fenton, 1972. A ninth-order solution for the solitary wave. J. of Fluid Mech., Vol. 53, No.2, pp.257-271) is used to generate solitary waves and implemented in original CADMAS-SURF code. Numerical results of the wave profiles and forces are in good agreements with available experimental data. Using the numerical results, the regression curves determined from the least-square analysis are proposed, which can be used to determine the maximum wave run-up height and force on a vertical revetment due to tsunami. In addition, the capability of CADMAS-SURF is demonstrated for tsunami wave forces acting on an onshore structure using various configuration computations including the variations of the crown heights of the vertical wall and the position of the onshore structure. Based on the numerical results such as water level, velocity field and wave force, the direct effects of tsunami on an onshore structure are discussed.

• 한국해양학회지
• /
• 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.4C
• /
• pp.231-238
• /
• 2008

A column testing device capable of measuring the electrical resistivity of soil at 3 different locations was developed to verify applicability of bulk electrical conductivity (BEC) breakthrough curves in monitoring contaminant transport. Tracer injection tests were conducted with three different types of saturated sands to obtain average linear velocities and longitudinal hydrodynamic dispersion coefficients based on BEC breakthrough curves and effluent solute breakthrough curves. Comparative analysis of transport parameters obtained from curve fitting the results into the analytical solutions confirmed the validity of resistance measurements in estimating time-continuous resident solute concentration. Under the assumption that a linear relationship exists between ${\sigma}_{sat}-{\sigma}_w-C$, the BEC breakthrough curves are able to effectively reduce the laborious and time-consuming processes involved in the conventional method of sampling and analysis. In order to reduce possible uncertainties in analyzing the BEC breakthrough curves, it was recommended that resistance measurements take place nearby the effluent boundary. In addition, a sufficient electrical contrast or difference in the electrical conductivity of the influent and the saturating solution is required to conduct reliable analysis.