Solutions Of Differential Equations Research Articles

Reduced Desensitization Formulation for Optimal Control Problems

Solutions to nonlinear optimal control problems (OCPs) exhibit dispersions under model uncertainties and it is desirable to generate optimal solutions that exhibit less sensitivity to model uncertainties. We propose a novel solution desensitization method dubbed “Reduced Desensitization Formulation (RDF)” by leveraging non-uniqueness of the solution of the costate differential equations when a hybrid indirect-direct optimization method is used. A key property of the RDF method is a significant reduction in the number of differential equations needed for generating desensitized solutions. This feature facilitates the generation of open-loop desensitized trajectories and makes the methodology applicable to OCPs with a larger number of uncertain parameters. To demonstrate the utility of the RDF method, three important classes of trajectory optimization problems are considered with uncertainty in the thrust magnitude of the propulsion system: (1) minimum-fuel low-thrust interplanetary rendezvous maneuvers, (2) low-thrust orbit-raising maneuvers, and (3) minimum-fuel high-thrust rocket-landing problems. For the considered problems with bang-bang control profiles, an analysis is presented on the change in the number of control switches between sensitive and desensitized optimal solutions. Numerical results demonstrate desensitization of the considered performance indices with respect to the thrust magnitude of the propulsion system.

Thermal analysis of a viscoelastic Maxwell hybrid nanofluid with graphene and polythiophene nanoparticles: Insights from an artificial neural network model

The utilization of solar radiation by converting them into thermal energy is discussed in this paper. Nanoparticles improve the ability of heat transfer therefore, it is beneficial in the use of solar thermal systems and energy storage devices. The novel mixture of nanoparticles Graphene and Polythiophene in base fluid, which has high thermodynamic properties for the improvement of thermal effect with electromagnetic effect by using Maxwell fluid model is discussed. Polyvinyl alcohol water is taken as base fluid flowing through a moveable flat plat. The governing partial differential equations are transformed into ordinary differential equations. The semi-analytical technique, homotopy analysis method is used to obtain the solution of the ordinary differential equations. The velocity is enhanced with magnetic and electric field strength. The increase of the Prandtl number, Eckert number and chemical reaction parameter, exceeds the thermal effect which produces more entropy generation and heat enhancement. The results show that the hybrid nanofluid with this Novel mixture is highly thermodynamic with higher entropy and rapid thermal augmentation which can be used in energy production and energy storage devices. A novel intelligent numerical computing technique multi-layer perceptron with feed-forward back-propagation, an artificial neural networking method with the Levenberg-Marquard algorithm is used in this model. The data is gathered for the neural networking method training, validation, and testing. The efficiency of the model is obtained and mean square error is obtained by artificial neural networking.

A Chebyshev polynomial based block integrator for the direct numerical solution of fourth order ordinary differential equations

This paper introduces an innovative method for numerically integrating fourth-order initial value problems by utilizing Chebyshev polynomials as the fundamental basis function. The block integrator base on Chebyshev polynomial demonstrates significant improvements in accuracy and stability, rendering it a valuable tool across various scientific and engineering fields. By leveraging the characteristics of Chebyshev polynomials, this approach accurately estimates solutions for fourth-order differential equations without reducing it to a system of first order Ordinary Differential Equations while at the same time effectively managing error accumulation within a block integration framework and thereby enhancing its accuracy over extended intervals. Through rigorous numerical experiments, the effectiveness and reliability of the new integrator are demonstrated and compared with existing methods. The new method is consistent, zero stable and convergent. The method also shows an appreciable error constants. The new method performed better in terms of accuracy than the existing methods in the literature in both linear and nonlinear problems.

Explicit optical solitons of a perturbed Biswas–Milovic equation having parabolic-law nonlinearity and spatio-temporal dispersion

