• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.245-255
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume, $$lt;S_{ij}'(t',l')$gt;$ and plume centroid variances, $$lt;R_{ij}'(t',l')$gt;$ were simulated with 3200 Monte Carlo runs for three variances of log K, $\omega^2_y1.0,,2.5,$ and 5.0, and three dimensionless lengths of line plume sources ( l=,5 and 10) normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are not fit well with the first order results. The first order theoretical results definitely underestimated the simulated transverse second spatial moments for the aquifers of large u: and small initial plume sources. The ergodic condition for the second spatial moments is far from reaching, and the first order theoretical results of the transverse second spatial moment of the ergodic plume slightly underestimated the simulated moments.

• The Journal of Engineering Geology
• /
• v.5 no.3
• /
• pp.289-300
• /
• 1995

At the initial stage of the underground reservoir design one should thoroughly consider surface and subsurface hydrology, hydrogeologic characteristics of aquifer system, and the function of cut - off wall because it is linked to the effective management. In this study, three dimensional finite difference model was applied to analyse the function of Ian underground reservoir at Kyungbuk Province. The steady and unsteady state conditions after construction of the underground dam were simulated through the model, and from these results the groundwater budget and the safe yield were determined. The model simulation indicates the infiltration of irrigation water to be one of the major factors of seasonal fluctuation of groundwater level. The recharge rates of irrigation water were estimated as 4.3mm/d during May and June, and 1.7mm/d during July and Agust. Groundwater recharge from the watershed area estimated to about $0.04m^3/s$, almost consistent through the year. In 1984, groundwater discharge through the transverse section of the dam was $0.002m^3/s$ and the optimum yield for two momths(July and Aguest)was $254000m^3$, however, the discharge became $0.013m^3/s$ in1993, implying the failure of cut -off function. without appropaiate of the cut - off wall, optiumum yield during the irrigaton period would be $93, 000m^3$.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.224-232
• /
• 1999

Latent heat storage system using micro-encapsuled phase change material is effective method for floor heating of house and building. The temperature profile in capsule block and flow rate of hot water are important parameters for the development of heat storage system. In the present study, a mathematical model based on 3-D, non-steady state, Navier-Stokes equations, scalar conservation equations and turbulence model ($\kappa$-$\varepsilon$), is used to predict the temperature profiles in capsule and the velocity vectors in hot water pipe. The multi-block grids and fine grids embedding are used to join the circle in hot water pipe and square in capsule block. The phase change process of the capsule is quite complex not only because the size of phase change material is very small, but also because phase change material is mixed with the cement to form thermal storage block. In calculation, it's assumed that the phenomena of phase change is limited only the thermal properties of phase change material and the change of boundary is not happened in capsule. The purpose of this study is to calculate the temperature profiles in capsule block and velocity vectors in hot water pipe using the numerical calculation. Two kinds of thermal boundary condition were considered, the first (case 1) is the adiabatic condition for the both outside surfaces of the wall, the second (case 2) is the case in which one surface is natural convection with atmosphere and another surface is adaibatic. Calculation results are shown that the temperature profile in capsule block for case 1 is higher than that for case 2 due to less heat loss in adaibatic surface. Specially, in the domain of near Y=0, the difference of temperature is greater in case 1 than in case 2. The detailed experimental data of capsule block on the temperature profile and the thermal properties such as specific heat and coefficient of heat transfer with the various temperature are required to predict more exact phenomena of heat transfer.

• Journal of Korea Soil Environment Society
• /
• v.4 no.3
• /
• pp.45-56
• /
• 1999

A biofilter was established to remove the ammonia, which is representative nitrogen-contained malodorous gas. in a compost factory. Removal efficiency of ammonia and hydrogen sulfide also was investigated. A quantity of malodor gas produced in a compost factory was affected greatly by the weather. compost states and working condition of a fertilizing mixer, and the produced gas concentrations doubled by above various parameters. By operating a water scrubbing system for removing water-soluble malodorous gases effectively. we could improve the removal efficiency over three times. We investigated long-term stability of biofilter under continuous gas flow(SV=500h-1) for 100 days. The results showed 30 days of microbial retention time. After the days, deodorization efficiency of biofilter was kept steady state. and the removal efficiency was kept over 95% for ammonia and 97% for hydrogen so]fide. respectively. The electric consumption of the biofilter, which could treat malodorous gas of 100$\textrm{m}^3$/min, applied in the compost factory was evaluated about 80u0day and water consumption was 80~100$\ell$/day. These results concluded that the biofilter is a excellent deodorization technology as well as cost-effective for removing malodorous gas produced in a compost factory.

• Tunnel and Underground Space
• /
• v.11 no.4
• /
• pp.319-327
• /
• 2001

