• Journal of applied mathematics & informatics
• /
• 제29권3_4호
• /
• pp.813-829
• /
• 2011

In this article, we present a numerical scheme for solving singularly perturbed (i.e. highest -order derivative term multiplied by small parameter) Burgers-Huxley equation with appropriate initial and boundary conditions. Most of the traditional methods fail to capture the effect of layer behavior when small parameter tends to zero. The presence of perturbation parameter and nonlinearity in the problem leads to severe difficulties in the solution approximation. To overcome such difficulties the present numerical scheme is constructed. In construction of the numerical scheme, the first step is the dicretization of the time variable using forward difference formula with constant step length. Then, the resulting non linear singularly perturbed semidiscrete problem is linearized using quasi-linearization process. Finally, differential quadrature method is used for space discretization. The error estimate and convergence of the numerical scheme is discussed. A set of numerical experiment is carried out in support of the developed scheme.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2008년도 정기총회 및 학술발표대회
• /
• pp.431-434
• /
• 2008

Pratical dispersion-correction scheme is applicated to simulate the distant propagation of tsunami. This scheme is based on the leap-frog finite difference scheme for the linear shallow-water equations. The new scheme has the advantage of using the constant spatial grid size and time step size even in area of variable depths. And this new model constructed by using the 2nd upwind scheme, dynamic linking method, and staggered grid system. This model is simulated to near Sokcho harbor about The Central East Sea Tsunami in 1983. And this result is compared to tide gage and result of former model.

• 정보보호학회논문지
• /
• 제10권3호
• /
• pp.37-48
• /
• 2000

서울Visual cryptography scheme is a simple method in which can be directly decoded the secret information in human visual system without performing any cryptographic computations. This scheme is a kind of secret sharing scheme in which if a secret of image type is scattered to n random images(slides) and any threshold (or more) of them are stacked together the original image will become visible. In this paper we consider (2, n) visual cryptography scheme and propose a new construction method in which the number of expanded pixels can be reduced by using the sample matrix. The proposed scheme can futhermore distribute the multiple secret image to each group according to the difference of relative contrast.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.197-207
• /
• 2013

The Cahn-Hilliard equation was proposed as a phenomenological model for describing the process of phase separation of a binary alloy. The equation has been applied to many physical applications such as amorphological instability caused by elastic non-equilibrium, image inpainting, two- and three-phase fluid flow, phase separation, flow visualization and the formation of the quantum dots. To solve the Cahn-Hillard equation, many numerical methods have been proposed such as the explicit Euler's, the implicit Euler's, the Crank-Nicolson, the semi-implicit Euler's, the linearly stabilized splitting and the non-linearly stabilized splitting schemes. In this paper, we investigate each scheme in finite-difference schemes by comparing their performances, especially stability and efficiency. Except the explicit Euler's method, we use the fast solver which is called a multigrid method. Our numerical investigation shows that the linearly stabilized stabilized splitting scheme is not unconditionally gradient stable in time unlike the known result. And the Crank-Nicolson scheme is accurate but unstable in time, whereas the non-linearly stabilized splitting scheme has advantage over other schemes on the time step restriction.

• 한국자동차공학회논문집
• /
• 제6권5호
• /
• pp.86-96
• /
• 1998

A numerical study has been carried out of three-dimensional turbulent flows around a MIRA reference vehicle model both with and without wheels in computation. Two convective difference schemes with two k-$\varepsilon$ turbulence models are evaluated for the performance such as drag coefficient, velocity and pressure fields. Pressure coefficients along the surfaces of the model are compared with experimental data. The drag coefficient, the velocity and pressure fields are found to change considerably with the adopted finite difference schemes. Drag forces computed in the various regions of the model indicate that design change decisions should not rely just on the total drag and that local flow structures are important. The results also indicate that the RNG model with the QUICK scheme predicts fairly well the tendency of velocity and pressure fields and gives more reliable drag coefficient rather than the other cases.

• 대한수학회지
• /
• 제44권6호
• /
• pp.1397-1415
• /
• 2007

We consider an inverse heat conduction problem(IHCP) in a quarter plane which appears in some applied subjects. We want to determine the heat flux on the surface of a body from a measured temperature history at a fixed location inside the body. This is a severely ill-posed problem in the sense that arbitrarily "small" differences in the input temperature data may lead to arbitrarily "large" differences in the surface flux. A semi-discrete central difference scheme in time is employed to deal with the ill posed problem. We obtain some error estimates which also give the information about how to choose the step length in time. Some numerical examples illustrate the effects of the proposed method.

• 대한수학회보
• /
• 제54권6호
• /
• pp.1893-1912
• /
• 2017

In this paper, we aim to construct a nonstandard finite difference (NSFD) scheme to solve numerically a mathematical model for cholera epidemic dynamics. We first show that if the basic reproduction number is less than unity, the disease-free equilibrium (DFE) is locally asymptotically stable. Moreover, we mainly establish the global stability analysis of the DFE and endemic equilibrium by using suitable Lyapunov functionals regardless of the time step size. Finally, numerical simulations with different time step sizes and initial conditions are carried out and comparisons are made with other well-known methods to illustrate the main theoretical results.

• Journal of applied mathematics & informatics
• /
• 제27권5_6호
• /
• pp.1001-1015
• /
• 2009

In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 추계 학술대회논문집
• /
• pp.48-58
• /
• 1999

The Characteristic Flux Difference Splitting (CFDS) scheme designed to adapt the characteristic boundary conditions at the wall and inflow/outflow boundary planes satisfies Roe's property U, although the CFDS Jacobian matrix is decomposed by a product of elaborate transformation matrices and explicit eigenvalue matrix. When the CFDS algorithm, thus a variant of Roe's scheme, is applied straightforwardly to hypersonic flows over a blunt body, the strong bow shock gradually breaks down near the stagnation point. This numerical instability is widely observed by many researchers employing flux-difference method, known in the literature as the carbuncle phenomenon. Many remedies have been proposed and resulted in partial cures. When the idea of Sanders et al. which identifies the minimum eigenvalues near the discontinuity present is applied to CFDS method, it is shown that the instability problem can be controlled successfully. A few flux splitting methods have also been tested and results are compared against the Nakamori's Mach 8 blunt body flow.

• 대한기계학회논문집A
• /
• 제27권4호
• /
• pp.585-592
• /
• 2003

Among various sensitivity evaluation techniques, semi-analytical method(SAM) is quite popular since this method is more advantageous than analytical method(AM) and global finite difference method(FDM). However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified fur individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, an iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes and the error of SAM caused by numerical difference scheme is alleviated by using a Von Neumann series approximation considering the higher order terms for the sensitivity derivatives.