# Forward Roll Coating of a Williamson’s Material onto a Moving Web: A Theoretical Study

Authors: M. Zahid; I. Siddique; S. Saleem; A. Al-Zubaidi; M. A. Rana; M. Zafar
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

This paper presents a mathematical model for the thin film roll coating process of an incompressible Williamson material, passing through a closed passage between a rotating roll and a web. In light of lubrication approximation theory, the flow equations are nondimensionalized. The regular perturbation approach is used to provide solutions for the velocity profile, pressure gradient, flow rate per unit width, and shear stress at the roll surface. Important engineering quantities such as coating thickness, maximum pressure, separation point, roll/sheet separating force, and roll-transmitted power to the fluid are also obtained. The effects of several factors are graphically projected. The study shows that the material factors that are involved determine the operating variables. Coating thickness and separation point are controlled by Weissenberg’s number, therefore acting as a controlling parameter for the rate of flow, thickness in coating, power contribution, pressure, roll separating force, and separation point. In comparison to the existing results in the literature, the current results are broader and zero-order results are more accurate.


---

## Body

## 1. Introduction

Deposition of fluid to the substrate under consideration is performed through roll coating in an industrial process. The phenomenon had gained a healthy reputation during the past few decades. As an example, the paper coating, decoration, or protection of fabrics or metal with such materials, films, coated products and magnetic recording, thin uniform liquid coating are the examples of coating in industry. To carry out these operations, a variety of apparatuses are used, in which the most common are roll coaters. The flow of non-Newtonian materials has several industrial and physiological applications in many coatings processes. Moreover, such flows have caused various challenges for mathematical scientists and model developers in the development of procedures suitable for the calculation of flows. This is due to the complexity and higher-order nature of leading flow expressions of non-Newtonian fluids than the Navier–Stokes equations.Due to the complex nature of flows, numerous leading expressions for non-Newtonian flows are suggested in the literature. Abegunrin et al. [1] analyzed the boundary layer flow of two non-Newtonian fluids over an upper horizontal surface of a paraboloid of revolution. Abegunrin and Animasaun [2] studied the motion of non-Newtonian Williamson fluid over an upper horizontal surface of a paraboloid of revolution due to partial slip and buoyancy. Shah et al. [3] signified the water base nanoparticle’s radius, heat flux due to concentration gradient, and mass flux due to the temperature gradient. Animasaun and Omowaye [4] studied the upper-convected Maxwell fluid flow with variable thermophysical properties over a melting surface situated in a hot environment subject to thermal stratification. Stability analysis of convection non-Newtonian vertical fluid layer in the presence of gold nanoparticles is performed by Mekheimer et al. [5]. Zaher et al. [6] studied the residual time of sinusoidal metachronal ciliary flow of non-Newtonian fluid through ciliated walls. Koumy et al. [7] investigated the hall and porous boundaries effects on peristaltic transport through a porous medium of a Maxwell model. Abdelsalam et al. [8] theoretically investigated the physical traits of electromagnetohydrodynamics (EMHD) of the bloodstream in the presence of electroosmotic forces via arteries having both aneurysm and stenosis for a hybrid fractional second-grade nanofluid. Raza et al. [9] investigated the thermal transport of radiative Williamson fluid over the stretchable curved surface. Pseudoplastic fluids are considered the most common in non-Newtonian fluids. Studies of pseudoplastic fluids have gained excessive attention for their utmost utility in industrial processes like extruding of polymer films, melts, and solutions with high molecular weight polymers and films coated with emulsion such as photographic films. Due to the pseudoplastic nature of Williamson fluid [10–12], its model is being focused on. However, less consideration has been paid to the Williamson fluid model, which describes the properties of pseudoplastic fluids. Nadeem and Akbar [13] debated the peristaltic flow of a Williamson fluid. Cramer and Marchello [14] demonstrated with experimental evidence that this model is better suited to polymer solutions and particle suspensions than other models.For better understanding the procedure and to foresee the complications faced during the development, many hypothetical surveys have been undertaken on the flows under discussion. Model for flow conduct inside toner of the printing press by considering the concept of lubrication was developed by Taylor and Zeltlemoyer [15]. The values for force and pressure distribution were calculated by them. Hintermaier and White [16] developed the model for the water flow between two rolls by using lubrication theory and observed the results compatible with their conclusions carried out experimentally. Moreover, the flow theory of roll coating of Newtonian Liquid, assuming the traditional lubrication concepts and some general physical concepts was developed by Greener and Middleman [17]. Their model was modified by Benkreira et al. [18] to a general situation of two rollers of equal or unequalized rotating at the same speed or uneven one, consequently establishing a model related to the general situation. Their conclusions were found close to experimental results. Numerical outcomes on the mentioned problem over a moving flat subs rate were provided by Sofou [19] undertaking the lubrication approximation theory. They manipulated the Herschel–Bulkley model of viscoplasticity, thereby reducing it to the Newtonian models under suitable modifications, Bingham, and power-law. Recently, Zahid et al. [20, 21] undertook the same situation with the viscoelastic fluid being incompressible in which the roll, as well as web, was assumed to move with the same rates. To simplify expressions for motion, they performed the above-mentioned theory. Using the regular perturbation method, solutions for main flow, pressure gradient, flow rate per unit width, and shear stress at the roll surface are computed.To the best of our knowledge, no study has been highlighted on the said fluid under roll coating in the mentioned works. The focus of the present study is to establish the flow technique for the model under consideration and to examine the influence of fluid physical properties on the coating process. Among physical properties, the Weissenberg numberWe is used in the study of viscoelastic flows and it usually gives the relation of stress relaxation time of the fluid and a specific process time. Some interesting studies regarding Weissenberg numbers are given in [22, 23]. The paper has been ordered as below. In the upcoming section, the flow equations and the construction of the problem are defined. The succeeding part discusses the exact solution of the flow parameters. Lastly, outcomes, discussion, and deductions are offered.


