• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.4
• /
• pp.201-202
• /
• 2008

Because of the importance of suction or injection in the fields of aerodynamics, space science and many other industrial applications, our present study is motivated. The effect of natural convection on MHD flow past a vertical plate with suction or injection is studied. We have tried to solve the dimensionless governing equations by using finite difference scheme. To ensure the validity of our numerical solutions, we have compared our numerical solutions for temperature and velocity for the case of suction and injection for unit Prandtl number with the available exact solutions in the literature. The corresponding codes were written in Mathematica 5.0 for calculating numerical solutions for temperature and velocity and the comparison between the exact and numerical solutions. For the purpose of discussing the results some numerical calculations are carried out for non-dimensional temperature T, velocity u, skin friction ${\tau}$ and the Nusselt number $N_u$, by making use of it, the rate of heat transfer is studied.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.905-908
• /
• 2009

In the management of groundwater in coastal areas, saltwater intrusion associated with extensive groundwater pumping, is an important problem. The groundwater optimization model is an advanced method to study the aquifer and decide the optimal pumping rates or optimal well locations. Cheng and Park gave the analytical solutions to the optimization problems basing on Strack's analytical solution. However, the analytical solutions have some limitations of the property of aquifer, boundary conditions, and so on. A simulation-optimization numerical method presented in this study can deal with non-homogenous aquifers and various complex boundary conditions. This simulation-optimization model includes the sharp interface solution which solves the same governing equation with Strack's analytical solution, therefore, the freshwater head and saltwater thickness should be in the same conditions, that can lead to the comparable results in optimal pumping rates and optimal well locations for both of the solutions. It is noticed that the analytical solutions can only be applied on the infinite domain aquifer, while it is impossible to get a numerical model with infinite domain. To compare the numerical model with the analytical solutions, calculation of the equivalent boundary flux was planted into the numerical model so that the numerical model can have the same conditions in steady state with analytical solutions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.1
• /
• pp.49-66
• /
• 2022

The G-Euler process has been proposed to overcome the difficulties of the calculation of the exponential function of the Jacobian. It is an explicit method that uses the exponential function of the scalar skew-symmetric matrix. We define the moving shapes of true solutions and the moving shapes of numerical solutions. It is discussed whether the moving shape of the numerical solution matches the moving shape of the true solution. The match rates of these two kinds of moving shapes are sequentially calculated by the G-Euler process without using the true solution. It is shown that the closer the minimum match rate is to 100%, the more closely the numerical solutions follow the true solutions to the end. The minimum match rate indicates the reliability of the numerical solution calculated by the G-Euler process. The graphs of the Lorenz system in Perko [1] are different from those drawn by the G-Euler process. By the way, there is no basis for claiming that the Perko's graphs are reliable.

• Journal of applied mathematics & informatics
• /
• v.42 no.2
• /
• pp.321-331
• /
• 2024

We consider the numerical treatment of blow-up problems having non-monotonic singular solutions that tend to infinity at some point in the domain. The use of standard numerical methods for solving problems with blow-up solutions can lead to significant errors. The reason being that solutions of such problems have singularities whose positions are unknown in advance. To be able to integrate such non-monotonic blow-up problems, we describe and use a method of non-local transformations. To show the efficiency of the method, we present a comparison of exact and numerical solutions in addition to some comparison with the work of other authors.

• Journal of the Korean earth science society
• /
• v.34 no.1
• /
• pp.51-58
• /
• 2013

An appropriate determination of time step is one of the important problems in atmospheric modeling. In this study, we investigate the sensitivity of numerical solutions to time step in a nonlinear atmospheric model. For this purpose, a simple nondimensional dynamical model is employed, and numerical experiments are performed with various time steps and nonlinearity factors. Results show that numerical solutions are not sensitive to time step when the nonlinearity factor is not influentially large and truncation error is negligible. On the other hand, when the nonlinearity factor is large (i.e., in a highly nonlinear regime), numerical solutions are found to be sensitive to time step. In this situation, smaller time step increases the intensity of the spatial filter, which makes small-scale phenomena weaken. This conflicts with the fact that smaller time step generally results in more accurate numerical solutions owing to reduced truncation error. This conflict is inevitable because the spatial filter is necessary to stabilize the numerical solutions of the nonlinear model.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.6
• /
• pp.1225-1234
• /
• 2010

Numerical solutions for the fractional differential dispersion equations with nonlinear forcing terms are considered. The backward Euler finite difference scheme is applied in order to obtain numerical solutions for the equation. Existence and stability of the approximate solutions are carried out by using the right shifted Grunwald formula for the fractional derivative term in the spatial direction. Error estimate of order $O({\Delta}x+{\Delta}t)$ is obtained in the discrete $L_2$ norm. The method is applied to a linear fractional dispersion equations in order to see the theoretical order of convergence. Numerical results for a nonlinear problem show that the numerical solution approach the solution of classical diffusion equation as fractional order approaches 2.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

• Kyungpook Mathematical Journal
• /
• v.57 no.4
• /
• pp.651-665
• /
• 2017

Some computational fluid dynamics problems concerning the thin films flow of viscous fluid with a free surface and draining or coating fluid-flow problems can be delineated by third-order ordinary differential equations. In this paper, the aim is to introduce the numerical solutions of the boundary value problems of such equations by variational iteration method. In this paper, it is shown that the third-order boundary value problems can be written as a system of integral equations, which can be solved by using the variational iteration method. These solutions are gleaned in terms of convergent series. Numerical examples are given to depict the method and their convergence.

• Geomechanics and Engineering
• /
• v.25 no.4
• /
• pp.267-273
• /
• 2021

In this work, we compare the analytical solutions with the numerical solutions for photothermal interactions in semiconductor medium containing cylindrical cavity. This paper is devoted to a study of the photothermal interactions in semiconductor medium in the context of the coupled photo-thermal model. The basic equations are formulated in the domain of Laplace transform and the eigenvalue scheme are used to get the analytical solutions. The numerical solution is obtained by using the implicit finite difference method (IFDM). A comparison between the analytical solution and the numerical solutions are obtained. It is found that the implicit finite difference method (IFDM) is applicable, simple and efficient for such problems.

• Journal of applied mathematics & informatics
• /
• v.26 no.3_4
• /
• pp.793-809
• /
• 2008

The aim of this article is focused on providing numerical solutions for system of second order robot arm problem using the RK-eight stage seventh order limiting formulas. The parameters governing the arm model of a robot control problem have also been discussed through RK-eight stage seventh order limiting algorithm. The precised solution of the system of equations representing the arm model of a robot has been compared with the corresponding approximate solutions at different time intervals. Results and comparison show the efficiency of the numerical integration algorithm based on the absolute error between the exact and approximate solutions. Based on the numerical results a thorough comparison is carried out between the numerical algorithms.