• Communications of the Korean Mathematical Society
• /
• v.24 no.3
• /
• pp.459-466
• /
• 2009

In this paper, a numerical scheme is introduced to solve the Cauchy problem for one-dimensional hyperbolic equations. The mesh points of the proposed scheme are distributed along characteristics so that the solution on the stencil can be easily and accurately computed. This is very important in reducing errors of the scheme because many numerical errors are generated when the solution is estimated over grid points. In addition, when characteristics intersect, the proposed scheme combines corresponding grid points into one and assigns new characteristic to the point in order to improve computational efficiency. Numerical experiments on the inviscid Burgers' equation have been presented.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1560-1565
• /
• 2003

Thermocapillary convection in an open cylindrical annulus heated from the inside wall is investigated by two-dimensional numerical simulations. The deformable free surface is obtained as a solution of the coupled transport equations at fixed Prandtl and aspect ratio. Only steady convection can be realized in this axisymmetric computations with either non-deformable or deformable surfaces. Dynamic free-surface deformations do not induce transitions to oscillatory convection even at large Reynolds numbers. Free surfaces are convex near the cold wall due to the stagnation point, and concave near the hot wall. Free surface deformation increases with increasing Ca at a fixed Re. Two peaks appear at the free surface with low Re, while additional ripples, four peaks, occur at larger Re. Thermocapillary convection in the open annulus interior is insensitive to variations in Ca.

• Kyungpook Mathematical Journal
• /
• v.59 no.1
• /
• pp.73-82
• /
• 2019

In this work, we give necessary and sufficient conditions under which every solution of a class of first-order neutral differential equations of the form $$(x(t)+p(t)x({\tau}(t)))^{\prime}+q(t)Hx({\sigma}(t)))=0$$ either oscillates or converges to zero as $t{\rightarrow}{\infty}$ for various ranges of the neutral coefficient p. Our main tools are the Knaster-Tarski fixed point theorem and the Banach's contraction mapping principle.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.