• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.76-83
• /
• 1995

Estimation of blasting behavior of explosives is prerequisite in the numerical analysis of blasting works. In this study, blasting pressure is estimated by the finite difference method using the Flux-Corrected Transport Algorithm. To formulate the behavior of blasting gas, the mass conservation equation, the moment conservation equation, the energy conservation equation and the ideal gas state equation are used. The simplified species conservation equation is included to simulate the behavior of reacting explosives. To verify the calculation, the Sod's shock tube problem, the strong shock problem and the reacting problem we used. Numerical results show that the shock wave can be captured by means of the FCT algorithm in the reacting and nonreacting states.

• Computers and Concrete
• /
• v.13 no.2
• /
• pp.149-171
• /
• 2014

This study presents the mathematical model for predicting chloride penetration into saturated concrete under non-isothermal condition. The model considers not only diffusion mechanism but also migration process of chloride ions and other chemical species in concrete pore solution such as sodium, potassium, and hydroxyl ions. The coupled multi-ionic transport in concrete is described by the Nernst-Planck equation associated with electro-neutrality condition. The coupling parameter taken into account the effect of temperature on ion diffusion obtained from available test data is proposed and explicitly incorporated in the governing equations. The coupled transport equations are solved using the finite element method. The numerical results are validated with available experimental data and the comparison shows a good agreement.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.7
• /
• pp.729-735
• /
• 2005

This paper has dealt with the characteristics of mass and heat transfer of vertical flat plate with free convection. The theory of similarity transformations applied to the momentum and energy equations for free convection. To derive the similarity equation of mass transfer. the equation for conservation of species was added to the continuity. momentum and energy equations. The momentum, energy and species equations set numerically to obtain the velocity, temperature and mass fraction of species as dimensionless. For cases where momentum transport dominates, the thermal boundary layers are shorter than the momentum boundary layer. The relationships between momentum, energy and species were clarified from this study.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.6
• /
• pp.609-617
• /
• 2012

The impact, spreading and recoil processes of a nanoparticle-laden droplet impacting on a horizontal solid surface are numerically investigated by solving the conservation equations for mass, momentum, energy and mass fraction. The liquid-air interface is tracked using a level-set method that is modified to include the effect of contact angle hysteresis at the wall. The species transport equation including a thermal diffusion term is additionaly solved to determine the nanoparticle distribution in the droplet. The effect of nanoparticle concentration and contact angle are also studied.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.9-16
• /
• 2013

The heat transfer enhancement of heat sink with mist flow is studied numerically by solving the conservation equations for mass, momentum and energy in the continuous and dispersed phases. A Lagrangian method is used for tracing dispersed water droplets in the heat sink and an Eulerian species transport model for air and steam mixture. The continuous and dispersed phases are interacted with the drag and evaporation source terms. The computed results show that addition of evaporating mist droplets enhances the cooling performance of heat sink significantly.

• Computers and Concrete
• /
• v.11 no.3
• /
• pp.201-222
• /
• 2013

A comprehensive simulation model for the transport process of fully coupled moisture and multi-species in non-saturated concrete structures is proposed. The governing equations of moisture and ion diffusion are formulated based on Fick's law and the Nernst-Planck equation, respectively. The governing equations are modified by explicitly including the coupling terms corresponding to the coupled mechanisms. The ionic interaction-induced electrostatic potential is described by electroneutrality condition. The model takes into account the two-way coupled effect of moisture diffusion and ion transport in concrete. The coupling parameters are evaluated based on the available experimental data and incorporated in the governing equations. Differing from previous researches, the material parameters related to moisture diffusion and ion transport in concrete are considered not to be constant numbers and characterized by the material models that account for the concrete mix design parameters and age of concrete. Then, the material models are included in the numerical analysis and the governing equations are solved by using finite element method. The numerical results obtained from the present model agree very well with available test data. Thus, the model can predict satisfactorily the ingress of deicing salts into non-saturated concrete.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.320-327
• /
• 2006

The computations of chemically reacting laminar and turbulent flows are performed using the preconditioned Navier-Stokes solver coupled with turbulent transport and multi-species equations. A low-Reynolds number $k-\varepsilon$ turbulence model proposed by Chien is used. The presence of the turbulent kinetic energy tenn in the momentum equation can materially affect the overall stability of the fluids-turbulence system. Because of this coupling effect, a fully coupled formulation is desirable and this approach is taken in the present study. Choi and Merkle's preconditioning technique is used to overcome the convergence difficulties occurred at low speed flows. The numerical scheme used for the present study is based on the implicit upwind ADI algorithm and is validated through the comparisons of computational and experimental results for laminar methane-air diffusion flame and $ H_2/O_2$ reacting turbulent shear flow. Preconditioning formulation shows better convergence characteristics than that of non-preconditioned system by approximately five times as much.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.10
• /
• pp.1279-1287
• /
• 2004

The vaporization characteristics of a liquid heptane droplet in a supercritical nitrogen flow are numerically studied. The transient conservation equations of mass, momentum, energy, and species are expressed in an axisymmetric coordinate system. The governing equations are solved time marching method with preconditioning scheme. The modified Soave-Redlich-Kwong equation of state is employed for taking account of real gas effects such as thermodynamic non-ideality and transport anomaly. Changing the convective velocity and ambient pressure, several parametric studies are conducted. The numerical results show that the two parameters, Reynolds number and dimensionless combined parameter(${\mu}$s/${\mu}$d)(equation omitted), have influence on supercritical droplet vaporization.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.

• Computers and Concrete
• /
• v.3 no.6
• /
• pp.401-422
• /
• 2006

A multi-species model based on the Nernst-Planck equation has been developed by using a finite volume method. The model makes it possible to simulate transport due to an electrical field or by diffusion and to predict chloride penetration through water saturated concrete. The model is used in this paper to assess and analyse chloride diffusion coefficients and chloride binding isotherms. The experimental assessment of the effective chloride diffusion coefficient consists in measuring the chloride penetration depth by using a colorimetric method. The effective diffusion coefficient determined numerically allows to correctly reproduce the chloride penetration depth measured experimentally. Then, a new approach for the determination of chloride binding, based on non-steady state diffusion tests, is proposed. The binding isotherm is identified by a numerical inverse method from a single experimental total chloride concentration profile obtained at a given exposure time and from Freundlich's formula. In order to determine the initial pore solution composition (required as initial conditions for the model), the method of Taylor that describes the release of alkalis from cement and alkali sorption by the hydration products is used here. Finally, with these input data, prediction of total and water-soluble chloride concentration profiles has been performed. The method is validated by comparing the results of numerical simulations to experimental results obtained on various types of concretes and under different exposure conditions.