• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6C
• /
• pp.267-274
• /
• 2012

The conventional approach to evaluate the contaminant transport in soils adopts the macro-scale implementation while the pore configuration and network is a dominant factor to determine the fate of contaminant. However, the observation of fate and transport at pore scale may not be readily approachable because of the computational expenses to solve Navier-Stokes equation. We herein present the 2D Lattice-Boltzmann method that enables to assess the local fluid velocity and density efficiently for the case of single phase and multi-components. The solute fate spatio-temperal space is explicitly determined by the advection of fluid flow. Two different types of idealized pore space provides the path of fluid. Also, solute transport, the velocity field and average concentration of solute are computed in steady state. Results show that the pore geometry such as tortuosity mainly affect the solute fate. It highlights the significance of the pore configuration and shape in granular soils and rock discontinuity in spite of the equivalent porosity.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.9 no.2
• /
• pp.251-259
• /
• 1999

A computer simulation was conducted for heat and momentum transfer in a roller kiln. Time-averaged Navier-Stokes equation conjugated with energy balance equation was numerically solved to predict the temperature distribution and fluid flow field in the roller kiln. A computer simulation was performed for a roller kiln for three cases. Firstly, when there are no ceramic materials in the roller kiln, the effect of natural convection was studied on the temperature distribution and fluid flow field. From the result, it was observed that air takes the heat of wall away from the roller kiln by natural convection and the heat was not transferred effectively. Secondly, with ceramic materials temperature difference of ceramic material from the borrom to the top of a ceramic material was about 255K in 5th zone and this is because the heat is transferred from the surface of a ceramic material to flowing air with relatively low temperature. Finally, we considered effect of radiation heat transfer. Temperature difference of ceramic material in 5th zone was about 300 K, due to radiation heat transfer on the ceramic material surfaces.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.61-70
• /
• 2016

Dong-tan Station, shared by high-speed railway and urban express railway, is a very complicated underground station having 6 tracks together with barrier and shafts between them, therefore it seems very hard to investigate the aerodynamic effects including the pressure variation and train gust in the station when a high-speed train runs through it. In this study, the aerodynamic effects on the structures and platform passengers when a high-speed train runs through an underground station have been studied using Computational Fluid Dynamics. STAR-CCM+ has been employed for numerical simulation based on Navier-Stokes equation and 2-equation turbulence model and moving mesh scheme supported by STAR-CCM+ has also been used to represent the relative motion between a train and station. Based on the simulation results, the unsteady flow fields in the underground station induced by the high-speed train have been analyzed and the pressures on the PSDs and pressure variation at the platform have quantitatively assessed.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.06a
• /
• pp.373-373
• /
• 2009

Relative humidity, membrane conductivity and water activity are critical parameters of polymer electrolyte membrane fuel cells (PEMFC) for high performance and reliability. These parameters are closely related with temperature. Moreover, the ideal values of these parameters are not always identical along the channels. Therefore, the cooling channel design and its operating condition should be well optimized along the all location of the channels. In the present study, we have performed a numerical investigation on the effects of cooling channels on performance of a PEMFC. Three-dimensional Navier-Stokes equations are solved with the energy equation including heat generated by the electrochemical reactions in the fuel cell. The present numerical model includes the gas diffusion layers (GDL) and serpentine channels for both anode and cathode gas flows, as well as cooling channels. To accurately predict the water transport across the membrane, the distribution of water content in the membrane is calculated by solving a nonlinear differential equation with a nonlinear coefficient, i.e., the water diffusivity which is a function of water content as well as temperature. Main emphasis is placed on the heat transfer between the solid bipolar plate and coolant flow. The present results show that local current density is affected by cooling channels due to the change of the oxygen concentration and the membrane conductivity as well as the water content. It is also found that the relative humidity is influenced by the generated water and the gas temperature and thus it affects the distribution of fuel concentration and the conductivity of the membrane, ultimately fuel cell performance. Unit-cell experiments are also carried out to validate the numerical models. The performance curves between the models and experiments show reasonable results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.72-78
• /
• 2008

