Next Article in Journal
Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes
Next Article in Special Issue
First Solution of Fractional Bioconvection with Power Law Kernel for a Vertical Surface
Previous Article in Journal
Identification of Inertial Parameters for Position and Force Control of Surgical Assistance Robots
Previous Article in Special Issue
Nonlocal Sequential Boundary Value Problems for Hilfer Type Fractional Integro-Differential Equations and Inclusions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth

1
Department of Mathematics, Indian Institute of Technology Patna, Patna 801106, Bihar, India
2
Department of Mathematics, Texas A & M University-Kingsville, 700 University Blvd., MSC 172, Kingsville, TX 78363-8202, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(7), 774; https://doi.org/10.3390/math9070774
Submission received: 15 March 2021 / Revised: 24 March 2021 / Accepted: 26 March 2021 / Published: 2 April 2021
(This article belongs to the Special Issue New Trends on Boundary Value Problems)

Abstract

:
In this work, the existence and nonexistence of stationary radial solutions to the elliptic partial differential equation arising in the molecular beam epitaxy are studied. Since we are interested in radial solutions, we focus on the fourth-order singular ordinary differential equation. It is non-self adjoint, it does not have exact solutions, and it admits multiple solutions. Here, λ R measures the intensity of the flux and G is stationary flux. The solution depends on the size of the parameter λ . We use a monotone iterative technique and integral equations along with upper and lower solutions to prove that solutions exist. We establish the qualitative properties of the solutions and provide bounds for the values of the parameter λ , which help us to separate existence from nonexistence. These results complement some existing results in the literature. To verify the analytical results, we also propose a new computational iterative technique and use it to verify the bounds on λ and the dependence of solutions for these computed bounds on λ .

1. Introduction

Epitaxy means the growth of a single thin film on top of a crystalline substrate. It is crucial for semiconductor thin film technology, hard and soft coatings, protective coatings, optical coatings, etc. The epitaxial growth technique is used to produce the growth of semiconductor films and multilayer structures under high vacuum conditions [1]. The major advantages of epitaxial growth are reducing the growth time, better structural and superior electrical properties, eliminating waste caused during growth, wafering cost, cutting, polishing, etc. Several types of epitaxial growth techniques, such as hybrid vapor phase epitaxy [2], chemical beam epitaxy [3], and molecular beam epitaxy (MBE), have been used for the growth of compound semiconductors and other materials. In this work, we strictly focus on MBE, and we restrict our attention to the differential equation model, which was proposed by Escudero et al. [4,5,6,7]. The mathematical description of epitaxial growth is carried out by means of a function σ defined as
σ : Ω R 2 × R + R ,
which describes the height of the growing interface in the spatial point x Ω R 2 at time t R + . The authors in [4,5,6,7] show that the function σ obeys the fourth-order partial differential equation
t σ + Δ 2 σ = det ( D 2 σ ) + λ η ( x , t ) , x Ω R 2 ,
where η ( x , t ) models the incoming mass entering the system through epitaxial deposition, λ measures the intensity of this flux, and the determinant of Hessian matrix is
det ( D 2 σ ) = 2 σ x 1 2 · 2 σ x 2 2 2 σ x 1 x 2 2 .
The stationary counterpart of the partial differential Equation (1) subject to the homogeneous Dirichlet boundary condition (4) and homogeneous Navier boundary condition (5) is defined as (see [6])
Δ 2 σ = det ( D 2 σ ) + λ G ( x ) , x Ω R 2 ,
σ = 0 , σ n = 0 on Ω ,
σ = 0 , Δ σ = 0 on Ω ,
where η ( x , t ) G ( x ) is a stationary flux, and n is the unit out drawn normal to Ω .
Using the transformation r = | x | and σ ( x ) = ϕ ( | x | ) , as a result of symmetry, the above set of equations are transformed into the following set of equations:
1 r r 1 r ( r ϕ ) = 1 r ϕ ϕ + λ G ( r ) ,
ϕ ( 0 ) = 0 , ϕ ( 1 ) = 0 , ϕ ( 1 ) = 0 , lim r 0 r ϕ ( r ) = 0 ,
ϕ ( 0 ) = 0 , ϕ ( 1 ) = 0 , ϕ ( 1 ) + ϕ ( 1 ) = 0 , lim r 0 r ϕ ( r ) = 0 ,
where = d d r .
In this paper, we also impose the following boundary conditions which complements the work in [6]:
ϕ ( 0 ) = 0 , ϕ ( 1 ) = 0 , ϕ ( 1 ) = 0 , lim r 0 r ϕ ( r ) = 0 .
For simplicity, we take G ( r ) = 1 , which physically means that the new material is being deposited uniformly on the unit disc.
Now, using lim r 0 r ϕ ( r ) = 0 , w = r ϕ and integrating parts from Equation (6), we have
r 2 w r w = 1 2 w 2 + 1 2 λ r 4 .
Using the transformation t = r 2 2 and u ( t ) = w ( r ) , it is possible to reduce Equation (10) into the following equation:
u = u 2 8 t 2 + λ 2 , for t 0 , 1 2 .
Corresponding to (11), we define the following three boundary value problems:
Problem 1 : u = u 2 8 t 2 + λ 2 , for t 0 , 1 2 lim t 0 + t u ( t ) = 0 , u 1 2 = 0 ,
Problem 2 : u = u 2 8 t 2 + λ 2 , for t 0 , 1 2 lim t 0 + t u ( t ) = 0 , u 1 2 = 0 ,
Problem 3 : u = u 2 8 t 2 + λ 2 , for t 0 , 1 2 lim t 0 + t u ( t ) = 0 , u 1 2 = u 1 2 .
The BVPs (12), (13) and (14) can be equivalently described as the following integral equations (IE):
  • IE corresponding to Problem 1:
    u ( t ) = 1 2 t 0 t u ( s ) 2 4 s d s + t t 1 2 u ( s ) 2 4 s 2 1 2 s d s + λ 4 t 1 2 t ,
  • IE corresponding to Problem 2:
    u ( t ) = 0 t u ( s ) 2 8 s d s + t t 1 2 u ( s ) 2 8 s 2 d s + λ 4 t 1 t ,
  • IE corresponding to Problem 3:
    u ( t ) = t + 1 2 0 t u ( s ) 2 4 s d s + t t 1 2 u ( s ) 2 4 s 2 s + 1 2 d s + λ 4 t 3 2 t .
