Genetic diversity of <i>Salmonella enterica</i> subsp. <i>enterica</i> serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles

Alexey V. Rakov; Natalya A. Kuznetsova; Anatoly A. Yakovlev; Alexey V. Rakov; Natalya A. Kuznetsova; Anatoly A. Yakovlev

doi:10.3934/microbiol.2020007

AIMS Microbiology

2020, Volume 6, Issue 2: 106-120. doi: 10.3934/microbiol.2020007

Previous Article Next Article

Review

Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles

Laboratory of Molecular Epidemiology and Ecology of Pathogenic Bacteria, Somov Institute of Epidemiology and Microbiology, Vladivostok, Russia

Received: 18 March 2020 Accepted: 06 May 2020 Published: 08 May 2020

For the first time, in the literature review we presents the molecular genetic structure of Salmonella Enteritidis populations in Russia, and particularly, in Siberia and the Far East of the country. Pathogen population in Russia has been compared with Salmonella populations circulating in different countries of the world. It has been shown that the microbial population is heterogeneous, but it is possible to identify the dominant and main genotypes, which determine up to 90% of the total population morbidity. The data were obtained as a result of a 30-year monitoring (1988 to 2018) by studying the microbial plasmid profiles. It was shown that the same S. Enteritidis clones circulate throughout Russia, however, their significance in the population morbidity may vary depending on geographic and temporal characteristics. Population is characterized by heterogeneity and relative stability of the plasmid types' structure. At the same time, the population is also specified by variability, reflected as a simple change of the annual number of constantly detected plasmid types, and the appearance of new ones that can play a significant role in the etiology of Salmonella infection.

Keywords:

Citation: Alexey V. Rakov, Natalya A. Kuznetsova, Anatoly A. Yakovlev. Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles[J]. AIMS Microbiology, 2020, 6(2): 106-120. doi: 10.3934/microbiol.2020007

Related Papers:

[1]	Manoj Singh, Ahmed Hussein, Msmali, Mohammad Tamsir, Abdullah Ali H. Ahmadini . An analytical approach of multi-dimensional Navier-Stokes equation in the framework of natural transform. AIMS Mathematics, 2024, 9(4): 8776-8802. doi: 10.3934/math.2024426
[2]	Qasim Khan, Anthony Suen, Hassan Khan, Poom Kumam . Comparative analysis of fractional dynamical systems with various operators. AIMS Mathematics, 2023, 8(6): 13943-13983. doi: 10.3934/math.2023714
[3]	Krunal B. Kachhia, Jyotindra C. Prajapati . Generalized iterative method for the solution of linear and nonlinear fractional differential equations with composite fractional derivative operator. AIMS Mathematics, 2020, 5(4): 2888-2898. doi: 10.3934/math.2020186
[4]	Qingkun Xiao, Jianzhu Sun, Tong Tang . Uniform regularity of the isentropic Navier-Stokes-Maxwell system. AIMS Mathematics, 2022, 7(4): 6694-6701. doi: 10.3934/math.2022373
[5]	Muhammad Akram, Ghulam Muhammad, Tofigh Allahviranloo, Ghada Ali . A solving method for two-dimensional homogeneous system of fuzzy fractional differential equations. AIMS Mathematics, 2023, 8(1): 228-263. doi: 10.3934/math.2023011
[6]	Zeliha Korpinar, Mustafa Inc, Dumitru Baleanu . On the fractional model of Fokker-Planck equations with two different operator. AIMS Mathematics, 2020, 5(1): 236-248. doi: 10.3934/math.2020015
[7]	Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen . New solutions of time-space fractional coupled Schrödinger systems. AIMS Mathematics, 2023, 8(11): 27033-27051. doi: 10.3934/math.20231383
[8]	Yasir Nadeem Anjam . The qualitative analysis of solution of the Stokes and Navier-Stokes system in non-smooth domains with weighted Sobolev spaces. AIMS Mathematics, 2021, 6(6): 5647-5674. doi: 10.3934/math.2021334
[9]	Shabir Ahmad, Aman Ullah, Ali Akgül, Fahd Jarad . A hybrid analytical technique for solving nonlinear fractional order PDEs of power law kernel: Application to KdV and Fornberg-Witham equations. AIMS Mathematics, 2022, 7(5): 9389-9404. doi: 10.3934/math.2022521
[10]	Tareq Eriqat, Rania Saadeh, Ahmad El-Ajou, Ahmad Qazza, Moa'ath N. Oqielat, Ahmad Ghazal . A new analytical algorithm for uncertain fractional differential equations in the fuzzy conformable sense. AIMS Mathematics, 2024, 9(4): 9641-9681. doi: 10.3934/math.2024472

Abstract

1. Introduction

In the past decade, many researchers have studied many fractional partial differential equations (FPDEs) types. Since the beginning of fractional calculus history in 1695, when L'Hospital raised the question: what is the meaning of $\frac{{d}^{n}y}{{dx}^{n}}$ if $n = \frac{1}{2}$ ? That is, what if n is fractional? ^[1] Even for new researchers, the fractional derivatives were complicated-although it appears in many parts of sciences such as physics, engineering, bioengineering, COVID-19 studies, and many other branches of sciences ^{[2,3,4,5,6,7,8,9]}. In addition, many definitions of fractional derivatives have been given ^[10,11]. Fractional order derivatives of a given function involve the entire function history where the following state of a fractional order system is dependent on its current state and all its historical states ^{[4,5,6,7,8,9,10,11]}.

There are several analytical and numerical techniques for handling fractional problems, B-spline functions, Bernoulli polynomials, Adomian decomposition, variational iteration, Homotopy analysis, and many others ^{[12,13,14,15,16,17]}. On top of that, some applicable analytical methods are developed to address nonlinear problems with fractional derivatives. One of these approaches is the Laplace residual power series method (LRPSM) ^{[18,19,20,21,22,23,24,25,26,27,28,29]}, which shows its efficiency and applicability in solving nonlinear problems.

The LRPSM is a very modern technique, and it is a hybrid method of two approaches, the Laplace transform (LT) and the idea of the residual power series method (RPSM) ^{[30,31,32,33]}. In 2020, the authors in the article ^[18] were able to adapt the LT to solve nonlinear neutral fractional pantograph equations using the residual function and the RPSM idea. The LT, usually, is implemented to solve linear equations only, but the LRPSM can overcome this disadvantage and thus adapts it to solve nonlinear equations of different types. The LRPSM presents an approximate analytical solution with a series form using the concept of the Laurent series and the power series ^{[18,19,20,21,22,23,24,25]}. What distinguishes LRPSM from RPSM is the use of the idea of limit at infinity in getting the coefficients of a series solution rather than the concept of a fractional derivative as in RPSM. Many articles used the proposed method to treat several types of differential equations of fractional orders. In 2021, El-Ajou adapted LRPSM to establish solitary solutions of nonlinear dispersive FPDEs ^[19] and to present series solutions for systems of Caputo FPDEs with variable coefficients ^[20]. Newly, the LRPSM is used for solving Fuzzy Quadratic Riccati Differential Equations ^[21], time-fractional nonlinear water wave PDE ^[22], fractional Lane-Emden equations ^[23], Fisher's equation and logistic system model ^[24], and nonlinear fractional reaction-diffusion for bacteria growth model ^[25].

Claude Louis Navier and Gabriel Stokes have created the so-called Navier-Stokes equations (NSEs). A French mechanical engineer Claude was affiliated in continuum mechanics with a physicist specializing and the French government, whose main contribution was the Navier-Stokes equations (1822). This famous equation made his name among the several names incised on the Eiffel Tower. Moreover, Newton's second law for fluid substance, which is central to fluid mechanics has been used in describing many physical phenomena in many applied sciences ^{[34,35,36,37]}. For example, the study of airflow around a wing and water flow in pipes and used as one of the continuity equations needed to build microscopic models in 1985 and also as a special case considered to establish the relationship between external and pressure forces on the fluid to the responses of fluid flow ^[38].

The motivation of this work is to adapt the LRPSM to provide analytical solutions for a multi-dimensional time-fractional Navier-Stokes (M-DT-FNS) system which takes the following form ^[39]:

${D}_{t}^{\alpha }u+\left(u\cdot \nabla \right)u = \upsilon {\nabla }^{2}u-\frac{1}{\rho }\nabla p, \mathrm{ }\mathrm{ }\mathrm{ }\mathrm{ }0 < \alpha \le 1,$

(1.1)

where ${D}_{t}^{\alpha }$ is the Caputo fractional-derivative operator of order $\alpha$ , $p = p\left(\chi, \varsigma, \zeta, t\right)$ is the pressure, $\rho$ is the density, $u$ is a vector field that represents the flow velocity vector, $\upsilon = \frac{\mu }{\rho }$ is the kinematic viscosity ( $\mu$ is the dynamic viscosity), and $\nabla$ & ${\nabla }^{2}$ are the gradient and Laplacian operators, respectively, subject to the initial conditions (ICs) at the initial velocity:

$u = \varphi .$

(1.2)

If the density is constant throughout the fluid domain, then the vector Eq (1.1) is an incompressible NSEs.

The vector Eqs (1.1) and (1.2) can be separated in a system form as follows ^[39,40]:

$\begin{array}{l} {D}_{t}^{\alpha }{u}_{1}+{u}_{1}\frac{\partial {u}_{1}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{1}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{1}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial p}{\partial \chi }, \\ {D}_{t}^{\alpha }{u}_{2}+{u}_{1}\frac{\partial {u}_{2}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{2}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{2}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{2}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial p}{\partial \varsigma }, \\ {D}_{t}^{\alpha }{u}_{3}+{u}_{1}\frac{\partial {u}_{3}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{3}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{3}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{3}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial p}{\partial \zeta }, \end{array}$

(1.3)

subject to the initial conditions (ICs):

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , \zeta , 0\right) = f\left(\chi , \varsigma , \zeta \right) = f, \\ {u}_{2}\left(\chi , \varsigma , \zeta , 0\right) = h\left(\chi , \varsigma , \zeta \right) = h, \\ {u}_{3}\left(\chi , \varsigma , \zeta , 0\right) = g\left(\chi , \varsigma , \zeta \right) = g, \end{array}$

(1.4)

where $u = \left\langle{{u}_{1}, {u}_{2}, {u}_{3}}\right\rangle$ and $\varphi = \left\langle{f, h, g}\right\rangle$ such that ${u}_{1}$ , ${u}_{2}$ , ${u}_{3}$ , and $p$ are analytical functions of four variables $\chi, \varsigma, \zeta \ \&\ t$ .

In this equation, the solution represents the fluid velocity and pressure. It is commonly used to describe the motion of fluids in models relevant to weather, ocean currents, water flow in pipes, etc.

The novelty of this study is obvious in the proposed method when dealing with the Navier-Stokes problem, we show the simplicity and the applicability of the method, and we mention also that the method needs no differentiation, linearization, or discretization, the only mathematical step we need after taking the LT and defining the residual functions, is taking the limit at infinity which is much easier compared to other analytical techniques. Moreover, in this research, we obtain a general formula of the solution that neither researcher has, allowing us to compute as many possible terms of the series solution directly.

This study is prepared as follows: After the introduction section, a few fundamental principles and theories are reviewed for constructing an analytic series solution to the M-DT-FNS system using LRPSM. In Section 3, we constructed a Laplace residual power series (LRPS) solution to the goal problem. Three interesting examples are presented to explain the technique's simplicity and accuracy, which are displayed in Section 4. Finally, some conclusions are made about the features of the method used and its applicability in solving other types of problems.

2. Basic concepts

This part presents fundamental definitions and properties of fractional operators and power series.

Definition 2.1. ^[1] The time Caputo fractional-derivative of order $\mathrm{\alpha }$ of the multivariable function $u\left(\chi, \varsigma, \zeta, t\right)$ , is defined by

${D}_{t}^{\alpha }\left[u\left(\chi , \varsigma , \zeta , t\right)\right] = \left\{\begin{array}{c}\frac{1}{\mathrm{\Gamma }(n-\alpha )}{\int }_{0}^{t}\frac{\frac{{\partial }^{n}}{\partial {\xi }^{n}}u\left(\chi , \varsigma , \zeta , \xi \right)}{{\left(t-\xi \right)}^{\alpha +1-n}}d\xi ,\ \ n-1 < \alpha < n, \\ \frac{{\partial }^{n}}{\partial {t}^{n}}u\left(\chi , \varsigma , \zeta , t\right), \ \ \alpha = n.\end{array}\right.$

(2.1)

Lemma 2.1. ^[1] The operator ${D}_{t}^{\alpha }$ , $t\ge 0$ , $n-1 < \alpha \le n$ satisfies the following properties:

1) ${D}_{t}^{\alpha }c = 0$ , $c\in \mathbb{R}.$

2) ${D}_{t}^{\alpha }{t}^{q} = \left\{\begin{array}{c}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 0\ \ \ \ \ \ \ \ \ \ , \ \ q < \alpha, q\ {\rm{is}} \ {\rm{integer}}\\ \frac{\mathrm{\Gamma }\left(q+1\right)}{\mathrm{\Gamma }\left(q-\alpha +1\right)}{t}^{q-\alpha },\ \ {\rm{otherwise}}\end{array}, \right.$ for $q > -1.$

Definition 2.2. ^[18] The time-LT for the multivariable function $u(\chi, \varsigma, \zeta, t)$ is defined by:

$U\left(\chi , \varsigma , \zeta , s\right) = \mathcal{L}\left[u\left(\chi , \varsigma , \zeta , t\right)\right] = \int_{0}^{\infty } {e}^{-st}u\left(\chi , \varsigma , \zeta , t\right)dt, s > \delta .$

(2.2)

We denote the inverse LT of the function $U(\chi, \upsilon, \zeta, s)$ and define it as

$u\left(\chi , \varsigma , \zeta , t\right) = {\mathcal{L}}^{-1}\left[U\left(\chi , \varsigma , \zeta , s\right)\right] = \int_{{z-i\mathrm{\infty }}}^{{{z+i\mathrm{\infty }}}}{e}^{st}U\left(\chi , \varsigma , \zeta , s\right)ds;z = Re\left(s\right) > {z}_{0}.$

(2.3)

The most popular properties of the LT are mentioned below.

Lemma 2.2. ^[18,19] Assume that $U\left(\chi, \varsigma, \zeta, s\right) = \mathcal{L}\left[u\left(\chi, \varsigma, \zeta, t\right)\right]$ and $V\left(\chi, \varsigma, \zeta, s\right) = \mathcal{L}\left[v\left(\chi, \varsigma, \zeta, t\right)\right]$ . Then

1) $\underset{s\to \mathrm{\infty }}{\mathrm{lim}}sU\left(\chi, \varsigma, \zeta, s\right) = u\left(\chi, \varsigma, \zeta, 0\right).$

2) $\mathcal{L}\left[a\frac{{t}^{i\alpha }}{\mathrm{\Gamma }\left(i\alpha +1\right)}\right] = \frac{a}{{s}^{i\alpha +1}}, a\in \mathbb{R}, i = 0, \mathrm{ }1, \mathrm{ }\cdots.$

3) $\mathcal{L}\left[{D}_{t}^{m\alpha }u\left(\chi, \varsigma, \zeta, t\right)\right] = {s}^{m\alpha }U\left(\chi, \varsigma, \zeta, s\right)-\sum _{n = 0}^{m-1}{s}^{\left(m-n\right)\alpha -1}{D}_{t}^{n\alpha }u\left(\chi, \varsigma, \zeta, 0\right), \mathrm{ }0 < \alpha \le 1,$

where ${D}_{t}^{m\alpha } = {D}_{t}^{\alpha }{D}_{t}^{\alpha }\dots {D}_{t}^{\alpha }$ ( $m$ -times).

Theorem 2.1. ^[19] The power series expansion of the multivariable function $u(\chi, \varsigma, \zeta, t)$ can be expressed as

$u\left(\chi , \varsigma , \zeta , t\right) = \sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){t}^{i\alpha }}{\mathrm{\Gamma }\left(i\alpha +1\right)},$

(2.4)

where ${f}_{i}\left(\chi, \varsigma, \zeta \right) = {D}_{t}^{i\alpha }u(\chi, \varsigma, \zeta, 0)$ , $i = 0, \mathrm{ }1, \mathrm{ }\cdots$ .

Thus, the LT of Eq (2.4) is a Laurent expansion in the Laplace space of the following form:

$u\left(\chi , \varsigma , \zeta , t\right) = \sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){t}^{i\alpha }}{\mathrm{\Gamma }\left(i\alpha +1\right)}.$

(2.5)

Theorem 2.2. ^[19] Assume that

$\left|s\mathcal{L}\left[{{D}_{t}}^{(m+1)\alpha }u\left(\chi , \varsigma , \zeta , t\right)\right]\right|\le M,$

on $0\le s\le q, \mathrm{ }0 < \alpha \le 1$ and $M = M\left(\chi, \varsigma, \zeta \right)$ for some $\chi$ , $\varsigma$ , and $\zeta \in I$ . Then the remainder ${R}_{m}(\chi, \varsigma, \zeta, s)$ of the new fractional Laurent series (2.5) satisfies the following inequality

$\left|{R}_{m}(\chi , \varsigma , \zeta , s)\right|\le \frac{M}{{s}^{\left(m+1\right)\alpha +1}}.$

It is known that the LT cannot be distributed in the case of multiplication. Therefore, the following Lemma is introduced to simplify the calculations at the application of LRPSM, based on the characteristics of the powers of the power series.

Lemma 2.3. Assume that $U\left(\chi, \varsigma, \zeta, s\right) = \mathcal{L}\left[u\left(\chi, \varsigma, \zeta, t\right)\right]\left(s\right)$ and $V\left(\chi, \varsigma, \zeta, s\right) = \mathcal{L}\left[v\left(\chi, \varsigma, \zeta, t\right)\right]\left(s\right)$ . Assume that the functions $U\left(\chi, \varsigma, \zeta, s\right)$ and $V\left(\chi, \varsigma, \zeta, s\right)$ have Laurent expansions as:

$\begin{array}{l} \left(\chi , \upsilon , \zeta , s\right) = \sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \upsilon , \zeta \right)}{{s}^{i\alpha +1}}, \\ V\left(\chi , \varsigma , \zeta , s\right) = \sum _{j = 0}^{\mathrm{\infty }}\frac{{g}_{j}\left(\chi , \varsigma , \zeta \right)}{{s}^{j\alpha +1}}. \end{array}$

(2.6)

Then, $\mathcal{L}\left[u\left(\chi, \varsigma, \zeta, t\right)v\left(\chi, \varsigma, \zeta, t\right)\right]\left(s\right)$ can be expanded in a Laurent series form, as follows

$\mathcal{L}\left[u\left(\chi , \varsigma , \zeta , t\right)v\left(\chi , \varsigma , \zeta , t\right)\right]\left(s\right) = \sum _{i = 0}^{\mathrm{\infty }}\sum _{j = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){g}_{j}\left(\chi , \varsigma , \zeta \right)\mathrm{\Gamma }(\left(i+j\right)\alpha +1)}{\mathrm{\Gamma }(i\alpha +1)\mathrm{\Gamma }(j\alpha +1){s}^{(i+j)\alpha +1}}.$

(2.7)

Proof. Directly, we can prove the Lemma as follows:

$\mathcal{L}\left[u\left(\chi , \varsigma , \zeta , t\right)v\left(\chi , \varsigma , \zeta , t\right)\right]\left(s\right) = \mathcal{L}\left[{\mathcal{L}}^{-1}\left[U\left(\chi , \varsigma , \zeta , s\right)\right]{\mathcal{L}}^{-1}\left[V\left(\chi , \varsigma , \zeta , s\right)\right]\right]\left(s\right)$

$= \mathcal{L}\left[{\mathcal{L}}^{-1}\left[\sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right)}{{s}^{i\alpha +1}}\right]{\mathcal{L}}^{-1}\left[\sum _{j = 0}^{\mathrm{\infty }}\frac{{g}_{i}\left(\chi , \varsigma , \zeta \right)}{{s}^{j\alpha +1}}\right]\right]\left(s\right)$

$= \mathcal{L}\left[\sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){t}^{i\alpha }}{\mathrm{\Gamma }\left(i\alpha +1\right)}\sum _{j = 0}^{\mathrm{\infty }}\frac{{g}_{i}\left(\chi , \varsigma , \zeta \right){t}^{j\alpha }}{\mathrm{\Gamma }\left(j\alpha +1\right)}\right]\left(s\right)$

$= \mathcal{L}\left[\sum _{i = 0}^{\mathrm{\infty }}\sum _{j = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){g}_{j}\left(\chi , \varsigma , \zeta \right){t}^{\left(i+j\right)\mathrm{\alpha }}}{\mathrm{\Gamma }\left(i\alpha +1\right)\mathrm{\Gamma }\left(j\alpha +1\right)}\right]\left(s\right)$

$= \sum _{i = 0}^{\mathrm{\infty }}\sum _{j = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){g}_{j}\left(\chi , \varsigma , \zeta \right)\mathrm{\Gamma }(\left(i+j\right)\alpha +1)}{\mathrm{\Gamma }(i\alpha +1)\mathrm{\Gamma }(j\alpha +1){s}^{(i+j)\alpha +1}}.$

It is noteworthy, that we can express the $k$ th-truncated series of the Laurent series (2.7) as follows:

$\mathcal{L}\left[{u}_{k}\left(\chi , \varsigma , \zeta , t\right){v}_{k}\left(\chi , \varsigma , \zeta , t\right)\right]\left(s\right) = \sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right){g}_{j}\left(\chi , \varsigma , \zeta \right)\mathrm{\Gamma }(\left(i+j\right)\alpha +1)}{\mathrm{\Gamma }(i\alpha +1)\mathrm{\Gamma }(j\alpha +1){s}^{(i+j)\alpha +1}},$

(2.8)

which we will use extensively throughout our work on the next pages.

3. The Laplace residual power series method

We employ the LRPSM to establish a series solution for the M-DT-FNS system (1.1) in this part of this article. This technique is mainly based on applying the LT on the target equations, assuming solutions of the generated equations have Laurent expansions, and then using the idea of the limit at infinity with the residual functions to get the unknown coefficients in expansions. Finally, we run the inverse LT to obtain the solution of the given equations in the original space.

To get the LRPS solution of the system (1.1), we first apply the LT to each equation in the system (1.1) and use the third part of Lemma 2.2 with the ICs (1.2). Then, after some simplification, we get the following algebraic system in Laplace:

$\begin{array}{l} {U}_{1}-\frac{f}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{1}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{1}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{1}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{1}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \chi } = 0, \\ {U}_{2}-\frac{h}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{2}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{2}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{2}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{2}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \varsigma } = 0, \\ {U}_{3}-\frac{g}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{3}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{3}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{3}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{3}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \zeta } = 0, \end{array}$

(3.1)

where ${U}_{1} = \mathcal{L}\left[{u}_{1}\right]$ , ${U}_{2} = \mathcal{L}\left[{u}_{2}\right]$ , ${U}_{3} = \mathcal{L}\left[{u}_{3}\right]$ and $P = \mathcal{L}\left[p\right]$ .

Suppose that the solution of the system (3.1), ${U}_{1}$ , ${U}_{2}$ , and ${U}_{3}$ , has the following fractional Laurent expansions as follows:

$\begin{array}{l} {U}_{1} = \sum _{i = 0}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right)}{{s}^{i\alpha +1}}, \\ {U}_{2} = \sum _{j = 0}^{\mathrm{\infty }}\frac{{h}_{j}\left(\chi , \varsigma , \zeta \right)}{{s}^{j\alpha +1}}, \\ {U}_{3} = \sum _{r = 0}^{\mathrm{\infty }}\frac{{g}_{r}\left(\chi , \varsigma , \zeta \right)}{{s}^{r\alpha +1}}. \end{array}$

(3.2)

Assume that also $P$ has the following Laurent expansion:

$P = \frac{\varphi }{s}+\sum _{i = 1}^{\mathrm{\infty }}\frac{{\varphi }_{i}}{{s}^{i\alpha +1}}.$

(3.3)

According to Theorem 2.1 and using the ICs (1.2), the expansions in Eq (3.2) can be rewritten as follows:

$\begin{array}{l} {U}_{1} = \frac{f}{s}+\sum _{i = 1}^{\mathrm{\infty }}\frac{{f}_{i}\left(\chi , \varsigma , \zeta \right)}{{s}^{i\alpha +1}}, \\ {U}_{2} = \frac{h}{s}+\sum _{i = 1}^{\mathrm{\infty }}\frac{{h}_{i}\left(\chi , \varsigma , \zeta \right)}{{s}^{i\alpha +1}}, \\ {U}_{3} = \frac{g}{s}+\sum _{r = 1}^{\mathrm{\infty }}\frac{{g}_{r}\left(\chi , \varsigma , \zeta \right)}{{s}^{r\alpha +1}}. \end{array}$

(3.4)

The $k$ th-truncated series of the expansions in Eqs (3.3) and (3.4) are given by:

$\begin{array}{l} {U}_{1, k} = \frac{f}{s}+\sum _{i = 1}^{k}\frac{{f}_{i}}{{s}^{i\alpha +1}}, \\ {U}_{2, k} = \frac{h}{s}+\sum _{j = 1}^{k}\frac{{h}_{j}}{{s}^{j\alpha +1}}, \\ {U}_{3, k} = \frac{g}{s}+\sum _{r = 1}^{k}\frac{{g}_{r}}{{s}^{r\alpha +1}}, \\ {P}_{k} = \frac{\varphi }{s}+\sum _{i = 1}^{k}\frac{{\varphi }_{i}}{{s}^{i\alpha +1}}. \end{array}$

(3.5)

To find the coefficients in the series expansions of Eq (3.4), we establish the Laplace residual functions ( $\mathcal{L}$ RF) of the equations in the system (3.1) as follows:

$\begin{array}{l} \mathcal{L}Re{s}_{1}\left(s\right) = {U}_{1}-\frac{f}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{1}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{1}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{1}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{1}}{\partial \chi }+\frac{{\partial }^{2}{U}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{1}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \chi }, \\ \mathcal{L}Re{s}_{2}\left(s\right) = {U}_{2}-\frac{h}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{2}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{2}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{2}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{2}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \varsigma }, \\ \mathcal{L}Re{s}_{3}\left(s\right) = {U}_{3}-\frac{g}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1}\frac{\partial {u}_{3}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2}\frac{\partial {u}_{3}}{\partial \varsigma }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3}\frac{\partial {u}_{3}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{3}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial P}{\partial \zeta }, \end{array}$

(3.6)

and the kth- $\mathcal{L}$ RF for $k = \mathrm{1, 2}, \dots$ , as follows:

$\begin{array}{l} \mathcal{L}Re{s}_{1, k}\left(s\right) = {U}_{1, k}-\frac{f}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1, k}\frac{\partial {u}_{1, k}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2, k}\frac{\partial {u}_{1, k}}{\partial \varsigma }\right] \\ \ \ \ \ \ \ \ \ \ \ \ +\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3, k}\frac{\partial {u}_{1, k}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{1, k}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{1, k}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{1, k}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial {P}_{k}}{\partial \chi }, \\ \mathcal{L}Re{s}_{2, k}\left(s\right) = {U}_{2, k}-\frac{h}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1, k}\frac{\partial {u}_{2, k}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2, k}\frac{\partial {u}_{2, k}}{\partial \varsigma }\right]\\ \ \ \ \ \ \ \ \ \ \ \ +\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3, k}\frac{\partial {u}_{2, k}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{2, k}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{2, k}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{2, k}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial {P}_{k}}{\partial \varsigma }, \\ \mathcal{L}Re{s}_{3, k}\left(s\right) = {U}_{3, k}-\frac{g}{s}+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{1, k}\frac{\partial {u}_{3, k}}{\partial \chi }\right]+\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{2, k}\frac{\partial {u}_{3, k}}{\partial \varsigma }\right] \\ \ \ \ \ \ \ \ \ \ \ \ +\frac{1}{{s}^{\alpha }}\mathcal{L}\left[{u}_{3, k}\frac{\partial {u}_{3, k}}{\partial \zeta }\right]-\frac{\upsilon }{{s}^{\alpha }}\left(\frac{{\partial }^{2}{U}_{3, k}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{U}_{3, k}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{U}_{3, k}}{\partial {\zeta }^{2}}\right)+\frac{1}{{\rho s}^{\alpha }}\frac{\partial {P}_{k}}{\partial \zeta }. \end{array}$

(3.7)

Using Lemma 2.3 and substituting the expansions in Eq (3.5) into Eq (3.7) gives the Laurent series form to the kth- $\mathcal{L}$ RFs as follows:

$\begin{array}{l} \mathcal{L}Re{s}_{1, \mathrm{k}}\left(s\right) = \sum _{i = 1}^{k}\frac{{f}_{i}}{{s}^{i\alpha +1}}+\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){f}_{i}\frac{\partial {f}_{j}}{\partial \chi }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}}\\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){h}_{i}\frac{\partial {f}_{j}}{\partial \varsigma }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}} \\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){g}_{i}\frac{\partial {f}_{j}}{\partial \zeta }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}f}{\partial {\chi }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{f}_{i}}{\partial {\chi }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}\\ -\frac{\upsilon \frac{{\partial }^{2}f}{\partial {\varsigma }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{f}_{i}}{\partial {\varsigma }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}f}{\partial {\zeta }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{f}_{i}}{\partial {\zeta }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}+\frac{\frac{\partial \varphi }{\partial \chi }}{{\rho s}^{\alpha +1}}\\ +\sum _{i = 1}^{k}\frac{\frac{\partial {\varphi }_{i}}{\partial \chi }}{{\rho s}^{\left(i+1\right)\alpha +1}}, \\ \mathcal{L}Re{s}_{2, \mathrm{k}}\left(s\right) = \sum _{i = 1}^{k}\frac{{h}_{i}}{{s}^{i\alpha +1}}+\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){f}_{i}\frac{\partial {h}_{j}}{\partial \chi }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}} \\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){h}_{i}\frac{\partial {h}_{j}}{\partial \varsigma }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}} \\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){g}_{i}\frac{\partial {h}_{j}}{\partial \zeta }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}h}{\partial {\chi }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{h}_{i}}{\partial {\chi }^{2}}}{{s}^{\left(i+1\right)\alpha +1}} \\ -\frac{\upsilon \frac{{\partial }^{2}h}{\partial {\varsigma }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{h}_{i}}{\partial {\varsigma }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}h}{\partial {\zeta }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{h}_{i}}{\partial {\zeta }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}+\frac{\frac{\partial }{\partial \varsigma }\varphi }{{\rho s}^{\alpha +1}} \\ +\sum _{i = 1}^{k}\frac{\frac{\partial {\varphi }_{i}}{\partial \varsigma }}{{\rho s}^{\left(i+1\right)\alpha +1}}, \\ \mathcal{L}Re{s}_{3, \mathrm{k}}\left(s\right) = \sum _{i = 1}^{k}\frac{{g}_{i}}{{s}^{i\alpha +1}}+\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){f}_{i}\frac{\partial {g}_{j}}{\partial \chi }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}}\\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){h}_{i}\frac{\partial {g}_{j}}{\partial \varsigma }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}} \\ +\sum _{i = 0}^{k}\sum _{j = 0}^{k}\frac{\mathit{\Gamma } \left(\left(i+j\right)\alpha +1\right){g}_{i}\frac{\partial {g}_{j}}{\partial \zeta }}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(j\alpha +1\right){s}^{\left(i+j+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}h}{\partial {\chi }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{g}_{i}}{\partial {\chi }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}\\ -\frac{\upsilon \frac{{\partial }^{2}g}{\partial {\varsigma }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{g}_{i}}{\partial {\varsigma }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}-\frac{\upsilon \frac{{\partial }^{2}g}{\partial {\zeta }^{2}}}{{s}^{\alpha +1}}-\sum _{i = 1}^{k}\frac{\upsilon \frac{{\partial }^{2}{g}_{i}}{\partial {\zeta }^{2}}}{{s}^{\left(i+1\right)\alpha +1}}+\frac{\frac{\partial \varphi }{\partial \zeta }}{{\rho s}^{\alpha +1}}\\ +\sum _{i = 1}^{k}\frac{\frac{\partial {\varphi }_{i}}{\partial \zeta }}{{\rho s}^{(i+1)\alpha +1}}. \end{array}$

(3.8)

El-Ajou in ^[19] proved the validity of the following formula, which is considered the essential tool of LRPSM:

$\underset{\mathrm{s}\to \mathrm{\infty }}{\mathrm{lim}}{s}^{k\alpha +1}\mathcal{L}Re\mathrm{s}\left(s\right) = \underset{\mathrm{s}\to \mathrm{\infty }}{\mathrm{lim}}{s}^{k\alpha +1}\mathcal{L}Re{s}_{\mathrm{k}}\left(s\right) = 0, k = \mathrm{1, 2}, 3, \dots .$

(3.9)

To find the first coefficients in the expansions in Eq (3.4), consider the first truncated series of the Laurent expansions of Eq (3.5) as follows:

$\begin{array}{l} {U}_{\mathrm{1, 1}} = \frac{f}{s}+\frac{{f}_{1}}{{s}^{\alpha +1}}, \\ {U}_{\mathrm{2, 1}} = \frac{h}{s}+\frac{{h}_{1}}{{s}^{\alpha +1}}, \\ {U}_{\mathrm{3, 1}} = \frac{g}{s}+\frac{{g}_{1}}{{s}^{\alpha +1}}, \\ {P}_{1} = \frac{\varphi }{s}+\frac{{\varphi }_{1}}{{s}^{\alpha +1}}. \end{array}$

(3.10)

Substituting Eq (3.10) into the first- $\mathcal{L}$ Res, multiplying each equation by ${s}^{\alpha +1}$ , and taking the limit as $s\to \mathrm{\infty },$ we have truth in Eq (3.9):

$\begin{array}{l} \underset{\mathrm{s}\to \mathrm{\infty }}{\mathrm{lim}}{s}^{\alpha +1}\mathcal{L}Re{s}_{\mathrm{1, 1}}\left(s\right) = {f}_{1}+f\frac{\partial f}{\partial \chi }+h\frac{\partial f}{\partial \varsigma }+g\frac{\partial f}{\partial \zeta }-\upsilon \frac{{\partial }^{2}f}{\partial {\chi }^{2}}-\upsilon \frac{{\partial }^{2}f}{\partial {\varsigma }^{2}}-\upsilon \frac{{\partial }^{2}f}{\partial {\zeta }^{2}}+\frac{1}{\rho }\frac{\partial \varphi }{\partial \chi } = 0, \\ \underset{s\to \mathrm{\infty }}{\mathrm{lim}}{s}^{\alpha +1}\mathcal{L}Re{s}_{\mathrm{2, 1}}\left(s\right) = {h}_{1}+f\frac{\partial h}{\partial \chi }+h\frac{\partial h}{\partial \varsigma }+g\frac{\partial h}{\partial \zeta }-\upsilon \frac{{\partial }^{2}h}{\partial {\chi }^{2}}-\upsilon \frac{{\partial }^{2}h}{\partial {\varsigma }^{2}}-\upsilon \frac{{\partial }^{2}h}{\partial {\zeta }^{2}}+\frac{1}{\rho }\frac{\partial \varphi }{\partial \varsigma } = 0, \\ \underset{s\to \mathrm{\infty }}{\mathrm{lim}}{s}^{\alpha +1}\mathcal{L}Re{s}_{\mathrm{3, 1}}\left(\mathrm{s}\right) = {g}_{1}+f\frac{\partial g}{\partial \chi }+h\frac{\partial g}{\partial \varsigma }+g\frac{\partial g}{\partial \zeta }-\rho \frac{{\partial }^{2}g}{\partial {\chi }^{2}}-\rho \frac{{\partial }^{2}g}{\partial {\varsigma }^{2}}-\rho \frac{{\partial }^{2}g}{\partial {\zeta }^{2}}+\frac{1}{\rho }\frac{\partial \varphi }{\partial \zeta } = 0. \end{array}$

(3.11)

Therefore, the first coefficients of the expansions of Eq (3.4) become known and take the following forms:

$\begin{array}{l} {f}_{1} = \upsilon \left(\frac{\partial f}{\partial {\chi }^{2}}+\frac{\partial f}{\partial {\varsigma }^{2}}+\frac{\partial f}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial \varphi }{\partial \chi }-\left(f\frac{\partial f}{\partial \chi }+h\frac{\partial f}{\partial \varsigma }+g\frac{\partial f}{\partial \zeta }\right), \\ {h}_{1} = \upsilon \left(\frac{\partial h}{\partial {\chi }^{2}}+\frac{\partial h}{\partial {\varsigma }^{2}}+\frac{\partial h}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial \varphi }{\partial \varsigma }-\left(f\frac{\partial h}{\partial \chi }+h\frac{\partial h}{\partial \varsigma }+g\frac{\partial h}{\partial \zeta }\right), \\ {g}_{1} = \upsilon \left(\frac{\partial g}{\partial {\chi }^{2}}+\frac{\partial g}{\partial {\varsigma }^{2}}+\frac{\partial g}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial \varphi }{\partial \zeta }-\left(f\frac{\partial g}{\partial \chi }+h\frac{\partial g}{\partial \varsigma }+g\frac{\partial g}{\partial \zeta }\right). \end{array}$

(3.12)

Now, to obtain the second coefficients of the series expansions in Eq (3.4), we substitute the following second truncated series:

$\begin{array}{l} {U}_{\mathrm{1, 2}} = \frac{f}{s}+\frac{{f}_{1}}{{s}^{\alpha +1}}+\frac{{f}_{2}}{{s}^{2\alpha +1}}, \\ {U}_{\mathrm{2, 2}} = \frac{h}{s}+\frac{{h}_{1}}{{s}^{\alpha +1}}+\frac{{h}_{2}}{{s}^{2\alpha +1}}, \\ {U}_{\mathrm{3, 2}} = \frac{g}{s}+\frac{{g}_{1}}{{s}^{\alpha +1}}+\frac{{g}_{2}}{{s}^{2\alpha +1}}, \\ {P}_{2} = \frac{\varphi }{s}+\frac{{\varphi }_{1}}{{s}^{\alpha +1}}+\frac{{\varphi }_{2}}{{s}^{2\alpha +1}}, \end{array}$

(3.13)

into the 2nd- $\mathcal{L}$ RF, multiply each equation by ${s}^{2\alpha +1}$ , and use Eq (3.9), we obtain the forms of required coefficients as follows:

$\begin{array}{l} {f}_{2} = \upsilon \left(\frac{{\partial }^{2}{f}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{f}_{1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{1}}{\partial \chi }-\left({f}_{1}\frac{\partial f}{\partial \chi }+f\frac{\partial {f}_{1}}{\partial \chi }+{h}_{1}\frac{\partial f}{\partial \varsigma }+h\frac{\partial {f}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial f}{\partial \zeta }+g\frac{\partial {f}_{1}}{\partial \zeta }\right), \\ {h}_{2} = \upsilon \left(\frac{{\partial }^{2}{h}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{h}_{1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{1}}{\partial \varsigma }-\left({f}_{1}\frac{\partial h}{\partial \chi }+f\frac{\partial {h}_{1}}{\partial \chi }+{h}_{1}\frac{\partial h}{\partial \varsigma }+h\frac{\partial {h}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial h}{\partial \zeta }+g\frac{\partial {h}_{1}}{\partial \zeta }\right), \\ {g}_{2} = \upsilon \left(\frac{{\partial }^{2}{g}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{g}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{g}_{1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{1}}{\partial \zeta }-\left({f}_{1}\frac{\partial g}{\partial \chi }+f\frac{\partial {g}_{1}}{\partial \chi }+{h}_{1}\frac{\partial g}{\partial \varsigma }+h\frac{\partial {g}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial g}{\partial \zeta }+g\frac{\partial {g}_{1}}{\partial \zeta }\right). \end{array}$

(3.14)

We can repeat the previous steps to find the third and the fourth coefficients of the expansions (3.4) as:

$\begin{array}{l} {f}_{3} = \upsilon \left(\frac{{\partial }^{2}{f}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{f}_{2}}{\partial {\zeta }^{2}}\right)+\frac{1}{\rho }\frac{\partial {\varphi }_{2}}{\partial \chi }\\ \ \ \ \ \ -\left({f}_{2}\frac{\partial f}{\partial \chi }+f\frac{\partial {f}_{2}}{\partial \chi }+{h}_{2}\frac{\partial f}{\partial \varsigma }+h\frac{\partial {f}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial f}{\partial \zeta }+g\frac{\partial {f}_{2}}{\partial \zeta }\right) \\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }}^{2}\left(\alpha +1\right)}\left({f}_{1}\frac{\partial {f}_{1}}{\partial \chi }+{h}_{1}\frac{\partial {f}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial {f}_{1}}{\partial \zeta }\right), \\ {h}_{3} = \upsilon \left(\frac{{\partial }^{2}{h}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{h}_{2}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{2}}{\partial \varsigma } \\ \ \ \ \ \ -\left({f}_{2}\frac{\partial h}{\partial \chi }+f\frac{\partial {h}_{2}}{\partial \chi }+{h}_{2}\frac{\partial h}{\partial \varsigma }+h\frac{\partial {h}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial h}{\partial \zeta }+g\frac{\partial {h}_{2}}{\partial \zeta }\right)\\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }}^{2}\left(\alpha +1\right)}\left({f}_{1}\frac{\partial {h}_{1}}{\partial \chi }+{h}_{1}\frac{\partial {h}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial {h}_{1}}{\partial \zeta }\right), \\ {g}_{3} = \upsilon \left(\frac{{\partial }^{2}{g}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{g}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{g}_{2}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{2}}{\partial \zeta } \\ \ \ \ \ \ -\left({f}_{2}\frac{\partial g}{\partial \chi }+f\frac{\partial {g}_{2}}{\partial \chi }+{h}_{2}\frac{\partial g}{\partial \varsigma }+h\frac{\partial {g}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial g}{\partial \zeta }+g\frac{\partial {g}_{2}}{\partial \zeta }\right) \\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }}^{2}(\alpha +1)}\left({f}_{1}\frac{\partial {g}_{1}}{\partial \chi }+{h}_{1}\frac{\partial {g}_{1}}{\partial \varsigma }+{g}_{1}\frac{\partial {g}_{1}}{\partial \zeta }\right). \end{array}$

(3.15)

Also, one can get

$\begin{array}{l} {f}_{4} = \upsilon \left(\frac{{\partial }^{2}{f}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{f}_{3}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{3}}{\partial \chi } \\ \ \ \ \ \ -\left({f}_{3}\frac{\partial f}{\partial \chi }+f\frac{\partial {f}_{3}}{\partial \chi }+{h}_{3}\frac{\partial f}{\partial \varsigma }+h\frac{\partial {f}_{3}}{\partial \varsigma }+{g}_{3}\frac{\partial f}{\partial \zeta }+g\frac{\partial {f}_{3}}{\partial \zeta }\right) \\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(3\alpha +1\right)}{\mathrm{\Gamma }\left(2\mathrm{\alpha }+1\right)\mathrm{\Gamma }\left(\alpha +1\right)}\left({f}_{2}\frac{\partial {f}_{1}}{\partial \chi }+{f}_{1}\frac{\partial {f}_{2}}{\partial \chi }+{h}_{2}\frac{\partial {f}_{1}}{\partial \varsigma }+{h}_{1}\frac{\partial {f}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial {f}_{1}}{\partial \zeta }+{g}_{1}\frac{\partial {f}_{2}}{\partial \zeta }\right), \\ {h}_{4} = \upsilon \left(\frac{{\partial }^{2}{h}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{h}_{3}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{3}}{\partial \varsigma } \\ \ \ \ \ \ -\left({f}_{3}\frac{\partial h}{\partial \chi }+f\frac{\partial {h}_{3}}{\partial \chi }+{h}_{3}\frac{\partial h}{\partial \varsigma }+h\frac{\partial {h}_{3}}{\partial \varsigma }+{g}_{3}\frac{\partial h}{\partial \zeta }+g\frac{\partial {h}_{3}}{\partial \zeta }\right) \\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(3\alpha +1\right)}{\mathrm{\Gamma }\left(2\mathrm{\alpha }+1\right)\mathrm{\Gamma }\left(\alpha +1\right)}\left({f}_{2}\frac{\partial {h}_{1}}{\partial \chi }+{f}_{1}\frac{\partial {h}_{2}}{\partial \chi }+{h}_{2}\frac{\partial {h}_{1}}{\partial \varsigma }+{h}_{1}\frac{\partial {h}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial {h}_{1}}{\partial \zeta }+{g}_{1}\frac{\partial {h}_{2}}{\partial \zeta }\right), \\ {g}_{4} = \upsilon \left(\frac{{\partial }^{2}{g}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{g}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{g}_{3}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{3}}{\partial \zeta } \\ \ \ \ \ \ -\left({f}_{3}\frac{\partial g}{\partial \chi }+f\frac{\partial {g}_{3}}{\partial \chi }+{h}_{3}\frac{\partial g}{\partial \varsigma }+h\frac{\partial {g}_{3}}{\partial \varsigma }+{g}_{3}\frac{\partial g}{\partial \zeta }+g\frac{\partial {g}_{3}}{\partial \zeta }\right) \\ \ \ \ \ \ -\frac{\mathrm{\Gamma }\left(3\alpha +1\right)}{\mathrm{\Gamma }\left(2\mathrm{\alpha }+1\right)\mathrm{\Gamma }\left(\alpha +1\right)}\left({f}_{2}\frac{\partial {g}_{1}}{\partial \chi }+{f}_{1}\frac{\partial {g}_{2}}{\partial \chi }+{h}_{2}\frac{\partial {g}_{1}}{\partial \varsigma }+{h}_{1}\frac{\partial {g}_{2}}{\partial \varsigma }+{g}_{2}\frac{\partial {g}_{1}}{\partial \zeta }+{g}_{1}\frac{\partial {g}_{2}}{\partial \zeta }\right) \end{array}$

(3.16)

Now, we can get the form of the unknown $k$ th-coefficients in the expansions of Eq (3.4) as follows:

$\begin{array}{l} {f}_{k} = \upsilon \left(\frac{{\partial }^{2}{f}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{k-1}}{\partial \chi } \\ \ \ \ \ \ \ -\sum _{i = 0}^{k-1}\frac{\mathrm{\Gamma }\left(\left(k-1\right)\alpha +1\right)}{\mathrm{\Gamma }\left(i\alpha +1\right)\mathrm{\Gamma }\left(\left(k-i-1\right)\alpha +1\right)}\left({f}_{i}\frac{\partial {f}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {f}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {f}_{k-i-1}}{\partial \zeta }\right), k = \mathrm{1, 2}, \cdots , \\ {h}_{k} = \upsilon \left(\frac{{\partial }^{2}{h}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{k-1}}{\partial \varsigma } \\ \ \ \ \ \ \ -\sum _{i = 0}^{k-1}\frac{\mathrm{\Gamma }\left(\left(k-1\right)\alpha +1\right)}{\mathrm{\Gamma }\left(i\alpha +1\right)\mathrm{\Gamma }\left(\left(k-i-1\right)\alpha +1\right)}\left({f}_{i}\frac{\partial {h}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {h}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {h}_{\mathrm{k}-i-1}}{\partial \zeta }\right), k = \mathrm{1, 2}, \cdots , \\ {g}_{k} = \upsilon \left(\frac{{\partial }^{2}{g}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{g}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{g}_{k-1}}{\partial {\zeta }^{2}}\right)-\frac{1}{\rho }\frac{\partial {\varphi }_{k-1}}{\partial \zeta }\\ \ \ \ \ \ \ -\sum _{i = 0}^{k-1}\frac{\mathrm{\Gamma }(\left(k-1\right)\alpha +1)}{\mathrm{\Gamma }\left(i\alpha +1\right)\mathrm{\Gamma }\left(\left(k-i-1\right)\alpha +1\right)}\left({f}_{i}\frac{\partial {g}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {g}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {g}_{k-i-1}}{\partial \zeta }\right), k = \mathrm{1, 2}, \cdots . \end{array}$

(3.17)

Thus, the solution of the algebraic system (3.1) can be expressed as:

$\begin{array}{l} {U}_{1}\left(\chi , \varsigma , \zeta , s\right) = \frac{f}{s}+\frac{{f}_{1}}{{s}^{\alpha +1}}+\frac{{f}_{2}}{{s}^{2\alpha +1}}+\dots , \\ {U}_{1}\left(\chi , \varsigma , \zeta , s\right) = \frac{h}{s}+\frac{{h}_{1}}{{s}^{\alpha +1}}+\frac{{h}_{2}}{{s}^{2\alpha +1}}+\dots , \\ {U}_{3}\left(\chi , \varsigma , \zeta , s\right) = \frac{g}{s}+\frac{{g}_{1}}{{s}^{\alpha +1}}+\frac{{g}_{2}}{{s}^{2\alpha +1}}+\dots . \end{array}$

(3.18)

Finally, to get the LRPS solution of the M-DT-FNS systems (1.1) and (1.2) in the original space, we apply the inverse LT on the solution in Eq (3.18), to get

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , \zeta , t\right) = f+\frac{{f}_{1}{t}^{\alpha }}{\mathrm{\Gamma }(\alpha +1)}+\frac{{f}_{2}{t}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}+\dots , \\ {u}_{2}\left(\chi , \varsigma , \zeta , t\right) = h+\frac{{h}_{1}{t}^{\alpha }}{\mathrm{\Gamma }(\alpha +1)}+\frac{{h}_{2}{t}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}+\dots , \\ {u}_{3}\left(\chi , \varsigma , \zeta , t\right) = g+\frac{{g}_{1}{t}^{\alpha }}{\mathrm{\Gamma }(\alpha +1)}+\frac{{g}_{2}{t}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}+\dots . \end{array}$

(3.19)

4. Numerical application

In this section, we present some numerical examples that explain the working mechanism of the LRPSM. Comparisons and graphical illustrations are made to demonstrate the accuracy and efficiency of the technique.

Example 4.1. ^[35,36] Consider the following two-DT-FNSEs:

$\begin{array}{l} {D}_{t}^{\alpha }{u}_{1}+{u}_{1}\frac{\partial {u}_{1}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{1}}{\partial \varsigma } = \upsilon \left(\frac{{\partial }^{2}{u}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{1}}{\partial {\varsigma }^{2}}\right), \\ {D}_{t}^{\alpha }{u}_{2}+{u}_{1}\frac{\partial {u}_{2}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{2}}{\partial \varsigma } = \upsilon \left(\frac{{\partial }^{2}{u}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{2}}{\partial {\varsigma }^{2}}\right), \end{array}$

(4.1)

with the ICs:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , 0\right) = -\mathrm{sin}\left(\chi +\varsigma \right), \\ {u}_{2}\left(\chi , \varsigma , 0\right) = \mathrm{sin}\left(\chi +\varsigma \right), \end{array}$

(4.2)

where $\upsilon \in \mathbb{R}$ , and ${u}_{1}$ and ${u}_{2}$ are two functions of three variables $\chi, \varsigma$ , and $t$ .

Note that, when $\alpha = 1$ the exact solution of the systems (4.1) and (4.2) is

$\begin{array}{l} {u}_{1} = -{e}^{-2\upsilon t}\mathrm{sin}\left(\chi +\varsigma \right), \\ {u}_{2} = {e}^{-2\upsilon t}\mathrm{sin}\left(\chi +\varsigma \right). \end{array}$

(4.3)

Based on the algorithm of the solution obtained in Section 3 and the result in Eq (3.19), we can obtain the LRPS solution of the systems (4.1) and (4.2) as follows:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , t\right) = f+\frac{{f}_{1}{t}^{\alpha }}{\mathrm{\Gamma }(\alpha +1)}+\frac{{f}_{2}{t}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}+\dots , \\ {u}_{2}\left(\chi , \varsigma , t\right) = h+\frac{{h}_{1}{t}^{\alpha }}{\mathrm{\Gamma }(\alpha +1)}+\frac{{h}_{2}{t}^{2\alpha }}{\mathrm{\Gamma }\left(2\alpha +1\right)}+\dots , \end{array}$

(4.4)

where the coefficients ${f}_{k}$ , ${h}_{k}$ , $k = \mathrm{1, 2}, \dots$ , are given by the general formula obtained in Eq (3.17) which in the two dimensions will be as follows:

$\begin{array}{l} {f}_{k} = \upsilon \left(\frac{{\partial }^{2}{f}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial \varsigma }\right)-\sum _{i = 0}^{k-1}\frac{\mathit{\Gamma } \left(\left(k-1\right)\alpha +1\right)\left({f}_{i}\frac{\partial {f}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {f}_{k-i-1}}{\partial \varsigma }\right)}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(\left(k-i-1\right)\alpha +1\right)}, \\ {h}_{k} = \upsilon \left(\frac{{\partial }^{2}{h}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\varsigma }^{2}}\right)-\sum _{i = 0}^{k-1}\frac{\mathit{\Gamma } \left(\left(k-1\right)\alpha +1\right)\left({f}_{i}\frac{\partial {h}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {fh}_{k-i-1}}{\partial \varsigma }\right)}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(\left(k-i-1\right)\alpha +1\right)}. \end{array}$

(4.5)

Thus, we conclude that

