• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.4
• /
• pp.209-216
• /
• 1991

Characteristics of suspended sediment concentration of equilibrium and non-equilibrium state caused by waves and currents are investigated by conducting a movable bed experiments in wave tanks. In the region where a downward flux of suspended sediment is larger than a upward flux, time-averaged vertical distribution of suspended sediment does not indicate logarithmic distribution. A new numerical procedure for predicting time-averaged suspended sediment concentration is also proposed based on two-dimensional advective diffusion equation by applying a split-operator approach. It is found that the unposed procedure can predict measured distribution of suspended sediment satisfactorily.

• Journal of Radiation Protection and Research
• /
• v.12 no.2
• /
• pp.28-36
• /
• 1987

Gaussian plume model is used to assess environmental dose for abnormal radioactive release in nuclear facility, but there has a problem to use it for complex terrain. In this report, MATTEW and WIND04 Codes which had been verified were used to calculate wind field in the complex terrain. Under the base of these codes principle, wind fields were obtained from the calculation of the finite difference approximation for advection-diffusion equations which satisfy the mass-conservative law. Particle concentrations and external doses were calculated by using PIC model which approximate the particle to radioactive cloud, and atmospheric diffusion of the particles from the random walk method. The results show that the adjusted wind fields and the distributions of the exposure dose vary with the topography of the complex terrain.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.10 no.4
• /
• pp.237-246
• /
• 2012

In this study, numerical model for transport of radionuclide and colloid was developed. In order to solve reaction-migration governing equation for colloid and radionuclide, Strang-splitting Sequential Non-Iterative (SNI), which is one of Operator Splitting Method, was used for numerical method and this was coded by MATLAB. From the verification by comparing the simulation results with analytical solution considering only solute transport and rock diffusion, the Pearson's correlation coefficient was greater than 0.99 which demonstrates the accuracy of the model.

• Water for future
• /
• v.24 no.2
• /
• pp.81-95
• /
• 1991

For investigation of the interaction of tide and river flow, the derived equations are solved analytically using the approximation method of perturbation. The convective inertia and nonlinear bottom friction terms are included in the derivations. The harmonic analysis is applied to decompose the complicated interaction of the freshwater discharge with various constituents of tide into its individual interaction with each constituent. In this study, four main constituents(M2, S2, Kl, 01) are included. The relations of dimensionless parameters of the tide, especially the dimensionless damping modulus, are then determined for each solution. The results show that analytical solution of dimensionless damping modulus underestimates the measured value obtained from harmonic analysis. Results of water level obtained by applying the analytical model to a tidal reach of the Keum River in the years 1981 and 1982 show very good agreement with those obtained from the harmonic analysis.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.208-215
• /
• 2003

Convective-diffusion equations appear in various disciplines such as hydrology, chemical engineering and oceanography dealing with the transport problem of scalar quantities. Since it is nonlinear, numerical methods are generally used to obtain its solution. Very limited number of analytical solutions are available usually in cases when the convective velocity is constant or has a simple functional form (for some collection of the solutions, see Noye, 1987). There is however a continuing need to develop analytical solutions because of its practical importance. Analytical solutions of the convection-diffusion equation are valuable not only for the better understanding on the transport process but the verification of numerical schemes. (omitted)

• Journal of KIISE:Computer Systems and Theory
• /
• v.32 no.11_12
• /
• pp.692-704
• /
• 2005

With the improvement of computer hardware, GPUs(Graphics Processor Units) have tremendous memory bandwidth and computation power. This leads GPUs to use in general purpose computation. Especially, GPU implementation of compute-intensive physics based simulations is actively studied. In the solution of differential equations which are base of physics simulations, tridiagonal matrix systems occur repeatedly by finite-difference approximation. From the point of view of physics based simulations, fast solution of tridiagonal matrix system is important research field. We propose a fast GPU implementation for the solution of tridiagonal matrix systems. In this paper, we implement the cyclic reduction(also known as odd-even reduction) algorithm which is a popular choice for vector processors. We obtained a considerable performance improvement for solving tridiagonal matrix systems over Thomas method and conjugate gradient method. Thomas method is well known as a method for solving tridiagonal matrix systems on CPU and conjugate gradient method has shown good results on GPU. We experimented our proposed method by applying it to heat conduction, advection-diffusion, and shallow water simulations. The results of these simulations have shown a remarkable performance of over 35 frame-per-second on the 1024x1024 grid.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.2 no.2
• /
• pp.20-30
• /
• 1999