This paper deals with a new variant of the Biswas–Milovic equation, referred to as the perturbed Biswas–Milovic equation with parabolic-law nonlinearity in spatio-temporal dispersion. To our best knowledge, the considered equation which models the pulse propagation in optical fiber is studied for the first time, and the abundant optical solitons are successfully obtained utilizing the auxiliary equation method. Utilizing a wave transformation technique on the considered Biswas–Milovic equation, and by distinguishing its real and imaginary components, we have been able to restructure the considered equation into a set of nonlinear ordinary differential equations. The solutions for these ordinary differential equations, suggested by the auxiliary equation method, include certain undetermined parameters. These solutions are then incorporated into the nonlinear ordinary differential equation, leading to the formation of an algebraic equation system by collecting like terms of the unknown function and setting their coefficients to zero. The undetermined parameters, and consequently the solutions to the Biswas–Milovic equation, are derived by resolving this system. 3D, 2D, and contour graphs of the solution functions are plotted and interpreted to understand the physical behavior of the model. Furthermore, we also investigate the impact of the parameters such as the spatio-temporal dispersion and the parabolic nonlinearity on the behavior of the soliton. The new model and findings may contribute to the understanding and characterization of the nonlinear behavior of pulse propagation in optical fibers, which is crucial for the development of optical communication systems.

An Interpretable Approach to the Solutions of High-Dimensional Partial Differential Equations

In recent years, machine learning algorithms, especially deep learning, have shown promising prospects in solving Partial Differential Equations (PDEs). However, as the dimension increases, the relationship and interaction between variables become more complex, and existing methods are difficult to provide fast and interpretable solutions for high-dimensional PDEs. To address this issue, we propose a genetic programming symbolic regression algorithm based on transfer learning and automatic differentiation to solve PDEs. This method uses genetic programming to search for a mathematically understandable expression and combines automatic differentiation to determine whether the search result satisfies the PDE and boundary conditions to be solved. To overcome the problem of slow solution speed caused by large search space, we propose a transfer learning mechanism that transfers the structure of one-dimensional PDE analytical solution to the form of high-dimensional PDE solution. We tested three representative types of PDEs, and the results showed that our proposed method can obtain reliable and human-understandable real solutions or algebraic equivalent solutions of PDEs, and the convergence speed is better than the compared methods. Code of this project is at https://github.com/grassdeerdeer/HD-TLGP.

Oscillator Simulation with Deep Neural Networks

The motivation behind this study is to overcome the complex mathematical formulation and time-consuming nature of traditional numerical methods used in solving differential equations. It seeks an alternative approach for more efficient and simplified solutions. A Deep Neural Network (DNN) is utilized to understand the intricate correlations between the oscillator’s variables and to precisely capture their dynamics by being trained on a dataset of known oscillator behaviors. In this work, we discuss the main challenge of predicting the behavior of oscillators without depending on complex strategies or time-consuming simulations. The present work proposes a favorable modified form of neural structure to improve the strategy for simulating linear and nonlinear harmonic oscillators from mechanical systems by formulating an ANN as a DNN via an appropriate oscillating activation function. The proposed methodology provides the solutions of linear and nonlinear differential equations (DEs) in differentiable form and is a more accurate approximation as compared to the traditional numerical method. The Van der Pol equation with parametric damping and the Mathieu equation are adopted as illustrations. Experimental analysis shows that our proposed scheme outperforms other numerical methods in terms of accuracy and computational cost. We provide a comparative analysis of the outcomes obtained through our proposed approach and those derived from the LSODA algorithm, utilizing numerical techniques, Adams–Bashforth, and the Backward Differentiation Formula (BDF). The results of this research provide insightful information for engineering applications, facilitating improvements in energy efficiency, and scientific innovation.

Numerical solution of steady nonlinear differential equations for compressible flow through a spinning convergent divergent nozzle

Many projectiles tend to spin about their longitudinal axis while progressing in the forward direction. It helps in providing stability and a reference direction for guidance during their run. Many different projectiles employ a supersonic convergent-divergent nozzle to produce thrust for their forward motion; hence, the nozzle and overall whole propulsion system tend to spin about its axis of rotation. The main aim of this study is to observe the effect of spin on the nozzle. In this research, a converging bell-shaped diverging nozzle is numerically designed using a method of characteristics (MOC) for exit Mach number 3.21. Viscous simulations are performed for both two- and three-dimensional cases. The analysis is then performed with nozzle spinning about its axis of symmetry with a constant angular velocity of 10 revolutions per second. The analysis is repeated for the value of constant angular velocities to be 15 and 20 revolutions per second, and the behavior of flow with increasing angular velocity is examined. It has been observed that the exit Mach number and velocity decrease due to the radial protrusion of the boundary layer, and it has a negative impact on the performance of the nozzle. Moreover, the decrease of exit Mach number is in direct relation to increasing angular velocity.