${f}_{1} = 2\upsilon \mathrm{sin}\left(\chi +\varsigma \right), \ \ \ \ \ \ \ {h}_{1} = -2\upsilon \mathrm{sin}\left(\chi +\varsigma \right)$

${f}_{2} = -{\left(2\upsilon \right)}^{2}\mathrm{sin}\left(\chi +\varsigma \right), \ \ \ \ \ \ \ {h}_{2} = {\left(2\upsilon \right)}^{2}\mathrm{sin}\left(\chi +\varsigma \right)$

${f}_{3} = {\left(2\upsilon \right)}^{3}\mathrm{sin}\left(\chi +\varsigma \right), \ \ \ \ \ \ \ {h}_{3} = {-\left(2\upsilon \right)}^{3}\mathrm{sin}\left(\chi +\varsigma \right)$

${f}_{4} = -{\left(2\upsilon \right)}^{4}\mathrm{sin}\left(\chi +\varsigma \right), \ \ \ \ \ \ \ {h}_{4} = {\left(2\upsilon \right)}^{4}\mathrm{sin}\left(\chi +\varsigma \right)$

$\ \ \ \ \ \ \ \ \ \ \ ⋮\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ ⋮$

${f}_{k} = {\left(-1\right)}^{k+1}{\left(2\upsilon \right)}^{k}\mathrm{sin}\left(\chi +\varsigma \right), \ \ \ \ \ \ \ {h}_{k} = {\left(-2\upsilon \right)}^{k}\mathrm{sin}\left(\chi +\varsigma \right),\ \ \ \ \ \ \ k = \mathrm{1, 2}, \cdots .$

Hence, we can express the LRPS solution of Eqs (4.1) and (4.2) as:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , t\right) = -\mathrm{sin}\left(\chi +\varsigma \right)\sum _{n = 0}^{\mathrm{\infty }}\frac{{\left(-2\upsilon \right)}^{n}{t}^{n\alpha }}{\mathrm{\Gamma }\left(n\alpha +1\right)}, \\ {u}_{2}\left(\chi , \varsigma , t\right) = \mathrm{sin}\left(\chi +\varsigma \right)\sum _{n = 0}^{\mathrm{\infty }}\frac{{\left(-2\upsilon \right)}^{n}{t}^{n\alpha }}{\mathrm{\Gamma }\left(n\alpha +1\right)}. \end{array}$

(4.6)

shows the stability of the results obtained, where the solution is stable with various values of $t$ .

Table 1. Numerical results for Example 1 at

$\alpha = 1, \chi = 0.1, \varsigma = 0.5$ .

$t$	Numerical results of ${u}_{1}\left(\chi, \varsigma, t\right)$
	${u}_{\mathrm{1, 6}}\left(\chi, \varsigma, t\right)$	${u}_{\mathrm{1, 8}}\left(\chi, \varsigma, t\right)$	${u}_{\mathrm{1, 10}}\left(\chi, \varsigma, t\right)$	${u}_{1}\left(\chi, \varsigma, t\right)$	$\left\|{u}_{1}\left(\chi, \varsigma, t\right)-{u}_{1, 10}\left(\chi, \varsigma, t\right)\right\|$
$0.16$	$\begin{array}{c}-0.481157\end{array}$	$\begin{array}{c}-0.481157\end{array}$	$\begin{array}{c}-0.481157\end{array}$	$-0.481157$	$5.551115\times {10}^{-17}$
$0.32$	$\begin{array}{c}-0.410015\end{array}$	$\begin{array}{c}-0.410015\end{array}$	$\begin{array}{c}-0.410015\end{array}$	$-0.410015$	$4.957146\times {10}^{-14}$
$0.48$	$\begin{array}{c}-0.349392\end{array}$	$\begin{array}{c}-0.349391\end{array}$	$\begin{array}{c}-0.349391\end{array}$	$-0.349391$	$4.238276\times {10}^{-12}$
$0.64$	$\begin{array}{c}-0.297736\end{array}$	$\begin{array}{c}-0.297732\end{array}$	$\begin{array}{c}-0.297731\end{array}$	$-0.297732$	$9.907064\times {10}^{-11}$
$0.80$	$\begin{array}{c}-0.253732\end{array}$	$\begin{array}{c}-0.253710\end{array}$	$\begin{array}{c}-0.253710\end{array}$	$-0.253710$	$1.138797\times {10}^{-9}$
$0.96$	$\begin{array}{c}-0.216273\end{array}$	$\begin{array}{c}-0.216199\end{array}$	$\begin{array}{c}-0.216198\end{array}$	$-0.216198$	$8.355848\times {10}^{-9}$
$t$	Numerical results of ${u}_{2}\left(\chi, \varsigma, t\right)$
	${u}_{\mathrm{2, 6}}\left(\chi, \varsigma, t\right)$	${u}_{\mathrm{2, 8}}\left(\chi, \varsigma, t\right)$	${u}_{\mathrm{2, 10}}\left(\chi, \varsigma, t\right)$	${u}_{2}\left(\chi, \varsigma, t\right)$	$\left\|{u}_{2}\left(\chi, \varsigma, t\right)-{u}_{2, 10}\left(\chi, \varsigma, t\right)\right\|$
$0.16$	$\begin{array}{c}0.481157\end{array}$	$\begin{array}{c}0.481157\end{array}$	$\begin{array}{c}0.481157\end{array}$	$0.481157$	$5.551115\times {10}^{-17}$
$0.32$	$\begin{array}{c}0.410015\end{array}$	$\begin{array}{c}0.410015\end{array}$	$\begin{array}{c}0.410015\end{array}$	$0.410015$	$4.957146\times {10}^{-14}$
$0.48$	$\begin{array}{c}0.349392\end{array}$	$\begin{array}{c}0.349391\end{array}$	$\begin{array}{c}0.349391\end{array}$	$0.349391$	$4.238276\times {10}^{-12}$
$0.64$	$\begin{array}{c}0.297736\end{array}$	$\begin{array}{c}0.297732\end{array}$	$\begin{array}{c}0.297731\end{array}$	$0.297732$	$9.907064\times {10}^{-11}$
$0.80$	$\begin{array}{c}0.253732\end{array}$	$\begin{array}{c}0.253710\end{array}$	$\begin{array}{c}0.253710\end{array}$	$0.253710$	$1.138797\times {10}^{-9}$
$0.96$	$\begin{array}{c}0.216273\end{array}$	$\begin{array}{c}0.216199\end{array}$	$\begin{array}{c}0.216198\end{array}$	$0.216198$	$8.355848\times {10}^{-9}$

| Show Table

DownLoad: CSV

This solution is the same as that obtained by the Laplace decomposition method ^[36] and the variational iteration transform method ^[35]. In a special case, taking $\alpha = 1$ gives the exact solution in terms of elementary functions as follows:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , t\right) = -{e}^{-2\upsilon t}\mathrm{sin}\left(\chi +\varsigma \right), \\ {u}_{2}\left(\chi , \varsigma , t\right) = {e}^{-2\upsilon t}\mathrm{sin}\left(\chi +\varsigma \right). \end{array}$

(4.7)

The behavior of the velocity field of the two-DT-FNSEs (4.1) and (4.2) is depicted in for various values of $\alpha$ at $t = 0.5$ and $\upsilon = 0.5$ . The 10th-truncated series of Eq (4.6) is plotted in - for $\alpha = 0.6$ , $\alpha = 0.8$ , and $\alpha = 1$ , respectively, whereas, the exact solution at $\alpha = 1$ is plotted in . The graphics indicate consistency in the behavior of the solution at various values of $\alpha$ , as well as the convention of the exact solution with the obtained solution in Figure 1(c, d).

Figure 1. The 3D surface plot of the 10th approximate solutions of

${u}_{1}$ and

${u}_{2}$ at different values of

$\alpha$ and

$t = 0.5$ &

$\upsilon = 0.5$ for the problem in Example 4.1. (a)

$\alpha = 0.6$ , (b)

$\alpha = 0.8$ , (c)

$\alpha = 1$ , (d)

$\alpha = 1$ (Exact solutions).

DownLoad: Full-Size Img PowerPoint

shows the action of the 10th approximate analytical solution of the initial value problems (IVP) (4.1) and (4.2) along the line $\varsigma = \chi$ and in the region $D = \{\left(\chi, t\right):-3\le \chi \le 3, \mathrm{ }0\le t < 1\}$ for distinct values of $\alpha$ , and at $\upsilon = 0.5$ . The 10th approximate solution is plotted in - for $\alpha = 0.6$ , $\alpha = 0.8$ , and $\alpha = 1$ , respectively, whereas, the exact solution at $\alpha = 1$ is plotted in (d). Also, the graphics indicate consistency in the action of the solution at distinct values of $\alpha$ , the accord nation of the exact solution with the approximate solution in Figure 2(c, d) as well as the region of convergence of the series solution.

Figure 2. The graph of the 3D surface of the 10th approximate solutions of

${u}_{1}$ and

${u}_{2}$ along the line

$\varsigma = \chi$ and at various values of

$\alpha$ and

$\upsilon = 0.5$ for the problem in Example 4.1. (a)

$\alpha = 0.6$ , (b)

$\alpha = 0.8$ , (c)

$\alpha = 1$ , (d)

$\alpha = 1$ (Exact).

DownLoad: Full-Size Img PowerPoint

Example 4.2. ^[35,36] Consider the following two-DT-FNSEs:

(4.8)

with the ICs:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , 0\right) = -{e}^{\chi +\varsigma }, \\{u}_{2}\left(\chi , \varsigma , 0\right) = {e}^{\chi +\varsigma }, \end{array}$

(4.9)

where $\upsilon \in \mathbb{R}$ and ${u}_{1}$ and ${u}_{2}$ are two functions of three variables $\chi, \varsigma,$ and $t$ .

The exact solution to problems (4.8) and (4.9) can be obtained, when putting $\alpha = 1$ , to be ${u}_{1} = -{e}^{-2\upsilon t+\chi +\varsigma }$ and ${u}_{2} = {e}^{-2\upsilon t+\chi +\varsigma }$ . Applying the same procedure in Example 4.1, one can obtain the following recurrence relations:

$\begin{array}{l} {f}_{k} = \upsilon \left(\frac{{\partial }^{2}{f}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial {\varsigma }^{2}}\right)-\sum _{i = 0}^{k-1}\frac{\mathit{\Gamma } \left((k-1)\alpha +1\right)\left({f}_{i}\frac{\partial {f}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {f}_{k-i-1}}{\partial \varsigma }\right)}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(\left(k-i-1\right)\alpha +1\right)}, \\{h}_{k} = \upsilon \left(\frac{{\partial }^{2}{h}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\varsigma }^{2}}\right)-\sum _{i = 0}^{k-1}\frac{\mathit{\Gamma } \left(\left(k-1\right)\alpha +1\right)\left({f}_{i}\frac{\partial {h}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {fh}_{k-i-1}}{\partial \varsigma }\right)}{\mathit{\Gamma } \left(i\alpha +1\right)\mathit{\Gamma } \left(\left(k-i-1\right)\alpha +1\right)}. \end{array}$

(4.10)

Thus, we conclude that

${f}_{1} = -2\upsilon {e}^{\chi +\varsigma },\ \ \ \ \ \ {h}_{1} = 2\upsilon {e}^{\chi +\varsigma }$

${f}_{2} = -{\left(2\upsilon \right)}^{2}{e}^{\chi +\varsigma }, \ \ \ \ \ \ {h}_{2} = {\left(2\upsilon \right)}^{2}{e}^{\chi +\varsigma }$

${f}_{3} = -{\left(2\upsilon \right)}^{3}{e}^{\chi +\varsigma }, \ \ \ \ \ \ {h}_{3} = {\left(2\upsilon \right)}^{3}{e}^{\chi +\varsigma },$

${f}_{4} = -{\left(2\upsilon \right)}^{4}{e}^{\chi +\varsigma },\ \ \ \ \ \ {h}_{3} = {\left(2\upsilon \right)}^{3}{e}^{\chi +\varsigma }$

$\ \ \ \ \ \ \ \ \ \ \ \ ⋮\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ ⋮$

${f}_{k} = -{\left(2\upsilon \right)}^{k}{e}^{\chi +\varsigma }, \ \ \ \ \ \ {h}_{k} = {\left(2\upsilon \right)}^{k}{e}^{\chi +\varsigma }, \ \ \ \ \ \ k = \mathrm{1, 2}, \cdots .$

Thus, the LRPS solution of the systems (4.8) and (4.9) can be expressed as follows:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , t\right) = -{e}^{\chi +\varsigma }\sum _{i = 0}^{\mathrm{\infty }}\frac{{\left(-2\upsilon \right)}^{i}{t}^{i\alpha }}{\mathrm{\Gamma }\left(i\alpha +1\right)}, \\ {u}_{2}\left(\chi , \varsigma , t\right) = {e}^{\chi +\varsigma }\sum _{j = 0}^{\mathrm{\infty }}\frac{{\left(-2\upsilon \right)}^{j}{t}^{j\alpha }}{\mathrm{\Gamma }\left(j\alpha +1\right)}. \end{array}$

(4.11)

For $\alpha = 1$ , the solution in Eq (4.11) has the form:

${u}_{1}\left(\chi , \varsigma , t\right) = -{e}^{-2\upsilon t+\chi +\varsigma },$

${u}_{2}\left(\chi , \varsigma , t\right) = {e}^{-2\upsilon t+\chi +\varsigma }.$

shows the velocity field behavior of the two-DT-FNSEs (4.8) and (4.9) for distinct values of $\alpha$ at $t = 0.5$ and $\upsilon = 0.5$ . The 10th LRPS approximate analytical solution of the IVP (4.8) and (4.8) plotted in - for $\alpha = 0.6$ , $\alpha = 0.8,$ and $\alpha = 1$ respectively, while the exact solution at $\alpha = 1$ is plotted in (d). The graphics indicate the consistency in the solution behavior at various values of $\alpha$ , as well as the exact solution agreement with the proposed analytical solution in Figure 3(c, d).

Figure 3. The 3D surface plot of the 10th approximate solutions of

${u}_{1}$ and

${u}_{2}$ at distinct values of

$\alpha$ and

$t = 0.5$ &

$\upsilon = 0.5$ for the problem in Example 4.2. (a)

$\alpha = 0.6$ , (b)

$\alpha = 0.8$ , (c)

$\alpha = 1$ , (d)

$\alpha = 1$ (Exact solutions).

DownLoad: Full-Size Img PowerPoint

illustrates the behavior of the 10th approximate solution of the IVP (4.8) and (4.9) along the line $\varsigma = \chi$ and in the region $D = \{\left(\chi, t\right):-1\le \chi \le \mathrm{1, 0}\le t < 1\}$ for different values of $\alpha$ , and at $\upsilon = 0.5$ . The 10th approximate solution is plotted in - for $\alpha = 0.6$ , $\alpha = 0.8$ , and $\alpha = 1$ , respectively, whereas, the exact solution at $\alpha = 1$ is plotted in (d). Also, the graphics indicate consistency in the action of the solution at distinct values of $\alpha$ , the coordination between the exact solution and the approximate analytical solution as illustrated in Figure 4(c, d) as well as the determination of the region of convergence for the series solution, is clear.

Figure 4.

${u}_{1}$ and

${u}_{2}$ along the line

$\varsigma = \chi$ and at various of

$\alpha$ and

$\upsilon = 0.5$ for the problem in Example 4.2. (a)

$\alpha = 0.6$ , (b)

$\alpha = 0.8$ , (c)

$\alpha = 1$ , (d)

$\alpha = 1$ (Exact).

DownLoad: Full-Size Img PowerPoint

Example 4.3. ^[36,41] Consider the following three-DT-FNSEs:

$\begin{array}{l} {D}_{t}^{\alpha }{u}_{1}+{u}_{1}\frac{\partial {u}_{1}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{1}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{1}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{1}}{\partial {\zeta }^{2}}\right), \\ {D}_{t}^{\alpha }{u}_{2}+{u}_{1}\frac{\partial {u}_{2}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{2}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{2}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{2}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{2}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{2}}{\partial {\zeta }^{2}}\right), \\ {D}_{t}^{\alpha }{u}_{3}+{u}_{1}\frac{\partial {u}_{3}}{\partial \chi }+{u}_{2}\frac{\partial {u}_{3}}{\partial \varsigma }+{u}_{3}\frac{\partial {u}_{3}}{\partial \zeta } = \upsilon \left(\frac{{\partial }^{2}{u}_{3}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{u}_{3}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{u}_{3}}{\partial {\zeta }^{2}}\right), \end{array}$

(4.12)

where $v\in \mathbb{R}$ , and with the ICs:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , \zeta , 0\right) = -0.5\chi +\varsigma +\zeta , \\ {u}_{2}\left(\chi , \varsigma , \zeta , 0\right) = \chi -0.5\varsigma +\zeta , \\ {u}_{3}\left(\chi , \varsigma , \zeta , 0\right) = \chi +\varsigma -0.5\zeta , \end{array}$

(4.13)

and the exact solution, when $\alpha = 1$ , is

${u}_{1} = \frac{-0.5\chi +\varsigma +\zeta -2.25\chi t}{1-2.25{t}^{2}}, {u}_{2} = \frac{\chi -0.5\varsigma +\zeta -2.25\varsigma t}{1-2.25{t}^{2}}, {u}_{3} = \frac{\chi +\varsigma -0.5\zeta -2.25\zeta t}{1-2.25{t}^{2}}.$

The general formula (3.12) for the systems (4.12) and (4.13) will be as follows:

${f}_{k} = \upsilon \left(\frac{{\partial }^{2}{f}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{f}_{k-1}}{\partial {\zeta }^{2}}\right)$

$-\sum\limits_{i = 0}^{k-1}\frac{\varGamma \left((k-1)\alpha +1\right)\left({f}_{i}\frac{\partial {f}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {f}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {f}_{k-i-1}}{\partial \zeta }\right)}{\varGamma \left(i\alpha +1\right)\varGamma \left(\left(k-i-1\right)\alpha +1\right)},$

${h}_{k} = \upsilon \left(\frac{{\partial }^{2}{h}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{h}_{k-1}}{\partial {\zeta }^{2}}\right)$

$-\sum\limits_{i = 0}^{k-1}\frac{\varGamma \left(\left(k-1\right)\alpha +1\right)\left({f}_{i}\frac{\partial {h}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {h}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {h}_{k-i-1}}{\partial \zeta }\right)}{\varGamma \left(i\alpha +1\right)\varGamma \left(\left(k-i-1\right)\alpha +1\right)},$

${g}_{k} = \upsilon \left(\frac{{\partial }^{2}{g}_{k-1}}{\partial {\chi }^{2}}+\frac{{\partial }^{2}{g}_{k-1}}{\partial {\varsigma }^{2}}+\frac{{\partial }^{2}{g}_{k-1}}{\partial {\zeta }^{2}}\right)$

$-\sum\limits_{i = 0}^{k-1}\frac{\varGamma \left(\left(k-1\right)\alpha +1\right)\left({f}_{i}\frac{\partial {g}_{k-i-1}}{\partial \chi }+{h}_{i}\frac{\partial {g}_{k-i-1}}{\partial \varsigma }+{g}_{i}\frac{\partial {g}_{k-i-1}}{\partial \zeta }\right)}{\varGamma \left(i\alpha +1\right)\varGamma \left(\left(k-i-1\right)\alpha +1\right)}.$

Thus, we conclude that

${f}_{1} = -2.25\chi , \ \ \ \ \ \ \ {h}_{1} = -2.25\varsigma ,\ \ \ \ \ \ \ {g}_{1} = -2.25\zeta ,$

${f}_{2} = -2.25\chi +4.5\varsigma +4.5\zeta ,\ \ \ \ \ \ \ {h}_{2} = 4.5\chi -2.25\varsigma +4.5\zeta ,\ \ \ \ \ \ \ {g}_{2} = 4.5\chi +4.5\varsigma -2.25\zeta ,$

${f}_{3} = -\chi \left(\frac{5.0625\varGamma \left(2\alpha +1\right)}{{\varGamma \left(\alpha +1\right)}^{2}}+20.25\right),$

${h}_{3} = -\varsigma \left(\frac{5.0625\varGamma \left(2\alpha +1\right)}{{\varGamma \left(\alpha +1\right)}^{2}}+20.25\right),$

${g}_{3} = -\zeta \left(\frac{5.0625\varGamma \left(2\alpha +1\right)}{{\varGamma \left(\alpha +1\right)}^{2}}+20.25\right),$

${f}_{4} = \frac{\varGamma \left(2\alpha +1\right)\left(-5.0625\chi +10.125\varsigma +10.125\zeta \right)}{{\varGamma \left(\alpha +1\right)}^{2}}+\frac{\varGamma \left(3\alpha +1\right)\left(-10.125\chi +20.25\varsigma +20.25\zeta \right)}{\varGamma \left(\alpha +1\right)\varGamma \left(2\alpha +1\right)}$

$-20.25\chi +40.5\varsigma +40.5\zeta ,$

${h}_{4} = \frac{\varGamma \left(2\alpha +1\right)\left(10.125\chi -5.0625\varsigma +10.125\zeta \right)}{{\varGamma \left(\alpha +1\right)}^{2}}+\frac{\varGamma \left(3\alpha +1\right)\left(20.25\chi -10.125\varsigma +20.25\zeta \right)}{\varGamma \left(\alpha +1\right)\varGamma \left(2\alpha +1\right)}$

$+40.5\chi -20.25\varsigma +40.5\zeta ,$

${g}_{4} = \frac{\mathrm{\Gamma }\left(2\alpha +1\right)\left(10.125\chi +10.125\varsigma -5.0625\zeta \right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{2}}+\frac{\mathrm{\Gamma }\left(3\alpha +1\right)\left(20.25\chi +20.25\varsigma -10.125\zeta \right)}{\mathrm{\Gamma }\left(\alpha +1\right)\mathrm{\Gamma }\left(2\alpha +1\right)}$

$+40.5\chi +40.5\varsigma -20.25\zeta , \text{}$

${f}_{5} = \chi \left(\begin{array}{l}-182.25-\frac{22.78125\mathrm{\Gamma }\left(4\alpha +1\right)\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{3}\mathrm{\Gamma }\left(3\alpha +1\right)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{2}}\\ +\frac{-\frac{91.125\mathrm{ }\mathrm{\Gamma }(3\alpha +1)}{\mathrm{\Gamma }(2\alpha +1)}-\frac{91.125\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{\mathrm{\Gamma }(3\alpha +1)}}{\mathrm{\Gamma }(\alpha +1)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{{\mathrm{\Gamma }(2\alpha +1)}^{2}}\end{array}\right),$

${h}_{5} = \varsigma \left(\begin{array}{l}-182.25-\frac{22.78125\mathrm{\Gamma }\left(4\alpha +1\right)\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{3}\mathrm{\Gamma }\left(3\alpha +1\right)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{2}}\\ +\frac{-\frac{91.125\mathrm{ }\mathrm{\Gamma }(3\alpha +1)}{\mathrm{\Gamma }(2\alpha +1)}-\frac{91.125\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{\mathrm{\Gamma }(3\alpha +1)}}{\mathrm{\Gamma }(\alpha +1)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{{\mathrm{\Gamma }(2\alpha +1)}^{2}}\end{array}\right),$

${g}_{5} = \zeta \left(\begin{array}{l}-182.25-\frac{22.78125\mathrm{\Gamma }\left(4\alpha +1\right)\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{3}\mathrm{\Gamma }\left(3\alpha +1\right)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }\left(2\alpha +1\right)}{{\mathrm{\Gamma }\left(\alpha +1\right)}^{2}}\\ +\frac{-\frac{91.125\mathrm{ }\mathrm{\Gamma }(3\alpha +1)}{\mathrm{\Gamma }(2\alpha +1)}-\frac{91.125\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{\mathrm{\Gamma }(3\alpha +1)}}{\mathrm{\Gamma }(\alpha +1)}-\frac{45.5625\mathrm{ }\mathrm{\Gamma }(4\alpha +1)}{{\mathrm{\Gamma }(2\alpha +1)}^{2}}\end{array}\right).$

So, the LRPS solution of the systems (4.12) and (4.13) has the following series form:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , \zeta , t\right) = -0.5\chi +\varsigma +\zeta -\frac{2.25}{\mathrm{\Gamma }\left(1+\alpha \right)}\chi {t}^{\alpha }+\frac{2\left(2.25\right)}{\mathrm{\Gamma }\left(1+2\alpha \right)}\left(-0.5\chi +\varsigma +\zeta \right){t}^{2\alpha } \\ \ \ \ \ \ \ \ \ \ \ \ \ -\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+3\alpha \right)}\left(4+\frac{\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right)\chi {t}^{3\alpha }+\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+4\alpha \right)}\left(8+\frac{2\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right. \\ \ \ \ \ \ \ \ \ \ \ \ \ \left.+\frac{4\mathrm{\Gamma }\left(1+3\alpha \right)}{\mathrm{\Gamma }\left(1+\alpha \right)\mathrm{\Gamma }\left(1+2\alpha \right)}\right)\left(-0.5\chi +\varsigma +\zeta \right){t}^{4\alpha }+\cdots \mathrm{ }, \\ {u}_{2}\left(\chi , \varsigma , \zeta , t\right) = \chi -0.5\varsigma +\zeta -\frac{2.25}{\mathrm{\Gamma }\left(1+\alpha \right)}\varsigma {t}^{\alpha }+\frac{2\left(2.25\right)}{\mathrm{\Gamma }\left(1+2\alpha \right)}\left(\chi -0.5\varsigma +\zeta \right){t}^{2\alpha }\\ \ \ \ \ \ \ \ \ \ \ \ \ -\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+3\alpha \right)}\left(4+\frac{\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right)\varsigma {t}^{3\alpha }+\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+4\alpha \right)}\left(8+\frac{2\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right. \\ \ \ \ \ \ \ \ \ \ \ \ \ \left.+\frac{4\mathrm{\Gamma }\left(1+3\alpha \right)}{\mathrm{\Gamma }\left(1+\alpha \right)\mathrm{\Gamma }\left(1+2\alpha \right)}\right)\left(\chi -0.5\varsigma +\zeta \right){t}^{4\alpha }+\cdots , \\ {u}_{3}\left(\chi , \varsigma , \zeta , t\right) = \chi +\varsigma -0.5\zeta -\frac{2.25}{\mathrm{\Gamma }\left(1+\alpha \right)}\zeta {t}^{\alpha }+\frac{2\left(2.25\right)}{\mathrm{\Gamma }\left(1+2\alpha \right)}\left(\chi +\varsigma -0.5\zeta \right){t}^{2\alpha }\\ \ \ \ \ \ \ \ \ \ \ \ \ -\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+3\alpha \right)}\left(4+\frac{\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right)\zeta {t}^{3\alpha }+\frac{{\left(2.25\right)}^{2}}{\mathrm{\Gamma }\left(1+4\alpha \right)}\left(8+\frac{2\mathrm{\Gamma }\left(1+2\alpha \right)}{{\mathrm{\Gamma }}^{2}\left(1+\alpha \right)}\right.\\ \ \ \ \ \ \ \ \ \ \ \ \ \left.+\frac{4\mathrm{\Gamma }\left(1+3\alpha \right)}{\mathrm{\Gamma }\left(1+\alpha \right)\mathrm{\Gamma }\left(1+2\alpha \right)}\right)\left(\chi +\varsigma -0.5\zeta \right){t}^{4\alpha }+\cdots . \end{array}$

(4.14)

The sum of the series in Eq (4.14) at $\alpha = 1$ has the following rational forms:

$\begin{array}{l} {u}_{1}\left(\chi , \varsigma , \zeta , t\right) = \frac{-0.5\chi +\varsigma +\zeta -2.25\chi t}{1-2.25{t}^{2}}, \\ {u}_{2}\left(\chi , \varsigma , \zeta , t\right) = \frac{\chi -0.5\varsigma +\zeta -2.25\varsigma t}{1-2.25{t}^{2}}, \\ {u}_{3}\left(\chi , \varsigma , \zeta , t\right) = \frac{\chi +\varsigma -0.5\zeta -2.25\zeta t}{1-2.25{t}^{2}}. \end{array}$

(4.15)

The behavior of the velocity field of the Three-DT-FNSEs (4.12) and (4.13) is depicted in Figure 5 for various values of at and. The th-truncated series of Eq (4.10) is plotted in Figure 5(a-c) for, and, respectively, whereas, the exact solution at is plotted in (d). The graphics indicate consistency in the behavior of the solution at various values of, as well as the agreement of the exact solution with the approximate solution in Figure 5(c, d).

Figure 5. The 3D surfaces plot of the 10th approximate solutions of

${u}_{1}, {u}_{2}$ , and

${u}_{3}$ at various values of

$\alpha$ and

$t = 0.5$ &

$\zeta = 3$ for the problem in Example 4.3. (a)

$\alpha = 0.6$ , (b)

$\alpha = 0.8$ , (c)

$\alpha = 1$ , (d)

$\alpha = 1$ (Exact solutions).

DownLoad: Full-Size Img PowerPoint

5. Mathematical reviews

This article presents the LRPSM in a new scheme. We proposed the method and used it to solve the M-DT-FNS system. In the following, we state the advantages of using the presented method and the disadvantages of treating the M-DT-FNS.

5.1. Advantages of the method

1) The method is simple to apply to solve linear and non-linear FPDE compared to other techniques, other power series methods are based on finding derivatives and the calculations are usually complex, but LRPSM mainly depends on computing the limit at infinity which is much easier.

