Order Ordinary Differential Equations Research Articles

Sur une équation définie par une réduction de la hiérarchie de Drinfeld et Sokolov

A seventh order ordinary differential equation (ODE) arising by reduction of the Drinfeld-Sokolov hierarchyis shown to be identical to a similarity reduction of an equationin the hierarchy of Sawada-Kotera.We also exhibit its link with a particular F-VI,a fourth order ODE isolated by Cosgrove which is likely to define a higher order Painlev\'e function.

New exact free vibration solution of two-dimensional decagonal quasicrystal cylindrical shells in the symplectic framework

An accurate free vibration analysis of two-dimensional decagonal quasicrystal (QC) cylindrical shells is performed under the framework of symplectic mechanics. Unlike the classical Lagrangain system, the governing equations are reduced into a set of lower order ordinary differential equations, and therefore exact solutions of displacements, rotation angles, internal forces and bending moments in the phonon and phason fields are simultaneously obtained and expressed in terms of symplectic vectors. The present results are compared with existing literature and finite element solutions, and excellent agreements are observed. Furthermore, the effects of key influencing factors on the free vibration characteristics are revealed also.

Development of a fully implicit ODE-solver for containment analysis code

The thermal–hydraulic dynamics in containment are governed by a system of stiff ordinary differential equations (ODEs). A fully implicit discretization scheme is adopted to discretize these ODEs in order to mitigate the effects of stiffness. In comparison with explicit or semi-implicit discretization schemes that are subject to Courant limits on time steps, the fully implicit discretization scheme is more suitable for a containment analysis code that focuses on predicting both short-term and long-term thermal–hydraulic parameters after an accident. This study introduces a general-purpose ODE solver for the containment analysis code. The outline of the solver is as follows: The fully implicit discrete equations lead to a large set of nonlinear equations that need to be solved using Newton’s iterative method. The partial derivative components in the Jacobi matrix are calculated by the perturbation method using finite difference approximation, which avoids the complicated derivation of partial derivatives. The scaling modification technique is incorporated into this ODE solver to deal with significant differences in unknown variable magnitudes, and the line search method is introduced to address the difficulty of obtaining an accurate root estimate with Newton’s method when the initial guess is far from the actual root. This proposed ODE solver was applied to two typical stiff ODE problems to test its stiffness-suppressed ability and to demonstrate that this proposed solver can perform calculations with a very large time step. Then, the CASSIA code, a containment analysis code developed by China Nuclear Power Technology Research Institute Co., Ltd (CNPRI), equipped with this ODE solver, was applied to the CSNI (Committee on the Safety of Nuclear Installations) benchmark problem and the Carolinas Virginia Tube Reactor (CVTR) test 3 problem to preliminarily demonstrate that the proposed ODE solver can perform containment thermal–hydraulic analysis correctly. This study could provide references for the development of a home-made containment analysis code.

Review of Adams-Bashforth method for numerical solution of first order ordinary differential equations

This project work presents a comprehensive review of the Adams-Bashforth method for the numerical solution of explicit first-order ordinary differential equations (ODEs). The study begins with a historical overview of the development of ordinary differential equations (ODEs), tracing back to the seminal works of Newton, Leibniz, and their contemporaries. The evolution of differential equations as a distinct mathematical discipline and their wide-ranging applications across various fields are discussed. Furthermore, the paper provides an in-depth analysis of the Adams-Bashforth method, a prominent numerical technique for solving ODEs. The method is derived and analyzed using Mathematica Software, demonstrating its efficacy and accuracy in approximating solutions to ODEs. A comparison with existing methods such as the Euler method and the Runge-Kutta method highlights the advantages of the Adams-Bashforth method, particularly in terms of computational efficiency and accuracy. Moreover, fundamental concepts related to ODEs, including the initial value problem, existence, and uniqueness of solutions, are explored in the context of numerical solution methods. The properties of linear multistep methods and their relevance to solving ODEs efficiently are also discussed. In conclusion, this project work provides valuable insights into the Adams-Bashforth method as a powerful tool for the numerical solution of first-order ODEs. By leveraging modern computational tools such as Mathematica Software, the method offers a practical and efficient approach to solving complex differential equations encountered in various scientific and engineering applications.

On Period Annuli and Induced Chaos

Nontrivial period annuli in the second order ordinary differential equation are continua of periodic trajectories that contain inside more than one critical point. They can appear in conservative equations, which are known to have no attractors. Nevertheless, according to some authors, their behavior may be done chaotic by adding a periodic external force. Is the period of the external force correlated with periods of solutions in period annuli? Is the chaotic behavior of a solution dependent on the initial value and, in turn, on a certain periodic annulus? These, and related questions are studied in the article.

