Nonlinear Ordinary Differential Equations Research Articles

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.

Significance of magnetic field and Darcy–Forchheimer law on dynamics of non-Newtonian hybrid nanofluid flow over a spinning disc with Arrhenius activation energy and shape factor

Purpose The purpose of the study is to explore the three-dimensional heat and mass transport dynamics of the magneto-hydrodynamic non-Newtonian (Casson fluid) hybrid nanofluid flow comprised of − as nanoparticles suspended in base liquid water as it passes through a flexible spinning disc. The influence of a magnetic field, rotation parameter, porosity, Darcy−Forchheimer, Arrhenius’s activation energy, chemical reaction, Schmidt number and nanoparticle shape effects are substantial physical features of the investigation. Furthermore, the influence of hybrid nanofluid on Brownian motion and thermophoresis features has been represented using the Buongiorno model. The novelty of the work is intended to contribute to a better understanding of Casson non-Newtonian fluid boundary layer flow. Design/methodology/approach The governing mathematical equations that explain the flow and heat and mass transport phenomena for fluid domains include the Navier−Stokes equation, the thermal energy equation and the solutal concentration equations. The governing equations are expressed as partial differential equations, which are then converted into a suitable set of non-linear ordinary differential equations by using the necessary similarity variables. The ordinary differential equations are computed by combining the shooting operation with the three-stage Lobatto BVP4c technique. Findings Graphs and tables are used in the process of analysing the characteristics of velocity distributions, temperature profiles and solutal curves at varying values of the parameters, along with friction drag, heat transfer rate and Sherwood number. It has been revealed that the radial and axial velocities decrease when the Casson parameter value increases and that the rate of heat transmission is higher in hybrid nanofluids with nanoparticles in the shape of a blade. The increase in Brownian motion and thermophoresis parameters causes a rise in the temperature profile. Also, an increase in the activation energy parameter improves the solutal curve. The use of nanoparticles was shown to improve extrusion properties, the rotary heat process and biofuel generation. Originality/value All results are presented graphically and all physical quantities are computed and tabulated. The current results are compared to previous investigations and found to agree significantly with them.

Combined fluid and ferrofluid buoyancy force, heat absorption and non linear ferro-viscosity effects on ferro- hydrodynamics fluid flow in orthogonal porous surface

Mixed convection incompressible ferro-hydrodynamic (FHD) boundary layer fluid flow on vertical porous curvilinear orthogonal surfaces is examined in our present study in presence of combined fluid and ferro-fluid buoyancy force, nonlinear ferro-viscosity parameter, and heat absorption/generation effects. Design/Methodology/Approach: Types of Prandtl momentum and thermal boundary layer partial differential equations in a curvilinear system are converted into non-dimensional ordinary nonlinear differential equations (ODE) by introducing dimensionless velocities, temperature functions and similarity variable. Missing initial conditions are found by shooting method and solving the system of nonlinear ODEs by Runge-Kutta integration scheme of order six. With the help of MATLAB, the numerical solutions are graphically displayed in the forms of primary velocity profile, secondary velocity profile, and temperature profile. Three types of magnetic nano-ferrofluid are considered such as water-based ferrofluid (Taiho W-40), Hydrocarbon based ferrofluid (EMG-901) and kerosene-based ferrofluid (C1-20B) whose Prandtl numbers are 44.3, 79.3 and 128 respectively shown in Table 1. For kerosene based ferro-fluid, Table 3 presents the coefficient of skin friction and heat flux. Finding: The effects of combined fluid and ferrofluid buoyancy forces, nonlinear ferro-viscosity parameter, suction/injection, and heat absorption/generation parameters on the velocity and temperature profiles are investigated. Our physical interest is to calculate the local shearing stresses and the local heat flux at the surface. Finally, we have made comparisons of the results which are highlighted the validity of our numerical calculations adopted for the presence study.

Boundary crossing problems and functional transformations for Ornstein–Uhlenbeck processes

