Second-order Ordinary Differential Equations Research Articles

Solving Nonlinear Second-Order ODEs via the Eisenhart Lift and Linearization

The linearization of nonlinear differential equations represents a robust approach to solution derivation, typically achieved through Lie symmetry analysis. This study adopts a geometric methodology grounded in the Eisenhart lift, revealing transformative techniques that linearize a set of second-order ordinary differential equations. The research underscores the effectiveness of this geometric approach in the linearization of a class of Newtonian systems that cannot be linearized through symmetry analysis.

Comparative analysis of thermal behaviour in functionally graded material fin under convective-radiative heat transfer: Homotopy perturbation method

In this study, the thermal characteristics of a fully wet, moving longitudinal fin made from linear Functionally Graded Material (FGM) is explored. The fin response to convective-radiative heat transfer is examined through three distinct FGM scenarios namely, a Homogeneous material (HM), Functionally Graded Material I (FGM I) having greater thermal grading near the base, and Functionally Graded Material II (FGM II) having increased thermal grading towards the tip. The derived energy equation is a second-order nonlinear ordinary differential equation is solved using the Homotopy Perturbation Method (HPM) compared with results found in the existing literature. The fin thermal profile has been graphically analysed for the thermal conductivity grading parameter (β), Peclet number (Pe) and other pertinent parameters. It is found that the temperature distribution along the fin is highest for the FGM II fin and lowest for the HM fin, with the FGM I fin displaying intermediate values. Efficiency of the fin (η) is discovered to be lowest in the case of HM fin structures followed by FGM I and FGM II fin structures at the base and vice versa at the tip. The results show significant potential of HPM which benefits the manufacture and design of FGM fin structures.

Mathematical model of a lightweight three-axle off-road vehicle construction for Arctic zone of Russia

Introduction. The model and results of calculating the smoothness of a light three-axle off-road vehicle for the Arctic zone of Russia are considered. The model on standard approaches and uses a system of assumptions that limits the number of degrees of freedom for the vehicle body to three, as well as one degree of freedom for the unsprung masses is based. The mathematical model is a system of ordinary differential equations and is supplemented with the necessary algebraic equations, as well as initial conditions. The system is integrated using the 4th order RungeKutta method, for which a program was written in C++. The calculations presented in the article demonstrate the possibility of conducting research on the smoothness of a vehicle under conditions of arbitrary terrain, typical for off-road conditions in the winter conditions of the Arctic zone. The dimensions and other parameters of the vehicle were taken from a full-scale model that was used in real expeditions in 2003-2019. Based on the model, suspension characteristics will be developed for a new all-terrain vehicle.Theory. When operating a wheeled vehicle in a wide range of conditions, even in northern regions, transverseangular vibrations are very often insignificant, so only vertical linear and longitudinal-angular vibrations of the frame can be considered. This problem enables to construct a system of equations of vehicle motion using selected degrees of freedom. From a mathematical point of view, these equations are classified as second-order ordinary differential equations with a variable structure of the right-hand sides, which reflects the non-linear nature of the behavior of the suspension in terms of its geometric constraints.Methods. The work uses numerical methods to solve the equations of the constructed model, which enables to gradually weaken the accepted assumptions and build more general calculation algorithms. The main integration method for ensuring the stability of solutions is the multi-step Adams method, which, with the correct choice of step, ensures the necessary stability of the solution over sufficiently long model times. However, in this work, the 4th order Runge-Kutta method was adopted, which turned out to be quite sufficient.Results and conclusions. The paper presents the results of a numerical study of the oscillatory processes of an off-road vehicle during uniform translational motion of the vehicle on a horizontal surface with a given profile of irregularities. The graphs show a transition process of oscillations, which ends with reaching a steady state. The shape of oscillations in a steady state can be irregular and significantly depends on the given speed of the all-terrain vehicle. Analysis of the dependencies presented in the figures shows that the shape of the oscillations of the allterrain vehicle’s frame, as well as the amplitude and frequency, significantly depend on the speed of the vehicle (at a constant road profile). Changing the road profile leads to corresponding changes in the forms and characteristics of forced vibrations of a vehicle on a suspension, which makes it possible to build the necessary amplitude-frequency characteristics, optimize the elastic and dissipative parameters of suspensions, optimize their number and location, and also monitor the movements of arbitrary points at which various units are located.

