Second-order Nonlinear Differential Equations Research Articles

Understanding the nonlinear reactive transport model in porous catalysts

This paper discusses the steady-state nonlinear reaction-diffusion processes in porous catalysts. This model is based on a nonlinear second-order differential equation that includes a nonlinear term associated with the Michaelis-Menten and non-Michaelis-Menten kinetics of the reaction. The nonlinear equations can be approximately solved using the Taylors series, modified Taylors series and Akabri-Ganji technique to obtain the concentration of dissolved species. The influence of the half-saturation parameter and the characteristic reaction rate on concentration is explored. Sensitivity analysis of parameters is discussed. Our analytical findings were compared with numerical solutions.

Stability and boundedness criteria for certain second-order nonlinear neutral stochastic functional differential equations

This paper presents stochastic stability and stochastic boundedness to certain second-order nonlinear neutral stochastic differential equations. The second-order differential equation is weakened to a neutral stochastic system of first-order equations and used together with a second-order quadratic function to obtain perfect Lyapunov-Krasovskii functional. This functional is adapted and applied to obtain criteria on the nonlinear functions to ensure novel results on stochastic stability and stochastic asymptotic stability of the zero solution. Furthermore, when the forcing term is nonzero, fresh results on stochastic boundedness and uniform stochastic boundedness of solutions are obtained. The results of this paper are original, new, essentially improving, complementing, and simplifying several related ones in the literature. Two special cases of the theoretical results are supplied to demonstrate the applicability of the hypothetical results.

Existenceof Heteroclinic Solutions in Nonlinear Differential Equations of the Second-Order Incorporating Generalized Impulse Effects with the Possibility of Application to Bird Population Growth

This work considers the existence of solutions of the heteroclinic type in nonlinear second-order differential equations with ϕ-Laplacians, incorporating generalized impulsive conditions on the real line. For the construction of the results, it was only imposed that ϕ be a homeomorphism, using Schauder’s fixed-point theorem, coupled with concepts of L1-Carathéodory sequences and functions along with impulsive points equiconvergence and equiconvergence at infinity. Finally, a practical part illustrates the main theorem and a possible application to bird population growth.

Periodic and subharmonic solutions in the motion of a bead on a rotating circular hoop

We establish the necessary conditions for the existence and multiplicity of periodic and subharmonic solutions to a second-order nonlinear ordinary differential equation (ODE). This ODE describes the motion of a bead on a rotating circular hoop subjected to a constant angular velocity ω and a T-periodic forcing. Our approach involves estimating bounds for the angular velocity and period using upper and lower solution methods.

Oscillatory behavior of solutions of second-order non-linear differential equations with mixed non-linear neutral terms

This study primarily seeks to expand upon these developments by encompassing neutral differential equations of mixed type, incorporating both delay and advanced terms, particularly in the case of the canonical operator. The presented results are derived from the application of the comparison method, Riccati transformation, and integral averaging technique. These methodologies lead to substantial improvements and extensions of existing results found in the literature. Additionally, illustrative examples are provided to demonstrate the practical implications of the developed criteria.

Minimum Principles for Sturm–Liouville Inequalities and Applications

A minimum principle for a Sturm–Liouville (S-L) inequality is obtained, which shows that the minimum value of a nonconstant solution of a S-L inequality never occurs in the interior of the domain (a closed interval) of the solution. The minimum principle is then applied to prove that any nonconstant solutions of S-L inequalities subject to separated inequality boundary conditions (IBCs) must be strictly positive in the interiors of their domains and are increasing or decreasing for some of these IBCs. These positivity results are used to prove the uniqueness of the solutions of linear S-L equations with separated BCs. All of these results hold for the corresponding second-order differential inequalities (or equations), which are special cases of S-L inequalities (or equations). These results are applied to two models arising from the source distribution of the human head and chemical reactor theory. The first model is governed by a nonlinear S-L equation, while the second one is governed by a nonlinear second-order differential equation. For the first model, the explicit solutions are not available, and there are no results on the existence of solutions of the first model. Our results show that all the nonconstant solutions are increasing and are strictly positive solutions. For the second model, many results on the uniqueness of the solutions and the existence of multiple solutions have been obtained before. Our results are applied to prove that all the nonconstant solutions are decreasing and strictly positive.

