Fourth-order Runge-Kutta Method Research Articles

Large sampling intervals for learning and predicting chaotic systems with reservoir computing

Reservoir computing (RC) is an efficient artificial neural network for model-free prediction and analysis of dynamical systems time series. As a data-based method, the capacity of RC is strongly affected by the time sampling interval of training data. In this paper, taking Lorenz system as an example, we explore the influence of this sampling interval on the performance of RC in predicting chaotic sequences. When the sampling interval increases, the prediction capacity of RC is first enhanced then weakened, presenting a bell-shaped curve. By slightly revising the calculation method of the output matrix, the prediction performance of RC with small sampling interval can be improved. Furthermore, RC can learn and reproduce the state of chaotic system with a large time interval, which is almost five times larger than that of the classic fourth-order Runge–Kutta method. Our results show the capacity of RC in the applications where the time sampling intervals are constrained and laid the foundation for building a fast algorithm with larger time iteration steps.

Peridynamics contact model: Application to healing using phase field theory

This paper presents a contact algorithm in non-ordinary state-based peridynamics (NOSB-PD) framework and its application in healing phase field (HPF) theory. The idea here is to exploit the nonlocality in PD to model the contact between two deformable bodies. The proposed method does not require artificial short-range force or additional degrees of freedom, e.g., Lagrange multiplier, which are often used in the conventional classical and PD contact models. In the regions where contact is likely to take place, neighbor search is performed at every time step and extra particles entering the horizon are identified. These extra particles form artificial bonds which transfer forces only in compression. This feature is essential for rebound simulation. Envisioning surface-to-surface bonding as the opposite process of cracking, we propose an HPF theory to model high velocity impact induced bonding between two material bodies. Here an HPF based bonding criterion is postulated and integrated with the proposed contact algorithm. The solution of the governing equations is obtained using the fourth order Runge-Kutta method. The efficacy of the proposed contact algorithm is demonstrated by simulations of impact and rebound of two plates, two cylinders, and cylinder with rectangular prism and comparing the energy time histories with finite element solutions obtained using ANSYS®. The potential of the proposed HPF theory is demonstrated through simulation of impact of two plates where the plates rebound when the impact velocity is below a critical value, and above that velocity, the plates join with each other.

Effects of Arrhenius Activation Energy on Thermally Radiant Williamson Nanofluid Flow Over a Permeable Stretching Sheet with Viscous Dissipation

This paper explores the role of viscous dissipation, Arrhenius activation energy, and thermal radiation of Williamson nanofluid flow over a permeable stretching sheet. The governing partial differential equations have been simplified through a transformation process, resulting in a set of non-linear differential equations. To find a solution for these equations, a numerical approach is employed, specifically the fourth order Runge-Kutta method. Additionally, a shooting technique is utilized to enhance the accuracy of the numerical solutions. Overall, the study involves reducing complex equations, solving them numerically, and refining the results through a combination of methods. The study investigates the impact of different physical parameters on key factors like velocity, temperature, skin friction coefficient, nano particle volume fraction, and rates of mass and heat transfer. This study exhibits that activation energy parameter enhances concentration profiles, whereas fitted rate constant shows opposite behavior. The activation energy into heat transfer model allows for the optimization of heat transfer systems utilizing Williamson nano fluids.

Analysis of SEIQV Epidemic Model on the Spread of Covid-19 in Jember Regency

Covid-19 is a disease that attacks the respiratory system caused by infection with SARSCoV-2. Efforts to prevent the spread of Covid-19 by vaccination. The spread of disease can be modeled into a mathematical equation. The model used in this study is the SEIQV Model. The disease spread model is analyzed by finding the equilibrium point and the stability of the equilibrium point as well as the basic reproduction number. In the Covid-19 distribution model, a bifurcation analysis is carried out which is needed to determine changes in stability and changes in the number of equilibrium points. Then perform a numerical solution using the fourth-order Runge-Kutta method which is simulated using MATLAB.

Efficient Numerical Solution of the Sharma-Tasso-Olver Equation Using Radial Basis Function-Pseudo Spectral Method with Comparative Analysis

In this study, we will solve the Sharma-Tasso-Olver equation by utilizing the radial basis function-pseudo spectral approach. Pseudo spectral techniques are renowned for their exceptional precision but have constraints regarding their ability to handle geometric variations. The integration of Radial Basis Functions (RBF) with pseudo spectral methods offers a solution to surmount this particular limitation. To estimate the differentials with respect to the variable , we will use the radial basis functions, including commonly used ones like multiquadric, inverse multiquadric, and inverse quadric, have been employed. Furthermore, the problem has transformed into the ordinary differential equation (ODE), and the solution to this ODE system has been acquired through the application of the fourth-order Runge-Kutta method. Additionally, a comparative analysis has been conducted between the solutions obtained through the suggested technique and the exact solutions.

