Order Differential Equations Research Articles

Nonlinear vibro-acoustic analysis of an axially moving submerged circular cylindrical shell

In this study, the nonlinear vibro-acoustic dynamics and stability of doubly-clamped axially moving cylindrical shell are investigated. The exterior surface of the shell is in contact with fluid and subjected to oblique incident plane sound wave. Donnell’s nonlinear shallow shell theory is used to derive the nonlinear partial differential equation of the cylindrical shell for the radial motion. The Galerkin method is employed to discretize the equations of motion into the set of coupled nonlinear nonhomogeneous ordinary second order differential equations. Considering both driven and companion modes, Multiple Scales Method is used to obtain the response of the system. The effects of sound level, incident angle and axially velocity on the frequency response of the system are studied. The results show that, with increase in the shell velocity transmission loss of the circular shell increases. In addition, at low shell axial velocities simultaneous excitation of the driven and companion modes occurs.

The dynamics and behavior of logarithmic type fuzzy difference equation of order two.

Fuzzy difference equations are becoming increasingly popular in fields like engineering, ecology, and social science. Difference equations find numerous applications in real-life problems. Our study demonstrates that the logarithmic-type fuzzy difference equation of order two possesses a nonnegative solution, and an equilibrium point, and exhibits asymptotic behavior. [Formula: see text] Where, xi represents the sequence of fuzzy numbers, and the parameters α, β, A, along with the initial conditions x-1 and x0, are positive fuzzy numbers. The characterization theorem is employed to convert each single logarithmic fuzzy difference equation into a set of two crisp logarithmic difference equations within a fuzzy environment. We evaluated the stability of the equilibrium point of the fuzzy system. Utilizing variational iteration techniques, the method of g-division, inequality skills, and a theory of comparison for logarithmic fuzzy difference equations, we investigated the governing equation dynamics, including its boundedness, existence, and both local and global stability analysis. Additionally, we provided some numerical solutions for the equation describing the system to verify our results.

On a Poincaré-Perron problem for high order differential equations

On a Poincaré-Perron problem for high order differential equations

Applications of Laguerre transform to solve Schrödinger-type equations and differential equations of order four

The finite Laguerre transform is applied to solve differential equations problems of order higher than two and a one-dimensional steady-state Schrödinger equation, by using elementary Linear Algebra methods.

Entropy generation of Williamson hydromagnetic non-linear radiative nanofluid flow near a stagnation point over a porous stretching sheet

Purpose This study aims to analyze the multi-slip effects of entropy generation in steady non-linear magnetohydrodynamics thermal radiation with Williamson nanofluid flow across a porous stretched sheet near a stagnation point. Also, the qualities of viscous dissipation, Cattaneo–Christove heat flux and Arrhenius activation energy are taken into account. Thermophoresis, Brownian motion and Joule heating are also considered. Design/methodology/approach The Navier–Stokes equation, the thermal energy equation and the Solutal concentration equations are the governing mathematical equations that describe the flow and heat and mass transfer phenomena for fluid domains. By using the proper similarity transformations, a set of ordinary differential equationss are retrieved from boundary flow equations. The classical Runge–Kutta fifth-order algorithm along with the shooting technique is implemented to solve the obtained first order differential equations. Findings The study concludes that the temperature distribution boosting for thermal radiation, magnetic field and Eckert number where as the velocity and entropy generation escalate for the Williamson parameter, diffusion parameter and Brinkman number. The skin-friction and heat and mass transfer rate increases with the fluid injection. In addition, tabulated values of friction drag and rate of heat and mass transfer for various values of constraints are provided. Originality/value The comparison of the present results is carried out with the published results and noted a good agreement.

Evaluation of thermal and concentration slip effects on heat and mass transmission of nanofluid over a moving wedge surface using Keller box scheme

The current research is based on the impact of thermal and solutal slip in the boundary layer nanofluid flow through a moving accelerating wedge. The present investigation is considered with the influence of Brownian motion and thermophoresis. Thermal insulation, geothermal engineering, crude oil extraction, and heat exchangers are very important applications of nanofluid movement over a wedge surface with thermal and concentration slip. The suggested mathematical analysis is expressed in terms of partial differential equations (PDEs). These PDEs are transformed into ordinary differential equations via similarity transformation. The Keller Box technique is used to integrate the resultant non-similar equations. The set of discretized and first order differential equations is formed with the help of central difference and the Newton–Raphson technique. The graphical and numerical results are extracted with the help of MATLAB. The numerical results with the influence of the Prandtl factor (Pr), constant moving factor (λ), thermal slip factor (S2), and concentration slip parameter (S2) are interpreted visually and numerically. Graphical representations of velocity, thermal, and mass concentration profiles are analyzed in depth. The solution for skin friction coefficient, heat transport rate, and mass transport rate is calculated. The moving velocity function increases as Pr increases. The rate of slip temperature and slip concentration rate is enhanced for a lower Prandtl factor. The maximum slip behavior in temperature function and fluid concentration slip is deduced for each value of thermal-slip and concentration-slip factors. For high Prandtl and Brownian motion factors, the rate of Nusselt number is enhanced significantly.