ABSTRACT We are interested in the law of the first passage time of an Ornstein–Uhlenbeck process to time-varying thresholds. We show that this problem is connected to the laws of the first passage time of the process to members of a two-parameter family of functional transformations of a time-varying boundary. For specific values of the parameters, these transformations appear in a realisation of a standard Ornstein–Uhlenbeck bridge. We provide three different proofs of this connection. The first is based on a similar result for Brownian motion, the second uses a generalisation of the so-called Gauss–Markov processes, and the third relies on the Lie group symmetry method. We investigate the properties of these transformations and study the algebraic and analytical properties of an involution operator which is used in constructing them. We also show that these transformations map the space of solutions of Sturm–Liouville equations into the space of solutions of the associated nonlinear ordinary differential equations. Lastly, we interpret our results through the method of images and give new examples of curves with explicit first passage time densities.

Classification of the non-null electrovacuum solution of Einstein–Maxwell equations with three-parameter abelian group of motions

The classification of the Stackel spaces of the electrovacuum of the type (3.0) has been done. These spaces are invariant under the action of the three-parameter abelian group of motions and belong to the first type Bianchi spaces. In the case of a non-zero cosmological term, the metrics and potentials contain solutions of a nonlinear ordinary differential equation of the second order. When the cosmological term equals zero, the metrics and the components of the electromagnetic field tensor are expressed through elementary functions. Thus the classification of the electrovacuum stackel spaces of all types is completed and complete list of these spaces is constructed.

Features of microorganism and two-phase nanofluid in a tangent hyperbolic Darcy-Forchhiemer flow induced by a stretching sheet with Lorentz forces

This work examines the two-dimensional tangent hyperbolic flow over a stretching sheet with a uniform magnetic field. The Buongiorno model is utilized to analyze and explain the spread of uneven coefficients in the presence of gyrotactic microorganisms. The concept of microorganisms and the resulting bioconvection enhance the stability of the nanoparticles. The impacts of thermal radiation, heat sources, convective heating, and chemical reactions are also evaluated. The suggested mathematical problem results in a nonlinear set of partial differential equations (PDEs), which are subsequently reduced to ordinary differential equations (ODEs) by applying the appropriate transformation. The resultant highly nonlinear ordinary differential equations (ODEs) are numerically solved using MATLAB's built-in package known as bvp4c. An in-depth investigation into the changes in the velocity field, the temperature profile, the concentration of nanoparticles profile, and the motile density profile is analyzed through graphs against various influencing parameters. Additionally, computations of engineering interest quantities such as skin friction, local Nusselt number, local Sherwood number, and molecular density are presented in both graphical and tabular formats for further examination. It has been explored that with greater values of the Weissenberg number, the fluid velocity upsurges when n<1, whereas the opposite behavior is noticed when n>1. It is also noted that an increment in the Peclet number decreases motile density for both dilatants n>1 and pseudoplastic n<1 fluids. The computed results are compared with existing literature in limiting cases and found good agreement.

Numerical simulation of an effective transform mechanism with convergence analysis of the fractional diffusion-wave equations

In the current study, we solve two very important mathematical models, such as the time fractional-order space-fractional telegraph and diffusion-wave equations using a reliable technique called the Adomian decomposition natural method (ADNM), which combines Adomian decomposition and natural transform. The diffusion wave equation describes the flood wave propagation, which is used in solving overland and open channel flow problems. For this reason, it is critical to fully understand and effectively solve the diffusion wave equations. Because telegraph equations are crucial for modeling and developing voltage or frequency transmission, they are widely used in physics and engineering. Furthermore, the designing process is greatly impacted by the uncertainty in the system parameters. For nonlinear ordinary differential equations based on the theorem of Banach fixed point, we provide existence and uniqueness theorem proofs. The present approach has been successfully used to explore exact solutions for time fractional-order and space fractional-order applications. The results show how effective and valuable the ADNM. This paper presents a methodology that will be used in future work to address similar nonlinear problems related to fractional calculus.

On the nonlinear dynamics of in-contact rigid bodies experiencing stick–slip and wear phenomena

