• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제12권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.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2009년도 학술발표회 초록집
• /
• 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
• /
• 제26권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
• /
• 제42권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.

• 한국지구과학회지
• /
• 제34권1호
• /
• pp.51-58
• /
• 2013

대기 모델링 연구에서 시간 간격을 적절하게 결정하는 것은 중요한 문제이다. 본 연구에서는 비선형 대기 모형에서 수치 해의 시간 간격에 대한 민감도를 조사하였다. 이를 위해 간단한 무차원화된 역학 모형을 사용하여 시간 간격과 비선형성 인자를 바꾸어가며 수치 실험을 수행하였다. 실험 결과, 비선형성 인자가 영향을 줄 만큼 크지 않고 절단 오차를 무시할 수 있는 경우에는 수치 해가 시간 간격에 민감하지 않았다. 그러나 비선형성 인자가 큰 경우에는 수치 해가 시간 간격에 민감한 것으로 밝혀졌다. 이 경우, 시간 간격이 감소할수록 공간 필터의 강도가 증가하여 작은 규모의 현상이 약하게 모의되었다. 이는 일반적으로 시간 간격이 감소하면 절단 오차가 감소하여 더 정확한 수치 해가 도출된다는 사실과 상충한다. 이러한 충돌은 비선형 모형의 수치 해를 안정하게 하기 위해 공간 필터가 반드시 필요하기 때문에 피할 수 없다.

• 대한수학회보
• /
• 제47권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.

• 대한기계학회논문집B
• /
• 제21권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
• /
• 제57권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
• /
• 제25권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
• /
• 제26권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.