ALGORITHM FOR ESTIMATING THE PARAMETRIC ROLL MOTIONS

For a ship, both synchronous and parametric roll motions can be modelled by nonlinear second-order differential equations with the roll angle as dependent variable and the nonlinearities coming from the damping and restoring moments. In the absence of exact solution techniques, approximate solutions for such equations can be obtained, in principle, analytically or numerically. In this paper, we focused on a fast and accurate recursive scheme based on the segmentation of the time domain into a fine grid of intervals of equal lengths. Due to the smallness of each time subinterval, a reset of the terms of the initial nonlinear equation of motion allows replacing it with a linear one in which the constant coefficients represent approximate values of the former variable coefficients in the middle of each subinterval. The linear equation is solved via the Laplace transform and the resulting solution together with its first derivatives is used to generate the iterative algorithm. This technique was applied to the case of two typical roll equations. The first describes the synchronous roll of a vehicle ferry model equipped or not with bilge keels while the second estimates the parametric roll of a container ship model. In both cases, the analyzed recursive algorithm not only generated results in full agreement with those provided by the ode45 solver in Matlab but also lead to a gain in computation time and offered more flexibility in imposing some restrictive conditions associated, for example, with exceeding some oscillation amplitudes or to the ship capsizing. Key words: synchronous and parametric ship rolling, iterative scheme

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

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

Attached Flows for Reaction–Diffusion Processes Described by a Generalized Dodd–Bullough–Mikhailov Equation

This paper uses the attached flow method for solving nonlinear second-order differential equations of the reaction–diffusion type. The key steps of the method consist of the following: (i) reducing the differentiability order by defining the first derivative of the variable as a new variable called the flow and (ii) a forced decomposition of the derivative-free term so that the flow appears explicitly in it. The resulting reduced equation is solved using specific balancing rules. Only step (i) would lead to an Abel-type equation with complicated integral solutions. Completed with (ii) and with the graduation procedure, the attached flow method used in the paper, without requiring such a great effort, allows for the obtaining of accurate analytical solutions. The method is applied here to a subclass of reaction–diffusion equations, the generalized Dodd–Bulough–Mikhailov equation, which includes a translation of the variable and nonlinearities up to order five. The equation is solved for each order of nonlinearity, and the solutions are discussed following the values of the parameters involved in the equation.

On the Oscillatory Behavior of Solutions of Second-Order Non-Linear Differential Equations with Damping Term

In this paper, we discuss the oscillatory behavior of solutions of two general classes of nonlinear second-order differential equations. New criteria are obtained using Riccati transformations and the integral averaging techniques. The obtained results improve and generalize some recent criteria in the literature. Moreover, a traditional condition is relaxed. Three examples are given to justify the results.

Self-similar solutions to a multidimensional singular heat equation with power nonlinearity

The paper deals with the construction of exact solutions to a singular heat equation with power nonlinearity in the case of numerous independent variables with spatial (e.g. axial or central) symmetry. A new class of self-similar solutions is proposed, which reduce to solving the Cauchy problem for a second-order nonlinear ordinary differential equation having singularities at the higher derivative with respect to the required function and/or the independent variable. The ordinary differential equation is studied in two ways: analytically and numerically. The analytical study uses a truncated Taylor series with recurrently computed coefficients, for which explicit formulas are obtained. The numerical solution to the problem uses an iteration algorithm based on the collocation method and radial basis functions. The numerical analysis shows the convergence of the proposed numerical algorithm and its sufficient accuracy enabling one to use the found self-similar solutions to verify approximate solutions to the original heat equation. Besides, the numerical analysis has allowed the radius of convergence of the constructed Taylor series to be evaluated. The form of the constructed self-similar solutions, namely their unboundedness near the symmetry center (axis), enables us to study the behavior and exactness of the numerical solutions to the nonlinear singular parabolic-type equation that have been obtained by the stepwise solution method proposed by us earlier and possess the same property.