We assume that u C l o c 2 0 , 1 2 ; R , where C l o c 2 0 , 1 2 ; R is defined as
u : 0 , 1 2 R | u C 2 [ a , b ] , R for everycompact set [ a , b ] 0 , 1 2 .
In [6], Escudero et al. proved the existence and nonexistence of solutions of Problems 1 and 3 using upper and lower solution techniques. Corresponding to Problems 1 and 3, they have also provided the rigorous bounds of the values of the parameter λ , which helps us to separate existence from nonexistence. In [8], Verma et al. provide numerical illustrations via VIM to verify the results of Escudero et al. [6]. To verify their numerical results, they provided other iterative schemes based on homotopy [9] and the Adomian decomposition method [10].
Equation (13) has not been investigated theoretically in the existing literature to the best of our knowledge. Moreover, many investigations are still pending relating to BVPs (12), (13), and (14). Here, we focus on both theoretical and numerical work. We derive the sign of the solution and prove its existence in continuous space. We also compute the bounds of the parameter λ . The results of this paper complements existing theoretical results. We also provide an iterative scheme based on Green’s function to compute the bounds and solutions to demonstrate the existence and nonexistence, which is dependent on λ .
To prove the existence of the solutions, we use the monotone iterative technique [11,12,13,14,15,16,17]. Recently, many researchers applied this technique on the initial value problem (IVP) for the nonlinear noninstantaneous impulsive differential equation (NIDE) [18], p-Laplacian boundary value problems with the right-handed Riemann–Liouville fractional derivative [19], etc. to prove the existence of the solution. Here, we also present numerical results to verify the theoretical results. To develop the iterative scheme based on Green’s function, we consider Equations (12)–(14). Recently, many authors have used numerical approximate methods like the VIM [8], the Adomian decomposition method (ADM), the homotopy perturbation method (HPM) etc. to find approximate solutions for different models involving differential equations [20,21], integral equations [22,23,24], fractional differential equations [25,26], the Stefan problem [27,28,29,30], system of integral equations [31], etc. Thereafter, Waleed Al Hayani [32] and Singh et al. [33] applied ADM with Green’s function to compute the approximate solution. Recently, Noeiaghdam et al. [34] proposed a technique based on ADM for solving Volterra integral equation with discontinuous kernels using the CESTAC method. To find out more about this method, please see [35,36]. They focused on the BVPs which have a unique solution. The major advantage of our proposed technique is its ability to capture multiple solutions together with a desired accuracy.
The remainder of the paper is focused on both theoretical and numerical results. We prove some of the basic properties of the BVPs in Section 2. The monotone iterative technique is presented in Section 3 to prove the existence of a solution. A wide range of λ for Equation (6), corresponding to different types of boundary conditions, is shown in Section 4. In Section 5, we apply our proposed technique to the integral equations and show a wide range of numerical results. Finally, in Section 6, we draw our main conclusions.

2. Preliminary Work

Corresponding to λ 0 , we prove some basic qualitative properties of the solution u C l o c 2 ( 0 , 1 2 ] , R , which satisfies the following inequality:
u u 2 8 t 2 + λ 2 , for t 0 , 1 2 .
Here, we omit the proof of Lemmas 1–3, Corollary 1, and Lemma 4, which were done by Escudero et al. in [6].
Lemma 1.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 + t u ( t ) = 0 and Equation (18), then lim t 0 u ( t ) = 0 .
Lemma 2.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 u ( t ) = 0 , u 1 2 = 0 and Equation (18), then u ( t ) 0 for all t 0 , 1 2 .
Lemma 3.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 u ( t ) = 0 , u 1 2 = u 1 2 and Equation (18), then u ( t ) 0 for all t 0 , 1 2 .
Corollary 1.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 u ( t ) = 0 , u ( t ) 0 and Equation (18), then lim t 0 + t u ( t ) = 0 if and only if lim t 0 + u ( t ) t = 0 .
Lemma 4.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 + u ( t ) t = 0 . Then, for every μ [ 0 , 1 ) , we have
lim t 0 + t 1 μ t 1 2 u 2 s 2 d s = 0 .
Lemma 5.
Let u C l o c 2 ( 0 , 1 2 ] , R satisfy lim t 0 u ( t ) = 0 , u 1 2 = 0 and Equation (18), then u ( t ) 0 for all t 0 , 1 2 .
Proof. 
First, we show that u 1 2 0 . Assume u 1 2 > 0 . Since lim t 0 u ( t ) = 0 , there exist a t 0 0 , 1 2 such that u ( t 0 ) < u 1 2 . Now, from (18), u ( t ) is an increasing function on 0 , 1 2 . Again, by mean value theorem, we have
u 1 2 u ( t 0 ) 1 2 t 0 = u ( ξ ) , min 1 2 , t 0 ξ max 1 2 , t 0 .
Since u 1 2 = 0 , we have u 1 2 u ( t 0 ) u 1 2 1 2 t 0 = 0 . Hence, we get u 1 2 u ( t 0 ) , which is a contradiction. Therefore, we have u 1 2 0 . Furthermore, u ( t ) is a convex function along with u 1 2 = 0 . Moreover, u ( t ) is increasing, which implies u ( t ) 0 . Again, u ( t ) is a decreasing function on 0 , 1 2 . Therefore, lim t 0 u ( t ) = 0 and u 1 2 0 lead to u ( t ) 0 on 0 , 1 2 . □
Lemma 6.
Let u C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 3, then u ( t ) satisfies the following integral equation:
u ( t ) = t + 1 2 0 t u ( s ) 2 4 s d s + t t 1 2 u ( s ) 2 4 s 2 s + 1 2 d s + λ 4 t 3 2 t ,
and
lim t 0 + | u ( t ) | t < + .
Proof. 
The Green’s function of Problem 3 can be written as
G ( s , t ) = 2 s t + 1 2 , 0 s t , 2 t s + 1 2 , t s 1 2 .
Therefore, from Equation (23) and Problem 3, we can easily deduce the integral Equation (21). Now, using the result of Lemma 1, we have
lim t 0 + | u ( t ) | t = 0 1 2 u ( s ) 2 4 s 2 s + 1 2 d s + 3 λ 8 .
Now, put
f ( t ) = u ( t ) 2 t , g ( t ) = 1 t μ and h ( t ) = 1 t 1 μ for t 0 , 1 2 .
Therefore, we get f g L ( 0 , 1 2 ] provided μ ( 0 , 1 ) . Consequently, we have
0 1 2 u ( s ) 2 s 2 d s < .
Hence, from Equation (24), we get Equation (22). □
Lemma 7.
Let u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 2, then u ( t ) can be written in the following form:
u ( t ) = 0 t u ( s ) 2 8 s d s + t t 1 2 u ( s ) 2 8 s 2 d s + λ 4 t 1 t ,
and also satisfies
lim t 0 + | u ( t ) | t < .
Proof. 
Using the boundary condition and properties of Green’s function, we have
G ( s , t ) = s , 0 s t , t , t s 1 2 .
From Equation (29) and Problem 2, we can easily derive Equation (27). Now, using the result of Lemma 1, we have
lim t 0 + | u ( t ) | t = 0 1 2 u ( s ) 2 8 s 2 d s + λ 4 .
Therefore, from Equation (30) and using a similar analysis as that in Lemma 6, we can prove the result (28). □
Lemma 8.
Let u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 1, then u ( t ) can be written in the following form:
u ( t ) = 1 2 t 0 t u ( s ) 2 4 s d s + t t 1 2 u ( s ) 2 4 s 2 1 2 s d s + λ 4 t 1 2 t ,
and satisfies
lim t 0 + | u ( t ) | t < .
Proof. 
The Green’s function of Problem 1 is given by
G ( s , t ) = 2 s 1 2 t , 0 s t , 2 t 1 2 s , t s 1 2 .
Again, from Equation (33) and Problem 3, we derive integral Equation (31). Furthermore, using the result of Lemma 1, we have
lim t 0 + | u ( t ) | t = 0 1 2 u ( s ) 2 4 s 2 1 2 s d s + λ 8 .
Again, using a similar analysis as that in Lemma 7, we get the inequality (32). □

