We consider a quasilinear singularly perturbed two-point boundary value problem:
where, is a sufficiently small positive perturbation coefficient, b and c are sufficiently smooth functions, and for any satisfies , and , respectively where , , and C is a generic positive constant, independent of the perturbation parameter . In this case, the above problem contains a unique solution u, and u has an exponential boundary layer at (see   ).
Linß considered the problem (1) with in . The problem was discretized using a simple upwind finite difference scheme on adapted meshes using grid equidistribution of monitor functions. Kopteva and Stynes considered the same case in . In , Jugal and Srinivasan considered a linear problem and proved the uniform first-order convergence of the numerical solution on adapted meshes by equidistributing the arc-length monitor function. In , Linß, Roos and Vulanović considered the problem (1), which was discretized using a nonstandard upwinded first-order difference scheme on generalized Shishkin-type mesh. Linß considered the problem (1) on Bakhvalov-type mesh in , on BS-type mesh in . Zheng, Li, and Gao proposed a class of hybrid difference schemes with variable weights on Bakhvalov?Shishkin mesh in  and the parameter-uniform second-order convergence of approximating to the solution and the derivative on BS-type mesh and that of nearly second-order on Shishkin mesh were proved.
In this paper, we consider the problem (1) in the general case. The problem is discretized using a simple upwind finite difference scheme on adapted meshes using grid equidistribution of the arc-length monitor function. We will prove the uniform first-order convergence and the numerical experiments verify our analysis.
2. Simple Upwind Scheme and a Posteriori Error Estimate
So, the problem (1) is equivalent to:
For any nodes on interval [0, 1], , we set:
The simple upwind scheme is:
is the solution computed on the mesh .
For any grid function , we define the norm:
Easily get: , and .
Define the linear case of the problem (1):
where v and w be two arbitrary functions.
To simplify some arguments, we assume that:
The operators L and T are related as follow:
Lemma 1: Suppose that operator L satisfies inequality (7), then equation (5) has a unique solution and .
Lemma 2: Let , and , , , then .
Lemma 3: There exists a constant C such that
From Lemma 3, we get:
The approach to prove the above lemmas is similar to that in , where with was derived.
3. First-Order Uniform Convergence
The adaptive mesh method based on the principle of equal arc-length distribution is applied here.
Define the discrete operator of operator L by
Let denote the operator adjoint to with respect to , i.e.,
Let denote the discrete Green’s function associated with the mesh node , i.e.,
We have , which can be proved by discrete comparison principle.
For any mesh function we have the representation
It’s easy to prove that is an M-matrix. It follows that is an increasing function of i, for . Similarly is a decreasing function of i, for .
Thus, for each we have
where is the arc length of discrete solution .
Set , and , , , when
For one has
Thus, from (10) (11), we have
Because the adaptive mesh is equally distributed according to the arc length of the solution, so the arc length from to is .
Theorem 4: Set is the solution of (1) under the grid that satisfies arc-length equidistribution, then: .
Proof: Because and satisfies arc-length equidistribution, so: .
According to Lemma 3 and (9) (10), we have
4. Numerical Experiment
Example 1 (see ): We study the following quasilinear singularly perturbed two-point boundary value problem:
The exact solution of this problem is:
For any grid function we define the following discrete maximum norm:
In order to verify the error estimate, the numerical convergence rate of the solution is calculated as follows:
where, we use the maximum node error approximately replaces , is the maximum nodal error of the numerical solution obtained from the initial mesh. Under the simple upwind difference scheme, the results of the adaptive method with the Bakhvalov-Shishkin grid as the initial grid are shown in Table 1.
From the numerical example, we can see that the adaptive mesh can converge faster and be more accurate when the number of grid points is small in Figure 1, and the convergence result reaches the first-order uniform convergence in Figure 2.
In this paper, we consider the general problem (1). A simple upwind difference scheme is used to discretize on adaptive meshs based on grid equidistribution of the arc-length monitor function. The uniform first-order convergence of the algorithm is proved, and numerical experiments verify our analysis.
Figure 1. The maximum node errors of numerical solutions in Example 1 with ε = 10−5.
Figure 2. Rate of Example 1 with ε = 10−5.
Table 1. The maximum node errors of numerical solutions in Example 1.
The authors thank the support of Natural Science Foundation of China (No. 11471019).
 Tollmien, W., Schlichting, H., Grtler, H. and Riegels, F.W. (1961) über Flüssigkeitsbewegung bei sehr kleiner Reibung. Ludwig Prandtl Gesammelte Abhandlungen. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-662-11836-8
 Kopteva, N.V. (1997) On the Uniform in Small Parameter Convergence of a Weighted Scheme for the One-Dimensional Time-Dependent Convection-Diffusion Equation. Computational Mathematics and Mathematical Physics, 37, 1173-1180.
 Kopteva, N. and Stynes, M. (2001) A Robust Adaptive Method for a Quasi-Linear One-Dimensional Convection-Diffusion Problem. SIAM Journal on Numerical Analysis, 39, 1446-1447. https://doi.org/10.1137/S003614290138471X
 Jugal, M. and Srinivasan, N. (2010) Parameter-Uniform Numerical Method for Global Solution and Global Normalized Flux of Singularly Perturbed Boundary Value Problems Using Grid Equi-distribution. Computers & Mathematics with Applications, 60, 1924-1939. https://doi.org/10.1016/j.camwa.2010.07.026
 Linβ, T., Roos, H.G. and Vulanovic, R. (2000) Uniform Pointwise Convergence on Shishkin-Type Meshes for Quasi-Linear Convection-Diffusion Problems. SIAM Journal on Numerical Analysis, 38, 897-912. https://doi.org/10.1137/S0036142999355957
 Zheng, Q., Li, X. and Gao, Y. (2015) Uniformly Convergent Hybrid Schemes for Solutions and Derivatives in Quasilinear Singularly Perturbed BVPs. Applied Numerical Mathematics, 91, 46-59. https://doi.org/10.1016/j.apnum.2014.12.010
 Vulanovic, R. and Hovhannisyan, G. (2006) A Posteriori Error Estimates for One- Dimensional Convection-Diffussion Problems. Computers & Mathematics with Applications, 51, 915-926. https://doi.org/10.1016/j.camwa.2005.11.027