In this paper, the dynamic behavior of one degree-of-freedom oscillator subject to stick–slip and wear phenomena at the contact interface with a rigid substrate is investigated. The motion of the oscillator, induced by a harmonic excitation, depends on the tangential contact forces, exchanged with the rigid soil, which are modeled through piecewise nonlinear constitutive laws, accounting for stick–slip phenomena due to friction as well as wear due to abrasion, already developed by the authors in a previous work. The nonlinear ordinary differential equations governing the problem are derived, whose solution is numerically obtained via a typical Runge–Kutta-based algorithm. The main target of this study is to analyze and discuss the strong nonlinear behavior, descending from the presence of stick–slip and wear phenomena, thus investigating the effect of the different interface modeling. In this framework, the analysis is carried out considering the whole evolution of non-smooth contact laws, starting from the virgin interface.

Nonlinear geometrically exact dynamics of hyperelastic pipes conveying fluid: Comparative study of different hyperelastic models

This work presents a comparative analysis of the nonlinear dynamic responses of fluid-conveying pipes modelled by four different phenomenological hyperelastic constitutive models, namely the Mooney-Rivlin model, the Neo-Hookean model, the Yeoh model, and the Ogden model. Based on the strain energy density functions and the exact geometric relation, the geometrically exact dynamic models corresponding to different hyperelastic contitutive models are established. Then the Galerkin's truncation is adopted to transfer the obtained nonlinear integral-partial differential equations into the nonlinear ordinary differential equations. Subsequently, linear dynamic study and finite difference method are taken to explore the natural frequencies and stabilities and nonlinear dynamic responses of the pipe, respectively. And the fitting coefficients for various models obtained from Treloar's experiment data are utilized in the numerical calculation. It is found the Hopf bifurcation and the limit cycle oscillations in the post-flutter status are influenced by the choice of the hyperelastic models, particularly under large flow velocity. However, hyperelastic constitutive models do not exert a significantly extreme influence on the results, which can be attributed to the relatively small strains even though large deformations occur in the pipe. Mathematically, it is proven that the symmetry observed in the pipe's displacement response stems from the symmetry of its cross-section and the inextensibility assumption of the centerline, with hyperelasticity having no effect on this symmetry.

Effects of Arrhenius activation energy and melting heat transfer on Carreau fluid through a stretching sheet

ABSTRACT This article investigates the effect of activation energy on Carreau fluid flow in the presence of variable thermal conductivity, inclined magnetic field, melting heat transfer, and Soret and Dufour phenomena. A set of suitable similarity transformations is applied for transforming the governing partial differential equations (PDEs) considered for the physical phenomena into non-linear ordinary differential equations (ODEs). We obtained the results by solving the non-linear ODEs numerically which describe the behavior of the velocity, temperature, and concentration profiles of the fluid against the governing parameters and illustrated these graphically. The velocity distribution decreases for angle of inclination while it is increasing against higher melting parameter. The temperature distribution enhances with higher modified Dufour parameter and Prandtl number while it declines for thermal conductivity, activation energy, and melting parameters. The concentration distribution increases for melting parameter, Soret number, and activation energy parameter while decreasing for temperature difference parameter, Schmidt number, and reaction rate parameter. Skin friction, Nusselt and Sherwood numbers are evaluated numerically against the various governing parameters. Accuracy of the present work is illustrated and established through comparison carried out for skin friction versus magnetic parameter with existing results in the literature.

Significant statistical model of heat transfer rate in radiative Carreau tri-hybrid nanofluid with entropy analysis using response surface methodology used in solar aircraft