A High-Resolution Exponential Spline Method and Its Convergence Analysis for Two-Point Mixed Boundary Value Problems

In this paper, we present a new numerical method for effectively solving general second-order ordinary differential equations with mixed boundary conditions. Our approach utilizes a quasi-variable mesh enabling us the flexibility to adapt the mesh density according to different boundary layer problems. By employing a discretization technique that incorporates the construction of exponential spline, we achieve third-order accuracy at internal grid points, while boundary points exhibit fourth-order accuracy. Since the method is based on off-step grid points, it eliminates the need to modify the method while employing it to singular problems. We provide computational results to various numerical problems which arise in different physical phenomena like Burger’s equation and Sturm–Liouville equation. A comparison with recent findings underscores the superior performance of our method.

Finding and implementing the numerical solution of an optimal control problem for oscillations in a coupled objects system

Finding and implementing the numerical solution of an optimal control problem for oscillations in a coupled objects system

On Periodic Solutions of a Second-Order Ordinary Differential Equation

On Periodic Solutions of a Second-Order Ordinary Differential Equation

One-step block scheme with optimal hybrid points for numerical integration of second-order ordinary differential equations

In this paper, a one-step block of optimized hybrid schemes for the numerical integration of second-order initial value problems (IVP) of ordinary differential equations (ODE) is constructed via collocation techniques. The developed scheme is obtained by considering two intra-step nodal points as hybrid points, which are chosen in order to achieve optimized errors of the main formulae approximating the solution such that 0 < v1 < v 2 < 1 where v1 and v2 are defined as hybrid points. The characteristics of the developed scheme are analyzed. Application of the new scheme on some second-order IVPs shows the accuracy and effectiveness of the scheme compared with some existing methods.

Double shooting method for FRCM reinforced systems in debonding problems

The debonding process in FRCM (Fiber Reinforced Cementitious Matrix) strengthened brittle substrates, such as concrete and masonry, is characterized by complex failure modes that include debonding at both the external and internal interfaces of the composite, as well as damage to the mortar matrix. This study addresses this problem classically as a mode-II fracture process by accounting for the one-dimensional non-linear interfacial shear stress-slip relationship combined with matrix axial non-linearity. The approach employed considers the longitudinal displacements of both the inner and outer layers of mortar, along with that of the fiber textile, as unknowns, and analytically solves the second-order ordinary differential equation (ODE) system that governs the debonding problem. Then, the boundary value problem (BVP) is transformed into an initial value problem (IVP) and a double shooting method is proposed to find displacements and stresses along the bonded length in the presence of non-linearity. The appropriate initial values of matrix layer displacements, satisfying the required boundary conditions, are efficiently determined through a two-dimensional bisection method working on rectangular and triangular patches. The softening behavior of both mortar layers and interfaces between mortar and fiber are taken into account using jagged constitutive relationships to enhance convergence. The numerical calculation speed, robustness, and key simulation parameters are thoroughly discussed. Moreover, the proposed numerical approach is compared with existing experimental data and models, on both concrete and masonry FRCM-reinforced specimens. The influence of mortar failure and interfacial strength is also explored. The model demonstrates its ability to effectively reproduce the global bond behavior observed in experimental studies, while capturing the local behavior and simulating various failure modes observed in the tests.

A Novel Approach to Linear and Nonlinear Time-History Analysis of Structures: Gauss–Lobatto–Hermite 4-Point (GLH-4P) Method

AbstractAn efficient time integration scheme has been devised for conducting time-history dynamic analysis of single-degree-of-freedom (SDOF) structural systems. The developed scheme belongs to Runge-Kuta family. The primary innovation of the new formulation lies in the constructive cooperation between the Hermite interpolation and the Gauss–Lobatto integration, linked to address the second-order ordinary differential equation of motion (DEOM). To put the idea into practice, several Hermite interpolators are generated to exploit the information of the internal points of Gauss-Lobatto integrator. The 4-point Gauss–Lobatto formula is employed to numerically integrate acceleration and velocity functions. This novel methodology has been designated the Gauss–Lobatto–Hermite 4-point (GLH-4P) method. GLH-4P furnishes a generally applicable implicit algorithm for conducting both linear and nonlinear analyses of structures subjected to seismic excitation. The advantage of the proposed algorithm over the conventional techniques is its higher accuracy level. Besides, it offers better stability compared with the conventional methods such as the Newmark-β and Wilson-θ methods. Unlike the other methods, GLH-4P can handle high-frequency systems without reducing the step size. Numerical examples reveal the efficacy of the GLH-4P method when juxtaposed against existing methodologies.