3. Existence of Solutions

In this section, we apply the monotone iterative technique coupled with lower and upper solutions to prove the existence of at least one solution for Problems 1–3. For this purpose, we need to prove some lemmas, which help us to prove the main results of this paper.

3.1. Construction of Green’s Function

To investigate the Problems 1–3, we consider the corresponding nonlinear singular boundary value problems, which are given by
Problem 1 ( a ) : u + k u = h ( t ) , for t 0 , 1 2 , lim t 0 + t u ( t ) = 0 , u 1 2 = b 1 ,
Problem 2 ( a ) : u + k u = h ( t ) , for t 0 , 1 2 , lim t 0 + t u ( t ) = 0 , u 1 2 = b 2 ,
Problem 3 ( a ) : u + k u = h ( t ) , for t 0 , 1 2 , lim t 0 + t u ( t ) = 0 , u 1 2 b 3 = u 1 2 ,
where h ( t ) = u 2 8 t 2 + λ 2 + k u , k R , b 1 , b 2 , b 3 0 , and λ R . Throughout the paper, we assume the following conditions:
  • H 0 = { k R : k < 0 } ;
  • H 1 = { k R : k < 0 | k | cosh | k | 2 sinh | k | 2 > 0 } ;
  • H 2 = { k R : 0 < k < 4 π 2 } ;
  • H 3 = { k R : 0 < k < π 2 } ;
  • H 4 = { k R : 0 < k < π 2 4 k cos k 2 sin k 2 > 0 } .
Lemma 9.
Let k satisfy H 0 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 1 ( a ) , then
u ( t ) = b 1 sinh | k | t sinh | k | 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = sinh | k | 2 | k | t sinh | k | s | k | sinh | k | 2 , 0 s t , sinh | k | 2 | k | s sinh | k | t | k | sinh | k | 2 , t s 1 2 ,
and G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
Using the boundary condition of Problem 1 ( a ) and the properties of Green’s function, we can easily prove Equation (39). Furthermore, we have G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 . □
Lemma 10.
Let k satisfy H 0 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 2 ( a ) , then
u ( t ) = b 2 sinh | k | t | k | cosh | k | 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = cosh | k | 2 | k | t sinh | k | s | k | cosh | k | 2 , 0 s t , cosh | k | 2 | k | s sinh | k | t | k | cosh | k | 2 , t s 1 2 ,
and G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
In a similar manner to that in Lemma 9, we can easily obtain Equation (41) and prove G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 . □
Lemma 11.
Let k satisfy H 1 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 3 ( a ) , then
u ( t ) = b 3 sinh | k | t | k | cosh | k | 2 sinh | k | 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = P ( t ) sinh | k | s | k | | k | cosh | k | 2 sinh | k | 2 , 0 s t , P ( s ) sinh | k | t | k | | k | cosh | k | 2 sinh | k | 2 , t s 1 2 ,
where P ( t ) = | k | cosh | k | 2 | k | t sinh | k | 2 | k | t . Moreover, G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
Again, using a similar analysis, we can easily derive the Green’s function. Now,
| k | cosh | k | 2 | k | t sinh | k | 2 | k | t = | k | cosh | k | 2 cosh | k | t | k | sinh | k | 2 sinh | k | t sinh | k | 2 cosh | k | t + cosh | k | 2 sinh | k | t = | k | cosh | k | 2 sinh | k | 2 cosh | k | t | k | sinh | k | 2 sinh | k | t + cosh | k | 2 sinh | k | t | k | cosh | k | 2 sinh | k | 2 cosh | k | t sinh | k | t 0 , since tanh | k | t 1 for all t 0 , 1 2 .
Hence, from (43), we have G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 . □
Lemma 12.
Let k satisfy H 2 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 1 ( a ) , then
u ( t ) = b 1 sin k t sin k 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = sin k 2 k t sin k s k sin k 2 , 0 s t , sin k 2 k s sin k t k sin k 2 , t s 1 2 ,
and G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
The proof is similar to that shown in Lemma 9. □
Lemma 13.
Let k satisfy H 3 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 2 ( a ) , then
u ( t ) = b 2 sin k t k cos k 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = cos k 2 k t sin k s k cos k 2 , 0 s t , cos k 2 k s sin k t k cos k 2 , t s 1 2 ,
and G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
The proof is similar to that shown in Lemma 10. □
Lemma 14.
Let k satisfy H 4 and u ( t ) C l o c 2 ( 0 , 1 2 ] , R be the solution of Problem 3 ( a ) , then
u ( t ) = b 3 sin k t k cos k 2 sin k 2 + 0 1 2 G ( s , t ) h ( s ) d s ,
where Green’s function G ( s , t ) is given by
G ( s , t ) = Q ( t ) sin k s k k cos k 2 sin k 2 , 0 s t , Q ( s ) sinh k t k k cos k 2 sin k 2 , t s 1 2 ,
where Q ( t ) = k cos k 2 k t sin k 2 k t . Moreover, G ( s , t ) 0 for all 0 s 1 2 and 0 t 1 2 .
Proof. 
The proof is similar to that shown in Lemma 11. □
Proposition 1.
Let k satisfy H 0 or H 2 (respectively, H 0 or H 3 and H 1 or H 4 ) and h ( t ) C l o c 2 ( 0 , 1 2 ] , R is such that h ( t ) 0 , then the solution of Problem 1 ( a ) (respectively, Problem 2 ( a ) and Problem 3 ( a ) ) is nonpositive.

3.2. Monotone Iterative Technique

Here, we define lower and upper solutions corresponding to Problems 1–3.
Definition 1
([37]). A function α C l o c 2 ( 0 , 1 2 ] , R is the upper solution of Problem 1 (respectively, Problem 2 and Problem 3) if
α α 2 8 t 2 + λ 2 , for t 0 , 1 2 ,
with lim t 0 α ( t ) t = 0 and α 1 2 0 (respectively, α 1 2 0 and α 1 2 α 1 2 ).
Definition 2
([37]). A function β C l o c 2 ( 0 , 1 2 ] , R is the lower solution of Problem 1 (respectively, Problem 2 and Problem 3) if
β β 2 8 t 2 + λ 2 , for t 0 , 1 2 ,
with lim t 0 β ( t ) t = 0 and β 1 2 0 (respectively, β 1 2 0 and β 1 2 β 1 2 ).
Now, we construct two sequences { α n } and { β n } corresponding to Problem 1 ( a ) (respectively, Problem 2 ( a ) and Problem 3 ( a ) ), which are defined by
α 0 = α ,
α n + 1 + k α n + 1 = α n 2 8 t 2 + λ 2 + k α n , for t 0 , 1 2 ,
lim t 0 α n + 1 ( t ) t = 0 and α n + 1 1 2 = 0 ,
respectively , α n + 1 1 2 = 0 and α n + 1 1 2 = α n + 1 1 2 and
β 0 = β ,
β n + 1 + k β n + 1 = β n 2 8 t 2 + λ 2 + k β n , for t 0 , 1 2 ,
lim t 0 β n + 1 ( t ) t = 0 and β n + 1 1 2 = 0 ,
respectively , β n + 1 1 2 = 0 and β n + 1 1 2 = β n + 1 1 2 . We assume the following properties:
  • P 1 : α 0 and β 0 satisfies
    lim t 0 | α 0 ( t ) | t < , lim t 0 α 0 ( t ) = 0 , α 0 ( t ) 0 ,
    and
    lim t 0 | β 0 ( t ) | t < , lim t 0 β 0 ( t ) = 0 ;
  • P 2 : h ( t , u ) is continuous on D 0 where
    D 0 = ( t , u ) 0 , 1 2 × R : β 0 = β u α 0 .
