System Of Nonlinear Differential Equations Research Articles

Impact of chemical reaction on the Cattaneo–Christov heat flux model for viscoelastic flow over an exponentially stretching sheet

In this article, the numerical solutions for the heat transfer flow of an upper-convected Maxwell fluid across an exponentially stretched sheet with a chemical reaction on the Cattaneo–Christov heat flux model have been investigated. Using similarity transformation, the controlling system of nonlinear partial differential equations was transformed into a system of ordinary differential equations. The resulting converted equations were solved numerically by a successive linearization method with the help of MATLAB software. A graphic representation was created to analyze the physical insights of the relevant flow characteristics. The findings were presented in the form of velocity, temperature, and concentration profiles. As the relaxation time parameter varied, the local Nusselt number increased. The thermal relaxation time was shown to have an inverse relationship with fluid temperature. Furthermore, the concentration boundary layer becomes thinner as the levels of the reaction rate parameter increase. The results of this model can be applicable in biological fluids and industrial situations. Excellent agreement exists between the analysis's findings and those of the previous studies.

A study using a semi-analytical approach to the reaction kinetics that describes the behavior of heterogeneous reaction–diffusion process with a glyoxyl-agarose immobilized enzyme (penicillin G acylase) in a batch reactor

This research article examines the reaction–diffusion process in a batch reactor by means of immobilized penicillin G acylase (PGA) with a glyoxyl-agarose matrix support. This study employs the Akbari-Ganji’s Method (AGM) when the system is under the internal (intra-particle) diffusional restrictions, and the external (film) diffusional restrictions is negligible. The semi-analytical approximations by the AGM technique are employed for computing the dimensionless steady-state solutions for a system of nonlinear differential equations to examine the impact of internal diffusional restrictions (IDR) as a function of catalyst enzyme loading and particle size. The model entails highly nonlinear components that adhere to the standard Michaelis-Menten kinetics. In addition to the proposed model, the dynamics for the mean integrated effectiveness factor (MIEF) of the dimensionless substrate concentration phenylglycine methyl ester (PGME) is presented. In light of this, there is a satisfactory level of agreement on the comparison of the semi-analytical result with the simulated numerical results across the whole concentration range where the dimensionless initial substrate and product concentrations are S1b=30,S2b=15, and P1b=60,P2b=100. The behavior of concentration profiles under the steady/unsteady state conditions for the effect of IDR on the overall reaction rates by immobilized biocatalyst PGA were examined. A sensitivity analysis is carried out to identify the effective parameter which can influence the diffusional restrictions on MIEF. An error analysis is performed on the dimensionless concentration equations to assess the amount of accuracy and the convergence of a solution.

A numerical inquiry of the MHD radiative thin film nanofluid flow and heat transfer augmentation for Ni and Ta nanoparticles of different shapes dispersed in C 12 H 26–C 15 H 32 over an unsteady, curved stretching surface

Thin films are a versatile technology used in various industries, including light-emitting diodes, magnetic recording media, electronic device manufacturing, and optical coating. They are crucial for studying quantum phenomena, superlattices, and multiferroic materials. Their small size allows for better thermal transport studies due to the combination of thin film flow and nanoparticles. Therefore, this study aims to investigate thin film flow and heat transfer with the effects of thermal radiations, magnetic field, and velocity slip in association with tantalum ( Ta ) and nickel ( Ni ) nanoparticles that emerged in kerosene oil ( C 12 H 26 − C 15 H 32 ) toward a curved stretching surface. A review of the pertinent literature indicates that no research has previously been done on the present problem of an unsteady thin film flow across a curved stretching surface. The problem under consideration is modeled by choosing the curvilinear coordinate system. The controlling nonlinear system of partial differential equations is transformed into the system of ordinary differential equations by using an appropriate similarity transformation and then solved numerically using the BVP 4 C technique in MATLAB . The impact of significant physical parameters on velocity, temperature, skin friction, and Nusselt numbers are depicted through graphs and radar charts. The investigation highlights the significant influence of magnetohydrodynamics ( MHD ) , thermal radiation, volumetric fraction, and velocity slip on the enhancement of nanofluid temperature. It is also noted that cylindrical-shaped Ta nanoparticles exhibited the maximum velocity, while blade-shaped Ni nanoparticles showed the lowest. Moreover, the maximum and minimum temperatures were observed for blade-shaped Ni nanoparticles and cylindrical Ta nanoparticles, respectively. Additionally, the maximum and minimum values of skin friction are found for cylindrical Ni nanoparticles and blade-shaped Ta nanoparticles, respectively. Similarly, the maximum and minimum Nusselt numbers are recorded for cylindrical Ta nanoparticles and blade-shaped Ni nanoparticles, respectively.