In this study, advanced (fluid-structure interaction (FSI)) analysis system has been developed in order to predict turbine cascade performance with blade deformation effect due to aerodynamic loads. Intereference effects due to the relative movement of the rotor cascade with respect to the stator cascade are also considered. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation k-ω SST turbulence models are solved to accurately predict fluid dynamic loads considering flow separation effects. A fully implicit time marching scheme based on the (coupled Newmark time-integration method) with high artificial damping is efficiently used to compute the complex fluid-structure interaction problem. Predicted aerodynamic performance considering structural deformation effect of the blade shows somewhat different results compared to the case of rigid blade model. Cascade performance evaluations for different elastic axis positions are importantly presented and its aeroelastic effects are investigated.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1675-1683
• /
• 1995

A numerical study of the fluid flow and heat transfer characteristics of the two-dimensional impingement jet with a confinement plate has been carried out. The fluid flow was calculated by solving the full Navier-Stokes equation. In doing that, the well known SIMPLER algorithm was used and the trouble making convection term was discretized according to QUICKER scheme. The energy equation was simply solved by using the SOR method. For the Reynolds number of 10000, two channel heights, say 1.5 and 3.0 times the jet exit width, and two thermal boundary conditions constant wall temperature and constant wall heat flux were considered. Discrete heat sources were flush mounted along the impingement plate at a distance of 0, 2, 3, 4, 5, 6, 10, 12, times the jet exit width from the stagnation point. The length of each heat source is 4 times the jet exit width long. The Nusselt number averaged over each heat source was compared with experiment. Comparison shows that both calculations and experiment have the secondary peak of Nusselt number at downstream of stagnation point, even though there is a little quantitative difference in between. The difference is believed due to abscure thermal boundary condition in experiment and also accuracy of turbulence model used. The secondary peak is shown to be caused by rigorous turbulent flow motion generated as the wall jet flow is retarded and developes into the channel flow without flow reversal.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.9
• /
• pp.915-922
• /
• 2011

In this paper, the unsteady incompressible Navier-Stokes equation coupled with energy equation was solved in order to investigate the thermal stimulation of continuous moxibustion using a commercial code (ANSYS-Fluent). In the simulations, various periods were selected for the continuous moxibustion, which was done by burning multiple disks successively. It has been found that the depth of the effective stimulation zone increases only when the replacing temperature is equal or larger than body temperature whereas the increase rate of the effective stimulation depth decreases as the number of disc increases. Further, it has been shown that the optimal period, for which the duration time of the effective stimulation zone is maximum, exists.

• Journal of Korea Water Resources Association
• /
• v.38 no.3 s.152
• /
• pp.189-198
• /
• 2005

This study showed that numerical simulation can be effectively used to analyze discharge capacity according to the form and arrangement of the lock gate of a tidal power plant. For the numerical simulation, FLOW-3D with Reynolds-averaged Navier-Stokes equation as a governing equation was used. This study found that improvement of apron length and approach angle of guide wall of the lock gate causes differences in discharge capacity of $10\%$ or more. In addition, there was a difference of discharge capacity caused by the connecting structures of the drainage gate and hydraulic turbine structure and the side slope at the end of apron. This study also showed that hydraulic investigation to enhance a discharge capacity is needed when the lock gate is designed and that numerical model experiments can be a useful analysis tool to design the drainage structure, as well as the hydraulic model experiment.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.245-250
• /
• 2002

This report describes the effects of following flaws due to ship's propeller on the fish farm structure when the ship's propeller is operated in full power. This study is applied an incompressible newtonian fluid theory, which is governed the Navier-Stokes equation. For the numerical solution, Neumann equation are applied as the boundary conditions. The result shows that the flow velocity near the fish farm is 1.0 m/sec. The actual measurement carries out by using propeller type velocimeter in order to measure the velocity of following flows and currents around the fish farm area. The result shows that the maximum velocity near the fish farm structure is 1.2 m/sec in depth of 1.5 m. This velocity is used for calculation of external force on the fish farm structure. The results of structural strength of the fish farm structures show that the actual maximum bending moment and bending stress are less than the damage strength of material. So the fish farm structure is not affected by the following flows and currents of ship's propeller.