Introduction of new design tools has been required to optimally design and operate the ventilation system of long vehicle tunnels.. The demand has led to wide spread use of the simulation technique throughout the would to analysis the dynamic relationship among the variables associated with vehicle tunnel ventilation. This paper aims at performing on-site study at local tunnels to test the applicability of NETVEN, a simulation model vehicle tunnel ventilation. The study was carried out at four urban as well as highway tunnels model of vehicle tunnel ventilation. The study was carried out at four urban as well as highway tunnels employing different ventilation systems as well as traffic methods. There were some discrepancies sound between the simulation output and measurements and the following four factors are considered to mainly cause those disagreement. (1) The real situation shows distinctive transient and retarding characteristics with respect to air flow and contaminant dispersion, while ventilation forces are not steady-state and in particular those traffic and climatic variables show significant instantaneous variation. (3) Near the exit portal, the CO levels show bigger differences. The general trend is that data with higher CO concentrations carry bigger discrepancies. Turbulent diffusion is though to be the main reason for it and also contribute to the fact hat the highest CO concentrations are found at the locations somewhat inward, not at the exit portals. (4) Higher traffic rate results in higher discrepancies of ventilation velocity. Along with the exhaust characteristics, the vehicle aerodynamic characteristics need to be studied continuously in order to reduce the velocity disagreement.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Journal of Korean Ophthalmic Optics Society
• /
• v.14 no.4
• /
• pp.33-37
• /
• 2009

Purpose: To examine the coloring condition of titanium spectacle frames with various colors by using anodization method. Methods: We made an anodization coater. Platinum plate with $3{\times}3cm^2$ was used for a cathode and titanium spectacle frame specimens was mounted on an anode in an electrolyte. An electric source device were designed to supply steady state current. The color of the coated spectacle frame specimens were measured by a spectrophotometer equipped with an integrating sphere. We use CIE $L^*a^*b$ color system as chromaticity coordinates. Results: The thickness of $TiO_2$ of titanium spectacle frame specimens was varied as controlling current flow time for electrodes. The specimens with various kinds of color as a walnut, a yellow brown, a navy blue, a blue, a light blue, a mung bean, a yellowish green, a light purple, a purple, a flower pink, a bluish green, an emerald, and a green color etc. were obtained. The values of CIE $L^*a^*b^*$ for these specimens were measured and analyzed to be changed clockwise in chromaticity coordinates as the thickness of $TiO_2$ increases. Conclusions: We identified the coloring mechanism by anodization method in titanium spectacle frame specimens.

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

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

• Journal of Korean Society of Water and Wastewater
• /
• v.9 no.3
• /
• pp.116-126
• /
• 1995

The purpose of this study was to investigate the removal of ammonium nitrogen by biological nitrification in raw water containing LAS using BAC. At batch teats, LAS removal by ozone followed the first order reaction, and the rate constants(k) by ozone dose 1, 3mg/min.L were $0.040min^{-1}$, $0.062min^{-1}$ respectively. Therefore, the more ozone was dosed, the higher LAS was removed The reaction between ozone and ammonium nitrogen also followed the first order, and rate constants(k) at pH7,8 and 9 were $8.9{\times}10^{-4}min-1$, $3.8{\times}10^{-3}min^{-1}$, and $2.9{\times}10^{-2}min^{-1}$ respectively at ozone dose of 3mg/min.L . Therefore, ammonium nitrogen was little removed by ozone under neutral pH of 7. The continuous flow apparatus had four sets composed of a ozone contacter and a GAC column. Through continuous filtration test for 50days, the following conclusions were derived; (1) LAS was removed 23%, 30% respectively by ozone dose 1, 3mg/L, and was not detected in all column effluents during the period of experiment. Therefore, it appeared that adsorption capacities of each column still remained. (2) Ammonium nitrogen concentration after ozone contact varied little in raw Water because pH of raw water was from 6 to 7, and was transfered to nitrite and nitrate within GAC columns as the result of staged nitrification. After 30days, nitrite was not detected in all column effluents due to biological equilbrium between nitro semonas and nitrobacter Average removals of ammonium nitrogen in each column after the lapse of 30days were the following; ${\cdot}$ column A (ozone dose 3mg/L, EBCT 9.5min): about 100% ${\cdot}$ column B (ozone dose 1mg/L, EBCT 9.5min): 91% ${\cdot}$ column C (ozone dose 3mg/L, EBCT 14.2min): about 100% ${\cdot}$ column D (ozone dose 0mg/L, EBCT 9.5min): 53% Though column A and C reached nitrification of about 100%, column C (longer EBCT than column A) was more stable than column A. (3) After backwash, nitrification reached steady state within 5 to 8 hours. Therefore, nitrification was not greatly affected by backwash. (4) According to the nitrification capacity in depth of column A, C, where 100% nitrification occured. LAS was removed within 20cm, while ammonium nitrogen required more depth to be removed by nitrification.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.15 no.2
• /
• pp.91-98
• /
• 2009

The long-term water-quality change of Asan Bay by the influx of polluted disposal water was predicted through a simulation with an Eco-hydrodynamic model. Eco-hydrodynamic model is composed of a multi-level hydrodynamic model to simulate the water flow and an ecosystem model to simulate water quality. The water quality simulation revealed that the COD(Chemical Oxygen Demand), dissolved inorganic nitrogen(DIN) and dissolved inorganic phosphorus(DIP) are increased at 5 stations for the subsequent 6 months after the influx of the effluent. COD, DIN and DIP showed gradual decreases in concentration during the period of one to two years after the increase of last 6 months and reached steady state for next three to ten years. Concentration levels of COD, DIN, and DIP showed the increase by the ranges of $11{\sim}67%$, $10{\sim}67%$, and $0.5{\sim}7%$, respectively, which represents that the COD and DIN are the most prevalent pollutants among substances in the effluent through the sewage treatment plant. The current water quality of Asan Bay based on the observed COD, TN and TP concentrations ranks into the class II of the Korean standards for marine water quality but the water quality would deteriorate into class III in case that the disposal water by the sewage plant is discharged into the Bay.