• Journal of the Korean Society of Visualization
• /
• v.20 no.3
• /
• pp.74-85
• /
• 2022

We propose a Cubic-Interpolated Pseudo-Particle Lattice Boltzmann method (CIP-LBM) for the convection-diffusion equation (CDE) based on the Bhatnagar-Gross-Krook (BGK) scheme equation. The CIP-LBM relies on an accurate numerical lattice equilibrium particle distribution function on the advection term and the use of a splitting technique to solve the Lattice Boltzmann equation. Different schemes of lattice spaces such as D1Q3, D2Q5, and D2Q9 have been used for simulating a variety of problems described by the CDE. All simulations were carried out using the BGK model, although another LB scheme based on a collision term like two-relation time or multi-relaxation time can be easily applied. To show quantitative agreement, the results of the proposed model are compared with an analytical solution.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.404-405
• /
• 2011

원전구조물의 방사성 폐기물 처분시설의 경우 지하수에 해수가 유입되어 콘크리트에 염소이온 침투가 발생할 수 있으며, 콘크리트 내부에 존재하는 인장철근의 부식에 의한 내구성 저하 및 수명 단축이 주된 문제가 된다. 본 논문에서는 원전콘크리트 구조물에서의 확산에 의한 염소이온 침투에 대한 수학적 모델을 제시하였다. 콘크리트 중의 염소이온의 침투는 콘크리트의 노출환경, 습윤상태에 따라 확산(Diffusion), 대류(Absorption), 전기적 이동(Migration)에 의해 발생한다. 이러한 조건을 모두 고려하여 제시한 방정식에 의해 염소이온의 침투를 예측할 수 있다.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.19-30
• /
• 2000

A finite-volume method is developed to solve turbulent flows around modern commercial hull forms with bow and stern bulbs. The RANS equations are solved. The cell-centered finite-volume method employs QUICK and central difference scheme for convective and diffusive flux discretization, respectively. The SIMPLEC method is adopted for the velocity-pressure coupling. The developed numerical methods are applied to calculate turbulent flow around KRISO 3600TEU container ship. Surface meshes are generated into five blocks: bow and stern bulbs, overhang, fore and afterbody. 3-D field grid system with O-H topology is generated using elliptic grid generation method. Surface friction lines and wake distribution at propeller plane is compared with experiment. The calculated results show that the present method can be used to predict flow around a modern commercial hull forms with bulbs.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.3
• /
• pp.528-536
• /
• 1987

In a turbulent free shear flow, the large scale motion is characterized by the intermittent flow which arises from the interaction between the turbulent fluid and the irrotational fluid of the environment through the mean velocity gradient. This large scale motion causes a bulk convection whose effect is similar to the spatial diffusion process. In this paper, the total diffusion process is proposed to be approximated by weighted sum of the bulk convection due to the large scale motion and the usual gradient diffusion due to small scale motion. The diffusion term in conventional .kappa.-.epsilon. model requires on more equation of the intermittency transport equation. A production term of this equation means mass entrainment from the irrotational fluid to the turbulent one. In order to test the validity of the proposed model, a plane jet is predicted by this method. Numerical results of this model is found to yield better agreement with experiment than the standard .kappa.-.epsilon. model and Byggstoyl & Kollmann's model(1986). Present hybrid diffusion model requires further tests for the check of universality of model and for the model constant fix.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.6
• /
• pp.772-778
• /
• 2002

A numerical analysis of turbulent gas-particle two-phase flow is performed in conjunction with the experiments of Fackrell ＆ Robins and Raupach & Legg that considered ground-level source and/or elevated source flat plate flow. K-$\omega$ turbulence model is used in order to analyze fully turbulent flow field and the concentration equation with settling velocity is adopted for the concentration field. The model of Einstein and Chien is applied that couples the velocity field and the concentration field. Turbulent eddy viscosity is re-evaluated in this model. The present numerical results have good agreement between the simulation and the experimental data for the mean flow velocities and particle concentrations. While the previous study shows about 27% error in the vicinity of the source of particle concentration, the .present study allows about 14% error. A new turbulent gas-particle flow model developed by this study is able to cut down error by 13% at a near source.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.162-171
• /
• 1993