Novel approach for solving higher-order differential equations with applications to the Van der Pol and Van der Pol–Duffing equations

BackgroundNumerical methods are used to solve differential equations, but few are effective for nonlinear ordinary differential equations (ODEs) of order higher than one. This paper proposes a new method for such ODEs, based on Taylor series expansion. The new method is a second-order method for second-order ODEs, and it is equivalent to the central difference method, a well-known method for solving differential equations. The new method is also simple to implement for higher-order differential equations. The proposed technique was applied to solve the Van der Pol and Van der Pol–Duffing equations. It is stable over a wide range of nonlinearity and produces accurate and reliable results. For the self-excitation Van der Pol equation, the proposed technique was applied with different values of nonlinear damping.ResultsThe results were compared with those obtained using the ODE15s solver in MATLAB. The two sets of results showed excellent agreement. For the forced Van der Pol–Duffing equation, the proposed technique was applied with different values of exciting force amplitude and frequency. It was found that for certain conditions, the solution obtained using the proposed technique differed from that obtained using ODE15s.ConclusionsThe solution obtained using the proposed technique showed good agreement with the solutions obtained using ODE45 and Runge–Kutta fourth order. The results show that the proposed approach is very simple to apply and produces acceptable error. It is a powerful and versatile tool for solving of high-order nonlinear differential equations accurately.

A systematic approach to obtain the analytical solution for coupled linear second-order ordinary differential equations: Part II

A systematic approach to obtain the analytical solution for coupled linear second-order ordinary differential equations: Part II

Integrable geodesic flows and metrisable second-order ordinary differential equations

It is well known that the system of ordinary differential equations (ODEs) describing geodesic flows of some Riemannian metrics on 2-surfaces admits a projection on a special class of second-order ODEs. In this paper we study in detail this special class of ODEs. We classify all such autonomous ODEs possessing autonomous first integrals that are fractional-quadratic in the first derivative. We construct all the families of integrable geodesic flows related to the ODEs under study. These families are parameterized by two arbitrary functions and contain metrics with superintegrable geodesic flows. We also study the inverse problem and find novel families of second-order ODEs that are related to integrable geodesic flows. We explicitly present first integrals of such ODEs. We find the general structure of metrisable second-order ODEs related to conformal metrics.

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.

Soret and dufour effects on unsteady non-linear mixed convection flow past a stretching sheet influenced by non-linear thermal radiation

This article presents a novel mathematical model to analyze the behavior of a time-dependent, incompressible, non-linear (quadratic) mixed convection flow across a stretching surface heated by convection. The model takes into account the influence of non-linear thermal radiation, along with the Soret and Dufour effects. The equations governing the flow characteristics in the boundary layer approximation are derived and transformed into ordinary differential equations (ODEs) using similarity variables. The second-order ODEs are transformed into an initial value problem, which is then numerically solved using the bvp4c solver in MATLAB. The novel findings of the study have been presented using graphs and tables. The intensification of non-linear convection parameters has been shown to lead to the development of velocity near the sheet while simultaneously reducing velocity toward the free stream. The non-linear convection parameters reduce the temperature distribution, while the temperature ratio and the non-linear thermal radiation parameters enhance it. The non-linear convection parameters ( 0.5 ≤ α 1 , α 2 ≤ 2.0 ) , temperature ratio ( 1.2 ≤ θ w ≤ 1.5 ) , and non-linear thermal radiation ( 1.0 ≤ Rd ≤ 1.5 ) , parameters are influential factors in decreasing the surface drag force (by about 15.4%, 28.5%, 4.4%, and 4.3%) and increasing the cooling rate (by about 0.4%, 0.5%, 19.5%, and 23%), respectively. The Biot ( 0.5 ≤ Bi ≤ 1.0 ) number is also found to significantly contribute toward enhancing the cooling rate of the system by about 51.6%.

A note on Laplace Adomian decomposition method on a linear problem

