Nonlinear Algebraic Equations Research Articles

A novel hybrid method with convergence analysis for approximation of HTLV-I dynamics model

This paper presents a novel numerical approach for approximating the solution of the model describing the infection of CD4+T\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$CD4^+T$$\\end{document}-cells by the human T-cell lymphotropic virus I (HTLV-I).The proposed method utilizes the operational matrix along with spectral method to convert the fractional model into a system of nonlinear algebraic equations. The Levenberg-Marquardt algorithm efficiently solves these equations. The study includes theoretical convergence analysis and error bounds to establish the validity of the proposed method. Through several test problems, we demonstrate the effectiveness and accuracy of the approach. We compare its performance and reliability to other existing methods in the literature. The results indicate that the proposed method is a reliable and efficient approach for solving the model.

Galerkin’s Method and Multivariate Newton’s Method for the Nonlinear Deformation of the Magneto-Electro-Elastic Bi-Layered Laminates

In this paper, mathematical modelling for the large deformation of a magneto-electro-elastic rectangular bi-layered laminate with general boundary conditions is presented. Constitutive equations involving the magneto-electro-elastic (MEE) material properties are introduced, Maxwell equations accounts for the electric and magnetic effects are also utilized. First-order shear deformation theory (FSDT) considering the von Karman nonlinear strain is adopted, and the plain strain/stress assumption applicable for thin plate analysis is used. A rather compact set of governing equations related to kinematical variables, electric/magnetic potentials and the Airy stress function is obtained as a consequence of the preliminary condensation for the electro-magnetic state to the plate kinematics. Semi-analytic solution for a bi-layered BaTiO3-CoFe2O4 laminate with specified boundary conditions subjected to various external applied loads is performed. By employing the Bubnov-Galerkin method, the set of nonlinear partial differential equations is transformed to a set of third-order nonlinear algebraic equations for the static deformation due to applied load. Numerical results are carried out by using the multivariate Newton's method with respect to various volume fractions indicating the volume ratio between piezoelectric BaTiO3 layer and piezomagnetic CoFe2O4. From the result, the nonlinearity of the von Karman strain appears to enhance system rigidity as smaller deformations will be detected when external load is applied. Also, some other interesting results are obtained which could be useful to future analysis and design of multiphase composite plates.

The Extended Weierstrass Transformation Method for the Biswas–Arshed Equation with Beta Time Derivative

In this article, exact solutions of the Biswas–Arshed equation are obtained using the extended Weierstrass transformation method (EWTM). This method is widely used in solid-state physics, electrodynamics, and mathematical physics, and it yields exact solution functions involving trigonometric, rational trigonometric, Weierstrass elliptic, wave, and rational functions. The process involves expanding the solution functions of an elliptic differential equation into finite series by transforming them into Weierstrass functions. Furthermore, it generates parametric solutions for nonlinear algebraic equation systems, which are particularly useful in mathematical physics. These solutions are derived using the Mathematica package program. To analyze the behavior of these determined solution functions, the article employs separate two- and three-dimensional graphs showing the real and imaginary components, along with contour and density graphs. These visuals aid in comprehending the physical characteristics exhibited by these solution functions.

Flow dynamics over cylinder in two‐dimensional benchmark problem: A finite element approximation

This research investigates 2D benchmark flow around a circular cylinder, utilizing the incompressible Navier–Stokes equations alongside the continuity and energy equations. The numerical solution is achieved through finite element discretization for space variable, combined with a second‐order Crank–Nicolson scheme for time integration. The computational results are derived using the FEATFLOW finite element based library package. Our study focuses on the dimensionless form of the flow equations and examines three key dimensionless parameters: drag, lift, and pressure drop. Upon applying finite element method (FEM) discretization, the system is converted into a set of linear or nonlinear ordinary differential equations or algebraic equations for steady‐state scenarios. We then apply the Newton–Raphson method as the outer nonlinear solver, and the multigrid method for efficiently resolving the linear subproblems. To ensure numerical accuracy, we evaluated the and errors, confirming that the experimental order of convergence matches the theoretical predictions. Flow profiles were both graphically represented and tabulated, offering a detailed understanding of the simulation results.

A Robust numerical technique based on the chromatic polynomials for the European options regulated by the time-fractional Black–Scholes equation

