• Communications of the Korean Mathematical Society
• /
• v.34 no.1
• /
• pp.231-238
• /
• 2019

We obtain ordinary differential equations related with some nonlocal Schr$Schr{\ddot{o}}dinger$dinger equations. As applications, we prove finite time blow-up or extinction of solutions.

• Kyungpook Mathematical Journal
• /
• v.48 no.4
• /
• pp.665-684
• /
• 2008

Runge-Kutta-Nystr$\"{o}$m (RKN) methods provide a popular way to solve the initial value problem (IVP) for a system of ordinary differential equations (ODEs). Users of software are typically asked to specify a tolerance ${\delta}$, that indicates in somewhat vague sense, the level of accuracy required. It is clearly important to understand the precise effect of changing ${\delta}$, and to derive the strongest possible results about the behaviour of the global error that will not have regular behaviour unless an appropriate stepsize selection formula and standard error control policy are used. Faced with this situation sufficient conditions on an algorithm that guarantee such behaviour for the global error to be asympotatically linear in ${\delta}$ as ${\delta}{\rightarrow}0$, that were first derived by Stetter. Here we extend the analysis to cover a certain class of ODEs with low-order derivative discontinuities, and the class of ODEs with constant delays. We show that standard error control techniques will be successful if discontinuities are handled correctly and delay terms are calculated with sufficient accurate interpolants. It is perhaps surprising that several delay ODE algorithms that have been proposed do not use sufficiently accurate interpolants to guarantee asymptotic proportionality. Our theoretical results are illustrated numerically.

• Steel and Composite Structures
• /
• v.52 no.5
• /
• pp.525-541
• /
• 2024

This study uses the state-space approach to study the bending behavior of Levy-type functionally graded (FG) plates sandwiched between two piezoelectric layers. The coupled governing equations are obtained using Hamilton's principle and Maxwell's equation based on the efficient four-variable refined plate theory. The partial differential equations (PDEs) are converted using Levy's solution technique to ordinary differential equations (ODEs). In the context of the state-space method, the higher-order ODEs are simplified to a system of first-order equations and then solved. The results are compared with those reported in available references and those obtained from Abaqus FE simulations, and good agreements between results confirm the accuracy and efficiency of the approach. Also, the effect of different parameters such as power-law index, aspect ratio, type of boundary conditions, thickness-to-side ratio, and piezoelectric thickness are studied.

• Communications of the Korean Mathematical Society
• /
• v.36 no.3
• /
• pp.549-555
• /
• 2021

We show the solvability of the ordinary differential equations describing curves on sphere or hyperbolic space. Making use of the geometric properties of the equations, we derive explicit solution formulae.

• Journal of applied mathematics & informatics
• /
• v.35 no.1_2
• /
• pp.191-204
• /
• 2017

A new fifth-order weighted Runge-Kutta algorithm based on heronian mean for solving initial value problem in ordinary differential equations is considered in this paper. Comparisons in terms of numerical accuracy and size of the stability region between new proposed Runge-Kutta(5,5) algorithm, Runge-Kutta (5,5) based on Harmonic Mean, Runge-Kutta(5,5) based on Contra Harmonic Mean and Runge-Kutta(5,5) based on Geometric Mean are carried out as well. The problems, methods and comparison criteria are specified very carefully. Numerical experiments show that the new algorithm performs better than other three methods in solving variety of initial value problems. The error analysis is discussed and stability polynomials and regions have also been presented.

• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.861-869
• /
• 2023

This article proposes a novel approach to measure the performance of Safety-Critical Systems (SCS). Such systems contain multiple processing nodes that communicate with each other is modeled by a Petri nets (PN). The paper uses the PN for the performance evaluation of SCS. A set of ordinary differential equations (ODEs) is derived from the Petri net model that represent the state of the system, and the solutions can be used to measure the system's performance. The proposed method can avoid the state space explosion problem and also introduces new metrics of performance, along with their measurement: deadlock, liveness, stability, boundedness, and steady state. The proposed technique is applied to Shutdown System (SDS) of Nuclear Power Plant (NPP). We obtained 99.887% accuracy of performance measurement, which proves the effectiveness of our approach.

• Structural Engineering and Mechanics
• /
• v.4 no.6
• /
• pp.589-600
• /
• 1996

An explicit 4th order time integration scheme for solving the convection-diffusion equation is discussed in this paper. A system of ordinary differential equations are derived first by discretizing the spatial derivatives of the relevant PDE using the finite difference method. The integration of the ODEs is then carried out using a 4th order scheme and a self-adaptive technique based on the spatial grid spacing. For a non-uniform spatial grid, different time step sizes are used for the integration of the ODEs defined at different spatial points, which improves the computational efficiency significantly. A numerical example is also discussed in the paper to demonstrate the implementation and effectiveness of the method.

• 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.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Journal of Biomedical Engineering Research
• /
• v.41 no.3
• /
• pp.128-137
• /
• 2020

Cardiac electrophysiology studies often use simulation to predict how cardiac will behave under various conditions. To observe the cardiac tissue movement, it needs to use the high--resolution heart mesh with a sophisticated and large number of nodes. The higher resolution mesh is, the more computation time is needed. To improve computation speed and performance, parallel processing using multi-core processes and network computing resources is performed. In this study, we compared the computational speeds of CPU parallelization and GPU parallelization in virtual heart simulation for efficiently calculating a series of ordinary differential equations (ODE) and partial differential equations (PDE) and determined the optimal CPU and GPU parallelization architecture. We used 2D tissue model and 3D ventricular model to compared the computation performance. Then, we measured the time required to the calculation of ODEs and PDEs, respectively. In conclusion, for the most efficient computation, using GPU parallelization rather than CPU parallelization can improve performance by 4.3 times and 2.3 times in calculations of ODEs and PDE, respectively. In CPU parallelization, it is best to use the number of processors just before the communication cost between each processor is incurred.