• 대한수학회논문집
• /
• 제24권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.

• Journal of applied mathematics & informatics
• /
• 제27권1_2호
• /
• pp.441-452
• /
• 2009

In this paper, a fifth order numerical method is presented for solving singularly perturbed differential-difference equations with negative shift. In recent papers the term negative shift has been using for delay. Similar boundary value problems are associated with expected first exit time problem of the membrane, potential in models for neuron and in variational problems in control theory. In the numerical treatment for such type of boundary value problems, first we use Taylor approximation to tackle terms containing small shifts which converts it to a boundary value problem for singularly perturbed differential equation. The two point boundary value problem is transformed into general first order ordinary differential equation system. A discrete approximation of a fifth order compact difference scheme is presented for the first order system and is solved using the boundary conditions. Several numerical examples are solved and compared with exact solution. It is observed that present method approximates the exact solution very well.

• 충청수학회지
• /
• 제26권3호
• /
• pp.501-516
• /
• 2013

The well-known Vlasov-Poisson equation describes plasma physics as nonlinear first-order partial differential equations. Because of the nonlinear condition from the self consistency of the Vlasov-Poisson equation, many problems occur: the existence, the numerical solution, the convergence of the numerical solution, and so on. To solve the problems, a viscosity term (a second-order partial differential equation) is added. In a viscosity term, the Vlasov-Poisson equation changes into a parabolic equation like the Fokker-Planck equation. Therefore, the Schauder fixed point theorem and the classical results on parabolic equations can be used for analyzing the Vlasov-Poisson equation. The sequence and the convergence results are obtained from linearizing the Vlasove-Poisson equation by using a fixed point theorem and Gronwall's inequality. In numerical experiments, an implicit first-order scheme is used. The numerical results are tested using the changed viscosity terms.

• Journal of Advanced Marine Engineering and Technology
• /
• 제28권4호
• /
• pp.593-602
• /
• 2004

This paper is a study on the model of simultaneous heat and mass transfer process in the absorption of refrigerant vapor into a lithium bromide solution of water-cooled vertical plate absorber. The model can predict temperature and concentration profiles as well as the effect of Reynolds number on them. Also. the variations of the absorption heat and mass fluxes. and the heat and mass transfer coefficients have been investigated. The numerical result shows that the interface temperature and concentration decrease as film Reynolds number does. The absorption heat and mass fluxes, and the heat and mass transfer coefficients get their maximum values adjacent to the inlet solution. Analyses on a constant wall temperature condition have been also carried out to exam the reliability of the present numerical method by comparing to previous investigations.

• Journal of applied mathematics & informatics
• /
• 제28권3_4호
• /
• pp.611-624
• /
• 2010

A Numerical scheme based on collocation method using quartic B-spline functions is designed for the numerical solution of one-dimensional modified equal width wave (MEW) wave equation. Using Von-Neumann approach the scheme is shown to be unconditionally stable. Performance of the method is validated through test problems including single wave, interaction of two waves and use of Maxwellian initial condition. Using error norms $L_2$ and $L_{\infty}$ and conservative properties of mass, momentum and energy, accuracy and efficiency of the suggested method is established through comparison with the existing numerical techniques.

• Steel and Composite Structures
• /
• 제19권3호
• /
• pp.679-695
• /
• 2015

A numerical solution using finite difference method to evaluate the thermal buckling of simply supported FGM plate with variable thickness is presented in this research. First, the governing differential equation of thermal stability under uniform temperature through the plate thickness is derived. Then, the governing equation has been solved using finite difference method. After validating the presented numerical method with the analytical solution, the finite difference formulation has been extended in order to include variable thickness. The accuracy of the finite difference method for variable thickness plate has been also compared with the literature where a good agreement has been found. Furthermore, a parametric study has been conducted to analyze the effect of material and geometric parameters on the thermal buckling resistance of the FGM plates. It was found that the thickness variation affects isotropic plates a bit more than FGM plates.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1997년도 봄 학술발표회 논문집
• /
• pp.3-10
• /
• 1997

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

• Structural Engineering and Mechanics
• /
• 제68권6호
• /
• pp.713-720
• /
• 2018

In this paper a new numerical method to determine the elastic curve of the simply supported beams of variable cross-section is demonstrated. In general case it needs to solve linear or small nonlinear second order differential equations with prescribed boundary conditions. For numerical solution the initial values of the slope and the deflection of the end cross-section of the beam is necessary. For obtaining the initial values a lively procedure is developed: it is a special application of the shooting method because boundary value problems can be transformed into initial value problems. As a result of these transformations the initial values of the differential equations are obtained with high accuracy. Procedure is applied for calculating of elastic curve of a simply supported beam of variable cross-section. Results of these numerical procedures, analytical solution of the linearized version and finite element method are compared. It is proved that the suggested procedure yields technically accurate results.

• 대한수학회지
• /
• 제58권4호
• /
• pp.835-847
• /
• 2021

We study behavior of numerical solutions for a nonlinear eigenvalue problem on ℝⁿ that is reduced from a dispersion managed nonlinear Schrödinger equation. The solution operator of the free Schrödinger equation in the eigenvalue problem is implemented via the finite difference scheme, and the primary nonlinear eigenvalue problem is numerically solved via Picard iteration. Through numerical simulations, the results known only theoretically, for example the number of eigenpairs for one dimensional problem, are verified. Furthermore several new characteristics of the eigenpairs, including the existence of eigenpairs inherent in zero average dispersion two dimensional problem, are observed and analyzed.