This paper deals with the anisotropy of the mushy region during solidification process of a binary mixture. A theoretical model which specifies a permeability tensor in terms of pricipal values is proposed. Also, the governing equations are modified into convenient forms for the numerical analysis with the existing algorithm. Some test computations are performed for soeidification of aqueous ammonium chloride solution contained in a square cavity. Results show that not only the present model is capable of resolving fundamental characteristics of the tranport phenomena, but also the anisotropy significantly affects the interdendritic flow structure, i.e., double-diffusive convection and macrosegregation patterns.

• Fire Science and Engineering
• /
• v.10 no.3
• /
• pp.10-18
• /
• 1996

In foreign country such as USA and Japan, considerable research has been done regarding the spread of smoke in room of fire involvement by using computer. but, in our country it has not been, so, this paper presents a detailed qualitative description of phenomena which occure during typical fire scenarios through numberical analysis. The governing equations are solved by using FVM method with non-staggered grid. The SIMPLE method for pressure-velocity couple and power-law scheme for convention terms are used. It shows that a plume is formed, hot plume gases impinge on the ceiling and they spread across it. then, it eventually reaches the bounding walls of the enclosure. It takes 20s for smoke to fill the enclosure.

• Journal of Radiation Protection and Research
• /
• v.17 no.2
• /
• pp.49-59
• /
• 1992

A new model was developed in order to investigate the effects of chelating agents on the migration of a radionuclide in the form of ion or chelate. The migration behavior of the chelated radionuclide was analyzed by formulating a convective-dispersion transport equation which included a degradation of chelating agent and chelated radionuclide. The mathematical model was analytically solved and checked with the existing retardation factor. The results show that the migration velocity of the chelated radionuclide was much faster than the ionic one due to the decreased retardation. Therefore, it was concluded that a new remedial action should be developed to reduce the generation and release of chelating agents from the nuclear power plant into the environment.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.29 no.6
• /
• pp.365-371
• /
• 2019

In this study, the flame temperature distribution of the diffusion flame burner for SiO₂ deposition was analyzed by the computational fluid analysis. This corresponds to the previous step for simulating the SiO₂ preform deposition process for manufacturing optical fibers using environmentally friendly raw materials. In order to model premixed combustion, heat flow, convection, and chemical reactions were considered, and Reynolds-averaged Navier-Stokes equations and k-ω models were used. As a result, the temperature distribution of the flame showed a tendency to increase the distance from the nozzle surface to the maximum temperature when the flow rate of the auxiliary oxygen increased. In addition, it was confirmed that the temperature distribution due to incomplete combustion was large in the combustion reaction with a large equivalence ratio of the mixed gas.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.290-297
• /
• 2009

A two-phase (gas and liquid) flow analysis solver, named CUPID, has been developed for a realistic simulation of transient two-phase flows in light water nuclear reactor components. In the CUPID solver, a two-fluid three-field model is adopted and the governing equations are solved on unstructured grids for flow analyses in complicated geometries. For the numerical solution scheme, the semi-implicit method of the RELAP5 code, which has been proved to be very stable and accurate for most practical applications of nuclear thermal hydraulics, was used with some modifications for an application to unstructured non-staggered grids. This paper is concerned with the effects of interpolation schemes on the simulation of two-phase flows. In order to stabilize a numerical solution and assure a high numerical accuracy, the second-order upwind scheme is implemented into the CUPID code in the present paper. Some numerical tests have been performed with the implemented scheme and the comparison results between the second-order and first-order upwind schemes are introduced in the present paper. The comparison results among the two interpolation schemes and either the exact solutions or the mesh convergence studies showed the reduced numerical diffusion with the second order scheme.