• Journal of computational fluids engineering
• /
• v.21 no.1
• /
• pp.43-49
• /
• 2016

Thermal cracking is commonly modeled as plug flow reaction, neglecting the lateral gradients present. In this paper, 2-dimensional computational fluid dynamics including turbulence model and molecular reaction scheme are carried out. This simulation is solved by means of coupled implicit scheme for stable convergence of solution. The reactor is modeled as an isothermal tube, whose length is 1.2 m and radius is 0.01 m, respectively. At first, The radial profile of velocity and temperature at each point are predicted in its condition. Then the bulk temperature and conversion curve along the axial direction are compared with other published data to identify the reason why discussed variations of properties are important to product yield. Finally, defining a new non-dimensional number, Effect of interaction with turbulence, heat transfer and chemical reaction are discussed for design of thermal cracking furnace.

• Journal of Korea Water Resources Association
• /
• v.30 no.2
• /
• pp.143-154
• /
• 1997

A numerical model for solving advection-diffusion equation is presented by splitoperator method combining the Holly-Preissmann scheme with a fifth-degree interpolating polynomial for advection operator and the explicit scheme porposed by Hobson et al. for diffusion operator. To examine the developed model, the obtained numerical solutions are compared with both the analytic solution and those from the existing models for the instantaneous source (Gaussian hill) and the continuous source (advanced front) at upstream boundary with constant velocity and diffusivity condition. For the various cases having different Courant and Peclet numbers, it is shown that the present study provides stable solutions even for Courant numbers exceeding one. The result obtained by the present study also agree well with existing analytical solutions for both cases. The proposed explicit scheme somewhat releases the conventional restriction of explicit schemes for determining the time step size and provides satisfactory results for relatively large time step size.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.88-88
• /
• 2020

수심평균 2차원 혼합모형은 하천환경에서 다양한 용존성 오염물질의 혼합현상을 모의하기 위해 널리 활용되어왔다. 2차원 혼합모형에서 분산계수는 하천의 전단 흐름에 의해 야기되는 오염물질의 퍼짐 현상을 표현하는 중요한 인자로서 작용하기 때문에 정교한 오염물질 혼합거동을 모의하기 위해서는 적합한 분산계수를 산정하는 것이 필수적이다. 분산계수를 실험적으로 산정하는 방법으로는 크게 모멘트법과 추적법으로 나뉘며, 비정상상태의 혼합거동에 대해 종방향 및 횡방향 분산계수를 동시에 산정할 수 있는 방법은 추적법 계열의 2차원 유관추적법(2D STRP)이 유일하다. 본 연구에서는 하천에 유입된 오염물질의 2차원 혼합해석을 위한 수치모형을 개발하였으며, 개발된 모형의 수치해를 바탕으로 다양한 Peclet 수의 범위에 대해 기존연구에서 제시된 2D STRP의 적용범위 및 성능을 정량적으로 분석하였다. 분석된 정보를 바탕으로 기존 2D STRP의 한계를 극복하기 위한 개선된 2차원 유관추적법(2D STRP-i)을 개발하고, 사행하천을 모형화한 실규모 하천실험시설에서 검증하였다. 기존 2D STRP의 성능평가 결과, Peclet 수가 낮은 조건일수록 농도분포의 예측 정확도가 감소하는 경향을 보였으며, 하안 경계에 도달하는 농도가 증가할수록 부정확한 결과를 초래하는 것으로 나타났다. 본 연구에서는 기존 2D STRP의 한계를 보완하여 더욱 정확한 분산계수를 산정하고자 하안 경계면 조건을 고려한 2차원 유관추적법(2D STRP-i)을 개발하였다. 2D STRP-i는 직교-곡선좌표계 기반의 2차원 이송-분산 방정식을 바탕으로 횡방향 유속분포 및 하안 경계조건을 고려할 수 있도록 개선되었다. 2D STRP-i는 공간적으로 상이한 이송효과 및 하안경계 조건을 적절히 반영함으로써 농도분포의 예측 정확도를 개선 시키는 것으로 평가되었으며, 하안경계면에서 농도가 증가하는 구간에서 기존 2D STRP의 결과와 비교하여 더욱 정확한 농도분포 및 분산계수를 제공하는 것으로 밝혀졌다.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.13 no.2
• /
• pp.64-76
• /
• 2009

Experimental and numerical studies are conducted on the characteristics of flame extinction and edge flame oscillation in counterflow diffusion flames. The characteristics of flame extinction and edge flame oscillation are well described varying burner diameter, separation distance between two burners, global strain rate, and velocity ratio. It is verified numerically and experimentally that radial conduction heat loss significantly contributes to flame extinction and edge flame oscillation at low strain rate flames in zero- and micro-gravity. It is also shown that for appropriately small burner diameters flame extinction modes are grouped into four and these are significantly attributed to excessive radial conduction heat loss. The edge flame oscillation can be characterized well by one curve with Strouhal number and Peclet number.

• Journal of Energy Engineering
• /
• v.10 no.2
• /
• pp.161-165
• /
• 2001

봉다발 내 온도장 해석을 위해 개발되어진 난류 Prand시 수 모델들을 중심으로 액체금속에 대한 비교연구를 수행하였다. VANTACY-II 코드에 사용된 Zeggel & Monir의 모델의 기초가 된 Jischa & Rieke 모델 및 상수형 모델(Pr$_{t}$=0.9)을 비교대상 모델로 선정하여 난류 Prandtl 수의 비등방성과 공간분포 및 분자 Prandtl 수의 영향을 고려한 본 연구모델과 P/D와 Peclet 수를 변화시키며 얻어진 Nusselt 수의 결과를 비교하였다. 비교결과 본 모델이 다른 모델에 비하여 봉다발 내 액체금속의 열전달 거동을 전반적으로 잘 예측하였다.