Asymptotic analysis of electrohydrodynamic flow through a swarm of porous cylindrical particles

The present article reveals the study of an electrohydrodynamic flow through a membrane composed of a swarm of porous layered cylindrical particles adopting a heat transfer approach. The configuration of the proposed theoretical model is segregated into two regions in which the region proximate to the solid core of the cylindrical particle is a porous region. However, a region surrounded by a porous region is a non-porous (clear fluid) region. The thermal equations are employed under steady-state conditions to establish the temperature distribution when heat conduction prevails over heat convection. The Brinkman and Stokes equations regulate fluid flow through a swarm of porous layered cylindrical particles in porous and non-porous regions, respectively. With the purpose of addressing an electric field in the fluid flow process through a swarm of porous layered cylindrical particles to understand the role of a Hartmann electric number, the momentum equation and the charge density are coupled and nonlinear. The nonlinear second-order differential equation governs the momentum equation and regulates fluid flow through a swarm of porous cylindrical particles. The solutions of the energy equations for both regions are analytically obtained. The asymptotic expansions of velocities for porous and non-porous regions have been derived using the perturbation technique for the small and large values of the nonlinearity parameter α. The effects of various parameters like Hartmann electric number, Grashof number, radiation parameter, viscosity ratio parameter, and porosity of the porous material on the hydrodynamical permeability, Kozeny constant of the membrane, and temperature are analyzed graphically. A noteworthy observation is that a rising Hartmann electric number, the ratio of electric force to the viscous force, enhances the velocity, which is relatively more significant for higher permeability and hence enhances the membrane permeability; however, decay in Kozeny constant is reported with a rising Hartmann electric number. Significant velocity and membrane permeability growth are described with a rising Grashof number, a ratio of thermal buoyancy and viscous forces. The observations from the present study hold promise for advancing our understanding of critical physical and biological applications, including wastewater treatment filtration processes, petroleum reservoir rocks, and blood flow through smooth muscle cells.

Periodic Solutions to Nonlinear Second-Order Difference Equations with Two-Dimensional Kernel

In this work, we provide conditions for the existence of periodic solutions to nonlinear, second-order difference equations of the form y(t+2)+by(t+1)+cy(t)=g(y(t)), where b and c are real parameters, c≠0, and g:R→R is continuous.

Non-linear electro-rheological model of a membrane immersed in Tanner-Power law fluids applied to outer hair cells: Shear-thinning mechanisms

Flexoelectric actuation employs an applied electric field to induce membrane curvature, which is the mechanism utilized by the outer hair cells (OHC) present in the inner ear. The model developed for this study, representing the OHC, integrates two key components: (i) an approximation of the flexoelectric membrane shape equation for circular membranes attached to the inner surface of a circular capillary, and (ii) the coupled capillary flow of contacting liquid viscoelastic phases characterized by the Tanner-Power law rheological equation of state. A second-order non-linear differential equation for average curvature has been derived, and a robust numerical method has been programmed. This model simplifies to a linear model used previously. The main challenge involves identifying and describing the enhancement in curvature change rate. It was observed that low symmetry, low viscosity, and soft membrane and shear-thickening behavior of the phases enhance the curvature change rate. Additionally, there exists a critical electric field frequency value that maximizes the curvature change rate (resonance effect). The current theory, model, and computational simulations add to the ongoing development comprehension of how biological membrane shape actuation through electromechanical couplings.

On the problem of solvability of nonlinear boundary value problems for shallow isotropic shells of Timoshenko type in isometric coordinates