Investigating the new perspectives of Caudrey–Dodd–Gibbon equation arising in quantum field theory

The main purpose of the paper is obtaining the analytical results for beta fractional Caudrey–Dodd–Gibbon equation which is used to resolve complex problems in fluid dynamics, chemical kinetics, plasma physics, quantum field theory, crystal dislocations, and nonlinear optics by using auxiliary method. Beta derivative is a useful fractional operator due to satisfying basic properties of integer order derivative and also, allows us using chain rule and wave transform to turn nonlinear fractional partial differential equations into integer order ordinary differential equations. By the way many analytical methods can be applied to these equations. In order to understand the physical features of the solutions, 3D and 2D graphical illustrations are given. Finally, authors expect that the obtained solutions may give a deep insight for the explanation of physical phenomena in the fluid dynamics, chemical kinetics, plasma physics, quantum field theory, crystal dislocations, and nonlinear optics.

A Chebyshev polynomial based block integrator for the direct numerical solution of fourth order ordinary differential equations

This paper introduces an innovative method for numerically integrating fourth-order initial value problems by utilizing Chebyshev polynomials as the fundamental basis function. The block integrator base on Chebyshev polynomial demonstrates significant improvements in accuracy and stability, rendering it a valuable tool across various scientific and engineering fields. By leveraging the characteristics of Chebyshev polynomials, this approach accurately estimates solutions for fourth-order differential equations without reducing it to a system of first order Ordinary Differential Equations while at the same time effectively managing error accumulation within a block integration framework and thereby enhancing its accuracy over extended intervals. Through rigorous numerical experiments, the effectiveness and reliability of the new integrator are demonstrated and compared with existing methods. The new method is consistent, zero stable and convergent. The method also shows an appreciable error constants. The new method performed better in terms of accuracy than the existing methods in the literature in both linear and nonlinear problems.

Explicit Runge-Kutta Method for Evaluating Ordinary Differential Equations of type v^vi = f(u, v, v′)

The initiative of this paper is to present the Runge Kutta Type technique for the development of mathematical solutions to the problems concerning to ordinary differential equation of order six of structure v vi = f(u, v, v′ ) denoted as RKSD with initial conditions. The three and four stage Runge-Kutta methods with order conditions up to order seven (RKSD7) have been designed to evaluate global and local truncated errors for the ordinary differential equation of order six. The framework and evaluation of equations with their results are well established to obtain the effectiveness of RK method towards implicit function satisfying the required initial conditions and for obtaining zero-stability of RKSD7 in terms of their accuracy with maximum precision under minimal processing.

Analysis and application of MLPG7 for diffusion equations with nonlinear reaction terms

This paper extends the recently proposed variant of meshless local Petrov Galerkin (MLPG) method, i.e., MLPG7, for solving time dependent PDEs. As test function, the method uses a novel modification of fundamental solution of Laplace operator that not only the test function itself but also its derivative vanish on boundary of local subdomains. Therefore, more stable local integral equations are obtained by replacing a boundary integral of derivatives of field variables with a domain integral of field variables itself. MLPG7 formulation converts the nonlinear governing equation into a system of first order ordinary differential equations (ODEs). For this system of ODEs, a theoretical stability analysis is carried out for the case of regularly spaced nodal points. Stability and convergence conditions for fully discretized equations by Euler (forward and backward) and Crank–Nicolson methods are investigated. The nonlinear term is treated iteratively. A numerical study investigates the convergence and stability of the method. A comparison is also done with some well-known methods and the results reveal the excellence of MLPG7. The convergence of the iterative approach is numerically verified by another illustrative test example. As nonlinear test problems, Allen–Cahn equations are studied.

Model order reduction, a novel method using krylov sub-spaces and genetic algorithm

Model Order Reduction (MOR) of complex and large systems in Electrical engineering, continuous to be an attractive field for Engineers and Scientists over the last few decades, this complexity of models makes the control designs and simulation using Computer Aided Design (CAD) more and more difficult and consuming a lot of time. There for, accurate, robust and fast algorithms for simulation are needed. The goal of MOR is to replace the original system by an appropriate reduced system which preserves the main properties of the original one such that stability and passivity. Several analytical MOR techniques have been proposed in the literature over the past few decades, to approximate high order linear dynamic systems like Krylov sub-space techniques and SVD (Singular Value Decomposition) techniques. However, most of these techniques lead to computationally demanding, time consuming, iterative procedures that usually result in non-robustly stable models with poor frequency response resemblance to the original high order model in some frequency ranges. Recently a set of new techniques based on Artificial Intelligence (AI) were proposed in [1] for MOR. This article considers the problem of model order reduction of Linear Time In varying (LTI) systems. It is described by first and second order ordinary differential equations model. A tow steps method for model order reduction of LTI systems is proposed here, which combined features of an analytic technique (Krylov approach) and an AI technique (Genetic Algorithm). In the first step, the size of the original model is reduced to an intermediate order, using an analytical technique based on Krylov sub-spaces. In the final step of the reduction process, an AI approach based on Genetic Algorithm (GA) is applied to obtain an optimized nominal model.

