Using Lie Symmetry Analysis to Solve a Problem That Models Mass Transfer from a Horizontal Flat Plate

We use Lie symmetry analysis to solve a boundary value problem that arises in chemical engineering, namely, mass transfer during the contact of a solid slab with an overhead flowing fluid. This problem was earlier tackled using Adomian decomposition method Fatoorehchi and Abolghasemi 2011 , leading to the Adomian series form of solution. It turns out that the application of Lie group analysis yields an elegant form of the solution. After introducing the governing mathematical model and some preliminaries of Lie symmetry analysis, we compute the Lie point symmetries admitted by the governing equation and use these to construct the desired solution as an invariant solution.


Introduction
A simplified model for mass transfer phenomenon from a horizontal flat plate fixed along a laminar fluid flow is considered 1 .Imagine a solid slab placed horizontally on the x-axis as shown in Figure 1.The free stream velocity is denoted by u ∞ m/s , C A0 is the initial concentration of molecules of species A mol/m 3 , C Ai is the concentration of molecules of species A at the plate interface mol/m 3 , δ is the momentum or hydrodynamic boundary layer thickness m , and δ c is the concentration boundary layer thickness m .Molecules of species A from the solid slab diffuse along the y-axis only to be swept downstream by the fluid flow in the hydrodynamic boundary layer.Subject to a number of assumptions given in 1 , mass balance for species A over an infinitesimal element based on the Cartesian coordinates leads to the PDE √ xC yy αyC x , 1.1 for the concentration C of molecules of species A, with the associated boundary conditions where δ x √ x/β, and α and β are constants.Precisely where D is the diffusion coefficient m 2 /s and v is the kinematic viscosity m 2 /s .The rest of the paper is organised as follows.In Section 2 we present preliminaries of Lie symmetry analysis.Section 3 is the thrust of the paper.In this section we determine Lie point symmetries admitted by the governing PDE and subsequently the basis generator of the one-dimensional Lie symmetry algebra admitted by the boundary value problem BVP .We then construct the solution to the BVP as an invariant solution.We give concluding remarks in Section 4.

Lie Symmetry Analysis of Differential Equations
There are many good introductions to Lie symmetry analysis of differential equations 2-6 .We will however indulge in a little introduction of the subject for completeness and to whet the appetite of the readers who could be new to the subject.Consider a scalar PDE of order k, with n independent variables, x x 1 , x 2 , . . ., x n , and one dependent variable u, where u j represents all the jth order partial derivatives of u with respect to x.An element of the set represented by u j is denoted by u i 1 i 2 ,...,i j ∂u/∂x i 1 ∂x i 2 , . . ., ∂x i j , i j 1, 2, . . ., n for j 1, 2, . . ., k.We want to define invariance of 2.1 under a one-parameter Lie group of transformations in the parameter ε 2.2 i 1, 2, . . ., n, acting on the x, u -space.The group 2.2 has as its infinitesimal generator where and the kth extension given by where i 1, 2, . . ., n and i j 1, 2, . . ., n for j 1, 2, . . ., k.The kth extension of the group 2.2 is generated by the following kth extended infinitesimal generator: with the explicit formulas for the extended infinitesimals given recursively by where D i is the total derivative operator defined by with summation over a repeated index.We say that 2.1 admits or is invariant under the Lie group of point transformations 2.2 if 2.1 has the same form in the new variables x x 1 , x 2 , . . ., x n and u, that is, When this is the case we sometimes loosely simply say that X the infinitesimal generator defined in 2.3 is a symmetry of 2.1 .Invariance of a differential equation under a given Lie group of transformations is neatly characterised by the infinitesimal criterion.Equation 2.1 is invariant under the Lie group of transformations 2.2 if and only if The infinitesimal criterion 2.10 provides a key to the explicit determination of symmetry groups admitted by differential equations.Using a straightforward algorithm based on 2.10 one obtains infinitesimals of the Lie group of point transformations that leaves a given differential equation invariant.The tedious algebraic calculations involved in this process are today done easily, often automatically, thanks to powerful Computer Algebra Systems CAS like Maple and Mathematica, Maxima and Reduce, and the many specific packages for performing symmetry analysis of differential equations 7-13 .An admissible Lie group characterises symmetry properties of a differential equation and is used for, among other things, complete integration in the case of ODEs or construction of special exact solutions of the differential equation see, e.g., 14-17 .For boundary initial value problems BIVPs there is a principle, the invariance principle, which states that if a BIVP is invariant under a given group, then one should seek the solution to the problem among the functions invariant under the admitted group.To state this principle concretely consider a BIVP F x, u, u 1 , u 2 , . . ., u k 0, 2.11 where s is a given manifold.Suppose that 2.11 admits m one-parameter symmetries X i ξ j x, u ∂ x j η x, u ∂ u , i 1, . . ., m.

