Nonlinear Algebraic Equations Research Articles

An intriguing numerical strategy for Zakharov–Kuznetsov equation through graph-theoretic polynomials

This paper explores graph-theoretic polynomials to find the approximate solution of the (2+1)D Time-fractional Zakharov-Kuznetsov(TF-Z-K) equation. The Zakharov-Kuznetsov equations govern the behavior of nonlinear acoustic waves in the plasma of hot isothermal electrons and cold ions in the presence of a homogeneous magnetic field. Independence polynomials of the Ladder-Rung graph serve as the polynomial approximation for the suggested Independence Polynomial Collocation Method (IPCM). The Caputo fractional derivatives are adopted to determine the fractional derivatives in the TF-Z-K equation. The TF-Z-K equation is converted into a system of nonlinear algebraic equations using the collocation points in IPCM. The Newton-Raphson approach yields the solution of the suggested method by solving the resulting system. We’ve compared a few scenarios with the tangible outcomes to validate the procedure. Quantitative outcomes match the current findings and validate the exactness of IPCM compared t o the recent numerical and semi-analytical approaches in the literature.

Semi-analytical framework for modelling of surface heating and its impact on nonlinear stability of nanotube reinforced doubly curved porous fiber composite panels

In elevated temperatures, the stiffness of the panel deteriorates and can lead to loss of stability of structures. Hence, the stability analysis of thin-walled structures is a critical investigation in the thermal environment. Finding such criticality, the present semi-analytical study investigated the nonlinear stability behaviour of porous doubly curved thin-walled panels subject to surface heating like dome heating and localised heating. To improve the stiffness of the panel, the matrix is reinforced with carbon nanotubes (CNTs). Improper mixing of such nanoparticles can lead to agglomeration or bundling effect, which may reduce the structure's stiffness. Still, its investigation in surface heating can be an essential aspect that is found to be untouched by researchers for porous shell panels. The effective material properties of the three-phase composite panel are determined using the Eshelby-Mori-Tanaka (E-M-T) approach and the Chamis homogenisation technique. Using the variational principle, governing equations are derived and further simplified to non-linear algebraic equations using the Galerkin technique. Gibson and Ashby foam model is used to model cell structure. It is observed that a complete agglomerated state results in higher post-buckling strength, a continuous deformation path due to biaxial compression and curved geometry, and porosity of cellular structure distribution and cell walls are the important parameters to study the stability of panel in a thermal environment are the major conclusions drawn from the current study. Current investigation can be an important contribution towards modelling aircraft structure in supersonic airflow and furnace walls of thermal power plants subjected to localised heating.

Fibonacci Wavelet Collocation Method for Solving Dengue Fever SIR Model

The main focus in this manuscript is to find a numerical solution of a dengue fever disease model by using the Fibonacci wavelet method. The operational matrix of integration has been obtained using Fibonacci wavelets. The proposed method is called Fibonacci wavelet collocation method (FWCM). This biological model has been transformed into a system of nonlinear algebraic equations by using the Fibonacci wavelet collocation scheme. Afterward, this system has been solved by using the Newton–Raphson method. Finally, we provide evidence that our results are better than those obtained by various current approaches, both numerically and graphically, demonstrating the method’s accuracy and efficiency.

A spectral collocation scheme for solving nonlinear delay distributed-order fractional equations

The basic purpose of this paper is to implement a variant of the Jacobi–Gauss collocation scheme for solving nonlinear delay (in difference format) fractional differential equations of distributed-order type. At the first stage of the solving procedure, the Legendre–Gauss quadrature rule is implemented to approximate the main problem by a difference delay fractional differential equation of multi-term type. So, one can apply the Jacobi–Gauss collocation approach for localizing the resulting multi-term fractional differential equation with the difference delay factor. Moreover, some special cases of Jacobi–Gauss quadrature rules can be used for approximating the involved integrals in the aforementioned problem. By considering the Lagrange interpolating polynomials, as the numerical solutions, the solution of the main problem can be approximated by solving the associated system of nonlinear algebraic equations in terms of unknown Lagrange multipliers. Convergence analysis of the problem is investigated rigorously under some mild conditions at a spectral rate. Extensive test problems are considered to justify the theoretical analysis.

An Effective Matrix Technique for the Numerical Solution of Second Order Differential Equations