Exponential stability results for second-order impulsive neutral stochastic differential equations

This paper is aimed to discuss the existence and exponential stability of mild solutions for second-order impulsive neutral stochastic differential equations (INSDEs). Firstly, we derive the existence of mild solutions for INSDEs by using the Banach fixed point theorem. Then, the exponential stability in the pth moment of mild solution for the considered INSDEs is proved by means of the generalized impulsive-integral inequality, which can effectively improve some known existing ones. Finally, as an application, an example is given to illustrate the efficiency of the obtained theoretical results.

On a computable Skorokhod's integral‐based estimator of the drift parameter in fractional SDE

AbstractThis paper deals with a Skorokhod's integral‐based least squares‐ (LS) type estimator of the drift parameter computed from multiple (possibly dependent) copies of the solution of a stochastic differential equation (SDE) driven by a fractional Brownian motion of Hurst index . On the one hand, some convergence results are established on our LS estimator when . On the other hand, when , Skorokhod's integral‐based estimators cannot be computed from data, but in this paper some convergence results are established on a computable approximation of our LS estimator.

The Lie point symmetry criteria and formation of exact analytical solutions for Kairat-II equation: Paul-Painlevé approach

The renowned non-linear Kairat-II equation, which expresses the surface curves, is examined meticulously in this work. There has never been a study before that discussed conserved quantities, sensitivity analysis, Hamiltonian function, Lie symmetry invariance criteria, and invariant solutions of Kairat-II equation The Lie invariance criteria are taken into consideration by the symmetry generators. The suggested technique leads to in a four-dimensional Lie algebra, where the translation point symmetries in space and time correlate with the conservation of mass and the energy, respectively, and remaining point symmetries are dilation and scaling. Firstly, we use symmetry reduction of Lie subalgebras to obtain closed-form invariant solutions. In specific reduction cases, we transform the Kairat-II equation into a spectrum of non-linear ordinary differential equations, which have the benefit of being able to offer an extensive selection of closed-form solitary wave solutions. For this equation, the Cauchy problem cannot be solved by the inverse scattering transform, therefore, traveling wave exact solutions are generated using the analytical Paul-Painlevé approach. In both two and three dimensions, the graphical behavior of specific solutions is presented for specific quantities of physical factors of examined equation. Utilizing conservation laws multipliers, the study culminates by determining an extensive set of local conservation laws for the non-linear Kairat-II equation that are applicable to arbitrary constant coefficients. Instead of focusing on the physical aspects of conservation laws, this study computed the conserved quantities using a mathematical perspective that can be used to identify potential symmetries of partial differential equations, guiding readers toward the solutions of partial differential equations. The existence of the Hamiltonian function is presented. The model’s sensitivity to various starting conditions is highlighted by the sensitivity analysis.

On the stability of periodic solutions of piecewise smooth periodic differential equations

On the stability of periodic solutions of piecewise smooth periodic differential equations

Plane Stress Problems for Isotropic Incompressible Hyperelastic Materials

AbstractThe analysis of plane stress problems has long been a topic of interest in linear elasticity. The corresponding problem for non-linearly elastic materials is considered here within the context of homogeneous incompressible isotropic elasticity. It is shown that when the problem is posed in terms of the Cauchy stress, a semi-inverse approach must be employed to obtain the displacement of a typical particle. If however the general plane stress problem is formulated in terms of the Piola-Kirchhoff stress, the deformation of a particle requires the solution of a non-linear partial differential equation for both simple tension and simple shear, the trivial solution of which yields a homogeneous deformation. It is also shown that the general plane stress problem can be solved for the special case of the neo-Hookean material.

Solving Nonlinear Time-Fractional Partial Differential Equations Using Conformable Fractional Reduced Differential Transform with Adomian Decomposition Method

