Search results

The purpose of this study is to solve two unique but difficult partial differential equations: the foam drainage equation and the nonlinear time-fractional fisher’s equation. Through our methods, we aim to provide accurate solutions and gain a deeper understanding of the intricate behaviors exhibited by these systems.

In this study, we use a dual technique that combines the Aboodh residual power series method and the Aboodh transform iteration method, both of which are combined with the Caputo operator.

We develop exact and efficient solutions by merging these unique methodologies. Our results, presented through illustrative figures and data, demonstrate the efficacy and versatility of the Aboodh methods in tackling such complex mathematical models.

Owing to their fractional derivatives and nonlinear behavior, these equations are crucial in modeling complex processes and confront analytical complications in various scientific and engineering contexts.

This study aims to investigate two newly developed (3 + 1)-dimensional Kairat-II and Kairat-X equations that illustrate relations with the differential geometry of curves and equivalence aspects.

The Painlevé analysis confirms the complete integrability of both Kairat-II and Kairat-X equations.

This study explores multiple soliton solutions for the two examined models. Moreover, the author showed that only Kairat-X give lump solutions and breather wave solutions.

The Hirota’s bilinear algorithm is used to furnish a variety of solitonic solutions with useful physical structures.

This study also furnishes a variety of numerous periodic solutions, kink solutions and singular solutions for Kairat-II equation. In addition, lump solutions and breather wave solutions were achieved from Kairat-X model.

The work formally furnishes algorithms for studying newly constructed systems that examine plasma physics, optical communications, oceans and seas and the differential geometry of curves, among others.

This paper presents an original work that presents two newly developed Painlev\'{e} integrable models with insightful findings.

The purpose of this paper is to investigate a variety of Painlevé integrable equations derived from a Hamiltonian equation.

The newly developed Painlevé integrable equations have been handled by using Hirota’s direct method. The authors obtain multiple soliton solutions and other kinds of solutions for these six models.

The developed Hamiltonian models exhibit complete integrability in analogy with the original equation.

The present study is to address these two main motivations: the study of the integrability features and solitons and other useful solutions for the developed equations.

The work introduces six Painlevé-integrable equations developed from a Hamiltonian model.

The work presents useful algorithms for constructing new integrable equations and for handling these equations.

The paper presents an original work with newly developed integrable equations and shows useful findings.

This study aims to discuss the numerical solutions of weakly singular Volterra and Fredholm integral equations, which are used to model the problems like heat conduction in engineering and the electrostatic potential theory, using the modified Lagrange polynomial interpolation technique combined with the biconjugate gradient stabilized method (BiCGSTAB). The framework for the existence of the unique solutions of the integral equations is provided in the context of the Banach contraction principle and Bielecki norm.

The authors have applied the modified Lagrange polynomial method to approximate the numerical solutions of the second kind of weakly singular Volterra and Fredholm integral equations.

Approaching the interpolation of the unknown function using the aforementioned method generates an algebraic system of equations that is solved by an appropriate classical technique. Furthermore, some theorems concerning the convergence of the method and error estimation are proved. Some numerical examples are provided which attest to the application, effectiveness and reliability of the method. Compared to the Fredholm integral equations of weakly singular type, the current technique works better for the Volterra integral equations of weakly singular type. Furthermore, illustrative examples and comparisons are provided to show the approach’s validity and practicality, which demonstrates that the present method works well in contrast to the referenced method. The computations were performed by MATLAB software.

The convergence of these methods is dependent on the smoothness of the solution, it is challenging to find the solution and approximate it computationally in various applications modelled by integral equations of non-smooth kernels. Traditional analytical techniques, such as projection methods, do not work well in these cases since the produced linear system is unconditioned and hard to address. Also, proving the convergence and estimating error might be difficult. They are frequently also expensive to implement.

There is a great need for fast, user-friendly numerical techniques for these types of equations. In addition, polynomials are the most frequently used mathematical tools because of their ease of expression, quick computation on modern computers and simple to define. As a result, they made substantial contributions for many years to the theories and analysis like approximation and numerical, respectively.

This work presents a useful method for handling weakly singular integral equations without involving any process of change of variables to eliminate the singularity of the solution.

To the best of the authors’ knowledge, the authors claim the originality and effectiveness of their work, highlighting its successful application in addressing weakly singular Volterra and Fredholm integral equations for the first time. Importantly, the approach acknowledges and preserves the possible singularity of the solution, a novel aspect yet to be explored by researchers in the field.