Now, we state our main existence theorems.
Theorem 1.
Assume H 0 (respectively, H 0 and H 1 ) is true, there exist α 0 , and β 0 C l o c 2 ( 0 , 1 2 ] , R are upper and lower solutions of Problem 1 (respectively, Problem 2 and Problem 3), which satisfy the properties P 1 and P 2 such that β 0 α 0 = 0 , then the Problem 1 (respectively, Problem 2 and Problem 3) has at least one solution in the region D 0 and the sequences { α n } , { β n } defined by (52)–(55) converge to solutions u, v uniformly and monotonically, respectively, such that
β u v α = 0 , t 0 , 1 2 .
Proof. 
We divide the proof into three parts. In the first part, we prove that
β n is a lower solution of problem 1 , β n β n + 1 , and β n + 1 α 0 n N { 0 } .
We apply mathematical induction on n. For n = 0 , from (54) and (55), we have
β 1 + k β 1 = β 0 2 8 t 2 + λ 2 + k β 0 , for t 0 , 1 2 ,
lim t 0 β 1 ( t ) t = 0 and β 1 1 2 = 0 .
Now, from Equation (51), we have
β 0 β 1 + k β 0 β 1 = β 0 2 8 t 2 λ 2 + β 0 0 , for t 0 , 1 2 ,
lim t 0 β 0 β 1 ( t ) t = 0 and β 0 β 1 1 2 0 .
Therefore, by Proposition 1, we have β 0 β 1 . Again from (50) and (60), we have
β 1 α 0 + k β 1 α 0 = β 0 2 8 t 2 + λ 2 + k β 0 α 0 α 0 , for t 0 , 1 2
β 0 + α 0 8 t + k t β 0 α 0 t .
Since β 0 α 0 , we have
β 1 α 0 + k β 1 α 0 0 , t 0 , 1 2 ,
lim t 0 ( β 1 α 0 ) ( t ) t = 0 and β 1 α 0 1 2 0 .
Hence, by Proposition 1, we have β 1 α 0 . Therefore, our assumptions are true for n = 0 . Let our assumptions be true up to n = m . Then, we find that
β n is a lower solution of problem 1 , β n β n + 1 , and β n + 1 α 0
for n = 1 , 2 , , m . Now, we want to show that our assumptions are true for n + 1 . Therefore, from Equation (54), we have
β n + 1 β n + 1 2 8 t 2 λ 2 = β n 2 β n + 1 2 8 t 2 + k β n β n + 1 , for t 0 , 1 2
β n + β n + 1 8 t + k t β n β n + 1 t .
Again, by using conditions (68), we have
β n + 1 β n + 1 2 8 t 2 + λ 2 , for t 0 , 1 2 ,
lim t 0 β n + 1 ( t ) t = 0 and β n + 1 1 2 0 .
Hence, β n + 1 is a lower solution of Problem 1. Now, from Equation (54) and (71), we have
β n + 1 β n + 2 + k β n + 1 β n + 2 = β n + 1 2 8 t 2 λ 2 + β n + 1 0 ,
lim t 0 ( β n + 1 β n + 2 ) ( t ) t = 0 and β n + 1 β n + 2 1 2 0 .
Therefore, by Proposition 1, we have β n + 1 β n + 2 . Again, from (50) and (54), we have
β n + 2 α 0 + k β n + 2 α 0 = β n + 1 2 8 t 2 + λ 2 + k β n + 1 α 0 α 0
β n + 1 + α 0 8 t + k t β n + 1 α 0 t .
Using a similar analysis, we have β n + 2 α 0 . Hence, by mathematical induction, we find that
β n is a lower solution of Problem 1 , β n β n + 1 and β n + 1 α 0 n N .
In the second part of the proof, we have to show that
α n is an upper solution of Problem 1 and α n + 1 α n n N .
Now, from (52) and (53), we have
α 1 + k α 1 = α 0 2 8 t 2 + λ 2 + k α 0 , for t 0 , 1 2 ,
lim t 0 α 1 ( t ) t = 0 and α 1 1 2 = 0 .
Therefore, by using (50), we have
α 1 α 0 + k α 1 α 0 = α 0 2 8 t 2 + λ 2 α 0 0 , for t 0 , 1 2 .
Again,
lim t 0 ( α 1 α 0 ) ( t ) t = 0 and ( α 1 α 0 ) 1 2 0 .
Hence, by Proposition 1, we have α 1 α 0 . Therefore, our assumptions are true for n = 0 . Let our assumptions be true up to n = m . Then, we find that
α n is an upper solution of Problem 1 and α n + 1 α n for n = 1 , 2 , , m .
Now, for n + 1 , we have
α n + 1 α n + 1 2 8 t 2 λ 2 = α n 2 α n + 1 2 8 t 2 + k α n α n + 1 , for t 0 , 1 2 ,
1 8 α n + α n + 1 t + k t α n α n + 1 t 0 .
Therefore,
α n + 1 α n + 1 2 8 t 2 + λ 2 , for t 0 , 1 2 ,
and
lim t 0 α n + 1 ( t ) t = 0 , α n + 1 1 2 0 .
Hence, α n + 1 is an upper solution of Problem 1. Therefore, by using (86), (52), and (53), we have
α n + 2 α n + 1 + k α n + 2 α n + 1 = α n + 1 2 8 t 2 + λ 2 α n + 1 0 ,
and
lim t 0 ( α n + 2 α n + 1 ) ( t ) t = 0 , ( α n + 2 α n + 1 ) 1 2 0 .
Therefore, by Proposition 1, α n + 2 α n + 1 . Hence, by mathematical induction, we conclude that
α n is an upper solution of Problem 1 and α n + 1 α n n N .
In the last part of the proof, we want to show β n α n for all n N . Again, from (71) and (86), we have
β n + 1 α n + 1 + k β n + 1 α n + 1 = β n 2 8 t 2 + λ 2 + k β n α n α n
β n + α n 8 t + k t β n α n t .
Since β n α n 0 , we have
β n + 1 α n + 1 + k β n + 1 α n + 1 0 ,
and
lim t 0 β n + 1 α n + 1 ( t ) t = 0 , β n + 1 α n + 1 1 2 0 .
Hence, by Proposition 1, β n + 1 α n + 1 . Finally, we have
β = β 0 β 1 β n α n α 1 α 0 = 0 .
Let t n 0 , 1 2 for n N such that
t n + 1 < t n for n N , lim n t n = 0 .
Therefore, for every n N , there exists a solution α n and β n to Equations (52) and (53), while (54) and (55) satisfy the inequality (95) on the interval [ t n , 1 2 ] . Since { α n } and { β n } are monotone and bounded, they converge to function u ( t ) and v ( t ) , respectively. Therefore, by Dini’s theorem, there exists u ( t ) and v ( t ) such that
lim n α n = u and lim n β n = v uniformly on every compact interval t n , 1 2
of 0 , 1 2 . Hence, from (52)–(55) and (38), there exists solutions v ( t ) C l o c 2 ( 0 , 1 2 ] , R and u ( t ) C l o c 2 ( 0 , 1 2 ] , R to Problem 1, satisfying
β u v α 0 = 0 , t 0 , 1 2 .
Hence, the proof is complete. □
Now, we assume the following conditions:
  • H 5 = k R : 0 < k < k , where k = min 4 π 2 , max t ( 0 , 1 2 ] α 0 2 t ,
  • H 6 = k R : 0 < k < k , where k = min π 2 , max t ( 0 , 1 2 ] α 0 2 t ,
  • H 7 = k R : 0 < k < k ,
    where k = min π 2 4 , max t ( 0 , 1 2 ] α 0 2 t k cos k 2 sin k 2 > 0
    .