AbstractRisk mitigation and control are critical for investors in the finance sector. Purchasing significant instruments that eliminate the risk of price fluctuation helps investors manage these risks. In theory and practice, option pricing is a substantial issue among many financial derivatives. In this scenario, most investors adopt the Black–Scholes model to describe the behavior of the underlying asset in option pricing. The exceptional memory effect prevalent in fractional derivatives makes it easy to understand and explain the approximation of financial options in terms of their inherited characteristics prompted by the given reason. Finding numerical solutions that are both successful and suitably precise is crucial when working with financial fractional differential equations. Hence, this paper proposes an innovative method, designated the Chromatic polynomial collocation method (CPM), for the theoretical study of the Time fractional Black–Scholes equation (TFBSE) that regulates European call options. The newly developed numerical algorithm CPM is on a functional basis of the Chromatic polynomials of Complete graphs (Kn) and operational matrices of the basis polynomials. The CPM transforms the TFBSE into a framework of nonlinear algebraic equations with the help of operational matrices and equispaced collocation points. The fractional orders in the PDE are concerned in the Caputo sense. The CPM findings further corroborate the results of the most recent numerical schemes to show the effectiveness of the suggested numerical algorithm.

A Meshless Radial Point Interpolation Method for Solving Fractional Navier–Stokes Equations

This paper aims to develop a meshless radial point interpolation (RPI) method for obtaining the numerical solution of fractional Navier–Stokes equations. The proposed RPI method discretizes differential equations into highly nonlinear algebraic equations, which are subsequently solved using a fixed-point method. Furthermore, a comprehensive analysis regarding the effects of spatial and temporal discretization, polynomial order, and fractional order is conducted. These factors’ impacts on the accuracy and efficiency of the solutions are discussed in detail. It can be shown that the meshless RPI method works quite well for solving some benchmark problems accurately.

Buckling disappearance via merging/divergence in a nonlinear three-d.o.f. system with linear constitutive law

The phenomenon of buckling disappearance, occurring in a parameter-dependent family of systems admitting a nontrivial fundamental path, is studied. Two different forms of disappearance are detected, namely: (i) the divergence, in which the critical load continuously tends to infinity, and (ii) the merging, in which two critical loads approach each other, coalesce, and then disappear at a finite value of the critical load. It is shown that the two phenomena can be exhibited by the same mechanical system, when a suitable elasto-geometric parameter is varied. More importantly, it is proved that merging continuously changes into divergence when a second parameter is changed. A paradigmatic system is chosen to illustrate the two forms of buckling, i.e., a three degree-of-freedom spherical pendulum, elastically constrained at the ground, loaded by a transverse force and/or a conservative couple, made of two longitudinal potential forces. The springs are taken elastically linear, to stress the fact that divergence not necessarily calls for introducing a nonlinear constitutive law, as also mentioned in literature. Only a linear bifurcation analysis is carried out here, aimed to find the bifurcation points along the nonlinear fundamental path. However, due to the presence of non-negligible prestrains, such a bifurcation problem is governed by nonlinear algebraic equations, whose number of roots cannot be predicted in advance.

Discretization methods and their extrapolations for two-dimensional nonlinear Volterra-Urysohn integral equations

In this paper, a class of nonlinear two-dimensional (2D) integral equations of Volterra type, i.e. Volterra-Urysohn integral equations, is studied. Following the ideas of [24], and assuming that the kernels of the integral equation are Lipschitz functions with respect to the dependent variable, the existence and uniqueness of a solution to the integral equation is shown by a technique based on the Picard iterative method. Then, the Euler and trapezoidal discretization methods are used to reduce the solution of the integral equation to the solution of a system of nonlinear algebraic equations. It is proved that the solution of the Euler method has first order convergence to the exact solution of the integral equation while the solution of the trapezoidal method has quadratic convergence. To prove the convergence of the trapezoidal method, a new Gronwall inequality is developed. Some numerical examples are given which confirm our theoretical results.

Numerical solution of fuzzy non-linear equations under generalized trapezoidal fuzziness by Newton–Raphson method

In the real world, we are faced with many problems that cannot be easily described by crisp values, as there is always some uncertainty involved in real-life situations and in those situations, we end up receiving fuzzy values rather than crisp ones in those situations. Some preliminary concepts of fuzzy sets and operations on trapezoidal fuzzy numbers incorporated with the membership function have been discussed here. This paper aims to demonstrate the advantages of the newly developed fuzzified newton-raphson method with the help of trapezoidal fuzzy numbers for solving algebraic nonlinear equations arising in fuzzy environments. The algorithm of newton-raphson techniques for solving the same has been described. This paper presents several examples of how the new fuzzified newton-raphson method works.