Solving an electrically conducting nanofluid over an impermeable stretching cylinder problem with a spectral reproducing kernel method

In this paper, a nonlinear mechanical system of ordinary differential equations (ODEs) with multi-point boundary conditions is considered by a novel type of reproducing kernel Hilbert space method (RKHSM). To begin, we define the unknown variables in terms of the reproducing kernel function. The roots of the Shifted Chebyshev polynomials (SCPs) are then utilized to collocate the resulting system. Finally, Newton’s iterative method is employed to find the unknown expansion coefficients. The solutions of this system of equations, which arise from the flow of an electrically conducting nanofluid over an impermeable stretching cylinder, are numerically analyzed, and convergence analysis is discussed to demonstrate the reliability of the presented method (PM). Tables and figures are provided to further discuss the solutions and assess the effectiveness of the method in comparison to other techniques in the literature.

Effect of slip velocity on Newtonian fluid flow induced by a stretching surface within a porous medium

This article aims to examine an unsteady 2-D laminar flow of magnetohydrodynamic fluid caused by an elastic surface immersed in a permeable medium under the influence of thermal radiation and extended heat flux. Thermal conductivity and viscosity both are supposed to vary with temperature. This flow model also includes velocity slip, heat source, and joule heating. The governing equations of the fluid, including momentum and energy equations, of the proposed problem are transfigured into a system of interconnected non-linear ordinary differential equations through similarity transformations. The resultant equations are solved efficiently by employing the shooting technique in combination with the fourth-order Runge-Kutta method. Numerical values and the effect of numerous governing factors on the flow field, temperature distribution, local skin friction coefficient, and Nusselt number are showcased via graphs and tables. The investigation reveals that velocity slip, heat source, and porosity parameters enhance the temperature field while diminishing the velocity field. Furthermore, the velocity slip parameter notably reduces both the coefficient of skin friction and the Nusselt number.

Mathematical modeling of unemployment dynamics with skills development and cyclical effects

In the present work, a new mathematical model is introduced to investigate the effects of skills development on the governance of cyclical unemployment. A system of non-linear differential equations is formed and analyzed. In the modeling approach, we categorize the considered workforce into three distinct classes: the fundamentally unemployed Us, employed individuals E and cyclically unemployed people Uc. We also make the assumption that some government programs can qualify individuals who have been out of work for a period to reintegrate and prevent them from becoming long-term unemployed. The qualitative analysis of differential equations is used to analyze the proposed model. It has been demonstrated that the suggested model has only one equilibrium that is globally asymptotically stable under particular conditions. The impact of key parameters is analyzed by examining their sensitivity. In order to verify the important theoretical results, numerical illustrations of the model have been performed. The obtained results demonstrate that the development of skills among the cyclically unemployed is effective in reducing unemployment.

Hybrid Nanofluid Flow over a Shrinking Darcy-Forchheimer Porous Medium with Shape Factor and Solar Radiation: A Stability Analysis