Theorem 2.
Let α 0 , β 0 C l o c 2 ( 0 , 1 2 ] , R be the upper and lower solutions of Problem 1 (respectively, Problem 2 and Problem 3), which satisfy the properties P 1 and P 2 such that β 0 α 0 . Assume H 5 (respectively, H 6 and H 7 ) is true and λ R . Then, Problem 1 (respectively, Problem 2 and Problem 3) has at least one solution in the region D 0 and the sequences { α n } , { β n } defined by (52)–(55) converge to solutions u, v uniformly and monotonically, respectively, such that
β u v α , t 0 , 1 2 .
Proof. 
The proof is same as that shown in Theorem 1. □

4. Estimations of λ

The objective of this section is to derive some qualitative bounds of the parameter λ , from which we can conclude about the nonexistence of solutions. Equation (11) can be written in the following form:
( t u u ) = u 2 8 t + λ t 2 , t 0 , 1 2 .
Put v ( t ) = u ( t ) t and integrating from 0 to t, Equation (100) becomes
v ( t ) = 1 8 t 2 0 t v 2 ( s ) s d s λ 4 , t 0 , 1 2 .
Therefore, we have
v ( t ) 0 , t 0 , 1 2 .
In view of the transformation, the boundary condition at r = 1 becomes
BC of Problem 1 : v 1 2 = 0 ,
BC of Problem 2 : v 1 2 = 1 2 v 1 2 ,
BC of Problem 3 : v 1 2 = v 1 2 .
Escudero et al. in [6] prove the following two lemmas:
Lemma 15.
The set of numbers λ 0 , for which there exists a solution u ( t ) C l o c 2 ( 0 , 1 2 ] , R of Equation (11) satisfying lim t 0 u ( t ) t = 0 and u ( t ) 0 , is nonempty and bounded from above.
Lemma 16.
If Problem 1, Problem 2, and Problem 3 are solvable for some λ 0 0 , then these are solvable for every 0 λ λ 0 .
We present the following results which complement the results proved by Escudero et al. [6].
Lemma 17.
Let there exist a function u C l o c 2 ( 0 , 1 2 ] , R satisfying Equations (101), (102), and (104), then
λ 384 11 34.91 .
Proof. 
Now from Equation (101), we have
v ( t ) 0 , t 0 , 1 2 .
Again, from Equation (101), we get
v ( t ) = 1 4 t 3 0 t v 2 ( s ) s d s v 2 ( t ) 8 t , t 0 , 1 2 .
Therefore, by using (107) and (102), from (108), we have
v ( t ) 0 , t 0 , 1 2 .
Therefore, v ( t ) is increasing in 0 , 1 2 . Now,
v ( t ) v 1 2 = 1 2 0 1 2 v 2 ( s ) s d s λ 4 , t 0 , 1 2 .
Therefore, we have
v ( t ) c , t 0 , 1 2 ,
where
c = 1 2 0 1 2 v 2 ( s ) s d s + λ 4 .
Now, integrating Equation (111) from 0 to t and by using Equation (104), we have
v ( t ) c ( 1 t ) , t 0 , 1 2 .
Therefore, from Equations (112) and (113), we get
11 384 c 2 c + λ 4 0 ,
which implies Equation (106). □
Lemma 18.
Let
0 λ 2 C a n d C 128 9 ,
then there exists a solution β C l o c 2 ( 0 , 1 2 ] , R that satisfies Equation (51), the assumption P 1 , β 1 2 = 0 , and β ( t ) 0 .
Proof. 
We put
β ( t ) = C t ( A 2 t ) , t 0 , 1 2 and C 0 .
Obviously, β ( t ) satisfies assumption P 1 . Now, β 1 2 = 0 implies A = 3 2 . Therefore, β ( t ) 0 is also fulfilled. Now, we have
β ( t ) β 2 ( t ) 8 t 2 λ 2
= 3 C 2 2 t C 2 t 2 3 2 2 t 2 8 t 2 λ 2
= C 2 t 3 2 λ 2 t 2 C C 2 3 2 2 t 2 8
C 2 t 3 2 2 t C 2 3 2 2 t 2 8 , since λ 2 C
= C 2 8 2 t 3 2 2 t 2 t 3 4 2 + 9 C + 128 16 C
0 , t 0 , 1 2 , since C 128 9 .
Hence, the inequality (51) is satisfied. □
Lemma 19.
Let
0 λ 3 C a n d C 48 ,
then there exists a solution β C l o c 2 ( 0 , 1 2 ] , R that satisfies Equation (51), the assumption P 1 , β 1 2 = 0 , and β ( t ) 0 .
Proof. 
We put
β ( t ) = C t ( A 2 t ) , t 0 , 1 2 and C 0 .
Again, β ( t ) satisfies the assumption P 1 . Now, β 1 2 = 0 implies A = 1 . Hence, β ( t ) 0 is also fulfilled. Now, we have
β ( t ) β 2 ( t ) 8 t 2 λ 2
= 3 C 2 2 t C 2 t 2 1 2 t 2 8 t 2 λ 2
= 3 C 2 2 t 1 λ 2 t 3 C C 2 1 2 t 2 8
3 C 2 2 t 1 2 t C 2 1 2 t 2 8 , since λ 3 C
= C 2 8 2 t 1 2 t 2 t 1 2 2 + C + 48 4 C
0 , t 0 , 1 2 , since C 48 .
This completes the proof. □
Lemma 20.
Let
0 λ 3 C 2 a n d C 6 ,
then there exists a solution β C l o c 2 ( 0 , 1 2 ] , R that satisfies Equation (51), the assumption P 1 , β 1 2 = β 1 2 , and β ( t ) 0 .
Proof. 
We put
β ( t ) = C t ( A 2 t ) , t 0 , 1 2 and C 0 .
Now, β ( t ) also satisfies assumption P 1 . Similarly, β 1 2 = β 1 2 implies A = 2 . Therefore, β ( t ) 0 is also fulfilled. Then, we have
β ( t ) β 2 ( t ) 8 t 2 λ 2
= 3 C 2 2 t C 2 t 2 2 2 t 2 8 t 2 λ 2
= 3 C 4 2 t 2 2 λ 2 t 3 C C 2 2 2 t 2 8
3 C 4 2 t 2 2 t C 2 2 2 t 2 8 , since λ 3 C 2
= C 2 8 2 t 2 2 t 2 t 1 2 + C + 6 C
0 , t 0 , 1 2 , since C 6 .
Hence, the proof is complete. □
Theorem 3.
Let λ 0 R + . If 0 λ < λ 0 , then Equation (10) corresponding to different types of boundary conditions are solvable. Moreover, there is no solution to these problems if λ > λ 0 . Furthermore, every solution w ( r ) of a governing equation corresponding to these three types of boundary condition satisfy
w ( r ) 0 , r ( 0 , 1 ] a n d lim r 0 + w ( r ) = 0 .
Proof. 
The proof of this can be deduced from Lemma 15, Lemma 16, Lemma 1, Lemma 2, Lemma 3, and Lemma 5. □
Proposition 2.
Corresponding to Equations (6) and (7), the value of λ 0 admits the estimates
144 λ 0 307 .
Proof. 
From Lemma 7.7 in [6] and Lemma 19, we get Equation (140). □
Proposition 3.
Corresponding to Equations (6) and (9), the value of λ 0 admits the estimates
256 9 λ 0 384 11 .
Proof. 
From Lemmas 17 and 18, we have Equation (141). □
Proposition 4.
Corresponding to Equations (6) and (8), the value of λ 0 admits the estimates
9 λ 0 11.63 .
Proof. 
By using Lemma 7.6 in [6] and Lemma 20, we have Equation (142). □