Energy harvesting for a vibration isolation system with negative stiffness under harmonic excitation

ABSTRACT In this paper, an energy harvesting system with a negative stiffness structure is established, whose main purpose is to achieve good vibration isolation, and at the same time, the system can also achieve the purpose of energy harvesting. In order to analyse the system model, the amplitude–frequency response equations for the displacement, induced voltage and electric power of the system when the system is in the main resonance state are obtained by using the multiscale method, and their amplitude–frequency response curves are plotted. In addition, we obtained the existence conditions for the system to satisfy Smale's horseshoe chaos using the Melnikov function. And the phase diagram and Poincaré mapping of the system were obtained using the fourth-order Runge–Kutta method. Numerical results show that the system has some energy harvesting performance along with vibration isolation. Variations in the parameters of the system affect the resonance frequency band and co-oscillation amplitude of the system as well as the nonlinear characteristics of the system itself.

Nonlocal strain gradient-based geometrically nonlinear vibration analysis of double curved shallow nanoshell containing functionally graded layers

It is a fact that while examining the mechanical behavior of nano/microscopic materials, size effects become apparent. This study aims at presenting a new approach to geometrically nonlinear vibration analysis of double curved shallow nanoshell on Kerr foundation containing functionally graded layers under the impacts of harmonic loading in the framework of nonlocal strain gradient theory and the classical shell theory. The three-parameter Kerr model is founded on the independence of the upper and lower elastic layers related to the shear layer. The motion equations of the double curved shallow nanoshell are deduced by introducing the Von Kármán strain-displacement relations. Governing equations are solved by displacement function method and Galerkin method to achieve the deflection differential equation, the fourth-order Runge-Kutta methods are implemented to solve the deflection differential equation of the dynamic system. In this study, the obtained outcomes are compared with other publications indicating that the theoretical results agree with the previously published ones. Afterward, the article investigates the effects of the material length scale parameter, the nonlocal parameter, the geometric parameters and the material properties on the geometrically nonlinear vibration analysis of the double curved shallow nanoshell containing functionally graded layers. This research paper can provide extensive references and valuable guidelines for the impact and application of double curved shallow nanoshell structures containing functionally graded layers.

SPECTRAL METHOD FOR ONE DIMENSIONAL BENJAMIN-BONA-MAHONY-BURGERS EQUATION USING THE TRANSFORMED GENERALIZED JACOBI POLYNOMIAL

The Benjamin-Bona-Mahony-Burgers equation (BBMBE) plays a fundemental role in many application scenarios. In this paper, we study a spectral method for the BBMBE with homogeneous boundary conditions. We propose a spectral scheme using the transformed generalized Jacobi polynomial in combination of the explicit fourth-order Runge-Kutta method in time. The boundedness, the generalized stability and the convergence of the proposed scheme are proved. The extensive numerical examples show the efficiency of the new proposed scheme and coincide well with the theoretical analysis. The advantages of our new approach are as follows: (i) the use of the transformed generalized Jacobi polynomial simplifies the theoretical analysis and brings a sparse discrete system; (ii) the numerical solution is spectral accuracy in space.

Innovative coupling of s-stage one-step and spectral methods for non-smooth solutions of nonlinear problems

The behavior of nonlinear dynamical systems arising in mathematical physics through numerical tools is a challenging task for researchers. In this context, an efficient semi-spectral method is proposed and applied to observe the robust solutions for the mathematical physics problems. Firstly, the space variable is approximated by the Vieta-Lucas polynomials and then the s-stage one-step method is applied to discretize the temporal variable which transfers the problem in the form Cn+1=Cn+Δtϕ(x,t,Cn,F(un)). Novel operational matrices of integer order are developed to replace the spatial derivative terms presented in the discussed problem. Related theorems are included in the study to validate the approach mathematically. The proposed semi-spectral schemes convert the considered nonlinear problem to a system of linear algebraic equations which is easier to tackle. We also accomplish an investigation on the error bound and convergence to confirm the mathematical formulation of the computational algorithm. To show the accuracy and effectiveness of the suggested computational method numerous test problems, such as the advection-diffusion problem, generalized Burger-Huxley, sine-Gordon, and modified KdV–Burgers’ equations are considered. An inclusive comparative examination demonstrates the currently suggested computational method in terms of credibility, accuracy, and reliability. Moreover, the coupling of the spectral method with the fourth-order Runge-Kutta method seems outstanding to handle the nonlinear problem to examine the precise smooth and non-smooth solutions of physical problems. The computational order of convergence (COC) is computed numerically through numerous simulations of the proposed schemes. It is found that the proposed schemes are in exponential order of convergence in the spatial direction and the COC in the temporal direction validates the studies in the literature.