In this paper, an effective technique for solving differential equations with initial conditions is presented. The method is based on the use of the Legendre matrix of derivatives defined on the close interval [-1,1]. Properties of the polynomial are outlined and further used to obtain the matrix of derivative which was used in transforming the differential equation into systems of linear and nonlinear algebraic equations. The systems of these algebraic equations were then solved using Gaussian elimination method to determine the unknown parameters required for approximating the solution of the differential equation. The advantage of this technique over other methods is that, it has less computational manipulations and complexities and also its availability for application on both linear and nonlinear second-order initial value problems is impressive. Other advantage of the algorithm is that high accurate approximate solutions are achieved by using a greater number of terms of the Legendre polynomial and once the operational matrix is obtained, it can be used to solve differential equations of higher order by introducing just a little manipulation on the operational matrix. Some existing sample problems from literature were solved and the results were compared to show the validity, simplicity and applicability of the proposed method. The results obtained validate the simplicity and applicability of the method and it also reveals that the method perform better than most existing methods.

A numerical study for non-linear multi-term fractional order differential equations

In this article, a numerical method is presented for solving non-linear multi-order fractional differential equations (FDEs) using orthonormal Boubaker polynomials as basis functions. The operational matrix is constructed for the product of basis functions. The fractional derivative is considered in the Caputo sense. Using this technique, non-linear multi-order FDEs converted into a system of non-linear algebraic equations, which can then be solved numerically. The proposed method is implemented using MATLAB for numerical analysis. Using numerical results, validity, and accuracy of the proposed method is given and compared the results with existing literature

Stagnant rings and uniform film analysis of Phan-Thien Tanner fluid film flow on a vertically upward moving tube

Analyzing Phan-Thien Tanner fluid film behavior on a vertically upward moving tube can help predictive models in engineering, notably in coating and lubrication operations. This paper provides a theoretical analysis of the dynamics of stagnant rings and uniform Phan-Thien Tanner fluid film adhered to a vertically upward moving tube. Formulated ordinary differential equations are solved to get exact analytic expressions for velocity, flow rate, average velocity, shear stress components, and stagnant rings. Highly nonlinear algebraic equations are solved using Newton's method in MAPLE to find the linear and exponential Phan-Thien Tanner film thicknesses. The uniform film thickness widens with increasing constant tube velocity, while it decreases with increasing Deborah number, Stokes number, and elongation parameter. The analysis delineates that stagnant rings shrink around the tube as the Stokes number, Deborah number, and elongation parameter increase. The minimal Stokes condition for the existence of a realistic stagnant ring is also determined. It has been established that whenever the Stokes number is less than the minimal Stokes condition, stagnant rings do not form. A comparison between the exponential Phan-Thien Tanner, linear Phan-Thien Tanner, upper convected Maxwell, and Newtonian fluids is also provided for stagnant rings and fluid film thickness. Following the application of an efficient approximation on tube geometry, plate geometry is approximated, and the outcomes are consistent with existing literature. The results of this research are significant for a wide range of biofluid applications, including agrochemical uses, paint and surface coating flow behavior, thin films on the cornea and lungs, and chemical and nuclear reactor design.

Modeling the Efficiency of Biogas Plants by Using an Interval Data Analysis Method

This article considers the task of developing mathematical models and their computer implementation that would establish the dependence of pH (acidity of the environment) on the volume and structure of raw materials for daily loading, as well as on the operating parameters of temperature and humidity based on the interval analysis of experimental data obtained during BGP research of a given type. In the process of research, based on the developed interval models, it was established that this indicator depends on the volume and structure of raw materials, as well as on the temperature and humidity of the substrate in the bioreactor. To build this mathematical model, it is proposed to use the method of interval data analysis and the method of identification of model parameters based on multidimensional optimization. The results of experimental studies for a specific type of biogas plant are given, and interval models with guaranteed prognostic properties that characterize the pH of the environment depending on the specific type of bio-raw material of solid and liquid fractions, temperature, and humidity are obtained. Based on the use of different types of raw materials, the developed models, based on experimental data, describe different configurations of the structure and volumes of raw materials for daily loading. The obtained mathematical models are an algebraic nonlinear equation that can be applied to control the level of pH of the environment in the bioreactor by determining the optimal volumes of raw materials of each type during the loading period depending on the temperature and humidity of the substrate in the bioreactor.

