New Seven-Step Numerical Method for Direct Solution of Fourth Order Ordinary Differential Equations

This paper presented a linear multistep method for solving fourth order initial value problems of ordinary differential equations. Collocation and interpolation methods are adopted in the derivation of the new numerical scheme which is further applied to finding direct solution of fourth order ordinary differentiation equations. This implementation strategy is more accurate and efficient than Adams–Bashforth Method solution. The newly derive scheme have better stabilities properties than that of the Adams-Bashforth Method. Numerical examples are included to illustrate the reliability and accuracy of the new methods.Keywords: Linear Multistep Methods, Region of Absolute Stability, Zero-Stability, Error Analysis, Collocation.

Computational methods for reaction–diffusion problems for fourth order ordinary differential equations with a small parameter at the highest derivative

This paper carries a different approach of collection and interpolation to develop a tenth block method for the numerical solution of linear or nonlinear ordinary differential equations of fourth order with initial conditions. The method has been implemented at the selected mesh points to generate a direct tenth block method through Taylor series. Some critical properties of this method such as zero stability, order of the method, and convergence have been analysed. Two numerical tests have taken to make a comparison of the approximate results with exact as well as results of other authors.

From the biharmonic equation of the plane problem in the polar coordinate system and taking into account the variable-separable form of the partial solutions, a homogeneous ordinary differential equation (ODE) of the fourth order is deduced. Our study is based on the investigation of the behavior of the coefficients of the above fourth order ODE, which are functions of the radial coordinate r. According to the proposed investigation additional terms, φ¯−m(r,θ)(1≤m≤n) other than the usually tabulated in the Michell solution (1899, “On the Direct Determination of Stress in an Elastic Solid, With Application to the Theory of Plates,” Proc. Lond. Math. Soc., 31, pp. 100–124) are found. Finally the stress and the displacement fields due to each one additional term of φ¯−m(r,θ) are determined.

Criteria for fourth-order ordinary differential equations to be linearizable by contact transformations

In this paper we present fifth order Runge-Kutta method (RK5) for solving initial value problems of fourth order ordinary differential equations. In this study RK5 method is quite efficient and practically well suited for solving boundary value problems. All mathematical calculation performed by MATLAB software for better accuracy and result. The result obtained, from numerical examples, shows that this method more efficient and accurate. These methods are preferable to some existing methods because of their simplicity, accuracy and less computational cost involved.

In this paper, the Elzaki transform method is used for solving two-dimensional (2D) elasticity problems in plane polar coordinates. Airy stress function was used to express the stress compatibility equation as a biharmonic equation. Elzaki transform was applied with respect to the radial coordinate to a modified form of the stress compatibility equation, and the biharmonic equation simplified to a fourth order ordinary differential equation (ODE). The general solution for the Airy stress potential function in the Elzaki transform space was obtained by solving the ODE. By inversion, the general solution for the Airy stress potential function was obtained in the physical domain space variables in terms of four unknown integration constants. Normal stresses and shear stress fields were also determined for the general case of 2D elasticity problems. The Flamant problem was solved as a particular illustration of 2D elasticity problems. The stress boundary conditions and the requirement of equilibrium of the internal stress resultants and the external forces were used simultaneously to determine the four constants of integration. The Airy stress potential function and the normal and shear stress fields were thus completely determined. The principle of superposition was used to obtain the elasticity solutions for the stress fields in the elastic half plane due to strip load of infinite extent, and solutions for horizontal stresses on smooth rigid retaining walls due to strip loads and parallel line load of infinite extent acting on the elastic half plane.

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.

To a system of second-order ordinary differential equations one can assign a canonical nonlinear connection that describes the geometry of the system. In this paper, we develop a geometric setting that also allows us to assign a canonical nonlinear connection to a system of higher-order ordinary differential equations (HODE). For this nonlinear connection we develop its geometry, and explicitly compute all curvature components of the corresponding Jacobi endomorphism. Using these curvature components we derive a Jacobi equation that describes the behavior of nearby geodesics to a HODE. We motivate the applicability of this nonlinear connection using examples from the equivalence problem, the inverse problem of the calculus of variations, and biharmonicity. For example, using components of the Jacobi endomorphism we express two Wuenschmann-type invariants that appear in the study of scalar third- or fourth-order ordinary differential equations.

The block hybrid collocation method with three off-step points is proposed for the direct solution of fourth order ordinary differential equations. The interpolation and collocation techniques are applied on basic polynomial to generate the main and additional methods. These methods are implemented in block form to obtain the approximation at seven points simultaneously. Numerical experiments are conducted to illustrate the efficiency of the method. The method is also applied to solve the fourth order problem from ship dynamics.

A boundary value technique for boundary value problems for singularly perturbed fourth-order ordinary differential equations

A novel and efficient algorithm for the concurrent numerical solution of second, third and fourth orders of ordinary differential equations is studied in this article. Collocation and interpolation technique was employed in the derivation of the method and power series approximate solution was used as the interpolating polynomial. The fourth derivative of the power series was collocated at the entire grid and off-grid points, while the fifth and sixth derivatives of the polynomial were collocated at the endpoint only. Proper investigation of the basic properties of the method was done. The results showed that the new block method applied on a nonlinear second, third and fourth orders of ordinary differential equations were better in terms of accuracy than the existing methods. The proposed method takes away the burden of developing a separate method for the solution of second, third and fourth order initial value problem of ordinary differential equations.

In this paper, the purpose is to investigate the stability and variation of natural frequency of a cracked cantilever beams subjected to follower force and tip mass. In addition, an analysis of the flutter instability(flutter critical follower force) of a cracked cantilever beam as slenderness ratio and crack severity is investigated. The governing differential equations of a Timoshenko beam subjected to an end tangential follower force is derived via Hamilton's principle. The two coupled governing differential equations are reduced to one fourth order ordinary differential equation in terms of the flexural displacement. Finally, the influence of the slenderness ratio and crack severity on the critical follower force, stability and the natural frequency of a beam are investigated.

We study natural convection from a downward pointing cone in a viscoelastic fluid embedded in a porous medium. The fluid properties are numerically computed for different viscoelastic, porosity, Prandtl and Eckert numbers. The governing partial differential equations are converted to a system of fourth order ordinary differential equations using the similarity transformations and then solved together by using the successive linearization method (SLM). Many studies have been carried out on natural convection from a cone but they did not consider a cone embedded in a porous medium with linear surface temperature. The results in this work are validated by the comparison with other authors.

In this paper we show that when the Froude number is less than but close to 1 and the Bond number is greater than but close to 1/3 there exists a new solitary wave solution for surface waves on water with surface tension. An approximate expression for the new solitary wave solution, which satisfies a fourth order ordinary differential equation and represents a wave of depression is presented.