Melting heat transfer and irreversibility analysis in Darcy-Forchheimer flow of Casson fluid modulated by EMHD over cone and wedge surfaces

The present mathematical model analyses the melting heat transfer and irreversibility in Darcy-Forchheimer flow of Casson fluids modulated by the electroosmosis and magnetohydrodynamic mechanisms over the wedge and cone surfaces. The effects of thermal radiation, thermal buoyancy and heat generation/absorption are taken into consideration. Appropriate similarity transformations are utilized to establish a system of ordinary differential equations which are non-linear. The model has been numerically simulated using shooting approach in conjunction with the fourth order Runge-Kutta method under the bvp4c tool of MATLAB. Numerical results are computed for fluid velocity, fluid temperature, Nusselt number and Bejan number for analysing the fluid flow behaviour and heat transfer characteristics for the Casson fluid. A comparative analysis is performed to determine the irreversibility and heat transfer in the absence and presence of melting parameters, respectively. Variations in skin friction coefficient, local Nusselt number and Bejan number have also been examined under the effects of key parameters. Key findings of the present study indicate enhanced velocity distribution with increasing zeta potential and electric field parameters, while it diminishes with higher Prandtl number and Forchheimer terms. Moreover, fluid temperature reduces with rising zeta potential, while the Bejan number escalates with greater thermal radiation in the core region of geometries. The findings of the model may be useful in various applications of thermal systems.

The Analysis and the Solution of Incubation Period in a Disease Model

This study deals with the analysis and the solution of incubation period in a disease model by adopting the mathematical model with incubation period of diseases and the mathematical model without the incubation period of diseases. In the model equations, we partitioned the population into Susceptible (S), Incubated (I), Infected (D) population. We have compared the model equations without incubation period with the model equation with incubation period by solving and incorporating the system of first order linear equations into fourth order Runge-kutta method which has better error accuracy for solving first order equations. Graphical results for incubation class show that the infectious diseases were fatal if immediate attention is not given to endemic villages and communities.Keywords: SID Model, Incubation period, Runge-kutta method, numerical simulation, transmission.

Linear and non-linear stability analysis of double-diffusive Hadley-Prats flow through horizontal porous layer with non-uniform thermal and solutal gradients

A scientific model is devised to understand the onset of double-diffusion Hadley-Prats flow through an infinite parallel permeable channel with non-uniform thermal, solutal gradients, and concentration based internal heat generation under the influence of convection conditions is considered. Darcy’s model is implemented in the context of a porous medium, which is characterized as isotropic and homogeneous. A linear and nonlinear stability analysis is conducted with the help of energy functional method and longitudinal roll perturbations are examined. The emerged dimensionless eigenvalue problem in the linear and nonlinear stability analysis is solved numerically by utilizing fourth order Runge-Kutta method (RK-4) method and shooting scheme. The critical values of wave and vertical thermal Rayleigh numbers are examined. A comprehensive analysis of solutions concerning the onset of convection mechanism is conducted. The produced results are discussed for the various values of emerged physical parameters. It is evident that, the presence of non-uniform thermal and solutal gradients and the internal heat generation mechanism causes the significant variations in the critical values of vertical thermal Rayleigh number in the channel regime.

ANALYSIS AND OPTIMAL CONTROL OF HIV GROWTH MODEL IN THE BODY WITH ANTIRETROVIRAL THERAPY

Human Immunodeficiency Virus (HIV) is a virus that affects the human immune system. HIV infection causes a decrease in the body's immunity because the virus attacks immune-building cells, especially T-CD4 cells. Currently, there is no treatment that can cure or eliminate HIV, but antiretroviral therapy can be done. This study discusses the growth model of HIV in the body that is given control in an effort to maximize healthy T-CD4 cells. In this model, the infection-free and infected equilibrium points are also discussed and their stability is analyzed. Then the optimal control is solved using the Pontryagin Maximum Principle method and solved numerically using the fourth-order Runge Kutta method. Based on the analysis and simulation results, the system is asymptotically stable around the infection-free equilibrium point and unstable around the infected equilibrium point. Simulation results show that with the control of antiretroviral therapy, the T-CD4 cell population grows significantly which can improve the quality of life of patients. And the growth of HIV in the body can be inhibited until it cannot reproduce itself.