The recent advancement of technologies focuses on the use of solar-based heat radiation and nanotechnology. The primary source of heat is solar energy, which is obtained by absorbing sunlight. In addition, solar thermal planes, photovoltaic cells, and sun-based hybrid nanofluids all help to harness solar energy. Therefore, main concern of the study is to explored the two-dimensional Engine oil-based Carreau tri-hybrid nanofluid (SiO2 (silica dioxide), MOS2 (Molybdenum disulfide), Fe3O4(black iron oxide)) flow via a stretching sheet within solar wings. The flow phenomena of the non-Newtonian fluid enrich for the effects of radiant energy and the external heat source. By adopting similarity transfiguration, the dimensional partial main flow equations are changed into nonlinear dimensionless ordinary differential equations. The most effective and powerful tool of the semi-analytical method ADM (Adomain Decomposition Method) is utilized to solve the nanofluid flow problem. The uniqueness of the proposed study arises due to the analysis of entropy generation obtained by various irreversibility processes. The graphical illustration of outcome of various governing parameters on velocity, temperature, entropy terms and engineering quantities have been presented and theoretically analyzed. The findings reveal that The Carreau ternary hybrid nanofluids enhance the system's thermal conductivity for optimal temperature regulation in solar aircraft and better heat dissipation. The use of Carreau ternary hybrid nanofluids significantly reduces entropy generation, which implies a more reversible process, enhancing the overall efficiency of the energy conversion systems in solar aircraft. The heat flow of the system is accurately reflected by the chosen variables and their interactions, as demonstrated by the constructed model's strong statistical significance. For real-world applications, this improves the model's predictability and dependability.

Uncovering Hidden Patterns: Approximate Resurgent Resummation from Truncated Series

We analyze truncated series generated as divergent formal solutions of non-linear ordinary differential equations. Motivating the study is a specific non-linear, first-order differential equation, which is the basis of the resurgent formulation of renormalized perturbation theory in quantum field theory. We use the Borel–Padé approximant and classical analysis to determine the analytic structure of the solution using the first few terms of its asymptotic series. Afterward, we build an approximant, consistent with the resurgent properties of the equation. The procedure gives an approximate expression for the Borel–Ecalle resummation of the solution useful for practical applications. Connections with other physical applications are also discussed.

Phase trajectories, chaotic behavior, and solitary wave solutions for (3+1)-dimensional integrable Kadomtsev–Petviashvili equation in fluid dynamics

This paper focuses on finding exact soliton solutions and further examining the qualitative characteristics of these solutions for a nonlinear partial differential equation known as the (3+1)-dimensional integrable Kadomtsev–Petviashvili equation that portrays a unique dispersion effect. The (3+1)-dimensional integrable Kadomtsev–Petviashvili equation models the evolution of weakly nonlinear dispersive waves in three spatial dimensions “x,y,z” and time “t”. It balances nonlinearity and dispersion, leading to solitary wave solutions that are stable and localized. This equation is crucial in understanding complex wave phenomena in fields like fluid dynamics, plasma physics, and optics. Initially, the traveling wave transformation is employed to convert the nonlinear partial differential equations into the non-linear ordinary differential equations. Then, we use the Riccati equation approach to find solitary wave solutions for the resulting ordinary differential equation. To better comprehend the physical significance of these derived solutions, we illustrate them using various visual representations. Diverse novel collections of solitary wave solutions are obtained such as bright solitons, dark solitons, and dark singular solitons. Second, we convert the nonlinear ordinary differential equation into a linear dynamical system using the planar dynamical transformation. Then, we undertake a qualitative analysis of the dynamical system to study its chaotic characteristics and bifurcation phenomena. Phase portraits of bifurcation are observed at the equilibrium points of the planner dynamical system. We discuss various tools to identify the chaos (random and unpredictable behavior) in dynamical systems. Graphical visualization of qualitative analysis of the system is also depicted in 3-D phase portraits, 2-D phase portraits, Poincaré mapping, and time series. Additionally, the Runge–Kutta technique is used to do a thorough sensitivity analysis of the dynamical system. Using this analytical procedure, it is verified that the stability of the solution has negligible impact while there is little change in the beginning circumstances. This paper aims to provide vital views to scientists and researchers seeking to advance their experimental operations. Furthermore, this research might be expanded by incorporating solitons and situations and investigating wave dynamics. The methods employed provide a direct and uncomplicated approximation of all solutions when compared to the previously established techniques. Various combinations and values of the physical parameters are utilized to explore the optical soliton solutions of the resultant system.

An exact solution of the lubrication equations for the Oldroyd-B model in a hyperbolic pipe

