On a generalized fractional boundary value problem based on the thermostat model and its numerical solutions via Bernstein polynomials

In this paper, we introduce a new structure of the generalized multi-point thermostat control model motivated by its standard model. By presenting integral solution of this boundary problem, the existence property along with the uniqueness property are investigated by means of a special version of contractions named μ-φ-contractions and the Banach contraction principle. Then, on the given nonlinear generalized BVP of thermostat, the Bernstein polynomials are introduced and numerical solutions obtained by them are presented. At the end, three different structures of nonlinear thermostat models are designed and the results are examined.


Introduction
Fractional calculus and the existing notions in it are of high interest in different aspects of applied sciences, and one can find some instances of applications like signal and image processing, control theory, economics, optical systems, thermal materials, aerodynamics, mechanical systems, biology, and bio-mathematics [1][2][3][4][5][6][7][8][9][10]. The starting point for such a topic can be seen in many published papers in which mathematicians deal with some properties such as the existence of solution, uniqueness property, the property of stability, positivity, etc., and establish these properties to various abstract boundary value problems. Such a diversity and importance led to the publication of many research papers in this field, which revealed the flexibility of fractional calculus theory in designing various mathematical models. The main methods conducted in these articles are by terms of fixed point techniques .
Along with the investigation of existence theory, the approximation of solution and numerical methods are also of interest. Al-Smadi et al. [34] utilized a method based on the homotopy analysis for finding approximate solutions of a SEIR epidemic model in the fractional settings. Dhage, Dhage, and Ntouyas in [35] applied two notions of partial compactness and continuity to develop the method of Kranoselskii theorem and found the approximation of solutions regarding a hybrid ODE. Chadha et al. published a paper on the Faedo-Galerkin approximation for solutions of a neutral nonlocal FDE in a separable Hilbert space [36]. Other applications of approximate techniques can be found in [37][38][39][40]. The Bernstein polynomials are one of the strongest numerical techniques which possess some important properties such as the unity partition and continuity on [0, 1] [41]. In recent years, due to the importance and accuracy of this technique, the numerical solutions of a wide range of linear\nonlinear BVPs have been obtained for Riccati type FDEs, Bessel FDE, Lane-Emden equations, etc., [11][12][13][42][43][44]. We here use these polynomials to find approximate solution of our given multi-point FBVP introduced in the sequel.
In 2006, two mathematicians, Infante and Webb, simulated a mathematical model of a mechanical instrument named thermostat in the form of a second-order boundary problem on the interval [0, 1] which is insulated at s = 0 via the controller for s = 1 [45]. This model has the following mathematical structure: where ζ ∈ J and the parameter μ > 0. Further, two other mathematicians, Nieto and Pimentel, converted the above problem to a similar version of arbitrary order [46] which takes the following structure: where c D p is the Caputo derivative and ζ ∈ J. In fact, at the time s = ζ and based on the existing temperature, the sensor detects that the thermostat discharges or adds heat. In this manuscript, we concentrate on this aim in which some existence and uniqueness aspects and numerical analysis of solutions for a generalized fractional boundary value problem (GFBVP) based on thermostat model are investigated. Indeed, we formulate the following structure of a generalized multi-point thermostat control model motivated by the standard model (1): in which p ∈ (1, 2], σ ∈ (0, 1), p -1 ∈ (0, 1], k, β, ρ > 0, 0 < ζ 1 < ζ 2 < · · · < ζ m < 1, m ∈ N, ε 1 , ε 2 ∈ R, and c D 1 = d ds . Along with these, the mapping g : J × (R ≥0 ) 3 → R ≥0 is continuous and c D and I display the derivation and integration operators of order ∈ {1, p, p -1, σ , ρ} in the sense of Caputo and Riemann-Liouville.
The novelty, motivation, and objective of this research are: -The multi-point nonlinear system (3) is a generalized form of the mathematical model of thermostat that by assuming p = 2, μ = k > 0, m = 1, ζ 1 = · · · = ζ m = ζ ∈ (0, 1), ε 1 = 0, ε 2 = 0, and β = 1, we obtain the second-order integro-differential FBVP (1); -The existence property of solutions to the generalized nonlinear GFBVP (3) is derived by terms of a special case of contractions entitled μϕ-contraction and μ-admissible maps; -The approximate solution of this generalized nonlinear GFBVP of thermostat is obtained via Bernstein polynomials; -The accuracy and absolute errors of the mentioned numerical technique are examined in different examples of thermostat model. The rest of the manuscript is provided as follows: The primitive notions are indicated and recalled in Sect. 2. The existence property along with the uniqueness property are investigated in Sect. 3. Bernstein polynomials and numerical solutions obtained by them are presented in Sect. 4. In Sect. 5, three different structures of nonlinear thermostat GFBVP are designed, and the results are examined. We end the manuscript by giving conclusions in Sect. 6.

