Modelling and analysis of HFMD with the effects of vaccination, contaminated environments and quarantine in mainland China

1.   Introduction

Stagnation point flow reveals in practically all flow branches of engineering and science. Free stagnation point or a line exists interior of the fluid domain while in some situations flow is stagnated by a solid wall.

Hiemenz ^[1] demonstrated that Navier-Stokes equations are able to describe the stagnation point flow. Plane stagnation point flow which is called the Hiemenz flow, firstly discussed by Hiemenz as follows:

with the boundaries:

Gersten et al. ^[2] generalized these results to the stagnation point flow in three dimensions. Homotopy analysis method is used by Liao ^[3] to the 2D laminar viscous flow over a semi-infinite plate. Existence and uniqueness consequences for nonlinear third-order ODEs appeared in the stagnation point flow of a hydromagnetic fluid are discussed by Van Gorder et al. ^[4]. The flow characteristics in an electrically conducting second grade fluid pervaded by a uniform magnetic field over a stretching sheet is discussed by Vajravelu et al. ^[5]. Steady 2D stagnation-point flow of an incompressible viscous electrically conducting fluid over a flat deformable sheet is studied by Mahapatra et al. ^[6]. Van Gorder et al. ^[7] established the existence and uniqueness consequences for fourth-order ODEs over the semi-infinite interval, arising in the hydromagnetic stagnation point flow of a second grade fluid over a stretching sheet.

The concept of Lie groups plays an significant rule in establishing some intense numerical approaches to integrate ODEs, which conserve the geometrical properties. Preserving the geometry construction under discretization, is essential in the improvement of a qualitatively true performance in the minimization of approximation errors. So, along with the geometric structure and invariance of the ODEs, novel approaches may be invented, which are more stable and effective than the common approximation techniques.

A novel geometric method based on $GL(6, \mathbb{R})$ is introduced for the magneto-hemodynamic flow in a semi-porous channel in ^[8]. Moreover, in ^[9], a nonlinear heat transfer equation is considered by $SL(2, \mathbb{R})$. Order of differential equations and different initial or boundary conditions conclude different geometric algorithms for solving differential equations. That is, numerical algorithms based on e.g. $GL(6, \mathbb{R})$ and $GL(4, \mathbb{R})$ is completely different. Combination of radial basis functions with GPS are discussed in ^[10]. Some other recent papers of geometric numerical approach, Lie group analysis, and exact solution can be found in ^{[11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33]}.

This paper is organized as follows. In the further section, we will touch only a few aspects of the Mathematical formulation for the hydromagnetic stagnation point flow. In Section 3, we derive an interesting algorithm to integrate the underlying nonlinear equation. Section 4 deals with the numerical results. Final remarks are provided in the conclusion section.

2.   Mathematical formulation

Let us investigate the 2D steady stagnation-point flow of an electrically conducting fluid in the existence of a uniform transverse magnetic field towards a flat surface coinciding with the plane $y = 0,$ the flow being enclosed to the region $y > 0.$ Two equivalent and converse forces are considered in direction of the $x-$axis such that the wall is stretched keeping the origin fixed, and a uniform magnetic field $B_0$ is forced along the $y-$axis (Figure 1).

In the usual notation, the MHD equations for steady 2D stagnation-point flow in the boundary layer along the stretching surface are ^[6,34]:

Figure 1.  A profile of the physical problem.

DownLoad: Full-Size Img PowerPoint

where $\rho, \ \sigma, \ B_0$ are density, electrical conductivity and magnetic field respectively. Also, $(u, v)$ are the fluid velocity components toward $x$ and $y$. $U(x)$ stands for the stagnation-point velocity in the inviscid free stream.

For this problem, the boundaries are

where $c, a\in \mathbb{R}^+$. Eqs (2.1) and (2.2) together with the conditions (2.3) and (2.4) have the similarities of the form

Substituting (2.5) into (2.2), concludes

where $M = B_0^2\left(\frac{\sigma}{\rho c}\right), \ C = \frac{a}{c}$ and $\sqrt{M}$ is the Hartmann number.