An exact analytical solution of the lubrication equations for the steady, isothermal, incompressible flow of a viscoelastic Oldroyd-B fluid in a hyperbolic cylindrical contracting pipe is derived. The solution is valid for values of the Deborah number, De, up to order unity (De is defined as the ratio of the longest relaxation time of the polymer to the characteristic residence time of the fluid in the pipe), all values of the ratio of the polymer viscosity to the total viscosity of the fluid, η, and typical values of the contraction ratio, Λ, encountered in experiments and practical applications. It is provided in terms of the streamfunction only and is used in the momentum balance to derive a strongly non-linear ordinary differential equation of second order with unknown a function which corresponds to a modified fluid velocity along the main flow direction. The final equation is solved semi-numerically using a fully spectral (Legendre)-Galerkin approach to resolve the unknown function almost down to machine accuracy. The exact solution for the polymer extra-stresses, which is emphasized is not the full solution of the complete lubrication equations, allows for the derivation of a variety of theoretical expressions for the average pressure-drop along the pipe. In all cases, a decrease in the pressure drop compared to the Newtonian value with increasing De, η and/or Λ is predicted. The differences between the corresponding analytical solution for the planar geometrical configuration are also identified and discussed.

Tangent hyperbolic MHD nanoliquids on non-isothermal stretched sheets: Analyzing the impact of transport parameters, variable fluid properties and convective boundary conditions

The key focus of this research is to look at how tangent hyperbolic magnetohydrodynamic (MHD) nano-liquids move in a non-isothermal coagulated stretched sheet. In this study, we intend to explore how fluid behavior responds to alternations in transport physical parameters. The coagulated sheet is subjected to two distinct boundary conditions to inspect the heat and mass transport rate of the nano liquid. The convective heat boundary condition (CBC) and mass-convective boundary condition (MCBC) are employed to specify the physical conditions at the boundaries of the problem domain. The boundary conditions, responsible for regulating heat flow and mass transfer rates, play a crucial role in shaping the liquid’s behavior. Using a similarity solution, the partial differential equation describing the flow of mass and heat within the system transforms into an interconnected network of non-linear ordinary differential equations. These mathematical equations are subsequently figured out using a numerical finite difference method. This study additionally examines the correlation between thermophoresis and Brownian motion, two fundamental concepts in the dynamics of colloidal suspensions. The study’s findings indicate that the temperature profile increases in all scenarios when the variable thermal conductivity and variable viscosity parameters are increased. In contrast, for the same parameters, the velocity profile is decreased. Further, increasing wall thickness reduces heat dissipation; consequently, the temperature profile similarly affects the velocity power index. The Biot number improves the rates of temperature transfer. These results underscore the significant influence of these parameters in predicting the behavior of MHD tangent hyperbolic nanofluids. This study elucidates the intricate interaction among different physical parameters in the dynamics of nano liquids, which holds significant implications for various industrial applications.

Convective MHD flow of Casson fluid through porous medium using the Darcy‐Forchheimer model under the impact of prescribed heat sources

AbstractBlood and other non‐Newtonian fluids are most appropriately described by a rheological model known as a Casson fluid. The presence of yield stress in Casson fluids makes them particularly relevant in the fields of biomechanics and polymer industries. In this analysis, a numerical simulation of flow of a Casson nanofluid over a stretching surface was estimated. The flow of fluid was subjected to a magnetic field of varying strength. The governing partial differential equations for momentum and heat transmission are converted into a set of nonlinear ordinary differential equations with the aid of similarity transformations. After that, numerical techniques are used to solve the equations. Figures and tables are utilized to illustrate how non‐dimensional parameters affect energy, concentration, and velocity trends.

Multilayer deep-learning intelligent computing for the numerical analysis of unsteady heat and mass transfer in MHD Carreau Nanofluid model