2.13
We say that 2.11 -2.12 is invariant under a symmetry for some constants ε i , provided that i the manifold s is invariant under X; ii the boundary initial condition u| s h x is invariant under X restricted to s.
We now state the invariance principle 15, 18 : if the BIVP 2.11 -2.12 admits oneparameter symmetries X μ , one should seek the solution of the problem among the functions invariant under X μ .Successful application of the invariance principle to solve BIVPs has been reported in a number of papers see, e.g., 17, 19-21 .

Solution of the BVP 1.1 -1.2
Using LIE 8, 13 we obtain that 1.1 admits an infinite dimensional Lie symmetry algebra spanned by the following operators: where φ x, y is an arbitrary solution of 1.1 .
In the light of the invariance principle, the starting point in the "search" for the solution of 1.1 -1.2 is the determination of the subalgebra of the Lie algebra spanned by the symmetries in 3.1 that leaves the boundary conditions 1.2 invariant.To do this we construct a special linear combination, of the symmetries in 3.1 , with the constants ε i 's suitably chosen so that the boundary conditions 1.2 are invariant under 3.3 , that is, we require that X C − C A0 0 when y δ x , X C − C Ai 0 when y 0.

3.4
From 3.4 we obtain that from which it follows that ε 3 ε 4 0. Therefore 3.3 reduces to the general form of the symmetry admitted by both the PDE 1.1 and the boundary conditions 1.2 , where κ is an arbitrary constant.We can now assume the existence of an invariant solution C φ x, y of the BVP 1.1 -1.2 , the form of which is determined by functions invariant under the group generated by 3.6 .The characteristic equations of 3.6 , Mathematical Problems in Engineering yield the integrals from the first and second equations, respectively.These integrals are the similarity variables of the BVP and the invariant solution is expressed in terms of them in the form where φ is an arbitrary function.Upon substituting 3.10 into 1.1 we obtain the following second-order linear ODE: 3zφ z 2 αz φ z 0.

3.11
The solution of 3.11 is where K 1 and K 2 are arbitrary constants and Γ a, z is the incomplete gamma function defined by Γ a, z ∞ z t a−1 e −t dt.

3.13
We now have that all invariant solutions of 1.1 arising from 3.6 are embedded in the family of solutions 3.10 .The values of K 1 and K 2 in 3.12 are determined from the boundary conditions 1.2 : 3.14 where Mathematical Problems in Engineering 7 We now solve the equations 3.14 simultaneously for K 1 and K 2 .We obtain that K 1 and K 2 are constants only if we set κ 0, 3.16 in which case

3.17
Finally from 3.10 , 3.12 , and 3.17 we have that the BVP 1.1 -1.2 is solved by

Concluding Remarks
When a PDE is richly endowed with Lie point symmetries i.e., when the PDE admits the Lie symmetry algebras sl 2, R s A 1 or sl 2, R s W, where W is the Heisenberg-Weyl 19 the symmetry analysis approach usually provides a neat solution when it exists to the associated BIVP.Algorithms for obtaining admitted Lie point symmetries and for using the symmetries to construct invariant solutions are well developed and quite routine.It is instructive to compare the solution obtained in this paper, 3.18 , with the series-form solution obtained by Fatoorehchi and Abolghasemi 1 via Adomian decomposition: where w i 's and λ are "suitable" constants.As reported in 1 one can also solve 1.1 -1.2 by the combination of variables method, the starting point of which is the assumption that C x, y ψ y n x m 4.2 for suitable numbers m and n and some function ψ.This approach also leads to the exact solution

4.3
We remark here that Fatoorehchi and Abolghasemi 1 only determine the constants w i 's and λ by comparing 4.1 with the exact solution 4.3 .These constants are determined to be λ 4.64/ 3 √ 9.28 Γ 4/3 Sc 1/3 , w i 2 i , 4.4 where Sc is the Schmidt number.Finally, it is noteworthy that in the combination of variables approach to the solution of 1.1 -1.2 it turns out that m −1/2 and n 1 1 .This makes 3.10 and 4.2 equivalent in the light of 3.16 .

Figure 1 :
Figure 1: Molecules of species A from the solid slab diffuse along the y-axis only to be swept downstream by fluid flow in the hydrodynamic boundary layer.