To improve water quality particularly for sea bathers along the Fylde coastal zone near Blackpool, North West England, waste water from a sewage outfall is studied using a mathematical model. The explicit second order accurate central scheme and the third order accurate QUICKEST scheme are used to represent the diffusion terms and the advection terms of the advective-diffusion equation, respectively. Hydrodynamic model is run for a coarse and fine grid, of 1km and 200m, respectively, obtaining good agreement with measured data. Water quality model is then used to predict faecal coliform levels in the region for four different scenarios, including discharges from: - (i) Fleetwood outfall, (ii)River Ribble for summer condition, (iii)River Ribble for winter condition, and (iv)combined sewer overflows for the Blackpool and Fleetwood communities. Main findings from the simulations are:- (i) Fleetwood outfall has a negligible impact on the beaches with respect to pathogen levels; (ii) Discharge from River Ribble for both summer and winter conditions is predicted in the range of coliform levels 10 -500 counts/100ml along the beach at Lytham St. Annes; and (iii) The CSO effluent discharges are predicted not to advect out into offshore by stronger tidal currents.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.123-136
• /
• 1996

A three-dimensional Galerkin-FEM model which can handle the temporal and spatial variation of density is presented. The hydrostatic approximation is used and density effects are included by means of conservation equation of heat and the equation of state. The finite difference grids are used in the horizontal plane and a set of linear-shape functions is used for the vertical expansion. The similarity transform is introduced to solve resultant matrix equations. The proposed model was first applied to the density-driven circulation in an idealized basin in the presence of the heat exchange between the air and the sea. The advection terms in the momentum equation were ignored, while the convection terms were retained in the heat equation. Coefficients of the vertical eddy viscosity and diffusivity were fixed to be constant. Calculation in a non-rotating idealized basin shows that the difference in heat capacity with depth gives rise to the horizontal gradient of temperature. Consequently, there is a steady new in the upper layer in the direction of increasing depth with compensatory counter flow .in the lower layer. With Coriolis force, geostrophic flow was predominant due to the balance between the pressure gradient and the Coriolis force. As a test in region of irregular topography, the model is applied to the Yellow Sea. Although the resultant flow was very complex, the character of the flow Showed to be geostrophic on the whole.

• Journal of agriculture & life science
• /
• v.45 no.6
• /
• pp.265-278
• /
• 2011

Effects of the open level of the side window were studied to control the temperature and relative humidity in the fog cooling greenhouse. The greenhouse was cooled by air atomizing spray nozzles of the air and water two-fluid process. The control process includes the measuring of environmental variables, setting and coding of the water balance equations and heat balance in greenhouse, calculating of the roof window open and spray water, and operating of the motor and pump. The target temperature and relative humidity were set at $28^{\circ}C$, 75%, respectively. The three modes of the side window open level were 0%, 50% and 100%. The average dry bulb temperatures of the inside air were 28.2, 27.2 and $26.3^{\circ}C$, respectively and their standard deviation was ranged from $0.4^{\circ}C$to $0.8^{\circ}C$. Also the relative humidity of the 0% mode was the best controlled one with the average of 76.3% and the standard deviation of 2.1%.

• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.921-932
• /
• 2010

In this study, the adjustment processes of the disturbed channels by sand or gravel mining were investigated by a two dimensional numerical model in the generalized coordinate system. As a numerical scheme, the CIP (cubic interpolated pseudoparticle method) method was used to calculate the advection term in the flow field and central difference method was used to the diffusion term in it. The pit of the channel was partially filled with sediment at the toe of the pit upstream. As time increased, the headcut erosion upstream in the pit was decreased due to the sediment inflow. The almost inflow sediment upstream was trapped into the pit and the sediment deposit wedge migrated downstream in the pit with the steep submerged angle of repose. The numerical model was reproduced well the evolution processes of the channel. The mining pit migrated with speed as the channel was steep, and the numerical results were in overall agreement with the experimental results.