The objective of this article is to investigate the mechanism of unsteady heat and mass transfer (HMT) in a magnetohydrodynamics Carreau nanofluid flow (MHD-CNFF). To accomplish this, we utilized a novel artificial intelligence-based (AI) method, deep learning predictive technique combined with Levenberg-Marquardt scheme (DLPT-LMS). The similarity transformation is applied to convert the partial differential equations (PDEs) into a set of non-linear ordinary differential equations (ODEs). We numerically solved these resulting ODEs using the Adam numerical approach in sophisticated “Mathematica” software. The influence of Brownian motion parameter, thermophoresis parameters, Prandtl number, Lewis number, and mass transfer parameters on the temperature, velocity, and concentration profile is analyzed using a synthetic dataset. This evaluation is conducted across different variation of MHD-CNFF by leveraging DLPT-LMS. The absolute error (AE) across various iteration of MHD-CNFF demonstrated a progressive drop, indicating the accuracy of DLPT-LMS. The outcomes of this investigation revealed that the AE ranged from 10-3 to 10-8 which shows a minimum error between actual and predicted value. Additionally, the MSE confirmed the accuracy of proposed DLPT-LMS in predicting the behavior of velocity, temperature, and concentration profiles, with lower values showing better predictions. The model accuracy was reinforced by performance results, including regression analysis and error histograms. Error values gradually decreased during training, validation, and testing phases, exhibiting a robust convergence and indicating the highly reliability of AI-based algorithm for investigating intricate fluid dynamics in MHD-CNFF systems.

Revealing the dynamics of stagnant rings of third-grade fluid film with heat transfer in the presence of surface tension

In many engineering applications, including coating and lubrication operations, analyzing the temperature behavior of thin film flows on a vertically upward-moving tube is crucial to improving predictive models. This paper examines a steady third-grade fluid film flow with a surface tension gradient on a vertical tube. The mechanisms responsible for the fluid motion are upward tube motion, gravity, and surface tension gradient. This analysis focuses on heat transfer and stagnant ring dynamics. The formulated highly nonlinear ordinary differential equations are solved using the Adomian decomposition method. The conditions for stagnant rings and uniform film thickness are attained and discussed. The inverse capillary number C, Stokes number St, Deborah number De, and Brinkman number Br emerged as flow control parameters. The temperature of the fluid film rises with an increase in the C, St, De, and Br, whereas it decreases with an increase in thermal diffusion rate. The radius of stagnant rings tends to shrink by the increase in C, St, and De. When the value of De is high, third-grade fluid behaves like solids; only free drainage happens with smaller radius stagnant rings and high temperatures. A comparison between Newtonian and third-grade fluids regarding surface tension, velocity, temperature, stationary rings, and fluid film thickness is also provided.

Conservative deep neural networks for modeling competition of ribosomes with extended length

We develop a network model that combines several ribosome flow models with extended objects (RFMEO) competing for the finite pool of ribosomes. This alleviates the need to systematically coarse-grain the mRNA molecules. The dynamics of the network is described by a system of non-linear ordinary differential equations. It is shown that the network always converges to a steady state for a fixed number of ribosomes. Our analysis shows that increasing any of the transition rates along an RFMEO increases its output rate and either the output rates of the other RFMEOs all increase or all decrease. Simulations also demonstrate a counterintuitive result that increasing the ribosomal footprint may sometimes lead to an increase in the network production rate. Next, we propose a conservative deep neural network (CDNN) framework to approximate the solution of the network. The proposed loss function also incorporates the term satisfying the first integral property of the network. Point-wise comparison of the solutions by CDNN is in good agreement with the Runge–Kutta based numerical solution. Also, the CDNN framework offers a closed-form solution of the RFMEONP as a function of free parameters, thus allowing evaluation of the solution at any parameter value without again simulating the system.

A Second-Order Numerical Method for a Class of Optimal Control Problems

The numerical solution of optimal control problems through second-order methods is examined in this paper. Controlled processes are described by a system of nonlinear ordinary differential equations. There are two specific characteristics of the class of control actions used. The first one is that controls are searched for in a given class of functions, which depend on unknown parameters to be found by minimizing an objective functional. The parameter values, in general, may be different at different time intervals. The second feature of the considered problem is that the boundaries of time intervals are also optimized with fixed values of the parameters of the control actions in each of the intervals. The special cases of the problem under study are relay control problems with optimized switching moments. In this work, formulas for the gradient and the Hessian matrix of the objective functional with respect to the optimized parameters are obtained. For this, the technique of fast differentiation is used. A comparison of numerical experiment results obtained with the use of first- and second-order optimization methods is presented.