The aim of the study is the question, whether printed electronics circuitry from low-cost printers can be improved in quality by certain ink combinations or after-treatments in order to achieve acceptable results in comparison with circuitry from higher-quality printers. For this purpose, the six samples circuitries from different printers (professional and semi-professional) and different ink combinations (PEDOT: PSS, Silver and Carbon) were subjected to an endurance test of 5 million switching cycles under varying climatic conditions. For this purpose test number of N = 5 experimental evaluations were carried out for each samples (in total of 30 experimental evaluations were done). The results show that, respectable results could be achieved with corresponding ink and post-treatment combinations. This opens up new possibilities for future developments in the field of printers, inks and post-treatments under the aspect of “low-cost”.
1.
Introduction
An important tool for describing the natural process is a system of fractional differential equations [1,2]. An ordinary system of fractional differential equations (FDEs) is described by
where a∈R, q:[0,a]→Rν is a ν-dimensional vector function and p:[0,a]×Rν→Rν. Here, C0Dαt is the Caputo fractional derivative of order α>0.
FDEs of the form (1.1) are used in an extensive amount of published papers for describing diverse evolutionary processes. Pollution model for lake [3], dynamical models of happiness [4] and Irving-Mullineux oscillaton [5] are among them.
A system of type (1.1) has at least n=⌈α⌉ (ceil of α) degree of freedom due to the presence of n the integer-order derivative. To have a unique solution more information of the system's states is required. In most applied problems this information on the boundary is known. These problems are referred to as boundary value problems. Initial value problems, terminal value problems (TVPs) and Sturm-Liouville problems are among boundary value problems. For initial value problems, the states of the system in the beginning time of the model are known. These states may include derivatives or more complex operators of the modeled process in the initial time. For initial value problems, the states of the system in the beginning time/point of the model are known. These states may include derivatives or more complex operators of the modeled process in the initial time. For terminal value problems, the system's states are known at the end terminal. However, for Sturm-Liouville problems states of the system are known in both initial and end boundaries. These systems have completely different dynamics and behaviors [6].
Adding boundaries information is not the only way to reduce the degree of freedom. For example, the information in the intermediate point/time can be used [7] for such a reduction. Such problems are known as intermediate value problems.
Recently, TVPs for fractional systems of the order less than one (0<α<1) have received some extensive attention and interest [8,9,10,11,12,13,14]. These problems can be described by
where Dα is fractional derivative and qa is a given ν-dimensional vector. At the first glance, these systems seem to be well-posed with regular source functions [15,16]. Surprisingly, Cong and Tuan [17] revealed some counterexamples. The papers [14,18] pointed out that the well-posedness depends on the terminal value a. For larger a we may not obtain a unique solution.
Surprisingly, TVPs for fractional differential equations of higher-order derivatives α≥1 are not well-studied. An order α∈(1,2] problem on an infinite interval have been studied in [19]. These problems have the form
where q′(∞)∈R is a given number and C0Dαt stands for Caputo derivative.
This paper's important aim is to study the existence results and regularity of the higher-order terminal value problem for systems of FDEs. The other important aim of this paper is to introduce an analyzed high-order numerical method for solving such problems. Thus, the piecewise polynomial collocation method (PPCM) as a numerical solver is introduced and analyzed in detail.
The PPCMs have some superior advantages in solving differential or integral equations. For example, if we use polynomials collocation methods (known as spectral methods) we may encounter Runge-phenomena. The piecewise characteristic of the PPCMs avoids such probable divergence. In comparison with Runge-Kutta methods, the coefficient and parameters of the given methods are obtained constructively. It is remarkable to mention that for ordinary differential equations PPCMs, for some piecewise polynomial spaces are equivalent to Runge-Kutta methods. Conclusively, having more parameters (collocation parameters, degree of polynomial space, step size, meshing method) for controlling error, complexity and order of convergence makes the PPCMs competitively superior category of methods.
This paper is organized as follows: In Sections 2 and 3 we obtain the existence and uniqueness of a mild solution for higher-order FDEs. The regularity for the linear case is studied in Section 4. A numerical method is introduced in Section 5 and analysand in Section 6. Supportive numerical experiments are given in Section 7.
2.
Higher order terminal value problem
Consider the system (1.1) with terminal values
where n−1<α<n and n∈N. The ν-dimensional vector function q=[q1,…,qν] is an unknown vector and q(i)a∈Rν (i=0,…,n−1) are known vectors.
Further, we impose the Lipschitz condition on components of the function p:[0,a]×Rν→Rν. Thus, we suppose
where Li (i=0,…,ν) are constants and are not depend on t. Let z=q(n−1), (z′=q(n)). Taking repeatedly integrals and using terminal values, we obtain
By definition of Caputo derivative, we have
Acting Riemann-Liouville integral 0Jα−n+1t on both sides of Eq (1.1), we obtain
Finally, putting t=a and using terminal conditions, we obtain
Remark 1. The solution of the coupled system of (2.2) and (2.4) can be regarded as a mild solution of the system (2.1).
3.
Existence of a unique solution
Since we use vector functions, we need a norm combined with vector norm and function norm. This combinations brings us to face complexity in calculating induced norm. In this respect, the max norm defined by
seems to be more easier. We establish well-posedness of the inverse problems (1.1) and (2.1). We define the operators P and Q:(C[0,a])2ν→(C[0,a])ν by the right hand sides of systems (2.2) and (2.4), i. e.,
and
Let Pi and Qi:(C[0,a])2ν→C[0,a] be the ith component of the operator P and Q, respectively. Setting w=[q,z] and defining the operator T:R2×ν→R2×ν by
the Eqs (2.2) and (2.4) can be compactly written as
It is straightforward to show that for w1 and w2∈(C[0,a])2ν. We can write
and
where LM=maxi=1,…,νLi. Now, the Eqs (3.4) and (3.5) give
Conclusively, we can estate the following theorem.
Theorem 3.1. Let each component of the continuous function p:[0,a]×Rν→Rν (ν∈N) satisfy the Lipschitz condition with the Lipschitz constant Li and let
Then, the problem (1.1) with terminal conditions (2.1) has a unique mild solution on (C[0,a])2ν if
Proof. Suppose w∈(C[0,a])2ν. Since all integral operators involving in the definition of T transform a continuous functions into a continuous functions, therefore, T(w)∈(C[0,a])2ν. The operator T is contractor by (3.6) and the space (C[0,a])2ν is a Banach space by the norm ‖.‖∞. Thus, by contraction mapping theorem w=T(w) has a unique solution which completes the proof.
The next corollary is an important result of Theorem 3.1 which shows the regularity of the solution on more smooth spaces.
Corollary 1. Suppose the hypotheses of Theorem 3.1 are satisfied. Then, the mild solution q obtained from solving systems (2.2) and (2.4) satisfies the system (1.1) with terminal conditions (2.1). Moreover,
Proof. Taking derivative from Eq (2.2) and using Theorem 3.1 results q(n−1)=z∈(C[0,a])ν. Thus, q∈(Cn−1[0,a])ν. Similarly, putting t=a in Eq (2.2) and its derivatives, proves that q satisfies terminal conditions. Replacing z=q(n−1) in (2.2) we obtain
where aJα−n+1t is the Riemann-Liouville fractional integral. Noting that α−n+1>0, we can take a fractional derivative of order α−n+1 from (3.8) to obtain
which shows that the component q of the mild solutions satisfies (1.1).
4.
Regularity for linear FDEs
The regularity of the solution is important for analyzing numerical methods, especially finite difference methods and methods based on projection [20]. The regularity speaks about the order of differentiability of solutions. Usually, the regularity is investigated in the space Cm[a,b], (see Corollary 1). But, some good functions such as √t∈C[a,b] do not belong to Cm[a,b]. The space Cm(a,b]∩C[a,b] contains such functions. This space is not Banach or complete. Therefore, we need to introduce another Banach space such that
In this paper, we study regularity on the following weighted space Cm,α(0,a] introduced in [13].
Definition 4.1. [13] Let 0≤α<1 and m∈N. Then q∈Cm,α(0,a] if there exist functions qi∈C[0,a] for i=0,⋯,m+1 such that q=tαq0+q1 and
Theorem 4.2. The space Cm,α(0,a] with the norm
is a Banach space.
Proof. Let qn∈Cm,α(0,a], n∈N be a Cauchy sequence with norm ‖.‖α,m. Thus, by the definition of Cm,α(0,a] there exists functions qn,i in C[0,a] for i=0,⋯,m+1 such that qn(t)=tαqn,0(t)+qn,1(t) and
From the definition of the norm ‖.‖α,m, the sequences qn,i for a fixed i are Cauchy sequence in the Banach space C[0,a] and thus have a unique limit say qi. Let q defined by q(t)=tαq0(t)+q1(t) on [0,a]. Following items show that q∈Cm,α(0,a].
● limn→∞qn(t)=limn→∞tαqn,0(t)+qn,1(t)=tαq0(t)+q1(t)=q for all t∈[a,b];
● q∈C[0,a], since q0∈C[0,a] and q1∈C[0,a];
● limn→∞Diqn(t)=limn→∞tα−iqn,i+1(t)=tα−iqi+1(t). For each t∈(a,b] we can find ϵ>0 such that Dϵ=[t−ϵ,t]⊂(a,b]. It is trivial that Diqn(t) converges uniformly to tα−iqi+1(t) on Dϵ. According to Theorem 7.17 of [21] Diq(t)=limn→∞Diqn(t) on Dϵ. Thus Diq(t)=tα−iqi+1(t) for all t∈Dϵ and hence for all t∈(a,b].
Let ν∈N. For a ν dimensional vector functions, p=[f1,…,fν] in (Cm,α(0,a]⊂(Cm(0,a])r, we use the max norm defined by
The system (1.1) with
is a system of linear FDEs. Here A is a given ν×ν dimensional matrix function and b∈(C[0,a])ν is a given source function. To study linear systems we introduce the operators W1 and W2:(C[0,a])2ν→(C[0,a])ν by
and
Let us also define the vector valued functions G2 and G1:[0,a]→Rν by
and
Let T=[W1,W2]T and G=[G1,G2]T. Then, the Eqs (2.2) and (2.4) can be written as an inhomogeneous system
where w=[q,z]T.
Theorem 4.3. Let a∈R and α>0. Assume each component of A and b are in Cm,α(0,a]) and the hypotheses of Theorem 3.1 are satisfied for related p:[0,a]×Rν→Rν (ν∈N).Then, the problem (1.1) with terminal conditions (2.1) has a unique mild solution on (Cm,α(0,a])2ν.
Proof. Since T is a combination of weakly singular integral operators, it is a compact linear operator on (C[0,a])2ν and (Cm,α(0,a])2ν. Also, it is clear that
System (25) has a unique solution on (C[0,a])2ν by Theorem 3.1. Therefore, the homogeneous system w−T(w)=0 has the trivial null space in (C[0,a])2ν and thus (Cm,α(0,a])2ν. Conclusively, alternative Fredholm theorem asserts that the system (4.7) has a unique solution on (Cm,α(0,a])2ν.
Remark 2. Regularity of the nonlinear case needs further investigations.
Example 1. Consider the system (1.1) with a=0.5, α=2.5 and p defined by (4.2)
Obviously, L1=0.5, L2=0.5. and thus LM=0.5. Therefore,
and Condition (3.7) holds. By Theorem 3.1 the terminal value problem (1.1) with arbitrary terminal value (2.1) has a unique solution on C[0,a]. Since all components of A and b belong to Cm,α(0,a], (m∈N) the mild solution belongs to (Cm,α(0,a])4, by Theorem 4.3. Furthermore, we obtain the regularity of the solution of the system (1.1) on the closed interval [0,a] by Corollary 1 and we have
5.
Numerical method
Considering the memory, one of the best methods for solving the coupled systems (2.2)–(2.4) is to use collocation methods on piecewise polynomial spaces. To implement such methods we partition the solution interval [0,a] into sub-intervals σi=[ti,ti+1], i=0,…,N−1, with length hi:=ti+1−ti where 0=t0<…<tN=a are nodes of a chosen mesh (uniform or graded mesh) and N∈N. A graded mesh with exponent r≥1 is described by
Let 0<c1<…<cm≤1 (m∈N) be collocation parameters, tn,i=tn+cihn (n=0,…,N−1) be collocation points and let ˆqN(t) and ˆzN(t) be approximate solutions. The restriction of approximate functions to the σk is fully determined by Lagrange polynomials interpolation formula
and
where Qk,j=ˆqN(tn,j), Zk,j=ˆzN(tk,j), tk,j=tk+hcj and Lj (j=1,…,m) are Lagrange polynomials of degree m−1. The integrals in the operators P and Q of systems (3.1) and (3.2) can be discretized as:
where
For t∈[tk,tk+1] (k=0,…,N−1) we have
where
and
Similarly, we discretize integrals of the operator Q. We have
By interpolating p(tl+shl,q(tl+shl)) on ci we obtain
and thus
where
and finally
where
The discretized mappings PN and QN:(C[0,a])2ν→(PC[0,a])ν (PC stands for piecewise continuous space) related to systems (3.1) and (3.2) can be defined by
and
on σk. Setting TN=[PN,QN], unknowns vectors Ql,j and Zl,j can be obtained by solving the nonlinear system
for k=0,…,N−1 and o=1,…,m. Taking into account that
the dense solution can be evaluated in all points of the desired interval by Eqs (5.1) and (5.2). We note that the six parameters γn,i,l,j, ηn,l,j, ˉηn,k,j, ϱn,l,j, ζn,l,j and ˉζn,k,j fully determines the method for each m. The equations of these six parameters for m=1 and m=2 are provided in Tables 1 and 2, respectively.
6.
Numerical analysis
For simplifying our analysis we recall some notations. Let S−1m−1(Ih) be the space of piecewise polynomials of degree less than m, (m∈N) on the mesh partitioning Ih={ti:i=0,⋯,N}. The projection operator Πm−1,N:(C[0,a])2ν→(S−1m−1(Ih))2ν is uniquely determined by interpolation on collocation points such that
Since ˆwN(t):=[ˆqN(t),ˆzN(t)]∈(S−1m−1(Ih))2ν we have
Thus, Eq (5.10) can be written as
A useful theorem regarding Πm−1,N that simplifies existence and convergence analysis is the following theorem.
Theorem 6.1. Let m,N∈N. Then, Πm−1,N is a bounded linear operator. Let
be induced norm of Πm−1,N and
Then, Λ≤Λ1 and
Proof. In each sub-interval σn of [0,a] by Lagrange interpolation formula, we have
The rest of the proof is straightforward by taking the max norm.
6.1. Existence results for the proposed numerical method
One of the most fundamental questions is whether the system (6.1) has a unique solution? The answer is affirmative. For ˆwN(t) in (S−1m−1(Ih))2ν the operators PN of (5.8) and P of (3.1) are equivalent and we have
Thus, the inequality (3.4) holds for PN and
for all w1,w2∈(S−1m−1(Ih))2ν. However, QN is different from Q in this space and we need further computing. Actually, QN can be described by
where q is the first ν element of w. More precisely
Therefore,
Theorem 6.2. Let
Then, the numerical solution of system (1.1) obtained by (6.1) exists and is unique.
Proof. The operator Πm−1,NTN:(S−1m−1(Ih))2ν→(S−1m−1(Ih))2ν is a contractor by (6.7). Thus, by using contraction mapping theorem, Πm−1,NTN admits a unique fixed-point.
Remark 3. Obviously, if
holds, then Eq (6.8) holds too. For m=1 we have Λ1=1, and Eq (3) matches with Eq (3.7). However, our estimate may not be optimal for higher degrees of approximations. Also, we guess using the convergence properties of Πm−1,N when N→∞, we may obtain a convergence result without dependency on Λ1.
6.2. Detailed implementation: Solving nonlinear system
An important question is how to solve the nonlinear Eq (6.1). There are many methods available in literature that we can employ. The Newton iteration method is one of them [14]. The advantage is that it is fast. The disadvantage is that it needs computation of Jacoby as well as a good initial value for an iteration. It increases competition cost for higher dimension (computation of ν2 function for each iteration) and complexity of the implementation. However, strict restriction of λΛ and analytic discussion of the previous section is constructive and suggests the use of an iterative method. Since the operator Πm−1,NTN is a contractor operator, beginning with an arbitrary initial value grantees the convergence of the related iterative method. In our algorithm, we initialize the iteration by
Then, we estimate the solution by the recursive formula
The computation of iteration (6.10) continues with the smallest κ such that
where Tol is a desired user tolerance. We report such κ as a number of iteration to achieve the given tolerance.
6.3. Convergence analysis
Theorem 6.3. Suppose ˆwN=[ˆwN,ˆwN] (N∈N) be a discretized collocation solution of (2.2) and (2.4) described by (6.1). Let the conditions of Theorems 3.1 and 6.2 be fulfilled. Then, the coupled systems (2.2) and (2.4) have a unique solution w on (C[0,a])2ν and the ˆwN converges to the w. Furthermore, supposing w∈(Cm,α(0,a])2ν implies
Proof. Let w be the solution of the coupled systems (2.2) and (2.4). Define the residue operator by
and let e=ˆwN−w. Subtracting Eq (3.3) from (6.1), we obtain
Taking maximum norm, we obtain
Applying Eq (6.8), we obtain
and finally
This is a good news since we can obtain the asymptotic behavior of ‖R(w)‖∞ for the known w. Thanks to Theorem 8 of [20], one can easily prove that
which completes the proof.
7.
Examples
All numerical experiments is implemented in MATLAB software. The tolerance number for solving related nonlinear equations is Tol=10−14. This tolerance is close to the floating-point relative accuracy 2.2204e−16 for data type of double-precision (64 bit) in standard machine. The numerical examples show theoretically the obtained order of convergence is sharp and can not be improved more.
Example 2. Let 1<α≤2 and
on [0,0.15]. We note that p(t,.)∈(C[0,0.25])2 but p(t,.)∉(C1[0,0.25])2. Since LM=0.1, one can check that the conditions of Theorem 3.1 for a=0.25 hold and the system (2.1) with boundary condition
and
has a unique solution on [0,a]. The solution by Corollary 1 belongs to (C1[0,a])2. We apply a numerical method based on collocation parameter c=[0.5,1]. Thus, the dense solution on each sub-interval σk (k=0,⋯,N−1) of a given partition can be described by
An estimated norm of the given method can be determined by (6.2) and we have
Taking into account (6.8), we obtain
Consequently, by using Theorem 6.2 the nonlinear Eq (6.1) has a unique solution and the iteration process (6.10) converges to that solution. Theorem 6.3 implies the numerical solution converges to the exact solution by following order (O)
In Tables 3 and 4, we provided the maximum error
estimated order of proposed method (ON) and the number of iteration for required tolerance. As we see ON=0.5,1,1.5,2 for r=1,2,3,4 respectively, which is in agreement with theoretical result (see Eq (7.2)).
Remark 4. Theorem 1 provides a sharp result when we consider all the components of the state of the system 1 (w=[q,z]). However, it can be improved for the desired state of the solution q. According to Corollary 1, q has better regularity for n≥1. Thus, we guess that the order of convergence for q component with uniform mesh r=1, can be greater than or equal to n−1. Thus, for that component, a greater graded mesh exponent is not necessary. The error Eq(r,N)=‖ˆqN−q‖∞ and the estimated order of convergence of Example 2 is reported in Table 5. The order 2 is achieved with r=1.5.
To show the sharpness of theoretical results, we add numerical experiments for the case m=1 with c=0.5 (for other choice of c on [0,1] we have similar results). The dense solution on each sub-interval σk can be described by
In this case Λ1=1 and by (6.8)
Consequently, by using Theorem 6.2 the nonlinear Eq (6.1) related to m=1 has a unique solution and the iteration process (6.10) converges to the exact solution. By Theorem 6.3 we expect
The sharpness of this result can be seen in Table 6.
Remark 5. The low order method have the advantage of supporting larger class of terminal value problem with large value of a. The reason is that with m=1, we have Λ1=1 which is three times smaller than Λ1 for the case m=2.
8.
Conclusions
TVPs for higher-order (greater than one) fractional differential equations are rarely studied in the literature. In this paper, we tried to have a comprehensive study with a simple analysis to cover all general interests. Many questions in this topic deserve to study with more scrutiny. The regularity of nonlinear cases as well as optimal bound for obtaining well-posed problems are among them. We think the terminal value problem is important in applied since it monitors the past evolution of a dynamical system. Also, it can be regarded as a control problem for having the desired future by finding suitable initial value. We think this topic will catch more interest similar to the initial value problem.
Conflict of interest
The authors declare that there are no conflicts of interest.