This research aimed to develop a numerical solution to analyze the effects of solar radiation and nanoparticle shape factors on the flow of a hybrid nanofluid past a shrinking Darcy-Forchheimer porous medium. The base fluid chosen for this study is water (H2O), and the hybrid nanofluid consists of nanoparticles of silver (Ag) and titanium dioxide (TiO2) in four different shapes: bricks, cylinders, platelets, and blades. To account for solar radiation, the energy model incorporated a radiative heat flux, while the momentum problem considers the influence of a magnetic field. The application of an appropriate similarity transformation method converts the partial differential equations (PDEs) model into a system of nonlinear ordinary differential equations (ODEs). The mathematical model is solved using the shooting technique method and the bvp4c solver. The obtained results, along with the effects of the nanoparticle shape factor, solar radiation parameter, shrinking parameter, Darcy-Forchheimer number, and nanofluid volume fraction, are visually presented through figures and tables. It is worth noting that, in our numerical results, we observed the presence of dual solutions when λ < 0. Our findings indicate that the thermal transmittance increases with an increase in the nanoparticle shape factor and solar radiative parameter. Additionally, we observed an escalation in the velocity distribution in relation to the shrinking parameter and nanofluid volume fraction. Before reaching the two solutions, a flow stability analysis revealed that the first branch appears to be the most stable. Overall, these findings provide valuable insights into the behaviour of hybrid nanofluid flow in the presence of solar radiation and porous media.

Long-Term Behavior of Positive Solutions of a Certain Nonlinear System of Difference Equations

Long-Term Behavior of Positive Solutions of a Certain Nonlinear System of Difference Equations

Similarity solution for the magnetogasdynamic shock wave in a self-gravitating and rotating ideal gas under the influence of radiation heat flux

The Lie invariance method is used to analyze the one-dimensional, unsteady flow of a cylindrical shock wave in a rotating, self-gravitating, radiating ideal gas under the influence of an axial or azimuthal magnetic field, with an emphasis on adiabatic conditions. The analysis assumes a stationary environment just ahead of the shock wave and considers variations in fluid velocity, magnetic field, and density within the perturbed medium just behind the shock front. In the governing equations, the impact of thermal radiation under an optically thin limit is integrated into the energy equation. Utilizing the Lie invariance method, the set of partial differential equations governing the flow in this medium is transformed into a system of nonlinear ordinary differential equations (ODEs) using similarity variables. Two distinct cases of similarity solutions are obtained by selecting different values for the arbitrary constants associated with the generators. Among these cases, one yields similarity solutions assuming a power-law shock path and the other an exponential-law shock path. For both cases, the resulting set of nonlinear ODEs are numerically solved using the 4th-order Runge–Kutta method in MATLAB software. The article thoroughly explores the influence of various parameters, including γ (adiabatic index of the gas), Ma−2 (Alfvén–Mach number), σ (ambient density exponent), l1 (rotational parameter), and G0 (gravitational parameter) on the flow properties. The findings are visually presented to offer a comprehensive insight into the effects of these parameters.

Implementation and Linearization of a Coupled Panel and Vortex Particle Method in State-Space Form

This paper describes the implementation and linearization of a coupled panel and vortex particle method in a state-variable form. More specifically, the coupled panel and vortex particle dynamics are formulated as a nonlinear system of ordinary differential equations in first-order form to be self-contained and inherently linearizable. Linearization of this coupled panel and vortex particle method is demonstrated via finite differencing and a novel analytical linearization technique to yield a linear time-invariant representation. The code is implemented in MATLAB® and validated against an open-source aerodynamic solver for a fixed wing in both stationary and unsteady conditions. The linearized models are verified against the nonlinear dynamics both in the time and frequency domains. Linearized models accurately represent the wake dynamics about the equilibrium condition for moderate amplitude inputs and for input frequencies covering the typical frequency range of flight dynamics and flight controls, i.e., 0.3–30 rad/s. Analytical linearization is shown to abate the cost of linearization over perturbation methods by O(n2), where n is the total number of states of the system. Linearized models of the coupled panel and vortex particle dynamics have applications in the flight dynamics of both fixed- and rotary-wing vehicles, where these dynamics can be used to augment the rigid-body dynamics. The coupled rigid-body and wake dynamics can be leveraged to assess the stability, response characteristics, and handling qualities of vehicles experiencing aerodynamic interactions between rotors, wings, and/or obstacles.

On the Role of Early Case Detection and Treatment Failure in Controlling Tuberculosis Transmission: A Mathematical Modeling Study