In this article, we use a new technique called conformable fractional reduced differential transform (CFRDT) with Adomian decomposition to estimate the solution of one and two-dimensional time-fractional partial linear and nonlinear differential equations with initial values. We explain the convergence analysis of this technique. The obtained results illustrate that the novel method is efficient and easy to use to find approximate solutions for the time-fractional partial differential equations (PDEs). Thus, the suggested method has a significant impact on how engineering, physics, and other disciplines solve fractional PDEs. Furthermore, we analyze the solution of problems with a 2D or 3D graphical representation by using Mathematica software.

Elzaki residual power series method to solve fractional diffusion equation.

The time-fractional order differential equations are used in many different contexts to analyse the integrated scientific phenomenon. Hence these equations are the point of interest of the researchers. In this work, the diffusion equation for a one-dimensional time-fractional order is solved using a combination of residual power series method with Elzaki transforms. The residual power series approach is a useful technique for finding approximate analytical solutions of fractional differential equations that needs the residual function's (n-1)α derivative. Since it is challenging to determine a function's fractional-order derivative, the traditional residual power series method's application is somewhat constrained. The Elzaki transform with residual power series method is an attempt to get over the limitations of the residual power series method. The obtained numerical solutions are compared with the exact solution of this equation to discuss the method's applicability and efficiency. The results are also graphically displayed to show how the fractional derivative influences the behaviour of the solutions to the suggested method.

Solving Linear and Nonlinear Fuzzy Fractional Volterra-Fredholm Integro Differential Equations Using Shehu Adomian Decomposition Method

Objectives: In applied sciences and engineering, fuzzy fractional differential equations (FFDEs) and fuzzy fractional integral equations (FFIEs) are a crucial topic. The main objective of this work is to discover an analytical approximate solution for the fuzzy fractional Volterra-Fredholm integro differential equations (FFVFIDE). In the Caputo concept, fractional derivatives are regarded. Methods: The Shehu transform is challenging to exist for nonlinear problems. So, the Shehu transform is combined with the Adomian decomposition method is called the Shehu Adomian decomposition method (SHADM) and has been proposed to solve both linear and nonlinear FFVFIDEs. Findings: Both linear and nonlinear FFVIFIDEs can be solved using this technique. For nonlinear terms, Adomian polynomials have been used. The main benefit of this approach is that it converges quickly to the exact solution. Figures and numerical examples demonstrate the expertise of the suggested approach. Novelty: The comparison between the exact solution and numerical solution is shown in figures for various values of fractional order . The numerical evolution demonstrates the efficiency and reliability of the proposed SHADM. The proposed approach is rapid, exact, and simple to apply and produce excellent outcomes. Keywords: Fractional calculus, fuzzy number, Mittag ­ Leffler function, Shehu Adomian decomposition method, fuzzy fractional Volterra­-Fredholm integro differential equation

Correlated dynamics of immune network and sl(3, R) symmetry algebra

Abstract We observed the existence of periodic orbits in immune network under transitive solvable Lie algebra. In this article, we focus to develop condition of maximal Lie algebra for immune network model and use that condition to construct a vector field of symmetry to study nonlinear pathogen model. We used two methods to obtain analytical structure of solution, namely normal generator and differential invariant function. Numerical simulation of analytical structure exhibits correlated periodic pattern growth under spatiotemporal symmetry, which is similar to the linear dynamical simulation result. We used Lie algebraic method to understand correlation between growth pattern and symmetry of dynamical system. We employ idea of using one parameter point group of transformation of variables under which linear manifold is retained. In procedure, we present the method of deriving Lie point symmetries, the calculation of the first integral and the invariant solution for the ordinary differential equation (ODE). We show the connection between symmetries and differential invariant solutions of the ODE. The analytical structure of the solution exhibits periodic behavior around attractor in local domain, same behavior obtained through dynamical analysis.

On the Solution of Blasius Boundary Layer Equations of Prandtl-Eyring Fluid Flow Past a Stretching Sheet