With the use of the transformation (2.5), the boundary conditions (2.3) and (2.4) read $f(0) = 0, \ f'(0) = 1, \ f'(\infty) = C.$

The self-similar solutions will normally be made to meet Eq (2.6) with additional nonlinear terms

in the governing differential equation, where $K\geq 0$ is the viscoelastic parameter. By the additional condition $f''(\infty) = 0$, we investigate the following fourth order differential equation:

with

where $M = \sigma_0 B_0^2/\rho B\geq 0$ is the magnetic parameter, $K = \omega B/\nu\geq 0$ is the viscoelastic parameter, $s\geq 0, \ C$ is the parameter due to the stagnation point flow, and $\chi$ is the stretching of the sheet.

3.   An $GL(4, \mathbb{R})$ Lie group method

There is no loss of generality in assuming

This supposition is possible. Because in the other case, we can obtain a constant $\gamma > 0$ such that

where

and $c_1$-$c_3$ are integration constants which we can assume all of them to be positive. Also, from this assumption we clearly have $|F(\xi)| > 0$. Using these transformations, Eq (2.8) converts into:

subject to

Now, we are willing to erect $GL(4, \mathbb{R})$ Lie group method. Suppose $y_1(\xi) = f(\xi), \ y_2(\xi) = f'(\xi), \ y_3(\xi) = f''(\xi)$ and $y_4(\xi) = f'''(\xi)$. Then equivalent system of (2.8) leaves

with boundary conditions:

Equivalent system of (3.7) is:

where

which is well-defined because of supposition (3.1). System (3.9), in spite of its nonlinearity, admits a Lie-group of $GL(4, \mathbb{R})$ and we have:

where

Now, we are able to erect a GPS to consider Eq (3.9) as follows:

where $\mathbf{Y}_n = \mathbf{Y}|_{\xi = \xi_n}$. Integrating (3.9) by (3.13) along with initial condition $\mathbf{Y}(0) = \mathbf{Y}_0$, results the demanded value $\mathbf{Y}(\xi)$. Suppose the GPS method:

where

and $\Delta\xi = 1/N$ is the GPS iteration stepsize. In fact, this approach approximates the solution of

Then

computes $\mathbf{Y}(\xi)|_{\xi = \eta_\infty}$, ^* at $\xi = \eta_\infty.$ From the closure property of the Lie groups, $\mathbb{G}: = \mathbb{G}_N(\Delta\xi)\cdots \mathbb{G}_1(\Delta\xi)\in GL(4, \mathbb{R})$. Therefore, we can develop a one-step transformation from $\mathbf{Y}_0$ to $\mathbf{Y}_f$ as:

^* By using the truncation rule for the differential equations over the infinite or semi-infinite intervals, we set $\xi^f = \eta_\infty$ as the final value of interval $[0, \infty)$.

From Eqs (3.11) and (3.18) we obtain

Along with (3.8) and generalized mid-point rule, let us assume

where $r\in[0, 1]$ has to be determined. Therefore Eq (3.19) may be rewritten as:

where

In this letter, $y_1^0 = s, \ y_2^0 = \chi, \ y_2^f = C, \ y_3^f = 0$ are known and $y_1^f, \ y_3^0, \ y_4^0, \ y_4^f$ are unknown values of the current problem, which after determining $y_3^0$ and $y_4^0$, the B.V.P. (2.8) converts to an I.V.P.

For the fixed $\mathcal{A}\in gl(4, \mathbb{R})$, ^† we have

^† $gl(4, \mathbb{R})$ is the Lie algebra associated with the Lie group $GL(4, \mathbb{R}).$

where

Therefore, each of this Lie group elements may be erected by a single parameter $\mathbb{G}(r),$ with $r\in [0, 1]$.

From (3.18) and (3.22) we have:

Thus, for a given $r$, we may get the values $y_1^f, \ y_3^0, \ y_4^0$ and $y_4^f$ by the further explained process.

Imposing the initial guesses of $y_3^0, \ y_4^0, \ y_1^f$ and $y_4^f$, and the given boundary values of $y_1^0 = s, \ y_2^0 = \chi, \ y_2^f = C$ and $y_3^0 = 0$, with a specified $r\in[0, 1]$, and from