Determination of the plastic J integral of ductile material using the XFEM with only Heaviside function and variable-node elements

In this paper, the extended finite element method (XFEM) with only Heaviside function is proposed for the elastic–plastic fracture mechanics (EPFM) modeling. The proposed method removes tip enrichment functions depending on the polar coordinates at the crack front, and a step function only determined by the level set function is utilized to track the crack front. To alleviate the volumetric locking phenomena caused by the plastic incompressibility, the B-bar method is incorporated into the three dimensional (3D) XFEM program. Therefore, the fully integration scheme can be chosen to ensure the accuracy when addressing large plastic deformation in EPFM analysis. Additionally, the material tangent stiffness matrix of Ramberg-Osgood constitutive is given, and the local refinement technique using variable-node elements is adopted to reduce the number of elements and nodes for efficient analysis. A Newton-Raphson iterative algorithm is developed to solve the nonlinear algebraic equations caused by material nonlinearity. Several numerical examples including the determination of crack opening displacement, and the fully plastic J integral in the ductile materials are presented to test the performance of the proposed method. Comparisons with the results from the existing methodologies show that the new enrichment scheme can save computational cost and obtain sufficient accuracy even in the case of 3D curved crack.

Differential-algebraic equation-based model predictive control for trajectory tracking of spacecraft-mounted continuum manipulators

Spacecraft-mounted continuum manipulators (SMCMs) exhibit great potential for performing dexterous operations in unstructured environments due to their inherent compliance and dexterity. However, the dynamic model of the rigid-flexible coupling SMCM is highly nonlinear and typically formulated as a set of implicit differential-algebraic equations (DAEs), posing significant challenges for precise trajectory tracking control. This paper proposes a novel model predictive control (MPC) framework specifically designed for generic DAEs to achieve precise trajectory tracking of the SMCM under uncertain disturbances and input limitations. The DAE model of the SMCM is discretized into a set of nonlinear algebraic equations. By performing implicit differentiation of these equations with respect to the system state, the state transition matrix (STM) for the DAE model is derived. The optimal control action for the SMCM can be further determined based on the derived STM. Additionally, nonlinear complementary functions are introduced to address the issue of input limitations, allowing the problem of determining the optimal control sequence to be equivalently transformed into a set of nonlinear algebraic equations for solving. Numerical simulations demonstrate that the proposed approach can achieve precise trajectory tracking of the SMCM while strictly adhering to input limitations.

A variational framework for Cahn–Hilliard-type diffusion coupled with Allen–Cahn-type multi-phase transformations in elastic and dissipative solids

This article presents a variational framework for coupled chemo-mechanical solids undergoing irreversible micro-structural changes at infinitesimal strains. The coupled problem is characterised by phenomena such as phase transitions, micro-structure coarsening and swelling. It is an extension of our previous work on variational inelasticity for a conserved chemo-mechanical setting to a unified conserved and non-conserved setting which include multi-phase transformations. The variational framework, again governed by continuous-time, discrete-time and discrete-space–time incremental variational principles, is outlined for coupled diffusion-phase transformation phenomena in elastic and dissipative solids. For the sake of simplicity, focus is restricted to isothermal conditions. It is shown that the governing macro- and micro-balance equations of the coupled problem appear as Euler equations of these minimisation and saddle point principles. In contrast to our previous work, extended variational principles (with the gradient of the chemical potential and phase fractions) are constructed that account for diffusion-phase transformation coupling. This is achieved by Legendre transformations. Note that the local–global solution strategy is still preserved and the resulting system of symmetric non-linear algebraic equations are solved by Newton–Raphson-type iterative methods. The applicability of the proposed framework is demonstrated by numerical simulations that qualitatively characterise lower bainitic micro-structure.

Energy harvesting from transverse galloping using a two-degree-of-freedom flow energy harvester