Haar Wavelet Approach for the Mathematical Model on Hepatitis B Virus

Abstract The Haar wavelet collocation method, a wavelet technique, is discussed in this article to examine the mathematical model of Hepatitis B virus infection. We took into account the HB virus, cytotoxic T lymphocytes (CTL) immune response, birth rate, death rate, and infected and uninfected hepatocytes to identify the dynamics of the hepatitis B virus infection. An ordinary differential equation (ODE) system that is nonlinear makes up this model. Using this method, the Hepatitis B Virus model can be solved by expressing each dependent variable as a Haar wavelet and then converting the system of ordinary differential equations into a system of nonlinear algebraic equations. The unknown coefficient values are thought to be extracted using the collocation procedure and the Newton–Raphson method. Tables and graphs are used to illustrate the characteristics of the Hepatitis B virus. The obtained results show that the current approach outperforms other approaches found in the literature in terms of accuracy. Mathematica software is utilized to obtain numerical results and nature.

Determining the influence of wheelset arrangement in the model 18-100 bogies on the level of steering efforts in the wheel-rail flange contacts

The object of this study is the process of directing the bogies of model 18-100 freight cars along a rail track, in particular, along curved sections of the track. The task to be solved was to determine the influence of wheelset arrangement on the level of steering forces in the flange contacts. A calculation diagram and a mathematical model of fitting the bogie into the curved section of the track have been built. The bogie loading scheme by external forces, including lateral rocking forces acting on the car in a curve, has been refined. In this case, the method of pseudo-statics of the mechanical system, which is a system of nonlinear algebraic equations, was applied. The calculation module Given-Find in the Mathcad software package (USA) was used to solve the mathematical model. It was established that the misalignment of wheelsets in the frames of model 18-100 bogies was of an accumulative nature. At the maximum operating angles of misalignment of wheelsets, the lateral steering forces in the flange contacts increase by 40–60 % compared to the rated setting. These angles can be up to 0.015 rad (0.85 degrees). The field of practical application of the results is railroad transportation, in particular, the system of maintenance and repair of freight cars on bogies of the 18-100 model. At the same time, the condition for the practical application of the research results is the expediency of introducing into the maintenance system the technological operation of controlling the deviation of wheelset arrangement in the bogie relative to the rated one. The current study will contribute to the construction of a measuring system for monitoring the deviation of wheelset arrangement in the bogie relative to the rated one. This proves the expediency of introducing into the trolley maintenance system the technological operation of controlling the deviation of wheelsets and designing a device for monitoring this parameter

An efficient approach for solving a class of fractional anomalous diffusion equation with convergence

This article presents a study on Fractional Anomalous Diffusion (FAD) and proposes a novel numerical algorithm for solving Caputo’s type fractional sub-diffusion equations to convert the fractional model into a set of nonlinear algebraic equations. These equations are efficiently solved using the Levenberg-Marquardt algorithm. The study provides the error analysis to validate the proposed method. The effectiveness and accuracy of the method are demonstrated through several test problems, and its performance and reliability are compared with other existing methods in the literature. The results indicate that the proposed method is a reliable and efficient technique for solving fractional sub-diffusion equations, with better accuracy and computational efficiency than other existing methods. The study’s findings could provide a valuable tool for solving FAD in various applications, including physics, chemistry, biology, and engineering.

Steady-state flow under surface recharge through unconfined aquifer with depth-decaying hydraulic conductivity: Analytical results

This paper presents some analytical results related to groundwater flow through an unconfined horizontal aquifer induced by areal recharge on the top and a decaying hydraulic conductivity. Under steady-state conditions, and in the context of Boussinesq equation, exact analytical solutions describing drainage into one or two penetrating channels at the aquifer boundaries have been derived in closed forms. The analytical expressions have been used to investigate the existence of water divide inside the unconfined aquifer. Based on the derived analytical solutions, a closed-form analytical formula is developed for calculating groundwater transit time within Dupuit-type flow system. The analytical transit time expression is simply expressed in terms of common mathematical functions. Analytical results are verified against a numerical model and show excellent agreements. The analytical transit time is validated against a numerical one obtained from a two-dimensional model representing a rectangular unconfined aquifer. Based on the explicit closed-form solution, a fully analytical inverse problem to estimate the decay parameter is proposed using a least-square objective function. The solution of the inverse problem simply corresponds to the solution of a nonlinear algebraic equation obtained from the definition of the objective function and the form of the analytical solution.