In last two decades, Laplace Adomian Decomposition Method (LADM) is vastly used to solve non-linear (or even fractional order) differential equations. The method approaches the solution with the partial sums of function series. However, it is not easy to show that the limit of the function series is the exact solution of the problem. In this article, we consider a simple problem such as homogenous second-order linear ordinary differential equation with constant coefficients. We prove analytically that the LADM gives the right exact solution to the considered problem.

Scalar quasi-normal modes of accelerating Kerr-Newman-AdS black holes

We study linear scalar perturbations of slowly accelerating Kerr-Newman-anti-de Sitter black holes using the method of isomonodromic deformations. The conformally coupled Klein-Gordon equation separates into two second-order ordinary differential equations with five singularities. Nevertheless, the angular equation can be transformed into a Heun equation, for which we provide an asymptotic expansion for the angular eigenvalues in the small acceleration and rotation limit. In the radial case, we recast the boundary value problem in terms of a set of initial conditions for the isomonodromic tau function of Fuchsian systems with five regular singular points. For the sake of illustration, we compute the quasi-normal modes frequencies.

Momentum transfer across an open-channel, turbulent flow

The distribution of stress generated by a turbulent flow matters for many natural phenomena, of which rivers are a prime example. Here, we use dimensional analysis to derive a linear, second-order ordinary differential equation for the distribution of stress across a straight, open channel, with an arbitrary cross-sectional shape. We show that this equation is a generic first-order correction to the shallow-water theory in a channel of large aspect ratio. It has two adjustable parameters – the dimensionless diffusion parameter, $\chi$ , and a local-shape parameter, $\alpha$ . By assuming that the momentum is carried across the stream primarily by eddies and recirculation cells with a size comparable to the flow depth, we estimate $\chi$ to be of the order of the inverse square root of the friction coefficient, $\chi \sim C_f^{-1/2}$ , and predict that $\alpha$ vanishes when the flow is highly turbulent. We examine the properties of this equation in detail and confirm its applicability by comparing it with flume experiments and field measurements from the literature. This theory can be a basis for finding the equilibrium shape of turbulent rivers that carry sediment.

On the Exact Solution of a Scalar Differential Equation via a Simple Analytical Approach

The existence of the advance parameter in a scalar differential equation prevents the application of the well-known standard methods used for solving classical ordinary differential equations. A simple procedure is introduced in this paper to remove the advance parameter from a special kind of first-order scalar differential equation. The suggested approach transforms the given first-order scalar differential equation to an equivalent second-order ordinary differential equation (ODE) without the advance parameter. Using this method, we are able to construct the exact solution of both the transformed model and the given original model. The exact solution is obtained in a wave form with specified amplitude and phase. Furthermore, several special cases are investigated at certain values/relationships of the involved parameters. It is shown that the exact solution in the absence of the advance parameter reduces to the corresponding solution in the literature. In addition, it is declared that the current model enjoys various kinds of solutions, such as constant solutions, polynomial solutions, and periodic solutions under certain constraints of the included parameters.

Implicit Block and Runge- Kutta type Methods for Solution of Second-Order Ordinary Differential Equations.

In this research, implicit discrete schemes which form our block integrators were developed for solving Initial Valued Problems of Ordinary Differential Equations. The equivalent second order Runge-Kutta type Methods (RKTM) were also constructed for the same purpose. Both methods were demonstrated on linear and nonlinear problems of Ordinary Differential Equations. Numerical results obtained from RKTM show that the method is competitive with the existing one.

Stability analysis of nonlinear convective boundary layer flow on a shrinking surface in the presence of viscous dissipation

This study investigates analytically the stability analysis of a nonlinear convective boundary layer (BL) flow of the nanofluids GO-EG and GO-W on a shrinking surface in the presence of a magnetic field and viscous dissipation. We provide a second-order nonlinear ordinary differential equation (NODE) for temperature distribution (TD) and a third-order NODE for velocity profile (VP) coupling based on the thermophysical characteristics of the base fluid and nanofluid as well as similarity transformations (STs) in the fundamental governing equations of momentum and energy. The analytical method HAM is used to answer the equations that have been gathered. In many different industries, such as manufacturing, automotive, microelectronics, and defense, cooling of various types of equipment and devices has very important challenges. To fulfill the demands of the engineering and industrial industries, it is hoped that this issue would significantly improve the heat transformation ratio.