• Resources Recycling
• /
• v.30 no.4
• /
• pp.20-26
• /
• 2021

Solar grade silicon (SoG-Si) has been commercially supplied mainly from off-grade high-purity silicon manufactured for electronic-grade Si (EG-Si). Therefore, for wider application of solar cells, the development of a refining process at a considerably lower cost is required. The most cost-effective and direct approach for producing SoG-Si is to purify and upgrade metallurgical-grade Si (MG-Si). In this study, directional solidification of molten MG-Si was conducted in a high-frequency induction furnace to remove iron from molten Si. The experimental conditions and results were also discussed with respect to the effective segregation coefficient, Scheil equation, and Peclet number. The study showed that when the descent velocity of the specimen decreased, the macro segregations of impurities and ingot purities increased. These results were derived from the decrease in the effective segregation coefficient with the decrease in the rate of descent of the specimen.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.6
• /
• pp.28-34
• /
• 2006

Laboratory column tests using $Cl^-$ tracer were conducted to study the correlation of soil particle distribution and hydrodynamic dispersion mechanism with three kinds of ununiformed soil samples, in which the ratio of gravel and sand versus silt and clay is 24.5 for S-1 soil, 4.48 for S-2 soil, and 0.4 for S-3 soil. Chloride breakthrough curves with time were fitted with gaussian functions. The relative concentrations of chloride were converged to 1.0 after 0.7 hours for S-1, 6.3 hours for S-2, and 389 hours for S-3. Average linear velocity, longitudinal dispersion coefficient, and longitudinal dispersivity were calculated by chloride breakthrough curves. Longitudinal dispersion coefficients were $1.20{\times}10^{-4}\;m^2/sec$ for S-1, $8.87{\times}10^{-7}\;m^2/sec$ for S-2, and $1.94{\times}10^{-9}\;m^2/sec$ for S-3. Peclet numbers calculated by the molecular diffusion coefficient of chloride and the mean grain diameters of soils were $2.59{\times}10^2$ for S-1, $6.27{\times}10^0$ for S-2, and $1.35{\times}10^{-4}$ for S-3. Mechanical dispersion was dominant for the hydrodynamic dispersion mechanism of S-1. Both mechanical dispersion and molecular diffusion were dominant for the hydrodynamic dispersion mechanism of S-2, but mechanical dispersion was ascendant over molecular diffusion. Hydrodynamic dispersion in S-3 was occurred mainly by molecular diffusion. When plotting three soils on the graph of $D_L/D_m$ versus Peclet number produced by Bijeljic and Blunt (2006), the values of $D_L/D_m$ for S-1 and S-2 were more than 2.0 order compared to their graph. S-3 was not plotted on their graph because the Peclet number was as small as $1.35{\times}10^{-4}$.

• 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 Soil Environment Society
• /
• v.2 no.2
• /
• pp.81-94
• /
• 1997

In many solute transport studies, either flux or resident concentration has been used. Choice of the concentration mode was dependent on the monitoring device in solute displacement experiments. It has been accepted that no priority exists in the selection of concentration mode in the study of solute transport. It would be questionable, however, to accept the equivalency in the solute transport parameters between flux and resident concentrations in structured soils exhibiting preferential movement of solute. In this study, we investigate how they differ in the monitored breakthrough curves (BTCs) and transport parameters for a given boundary and flow condition by performing solute displacement experiments on a number of undisturbed soil columns. Both flux and resident concentrations have been simultaneously obtained by monitoring the effluent and resistance of the horizontally-positioned TDR probes. Two different solute transport models namely, convection-dispersion equation (CDE) and convective lognormal transfer function (CLT) models, were fitted to the observed breakthrough data in order to quantify the difference between two concentration modes. The study reveals that soil columns having relatively high flux densities exhibited great differences in the degree of peak concentration and travel time of peak between flux and resident concentrations. The peak concentration in flux mode was several times higher than that in resident one. Accordingly, the estimated parameters of flux mode differed greatly from those of resident mode and the difference was more pronounced in CDE than CLT model. Especially in CDE model, the parameters of flux mode were much higher than those of resident mode. This was mainly due to the bypassing of solute through soil macropores and failure of the equilibrium CDE model to adequate description of solute transport in studied soils. In the domain of the relationship between the ratio of hydrodynamic dispersion to molecular diffusion and the peclet number, both concentrations fall on a zone of predominant mechanical dispersion. However, it appears that more molecular diffusion contributes to the solute spreading in the matrix region than the macropore region due to the nonliearity present in the pore water velocity and dispersion coefficient relationship.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.7
• /
• pp.975-980
• /
• 2001

The thermal analysis of the hot roller in a dry film laminator is studied numerically by steady-state two-dimensional heat transfer. In the laminating process for PDP glass or PCB, the temperature distributions in a hot roller are presented considering the effects of the roller rotation speed and the inner and outer radii of the roller. The results show that the temperature distributions are strongly dependent on Peclet number. If Pe number becomes larger, the iso-thermal lines are more concentric about the rotating axis and the temperature difference on the hot roller surface decreases exponentially. It also shows that if the contact angle between the roller and the film becomes smaller the temperature difference becomes smaller. However, the changes of the rollers inner or outer radius have little effect on the temperature difference.