we may produce new values of $y_1^f, \ y_3^0, \ y_4^0$ and $y_4^f$. If these values satisfy the convergence criterion:

then iterations stop.

For a fixed value of $r$, we are able to compute $y_3^0$ and $y_4^0$ from Eq (3.25) and then approximately integrate Eq (3.7) by GPS from $0$ to $\eta_\infty$, and compare the end value of $y_2^f$ and $y_3^f$ with the exact ones $y_2(\eta_\infty) = C$ and $y_3(\eta_\infty) = 0$ respectively. That is, we have to find the minimum value of

In the appendix, we give a brief algorithm of our proposed method to find the initial values of problem.

In the non-geometric numerical methods, investigation of the method stability is an important task which guarantees the validity of techniques. Without considering the convergence and error analysis, results may lead to solutions with no physical meaning or high-level of errors. This issue may arise from that, we do not know information about the geometric structure of original problem. However, in this paper we introduce a $GL(4, \mathbb{R})$ method which conserve the geometrical properties of problem. Preserving the geometry construction under discretization, is essential in the improvement of a qualitatively true performance in the minimization of approximation errors.

4.   Numerical results

We emphasize that the stepping approach developed for integrating the IVP needs the initial values $f_1(0) = f(0), \ f_2(0) = f'(0), \ f_3(0) = f''(0)$ and $f_4(0) = f'''(0)$ for the fourth-order ODEs. Both of the last initial values in Eq (2.8) are unknown. Therefore, approximating $f''(0)$ and $f'''(0)$ converts Eq (2.8) into an IVP. It is worth pointing out that every IVP can be integrated by GPS.

Indeed, our proposed method in this paper is split into two parts. In the first part, initial conditions of the problem are obtain by the $GL(4, \mathbb{R})$ method. Our method is implemented by the convergence criterion of (3.26) with $\epsilon = 10^{-8}$. For this value of convergence criterion, only four or five steps of algorithm, mentioned in the Appendix, are needed. So, very fast convergence occur in this part of the proposed method Second part of this technique starts when we obtain the initial values of BVP from the first part. Second part which is nominated as the GPS, is an explicit iteration method which is well-known as the weighted Euler method. Convergence of this part can be guaranteed from the geometric structure of method.

The results of this paper were announced with $\eta_\infty = 7$ and $\Delta\xi = 0.01$ for the GPS procedure. There is no loss of generality in assuming $y_3^0 = y_4^0 = y_1^f = y_4^f = 1$ as initial guesses.

Some mis-matching errors with various $M, \ C, \ s$, $\ K$ and $\chi$ are plotted in Figure 2(a)-(d). Desired values of $r = 0.5314, \ r = 0.4993, \ r = 0.814$ and $r = 0.9546$ are noticeable from these figures.

Figure 2.  Plots of mis-matching errors for some values of $M, C, s, K$ and $\chi$.

DownLoad: Full-Size Img PowerPoint

In Figure 3(a)-(d), we plot the profiles of $f(\xi), \ f'(\xi), \ f''(\xi)$ and $f'''(\xi)$ for various values of magnetic parameter (the Hartmann number) $M\ (\sqrt{M})$, and fixed values of the other parameters. It is observable that the profiles of $f(\xi)$ and $f'(\xi)$ decrease, when the magnetic parameter increases, while the reverse situation occur for $f''(\xi)$ and it decrease more rapidly toward zero as $\xi$ tends to $\eta_\infty.$

Figure 3.  Profiles of $(a):f(\xi), \ (b):f'(\xi), \ (c):f''(\xi)$ and $(d):f'''(\xi)$, for various $M$ when $C = s = K = 1$ and $\chi = 3$.

DownLoad: Full-Size Img PowerPoint

Figure 4(a)-(d) show the absorbing dependence of the solutions on the various viscoelastic parameter $K.$ Uniform decrease in the magnitudes of $f(\xi)$ and $f'(\xi)$ are observable for increasing values in $K.$ Inverse results occur for the profiles of $f(\xi)$ and $f'(\xi)$ when $\chi$ is greater than the boundary value $C$; see Figure 5(a)-(d). Notice that the plots are in a good agreement with the boundary conditions $f'(\eta_\infty) = C$ and $f''(\eta_\infty) = 0.$