Curvature form of Raychaudhuri equation and its consequences: A geometric approach

The paper aims at deriving a curvature form of the famous Raychaudhuri equation (RE) and the associated criteria for focusing of a hyper-surface orthogonal congruence of time-like geodesic. Moreover, the paper identifies a transformation of variable related to the metric scalar of the hyper-surface which converts the first order RE into a second order differential equation that resembles an equation of a Harmonic oscillator and also gives a first integral that yields the analytic solution of the RE and Lagrangian of the dynamical system representing the congruence.

Three-dimensional elasticity solutions for doubly-curved composite shells by extended differential quadrature method

This paper presents the three-dimensional (3D) stress analysis of doubly-curved composite shells with general boundary conditions using the strong sampling surfaces (SaS) formulation. The SaS method is based on the choice of SaS parallel to the middle surface and located at Chebyshev polynomial nodes to introduce the displacements of these surfaces as unknown functions that leads to a non-conventional shell formulation, in which strain–displacement and stress–strain relationships are represented in terms of SaS displacements. This is due to the use of Lagrange polynomials in approximating displacements, strains and stresses in the thickness direction. The outer surfaces are not included into a set of SaS that makes it possible to uniformly minimize the error ca used by Lagrange interpolation. The strong SaS formulation, based on direct integration of elasticity equilibrium equations in the thickness direction by the extended differential quadrature (EDQ) method, can be applied efficiently for high-precision calculations of doubly-curved composite shells with clamped and free edges. This is because in the SaS/EDQ formulation, displacements, strains and stresses of SaS are interpolated in a rectangular domain specified in a curvilinear coordinate system using a Chebyshev-Gauss-Lobatto grid and Lagrange polynomials are also used as basis functions. The proposed approach deals with equilibrium equations in terms of SaS stresses, avoiding the integration of second order differential equations in terms of SaS displacements, that greatly simplifies the implementation of the EDQ method.

A New Methodology for the Development of Efficient Multistep Methods for First–Order IVPs with Oscillating Solutions IV: The Case of the Backward Differentiation Formulae

A theory for the calculation of the phase–lag and amplification–factor for explicit and implicit multistep techniques for first–order differential equations was recently established by the author. His presentation also covered how the approaches’ efficacy is affected by the elimination of the phase–lag and amplification–factor derivatives. This paper will apply the theory for computing the phase–lag and amplification–factor, originally developed for implicit multistep methods, to a subset of implicit methods, called backward differentiation formulae (BDF), and will examine the impact of the phase–lag and amplification–factor derivatives on the efficiency of these strategies. Next, we will show you the stability zones of these brand-new approaches. Lastly, we will discuss the results of numerical experiments and draw some conclusions about the established approaches.

Design of integrated evolutionary finite differences for nonlinear electrohydrodynamics ion drag flow in cylindrical conduit model

This research implements an evolutionary optimized finite differences scheme (FDS) for nonlinear electrohydrodynamics ion drag flow dynamics in a cylindrical conduit (EHD-IDFCC). In the scheme, finite differences are exploited to discretize governing expressions of EHD-IDFCC, represented with a nonlinear singular seconder order differential equation, into a nonlinear equations-based system. The fitness function based on residual error is constructed for the FDS-based discretized EHD-IDFCC model. Its optimization is conducted with the global search capability of genetic algorithms (GAs) aided by the local search efficacy of the interior-point method (IPM), i.e., FDS-GA-IPM. The performance of the designed stochastic numerical solver FDS-GA-IPM is evaluated for a variant of the EHD-IDFCC model for different scenarios to measure the effect of axial flow velocity by varying the electric Hartmann numbers, as well as the nonlinearity factor. Statistical interpretations based on histogram illustrations, probability plots, and boxplots in terms of the cost function, mean absolute error (MAE), Thai inequality coefficients (TIC), and coefficient of determination metrics (R2) are used to endorse the accuracy, convergence, stability, and strength of the FDS-GA-IPM.

