• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.3
• /
• pp.201-215
• /
• 2018

In this paper we consider a class of singularly perturbed system of delay differential equations of convection diffusion type with integral boundary conditions. A finite difference scheme on an appropriate piecewise Shishkin type mesh is suggested to solve the problem. We prove that the method is of almost first order convergent. An error estimate is derived in the discrete maximum norm. Numerical experiments support our theoretical results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• Biotechnology and Bioprocess Engineering:BBE
• /
• v.5 no.5
• /
• pp.332-339
• /
• 2000

Experimental and theoretical works were performed for the separation of large polyelectrolytes such as DNA in the column packed with gel particles under an electric field. This paper shows how intraparticle convection effects the separation of DNAs in the column because DNAs quickly oriented through the pores in the field direction. Dimensionless transient mass balance equations were derived considering diffusion and electrophoretic convection. The separation criteria is theoretically studied using two different Peclet numbers in the fluid and solid phases and these criteria were verified uing two different DNAs by electrophoretic mobilities measured experimentally, showing how the separation position of DNAs varies in the column according to values of Pe(sub)f/Pe(sub)g of individual DNA. Governing equations are simultaneously solved by operator theoretic and characteristic methods to yield the column response.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.123-129
• /
• 1998

A numerical study has been conducted to investigate the effect of shock waves on the mixing and the recirculation zone of a hydrogen jet diffusion flame in a supersonic combustor. The general trends are compared with the experimental results obtained from the supersonic combustor at the University of Michigan. For the numerical simulation of supersonic diffusion flames, multi-species Navier-Stokes equations and detailed chemistry reaction equations of $H_2$-Air are considered. The $K-{\omega}/k-{\varepsilon}$ blended two equation turbulent model is used. Roe's FDS method and MUSCL method are used for convection fluxes in governing equations. Numerical results show that when slender wedges are mounted at the combustor wall the mixing and the combustion are enhanced and the size of recirculation zone is increased . The flame shape of supersonic flames is different in the flame-tip; it is not closed but open. The flame shape is shown to be greatly affected by shock waves.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.12
• /
• pp.1774-1783
• /
• 1998

The transition from steady laminar to chaotic convection in a glass melting furnace specified by upper surface temperature distribution has been studied by the direct numerical analysis of the two and three-dimensional time dependent Navier-Stokes equations. The thermal instability of convection roll may take place when modified Rayleigh number($Ra_m$) is larger than $9.71{\times}10^4$. It is shown that the basic flows in a glass melting furnace are steady laminar, unsteady periodic, quasi-periodic or chaotic flow. The dimensionless time scale of unsteady period is about the viscous diffusion time, ${\tau}_d=H^2/{\nu}_0$. Through primary and secondary instability analyses the fundamental unsteady feature in a glass melting furnace is well defined as the unsteady periodic or weak chaotic flow.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 2003.08a
• /
• pp.208-215
• /
• 2003

Convective-diffusion equations appear in various disciplines such as hydrology, chemical engineering and oceanography dealing with the transport problem of scalar quantities. Since it is nonlinear, numerical methods are generally used to obtain its solution. Very limited number of analytical solutions are available usually in cases when the convective velocity is constant or has a simple functional form (for some collection of the solutions, see Noye, 1987). There is however a continuing need to develop analytical solutions because of its practical importance. Analytical solutions of the convection-diffusion equation are valuable not only for the better understanding on the transport process but the verification of numerical schemes. (omitted)

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.2
• /
• pp.453-462
• /
• 1990

The interaction of forced convection-conduction with thermal radiation in laminar boundary layer over a plate fin is studied numerically. The analysis is based on complete solution whereby the heat conduction equation for the fin is solved simultaneously with the conservation equations for mass, momentum and energy in the fluid boundary layer adjacent to the fin. The fluid is a gray medium and diffusion(Rosseland) approximation is used to describe the radiative heat flux in the energy equation. The resulting boundary value problem are convection-conduction parameter N$_{c}$ and radiation-conduction parameter m, Prandtl number Pr. Numerical results are presented for gases with the Prandtl numbers of 0.7 & 5 with values of N$_{c}$ and M ranging from 0 to 10 respectively. The object of this study is to provide the first results on forced convection-radiation interaction in boundary layer flow over a semi-infinite flay plate which can be used for comparisons with future studies that will consider a more accurate expression for the radiative heat flux. The agreement of the results from the complete solution presented by E. M. Sparrow and those from this paper for the special case of M=0 is good. The overall rate of heat transfer from the fin considering radiative effect is higher than that from the fin neglecting radiative effect. The local heat transfer coefficient with radiative effect is higher than that without radiative effect. In the direction from tip to base, those coefficients decrease at first, attain minimum, and then increase. The larger values of N$_{c}$ M, Pr give rise to larger fin temperature variations and the fin temperature without radiative effect is always higher than that with radiative effect.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.31-48
• /
• 2010

In this article, we consider singularly perturbed Boundary Value Problems(BVPs) for second order Ordinary Differential Equations (ODEs) with Neumann boundary condition and discontinuous source term. A parameter-uniform error bound for the solution is established using the Streamline-Diffusion Finite Element Method (SDFEM) on a piecewise uniform meshes. We prove that the method is almost second order of convergence in the maximum norm, independently of the perturbation parameter. Further we derive superconvergence results for scaled derivatives of solution of the same problem. Numerical results are provided to substantiate the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.28 no.5_6
• /
• pp.1035-1054
• /
• 2010

We consider a mixed type singularly perturbed one dimensional elliptic problem with discontinuous source term. The domain under consideration is partitioned into two subdomains. A convection-diffusion and a reaction-diffusion type equations are posed on the first and second subdomains respectively. Two hybrid difference schemes on Shishkin mesh are constructed and we prove that the schemes are almost second order convergence in the maximum norm independent of the diffusion parameter. Error bounds for the numerical solution and its numerical derivative are established. Numerical results are presented which support the theoretical results.