• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Korean Journal of Mathematics
• /
• v.29 no.2
• /
• pp.227-237
• /
• 2021

Opial inequality and its consequences are useful in establishing existence and uniqueness of solutions of initial and boundary value problems for differential and difference equations. In this paper we analyze Opial-type inequalities for convex functions. We have studied different versions of these inequalities for a generalized integral operator. Further difference of Opial-type inequalities are utilized to obtain generalized mean value theorems, which further produce various interesting derivations for fractional and conformable integral operators.

• Environmental Engineering Research
• /
• v.11 no.3
• /
• pp.164-171
• /
• 2006

Numerical modeling of the flow velocity fields for the near corona wire electrohydrodynamic (EHD) flow was conducted. The steady, two-dimensional momentum equations have been computed for a wire-plate type electrostatic precipitator (ESP). The equations were solved in the conservative finite-difference form on a fine uniform rectilinear grid of sufficient resolution to accurately capture the momentum boundary layers. The numerical procedure for the differential equations was used by SIMPLEST algorithm. The Phoenics (Version 3.5.1) CFD code, coupled with Poisson's electric field, ion transport equations and the momentum equation with electric body force were used for the numerical simulation and the Chen-Kim ${\kappa}-{\varepsilon}$ turbulent model numerical results that an EHD secondary flow was clearly visible in the downstream regions of the corona wire despite the low Reynolds number for the electrode ($Re_{cw}=12.4$). Secondary flow vortices caused by the EHD increases with increasing discharge current or EHD number, hence pressure drop of ESP increases.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Journal of the Korean Mathematical Society
• /
• v.50 no.4
• /
• pp.715-725
• /
• 2013

In this paper, a finite difference scheme for a two-point boundary value problem of 2nth-order ordinary differential equations is presented. The convergence and uniqueness of the solution for the scheme are proved by means of theories on matrix eigenvalues and norm. Numerical examples show that our method is very simple and effective, and that this method can be used effectively for other types of boundary value problems.

• Communications of the Korean Mathematical Society
• /
• v.37 no.3
• /
• pp.939-955
• /
• 2022

In this paper, we study a first-order non-linear singularly perturbed Volterra integro-differential equation (SPVIDE). We discretize the problem by a uniform difference scheme on a Bakhvalov-Shishkin mesh. The scheme is constructed by the method of integral identities with exponential basis functions and integral terms are handled with interpolating quadrature rules with remainder terms. An effective quasi-linearization technique is employed for the algorithm. We establish the error estimates and demonstrate that the scheme on Bakhvalov-Shishkin mesh is O(N^-1) uniformly convergent, where N is the mesh parameter. The numerical results on a couple of examples are also provided to confirm the theoretical analysis.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.2-7
• /
• 2008

A highly efficient moving least squares finite difference method (MLS FDM) for heat transfer analysis of composite material with interface. In the MLS FDM, governing differential equations are directly discretized at each node. No grid structure is required in the solution procedure. The discretization of governing equations are done by Taylor expansion based on moving least squares method. A wedge function is designed for the modeling of the derivative jump across the interface. Numerical examples showed that the numerical scheme shows very good computational efficiency together with high aocuracy so that the scheme for heat transfer problem with different heat conductivities was successfully verified.

• Journal of KIISE:Computer Systems and Theory
• /
• v.32 no.9
• /
• pp.445-456
• /
• 2005

We propose a fast and accurate fluid solver of the Wavier-Stokes equations for the physics-based fluid simulations. Our method utilizes the solution of the Stokes equation as an initial guess for the velocity of the nonlinear term in the Wavier-Stokes equations. By guessing the initial velocity close to the exact solution of the given nonlinear differential equations, we can develop remarkably accurate and stable fluid solver. Our solver is based on the implicit scheme of finite difference methods, that makes it work well for large time steps. Since we employ the ADI method, our solver is also fast and has a uniform computation time. The experimental results show that our solver is excellent for fluids with high Reynolds numbers such as smoke and clouds.

• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.