Abstract A two degree-of-freedom (2-DOF) galloping piezoelectric energy harvesting system where both oscillators are excited by the incoming flow is studied in this paper. A coupled lumped-parameter model is developed to simulate the electro-fluid-structural coupling behaviour assuming a quasi-steady aerodynamic flow field. The differential equations governing the dynamics of the lumped system are converted into nonlinear algebraic equations employing the harmonic balance method (HBM) and then solved using Newton’s method. The approximate analytical solutions are compared to the solutions obtained by numerical integration, and the results agree, both qualitatively and quantitatively. The solutions were subsequently used to investigate the effects of the bluff body dimension, mechanical, electrical, aerodynamic, and electromechanical parameters on the performance of the harvester and to illustrate how to optimise some of these parameters to reduce the cut-in wind speed and maximise the power output of the energy harvester. It will be shown that the energy harvesting performance could be significantly improved by introducing piezoelectric transducers on both oscillators and including both masses in the incoming flow. A comparative study between the proposed 2-DOF flow energy harvesting system, a single-degree-of-freedom (SDOF) galloping piezoelectric energy harvester (GPEH), and other configurations of 2-DOF GPEH shows that the present 2-DOF GPEH generates the largest power peak.

A novel unified variational finite difference (VFD) solution method for optimal control problems

This article is aimed to propose a simple yet efficient unified numerical strategy for solving both linear and non-linear optimal control problems. To do so, the general form of quadratic performance index function and nonlinear state equations are considered first. Then, the idea of variational differential quadrature method is used to convert the integral/differential equations to the equivalent algebraic form. Since time is the only independent variable is this research, the finite difference method with an equally-spaced discretization scheme would be a more appropriate technique rather than the differential quadrature approach. So, the implemented numerical solution is called now as the variational finite difference method. The method of Lagrange multipliers is then utilized for minimization purpose and, as a result, the final set of nonlinear algebraic equations are obtained. Finally, the quadratic and triadic forms of non-linearity are considered and an explicit formulation is represented for the residual and Jacobian of the Newton-based iterative solution procedure. To demonstrate the accuracy and efficiency of the proposed approach in the quadratic optimal control area, several benchmark problems involving linear time-invariant, linear time-variant and nonlinear examples are successfully solved and the results are confirmed with those existed in the literature.

Assessment of data‐driven modeling approaches for chromatographic separation processes

AbstractChromatographic separation processes are described by nonlinear partial differential and algebraic equations, which may result in high computational cost, hindering further online applications. To decrease the computational burden, different data‐driven modeling approaches can be implemented. In this work, we investigate different strategies of data‐driven modeling for chromatographic processes, using artificial neural networks to predict pseudo‐dynamic elution profiles, without the use of explicit temporal information. We assess the performance of the surrogates trained on different dataset sizes, achieving good predictions with a minimum of 3400 data points. Different activation functions are used and evaluated against the original high‐fidelity model, using accuracy, interpolation, and simulation time as performance metrics. Based on these metrics, the best performing data‐driven models are implemented in a process optimization framework. The results indicate that data‐driven models can capture the nonlinear profile of the process and that can be considered as reliable surrogates used to aid process development.

Improved numerical schemes to solve general fractional diabetes models

In this article, we propose a new class of nonlinear fractional differential equations of diabetes disease based on the concept of Caputo fractional derivative. Two numerical techniques are introduced to analyse the solution of the general fractional diabetes model, that describes glucose homeostasis. The first method, which constructed using the nonstandard finite difference technique involves an asymptotically stable difference scheme. This method maintains important properties of the solutions of the studied system, such as the positivity and boundedness. The second method is the Jacobi-Gauss–Lobatto spectral collocation approach, known for its exponential accuracy. By employing this collocation method, the problem is transformed into a set of algebraic nonlinear equations, simplifying the overall task. Numerical simulations were conducted to compare the performance of these two techniques with other standard methods and the analytic solution in specific cases. Our findings show that the Jacobi-Gauss–Lobatto spectral collocation technique provides higher accuracy in solving the fractional diabetes model system, while the nonstandard finite difference approach requires lower computational duration.

A compact combination of second-kind Chebyshev polynomials for Robin boundary value problems and Bratu-type equations

AbstractThis paper presents a new algorithm for resolving linear and non-linear second-order Robin boundary value problems (BVPS) and the Bratu-type equations in one and two dimensions using spectral approaches. Basis functions according to second-kind shifted and modified shifted Chebyshev polynomials that comply with the Robin conditions are created. It has produced operational matrices for its derivatives. The provided solutions are the result of applying the collocation and tau approaches. These methods convert the problem dictated by its boundary conditions into a system of linear or non-linear algebraic equations that may be solved using any suitable numerical solver. Convergence analysis has been provided and it accords with the numerical results. Six numerical problems are provided to investigate and demonstrate the practical utility of the suggested method. The current results show that our method outperforms the previous methods in terms of accuracy which are presented in tables and figures.

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.