## 2. Mathematical Formulation

Consider a steady, laminar flow of an incompressible Williamson fluid between a flat plate and cylindrical roll of radiusR as represented in Figure 1. The plate moves with a constant velocity U in a positive x-direction, as the roll rotates counterclockwise with angular velocity ω. The plate and cylinder have a gap at the nip, H0. The polymer first bites the plane at the location x=−xb, as shown in Figure 2. The expressions leading the flow of a said fluid are as follows:(1)∇⋅V=0,(2)ρdVdt=−∇p+∇⋅τ,where V is the velocity, ρ is the density, p is the pressure, d/dt⋅=∂/∂t⋅+V⋅∇⋅ is the material derivative, and τ is the extra stress tensor, for a Williamson fluid [24] defined as(3)τ=η∞+η0−η∞1−Γγ˜−1A1,where η∞ and η0 are correspondingly the infinite and zero shear rate viscosities, respectively, A1=∇V+∇Vt is the first Rivlin–Ericksen tensor, Γ is the time constant, and the shear rate γ˜ [5] is defined as(4)γ˜=12∑iγ˜ij∑jγ˜ji=12π,where π=traceA12 is the second invariant stress tensor. Constitutive equation (3) takes the following form, when we consider the case for which η∞=0 and Γγ˜<1:(5)τ=η01+Γγ˜A1.Figure 1
Physical representation of forward roll coating.Figure 2
Geometry of the studied physical model.The above model reduces to Newtonian whenΓ=0. The curved channel length molded by the plane and the roll is much larger than the parting at the nip, that is, H0≪R, thereby causing the flow to become two-dimensional. The velocity profile consequently takes the following form:(6)V=u11x,y,v11x,y.In view of equation (6), the flow equations (1) and (2) in components form can be written as(7)∂u11∂x+∂v11∂y=0,(8)ρu11∂u11∂x+u11∂v11∂y=−∂p11∂y+∂σxx∂x+∂τyx∂y,(9)ρu11∂v11∂x+v11∂v11∂y=−∂p11∂y+∂τxy∂x+∂σyy∂y,where σxx, τyx, τxy, and σyy are the stress components which can be written from equation (5) as(10)σxx=2η01+Γγ˜∂u11∂x,τxy=τyx=2η01+Γγ˜∂u11∂x+∂v11∂y,σyy=2η01+Γγ˜∂u11∂y.We start with the LAT assumptions that the most significant dynamic phenomena happen in the region where the minimum gap exists, that is, the nip region. In that region, and extending to either side (i.e., the±x direction) by a distance of the order of x0, the roll surfaces are nearly parallel. Then, it is reasonable to assume that v11≪u11 and ∂/∂x≪∂/∂y. In order to obtain the characteristic scale for the velocity and pressure, we conduct in brevity an order of magnitude analysis, and we can identify the following scales for x,y, and u11,x∼Lc,y∼H0,u11∼U, u11∼U. Then, from equation (1), we obtain vc/U∼H0/Lc≪1, which shows that the order of magnitude of a transversal velocity vc is smaller than the longitudinal velocity, where the longitudinal characteristic length is given by Lc=2RH0. Based on the above discussion, we introduce the following dimensionless quantities:(11)x¯=xRH0,y¯=yH0,u¯=u11U,v¯=v11ϑ1U,p¯=p11H0η0UH0R,γ˜¯=γ˜H0U,ϑ1=H0R,into equations (7)–(10), and we get the following nondimensional equations (after dropping (−) for simplicity):(12)∂u∂x+∂v∂y=0,(13)Reϑ1u∂u∂x+v∂u∂y=−∂p∂x+ϑ1∂σxx∂x+∂τxy∂y,(14)Reϑ13u∂v∂x+v∂v∂y=−∂p∂y+ϑ1∂σyy∂y+ϑ12∂τxy∂y,(15)σxx=2ϑ11+Weγ˜∂u∂x,τxy=τyx=1+Weγ˜∂u∂y+ϑ12∂v∂x,σyy=2ϑ11+Weγ˜∂v∂x,where Re=ρUH0/η0,We=ΓU/H0, and γ˜=2ϑ12∂u/∂x2+∂u/∂y+ϑ12∂v/∂x2+2ϑ12∂v/∂y21/2.Neglecting the terms havingϑ1 being small quantity as ϑ1 is a square root of the ratio of H0 to R, we get(16)∂∂y∂u∂y1+We∂u∂y=∂p∂x,(17)∂p∂y=0.It is obvious from equation (17) that p is not a function of y, so it is the function of x only; that is, p=px. Thus, equation (16) becomes(18)dpdx=∂2u∂y2+We∂∂y∂u∂y2.The corresponding nondimensional boundary conditions are(19)u=1,aty=0,u=1,aty=hx,,where hx=1+x2/2 represents the dimensionless gap between roll and the web. Along with equation (18), the following nondimensional balanced mass equation is needed:(20)λ=∫0hxuydy,where λ=Q/UWH0,Q is the dimensionless flow rate per unit width.