5. Numerical Results and Discussion

Here, we present the numerical data to validate our derived theoretical results. In Section 5.1, we derive the numerical estimation of the bounds computed by ADM. In Section 5.3, we numerically show the existence of at least one solution.

5.1. ADM

To find the approximate solutions, we develop the iterative numerical schemes with the help of the Fredholm integral Equations (15), (16), and (17), respectively. Now, we decompose the solution u ( t ) of the form u ( t ) = i = 0 u i ( t ) , and approximate the nonlinear term in terms of Adomian’s polynomials [38], which is given by
N u ( t ) = 1 2 u 2 ( t ) = i = 0 A i ( u 0 , u 1 , , u i ) ,
where
A i = 1 i ! d i d β i N j = 0 i β j u j β = 0 , i = 0 , 1 , 2 , .
Therefore, from integral Equation (15), we define
Scheme of Problem 1 = u 0 ( t ) = c t λ 4 t 1 2 t , u n + 1 ( t ) = 0 t s 2 t 2 A n 2 s 2 d s , , and c = 0 1 2 i = 0 n A i 2 s 2 1 2 s d s .
We compute the arbitary constant c using the Mathematica 10.9 program. For better understanding, we present the algorithm of our proposed technique corresponding to Equation (15) below.
Residue Error:
Here, we define the residue error [39] corresponding to Equation (15) for error analysis, which is given by
R ( t ) = u ( t ) + 1 2 t 0 t u ( s ) 2 4 s d s + t t 1 2 u ( s ) 2 4 s 2 1 2 s d s + λ 4 t 1 2 t ,
where λ is the parameter. Therefore, the maximum absolute residue error can be defined as
L = max i = 0 , , 10 | R ( t 0 + i × 0.1 ) | , where t 0 = 0 .

5.2. Algorithm

Step 1 . Convert Fredholm integral Equation (15) into the Voltera integral equation.
Step 2 . Identify the constant term, and approximate the nonlinear term by Equation (143).
Step 3 . Consider u 0 ( t ) as in (145), and obtain u i ( t ) for i = 1 , 2 , , n + 1 .
Step 4 . Approximate the term u 2 4 s 2 by i = 0 n A i 2 s 2 in the equation c = 0 1 2 u 2 4 s 2 1 2 s d s .
Step 5 . Compute the values of the constant and the approximate solutions u ( t ) = i = 0 n + 1 u i .
Step 6 . Determine the residue error R ( t ) and set the stopping criteria L < ϵ , where ϵ is the tolerance.
Again, we apply the algorithm Section 5.2 to Equations (16) and (17), and we define the following iterative schemes:
Scheme of Problem 2 = u 0 ( t ) = c t λ 4 t 1 t , u n + 1 ( t ) = 0 t s t A n 4 s 2 d s , , and c = 0 1 2 i = 0 n A i 4 s 2 d s ,
and Scheme of Problem 3 = u 0 ( t ) = c t λ 4 t 3 2 t , u n + 1 ( t ) = 0 t s 2 t 2 A n 2 s 2 d s , , and c = 0 1 2 i = 0 n 1 2 + s A i 2 s 2 d s .
Approximate solutions for Equations (16) and (17) can be written as u ( t ) = i = 0 n + 1 u i ( t ) , provided the series is convergent for n . Recently, the convergence of ADM was established by Verma et al. in [9]. Now, by using the transformation t = r 2 2 , u ( t ) = w ( r ) , w ( r ) = r ϕ ( r ) , and ϕ ( 1 ) = 0 , we get the solutions of Equation (6). We arrive at two cases:
Case (a): λ 0 .
For λ = 0 , we get one trivial and one nontrivial solution. For 0 < λ λ critical , we always find two nontrivial solutions. We may refer to them as upper and lower solutions, respectively. Corresponding to Equations (9), (8), and (7), we find the critical values of λ , i.e., λ critical , are 31.94 , 11.34 , and 168.76 , respectively. For λ > λ critical , we do not find any numerical solutions, as the value of c become imaginary. In Section 5.2.1, we tabulate residual errors of the approximate solutions corresponding to some λ .
Case (b): λ < 0 .
In this case, we always have two nontrivial numerical solutions corresponding to three types of boundary conditions. One solution is negative (namely, the negative solution) and the other solution is positive (namely, the positive solution). We do not find any negative critical λ . Please refer to Section 5.2.1 as regards residue errors.

5.2.1. Tables

Here, we have placed some numerical data of approximate solutions of ϕ ( r ) corresponding to different types of boundary conditions below. If we are increasing the value of λ , we see that the residue error of the lower solution is increasing and the residue error of the upper solution is decreasing (see: Table 1). Similarly, if we are decreasing the value of negative λ , we see that the residue error of both positive and negative solutions are decreasing (see: Table 2). The same can be seen in Table 3, Table 4, Table 5 and Table 6.

5.3. Monotone Iterative Method