Objective: This paper investigates velocity profile for two-dimensional, incompressible, laminar forced convection flow of the fluid model for Prandtl-Eyring fluid past a stretching sheet in the presence of fluid parameters. Methods: The governing partial differential equation for the flow was transformed into non-linear ordinary differential equation by using the deductive one parameter group theoretic method and numerical solution of non-linear ordinary differential equation (ODE) is solved by MATLAB bvp4c solver. Findings: The solution of velocity profile obtained as a function of parameter and . The effect of the fluid parameter was discussed graphically. Novelty: The main goal of this article is to analyze boundary layer flow of Prandtl-Eyring fluid over a stretching surface. The conservation equations of mass, momentum are converted into non-linear ordinary differential equations along with boundary conditions using deductive one parameter group theoretic method and solved by MATLAB ODE solver. Comparisons with previously published works are made, and results show a high level of agreement. This type of research is applicable to extrusion, paper production, fiber glass production, hot rolling, condensation process, crystal growing, polymer sheets etc. Keywords: Boundary layer, laminar flow, Deductive one parameter Group theoretic method, Absolute invariant, Stretching Sheet, Prandtl-Eyring fluid

Tensor networks for solving the time-independent Boltzmann neutron transport equation

Tensor network techniques, known for their low-rank approximation ability that breaks the curse of dimensionality, are emerging as a foundation of new mathematical methods for ultra-fast numerical solutions of high-dimensional Partial Differential Equations (PDEs). Here, we present a mixed Tensor Train (TT)/Quantized Tensor Train (QTT) approach for the numerical solution of time-independent Boltzmann Neutron Transport equations (BNTEs) in Cartesian geometry. Discretizing a realistic three-dimensional (3D) BNTE by (i) diamond differencing in space, (ii) multigroup-in-energy, and (iii) discrete ordinates collocation in angle leads to large generalized eigenvalue problems that generally require a matrix-free approach and large computer clusters. Starting from this discretization, we construct a TT representation of the PDE fields and discrete operators, followed by a QTT representation of the TT cores. We then solve the tensorized generalized eigenvalue problem using a fixed-point scheme with tensor network optimization techniques. We validate our approach by applying the method to two examples of 3D neutron transport problems, currently solved by the Los Alamos National Laboratory PARallel TIme-dependent SN (PARTISN) solver.1 We demonstrate that our TT/QTT method, executed on a standard desktop computer, leads to large compression. This allows for the storage of terrabyte-sized neutron angular flux eigenvectors in megabytes. Additionally, we create megabyte-sized full access TT representations of yottabyte-sized transport matrix operators. By leveraging the TT operators and solution methods, we obtain a 7500 times speedup when compared to the PARTISN solution time with an error of less than 10−5.

Existence of basic solutions of first order linear homogeneous set-valued differential equations

The paper presents various derivatives of set-valued mappings,their main properties and how they are related to each other.Next, we consider Cauchy problems with linear homogeneousset-valued differential equations with different types ofderivatives (Hukuhara derivative, PS-derivative andBG-derivative). It is known that such initial value problems withPS-derivative and BG-derivative have infinitely many solutions.Two of these solutions are called basic. These are solutions suchthat the diameter function of the solution section is amonotonically increasing (the first basic solution) or monotonicallydecreasing (the second basic solution) function. However, the secondbasic solution does not always exist. We provideconditions for the existence of basic solutions of such initialvalue problems. It is shown that their existence depends on thetype of derivative, the matrix of coefficients on the right-handand the type of the initial set. Model examples are considered.

Approximate solutions for fuzzy stochastic differential equations with Markovian switching

The fuzzy stochastic differential equations with Markovian switching are considered. First, under the Carath e ´ odory assumptions, the existence and uniqueness theorem for the aforementioned equations is given by means of stopping time techniques, discretization method, and the Gronwall-Bellman inequality. Then the boundedness of the second moment of the solutions for such equations is established. Furthermore, by using the averaging method and stochastic analysis, as well as the Bihari's inequality, this paper mainly demonstrates the asymptotic relationship between the solutions of the original equations and the corresponding averaged equations in terms of mean square and probability. Finally, the theoretical result of the averaging method is illustrated by an example.