Basic notions
In this section, we provide some general basic tools and results of fractional calculus that allow us to achieve our desired results. For more details on this subject, we advise the authors to consult, for example, [47,48].

Definition 2.1
The Riemann-Liouville fractional integral (FRL-integral) of order ν > 0 for a continuous function f : R ≥0 → R is defined by such that integral (4) converges.

The existence property
It is known that X = {u : u, c D σ u ∈ C R (J)} is a Banach space with the sup norm u X := sup s∈J |u(s)| + sup s∈J | c D σ u(s)|, where C R (J) denotes the collection of all continuous realvalued functions on J. In the following, we characterize the structure of the solutions for given GFBVP caused by thermostat model (3) which plays a key role in our required method. Before it, we introduce some notations for simplicity: Proposition 3.1 Let p ∈ (1, 2], p -1 ∈ (0, 1], k > 0, 0 < ζ 1 < ζ 2 < · · · < ζ m < 1, m ∈ N, ε 1 , ε 2 ∈ R, and h ∈ C R (J). Then the solution of the linear thermostat GFBVP is given by where A, G ∈ C R (J) are introduced as Proof We assume that u * satisfies the linear thermostat GFBVP (7). Then c D p u * (s) = h(s). By integrating of order 1 < p ≤ 2 on both sides of it, we get where we try to obtain the constant values of the coefficients c 0 , c 1 ∈ R. On the other hand, we have and Now, in view of notations (6) and by using boundary conditions (7) and by invoking relations (11), (12), and (13), we reach and Eventually, by (14) and (15), we substitute the obtained values for the coefficients c 0 and c 1 in (10) and it becomes which confirms that u * satisfies (8), and accordingly, the proof is finished.
To follow the procedure of the paper, we introduce the operator K : X → X associated with the nonlinear thermostat GFBVP which takes the form where the functions A, G ∈ C R (J) are introduced by (9). Before presenting our main theorems, we equip the space X with the metric d formulated as d(x, y) = xy X .
It is well known that (X, d) is a complete metric space (see [51]).
Theorem 3.2 Consider a continuous function g : J × R 3 → R and assume that the following assumptions hold: (ASS1) There are a map ϕ ∈ ( is the family defined in Section 2) and a function w : where ϑ 1 and ϑ 2 are two positive real constants which satisfy the following inequalities: and (ASS4) For each sequence x n ∈ X which converges to x in X and w(x n (s), x n+1 (s)) ≥ 0, ∀s ∈ J and ∀n ∈ N, we have w(x n (s), x(s)) ≥ 0.

Then problem (3) has a solution.
Proof Let us define a map μ : For all x, y ∈ X and w(x(s), y(s)) ≥ 0 for each s ∈ J, we have and Therefore, from (17) and (18) This means that d(Kx, Ky) ≤ ϕ (d(x, y)). Consequently, from the definition of the map μ it follows that d(x, y) , ∀x, y ∈ X, which means that K is a μϕ-contraction. Furthermore, in view of the definition of the map μ and assumption (ASS3), we can easily verify that K is μ-admissible.
Let now x n be a sequence in X which approaches to x in X and satisfies μ(x n , x n+1 ) ≥ 1, ∀n ∈ N and ω(x n (s), x n+1 (s)) ≥ 0. Then from the definition of the map μ together with assumption (ASS5), we can directly verify that μ(x n , x) ≥ 1. At this time, all the assumptions of Theorem 2.6 are fulfilled. Consequently, the operator K admits a fixed point which is solution of our nonlinear thermostat GFBVP (3). (ASS6) ∃ R > 0, s.t. ∀s ∈ J, and ∀x, x, y, y, z, z ∈ R, The constants β and ρ are linked by the relation β + 1 (ρ+1) < 1 and Then the nonlinear thermostat GFBVP (3) has exactly one solution.
Proof By following the same arguments of the calculations used in Theorem 3.2 together with the hypotheses of Theorem 3.3, we write Kx(s) -Ky(s) Similarly, we obtain Thus, It yields Kx -Ky X ≤ γ xy X , ∀x, y ∈ X. Now, from the assumption γ < 1 and the Banach contraction principle, we conclude that K admits a unique fixed point which represents the unique solution of our nonlinear thermostat GFBVP (3).