Figure 4.  Variation of $(a):f(\xi), \ (b):f'(\xi), \ (c):f''(\xi)$ and $(d):f'''(\xi)$, for various viscoelastic parameter $K$ when $C = s = M = 1,$ and $\chi = 0$.

DownLoad: Full-Size Img PowerPoint

Figure 5.  Profiles of $(a):f(\xi), \ (b):f'(\xi), \ (c):f''(\xi)$ and $(d):f'''(\xi)$, for various viscoelastic parameter $K$ when $C = s = M = 1,$ and $\chi = 3$.

DownLoad: Full-Size Img PowerPoint

It may be seen from Figure 6(a)-(d) that for fixed values of $K, s, M$ and $\chi,$ the transverse velocity increases with increase in the stagnation point flow $C$. Notice the behavior of solutions in the profiles of $f(\xi), \ f'(\xi)$ and $f''(\xi)$ for different cases $\chi < C$ and $\chi > C$ which are against to the profile of $f'''(\xi)$.

Figure 6.  Variation of $(a):f(\xi), \ (b):f'(\xi), \ (c):f''(\xi)$ and $(d):f'''(\xi)$, for various stagnation point flow $C$ when $K = s = M = 1,$ and $\chi = 1$.

DownLoad: Full-Size Img PowerPoint

From the Figure 7(a)-(d), an increase in the stretching sheet parameter $\chi$ results in a increase in the profiles of $f(\xi), \ f'(\xi)$ and $f'''(\xi)$. Also, an increase in $\chi$ consequences in a decrease in the magnitude of the shear stress at the wall, $f''(0).$

Figure 7.  Profiles of $(a):f(\xi), \ (b):f'(\xi), \ (c):f''(\xi)$ and $(d):f'''(\xi)$, for various stretching sheet parameter $\chi$ when $K = s = M = 1,$ and $C = 1$.

DownLoad: Full-Size Img PowerPoint

The double mesh principle is applied to estimate the error

where $N$ is the number of GPS iterations. Residual errors of the present method applied for Eq (2.8) with $N = 10000$ are plotted in Figure 8 for various parameters of equation. This figure shows that our method is reliable for current types of equations. Besides, obtained solutions and results are in good agreement with the results of ^[7].

Figure 8.  Plot of estimate errors of $f(\xi)$ for some values of $M, C, K, s$ and $\chi$.

DownLoad: Full-Size Img PowerPoint

5.   Conclusions

In this paper, a novel iterative technique based on the Lie group of $GL(4, \mathbb{R})$, as the first time, is proposed to solve the hydromagnetic stagnation point flow of a second grade fluid over a stretching sheet. Structures of the Lie group elements in $GL(4, \mathbb{R})$ and corresponding Lie algebra elements in $gl(4, \mathbb{R})$ are completely innovative for the underlying problem in this paper. Against the non-geometric numerical techniques, the proposed method uses the geometric structure of original problem and this guarantees that obtained results are reliable and the proposed method avoid the ghost solutions. Effects of equation parameters are graphically considered. Point-wise errors for various parameters show the power and effectiveness of proposed method.

Conflict of interest

No potential competing interest was reported by the authors.

Appendix

Algorithm of finding initial values of the problem:

(i) Take arbitrary initial guesses for $y_3^0, \ y_4^0, \ y_1^f$ and $y_4^f$.

(ii) Specify an $r\in[0, 1]$ and a small $\epsilon\in\mathbb{R}^+$.

(iii) Calculate $\mathbb{G}(r)$ and $\omega$ from (3.22) and (3.23), respectively.

(iv) Compute updated values of $y_3^0, \ y_4^0, \ y_1^f$ and $y_4^f$ from (3.25).

(v) If convergence criterion (3.26) is valid for given $\epsilon$, then $y_3^0$ and $y_4^0$ are desired initial values of the problem. Otherwise, go to Step (ii).