Tuberculosis (TB) remains a pressing global health concern, demanding urgent attention to mitigate its spread and impact. In this study, we present a rigorous mathematical model of TB transmission that incorporates early case detection and addresses the critical issue of treatment failure. Through the development of a system of nonlinear ordinary differential equations, we conduct comprehensive analyses to assess the dynamics of TB transmission and the efficacy of intervention strategies. Our findings underscore the urgent need for effective TB control measures. Mathematical analyses reveal that the model exhibits a TB-free equilibrium, which is globally asymptotically stable only if the control reproduction number falls below one. However, we identify a concerning phenomenon: the model demonstrates a forward bifurcation when the control reproduction number equals one, suggesting that the disease-free equilibrium loses its stability, while simultaneously, the stable unique endemic equilibrium begins to emerge. Moreover, sensitivity analysis highlights the complex interplay between case detection rates, treatment failure probabilities, and TB transmission dynamics. Contrary to expectations, increasing case detection rates and minimizing treatment failure probabilities may not consistently reduce the basic reproduction number or the size of the infected population. Instead, there exists a critical threshold for intervention effectiveness, beyond which TB transmission can be significantly curtailed. Biologically, this phenomenon may occur if there is no balance between case detection and treatment efforts. If treatment quality does not improve, then case detection will not have a significant impact, and in the worst case scenario, it can exacerbate the intervention’s negative effects. These findings underscore the urgency of implementing targeted intervention strategies to combat TB transmission effectively. Failure to meet the critical intervention threshold risks undermining TB elimination efforts and exacerbating the global TB burden. Through numerical simulations, we elucidate potential intervention scenarios necessary for achieving TB elimination goals in human populations. In conclusion, our study highlights the urgent imperative for coordinated action to control TB transmission effectively. By elucidating the dynamics of TB spread and intervention efficacy, we provide valuable insights to inform evidence-based policy decisions and accelerate progress towards TB elimination on a global scale.

Flow analysis of temperature-dependent variable viscosity Phan Thien Tanner fluid thin film over a horizontally moving heated plate

Comprehending the temperature-dependent variable viscosity Phan Thien Tanner fluid film flow over a horizontal heated plate with surface tension is essential to enhancing coating and lubrication process prediction models. This paper accords with the flow analysis of a temperature-dependent variable viscosity Phan Thien Tanner fluid film over a horizontal heated plate. The flow over a heated plate occurs as a result of the plate’s motion and the surface tension gradient. The Adomian decomposition method is utilized to solve a system of linear and nonlinear ordinary differential equations, resulting in series-form solutions. The series-form expressions for flow variables such as velocity, temperature, volume flow rate, and surface tension are derived. Moreover, the positions of stationary points are computed using MATHEMATICA. The analysis delineated that when the inverse capillary number, variable viscosity parameter, Deborah number, elongation parameter, and Brinkman number increase, the stationary points move closer to the heated plate. The temperature also rises with an increase in these parameters. The temperature rises with increasing viscous dissipation while it lowers with increasing thermal diffusion. When the Deborah number is high, the Phan Thien Tanner fluid behaves like a solid and the flow is only driven by the motion of the plate. A comparison between Newtonian and Phan Thien Tanner fluids is made for velocity, temperature, and stationary points as a special case.

Mathematical model for applying electromagnetic (EM) on a Carreau fluid in a tube: Stimulation to avoid sediment accumulation in oil tanks

AbstractThe purpose of this article is to study the effect of electromagnetic (EM) stimulation on the pipeline, which has an electrical and thermal effect in addition to the chemical reaction on the crude oil and makes a sinusoidal wave on the wall. Modeling the crude oil as Carreau fluid is done. EM stimulation is an effective and safe technology that can be used to improve fluid movement in a variety of industrial applications. The flow analysis by applying EM may avoid blocking the crude oil pipeline which leads to a loss of production and capital investment. The basic partial differential equations of momentum, temperature and concentration are reduced to a system of nonlinear partial differential equations, which is solved numerically by using the Rung–Kutta–Merson method with Newton iteration in a shooting and matching technique under the assumption of long wavelength and the effect of physical implanted parameters is represented through charts for velocity, temperature, and concentration and numerical application.

Time Optimal and PID Controller for Armed Manipulator Robots

This paper is fourfold. First, three different time-optimal control problems for simulating the manipulator robots with one, two, and three arms are mathematically formulated. Corresponding to the related dynamical systems, the nonlinear system of ordinary differential equations is derived. It has been found that these problems are time-free and have distinct initial and boundary conditions, making them hard to solve. To find the minimum time with different payloads, a successful numerical method based on the finite difference and the three-stage Lobatto formula is proposed. Secondly, the related torquecontrolling problems are simulated, and then for one, two, and three armed manipulator robots, they are solved using the PID controller method. Thirdly, it is shown that, compared to the time-optimal controlling problem, the optimal PID torque controller solution takes more time to do the job than was expected. However, the solution in the PID controller method shows less oscillation than the time-optimal control problem. Fourthly, mathematical theories are used, and the numerical results for both methods and different payloads are compared.

Computational investigation of methanol-based hybrid nanofluid flow over a stretching cylinder with Cattaneo–Christov heat flux

Abstract This study investigates heat transfer rates in (AA7075-AA7072/Methanol) hybrid nanofluid flows, considering non-uniform heat sources and Cattaneo–Christov heat flux, with significant implications for aerospace engineering by enhancing thermal management in aircraft engines. The findings could revolutionize automotive cooling system efficiency, optimize heat dissipation in electronic devices, and advance the design of renewable energy systems such as concentrated solar power plants. The study aims to conduct a comparative analysis of (AA7075/Methanol) nanofluid and (AA7075-AA7072/Methanol) hybrid nanofluid flow, examining heat transfer rates, non-uniform heat sources, and Cattaneo–Christov heat flux theory around a stretching cylinder. Thermal radiation and the Biot number are also evaluated. Two different nanoparticles, AA7072 and AA7075, are used with methanol to create AA7075/Methanol nanofluid and AA7075-AA7072/Methanol hybrid nanofluid. The study compresses the resultant non-linear partial differential equation system and applies suitable similarity transformations to reduce the governing partial differential equations with boundary conditions to dimensionless form. The BVP4C shooting method in MATLAB is employed to numerically and graphically solve these dimensionless ordinary differential equations. The results indicate that higher curvature parameter values correlate with increased velocity and temperature distribution profiles. A rise in nanoparticle volume fraction reduces the radial velocity profile but increases the temperature profile. Temperature distribution profiles increase with higher thermal radiation parameter and Biot number values, while higher thermal relaxation parameter values decrease temperature. Additionally, thermal distribution profiles rise with increasing values of both the time-dependent heat source constant and space-dependent heat source parameter.

The Effect of a Riga Plate On Casson Hybrid Nanofluid Flow Along a Stretching Cylinder with an Exponential Heat Source and Thermal Radiation

This study investigates the effect of a Riga plate on the flow characteristics of a Casson hybrid nanofluid through a stretching cylinder embedded in a porous medium in the presence of an exponential heat source and thermal radiation. This model is used toexplore the potential applications of this analysis in the fields of cancer treatment and wound healing. The governing partial differential equations are converted into a system of nonlinear ordinary differential equations using suitable similarity transformations. The governing partial differential equations (PDEs) were reduced to ordinary differential equations (ODEs) using similarity variables, and the heat transfer phenomena, fluid flow dynamics, and nanoparticle behavior in the base fluid were captured through numerical analysis. The Casson fluid model was used to capture the non-Newtonian characteristics of blood-like fluids in the circulatory system, and stretching was employed to mimic the blood vessel. The modeling involved incorporating thermal radiation (Ra) and an exponential heat source, which are encountered in cancer treatment. The Runge-Kutta order 4 with shooting technique was used to obtain numerical solutions of the governing equations, and the effects of the Casson fluid parameter, curvature parameter, nanoparticle volume fraction, Eckert number, and radiation parameter were analyzed. The results showed that the Riga plate affected the flow patterns by increasing the temperature, and the temperature distribution was improved by increasing the radiation parameter and nanoparticle concentration. The insights from this study could be used to optimize heat-based treatments for cancer and wound healing, where control of fluid dynamics and heat transfer processes is critical.

Mathematical Model of Basal Stem Rot (BSR) Disease Spread in Oil Palm Plants

This study discusses the stability analysis of mathematical models for the spread of Basar Stem Rot(BSR) in oil palm plants. In developing this mathematical model, several assumptions are taken to obtain a model that is suitable for the spread of BSR disease. The resulting model is a system of first-order nonlinear differential equations with three variables. This research includes both analytical and numerical analysis. Analytical analysis includes determination of equilibrium points and local stability analysis, while numerical analysis is conducted using Microsoft Excel application. From this study, two equilibrium points were found with stability conditions that depend on the fulfillment of certain conditions. One important result obtained is that the equilibrium point will be locally stable if and only if α > μ and b > √D, where D is the discriminant of a quadratic equation. After analyzing analytically, the study continued with numerical simulations to illustrate and test the analytical results. Numerical results in the form of graphs show that the solution of the system is stable, which indicates that the disease will be endemic under certain conditions and time. This research provides a deeper understanding of the dynamics of the spread of BSR and the conditions that affect the stability of the spread of the disease. With this analysis, it is expected to contribute to efforts to control BSR disease in oil palm plants. In addition, this research also opens opportunities for the development of similar mathematical models for other plant diseases.

Mathematical Assessment of the Role of Human Behavior Changes on SARS-CoV-2 Transmission Dynamics in the United States.

The COVID-19 pandemic has not only presented a major global public health and socio-economic crisis, but has also significantly impacted human behavior towards adherence (or lack thereof) to public health intervention and mitigation measures implemented in communities worldwide. This study is based on the use of mathematical modeling approaches to assess the extent to which SARS-CoV-2 transmission dynamics is impacted by population-level changes of human behavior due to factors such as (a) the severity of transmission (such as disease-induced mortality and level of symptomatic transmission), (b) fatigue due to the implementation of mitigation interventions measures (e.g., lockdowns) over a long (extended) period of time, (c) social peer-pressure, among others. A novel behavior-epidemiology model, which takes the form of a deterministic system of nonlinear differential equations, is developed and fitted using observed cumulative SARS-CoV-2 mortality data during the first wave in the United States. The model fits the observed data, as well as makes a more accurate prediction of the observed daily SARS-CoV-2 mortality during the first wave (March 2020-June 2020), in comparison to the equivalent model which does not explicitly account for changes in human behavior. This study suggests that, as more newly-infected individuals become asymptomatically-infectious, the overall level of positive behavior change can be expected to significantly decrease (while new cases may rise, particularly if asymptomatic individuals have higher contact rate, in comparison to symptomatic individuals).

Stability of Periodic solutions by Krasnoselskii fixed point theorem of neutral nonlinear system of dynamical equation with Variable Delays

The fixed point theorem is used in this study to provide stability results for the zero solution of a nonlinear neutral system of differential equations with functional delay.

Semi-analytical solution of MHD free convective flow of a nanofluid over an exponentially stretching porous sheet in the presence of heat source/sink

In this study, the magnetohydrodynamic flow and heat transfer of electrically conducting nanofluid over an exponentially stretching porous sheet is investigated. By introducing similarity variables, the non-linear set of partial differential governing equations, are transformed into a system of ordinary nonlinear differential equations. The equations are solved semi-analytically by the Optimal Collocation Method (OCM) which is a powerful method to solve equations containing infinite boundary conditions. The accuracy of the OCM results is examined by the Richardson method. The effects of the different physical parameters such as magnetic field strength, nanofluid volume fraction, suction strength, and the heat sink/source value are discussed on the flow and heat transfer characteristics of the problem. The obtained results show that the velocity and temperature of nanofluid increase with the increase in volume fraction. Also, when suction increases, the thickness of the boundary layers decreases. As the strength of the magnetic field (Hartmann) increases, the temperature slightly increases but the flow rate of the nanofluid decreases. The skin friction coefficient and Nusselt number increase when the values of suction parameters, volume fraction, magnetic field, and the field angle increase.