The solvability of a boundary value problem for a system, which describes the equilibrium state of elastic shallow inhomogeneous isotropic shells with loose edges referred to isometric coordinates in the Timoshenko shear model and consists of five non-linear second-order partial differential equations under given non-linear boundary conditions, is studied. The boundary value problem is reduced to a nonlinear operator equation for generalized displacements in Sobolev space, the solvability of this equation is established with the help of the contraction mapping principle.

New Conditions for Testing the Oscillation of Solutions of Second-Order Nonlinear Differential Equations with Damped Term

This paper deals with the oscillatory behavior of solutions of a new class of second-order nonlinear differential equations. In contrast to most of the previous results in the literature, we establish some new criteria that guarantee the oscillation of all solutions of the studied equation without additional restrictions. Our approach improves the standard integral averaging technique to obtain simpler oscillation theorems for new classes of nonlinear differential equations. Two examples are presented to illustrate the importance of our findings.

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

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

Dynamic Model of Lower Limb Motion in the Sagittal Plane during the Gait Cycle

Context: This work presents the development of a dynamic model for human lower limb motion in the sagittal plane during the gait cycle. The primary objective of this model is to serve as a powerful tool for the design of rehabilitation and assistive devices, such as exoskeletons, prostheses, and orthoses. It achieves this by facilitating the estimation of joint torques, the detailed analysis of kinematic variables, optimal actuator selection, and the exploration of advanced control techniques. Method: The dynamic model consists of two primary components: (1) the plant model and (2) a closed-loop controller. The plant model represents the forward dynamics of human gait and is based on a multi-mass pendulum composed of three segments of the lower limb (thigh, lower leg, and foot) and three joints (hip, knee, and ankle). It is analyzed using the Euler-Lagrange formulation and the nonlinear second-order differential equations are implemented in MATLAB’s Simulink. To reproduce reference human gait trajectories and simulate the functioning of the neuromusculoskeletal system and the central nervous system, a closed-loop PID controller is incorporated into the plant model. It is noteworthy that the scope of this dynamic model is specifically confined to the sagittal plane. Results: The dynamic model is evaluated in terms of angular displacement tracking using the relative maximum error (RME) and the root mean square error (RMSE) for reference trajectories of healthy adult male human gait as reported in the literature. The model demonstrates tracking with errors below 2.2 [°] in magnitude and 3,5% for all three considered segments (thigh, lower leg, and foot). Conclusions: The quantitative results show that the dynamic model developed in this work is reliable and allows for a precise reproduction of human gait trajectories.

Matrix factorizations for the generalized Charlier and Meixner orthogonal polynomials

The Cholesky factorization of the moment matrix is considered for the generalized Charlier and generalized Meixner discrete orthogonal polynomials. For the generalized Charlier, we present an alternative derivation of the Laguerre–Freud relations found by Smet and Van Assche. Third-order and second-order nonlinear ordinary differential equations are found for the recursion coefficient γn, that happen to be forms of the Painlevé deg-PV in disguise. New Laguerre–Freud relations are also found for the generalized Meixner case, which are compared with those of Smet and Van Assche.

Qualitative analysis on the electrohydrodynamic flow equation

<abstract><p>In this paper, we present a comprehensive analysis of the lower and upper bounds of solutions for a nonlinear second-order ordinary differential equation governing the electrohydrodynamic flow of a conducting fluid in cylindrical conduits. The equation describes the radial distribution of the flow velocity in an "ion drag" configuration, which is affected by an externally applied electric field. Our study involves the establishment of rigorous analytical bounds on the radial distribution, taking into account the Hartmann number $ H $ and a parameter $ \alpha. $ An analytic approximate solution is obtained as an improvement of the a priori estimates and it is found to exhibit strong agreement with numerical solutions, particularly when considering small Hartmann numbers. Further, estimates for the central velocity $ w(0) $ of the fluid occurring at the center of the cylindrical conduit were also established, and some interesting relationships were found between $ H $ and $ \alpha. $ These findings establish a framework that illuminates the potential range of values for the physical parameter within the conduit.</p></abstract>