The main objective of this paper is to investigate the response of human skin to an intense temperature drop at the surface. In addition, this paper aims to evaluate the efficiency of finite difference and finite volume methods in solving the highly nonlinear form of Pennes’ bioheat equation.

One-dimensional linear and nonlinear forms of Pennes’ bioheat equation with uniform grids were used to study the behavior of human skin. The specific heat capacity, thermal conductivity and blood perfusion rate were assumed to be linear functions of temperature. The nonlinear form of the bioheat equation was solved using the Newton linearization method for the finite difference method and the Picard linearization method for the finite volume method. The algorithms were validated by comparing the results from both methods.

The study demonstrated the capacity of both finite difference and finite volume methods to solve the one-dimensional and highly nonlinear form of the bioheat equation. The investigation of human skin’s thermal behavior indicated that thermal conductivity and blood perfusion rate are the most effective properties in mitigating a surface temperature drop, while specific heat capacity has a lesser impact and can be considered constant.

This paper modeled the transient heat distribution within human skin in a one-dimensional manner, using temperate-dependent physical properties. The nonlinear equation was solved with two numerical methods to ensure the validity of the results, despite the complexity of the formulation. The findings of this study can help in understanding the behavior of human skin under extreme temperature conditions, which can be beneficial in various fields, including medical and engineering.

In this paper, the authors study the nonlinear matrix equation Xp=Q±A(X-1+B)-1AT, that occurs in many applications such as in filtering, network systems, optimal control and control theory.

The authors present some theoretical results for the existence of the solution of this nonlinear matrix equation. Then the authors propose two iterative schemes without inversion to find the solution to the nonlinear matrix equation based on Newton's method and fixed-point iteration. Also the authors show that the proposed iterative schemes converge to the solution of the nonlinear matrix equation, under situations.

The efficiency indices of the proposed schemes are presented, and since the initial guesses of the proposed iterative schemes have a high cost, the authors reduce their cost by changing them. Therefore, compared to the previous scheme, the proposed schemes have superior efficiency indices.

Finally, the accuracy and effectiveness of the proposed schemes in comparison to an existing scheme are demonstrated by various numerical examples. Moreover, as an application, by using the proposed schemes, the authors can get the optimal controller state feedback of $x(t+1) = A x(t) + C v(t)$.

The purpose of the present article is to obtain the similarity solution for the shock wave generated by a piston propagating in a self-gravitating nonideal gas under the impact of azimuthal magnetic field for adiabatic and isothermal flows.

The Lie group theoretic method given by Sophus Lie is used to obtain the similarity solution in the present article.

Similarity solution with exponential law shock path is obtained for both ideal and nonideal gas cases. The effects on the flow variables, density ratio at the shock front and shock strength by the variation of the shock Cowling number, adiabatic index of the gas, gravitational parameter and nonidealness parameter are investigated. The shock strength decreases with an increase in the shock Cowling number, nonidealness parameter and adiabatic index, whereas the strength of the shock wave increases with an increase in gravitational parameter.

Propagation of shock wave with spherical geometry in a self-gravitating nonideal gas under the impact of azimuthal magnetic field for adiabatic and isothermal flows has not been studied by any author using the Lie group theoretic method.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

Two-dimensional (2D) problems are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are considered. The problems are transformed into a boundary-only integral equation which can be solved numerically using a standard boundary element method (BEM). Some examples are solved to show the validity of the analysis and examine the accuracy of the numerical method.

The 2D problems which are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are solved using a combined BEM and Laplace transform. The time–space dependent coefficient equation is reduced to a time-dependent coefficient equation using an analytical transformation. Then, the time-dependent coefficient equation is Laplace transformed to get a constant coefficient equation, which can be written as a boundary-only integral equation. By utilizing a BEM, this integral equation is solved to find numerical solutions to the problems in the frame of the Laplace transform. These solutions are then inversely transformed numerically to obtain solutions in the original time–space frame.

The main finding of this research is the derivation of a boundary-only integral equation for the solutions of initial-boundary value problems governed by a modified-Helmholtz equation of time–space dependent coefficients for anisotropic functionally graded materials with time-dependent properties.

The originality of the research lies on the time dependency of properties of the functionally graded material under consideration.