Geometrical description and Faddeev-Jackiw quantization of electrical networks

In lumped-element electrical circuit theory, the problem of solving Maxwell's equations in the presence of media is reduced to two sets of equations, the constitutive equations encapsulating local geometry and dynamics of a confined energy density, and the Kirchhoff equations enforcing conservation of charge and energy in a larger, topological, scale. We develop a new geometric and systematic description of the dynamics of general lumped-element electrical circuits as first order differential equations, derivable from a Lagrangian and a Rayleigh dissipation function. Through the Faddeev-Jackiw method we identify and classify the singularities that arise in the search for Hamiltonian descriptions of general networks. The core of our solution relies on the correct identification of the reduced manifold in which the circuit state is expressible, e.g., a mix of flux and charge degrees of freedom, including the presence of compact ones. We apply our fully programmable method to obtain (canonically quantizable) Hamiltonian descriptions of nonlinear and nonreciprocal circuits which would be cumbersome/singular if pure node-flux or loop-charge variables were used as a starting configuration space. We also propose a specific assignment of topology for the branch variables of energetic elements, that when used as input to the procedure gives results consistent with classical descriptions as well as with spectra of more involved quantum circuits. This work unifies diverse existent geometrical pictures of electrical network theory, and will prove useful, for instance, to automatize the computation of exact Hamiltonian descriptions of superconducting quantum chips.

Influence of pressure anisotropy on mass-radius relation and stability of millisecond pulsars in f(Q) gravity

In this study we explore the astrophysical implications of pressure anisotropy on the physical characteristics of millisecond pulsars within the framework of f(Q) gravity, in particular f(Q) = - α Q - β,where α and β are constants. Starting off with the field equations for anisotropic matter configurations, we adopt the physically salient Durgapal-Fuloria ansatz together with a well-motivated anisotropic factor for the interior matter distribution. This leads to a nonlinear second order differential equation which is integrated to give the complete gravitational and thermodynamical properties of the stellar object. The resulting model is subjected to rigorous tests to ensure that it qualifies as a physically viable compact object within the f(Q)-gravity framework. We study in detail the impact of anisotropy on the mass, radius and stability of the star. Our analyses indicate that our models are well-behaved, singularity-free and can account for the existence of a wide range of observed pulsars with masses ranging from 2.08 to 2.67 M ⊙, with the upper value being in the so-called mass gap regime observed in gravitational events such as GW190814. A comparison of the so-called Symmetric Teleparallel Equivalent to GR (STEGR) models with classical General Relativity (GR) models reveal that the anisotropy parameter and the sign of β impact on the predicted radii of pulsars. In particular, STEGR models have larger radii than their GR counterparts.

Dynamics of quasi-periodic, bifurcation, sensitivity and three-wave solutions for (n + 1)-dimensional generalized Kadomtsev-Petviashvili equation.

This study endeavors to examine the dynamics of the generalized Kadomtsev-Petviashvili (gKP) equation in (n + 1) dimensions. Based on the comprehensive three-wave methodology and the Hirota's bilinear technique, the gKP equation is meticulously examined. By means of symbolic computation, a number of three-wave solutions are derived. Applying the Lie symmetry approach to the governing equation enables the determination of symmetry reduction, which aids in the reduction of the dimensionality of the said equation. Using symmetry reduction, we obtain the second order differential equation. By means of applying symmetry reduction, the second order differential equation is derived. The second order differential equation undergoes Galilean transformation to obtain a system of first order differential equations. The present study presents an analysis of bifurcation and sensitivity for a given dynamical system. Additionally, when an external force impacts the underlying dynamic system, its behavior resembles quasi-periodic phenomena. The presence of quasi-periodic patterns are identified using chaos detecting tools. These findings represent a novel contribution to the studied equation and significantly advance our understanding of dynamics in nonlinear wave models.

Shifted Legendre Gauss‐Lobatto collocation scheme for solving nonlinear coupled space‐time fractional reaction‐advection diffusion equations

