• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Journal of Korean Society for Atmospheric Environment
• /
• v.22 no.1
• /
• pp.73-84
• /
• 2006

Hazardous Air Pollutants (HAPs) are characterized by being relatively heavier and denser than that of ambient air due to the various reasons such as higher molecular weight, low temperature and other complicated chemical transformations (Witlox, 1994). In an effort to investigate transport and diffusion from instantaneous emission of heavy gas, Lagrangian Particle Dispersion Model (LPDM) coupled with the RAMS output was employed. Both deposition process and buoyancy term were added on the atmospheric diffusion equations of LPDM, and the locations and concentrations of dense gas particle released from instantaneous single point source (emitting initially for 10 minutes only) were analyzed. The result overall shows that adding deposition process and buoyancy terms on the diffusion equation of LPDM has very small but detectable effect on the vertical and horizontal distribution of Lagrangian particles that especially transported for a fairly long traveling time. Also the slumping of dense gas can be found to be ignored horizontally compared to the advection by the horizontal wind suggesting that it was essential to couple the Lagrangian particle dispersion model coupled with the RAMS model in order to explain the dispersion of HAPs more accurately. However, during the initial time of instantaneous emission, buoyancy term play an important role on the vertical locations of dense particles for near surface atmosphere and around source area, indicating the importance of densities of HAPs in the beginning stage or short duration for the risk assessment of HAPs or management of heavy vapors during the explosive accidents.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.621-628
• /
• 1997

Excavation works on stream beds have been done for various reasons including aggregate collection, sediment dredging, bridge constructions, or laying pipes under the ground. These activities may cause significant loadings of SS (suspended solids) resulting in water pollution and other detrimental effects to the surrounding environment. This research investigates application potential of a fickian diffustion model, derived from two dimensional advection-diffusion equation through some simplifying assumptions, as a planning tool for the estimation of SS loadings from excavation works and evaluation fo pollution prevention measures in case that sophisticated numerical simulation models are not applicable due to various practical reasons. Through a case study of the Juncheon stream in the Donghae City on the Kangwondo Province, this study demonstrates applicability of the fickian diffustion model as a practical method for the preliminary estimation of Ss loadings from excavation works and evaluation of performance of fabrics made of synthetic fiber for the reduction of downstream SS concentration with deficient field data.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Korean Journal of Environmental Agriculture
• /
• v.8 no.1
• /
• pp.37-46
• /
• 1989

A 2-dimensional pollutant transport phenomenon in shallow reservoirs was analyzed by using a finite element method. The Galerkin's weighted residual method, based on linear interpolation, was used and a triangle was adopted as an element. The two dimensional Stock's equation and the advection-diffusion equation integrated over depth were used as governing equations. Also the Newton-Raphson method was introduced to solve the non-linear terms of the equation. The results calculated by the model are in good agreement with the analytical solution for a simplified channel where a known solution is avaiable. An actual application of the model is attempted for Sapkyo Lake with a consideration of the influx of the Sapkyo Chun, the Muhan Chun and kogkyo Chun. Further refined research is needed to evaluate the water quality in the other reservoirs.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.450-457
• /
• 2013

For the design guide of a vessel operating in cold region, numerical analysis was carried out to evaluate the weather-tight door which installed the heating cables by using ANSYS 13.0 Transient Thermal. The numerical analysis was performed by considering Advection-Diffusion equation. This study based on the experimental results of 'A study on Anti-Icing Technique for Weather-Tight Door of Ice-Strengthened Vessels'(Jeong, et al., 2011a) in KIOST. For validation of the numerical analysis results, the cold chamber experimental data measured by the heat sensors in certain location of the weather-tight door was used. The external environmental temperature which varies from $5^{\circ}C$ to $-55^{\circ}C$ was considered in numerical analysis. Also three different heating cables which have the heat capacity of 33W/m, 45W/m and 66W/m were adapted for the design parameters to be the most efficient and guidelines for anti-icing design of the weather tight door.

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

• Journal of the Korean Society of Groundwater Environment
• /
• v.3 no.2
• /
• pp.51-59
• /
• 1996

In case waste disposal site is to be constructed close to the underground facilities such as LPG storage cavern which is completely maintained by groundwater pressure, it is generally requested that the possibility on leachate contamination of cavern area be reviewed and the countermeasure, if it is estimated cavern area is severely affected by leachate, be taken into consideration. Prediction was performed and leachate control plan was made using by analytical and the numerical analysis on the leachate migration which is likely to happen at the area between the proposed waste disposal site and the underground LPG storage cavern located at the U petrochemical complex. Analytical solutions were obtained by the conservative mass advection-diffusion equation and the effect of advection and dispersion factor on the leachate migration was reviewed through peclet number calculation and the functional relationship between the factors and leachate transport velocity was established, which leads to enable us to predict the leachate transport velocity without difficulties when different parameters (factors) are used for analytical solution. Numerical solutions were obtained by FEM using AQUA2D which is for the simulation of groundwater flow and contaminant transport. 3-D discrete fracture models were simulated and fracture flow analysis was performed and feasibility study on the water-curtain system was conducted through the fracture connectivity analysis in rock mass. As results of those analyses, it was interpreted that the leachate would trespass on the LPG storage cavern area in 30 years from the proposed wate disposal site and the vertical water-curtain system was effective mathod for the prevention of leachate's migration further into the cavern area.

• Journal of Korea Water Resources Association
• /
• v.33 no.2
• /
• pp.181-193
• /
• 2000

This paper presents the flow and the suspended sediment movement over ripples for oscillatory flows. A new numerical model system is developed, and applied to a laboratory experimental condition of regular waves and a fictitious condition of irregular waves. The flow field is obtained from a programme proposed by Kim et. al.(1994), which is a modified version of SOLA based on SMAC scheme. The sub-model solves the continuity and Reynolds momentum equations in the x-z plane. The wave orbital velocities, shear stresses, and pressure are all reasonably reproduced by the model. The model results on the vertical velocity component show good agreement with the measurements. The suspended sediment transport sub-model is newly set up to solve the advection-diffusion equation of suspended sediment using a split method, and involving a special shear entrainment from the whole ripple surface. The calculated suspended sediment concentrations for regular waves show reasonable agreement with measurements at Deltaflume. The model results for random waves show that the suspended sediment concentration is higher than those for regular waves and that the sediment diffuses higher than for regular waves with the significant wave height and the peak wave period of the irregular waves.