## 3. Solution of the Problem

Most of the mathematical problems do not possess explicit solutions, so it is very important to know methods for approximating the solution of such type of problems. Asymptotic analysis like the perturbation technique provides a powerful procedure for obtaining approximate solutions to complex problems. Such problems contain a parameter that is very large or very small, and one wants to take advantage of this structure to get the best approximation. These techniques are very useful in converging/diverging geometries like coating flow problems. The closed-form solution of equation (18) is difficult being a nonlinear differential equation. Therefore, regular perturbation is used to find its solution. We take We≪1 (as a perturbation parameter) and expand velocity u, pressure p, and dimensionless film thickness in power series of We:(21)u=u0+Weu1+We2u2+OWe3,(22)dpdx=dp0dx+Wedp1dx+We2dp2dx+OWe3,(23)λ=λ0+Weλ1+We2λ2+OWe3.Substituting expressions (21)–(23) into equations (18)–(20) and equating like powers of We, we obtain the following problems.


### 3.1. Zeroth-Order Problem and Its Solution

Comparing terms free fromWe in equation (18), we have the following zero-order problem:(24)dp0dx=d2udy2,(25)λ0=∫0hxu0dy,with the agreeing zeroth-order boundary conditions:(26)u0=1,aty=0,u0=1,aty=hx.The solution for equation (24) with conditions (26) is given by(27)u0x,y=1+12dp0dxy2−hy.The pressure distribution in equation (27) can be found by utilizing the standard of preservation of mass and has the accompanying structure:(28)dp0dx=12hx−λ0hx3.By using equation (28) into equation (27), one can write the zero-order velocity distribution as(29)u0x,y=1+6h−λ0h3y2−hy.At the separation point located atx=xs and y=h¯xs/2, as shown in Figure 2, the velocity goes to zero. Hence, from equation (29), we can write(30)xs=3λ0−1.Integrating equation (28), write in the following form:(31)p0=12∫xsxhx−λ0hx3dx,where xs is the value of x at p=0. Since we are assuming that the liquid splits evenly to the web and the roll, we can write Q0=2Uhc and in dimensionless form λ0=2hc/H0. The pressure at the separation point may be found by the first integral of equation (28), using assumption xb⟶−∞, as(32)p0x=6−9/2λ0x1+x2/2−3λ0x1+x2/22+122−9λ02tan−1x2+6π21−3λ04.The simplest dynamic model of the separation region is based on the assumption that the film splits at the point, whereu0=0andp0=0; resultantly, the above equation becomes a transcendental equation in λ0 whose solution may be obtained by a numerical technique named as the modified regula falsi method with a predefined accuracy of 10−10. For the said accuracy, the zero-order flow rate is λ0=1.3232524, which is correct up to five decimal places, and the residual error in this case is 0.0000011085; this coating thickness has been achieved at the separation point xs=1.723298349. It is worth to mention that the above zero-order results are more accurate as calculated by Middleman [17].


### 3.2. First-Order Problem and Its Solution

Equating the coefficients of termsWe in equation (18), the first-order problem is(33)dp1dx=d2u1dy2+ddydu0dy2,(34)Q1=λ1=∫0hxu1dy,with the agreeing first-order boundary conditions:(35)u1=0,aty=0,u1=0,aty=hx.Adopting the same procedure as in Section3.1, applying equation (27) into equation (33), after twice integration and using conditions (35), we acquire the following first-order solution:(36)u1x,y=12dp1dxy2−hy−16dp0dx22y3−3hy2+h2y,(37)dp1dx=−12λ1h3.Making use of equations (28) and (37) into equation (36), the first-order velocity distribution takes the form(38)u1x,y=6λ1h3hy−y2−24h−λ0h322y3−3hy2+h2y,(39)px=−92λ12x2+22arctan1/2x2+2x3+20/3xx2+22.The film splits at the point where p1=0; therefore, from the above pressure expression, one can get λ1=0. From equation (39) and the separation point condition for the pressure, we can see clearly that the first-order flow rate is zero; therefore, in order to get more insight into Williamson’s model, we need to take one more term of the perturbation analysis.


### 3.3. Second-Order Problem and Its Solution

Equating the coefficients of termsWe in equation (18), the second-order problem is(40)dp2dx=d2u2dy2+2ddydu0dydu1dy2,(41)Q2=λ2=∫0hxu1dy,with the agreeing first-order boundary conditions:(42)u2=0,aty=0,u2=0,aty=hx.Adopting a similar procedure, using equations (27) and (36) into equation (40) after twice integration and using the boundary conditions (42), we get the following second-order solution:(43)u2x,y=24h3dp0dx33y4−6hy3+4h2y2+12dp2dxy2−hy−12576h3−λ03−1728h2λ0−hλ02h6,(44)dp2dx=−1255h4λ2+96h3−λ03−288h2λ0−hλ02h7.By substituting the first- and second-order pressure gradient distributions into equation (43), one can find the third-order velocity distribution. The third-order pressure distribution at the separation point may be found by integrating the above equation as xb⟶−∞ given by(45)p2x=2.252λ2−0.142−2tan−1x2−π−λ2−0.142x11+34/3λ2−0.142x9+152/3λ2−7.5x7+112λ2−76.74x5+368/3λ2−1.67x3+160/3λ2−49.15x6250000x2+26.From equation (43), a relation between the separation point and third-order flow rate can be found. Setting x=xs and y=h¯xs/2, in order to find the third-order separation point, the regula falsi method has been used with a predefined tolerance of 10−10 which results in λ2=0.315, and the resulting third-order separation point is xs=0.2475. Maple has a Numerical Analysis package that implements many of the techniques we will discuss. The results till the third-order approximation are tabulated in Table 1.Table 1
Effect of Weissenberg number on operating variables.
WeλH/H0xsFPw0.011.32320.66162.41710.02258−5.062430.11.32640.66322.44340.00302−4.706150.21.33580.66792.4874−0.03125−4.127050.31.35160.67582.5486−0.08286−3.152190.41.37360.68682.6263−0.15480−1.541700.51.40200.70102.7191−0.250301.026030.61.43660.71832.8248−0.372134.968890.71.47760.73882.9409−0.5224510.81580.81.52480.76243.0646−0.7024819.22180.91.57840.78923.1804−0.8903529.62930.991.63820.81913.2653−0.8921932.4523Combining the solutions at each order of approximation yield the solutions up to second order for velocity, pressure gradient, and pressure.


## 3.1. Zeroth-Order Problem and Its Solution

Comparing terms free fromWe in equation (18), we have the following zero-order problem:(24)dp0dx=d2udy2,(25)λ0=∫0hxu0dy,with the agreeing zeroth-order boundary conditions:(26)u0=1,aty=0,u0=1,aty=hx.The solution for equation (24) with conditions (26) is given by(27)u0x,y=1+12dp0dxy2−hy.The pressure distribution in equation (27) can be found by utilizing the standard of preservation of mass and has the accompanying structure:(28)dp0dx=12hx−λ0hx3.By using equation (28) into equation (27), one can write the zero-order velocity distribution as(29)u0x,y=1+6h−λ0h3y2−hy.At the separation point located atx=xs and y=h¯xs/2, as shown in Figure 2, the velocity goes to zero. Hence, from equation (29), we can write(30)xs=3λ0−1.Integrating equation (28), write in the following form:(31)p0=12∫xsxhx−λ0hx3dx,where xs is the value of x at p=0. Since we are assuming that the liquid splits evenly to the web and the roll, we can write Q0=2Uhc and in dimensionless form λ0=2hc/H0. The pressure at the separation point may be found by the first integral of equation (28), using assumption xb⟶−∞, as(32)p0x=6−9/2λ0x1+x2/2−3λ0x1+x2/22+122−9λ02tan−1x2+6π21−3λ04.The simplest dynamic model of the separation region is based on the assumption that the film splits at the point, whereu0=0andp0=0; resultantly, the above equation becomes a transcendental equation in λ0 whose solution may be obtained by a numerical technique named as the modified regula falsi method with a predefined accuracy of 10−10. For the said accuracy, the zero-order flow rate is λ0=1.3232524, which is correct up to five decimal places, and the residual error in this case is 0.0000011085; this coating thickness has been achieved at the separation point xs=1.723298349. It is worth to mention that the above zero-order results are more accurate as calculated by Middleman [17].


## 3.2. First-Order Problem and Its Solution

Equating the coefficients of termsWe in equation (18), the first-order problem is(33)dp1dx=d2u1dy2+ddydu0dy2,(34)Q1=λ1=∫0hxu1dy,with the agreeing first-order boundary conditions:(35)u1=0,aty=0,u1=0,aty=hx.Adopting the same procedure as in Section3.1, applying equation (27) into equation (33), after twice integration and using conditions (35), we acquire the following first-order solution:(36)u1x,y=12dp1dxy2−hy−16dp0dx22y3−3hy2+h2y,(37)dp1dx=−12λ1h3.Making use of equations (28) and (37) into equation (36), the first-order velocity distribution takes the form(38)u1x,y=6λ1h3hy−y2−24h−λ0h322y3−3hy2+h2y,(39)px=−92λ12x2+22arctan1/2x2+2x3+20/3xx2+22.The film splits at the point where p1=0; therefore, from the above pressure expression, one can get λ1=0. From equation (39) and the separation point condition for the pressure, we can see clearly that the first-order flow rate is zero; therefore, in order to get more insight into Williamson’s model, we need to take one more term of the perturbation analysis.


## 3.3. Second-Order Problem and Its Solution

Equating the coefficients of termsWe in equation (18), the second-order problem is(40)dp2dx=d2u2dy2+2ddydu0dydu1dy2,(41)Q2=λ2=∫0hxu1dy,with the agreeing first-order boundary conditions:(42)u2=0,aty=0,u2=0,aty=hx.Adopting a similar procedure, using equations (27) and (36) into equation (40) after twice integration and using the boundary conditions (42), we get the following second-order solution:(43)u2x,y=24h3dp0dx33y4−6hy3+4h2y2+12dp2dxy2−hy−12576h3−λ03−1728h2λ0−hλ02h6,(44)dp2dx=−1255h4λ2+96h3−λ03−288h2λ0−hλ02h7.By substituting the first- and second-order pressure gradient distributions into equation (43), one can find the third-order velocity distribution. The third-order pressure distribution at the separation point may be found by integrating the above equation as xb⟶−∞ given by(45)p2x=2.252λ2−0.142−2tan−1x2−π−λ2−0.142x11+34/3λ2−0.142x9+152/3λ2−7.5x7+112λ2−76.74x5+368/3λ2−1.67x3+160/3λ2−49.15x6250000x2+26.From equation (43), a relation between the separation point and third-order flow rate can be found. Setting x=xs and y=h¯xs/2, in order to find the third-order separation point, the regula falsi method has been used with a predefined tolerance of 10−10 which results in λ2=0.315, and the resulting third-order separation point is xs=0.2475. Maple has a Numerical Analysis package that implements many of the techniques we will discuss. The results till the third-order approximation are tabulated in Table 1.Table 1
Effect of Weissenberg number on operating variables.
WeλH/H0xsFPw0.011.32320.66162.41710.02258−5.062430.11.32640.66322.44340.00302−4.706150.21.33580.66792.4874−0.03125−4.127050.31.35160.67582.5486−0.08286−3.152190.41.37360.68682.6263−0.15480−1.541700.51.40200.70102.7191−0.250301.026030.61.43660.71832.8248−0.372134.968890.71.47760.73882.9409−0.5224510.81580.81.52480.76243.0646−0.7024819.22180.91.57840.78923.1804−0.8903529.62930.991.63820.81913.2653−0.8921932.4523Combining the solutions at each order of approximation yield the solutions up to second order for velocity, pressure gradient, and pressure.


## 4. Operating Variables

Once the velocity, pressure distribution, and pressure gradient are found, then all other interesting engineering quantities are readily available. The operating variables are computed in the following manner.


### 4.1. Coating Thickness

The substrate thicknessHb at the entering is given by(46)Hb=1+xb22.It is a very large number, in the case of an infinite reservoir.Since we are assuming that fluid separates uniformly to coat both the roll and the plane, we can write volumetric flow rate in the dimensional form asQ=2UH and(47)HH0=λ2.This gives the resulting leaving coating thickness in a dimensionless form.The examination is yet incomplete at this point as the relationship between the nondimensional flow rateλ and the separation point xs needs to be evaluated. Figure 2 indicates that the web separates equally, and, consequently, the segregation point is xs,hxs/2; here, the pressure and velocity profile approach zero. This is assumed as the roll and sheet has equal velocity. It causes segregation point xs depending upon the dimensionless coating thickness λ.In case of Newtonian fluid, the relationship between segregation pointxs, nondimensional flow rate, and the rheological parameters admits the value of one of these parameters depending on others. As in case of Williamson fluid, it is difficult to calculate explicitly xs in terms of λ. To find the segregation point xs, some numerical system is needed.


### 4.2. Roll Separating Force

The roll separating forceF is given by(48)F=∫−∞xspxdx,where F=F¯H0/η0URW,F¯ is the dimensional roll separating force per unit width W.


### 4.3. Power Contribution

The power transmitted to the fluid by the roll is calculated by integrating the product of shear stress and the roll surface speed over the roll surface which is obtained by settingy¯=H0, as(49)Pw=∫−∞xsτxyx,1dx.Here,Pw=Pw¯/η0WU2 is the nondimensional power and τxy=τ¯xyH0/η0U the nondimensional stress tensor given by equation (14).


### 4.4. Adiabatic Temperature

The power contribution has the capability to increase the temperature of fluid to the maximum, given by rise in adiabatic temperatureΔTave:(50)ΔTave=PwQρCp,where Cp is the melt heat capacity at constant pressure.


## 4.1. Coating Thickness

The substrate thicknessHb at the entering is given by(46)Hb=1+xb22.It is a very large number, in the case of an infinite reservoir.Since we are assuming that fluid separates uniformly to coat both the roll and the plane, we can write volumetric flow rate in the dimensional form asQ=2UH and(47)HH0=λ2.This gives the resulting leaving coating thickness in a dimensionless form.The examination is yet incomplete at this point as the relationship between the nondimensional flow rateλ and the separation point xs needs to be evaluated. Figure 2 indicates that the web separates equally, and, consequently, the segregation point is xs,hxs/2; here, the pressure and velocity profile approach zero. This is assumed as the roll and sheet has equal velocity. It causes segregation point xs depending upon the dimensionless coating thickness λ.In case of Newtonian fluid, the relationship between segregation pointxs, nondimensional flow rate, and the rheological parameters admits the value of one of these parameters depending on others. As in case of Williamson fluid, it is difficult to calculate explicitly xs in terms of λ. To find the segregation point xs, some numerical system is needed.


## 4.2. Roll Separating Force

The roll separating forceF is given by(48)F=∫−∞xspxdx,where F=F¯H0/η0URW,F¯ is the dimensional roll separating force per unit width W.


## 4.3. Power Contribution

The power transmitted to the fluid by the roll is calculated by integrating the product of shear stress and the roll surface speed over the roll surface which is obtained by settingy¯=H0, as(49)Pw=∫−∞xsτxyx,1dx.Here,Pw=Pw¯/η0WU2 is the nondimensional power and τxy=τ¯xyH0/η0U the nondimensional stress tensor given by equation (14).


## 4.4. Adiabatic Temperature

The power contribution has the capability to increase the temperature of fluid to the maximum, given by rise in adiabatic temperatureΔTave:(50)ΔTave=PwQρCp,where Cp is the melt heat capacity at constant pressure.


## 5. Results and Discussions

This article investigates the roll coating procedure for an incompressible Williamson fluid. For the simplification of motion equations, the LAT is applied. Maple 15 software has been used to generate numerical results for the volumetric flow rateλ, the separation point xs, the exit sheet thickness H/H0, the power contribution, and the roll separating force presented in Table 1 for various values of the Weissenberg number We. The maximum coating thickness has been observed as 0.8200 by approaching We⟶1, and for this coating thickness, the optimum point of separation is detected. It is noticed that an increase in We causes an increase in the coating thickness. The minimum coating thickness when We⟶0 has been observed as 0.6610. The coating thickness, volume flow rate, and separation point are increasing functions of We (forces of inertia are controlling viscous force).The outcome for the dimensionless velocity profiles is projected in Figures3–10. Through the variation of We, velocity profiles at various locations of the roll coating procedure are presented. Figures 3 and 4 are plotted at positions x = 0 and x = 0.25, respectively, indicating that the increase in We first causes the velocity of the fluid to decrease in 0 < y < 0.602 and then increase monotonically beyond 0.6.Figure 3
Effect of the Weissenberg number on velocity distribution at nip region (x = 0).Figure 4
Effect of the Weissenberg number on velocity distribution atx = 0.25.Figure 5
Effect of the Weissenberg number on velocity distribution atx = 0.5.Figure 6
Effect of the Weissenberg number on velocity distribution atx = 0.75.Figure 7
Effect of the Weissenberg number on velocity distribution atx = 1.Figure 8
Effect of the Weissenberg number on velocity distribution atx = 1.5.Figure 9
Effect of the Weissenberg number on velocity distribution atx = 2.Figure 10
Effect of the Weissenberg number on velocity distribution atx = 2.5.It is interesting to see from Figure3 that, at the nip region, the reverse flow has been observed for the higher values of We, and as the material moves away from the nip in the direction of separation point, the backward flow decreases as shown graphically in Figures 4–6 and the fluid becomes more developed. Figure 4 indicates that all the velocity curves have the same point of intersection, that is, (y = 0.6, u = 1.4263) for different values of the Weissenberg number.Figures5 and 6 have been depicted at positions x = 0.5 and 0.75, respectively; here, one can observe that an increase in We causes increase in the velocity profile as We is the ratio of viscous forces to the elastic forces. It indicates that the viscous forces become controlling forces once the web reaches the position x > 0.25. It is also observed from these graphs that fluid becomes fully developed when one moves toward the separation point.Comparative to Figures5 and 6, a symmetric but opposite behavior has been observed in Figures 7–10, which are sketched at x = 1, x = 1.5, x = 2, and x = 2.5, respectively. With an increase in the Weissenberg number, the velocity distribution increases. It is also highlighted that, as the fluid approaches the separation point, the fluid viscosity increases and dominance of viscous force over elastic force will cause coating on the web; this behavior can clearly been seen from Figures 7–10, Furthermore, from these figures, one can analyze that the gap between the velocity lines near to the web is smaller as compared to the gab between the velocity line near to the roll.The result for the dimensionless pressure gradient distributions is presented in Figure11, whereas the outcomes for the dimensionless pressure profile are projected in Figure 12. It is observed that the pressure gradient rises by increasing the Weissenberg number, whereas the pressure distribution decreases by increasing the Weissenberg number. Starting from zero pressure at the attachment point, pressure starts increasing and attains its extreme values right at x=±0.86.Figure 11
Effect of the Weissenberg number on pressure gradient distribution.Figure 12
Effect of the Weissenberg number on pressure distribution.


## 6. Conclusion

Using lubrication approximation theory, the system of nonlinear partial differential equations is reduced to ordinary differential equations. Then, the obtained nonlinear ODEs are linearized using the regular perturbation technique. The influence of various evolving physical parameters on velocity components, flow rate, pressure profile, roll separating force, power contribution, separation point, and most important coating thickness have been studied theoretically. Some of the results are tabulated numerically and some are presented graphically.The key findings of the present study are as follows:(i)
For the pseudoplastic fluids, the results by Middleman [17] were verified when We⟶0 and extended(ii)
Weissenberg’s number acts as a controlling parameter for the rate of flow, thickness in coating, power contribution, pressure, roll separating force, and separation point(iii)
The extreme values for the pressure distribution have been found atx=±0.86(iv)
Viscous force plays a vital and dominating role in coating thickness, force of separation, power, and pressure profile(v)
The reverse flow has been observed at the nip region for the higher values ofWe, and as the material moves away from the nip in the direction of separation point, the backward flow decreases(vi)
It is worth mentioning that the above zero-order results are more accurate as calculated by Middleman



---

*Source: 3431885-2021-08-23.xml*