AbstractThis work aims to develop and apply an algorithm for solving space‐time fractional reaction‐advection diffusion (STFRAD) nonlinear coupled equations with specified initial and boundary conditions by employing a Shifted Legendre Gauss‐Lobatto collocation (SLGLC) scheme. The fractional derivatives of time and space are defined in Caputo sense. The approximate solutions have been constructed with the help of shifted Legendre polynomials. The proposed numerical scheme reduces the coupled fractional order differential equations into a system of algebraic equations which can then be easily solved. The error analysis through the two test problems is carried out to validate the efficiency and efficacy of the method. The effect of advection and reaction terms, and various space‐time fractional order derivatives on the solution profiles have been studied and are shown graphically for different particular cases. It has been observed that the species whose transport is simulated using STFRADE shows higher concentration at a specific distance in porous media in case of integer order differential equations as compared to fractional order equations demonstrating overestimation of its impact in case of integer order differential equation. This is extremely important observation with respect to applications in the field of groundwater pollution management. The salient feature of the article is the stability and convergence analyses of the proposed scheme on the concerned model, which demonstrates the effectiveness and capabilities of the developed method.

Viscosity solutions of obstacle equations by inhomogeneous convex envelopes

We characterize the convex envelope of a continuous function as a viscosity solution of a second order differential equation of obstacle style. The notion of convexity that we investigate is an inhomogeneous generalization of the classical concept of convexity motivated by problems of singular control in finance.

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

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

A novel method to solve analytically the non-linear Poisson equation in the inversion layer of a MOSFET

Abstract Despite the enormous importance that the metal oxide semiconductors (MOS) and the field effect transistors (MOSFETs) have in the actual semiconductor technology, the task of finding the analytical solution of the Poisson equation in the inversion layer was not fully accomplished for more than half a century. It is a well-known fact that due to the non-linear nature of this equation, the attempts to solve it got stuck after the first integration. Nevertheless, experimental and applied researchers found ways to characterize, control, and develop MOSFET devices based on approximate models, and numerical calculations. Here I present a new method to analytically solve the nonlinear Poisson equation, in principle, for any charge distribution in the inversion layer of a MOS, and for the charge density of the original Shockley model, in particular. To this purpose, a physical argument related to the charge and field energies in the transient population-inversion process is introduced and a new nonlinear but solvable second order differential equation is obtained, whose solution also solves the original Poisson equation. The analytical results presented here allow us to derive explicit expressions for the electrical potential distribution and, very importantly, for the charge distribution, the inversion layer width and the effective impurity concentration in the depletion layer. The quantum Hall effect and other physical phenomena were discovered in the inversion layer of a MOS. The Hall effect was explained under the assumption that the charge distribution is a two-dimensional electron gas, we show here however that the 2D limit is reached only at high gate voltages. When applied to MOSFET structures, we obtain new expressions for the drain currents in the inversion and depletion layers, as functions of the impurity concentration and the gate voltage for single and double-gate MOSFETs, and we give an insight for multi-gate MOSFETs.

Polynomial solutions to the first order difference equations in the bivariate difference field

The bivariate difference field provides an algebraic framework to study sequences satisfying a recurrence of order two, and it can be used to transform summations involving such sequences into first order difference equations over the bivariate difference field. Based on this, we present an algorithm for finding a series of polynomial solutions of such equations in the bivariate difference field, and show an upper bound on the degree of any possible polynomial solutions, which in turn is sufficient to compute all polynomial solutions by using the method of undetermined coefficients.

Stability of equilibrium states for a stochastic difference equation of exponential form

It is well know that many process and problems in different fields of sciences and engineering can be modelled by linear and nonlinear difference equations. Particularly, in case of systems biology, ecology, biochemistry, genetics and physiology dynamics, many population models are governed by exponential difference equations, and lot of papers were published in this matter. The goal of this work is to study a nonlinear second order difference equation of exponential form: Xn+1 = bXn + cXne−σXn−1, n ≥ 0 where the parameters b, c have arbitrary values, σ > 0 and the initial values, X0, X−1 are arbitrary positive numbers. Our equation has an equilibrium points and is exposed to additive stochastic perturbations that are assumed a sequence of independent random variables with zero mean, unit variance and these perturbations are proportional to the deviation of the system state Xn from equilibrium points, and the result is stochastic difference equation. We cannot find the solutions of nonlinear stochastic difference equations of exponential form in all cases. So, one can study the behavior of solutions by asymptotic stability of equilibrium points. Additionally, with the general method of Lyapunov functional constructions for stochastic difference equations with discrete time, we provide the necessary and sufficient conditions for the asymptotic mean square stability of the two equilibrium points (zero and positive) within the linear approximation of the stochastic difference equation. These conditions also serve as sufficient conditions for the stability in probability of the equilibrium points of the initial nonlinear equation. Finally, to demonstrate the validity of the obtained results, some numerical examples and simulations of equations solutions are made and numerous graphical illustrations of stability trajectories of solutions are plotted.

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

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