2) The proposed method is applicable in finding approximate solutions for physical applications, and in finding many terms of the analytical series solutions.

3) The method is accurate and gives approximate solutions close to the exact ones.

5.2. Disadvantages of the method

LRPSM needs first to find the LT of the target equations and finally to run the inverse LT to obtain the solution in the original space. So, if we have nonhomogeneous equations, the source functions need to be piecewise continuous and of exponential order, and after the computations, the inverse LT must exist.

6. Conclusions

In this article, we have introduced the LRPSM in a new scheme and simplified the technique to present series solutions for the M-DT-FNS system in the sense of the Caputo derivative. It is worth noting here that we obtained a general formula for an analytic solution of M-DT-FN, which other researchers have not previously obtained by other methods. We tested three examples by solving them in the proposed technique and then analyzing the results. In the future, we will use LRPSM to solve more problems and make new modifications to address the flaws of the presented technique.

Acknowledgments

The authors express their gratitude to the dear referees, who wish to remain anonymous, and the editor for their helpful suggestions, which improved the final version of this paper.

Conflict of interest

The authors declare no conflicts of interest.

Acknowledgments

This work was supported by the Ministry of Science and Higher Education of the Russian Federation (2019/# 0545-2019-0011).

Conflict of interest

All authors declare no conflicts of interest in this paper.

References

[1]	Ryan MP, O'Dwyer J, Adley CC (2017) Evaluation of the complex nomenclature of the clinically and veterinary significant pathogen Salmonella. BioMed Res Int 2017: 3782182.
[2]	Lamas A, Miranda JM, Regal P, et al. (2018) A comprehensive review of non-enterica subspecies of Salmonella enterica. Microbiol Res 206: 60-73. doi: 10.1016/j.micres.2017.09.010
[3]	Shubin FN, Rakov AV, Kuznetsova NA, et al. (2015) Zoonotic salmonellosis in the Far East: main aspects of a problem. Far-Eastern J Infect Pathol 29: 72-79.
[4]	EFSA and ECDC (European Food Safety Authority and European Centre for Disease Prevention and Control) (2018) The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2017. EFSA J 16: 5500.
[5]	Hendriksen RS, Vieira AR, Karlsmose S, et al. (2011) Global monitoring of Salmonella serovar distribution from the world health organization global foodborne infections network country data bank: results of quality assured laboratories from 2001 to 2007. Foodborne Pathog Dis 8: 887-900. doi: 10.1089/fpd.2010.0787
[6]	Shubin FN (2015) Zoonotic salmonellosis in Russia: main aspects of a problem. Epidemiol Vaccinal Prev 80: 28-30. doi: 10.31631/2073-3046-2015-14-1-28-30
[7]	Rakov AV, Mastriani E, Liu SL, et al. (2019) Association of Salmonella virulence factor alleles with intestinal and invasive serovars. BMC Genomics 20: 429. doi: 10.1186/s12864-019-5809-8
[8]	Dagan T, Blekhman R, Graur D (2006) The ‘domino theory’ of gene death: gradual and mass gene extinction events in three lineages of obligate symbiotic bacterial pathogens. Mol Biol Evol 23: 310-316. doi: 10.1093/molbev/msj036
[9]	Nuccio SP, Bäumler AJ (2014) Comparative analysis of Salmonella genomes identifies a metabolic network for escalating growth in the inflamed gut. mBio 5: e00929-14. doi: 10.1128/mBio.00929-14
[10]	Valenzuela LM, Hidalgo AA, Rodríguez L, et al. (2015) Pseudogenization of sopA and sopE2 is functionally linked and contributes to virulence of Salmonella enterica serovar Typhi. Infect Genet Evol 33: 131-142. doi: 10.1016/j.meegid.2015.04.021
[11]	Ortega AP, Villagra NA, Urrutia IM, et al. (2016) Lose to win: marT pseudogenization in Salmonella enterica serovar Typhi contributed to the surV-dependent survival to H₂O₂, and inside human macrophage-like cells. Infect Genet Evol 45: 111-121. doi: 10.1016/j.meegid.2016.08.029
[12]	Call DR, Kang MS, Daniels J, et al. (2006) Assessing genetic diversity in plasmids from Escherichia coli and Salmonella enterica using a mixed-plasmid microarray. J Appl Microbiol 100: 15-28. doi: 10.1111/j.1365-2672.2005.02775.x
[13]	Shubin FN, Kovalchuk NI, Kuznetsova NA, et al. (2002) Microbiological monitoring for Salmonella enteritidis in Primorye Region. Phenotypical and plasmid characterization of the pathogen. Epidemiol Infect Dis 1: 36-40.
[14]	Rakov AV, Shubin FN (2019) Comparative genomics analysis of Salmonella enterica subsp. enterica serotype Enteritidis virulence plasmid. Russ J Genet 55: 144-153. doi: 10.1134/S102279541902011X
[15]	Orskov F, Orskov I (1983) Summary of a workshop on the clone concept in the epidemiology, taxonomy and evolution of the Enterobacteriaceae and other bacteria. J Infect Dis 148: 346-357. doi: 10.1093/infdis/148.2.346
[16]	Selander RK, Caugant DA, Whittman TS (1987) Genetic structure and variation in natural populations of Escherichia coli. Cell Mol Biol 2: 1625-1648.
[17]	Kado CI, Liu ST (1981) Rapid procedure for detection and isolation of large and small plasmids. J Bacteriol 145: 1365-1373. doi: 10.1128/JB.145.3.1365-1373.1981
[18]	Jones GW, Rabert DK, Svinarich DM, et al. (1982) Association of adhesive, invasive, and virulent phenotypes of Salmonella Typhimurium with autonomous 60-megadalton plasmids. Infect Immun 38: 476-486. doi: 10.1128/IAI.38.2.476-486.1982
[19]	Nakamura M, Sato S, Ohya T, et al. (1985) Possible relationship of a 36-megadalton Salmonella Enteritidis plasmid to virulence in mice. Infect Immun 47: 831-833. doi: 10.1128/IAI.47.3.831-833.1985
[20]	Rivera MJ, Rivera N, Castillo J, et al. (1991) Molecular and epidemiological study of Salmonella clinical isolates. J Clin Microbiol 29: 927-932. doi: 10.1128/JCM.29.5.927-932.1991
[21]	Rodrigue DC, Cameron DN, Puhr ND, et al. (1992) Comparison of plasmid profiles, phage types, and antimicrobial resistance patterns of Salmonella Enteritidis isolates in the United States. J Clin Microbiol 30: 854-857. doi: 10.1128/JCM.30.4.854-857.1992
[22]	Holmberg SD, Wachsmuth IK, Hickman-Brenner FW, et al. (1984) Comparison of plasmid profile analysis, phage typing, and antimicrobial susceptibility testing in characterizing Salmonella Typhimurium isolates from outbreaks. J Clin Microbiol 19: 100-104. doi: 10.1128/JCM.19.2.100-104.1984
[23]	Farrar WE (1983) Molecular analysis of plasmids in epidemiologic investigation. J Infect Dis 148: 1-6. doi: 10.1093/infdis/148.1.1
[24]	Southern EM (1975) Detection of specific sequences among DNA fragments separated by gel electrophoresis. J Mol Biol 98: 503-517. doi: 10.1016/S0022-2836(75)80083-0
[25]	Olsen JE, Skov MN, Threlfall EJ, et al. (1994) Clonal lines of Salmonella enterica serotype Enteritidis documented by IS200-, ribo-, pulsed-field gel electrophoresis and RFLP typing. J Med Microbiol 40: 15-22. doi: 10.1099/00222615-40-1-15
[26]	Guerra B, Landeras E, Gonzalez-Hevia MA, et al. (1997) A three-way ribotyping scheme for Salmonella serotype Typhimurium and its usefulness for phylogenetic and epidemiological purposes. J Med Microbiol 46: 307-313. doi: 10.1099/00222615-46-4-307
[27]	Stanley J, Baquar N, Threlfall EJ (1993) Genotypes and phylogenic relationships of Salmonella Typhimurium are defined by molecular fingerprinting of IS200 and 16S rrn loci. J Gen Microbiol 139: 1133-1140. doi: 10.1099/00221287-139-6-1133
[28]	Hulton CS, Higgins CF, Sharp PM (1991) ERIC sequences: a novel family of repetitive elements in the genomes of Escherichia coli, Salmonella Typhimurium and other enterobacteria. Mol Microbiol 5: 825-834. doi: 10.1111/j.1365-2958.1991.tb00755.x
[29]	Hilton AC, Banks JG, Penn CW (1996) Random application of polymorphic DNA (RAPD) of Salmonella: strain differentiation and characterization of amplified sequences. J Appl Bacteriol 81: 575-584.
[30]	Schwartz DC, Cantor CR (1984) Separation of yeast chromosome-sized DNAs by pulsed gel electrophoresis. Cell 37: 67-75. doi: 10.1016/0092-8674(84)90301-5
[31]	Swaminathan B, Barrett TJ, Hunter SB, et al. (2001) PulseNet: the molecular subtyping network for foodborne bacterial disease surveillance, United States. Emerg Infect Dis 7: 382-389. doi: 10.3201/eid0703.017303
[32]	Ridley AM, Threlfall EJ, Rowe B (1998) Genotypic characterization of Salmonella Enteritidis phage types by plasmid analysis, ribotyping, and pulsed-field gel electrophoresis. J Clin Microbiol 36: 2314-2321. doi: 10.1128/JCM.36.8.2314-2321.1998
[33]	Ramisse V, Houssu P, Hernandez E, et al. (2004) Variable number of tandem repeats in Salmonella enterica subsp. enterica for typing purposes. J Clin Microbiol 42: 5722-5730. doi: 10.1128/JCM.42.12.5722-5730.2004
[34]	Boxrud D, Pederson-Gulrud K, Wotton J, et al. (2007) Comparison of multiple-locus variable-number tandem repeat analysis, pulsed-field gel electrophoresis, and phage typing for subtype analysis of Salmonella enterica serotype Enteritidis. J Clin Microbiol 45: 536-543. doi: 10.1128/JCM.01595-06
[35]	Liu Y, Shi X, Li Y, et al. (2016) The evaluation and application of multilocus variable number tandem repeat analysis (MLVA) for the molecular epidemiological study of Salmonella enterica subsp. enterica serovar Enteritidis infection. Ann Clin Microbiol Antimicrob 15: 4. doi: 10.1186/s12941-016-0119-3
[36]	Kotetishvili M, Stine OC, Kreger A, et al. (2002) Multilocus sequence typing for characterization of clinical and environmental Salmonella strains. J Clin Microbiol 40: 1626-1635. doi: 10.1128/JCM.40.5.1626-1635.2002
[37]	Achtman M, Wain J, Weill FX, et al. (2012) Multilocus sequence typing as a replacement for serotyping in Salmonella enterica. PLoS Pathog 8: e1002776. doi: 10.1371/journal.ppat.1002776
[38]	Feil EJ, Li BC, Aanensen DM, et al. (2004) eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol 186: 1518-1530. doi: 10.1128/JB.186.5.1518-1530.2004
[39]	Deng X, Shariat N, Driebe EM, et al. (2015) Comparative analysis of subtyping methods against a whole-genome-sequencing standard for Salmonella enterica serotype Enteritidis. J Clin Microbiol 53: 212-218. doi: 10.1128/JCM.02332-14
[40]	Alikhan NF, Zhou Z, Sergeant MJ, et al. (2018) A genomic overview of the population structure of Salmonella. PLoS Genet 14: e1007261. doi: 10.1371/journal.pgen.1007261
[41]	Rivera MJ, Rivera A, Castillo J, et al. (1993) Plasmid profile in epidemiological studies of human Salmonella infections. J Chemother 5: 288-290.
[42]	Rychlik I, Karpiskova R, Faldynova M, et al. (1998) Computer-assisted restriction endonuclease analysis of plasmid DNA in field strains of Salmonella Enteritidis. Can J Microbiol 44: 1183-1185. doi: 10.1139/w98-112
[43]	Rychlik I, Svestkova A, Karpiskova R (2000) Subdivision of Salmonella enterica serovar enteritidis phage types PT14b and PT21 by plasmid profiling. Vet Microbiol 74: 217-225. doi: 10.1016/S0378-1135(00)00185-1
[44]	Cieslik A, Brown D, Paciorek J, et al. (2001) Phage types, plasmid profiles and chromosomal restriction profiles of Salmonella enterica subsp. enterica ser. Enteritidis (S. Enteritidis) isolated in Poland in 1999–2000. Med Dosw Microbiol 53: 185-196.
[45]	Desai M, Threlfall EJ, Stanley J (2001) Fluorescent amplified-fragment length polymorphism subtyping of the Salmonella enterica serovar Enteritidis phage type 4 clone complex. J Clin Microbiol 39: 201-206. doi: 10.1128/JCM.39.1.201-206.2001
[46]	Nauerby B, Pedersen K, Dietz HH, et al. (2000) Comparison of Danish isolates of Salmonella enterica serovar enteritidis PT9a and PT11 from hedgehogs (Erinaceus europaeus) and humans by plasmid profiling and pulsed-field gel electrophoresis. J Clin Microbiol 38: 3631-3635. doi: 10.1128/JCM.38.10.3631-3635.2000
[47]	Ling JM, Koo IC, Kam KM, et al. (1998) Antimicrobial susceptibilities and molecular epidemiology of Salmonella enterica serotype Enteritidis strains isolated in Hong Kong from 1986 to 1996. J Clin Microbiol 36: 1693-1699. doi: 10.1128/JCM.36.6.1693-1699.1998
[48]	Su LH, Chiu CH, Wu TL, et al. (2002) Molecular epidemiology of Salmonella enterica serovar Enteritidis isolated in Taiwan. Microbiol Immunol 46: 833-840. doi: 10.1111/j.1348-0421.2002.tb02770.x
[49]	Shaginyan IA (2000) Role and significance of molecular methods in epidemiological analysis of nosocomial infections. Clin Microbiol Antimicrob Chemother 3: 82-95.
[50]	Khazenson LB, Poplavskaia ZhV, Kariagina EI, et al. (1996) Epidemiological data on salmonellosis due to Salmonella enteritidis in some areas of the Russian Federation. J Microbiol Epidemiol Immunobiol 4: 53-57.
[51]	Rodrigue DC, Tauxe RV, Rowe B (1990) International increase in Salmonella Enteritidis: a new pandemic? Epidemiol Infect 105: 21-27. doi: 10.1017/S0950268800047609
[52]	Bäumler AJ, Hargis BM, Tsolis RM (2000) Tracing the origins of Salmonella outbreaks. Science 287: 50-52. doi: 10.1126/science.287.5450.50
[53]	Rakov AV, Shubin FN, Ivanis VA, et al. (2001) Comparative characterization of infections caused by Salmonella enteritidis of different plasmid profiles. Epidemiol Infect Dis 5: 50-54.
[54]	Shubin FN, Rakov AV, Kuznetsova NA, et al. (2006) Structure of Salmonella enteritidis population in Primorye Region on the plasmid analysis data. J Microbiol Epidemiol Immunobiol 3: 28-32.
[55]	Shubin FN, Rakov AV, Kuznetsova NA (2011) Microbiological molecular genetic monitoring of enteric infection pathogens as an element of epidemiological surveillance. Bulletin SB RAMS 31: 99-105.
[56]	Rakov AV, Shubin FN, Kuznetsova NA (2013) Heterogeneity of 1.4 MDa plasmids in Salmonella enteritidis strains. Bulletin SB RAMS 33: 10-15.
[57]	Shubin FN, Rakov AV, Kuznetsova NA, et al. (2017) Formation of population morbidity with salmonellosis caused by Salmonella enteritidis in regions with incomplete supply of local poultry products. J Microbiol Epidemiol Immunobiol 1: 61-67. doi: 10.36233/0372-9311-2017-1-61-67
[58]	Shubin FN, Kuznetsova NA, Rakov AV, et al. (2018) Specific features of molecular epidemiology of imported morbidity caused by Salmonella Enteritidis introduced strains with specific plasmid types. Epidemiol Infect Dis 23: 73-76.
[59]	Rakov AV, Kuznetsova NA, Solovyeva AS, et al. (2018) Cluster analysis of Salmonella Enteritidis isolated in the Siberia and Far East of Russia. Pacific Med J 4: 23-26. doi: 10.17238/PmJ1609-1175.2018.4.23-26
[60]	Rakov AV, Shubin FN, Kuznetsova NA, et al. (2019) Heterogeneity of 2.3 MDa plasmids in Salmonella Enteritidis strains. Siberian Sci Med J 39: 40-45.
[61]	Kuznetsova NA, Solovyeva AS, Rakov AV (2018) Antibiotic resistance of Salmonella Enteritidis strains, circulated in territory of the Siberia and Far East, at multi-year monitoring. Health Med Ecol Sci 3: 50-58.
[62]	Rakov AV, Shubin FN, Kuznetsova NA, et al. (2016) Estimated microbiological and epidemiological of the outbreak of salmonellosis in the child care institution in Murmansk. Health Med Ecol Sci 4: 26-30. doi: 10.4103/2321-4848.183351
[63]	Döpfer D, Buist W, Soyer Y, et al. (2008) Assessing genetic heterogeneity within bacterial species isolated from gastrointestinal and environmental samples: how many isolates does it take? Appl Environ Microbiol 74: 3490-3496. doi: 10.1128/AEM.02789-07
[64]	Timchenko NF, Rakov AV, Terentyeva NA, et al. (2019) Characteristics of the mixed bacteria of the Enterobacteriaceae family Yersinia pseudotuberculosis and Salmonella Enteritidis in vitro. Health Med Ecol Sci 1: 19-22.
[65]	Yakovlev AA, Pozdeeva ES (2018) Possible mechanisms of self-regulation of parasitic systems in the biogeocenosis. Vestnik RAMS 73: 195-205.

This article has been cited by:

1.	Rania Saadeh, A Reliable Algorithm for Solving System of Multi-Pantograph Equations, 2022, 21, 2224-2880, 792, 10.37394/23206.2022.21.91
2.	Ahmad Qazza, Rania Saadeh, Said El Kafhali, On the Analytical Solution of Fractional SIR Epidemic Model, 2023, 2023, 1687-9732, 1, 10.1155/2023/6973734
3.	Manoj Singh, Ahmed Hussein, Mohammad Tamsir, Abdullah Ali H. Ahmadini, An analytical approach of multi-dimensional Navier-Stokes equation in the framework of natural transform, 2024, 9, 2473-6988, 8776, 10.3934/math.2024426
4.	Ahmad Qazza, Rania Saadeh, Osama Ala'yed, Ahmad El-Ajou, Effective transform-expansions algorithm for solving non-linear fractional multi-pantograph system, 2023, 8, 2473-6988, 19950, 10.3934/math.20231017
5.	Ahmad El-Ajou, Mohammed Shqair, Ibrahim Ghabar, Aliaa Burqan, Rania Saadeh, A solution for the neutron diffusion equation in the spherical and hemispherical reactors using the residual power series, 2023, 11, 2296-424X, 10.3389/fphy.2023.1229142
6.	AHMAD EL-AJOU, RANIA SAADEH, ALIAA BURQAN, MAHMOUD ABDEL-ATY, A MODERN TRAVELING WAVE SOLUTION FOR CAPUTO-FRACTIONAL KLEIN–GORDON EQUATIONS, 2024, 32, 0218-348X, 10.1142/S0218348X24500841
7.	Nada Sabeeh Mohammed, Emad A. Kuffi, Rania Saadeh, 2023, Implementation of “Emad-Falih Transform” to Resolve Electrical Circuits, 979-8-3503-7336-3, 1, 10.1109/EICEEAI60672.2023.10590438
8.	Rania Saadeh, Mohamed A. Abdoon, Ahmad Qazza, Mohammed Berir, A Numerical Solution of Generalized Caputo Fractional Initial Value Problems, 2023, 7, 2504-3110, 332, 10.3390/fractalfract7040332
9.	Amjad E. Hamza, Abdelilah K. Sedeeg, Rania Saadeh, Ahmad Qazza, Raed Khalil, A New Approach in Solving Regular and Singular Conformable Fractional Coupled Burger’s Equations, 2023, 22, 2224-2880, 298, 10.37394/23206.2023.22.36
10.	Osama Ala'yed, Rania Saadeh, Ahmad Qazza, Numerical solution for the system of Lane-Emden type equations using cubic B-spline method arising in engineering, 2023, 8, 2473-6988, 14747, 10.3934/math.2023754
11.	Mohamed A. Abdoon, Rania Saadeh, Mohammed Berir, Fathelrhman EL Guma, Mawada ali, Analysis, modeling and simulation of a fractional-order influenza model, 2023, 74, 11100168, 231, 10.1016/j.aej.2023.05.011
12.	Aliaa Burqan, Rania Saadeh, Ahmad Qazza, Ahmad El-Ajou, Trustworthy Analytical Technique for Generating Multiple Solutions to Fractional Boundary Value Problems, 2023, 9, 2349-5103, 10.1007/s40819-023-01554-y
13.	Faeza Hasan, Mohamed Abdoon, Rania Saadeh, Mohammed Berir, Ahmad Qazza, A New Perspective on the Stochastic Fractional Order Materialized by the Exact Solutions of Allen-Cahn Equation, 2023, 8, 2455-7749, 912, 10.33889/IJMEMS.2023.8.5.052
14.	Shuguang Li, Rania Saadeh, J.K. Madhukesh, Umair Khan, G.K. Ramesh, Aurang Zaib, B.C. Prasannakumara, Raman Kumar, Anuar Ishak, El-Sayed M. Sherif, Aspects of an induced magnetic field utilization for heat and mass transfer ferromagnetic hybrid nanofluid flow driven by pollutant concentration, 2024, 53, 2214157X, 103892, 10.1016/j.csite.2023.103892
15.	P. Dunnimit, W. Sawangtong, P. Sawangtong, An approximate analytical solution of the time-fractional Navier–Stokes equations by the generalized Laplace residual power series method, 2024, 9, 26668181, 100629, 10.1016/j.padiff.2024.100629
16.	Ahmad El-Ajou, Rania Saadeh, Moawaih Akhu Dunia, Ahmad Qazza, Zeyad Al-Zhour, A new approach in handling one-dimensional time-fractional Schrödinger equations, 2024, 9, 2473-6988, 10536, 10.3934/math.2024515
17.	Ahmad El-Ajou, Haneen Al-ghananeem, Rania Saadeh, Ahmad Qazza, Moa’ath N. Oqielat, A modern analytic method to solve singular and non-singular linear and non-linear differential equations, 2023, 11, 2296-424X, 10.3389/fphy.2023.1167797
18.	Mubashir Qayyum, Efaza Ahmad, Syed Tauseef Saeed, Ali Akgül, Sayed M. El Din, New solutions of fractional 4D chaotic financial model with optimal control via He-Laplace algorithm, 2024, 15, 20904479, 102503, 10.1016/j.asej.2023.102503
19.	Jayaram Prakash K, Saratha S R, Ahmet Yildirim, A novel approach to fractional differential equations via JSN transform coupled fractional residual power series method, 2024, 99, 0031-8949, 125285, 10.1088/1402-4896/ad92bb

Reader Comments

Your name:*

Email:*
© 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Microbiology

2.7 7.0

Metrics

Article views(4846) PDF downloads(439) Cited by(6)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(3)

AIMS Microbiology

Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles

Related Papers:

Abstract

1. Introduction

2. Basic concepts

3. The Laplace residual power series method

4. Numerical application

5. Mathematical reviews

5.1. Advantages of the method

5.2. Disadvantages of the method

6. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Microbiology

Genetic diversity of Salmonella enterica subsp. enterica serovar Enteritidis in the Siberia and Far East of Russia based on plasmid profiles

Related Papers:

Abstract

1. Introduction

2. Basic concepts

3. The Laplace residual power series method

4. Numerical application

5. Mathematical reviews

5.1. Advantages of the method

5.2. Disadvantages of the method

6. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog