International Journal of Engineering MathematicsVolume 2014 (2014), Article ID 485431, 12 pageshttp://dx.doi.org/10.1155/2014/485431

Research Article

Process Parameter Identification in Thin Film Flows Driven by a Stretching Surface

Satyananda Panda,¹ Mathieu Sellier,² M. C. S. Fernando,³ and M. K. Abeyratne³

¹Department of Mathematics, National Institute of Technology Calicut, Kerala 673601, India²Department of Mechanical Engineering, University of Canterbury, Private Bag 4800, Christchurch 8140, New Zealand³Department of Mathematics, University of Ruhuna, 81000 Matara, Sri Lanka

Received 25 January 2014; Revised 10 June 2014; Accepted 12 June 2014; Published 21 July 2014

Academic Editor: Francisco Chinesta

Copyright © 2014 Satyananda Panda et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The flow of a thin liquid film over a heated stretching surface is considered in this study. Due to a potential nonuniform temperature distribution on the stretching sheet, a temperature gradient occurs in the fluid which produces surface tension gradient at the free surface of the thin film. As a result, the free surface deforms and these deformations are advected by the flow in the stretching direction. This work focuses on the inverse problem of reconstructing the sheet temperature distribution and the sheet stretch rate from observed free surface variations. This work builds on the analysis of Santra and Dandapat (2009) who, based on the long-wave expansion of the Navier-Stokes equations, formulate a partial differential equation which describes the evolution of the thickness of a film over a nonisothermal stretched surface. In this work, we show that after algebraic manipulation of a discrete form of the governing equations, it is possible to reconstruct either the unknown temperature field on the sheet and hence the resulting heat transfer or the stretching rate of the underlying surface. We illustrate the proposed methodology and test its applicability on a range of test problems.

1. Introduction

The analysis of thin film flow and heat transfer over a stretching surface has been a subject of fundamental importance as it is relevant to several industrial applications such as metal and polymer extrusion, continuous casting, drawing of plastic sheets, or cable coatings to name a few. This industrial context has drawn fluid dynamists and applied mathematicians alike to study this problem from a more canonical angle. The first important contribution to the understanding of this problem is the work of Wang [1] who formulated a mathematical model and developed a solution strategy based on the homotopy analysis method (HAM). The problem has been revisited several times since this seminal work with the inclusion of additional physics or more complex rheology. Andersson et al., for example, extended Wang’s contribution by analyzing the associated heat transfer problem [2, 3], while Noor and Hashim built on the work of Dandapat et al. [4, 5] to consider the thermocapillary and magnetic field effects [6] and Aziz and Hashim that of viscous dissipation [7]. Khan et al. focused on the effect of the temperature-dependency on the viscosity and thermal conductivity on the flow in the film [8] and Andersson et al. extended the standard formulation to power-law fluids [9]. A common feature of the literature cited above is that it is implicitly assumed that the film thickness is uniform in the domain, a required assumption to enable the similarity transformation which reduces the set of partial differential equations to a more tractable one of ordinary differential equations.

Recognizing the restrictions of the plane interface assumption, Dandapat and co-workers were the first to extend the formulation to account for local deformation of the free surface in [10, 11]. The authors exploit the slenderness of the flow domain to derive a long-wave approximation of the Navier-Stokes equations and solve the resulting governing equation using the matched asymptotic method. Lately, this work was extended to include the heat transfer problem [12]. In this work, the nonuniform temperature distribution at the stretching sheet induces an inhomogeneous temperature field in the film. Consequently, a surface temperature gradient develops at the film free surface. As a result of the surface tension gradients, the film thickness varies along the flow and these deformations are advected in the stretching direction.

This work focuses on the flow of a thin liquid film over a heated stretching surface. The aforementioned literature provides solid modeling foundations upon which one can build to indirectly infer process parameters difficult to measure in practice. More specifically, since the developed models correlate the film thickness to the sheet stretch rate and/or temperature through a set of differential equations, it is only natural to wonder whether the knowledge of the film thickness variation allows the reconstruction of either the sheet temperature or the sheet stretch rate. It is precisely this problem which this paper tackles. The motivation for this work is to demonstrate theoretically that it is possible to infer either the stretching rate or the surface temperature from the knowledge of the free surface evolution. This could have interesting practical outcome such as indirectly inferring the heat transfer involved in the process. The solution strategy to solve this inverse problem is inspired from the recent work of Sellier and Panda [13–15], which proposed a strategy to reconstruct an unknown field from free surface data in thin film flows. The general idea is to algebraically manipulate the field equation in order to obtain an explicit differential equation governing the inverse problem. This solution can then be solved using standard numerical techniques. This concept is applied for the first time here to a transient problem. The simplicity of the proposed approach contrasts with the traditional way of dealing with such inverse problems which rely on pde-constrained optimization framework, [16].

The following section briefly describes the mathematical model and numerical solution procedure used to compute the free surface evolution. This section is followed by the development of the inverse problem solution methodology. Some examples illustrating the success of the proposed approach are presented in the penultimate section and finally concluding remarks are drawn in the final section.

2. Description of the Direct Problem

2.1. Mathematical Model

We consider here an unsteady liquid film which lies over a horizontal plane heated elastic sheet as shown in Figure 1. Because surface tension tends to be a monotonically decreasing function of temperature, temperature gradients at the stretching sheet surface induce surface tension variations at the free surface. As regions of high surface tension tend to “pull” on regions of lower surface tension, convective flow cells develop which deform the free surface. This deformation is advected by the flow in the stretching direction. The question this work addresses is how the temperature field at the sheet surface and the stretching velocity can be reconstructed from a known deformation of the free surface.

Figure 1: Sketch of flow geometry.

This study is based on the model developed by Santra and Dandapat [12] using the long wave theory. This model is briefly described here for completeness. The elastic sheet lies at and the liquid gas interface lies at , where the axis is directed along the stretching sheet and the axis is normal to the sheet in the outward direction toward the fluid. Gravity acts along negative direction. Further, the surface at starts stretching from rest and within a very short time attains the stretching velocity . The elastic surface is heated with its temperature a function of alone and the ambient gas phase is at constant temperature .

The fluid is assumed to be incompressible and Newtonian with constant kinematic viscosity , density , specific heat , and thermal conductivity . The surface tension of the liquid-gas interface decreases linearly with temperature according to where is the inter-facial temperature, is the surface tension at , and is a positive constant specific to the fluid.

Moreover, Newton’s law of cooling describes the convective heat flux at the interface by where is the temperature in the fluid, is the unit normal vector on the interface, and is the rate of heat transfer from the liquid to the ambient gas phase.

Following the analysis of Santra and Dandapat [12] based on the the lubrication equation, the dimensionless evolution equation for film thickness is given by where is the flow discharge given by Here the subscript means partial differentiation with respect to . The dimensionless numbers are(i)the Froude number that expresses the ratio of inertia to body forces;(ii)the Prandtl number that represents the ratio of momentum to thermal diffusivity;(iii)the Marangoni number that characterizes the relation between the temperature dependent surface tension and viscous forces;(iv)the Biot number that compares the relative magnitudes of resistances to internal conduction and surface convection;(v)the dimensionless number is known as surface tension parameter and is the aspect ratio.In the above expressions, and are characteristic length scales in the vertical and horizontal directions, respectively. The symbol stands for the sheet temperature at the origin.

At the origin, we apply the symmetry conditions and and at the other end of the domain; we assume that the same sheet temperature profiles and the sheet stretching rate continues beyond the computed domain. We also assume that the gradient of the free surface extends out of the computational domain. These boundary conditions are consistent with those mentioned in [12].

2.2. Numerical Solution Procedure

There are many different numerical methods that have been presented in the past for this kind of analysis [12, 17] and in this paper, we mostly follow the finite volume technique described in Sellier and Panda [14] on a uniform grid system with implicit flux discretization.

We discretize our flow domain into the uniform grid and the flow variable , the steady stretching velocity , and sheet temperature are located at the cell centers as shown in Figure 2. Let spatial grid points be equally distributed over the domain with spatial increment , . The numerical solution is sought at the discrete time levels , with time step .

Figure 2: Typical grid used for the finite volume discretizations. Blacked dots and dashed vertical lines represent the nodal points and the cell faces, respectively.

To set up the discrete equation, the function is approximated over the cell . That is the cell average of over is denoted as

For the given solution , the solution at the next time level is obtained by integrating (3) over the space and time intervals . This yields the following discrete equation: for nodes , where the discrete flux function is given by where face values are evaluated using linear interpolation from nodal values and gradients using forward differences. For example, Similar expressions can easily be derived for other terms. The following second-order accurate approximation of the third-order derivative is used:

The first and the last nodes are located at the boundary and special treatment is required there. Therefore, we apply the boundary conditions in the first node and at the last nodes. In brief, at the first node we have where the value of the functions , , , and are known at the cell face . Similarly, we approximate the discrete values at the last node.

Equation (6) describes an implicit time discretization scheme. Since the governing equation is nonlinear, a system of nonlinear algebraic equations needs to be solved at each time step. We use fsolve in MATLAB for this purpose. A good initial starting guess is required to solve the nonlinear equations. A reasonable initial guess for the free surface is chosen to be unity throughout the discrete domain at the first time step. The solution from the previous time step can be used otherwise. Convergence is usually achieved in less than 10 iterations and the convergence criterion is that the norm of the residuals should be less than .

To demonstrate the successful implementation of the proposed finite volume algorithm, the numerical results obtained with the proposed algorithm are compared to those obtained by Santra and Dandapat, [12]. Specifically, the results of Figure 7 in [12] are reproduced here for validation purpose. For this benchmark case, the following temperature distribution and stretching velocity profiles are imposed: The temperature profile is illustrated in Figure 3(a). The other nondimensional parameters used in the simulation are reported in the figure caption. The simulation was performed using 150 grid points and time steps of (dimensionless units).

(a)

(b)

Figure 3: Results for the temperature and stretching velocity profile given by (11) and with , , , , , and : (a) temperature profile; (b) evolution of free surface profile at different time from to with the increment of time step 1.

Figure 3(b) illustrates the evolution of free surface at different times. The observed results replicate well those obtained by Santra and Dandapat in [12], and thereby validating the formulation and implementation of the proposed numerical scheme. It is clear from Figure 3(b) that thermocapillary deforms the free surface during the early stages and this deformation is advected downstream by the stretching sheet. To clearly see the motion of the free surface disturbance, the local maximum on each curve is identified with star.

The effect of stretching velocity distribution on film height is demonstrated next. The following Gaussian temperature profile is imposed at the sheet and three different stretching velocity distributions are considered: The free surface is plotted in Figure 4(b) for the different stretching velocities (illustrated in Figure 4(a)) at a fixed time of and the parameters reported in the figure caption. Visual inspection confirms that the computed free surface profiles match well with those computed by Santra and Dandapat in Figure 3 of [12] further demonstrating the accuracy of the proposed numerical approach.

(a)

(b)

Figure 4: Results for the different stretching velocity profile with temperature profile given by (12) and with , , , , , and : (a) stretching velocity profile; (b) free surface profile.

3. The Inverse Problem

The evolution of the free surface provides valuable information on the flow process parameters. Such information can be used to infer unknown conditions otherwise difficult to measure. In the following, we will consider two inverse problems which can be solved using the same conceptual idea.

3.1. Temperature Reconstruction from Two Free Surface Snapshots and Prescribed Stretching Velocity

This first inverse problem consists in inferring the temperature of the stretching sheet from known free surface deformation data and a known stretching velocity distribution over the sheet. The inverse counterpart of the above problem is discussed assuming the transient solution for exist. More specifically, it is assumed that two free surface profiles at two successive time steps are available. In the following, we reconstruct the sheet temperature distribution from the film thickness profiles obtained by solving the direct problem up to a specific time level with a prescribed stretching velocity. The main idea behind the reconstruction strategy is that the governing equations, (3) and (4), hold equally well for the forward or the inverse problem. Thus, we can obtain an explicit partial differential equation which “governs” what the temperature distribution needs to be to result in a given film thickness distribution variation. Therefore, unlike traditional inverse problem methods which tend to rely on recasting the reconstruction problem into an optimization one whereby one has to infer the “optimal” temperature distribution which minimizes a given functional, the proposed methodology only involves finding the solution of a partial differential equation. Like the direct problem, we make use of the the same implicit discretized equation (6). Recall the equation as This nonlinear discrete equation could equally be written as where the bold font indicates the vector of discrete values. The explicit form of the operator is lengthy and not reported here. With the governing equation in this form, it becomes apparent that this equation describes three possible problems. The usual forward problem, as it is usually described, is recovered if we assume that , , and are known and is sought. The two inverse problems of interest here are recovered if(i), , and are known and is sought;(ii), , and are known and is sought;Irrespective of the inverse problem considered, (15) is an implicit equation for the temperature distribution (given the stretching rate), or the stretching rate (given the temperature distribution). Boundary conditions similar to the ones used for the forward problem are implemented in the inverse problem. Equation (15) is solved implicitly in the scientific computing program MATLAB using the fsolve routine for the reconstruction of from the knowledge of the film thickness distribution at two successive time steps and the temperature at the end points of the flow domain. An initial guess of zero temperature distribution, that is, , is considered for the implicit solver.

For the first test case, the artificial “experimental” free surface data is obtained for the prescribed linear stretching velocity and the following sine temperature profile: Additional noise with a standard uniform distribution is added to the obtained film thickness profile to replicate the likely measurement uncertainty. The amplitude of the added noise was either 3% of the total film thickness variation or 5%. Without any prior processing of the noisy data, the reconstruction algorithm invariably failed. However, the algorithm was successful if a cubic spline was first fitted through the noisy film thickness variation. The MATLAB function csaps() was used for this purpose with a smoothing parameter set to 0.9, [18]. The nondimensional parameters used in the simulation are reported in the figure caption. The number of grid points is 150 and time steps of 0.03 are used in the simulation. The corresponding free surface profile with and without noise at is shown in Figure 5(a). Figure 5(b) shows the reconstructed and the actual temperature profiles. Without noise, these two profiles are indistinguishable thereby validating the conceptual reconstruction idea and the implementation of the algorithm. Given the smoothness of the actual temperature profile, it comes as little surprise that even with the addition of 3% or 5% noise on the free surface data, the reconstructed temperature profile still is in good agreement with the actual one. This highlights the robustness of the proposed methodology even in the presence of noise.

(a)

(b)

Figure 5: Results for the sine sheet temperature profile, that is, , and with , , , , , , and stretching velocity : (a) free surface profiles with and without noise; (b) actual and reconstructed temperature profiles.

The ability of the proposed algorithm to reconstruct temperature profiles with a much steeper variation is tested in the next test case. The artificial temperature profile given by a smoothed step temperature function, that is, As in the previous case, 3% and 5% noise was then added to the film thickness profiles. Cubic splines were then fitted through these free surface profiles and used as input to the reconstruction algorithm. The resulting free surface profiles with and without noise at can be seen in Figure 6(a). It is apparent from these profiles that the film thickness variations occur on a much shorter length scale which will lead to a more difficult test case for the reconstruction algorithm, as one would expect. Reconstruction results can be seen on Figure 6(b) without noise, for 2% added noise, and for 5% added noise. Without noise, it is apparent that the actual and reconstructed temperature profiles are strictly identical again confirming that it is possible to reconstruct the surface temperature from the knowledge of two successive free surface profiles and the stretching velocity. It is also clear that the actual and reconstructed temperature profiles are still in reasonable agreement when noise is added to the input data. The reconstruction algorithm only mildly amplifies the input noise. The quality of the reconstruction for this second case is not as good as for the first one and one can attribute this to the difficulty of fitting a good representative cubic spline through rapidly varying data. Hence, one would expect this reconstruction algorithm to work best when the unknown temperature profile does not vary too rapidly.

(a)

(b)

Figure 6: Results for the smooth step sheet temperature profile given in (11) and with , , , , , , and stretching velocity : (a) free surface profiles; (b) actual and reconstructed temperature profiles with or without added noise on the input data.

In the inevitable presence of noise in an experimental system, the natural procedure is to run several realizations of the same experiment and average the output data. This is precisely what we proceeded to do in order to check whether, as the number of realization increases, the agreement between the reconstructed and actual temperature profile improves, as one would expect. Results are reported in Figure 7 for the 5% noise case for which the agreement between the actual and reconstructed temperature profiles was the poorest. In effect, we ran numerous reconstruction with 5% noise added to the film thickness and averaged the reconstructed temperature profiles. Figure 7(a) shows two realizations of the reconstructed temperature profiles (Realizations 1 and 2) which are rather poor. However, averaging the reconstructed temperature profile over 50 realizations yields a very good agreement with the actual distribution, as anticipated. Figure 7(b) shows the norm of the difference between the actual and reconstructed temperature profiles. It is clear from this graph that the norm decreases very rapidly for the first few realizations with very marginal reduction of this norm between 10 to 50 realizations. This indicates that, in practice, few realizations of the experiment would be necessary in order to obtain a reliable reconstruction of the temperature distribution.

(a)

(b)

Figure 7: Effect of the number of realizations on the reconstructed temperature profile: (a) individual and averaged reconstructed temperature profiles; (b) norm of the error between the actual temperature profiles and the one reconstructed using multiple realizations.

To explore the effects of the Prandtl number (), the Biot number (), and the surface tension parameter () on the reconstruction algorithm the algorithm was run for four different scenarios where , , and are fixed and the Biot number, Prandtl number and surface tension parameter vary, that is, , , and . A nonlinear parabolic (concave) profile is considered for the stretching velocity that is, for all four cases. Considering Figures 8(b), 9(b), and 10(b), it can be seen that the reconstruction algorithm recovers the temperature profile extremely well without noise even if the imposed stretching velocity is nonlinear. For and added noise, the quality of the reconstruction is still reasonable. The corresponding free surface profiles can be seen in Figures 8(a), 9(a), and 10(a). Slightly larger discrepancies can be observed when the surface tension parameter is increased from to , see Figure 11(b).

(a)

(b)

Figure 8: Results for the temperature profile at time with , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed temperature profiles.

(a)

(b)

Figure 9: Results for the temperature profile at time with , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed temperature profiles.

(a)

(b)

Figure 10: Results for the temperature profile at time with , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed temperature profiles.

(a)

(b)

Figure 11: Results for the temperature profile at time with , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed temperature profiles.

3.2. Stretching Velocity Reconstruction from Two Free Surface Snapshots and Prescribed Temperature Distribution

It is clear from the work of [10–12] that the stretching velocity plays an important role in the evolution of the thin film height. In addition to the Marangoni effect, the type of stretching applied at the bottom sheet mainly influences the film profile. In the following analysis, we discuss how to reconstruct the stretching velocity from given free surface data and temperature distribution at the surface. For this purpose, we first solve the forward problem for the film thickness profile up to a particular time for a given temperature distribution and stretching sheet velocity. To reconstruct the stretching velocity profile, we consider the same governing equations: (3) in continuous form and (6) in discretized form. In order to solve the nonlinear algebraic equation (6) for the stretching velocity from the known two successive free surface profiles and the temperature function , we first prescribe the stretching sheet velocity at the boundary points and then solve the equation using MATLAB fsolve routine. It is observed that for any nonzero constant initial guess, the MATLAB solver converges within less than 10 iterations.

For the first case, the film thickness is obtained up to a particular time of for the stretching velocity , and temperature profile is given by (12). The other flow parameters used in the simulation are reported in the figure caption. The transient solution is obtained by solving (6), for 150 grid points and time steps of . The film thickness at is shown in Figure 12(a) without noise and with 3% and 5% noise. Figure 12(b) shows the reconstructed and actual stretching velocity profiles. In the absence of noise, the perfect agreement between the two profiles confirms that the reconstruction algorithm is able to recover the stretching velocity equally well as it could recover the substrate temperature profile. The addition of up to 5% noise still lead to a reasonable reconstruction.

(a)

(b)

Figure 12: Results for the nonlinear stretching velocity , temperature profile given by (12) and with , , , , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed stretching velocity profiles with or without added noise on the input data.

For the second case, the film thickness is obtained up to a particular time of for the nonmonotonic stretching velocity and temperature profile . The other flow parameters used in the simulation are reported in the figure caption. The simulation results for reconstructed surface velocity without surface tension () and with surface tension () are given in Figures 13(b) and 14(b), respectively. In the absence of surface tension the reconstruction algorithm is able to recover the stretching velocity profile. It can be seen that with the inclusion of surface tension the algorithm perfectly reconstruct the sheet velocity in the absence of noise (Figure 14(b)). With the addition of noise the reconstructed results differ mildly from the actual one.

(a)

(b)

Figure 13: Results for the nonmonotonic stretching velocity with temperature profile and , , , , , : (a) free surface profiles with and without noise; (b) actual and reconstructed stretching velocity profiles with or without added noise on the input data.

(a)

(b)

Figure 14: Results for the nonmonotonic stretching velocity with temperature profile , , , , , , and : (a) free surface profiles with and without noise; (b) actual and reconstructed stretching velocity profiles with or without added noise on the input data.

4. Concluding Remarks

A solution strategy is presented in this work to reconstruct the surface temperature profile or the surface stretching rate from the knowledge of the thickness at a given time and its rate of change. The main idea is simply based on a rearrangement of the equations derived by Santra and Dandapat in [12] for nonisothermal thin liquid films on stretching surfaces. This rearrangement of the governing equations leads to an explicit partial differential equation which governs the inverse problem. This allows to solve the inverse problem in “one shot” which is a considerable advantage compared to the traditional way of approaching such inverse problems based on pde-constrained optimization theory. We have demonstrated that the algorithm is quite robust in the sense that noise in the input data is not too amplified provided the input data is smoothed using a cubic spline. Whilst the focus of this paper is on the presentation of the conceptual idea, future work will assess the feasibility and practicality of the suggested approach using true experimental data.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

B. Wang, “Liquid film on an unsteady stretching surface,” Quarterly of Applied Mathematics, vol. 48, no. 4, pp. 601–610, 1990.

H. I. Andersson, J. B. Aarseth, and B. S. Dandapat, “Heat transfer in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 43, no. 1, pp. 69–74, 2000.

I.-C. Liu and H. I. Andersson, “Heat transfer in a liquid film on an unsteady stretching sheet,” International Journal of Thermal Sciences, vol. 47, no. 6, pp. 766–772, 2008.

B. S. Dandapat, B. Santra, and H. I. Andersson, “Thermocapillarity in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 46, no. 16, pp. 3009–3015, 2003.

B. S. Dandapat, B. Santra, and K. Vajravelu, “The effects of variable fluid properties and thermocapillarity on the flow of a thin film on an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 991–996, 2007.

N. F. M. Noor and I. Hashim, “Thermocapillarity and magnetic field effects in a thin liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 53, no. 9-10, pp. 2044–2051, 2010.

R. C. Aziz and I. Hashim, “Liquid film on unsteady stretching sheet with general surface temperature and viscous dissipation,” Chinese Physics Letters, vol. 27, no. 11, Article ID 110202, 2010.

Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers & Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011.

H. I. Andersson, J. B. Aarseth, N. Braud, and B. S. Dandapat, “Flow of a power-law fluid film on an unsteady stretching surface,” Journal of Non-Newtonian Fluid Mechanics, vol. 62, no. 1, pp. 1–8, 1996.

B. S. Dandapat, A. Kitamura, and B. Santra, “Transient film profile of thin liquid film flow on a stretching surface,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 57, no. 4, pp. 623–635, 2006.

B. S. Dandapat and S. Maity, “Flow of a thin liquid film on an unsteady stretching sheet,” Physics of Fluids, vol. 18, no. 10, Article ID 102101, 7 pages, 2006.

B. Santra and B. S. Dandapat, “Unsteady thin-film flow over a heated stretching sheet,” International Journal of Heat and Mass Transfer, vol. 52, pp. 1965–1970, 2009.

M. Sellier, “Substrate design or reconstruction from free surface data for thin film flows,” Physics of Fluids, vol. 20, no. 6, Article ID 062106, 2008.

M. Sellier and S. Panda, “Surface temperature reconstruction based on the thermocapillary effect,” The ANZIAM Journal, vol. 52, no. 2, pp. 146–159, 2010.

M. Sellier and S. Panda, “Inverse temperature reconstruction in thermocapillary-driven thin liquid films,” International Journal of Numerical Analysis and Modeling B, vol. 3, no. 3, pp. 285–296, 2012.

O. Pironneau, Optimal Shape Design for Elliptic Systems, Springer, New York, NY, USA, 1984.

Y. Ha, Y.-J. Kim, and T. G. Myers, “On the numerical solution of a driven thin film equation,” Journal of Computational Physics, vol. 227, no. 15, pp. 7246–7263, 2008.

The MathWorks Inc., MATLAB R2011b Documentation, The MathWorks Inc., 2011.