Here, we compute the monotone iterations using Equations (52)–(55) corresponding to three types of boundary condition.
Corresponding to problem (12): By using Lemma 19, we chose the lower and upper iterations
β 0 ( t ) = 3 t 1 2 t α 0 ( t ) = 0 , t [ 0 , 1 ] .
Therefore, it is easy to show that α 0 and β 0 both satisfy the inequalities (50), (51), (56), and (57), such that β 0 α 0 . We consider that k = 1 . Hence, by using Theorem 1, we have a monotonically and uniformly convergent sequence { β n } and { α n } , which are converging to the solution v and u of the problem (12). We denote u A D M as the approximation of the solution of (12) computed by ADM. By using the transformation t = r 2 2 , w ( r ) = u ( t ) , and w ( r ) = r ϕ ( r ) , we have the fourth-order iterations corresponding to the monotone iteration α i and β i . ϕ α i and ϕ β i are the fourth-order solutions corresponding to the monotone iteration α i and β i for i = 0 , 1 , , respectively.
In Figure 1a, we plotted β 0 , β 1 , u A D M , α 1 , and α 0 corresponding to problem (12) for λ = 2 . We have seen that the lower sequence { β n } and upper sequence { α n } always satisfies the inequality β n α n . In Figure 1b, we have placed the monotone iterations of fourth-order SBVP corresponding to problem (12) for λ = 2 . Here, we also observed the existence of at least one solution for a fourth-order SBVP corresponding to the Dirichlet boundary condition.
Corresponding to problem (13): From Lemma 18, we chose the initial monotone iterations as follows:
β 0 ( t ) = t 3 2 2 t α 0 ( t ) = 0 , t [ 0 , 1 ] .
The same remarks follow as discussed above.
In Figure 2a,b, we present the numerical results for k = 1 and λ = 2 . We noticed that the approximate solution u A D M and ϕ A D M always lies between the lower sequence { β n } and upper sequence { α n } .
Corresponding to problem (14): Here, we consider the initial monotone iterations
β 0 ( t ) = 2 3 t 2 2 t α 0 ( t ) = 0 , t [ 0 , 1 ] .
We also included the same remarks as those stated above.
Monotone lower and upper iterations corresponding to the second-order and the fourth-order differential equation are plotted in Figure 3.

6. Conclusions

In this work, we derived some qualitative properties of the singular boundary value problems that arise in the theory of epitaxial growth. Moreover, we proved the existence of a solution and discovered a range of parameter k, for which the nonlinear problem has multiple solutions in the region D 0 . We established the bounds of the parameter λ , from which we confirmed the nonexistence of solutions. Furthermore, the boundary value problems have multiple solutions, therefore it is challenging for researchers to obtain a suitable scheme to capture both solutions with the desired accuracy. However, we successfully developed iterative schemes and captured both solutions with a high accuracy. From Table 1, Table 2, Table 3 and Table 4, we can see that the approximate solutions computed by our proposed method converge to the exact solutions very quickly. Corresponding to the boundary conditions (7), we notice that the positive approximate solution converges to the exact positive solution very slowly (See Table 6). We verified that our numerical results matched well with our theoretical results as well as the existing numerical results [9]. We conclude that our proposed technique is relatively powerful and efficient. Furthermore, this technique is an effective tool to solve BVPs, which have multiple solutions.

Author Contributions

Conceptualization, A.K.V. and B.P.; validation, R.P.A. and B.P.; writing—original draft preparation, A.K.V. and B.P.; writing—review and editing, A.K.V., B.P. and R.P.A.; visualization, A.K.V., B.P. and R.P.A.; supervision, R.P.A.; project administration, A.K.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by grant provided by DST project, file name: SB/S4/MS/805/12 and INSPIRE Program Division, Department of Science & Technology, file no: IF160984 New Delhi, India-110016.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barabasi, A.L.; Stanley, H.E. Fractal Concepts in Surface Growth; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  2. Lourdudoss, S.; Kjebon, O. Hybrid vapor phase epitaxy revisited. IEEE J. Sel. Top. Quantum Electron. 1997, 3, 749–767. [Google Scholar] [CrossRef]
  3. Foord, J.S.; Davies, G.J.; Tsang, W.T. Chemical Beam Epitaxy and Related Techniques; John Wiley and Sons Ltd.: Chichester, UK, 1997. [Google Scholar]
  4. Escudero, C. Geometric principles of surface growth. Phys. Rev. Lett. 2008, 101, 1–4. [Google Scholar] [CrossRef] [Green Version]
  5. Escudero, C.; Hakl, R.; Peral, I.; Torres, P.J. On radial stationary solutions to a model of non equilibrium growth. Eur. J. Appl. Math. 2013, 24, 437–453. [Google Scholar] [CrossRef] [Green Version]
  6. Escudero, C.; Hakl, R.; Peral, I.; Torres, P.J. Existence and nonexistence result for a singular boundary value problem arising in the theory of epitaxial growth. Math. Methods Appl. Sci. 2014, 37, 793–807. [Google Scholar] [CrossRef] [Green Version]
  7. Escudero, C.; Korutcheva, E. Origins of scaling relations in non equilibrium growth. J. Phys. Math. Theor. 2012, 45, 1–14. [Google Scholar]
  8. Verma, A.K.; Pandit, B.; Escudero, C. Numerical solutions for a class of singular boundary value problems arising in the theory of epitaxial growth. Eng. Comput. 2019, 37, 2539–2560. [Google Scholar] [CrossRef]
  9. Verma, A.K.; Pandit, B.; Agarwal, R.P. On approximate stationary radial solutions for a class of boundary value problems arising in epitaxial growth theory. J. Appl. Comput. Mech. 2019, 4, 713–734. [Google Scholar]
  10. Verma, A.K.; Pandit, B.; Agarwal, R.P. On multiple solutions for a fourth order nonlinear singular boundary value problems arising in epitaxial growth theory. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  11. Cabada, A.; Cid, J.A.; Sanchez, L. Positivity and lower and upper solutions for fourth order boundary value problems. Nonlinear Anal. 2007, 67, 1599–1612. [Google Scholar] [CrossRef]
  12. Ladde, G.S.; Lakshmikantham, V.; Vatsala, A.S. Monotone Iterative Techniques for Nonlinear Differential Equations; Pitman Advance Publishing Program: Boston, MA, USA, 1985. [Google Scholar]
  13. Pandey, R.K.; Verma, A.K. Existence-uniqueness results for a class of singular boundary value problems-ii. J. Math. Anal. Appl. 2008, 338, 1387–1396. [Google Scholar] [CrossRef] [Green Version]
  14. Pandey, R.K.; Verma, A.K. Existence-uniqueness results for a class of singular boundary value problems arising in physiology. Nonlinear Anal. 2008, 9, 40–52. [Google Scholar] [CrossRef]
  15. Pandey, R.K.; Verma, A.K. On solvability of derivative dependent doubly singular boundary value problems. J. Appl. Math. Comput. 2010, 33, 489–511. [Google Scholar] [CrossRef]
  16. Pandey, R.K.; Verma, A.K. Monotone method for singular bvp in the presence of upper and lower solutions. Appl. Math. Comput. 2010, 215, 3860–3867. [Google Scholar] [CrossRef]
  17. Wang, G.; Agarwal, R.P.; Cabada, A. Existence results and the monotone iterative technique for systems of nonlinear fractional differential equations. Appl. Math. Lett. 2012, 25, 1019–1024. [Google Scholar] [CrossRef] [Green Version]
  18. Agarwal, R.; O’Regan, D.; Hristova, S. Monotone iterative technique for the initial value problem for differential equations with non-instantaneous impulses. Appl. Math. Comput. 2017, 298, 45–56. [Google Scholar] [CrossRef]
  19. Xue, T.; Liu, W.; Shen, T. Extremal solutions for p-laplacian boundary value problems with the right-handed riemann-liouville fractional derivative. Math. Methods Appl. Sci. 2019, 42, 4394–4407. [Google Scholar] [CrossRef]
  20. Mittal, R.C.; Nigam, R. Solution of a class of singular boundary value problems. Numer. Algorithm 2008, 47, 169–179. [Google Scholar] [CrossRef]
  21. Odibat, Z.; Momani, S. A reliable treatment of homotopy perturbation method for klein-gordon equations. Phys. Lett. A 2007, 365, 351–357. [Google Scholar] [CrossRef]
  22. Abbasbandy, S. Numerical solutions of the integral equations: Homotopy perturbation method and adomian’s decomposition method. Appl. Math. Comput. 2006, 173, 493–500. [Google Scholar] [CrossRef]
  23. Babolian, E.; Davari, A. Numerical implementation of adomian decomposition method for linear volterra integral equations of the second kind. Appl. Math. Comput. 2005, 165, 223–227. [Google Scholar] [CrossRef]
  24. Maleknejad, K.; Hadizadeh, M. A new computational method for volterra-fredholm integral equations. Comput. Math. Appl. 1999, 37, 1–8. [Google Scholar] [CrossRef] [Green Version]
  25. Abdulaziz, O.; Hashima, I.; Momani, S. Application of homotopy perturbation method to fractional ivps. J. Comput. Appl. Math. 2008, 216, 574–584. [Google Scholar] [CrossRef] [Green Version]
  26. Hu, Y.; Luo, Y.; Lua, Z. Analytical solution of the linear fractional differential equation by adomian decomposition method. J. Comput. Appl. Math. 2008, 215, 220–229. [Google Scholar] [CrossRef] [Green Version]
  27. Grzymkowski, R.; Pleszczyński, M.; Słota, D. Comparing the Adomian decomposition methodand the Runge–Kutta method for solutions of theStefan problem. Int. J. Comput. Math. 2006, 83, 409–417. [Google Scholar] [CrossRef]
  28. Grzymkowski, R.; Słota, D. Stefan problem solved by Adomian decomposition method. Int. J. Comput. Math. 2005, 82, 851–856. [Google Scholar] [CrossRef]
  29. Grzymkowski, R.; Słota, D. One—Phase Inverse Stefan Problem Solved Adomian Decomposition Method. Comput. Math. Appl. 2006, 51, 33–40. [Google Scholar] [CrossRef] [Green Version]
  30. Grzymkowski, R.; Pleszczyński, M.; Słota, D. The Two-Phase Stefan Problem Solved by the Adomian Decomposition Method. In Proceedings of the 15th IASTED International Conference Applied Simulation and Modelling, Rhodes, Greece, 26–28 June 2006. [Google Scholar]
  31. Verma, A.K.; Pandit, B.; Agarwal, R.P. Existence and nonexistence results for a class of fourth-order coupled singular boundary value problems arising in the theory of epitaxial growth. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  32. Hayani, W.A. Adomian decomposition method with green’s function for sixth order boundary value problems. Comput. Math. Appl. 2011, 61, 1567–1575. [Google Scholar] [CrossRef] [Green Version]
  33. Singh, R.; Kumar, J. The adomian decomposition method with green’s function for solving nonlinear singular boundary value problems. J. Appl. Math. Comput. 2014, 44, 397–416. [Google Scholar] [CrossRef]
  34. Noeiaghdam, S.; Sidorov, D.; Wazwaz, A.M.; Sidorov, N.; Sizikov, V. The Numerical Validation of the Adomian Decomposition Method for Solving Volterra Integral Equation with Discontinuous Kernels Using the CESTAC Method. Mathematics 2021, 9, 260. [Google Scholar] [CrossRef]
  35. Verma, A.K.; Pandit, B.; Verma, L.; Agarwal, R.P. A review on a class of second order nonlinear singular BVPs. Mathematics 2020, 8, 1045. [Google Scholar] [CrossRef]
  36. Verma, A.K.; Pandit, B.; Agarwal, R.P. Continuous Galerkin Method and Lane-Emden Equations. Adv. Math. Sci. Appl. 2021, 30, 39–64. [Google Scholar]
  37. Cabada, A. An overview on the lower and upper solutions method with nonlinear boundary value problem. Bound. Value Probl. 2011, 2011, 893753. [Google Scholar] [CrossRef] [Green Version]
  38. Ghorbani, A. Beyond adomian polynomials: He polynomials. Chaos Solitons Fractals 2009, 39, 1486–1492. [Google Scholar] [CrossRef]
  39. Singh, R.; Wazwaz, A.M. An efficient algorithm for solving coupled Lane-Emden boundary value problems in Catalytic diffusion reactions: The homotopy analysis method. Commun. Math. Comput. Chem. 2019, 81, 785–800. [Google Scholar]