Fourier series method for the stability solution of simply supported thin beams on two-parameter elastic foundations of the Pasternak, Filonenko-Borodich, Hetenyi or Vlasov models

The analysis of stability problems of beams on two-parameter foundations (Bo2PFs) is an important part of their design for compressive loads. This work presents novel first principles derivation of the governing differential equations of elastic stability (GDES) of thin beams resting on two-parameter elastic foundations of the Pasternak, Filonenko-Borodich, Hetenyi or Vlasov models. The requirements of translational and rotational equilibrium of all the applied, reactive and internal forces on an infinitesimal segment of the Bo2PF and the laws of infinitesimal calculus were used to formulate the GDES as a fourth order ordinary differential equation (ODE) in terms of the transverse displacement function ux. The GDES is non-homogeneous in the presence of applied transverse load qx but homogeneous when qx vanishes. This study presents the Fourier series method (FSM) for solving the governing differential equation of stability (GDES) for the case of Dirichlet boundary conditions. The FSM has the advantage of amenability to differentiation, and integration due to the orthogonality properties of the sinusoidal functions. Implementation of the FSM by assuming the unknown function in the GDES as a Fourier series of infinite terms and the exploitation of orthogonalization simplifies the problem to an algebraic eigenvalue problem which is the characteristic buckling equation. The exact eigenvalues are found by algebraic solution of the buckling equation. The exact eigenvalues were used to find the exact buckling loads and the exact buckling load coefficients. The critical buckling load was found to correspond to the first buckling mode (n= 1), and is identical with previous solutions in the literature. Numerical calculations for the critical buckling load parameters Kcr were presented for the Bo2PF problem for values of the dimensionless foundation parameters k-1= 0, k-2= 0; k-1= 100, k-2= 0; k-1= 0, k-2= 1; k-1= 100, k-2= 100; k-1= 0, k-2= 2.5; k-1= 100, k-2= 2.5. The present solutions were compared with previous solutions for Kcr in the literature. The comparison shows that the present FSM results are identical with previous results obtained using various other methods such as Recursive Differentiation Method, Finite Element Method, Generalized Integral Transform Method (GITM) and Stodola-Vianello Iteration Method. The study has illustrated the effectiveness of the FSM for solving Bo2PFs.

Penyelesaian khusus persamaan diferensial biasa ordo dua linier tak homogen dengan koefisien konstan untuk fungsi bagian demi bagian

Spring mechanical vibration motion system with a damped degree of freedom and influenced by external forces, is mathematically expressed as an ordinary differential equation of order of two linear constant coefficients that are not homogeneous. If an external force acts on a stationary system expressed as a continuous function f(t) for any time t, then the system will experience mechanical vibrational motion which mathematically the equation of motion can be expressed as a superposition. The equation consists of as a solution to a homogeneous form with mechanical vibrations as a solution to a particular form. In terms of the particular solution this article will show a mathematical way when f(t) is a continuous function section by part which is defined at an interval, such that the mechanical vibration motion equation is at the same time a special solution of the equation the mechanical vibration system.
 Keywords: Vibration, Impulse Functions, ConvolutionMSC2020: 34A37

A systematic approach to obtain the analytical solution for linear second order ordinary differential equations: part I

A systematic approach to obtain the analytical solution for linear second order ordinary differential equations: part I

C PROGRAM AS A TOOL FOR THE TEACHING OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATION

Second order ordinary differential equation (ODE) has many applications in science and engineering. Undergraduate students in science and engineering departments must study this subject in order to understand other subjects which are related to real applications they will encounter later. There are many excellent textbooks on differential equations where the students can study the theory and solve the problems. However, a textbook cannot give a quick answer for a problem particularly when the problem is quite difficult. A good choice is to use software such as Maple or Mathematica. However, this software is not always available for the teaching and purchasing it is usually beyond the ability of a student or even a lecturer. On the other hand, lecturers who want to create problems by themselves will follow the theory of the ODE. While creating the problems may not be difficult, answering them are harder. Problems which are very easy to answer are less worthless because they will not increase the students’ knowledge. Here comes the solution. A C program has been created to help lecturers create problems and solve them quickly. The program is interactive and can be easily understood by anyone who has basic theory of ODE. No knowledge of programming is needed; a user just runs it and follows the instruction. Students can also use the program to sharpen their knowledge. They can compare the solution of a problem they have solved with the answer given by the program. While commercial software such as Maple and Mathematica is very powerful, they cannot be used without writing necessary commands to solve a problem.

Boundedness of Vector Linéard Equation with Multiple Variable Delays

In this article, we consider a system of ordinary differential equations (ODEs) of second order with two variable time delays. We obtain new conditions for uniform ultimate bounded (UUB) solutions of the considered system. The technique of the proof is based on the Lyapunov–Krasovskii functional (LKF) method using a new LKF. The main result of this article extends and improves a recent result for ODEs of second order with a constant delay to a more general system of ODEs of second order with two variable time delays. In this particular case, we also give a numerical example to verify the application of the main result of this article.

Estimation of the Impact of Vaccination Intervention on Recovered Coronavirus Patients

This work estimated the impact of vaccination intervention on coronavirus patients who have recovered from the disease and the vulnerability index of the recovered population due to the impact of vaccine was also investigated. This work adopted a numerical solution to study the continuous dynamical system of linear first order differential equations describing a SEIR (Susceptible, Exposed, Infected, Recovered) model on the spread of Coronavirus Disease – 2019 (COVID-19). To tackle this problem, MATLAB ordinary differential equation of order 45 (ODE45) numerical method was adopted for the analysis. The vulnerability index of the recovered population was low due to the impact of vaccine meaning that the recovered population will gain immunity and they will not be re-infected. The study recommended that coronavirus patients who have recovered from the disease should ensure that they have vaccination administered to them to avoid re-occurrence of the virus attack as an intervention strategy.

Initial Value Solvers for Direct Solution of Fourth Order Ordinary Differential Equations in a Block from Using Chebyshev Polynomial as Basis Function

The numerical computation of fourth order ordinary differential equations cannot be gloss over easily due to its significant and importance. There have been glowing needs to find an appropriate numerical method that will handle effectively fourth order ordinary differential equations without resolving such an equation to a system of first order ordinary differential equations. To this end, this presentation focuses on direct numerical computation to fourth order ordinary differential equations without resolving such equations to a system of first order ordinary differential equations. The method is not predictor – corrector one due to its limitation in the level of accuracy. The method is order wise christened “Block Method” which is a self-starting method. In order to achieve this objective, Chebyshev polynomial is hereby used as basis function.

Modeling the effect of shape deformation induced by gravity on the evaporation of pendant and sessile drops

Pendant and sessile drops form a spherical cap only in the absence of gravity. The effect of gravity on drop shape is often neglected on the basis of the assumption that the drop size is smaller than the capillary length [Lc=(σ/gρ)1/2], although the deformation may not be fully negligible even in those cases. This paper focuses on evaluation of the effect that deformation due to gravity has on the evaporation characteristics of pendant and sessile drops. The drop shape is described by the Bashforth–Adams equation, a non-linear second order ordinary differential equation, which is solved numerically using a Runge–Kutta method with variable time steps. Under quasi-steady approximation, the species and energy conservation equations in the gas phase have analytical solutions, even for temperature-dependent gas thermophysical properties, once the solution of a basic Laplace problem is known. The Laplace equation is solved in axial symmetric geometry by using COMSOL Multiphysics®, for a wide range of drop sizes and contact angles, yielding vapor distribution, vapor fluxes, and evaporation rates. Comparison with the results from drops of same size in microgravity (i.e., having a spherical cap shape) shows that the effect is also perceptible for drops with a size smaller than the capillary length and that it can become quite important for those with larger sizes. Complementary results are found for sessile and pendant drops with respect to wall wettability, suggesting that the phenomenon can be analyzed using a unitary approach.

On polynomial solutions of certain finite order ordinary differential equations

Given a finite order differential operator, some properties and relations satisfied by its polynomial eigenfunctions are studied. Under certain restrictions, such eigenfunctions are explicitly obtained, as well as the corresponding eigenvalues. Also, some linear transformations are applied to sequences of eigenfunctions and a necessary condition for this to be a sequence of eigenfunctions of a new differential operator is obtained. These results are applied to the particular case of classical Hermite polynomials.

On the Dynamic Characteristics of Orthotropic Rectangular Plates under the Influence of Moving Distributed Masses and Resting on a Variable Elastic Pasternak Foundation

This work studies the dynamic characteristics of orthotropic rectangular plates under the influence of moving distributed masses and resting on a variable elastic Pasternak foundation. The governing equation is a fourth order partial differential equation with variable and singular coefficients. The solutions to the problem are obtained by transforming the fourth order partial differential equation for the problem to a set of coupled second order ordinary differential equations using the technique of Shadnam et al, then simplified using asymptotic method of Struble [11]. The closed form solution is analyzed, resonance conditions are obtained and the results are depicted graphically for both cases of moving distributed mass and moving distributed force