Analytical method for systems of nonlinear singular boundary value problems

The standard Lane–Emden equations model several physical phenomena such as isotropic continuous media, thermal behaviour of a spherical cloud of gas, and isothermal gas spheres. Systems of Lane–Emden-type equations appear in the modelling of the concentration of carbon substrate and oxygen, catalytic diffusion reactions, the steady state concentration of carbon dioxide, dusty fluid models, and pattern formation. In solving singular boundary value problems, one faces challenges resulting from the divergence of the associated variable coefficients at the singular points. This research article explores the power series approach to analytically approximate solutions to a class of systems of strongly nonlinear singular boundary value problems of Lane–Emden-type. The nonlinear terms in the proposed problems are transformed into power series using the generalised Cauchy product before establishing explicit recursion formulae for the expansion coefficients of the system of series solutions. The initial conditions required in the proposed boundary value problems are assumed and determined from a set of nonlinear algebraic equations resulting from the given right boundary conditions. Three special systems of nonlinear singular boundary value problems of Lane–Emden type are presented to demonstrate the proposed method’s reliability, effectiveness, and accuracy. The obtained approximate solutions are compared with the exact solutions (where they are available), or with other existing results (where the exact solution is not readily available). The series solution of the first example has a slow convergent rate, while the results of the other two examples are in excellent agreement with the exact solutions and other published results.

Physical invariants-preserving compact difference schemes for the coupled nonlinear Schrödinger-KdV equations

In this paper, we develop efficient compact difference schemes for the coupled nonlinear Schrödinger-KdV (CNLS-KdV) equations to conserve all the physical invariants, namely, the energy of oscillations, the number of plasmon, the number of particle and the momentum. By combining the exponential scalar auxiliary variable (E-SAV) approach, we reconstruct the original CNLS-KdV equations and adopt the compact difference method and Crank-Nicolson method to develop energy stable scheme. The E-SAV compact difference scheme preserves the total energy and the number of particle. We further introduce two Lagrange multipliers for the E-SAV reformulation system to develop compact difference scheme, which preserves exactly the number of plasmon and the momentum. At each time step for the second scheme, we only need to solve linear systems with constant coefficients and nonlinear quadratic algebraic equations which can be efficiently solved by Newton's iteration. Numerical experiments are given to show the effectiveness, accuracy and performance of the proposed schemes.

Non-local second-order boundary value problems with exponential nonlinearity

In this article, the Numerical treatment (Impilict-Explicite) was applied, for solving problems- related to the nonlocal second-order boundary value problem with exponential nonlinearity. The finite difference method of treatment of the BVP is used as implicit technique within the smallest available subdomain, . A reduced system of nonlinear algebraic equations corresponding to standard BVP is generated within the smallest domain, . A reduction of the computational work is introduced through reducing the size of the nonlinear algebraic system. A marching explicit approach is reemployed to generate the solution in , the complement of the interval with respect to the overall domain. Simple numerical experiments confirming the applicability of the treatment are introduced.

Variable order fractional diabetes models: numerical treatment

ABSTRACT In this paper, we present two numerical approximation solutions for the updated nonlinear variable order fractional differential equations of diabetes disease. The variable order fractional differential operator is based on the Caputo fractional derivative concept. The first used method is developed using the nonstandard finite difference process, such that a new quick approximation is applied to tackle the nonlinear terms. By applying some reasonable assumptions to the given data, we can prove that this scheme maintains the positivity and boundedness of the discretized solutions. The second method is based on approximation with shifted Jacobi polynomials. The properties of Jacobi polynomials, together with the shifted Jacobi-Gauss-Lobatto nodes were utilized to reduce the proposed system into a set of algebraic nonlinear equations. Furthermore, some numerical experiments are conducted to compare the performance of the introduced schemes. Our numerical simulations show that the nonstandard finite difference approach requires less computational time compared to the Jacobi-Gauss-Lobatto spectral collocation. While the last one is more accurate when used for solving the variable order fractional diabetes model.

Application of generalized Haar wavelet technique on simultaneous delay differential equations

In this paper, the Haar wavelet technique is applied to a system of simultaneous linear and nonlinear proportional delay differential equations. By implementing this method, the delay differential equations are transformed into a series of linear and nonlinear algebraic equations, respectively. These equations are then solved using numerical code constructed in MATLAB. The obtained approximate solutions are compared with exact solutions, demonstrating the accuracy and efficacy of the Haar wavelet method. Using the fixed point theorem, the existence and uniqueness of the solutions is established. The reliability and efficiency of the proposed technique are illustrated through two numerical examples. Additionally, a comprehensive error analysis is presented and the convergence result is discussed, providing a thorough examination of the robustness of the method. This study not only shows that the Haar wavelet technique is a powerful tool for solving delay differential equations but also sets the groundwork for future research in this area.

A linearly implicit finite element full-discretization scheme for SPDEs with nonglobally Lipschitz coefficients

Abstract The present article deals with strong approximations of additive noise driven stochastic partial differential equations (SPDEs) with nonglobally Lipschitz nonlinearity in a bounded domain $ \mathcal{D} \in{\mathbb{R}}^{d}$, $ d \leq 3$. As the first contribution, we establish the well-posedness and regularity of the considered SPDEs in space dimension $d \le 3$, under more relaxed assumptions on the stochastic convolution. This improves relevant results in the literature and covers both the space-time white noise ($d=1$) and the trace-class noises ($\text{Tr} (Q) < \infty $) in multiple dimensions $d=2,3$. Such an improvement is achieved based on a key perturbation estimate for a perturbed PDE, with the aid of which we prove the convergence and uniform regularity of a spectral approximation of the SPDEs and thus get the improved regularity results. The second contribution of the paper is to propose and analyze a spatio-temporal discretization of the SPDEs, by incorporating a standard finite element method in space and a linearly implicit nonlinearity-tamed Euler method for the temporal discretization. The proposed time-stepping scheme is linearly implicit and does not suffer from solving nonlinear algebra equations as the backward Euler scheme does. Based on the improved regularity results, we recover the expected strong convergence rates of the fully discrete scheme and reveal how the convergence rates rely on the regularity of the noise process. In particular, a classical convergence rate of order $O(h^{2} +\tau )$ can be obtained even in high dimension $d=3$, as the driven noise is of trace class and satisfies certain regularity assumptions. The optimal error estimates turn out to be challenging and face some essential difficulties when the tamed time-stepping scheme meets the finite element spatial discretization, particularly in the context of low regularity and multiple dimensions $d \le 3$. Some highly nontrivial arguments are introduced to overcome the difficulties. Finally, numerical examples corroborate the claimed strong orders of convergence.

Advancements in numerical modeling of ideal multicomponent isotope separation Cascades: A comparative study

The complexity of operational conditions in stable isotope separation cascades necessitates the utilization of numerical methods for cascade design and performance analysis. This study focuses on solving a system of nonlinear algebraic equations derived from the mass and energy balance equations, which is a fundamental challenge in analyzing and designing cascades and developing numerical methods. Two existing numerical algorithms were selected from the literature and modified for quasi-ideal cascade systems. The performance of these algorithms was examined and compared, with a particular emphasis on influential parameters such as the number of stages, the number of isotopes, and the separation factor. The first algorithm solves the concentration and flow equations separately, while the second algorithm solves all equations simultaneously. The results indicate that the first algorithm exhibits significant superiority over the second algorithm. In most cases, the second algorithm either fails to converge or exhibits a slower convergence speed. Additionally, the results demonstrate that increasing the number of stages impacts the convergence time, while an increase in the separation factor leads to greater convergence complexity.

Efficient invariant-preserving scheme for the [formula omitted]-coupled nonlinear Schrödinger equations

The N-coupled nonlinear Schrödinger equations are reformulated into an expanded form by using a scalar auxiliary variable method, and then a scheme preserving exactly original mass and energy conservation laws is proposed based on discretizing the expanded form. The scheme is efficient as it consists of decoupled linear systems with constant coefficients, along with a nonlinear algebraic equation that can be solved with negligible computational cost. Some numerical experiments are carried out to demonstrate the behavior of wave solutions, the accuracy of solution and the preservation of physical invariants.