Figure 1. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (12) for k = −1 and λ = 2 .
Figure 1. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (12) for k = −1 and λ = 2 .
Mathematics 09 00774 g001
Figure 2. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (13) for k = −1 and λ = 2 .
Figure 2. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (13) for k = −1 and λ = 2 .
Mathematics 09 00774 g002
Figure 3. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (14) for k = −1 and λ = 1 .
Figure 3. Approximate monotone iterations of Equations (52)–(55) corresponding to problem (14) for k = −1 and λ = 1 .
Mathematics 09 00774 g003
Table 1. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (8).
Table 1. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (8).
Lower SolutionUpper Solution
λ 0 31.94 0 31.94
0 1.95399 × 10 14 8.3347 × 10 7 1.86517 × 10 14
Table 2. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (8).
Table 2. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (8).
Positive SolutionNegative Solution
λ 1 15 1 15
1.42341 × 10 6 0.001199979 1.37856 × 10 16 3.66374 × 10 15
Table 3. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (9).
Table 3. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (9).
Lower SolutionUpper Solution
λ 0 11.34 0 11.34
0 3.55271 × 10 15 3.55271 × 10 15 2.66454 × 10 15
Table 4. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (9).
Table 4. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (9).
Positive SolutionNegative Solution
λ 1 15 1 15
6.83897 × 10 14 4.59188 × 10 13 2.54571 × 10 16 7.77156 × 10 16
Table 5. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (7).
Table 5. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (7).
Lower SolutionUpper Solution
λ 0 168.76 0 168.76
0 6.1668 × 10 11 0.000557377 1.70296 × 10 10
Table 6. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (7).
Table 6. Maximum absolute residue error of approximate solutions ϕ ( r ) corresponding to boundary conditions (7).
Positive SolutionNegative Solution
λ 1 15 1 15
0.0005578320.000562261 7.29668 × 10 17 6.2797 × 10 16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Verma, A.K.; Pandit, B.; Agarwal, R.P. Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth. Mathematics 2021, 9, 774. https://doi.org/10.3390/math9070774

AMA Style

Verma AK, Pandit B, Agarwal RP. Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth. Mathematics. 2021; 9(7):774. https://doi.org/10.3390/math9070774

Chicago/Turabian Style

Verma, Amit K, Biswajit Pandit, and Ravi P. Agarwal. 2021. "Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth" Mathematics 9, no. 7: 774. https://doi.org/10.3390/math9070774

APA Style

Verma, A. K., Pandit, B., & Agarwal, R. P. (2021). Analysis and Computation of Solutions for a Class of Nonlinear SBVPs Arising in Epitaxial Growth. Mathematics, 9(7), 774. https://doi.org/10.3390/math9070774

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop