• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.7 no.4
• /
• pp.599-609
• /
• 1997

The phase behavior and dynamics of colloid suspensions and the resulting structures and properties of powder compacts were examined by a computer experimental method for cooperative packing processes. A wide range of properties and process conditions such as arbitrary particle size, medium densities, field strength, and temperature could be examined using the Peclet number (Pe). We demonstrated that an optimum range of Peclet number for the ordering of sediments was present and that the phenomena related to the ordering such as the onset of crystallization, the phase behavior, etc. strongly depend on process conditions. The present work appears to be useful to design the processing method of ceramic spherical submicron powders for the preparation of high-density green compacts.

• Journal of the Korean GEO-environmental Society
• /
• v.10 no.1
• /
• pp.5-14
• /
• 2009

This study is to analyze numerically the spatial behaviors of the solute transport in a spatially correlated variable-aperture fracture under the effective normal stress conditions. Numerical results show that the solute transport in a fracture is strongly affected by the spatial correlation length of apertures and applied effective normal stress. According to increasing spatial correlation length, the mean residence time of solute is decreased and the tortuosity and Peclet number (is a dimensionless number relating the rate of advection of a flow to its rate of diffusion) is also decreased. These results mean that the geometry of the aperture distribution is favorable to the solute transport as the spatial correlation length is increased. However, according to the applied effective normal stress is increased, the mean residence time and tortuosity have a tendency to increase but the Peclet number is decreased. The main reason that the Peclet number is decreased, is that the solute is displaced by one or two channels with relatively higher local flow rate due to the increment of contact areas by increasing effective normal stress. Moreover, based on numerical results of the solute transport in this study, the exponential-type correlation formulae between the mean residence time and the effective normal stress are proposed.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.179-189
• /
• 2012

To simulate the transport of pollutant, a numeric model for the advection-diffusion equation with the decay term is developed. This is finite-difference model using the implicit method (with the weight factor ${\alpha}$) and Gauss-Seidel SOR(successive over-relaxation). This model is compared to the analytical solutions (of simpler dimensional or boundary conditions), and in the condition of Peclet number < 5~20, the result shows stable condition, and Crank-Nicolson method (${\alpha}$=0.5) shows the more accurate results than fully-implicit method (${\alpha}$=1). The mass of advection, diffusion and decay is calculated and the error of mass balance is less than 3%. This model can evaluate the 3-D concentrations of the advection-diffusion and decay problems, but this model uses only the finite-difference method with the fixd grid system, so it can be effectively used in the problems with small Peclet numbers like the pollutant transport in groundwater.

• Polymer(Korea)
• /
• v.29 no.6
• /
• pp.549-556
• /
• 2005

Effects of the dimensionless variables on the heat transport phenomena in the extrusion process of a single screw extruder have been studied numerically. Based on the understanding of the solids conveying related to the geometrical structure and characteristics of the screw, the heat balance equation for the solids conveying zone was established and normalized. The finite volume method and power-law scheme were applied to derive a discretized equation and the equation was solved using the alternating direction iterative method with relaxation. Effects of the dimensionless parameters, Biot and Peclet numbers, that define the heat transfer characteristics of the solids conveying zone have been investigated with respect to the temperature of the feeding zone and the length of the solids conveying zone. As the Biot number is increased, the heat loss by cooling dominates to decrease the temperature of the barrel but it has little effects on the temperature of the solids bed and the length of the solids conveying zone. On the other hand, if the Peclet number is increased, the convection term dominates to decrease the temperature of the solids bed and it results in an increase in the length of the solids conveying zone.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• Economic and Environmental Geology
• /
• v.42 no.6
• /
• pp.609-618
• /
• 2009

High permeable faults are important geological structures for fluid flow, energy, and solute transport. Therefore, high permeable faults play an important role in the formation of hydrothermal fluid (or hot spring), high heat flow, and hydrothermal ore deposits. We conducted 2-D coupled thermal and hydraulic modeling to examine thermohydraulic behavior in fault zones with various permeabilities and geometric conditions. The results indicate discharge temperature in fault zones increases with increasing fault permeability. In addition, discharge temperature in fault zones is linearly correlated with Peclet number ($R^2=0.98$). If Peclet number is greater than 1, discharge temperature in fault zones can be higher than $32^{\circ}C$. In this case, convection is dominant against conduction for the heat transfer in fault zones.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2005.04a
• /
• pp.135-138
• /
• 2005

최근에 ELLAM 기법을 이용한 오염물 거동 문제를 많은 사람들이 다루어 오고 있다. ELLAM 기법은 기존의 Eulerian-Lagrangian 방식에서 일어나는 질량보존 문제점과 일반경계조건의 체계적인 적용 한계점을 극복하였다. 그러나 본 연구에서는 이 방식의 장단점을 네 개의 예제를 통하여 다른 모델들과 비교 검토하여 ELLAM의 수치적 고찰을 수행하고자 한다. 예제 수행 결과 Mesh Peclet Number가 무한대일때 ELLAM은 수치확산 및 수치진동과 같은 수치오차로 인해 음수의 농도 값을 갖거나 1 보다 큰 농도를 갖는 경향을 보인다. 그러나 Mesh Peclet Number 50 일때는 전체적으로 해석해와 잘 일치함을 볼 수 있다. 반면, LEZOOMPC(Lagrangian-Eulerian ZOOMing Peak and valley Capturing)는 항상 좋은 결과를 보여주고 있다. 따라서 위의 결과를 종합하여 볼 패 ELLAM의 단점은 LEZOOMPC의 성질을 이용하여 개선 및 보완될 수 있음을 간접적으로 시사해준다. 즉 LEZOOMPC에서 사용되는 선택적 국부 격자 세립화 과정을 이용하면 ELLAM에시 일어나는 다양한 수치오차를 줄일 수 있을 것이라고 판단된다.

• Journal of the Korean Society for Railway
• /
• v.12 no.3
• /
• pp.437-442
• /
• 2009

In the case of a Linear Induction Motor (LIM), numerical analysis method like Finite Element Method (FEM) has been mainly used to analyze the travelling magnetic field problem which includes the velocity-induced electromotive force. If the problem including the velocity-induced electromotive force is analyzed by FEM using the Galerkin method, the solution can be oscillated according to the Peclet Number, which is determined by conductivity, permeability, moving velocity and size of mesh. Consequently, the accuracy of the solution can be low and the vortexes of flux can be occurred at the secondary back-iron. These vortexes of the flux occurred at the secondary back-iron does not exist physically, but it can be occurred in the analysis. In this case, the vortexes of the flux can be generally removed by using Up-Wind method which is impossible to apply a conventional S/W tool (Maxwell 2D). Therefore, in this paper, authors examined the vortexes of the flux occurred at the secondary back-iron of the LIM according to variations of the Peclet Number, and analyzed whether these vortexes of the flux affect on the dynamic force characteristics of the LIM or not.

• Journal of Korea Soil Environment Society
• /
• v.5 no.1
• /
• pp.3-12
• /
• 2000

The multi-region model, to describe preferential flow, is an equation representing solute transport in soils by dividing soil into numerous pore groups and using the hydraulic properties of the soil. As the model partial differential equation (PDE) is solved numerically with finite difference methods. a modified equivalent partial differential equation(MEPDE) of the partial differential equation of the multi-region model is derived to analyze the accuracy and consistency of the solution of the model PDE and the Von Neumann method is used to analyze the stability of the finite difference scheme. The evaluation obtained from the MEPDE indicated that the finite difference scheme was found to be consistent with the model PDE and had the second order accuracy The stability analysis is performed to analyze the model PDE with the amplification ratio and the phase lag using the Von Neumann method. The amplification ratio of the finite difference scheme gave non-dissipative results with various Peclet numbers and yielded the most high values as the Peclet number was one. The phase lag showed that the frequency component of the finite difference scheme lagged the true solution. From the result of the stability analysis for the model PDE, it is analyzed that the model domain should be discretized in the range of Pe < 1.0 and Cr < 2.0 to obtain the more accurate solution.

• Journal of Korea Water Resources Association
• /
• v.31 no.5
• /
• pp.599-612
• /
• 1998

Numerical analysis of convection-dominated transport problems are challenging because of dual characteristics of the governing equation. In the finite element method, a strategy is to modify the test function to weight more in the upwind direction. This is called as the Petrov-Galerkin method. In this paper, both N+1 and N+2 Petrov-Galerkin methods are applied to transport problems at high grid Peclet number. Frequency fitting algorithm is used to obtain optimal levels of N+2 upwinding, and the results are discussed. Also, a new Petrov-Galerkin method, named as "Optimal Test Function Petrov-Galerkin Method," is proposed in this paper. The test function of this numerical method changes its shape depending upon relative strength of the convection to the diffusion. A numerical experiment is carried out to demonstrate the performance of the proposed method.