Numerical Investigation of Diffusion thermo and Thermal Diffusion on Mixed Convection flow of Williamson Nanofluid through a vertical cone with porous material in the presence Thermal radiation and Activation Energy

In this paper, analyze the impact of Activation energy, Diffusion thermo and thermal diffusion with heat and mass transfer inherent of over a vertical cone through Thermally Radiant Williamson nanofluid with saturated porous medium under the convective boundary condition in the presence of thermal radiation has been studied. The coefficients of Brownian and thermophoresis diffusions are also taken into consideration. The governing partial differential equations are reduced to a couple of nonlinear ordinary differential equations by using suitable transformation equations; these equations are then solved numerically with the use of the conventional fourth-order Runge Kutta method accompanied by the shooting technique. As a result, the effects of various physical parameters on the velocity, temperature, and nanoparticle concentration profiles as well as on the skin friction coefficient and rate of heat transfer are discussed with the aid of graphs and tables. This study has been directly applied in the pharmaceutical industry, microfluidic technology, microbial improved oil recovery, modelling oil and gas-bearing sedimentary basins, and many other fields. Further, to check the accuracy and validation of the present results, satisfactory concurrence is observed with the existing literature.

Ultra-low frequency vibration energy harvesting of piezoelectric vibration systems with an adjustable device

Energy harvesting from external Ultra-low frequency vibrations is currently attracting a lot of attention. However, for the same system, it is difficult to satisfy the need for more efficient energy harvesting in a variety of external environments. Accordingly, an adjustable gearing device connected to the energy harvesting system is used herein for adapting energy harvesting to a wide range of external environments. The device uses four symmetrical inclined springs. One end of the springs is connected to a gear and the other end to a mass element. The geometrical parameters of the system can be varied by rotating the gears, thus varying the dynamic stiffness of the system. The kinetic equations of the system are derived and their nonlinear terms are Taylor-expanded. The amplitude-frequency response curve of the system is obtained by solving using the harmonic balance method. Multi-color plots of RMS values of induced voltage and electric power are obtained using the fourth-order Runge-Kutta method. The results show that the device has excellent adjustable performance. Adjusting the rotation angle of the gears can change the soft and hard characteristics of the system's dynamic stiffness and the intrinsic frequency, thus improving the system's energy harvesting performance in different environments.

A general and efficient harmonic balance method for nonlinear dynamic simulation

This paper presents a general and efficient harmonic balance (GE-HB) method to expeditiously evaluate the periodic responses of nonlinear systems. By employing the fast Fourier transform (FFT) technique, it rapidly computes harmonic expansion coefficients of the nonlinear terms and their corresponding constant tensor matrices. Notably, the nonlinear part of the Jacobian matrix is just the product of the coefficients and the constant tensor matrices, loop computations are no longer requisite for calculating the Jacobian matrix. This, facilitates accelerated determination of the nonlinear portion of the Jacobian matrix, thus obtaining the periodic responses of the system in an efficient manner. Furthermore, the constant tensor matrices are exclusively contingent upon the harmonic bases of the system's periodic response, devoid of reliance on the equation's specific formulation, facilitates obtaining the Jacobian matrix in a general manner. Simulation results on a 20-degree-of-freedom nonlinear oscillator chain system demonstrate that the efficiency of the GE-HB method is about tenfold that of the harmonic balance-alternating frequency/time domain (HB-AFT) method and twenty-fivefold that of the fourth-order Runge-Kutta (RK4) method. Further, the simulation of a 284-degree-of-freedom nonlinear rotor system exhibits the GE-HB method's computational efficiency surpassing HB-AFT method by over 500 times and Newmark method by over 3000 times. Moreover, in these two examples, the GE-HB method captures all solution branches, including unstable periodic solutions, providing comprehensive insights into periodic responses of the systems. The GE-HB method proposed here has the potential to explore dynamic characteristics for high-dimensional complex nonlinear systems in engineering.

Design and optimization of a two-degrees-of-freedom single-sided vibro-impact nonlinear energy sink for transient vibration suppression of a thin plate

To rapidly suppress the transient vibration response of a thin plate, a two-degrees-of-freedom single-sided vibro-impact nonlinear energy sink (TSSVI NES) is proposed in this paper. Based on the Kirchhoff hypothesis, the coupled equations of motion between the plate and the TSSVI NES are established. The system is modeled by the Galerkin technique and numerically solved using the fourth order Runge-Kutta method. To obtain the optimal parameters of the TSSVI NES, a particle swarm optimization method is adopted. The system response with the TSSVI NES is compared with that of the cubic nonlinear energy sink (CNES) and the single-sided vibro-impact nonlinear energy sink (SSVI NES), showing that the vibration suppression ability of the TSSVI NES is better. The mechanism of rapid response suppression of the TSSVI NES is demonstrated by way of simulations. A parametric study is also conducted to determine the effects of damping, the coefficient of restitution, mass, stiffness, and installation position of TSSVI NES on the efficacy of the vibration control device. The numerical results demonstrate that the performance of the TSSVI NES is robust to the damping, coefficient of restitution, and installation position, but sensitive to mass difference and stiffness difference. Finally, the effect of the TSSVI NES is examined when the plate is excited by seismic loads, demonstrating a good response suppression performance.

Effect of Interstage Torsional Damping on Nonlinear Dynamic Characteristics of Two-Stage Planetary Gear System

Abstract This paper explores the influence of interstage torsional damping on the nonlinear dynamics of a two-stage planetary gear transmission system. A comprehensive nonlinear dynamic model is developed, incorporating interstage torsional damping, stiffness, time-varying meshing parameters, damping, and tooth side clearance. The nonlinear equations are derived and solved using the fourth-order Runge–Kutta method. The investigation reveals that the changes in torsional damping between stages significantly affect the nonlinear dynamic characteristics of the system, especially under load fluctuations. Bifurcation characteristics are analyzed using phase diagrams, Poincaré diagrams, time history diagrams, and frequency–domain spectrograms. The simulation results demonstrate that heightened torsional damping between stages distinctly affects the motion states across excitation frequencies. Notably, it mitigates the unstable motion state caused by high-frequency excitation when tooth side clearance is minimal.

Study of the Six-Compartment Nonlinear COVID-19 Model with the Homotopy Perturbation Method

The current study aims to utilize the homotopy perturbation method (HPM) to solve nonlinear dynamical models, with a particular focus on models related to predicting and controlling pandemics, such as the SIR model. Specifically, we apply this method to solve a six-compartment model for the novel coronavirus (COVID-19), which includes susceptible, exposed, asymptomatic infected, symptomatic infected, and recovered individuals, and the concentration of COVID-19 in the environment is indicated by S(t), E(t), A(t), I(t), R(t), and B(t), respectively. We present the series solution of this model by varying the controlling parameters and representing them graphically. Additionally, we verify the accuracy of the series solution (up to the (n−1)th-degree polynomial) that satisfies both the initial conditions and the model, with all coefficients correct at 18 decimal places. Furthermore, we have compared our results with the Runge–Kutta fourth-order method. Based on our findings, we conclude that the homotopy perturbation method is a promising approach to solve nonlinear dynamical models, particularly those associated with pandemics. This method provides valuable insight into how the control of various parameters can affect the model. We suggest that future studies can expand on our work by exploring additional models and assessing the applicability of other analytical methods.

Generalized Shannon entropy sparse wavelet packet transform for fault detection of traction motor bearings in high-speed trains

An effective structural health monitoring method of traction motor bearings is a powerful guarantee for the safety operation of high-speed trains. However, it is exceptionally difficult to detect bearing fault characteristics from the vibration signals of traction motor bearings operating at high rotational speeds. In this scenario, a generalized Shannon entropy sparse wavelet packet transform (GSWPT) for fault detection of motor bearings is proposed in this paper. Firstly, a generalized Shannon entropy sparse regularization method is proposed to obtain sparse wavelet reconstruction coefficients by extending the definition of the Shannon information entropy, and the non-convex sparse regularization function is minimized by synergistic swarm optimization algorithm. Then, the wavelet node coefficients are weighted according to the second-order cyclostationarity index of the wavelet packet node to further enhance the sparsity of the reconstructed signal. Moreover, the optimal decomposition level of GSWPT is adaptively selected by the maximum sparsity and cyclostationarity criterion. Particularly, in order to verify the bearing fault detection performance of GSWPT in practical engineering, a bearing fault dynamic model of traction motor in high-speed train was established based on Hertz contact theory and the fourth-order Runge-Kutta method to obtain simulated data under strong Gaussian white noise, and a corresponding test platform was constructed to collect experimental data under different operating conditions. Finally, the applications on the simulated and experimental signals of traction motor bearings in high-speed trains demonstrate that GSWPT significantly outperforms the conventional wavelet packet transform, dual-tree complex wavelet packet transform, blind deconvolution, modal decomposition, and Infogram methods to some extent for fault detection.