Numerical solutions via Bernstein polynomials
In the first place, we describe the basic formulation of Bernstein polynomials which are necessary to derive our developed results.
A recursive expression also can be used to formulate the Bernstein basis polynomials on the interval [0, 1] such that the Bernstein polynomials of (j, m)th degree can be rewritten as We can easily show that each of Bernstein basis polynomials is positive and also the sum of all Bernstein basis polynomials is equal to unity for any real z belonging to the interval In [52], Juttler has explicitly represented d j,m (z) by the following formula: where for 0 ≤ j, k ≤ m, Furthermore, where I ν denotes the ν th -FRL-integral and I (ν) represents the (m + 1) × (m + 1) operational matrix of the ν th -FRL-integral. Also, by using the binomial expansion of (1z) m-j , we can write As we can find in [52], an approach for the direct least squares approximation with the help of Bernstein polynomials is based on the construction of the basis {d 0,m (z), d 1,m (z), . . . , d m,m (z)} which represents the dual in Bernstein basis of mth-degree on [0, 1]. It is specified as where δ jk denotes the Kronecker symbol.

Fractional matrix of integration
The conclusion of the next theorem is useful for us.

Theorem 4.2 Let B(z) be the Bernstein vector introduced in (23) and I (ν) be the (m + 1) × (m + 1) operational matrix of the ν th -FRL-integral which is formulated by
Then we have where and λ jk are the coefficients expressed by (22).
Proof From expression (23) and by exploiting relationship (25) which connects Bernstein polynomials with their dual, we can write where D denotes the vector of dual polynomials of Bernstein polynomials and I ν B(z), D T stands for the (m + 1) × (m + 1)-matrix formulated as Here, we have By applying (24), we get Afterwards, in view of (21) we conclude that but we have Therefore, a combination of (28) and (29) ends the proof of our Theorem 4.2.

Fractional matrix of derivative
We can write the derivative of B(z) as where D (1) stands for the (m + 1) × (m + 1)-operational matrix of derivative given in the following format: For more details, we refer to [43,[53][54][55].
By applying (30), it is obvious that for each n ∈ N we have d n B(z) dz n = D (1) n B(z).
For ν > 0, the ν th -Caputo derivative of B(z) is given as in which c D (ν) stands for the (m + 1) × (m + 1)-operational matrix of the ν th -Caputo derivative which is given by where λ pk μ kj , λ pk is defined as (22) and

Some simulative examples
Before illustrating our theoretical results by some numerical examples, we present, in general, the principle of the Bernstein collocation method applied to our problem (3) to obtain an accurate numerical solution. For this fact, let us consider for all s ∈ J with the following conditions: Now, to determine an approximation of the exact solution u(s) by Bernstein polynomials, we utilize the FRL-integral of Bernstein polynomials with the operational matrix of the Caputo derivative used in [53]. We know that the approximate solution of u(s) by Bernstein polynomials is defined by such that C is an indeterminate vector. By replacing the approximate solution given by (36) in (34) and (35), we can write respectively and For convenience of computations, we consider equidistant points and the roots of the Legendre polynomial of degree (m-1) in [0, 1]. So, to obtain the solution u(s), we collocate equation (37) at (m -1) points together with equations (38) and (39). Therefore, we get (m + 1) equations with (m + 1) indeterminate coefficients. Consequently, the approximate solution can be determined.
At the moment, we are ready to illustrate the Bernstein collocation method with some simulative examples.

Conclusions
In this work, we introduced a new generalized version of the mathematical model of the thermostat in the form of the nonlinear GFBVP given as (3). The existence property for its solutions was established via a special form of contractions and μ-admissible maps. The uniqueness property was verified by the Banach principle. Further, we used the Bernstein operational matrix of the Caputo fractional derivative and the Bernstein operational matrix of FRL-integral which are necessary to obtain accurate numerical solutions to the nonlinear thermostat GFBVP (3) via Bernstein polynomials. We have designed two examples to illustrate the accuracy of the numerical method in finding the exact and approximate solutions. Then we checked the uniqueness property in the third example. For the next works, we will apply these methods on different mathematical models designed by nonsingular operators.