Research Article  Open Access
On the Study of Viscoelastic Walters' B Fluid in Boundary Layer Flows
Abstract
Viscoelastic Walters' B fluid flows for three problems, stagnationpoint flow, Blasius flow, and Sakiadis flow, have been investigated. In each problem, Cauchy equations are changed to a nondimensional differential equations using stream functions and with assumption of boundary layer flow. The fourthorder predictorcorrector finitedifference method for solving these nonlinear differential equations has been employed. The results that have been obtained using this method are compared with the results of the last studies, and it is clarified that this method is more accurate. It is also shown that the results of last study about Sakiadis flow of Walter's B fluid are not true. In addition, the effects of order of discretization in the boundaries are investigated. Moreover, it has been discussed about the valid region of Weissenberg numbers for the secondorder approximation of viscoelastic fluids in each case of study.
1. Introduction
Boundary layer flows of nonNewtonian fluids can be a beneficial approach for modeling several manufacturing processes of industry such as the aerodynamic extrusion of plastic sheets, the cooling of metallic plates, fabrication of adhesive tapes, application of coating and layers onto rigid substrates. Some of articles that have investigated these kinds of fluid flows are given in [1–24]. These types of phenomenon can be considered in two branches, one about the geometry of flow and it’s boundary conditions and the other about the viscoelastic model. If we take a glance at the literature, we find that there are many models for considering fluids that have both elastic and viscous properties. One of these models is Walters’ B that has been used in several studies and for modeling many phenomena [25–33]. Study of the flow problems of this class of fluids not only is important technologically, but is also challenging to engineers and applied mathematicians who are interested in obtaining accurate solutions.
Constitutive equations of viscoelastic fluids usually generate higherorder derivative terms in the momentum equations that make them more difficult to solve in comparison with Newtonian fluids. There are some analytical and numerical methods for investigating these kinds of fluids such as homotopy analysis method [29] and spectral method [34]. Because of the apparent nonavailability of extra boundary conditions, researchers tend to develop a regular perturbation solution of the problem, taking the solution for the Newtonian fluid as the primary solution and the firstorder perturbed solution as the secondary solution. A classical example of this technique is the analysis of the twodimensional flow of Walters’ B fluid about a stagnation point given by Beard and Walters [11]. Perturbation approach can lead to erroneous results not only quantitatively but also qualitatively which is shown by some researchers such as Ariel [35]. He proposed a predictorcorrector method for solving linear and nonlinear differential equations such as boundary layer stagnationpoint flow. This method can solve a differential equation by an explicit method with fourth order of approximation and without any restriction of convergence.
In this paper, similarity solutions of three viscoelastic Walters’ B boundary layer flow problems, stagnationpoint flow, Blasius flow and Sakiadis flow, are obtained and all of the equations are changed to nondimensional forms in the second section. The mentioned fourthorder predictorcorrector method for solving these nonlinear differential equations is presented in the third section, and the results are compared with last studies in the fourth section. It is obtained that this method is more accurate and straightforward. It is also shown that the results of last article about Sakiadis flow by Sadeghy and Sharifi [36] are not true. Moreover, valid regions of Weissenberg numbers for the secondorder model in each problem are obtained, and the effects of order of discretization in the boundaries are investigated.
2. Mathematical Formulation
In this section, nondimensional equations of motions for three different boundary layer flows, stagnationpoint flow, Blasius flow, and Sakiadis flow, are obtained. These equations are solved in Section 4 using a predictorcorrector method.
2.1. StagnationPoint Flow
Cauchy equations are employed for obtaining boundary layer equation of twodimensional stagnationpoint flow. Steady Cauchy equations are where is pressure, is density of fluid, and and are velocity components in and directions, respectively. Also, , , , and are components of stress matrix. These variables can be obtained for WaltersB viscoelastic fluid by following equation: which is the constitutive equation of Walters’ B fluid and is the dynamic viscosity, is the shortmemory coefficient, and is the deformation rate tensor. Also is upperconvected derivative of a tensor and can be obtained from this equation: So, the stress components in steady conditions are where can be obtained as follows: On the basis of the above relation, the stress components can be written as follows:
These equations are the relation of components of stress tensor and velocity components. Because of requirement for derivatives of these stress components, these derivatives are presented.
By substituting these relations in (2.1), Cauchy equations become
momentum:momentum: The and momentum equations can be simplified as follows after employing boundary layer approximation:
On the basis of viscous flow theory, the change in pressure across boundary layer is . Therefore, the pressure gradient term in (2.9) can be obtained from the flow just outside the boundary layer. So, following equation is obtained in this region and under the steady flow condition: where is the mainstream velocity. So, boundary layer equation for this fluid flow is obtained as follows: Because of the similarity of the velocity profiles in different crosssections, it is concluded that the only possible form of under these conditions for obtaining the solution of (2.11) is where is a constant parameter [37]. This procedure has been described in [37] completely.
The momentum equation can be nondimensionalized using stream function . So, the velocity components can be written as The stream function can be made dimensionless as follows: where is the dimensionless stream function, is the kinematic viscosity, and is the similarity variable, defined by By substituting (2.11)–(2.15) into (2.9), dimensionless equation can be obtained: where is the local Weissenberg number defined by On the basis of no slip condition, no permeability, and , three boundary conditions of (2.16) are
Order of dimensionless equation of motion, (2.16), is the fourth while the momentum equations are the third order. This is because of derivative operations that are applied to momentum equations for obtaining dimensionless equation (2.16).
In addition to velocities and stream functions that are common parameters for investigating boundary layer problems, shear stress as another important parameter should obtained in the wall. Therefore, the following equation for shear is presented: By substituting and on the basis of (2.13) and (2.14) on (2.19), shear stress on the solid boundary can be obtained using (2.18) as follows:
2.2. Blasius Flow
Equation of motion for Blasius flow can be obtained using (2.1) that leads to (2.9).
The momentum equation of (2.9) can be nondimensionalized using stream function . So, the velocity components can be written as The stream function can be made dimensionless as follows: where is the dimensionless stream function, is the kinematic viscosity, and is the similarity variable, defined by So, nondimensional equation of motion for Blasius flow can be obtained on the basis of as follows: where is the local Weissenberg number defined by Boundary conditions of (2.24) are By substituting and from (2.21) and (2.22) on (2.19) and employing boundary conditions that are presented in (2.26), shear stress on the solid boundary can be obtained as follows:
2.3. Sakiadis Flow
Sakiadis flow for viscoelastic Walters’ B fluid has been investigated by Sadeghy and Sharifi [36] at 2004. The results that are presented in [36] are not true and reliable which are based on false nondimensional equation of motion that is derived in [36]. In this paper, the true equation of motion is presented and reliable results are obtained.
Sakiadis flow is similar to Blasius flow and the only difference of these fluid flows is their boundary conditions that are presented in following relations:
In the next section, the predictorcorrector method for analysis of nondimensional equations of these flows is described.
In addition, shear stress on the solid boundary can be obtained similar to (2.27) and using boundary conditions that are presented in (2.28) as follows:
3. Method of Solution
In this section, fourthorder predictorcorrector method that was proposed by Ariel [35] has been described. It should be noted that stagnationpoint flow has been investigated by Ariel [35] using this method. So, solution of this flow is just obtained to show the true using of this approach.
3.1. TwoDimensional StagnationPoint Flow
It should be noted that because of the singularity of (2.16) and (2.19) at , some usual methods, such as shooting method, cannot be employed for solving these equations. Henceforth, the predictorcorrector method, which not only is an explicit method without requiring for solving any matrix but also converges without any convergence condition, is used in this study. For clarity, this method is described by applying to (2.16) that is the nondimensional equation of motion for stagnationpoint flow.
The auxiliary parameters should be used to rewritte (2.16) as follows: and boundary conditions (2.18) become
The key of the algorithm is to retrain the secondorder derivative in (3.2) and not to further decompose it into a pair of differential equations in and .
The following equations should be used for approximating the first and secondorder derivatives: where is the derivative operator and is the backward operator.
These operators are employed for discretizing (3.1) and (3.2) which lead to the following equations: Equation (3.5) can be explicitly solved for at the ()th point in a supposing mesh generation. Using the subscript “” for the predicted values, we can write These steps were about the predicting the values of . So, the following operators are necessary to obtain the correct values of these variables: On the basis of (3.7), (3.1) and (3.2) become
Therefore, we have a fourthorder predictorcorrector method for obtaining the values of and at a particular point based with respect to their values at the previous five mesh points. The values of and are given at . Since is given at infinity, (2.16) can be thought as a boundary value problem (BVP). Therefore, we can solve these equations by changing the initial guess of so that the boundary condition at infinity is satisfied.
Assume that the amount of the initial guess of is
The amounts of in at five extra nodes are needed to start the solution of the discretized equation. So, by applying the Taylor series around zero, we can find the amount of these parameters in . Taylor series for expansion of of around results in the following equation: where are constant and should be obtained using the boundary values of at . Henceforth, some of these constants are where are obtained on the basis of boundary values of parameters at and others are obtained using (2.16).
3.2. Blasius Flow
Nondimensional equation of Blasius flow (2.19) can be solved by the proposed fourthorder predictorcorrector method. So, the predictor and corrector relations are obtained as follows.
Predictor: Corrector: It is obvious that these equations can be solved explicitly.
Quantities of extra nodes that are necessary for solving above equations are obtained similar to those obtained for stagnationpoint flow using Taylor series. So, we have where where are obtained by employing the boundary values of parameters at and others are obtained using (2.19).
3.3. Sakiadis Flow
It is mentioned that the only difference between Sakiadis flow and Blasius flow is just in their boundary conditions. So, their predictor and corrector relations are same and we can use those of Blasius flow for solving (2.19) which is the nondimensional equation of motion of Sakiadis flow. On the other hand, the constants , and so forth, are different from Blasius flow that using (2.19) and (2.21) are obtained as In the next section, the predictorcorrector method that is described in this section is applied to the mentioned problems, stagnationpoint flow, Blasius flow, and Sakiadis flow. Therefore, the results are obtained and compared with the results of the other papers that are obtained by different methods. In addition, the results presented by Sadeghy and Sharifi [36] in 2004 are reformed.
4. Result and Discussion
4.1. StagnationPoint Flow
The results of the solution of stagnationpoint flow with the proposed fourthorder method are presented in Table 1.

On the basis of Table 1, it is obvious that step size () must be smaller than 0.05 in this problem to obtain the reliable values of which is the key point of the solution of this problem. In addition, larger values cause some errors in the results for smaller value of Weissenberg number which can be omitted using smaller values of which are chosen as 10 in these situations. Because of these errors, far from the wall, the numerical error becomes larger and makes an unstable numerical procedure. Thus, for we used , and for greater values of this viscoelastic parameter, has been set equal to 20 to obtain the suitable results.
Moreover, the comparison of values of shear stress at solid boundary can be beneficial in this problem. On the basis of (2.20), this parameter in a particular location of the domain depends linearly on . So, it is obvious that, for larger values of , shear stress on the solid boundary is larger than for smaller ones. Therefore, it can be concluded from Table 1 that this parameter increases as Weissenberg number increases. It should be noted that valid values of are based on step sizes smaller than 0.05.
The boundary layers in stagnationpoint flow for different Weissenberg numbers are compared in Figure 1. As Weissenberg number increases, the boundary layer thickness decreases. In addition, Figure 1 clarifies the differences between the quantities of nondimensional horizontal velocity at any particular nondimensional parameter . On the basis of this figure, values of nondimensional horizontal velocity increases as Weissenberg number increases.
As we see in the literature [38, 39], this ODE equation has been solved with other numerical scheme, perturbation method, kellerbox method, and so forth, nonphysical phenomena, which are overshoot in the boundary layer, and occurred for . It is clarified that Walters’ B model and its constitutive equation are not suitable for this range of Weissenberg number. Figure 2 shows stream lines for different values of Weissenberg number. The differences of boundary layer thicknesses for different values of Weissenberg numbers are clarified in this figure.
4.2. Blasius Flow
Quantities of for different Weissenberg numbers are obtained, for , as shown in Table 2. It is shown that increasing the amount of causes to reduce the amount of . Therefore, on the basis of (2.27), shear stress at the solid boundary is smaller for larger values of Weissenberg number. The amount of for all values of is opted 20, that is, similar to the solution of stagnationpoint flow.

Figure 3 shows the growth of boundary layers in Blasius flow for different Weissenberg numbers. As increases, the boundary layer thickness decreases. It should be noted that a nonphysical phenomenon occurred for larger values of Weissenberg number, , which clarifies that this model is useless for investigating this problem in this this range of Weissenberg number.
The differences between the values of nondimensional horizontal velocity at any particular nondimensional parameter is shown in Figure 3. So, it is obvious that this variable decreases as Weissenberg number increases.
4.3. Sakiadis Flow
Variation of in twodimensional Sakiadis flow is clarified in Table 3. So, it is obvious that values of become larger as Weissenberg number increases. Therefore, similar to Blasius flow and on the basis of (2.29), shear stress at the solid boundary is smaller for larger values of Weissenberg number with respect to smaller values.

The most important result of this study is shown in Figure 4. This figure that presents the growth of boundary layer of Sakiadis flow clarifies that results of Sadeghy and Sharifi [36] are unsatisfactory and clearly false. It is obvious that, for larger values of Weissenberg number, more than 0.3, Walters’ B model is not a reliable model for investigating the viscoelastic Sakiadis flow. But, Sadeghy and Sharifi [36] had concluded that this model is satisfactory for larger values even more than 1.
4.4. Order of Approximation of Boundary Conditions
One of the fundamental principles of numerical methods to obtain an arbitrary order of approximation is the sameness of order of discretization both in the entire domain and in the boundaries. The predictorcorrector method that is used in this study is a the fourth order method. So, boundary condition also should be approximated in the fourthorder. Amounts of for 2nd and 4th order of approximation for calculating boundary conditions in stagnation point, Blasius, and Sakiadis flows are clarified in Tables 4, 5, and 6. On the basis of these tables, the effect of the approximation of boundary conditions is clarified.



It is obvious that the amounts of for all three problems are more accurate by the 4th order of approximation in boundary conditions.
5. Conclusion
In this study, the fourthorder predictorcorrector method is employed for solving the three viscoelastic fluid flow problems. Equations of motion for these problems, stagnationpoint flow, Blasius flow, and Sakiadis flow, are obtained using Cauchy equation of motion, with assumption of boundary layer flow and turned into nondimensional forms. The discretized equations on the approach of the predictorcorrector method are solved, and the results are compared with the results of last studies, that is, clarified the high accuracy of the method. Also it has been shown that as elasticity increases, the stress on the solid boundary increases in stagnationflow problem and decreases in Blasius and Sakiadis flows. In addition, it is shown that the results of the last study about Sakiadis flow by Sadeghy and Sharifi [36] are not true. Moreover, it has been shown that the effect of order of discretizations in the boundaries is an important issue. Finally, it is shown that there are restrictions for valid region of Weissenberg numbers and the secondorder method of viscoelastic fluids failed out of these ranges.
References
 K. L. Hsiao, “Mixed convection with radiation effect over a nonlinearly stretching sheet,” Proceedings of World Academy of Science, Engineering and Technology, vol. 62, pp. 338–342, 2010. View at: Google Scholar
 H. I. Andersson and J. B. Aarseth, “Sakiadis flow with variable fluid properties revisited,” International Journal of Engineering Science, vol. 45, no. 2–8, pp. 554–561, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. M. Patil, S. Roy, and A. J. Chamkha, “Double diffusive mixed convection flow over a moving vertical plate in the presence of internal heat generation and a chemical reaction,” Turkish Journal of Engineering and Environmental Sciences, vol. 33, no. 3, pp. 193–205, 2009. View at: Publisher Site  Google Scholar
 M. A. A. Mahmoud, “Variable viscosity effects on hydromagnetic boundary layer flow along a continuously moving vertical plate in the presence of radiation,” Applied Mathematical Sciences. Journal for Theory and Applications, vol. 1, no. 17–20, pp. 799–814, 2007. View at: Google Scholar  Zentralblatt MATH
 A. Pantokratoras, “Asymptotic suction profiles for the Blasius and Sakiadis flow with constant and variable fluid properties,” Archive of Applied Mechanics, vol. 79, no. 5, pp. 469–478, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Pantokratoras, “The Blasius and Sakiadis flow with variable fluid properties,” Heat and Mass Transfer, vol. 44, no. 10, pp. 1187–1198, 2008. View at: Publisher Site  Google Scholar
 R. C. Bataller, “Numerical comparisons of blasius and sakiadis flows,” Mathematica, vol. 26, no. 2, pp. 187–196, 2010. View at: Google Scholar
 T. Fang and J. Zhang, “Closedform exact solutions of MHD viscous flow over a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2853–2857, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Y. Ghaly and E. M. E. Elbarbary, “Radiation effect on MHD freeconvection flow of a gas at a stretching surface with a uniform free stream,” Journal of Applied Mathematics, vol. 2, no. 2, pp. 93–103, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Cortell, “A numerical tackling on Sakiadis flow with thermal radiation,” Chinese Physics Letters, vol. 25, no. 4, pp. 1340–1342, 2008. View at: Publisher Site  Google Scholar
 D. W. Beard and K. Walters, “Elasticoviscous boundarylayer flows. I. Twodimensional flow near a stagnation point,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 60, pp. 667–674, 1964. View at: Google Scholar  Zentralblatt MATH
 A. Acrivos, M. J. Shah, and E. E. Peterson, “Momentum and heat transfer in laminar boundary layer flows of nonNewtonian fluids past external surfaces,” AIChE Journal, vol. 6, p. 312, 1960. View at: Publisher Site  Google Scholar
 A. Acrivos, M. J. Shah, and E. E. Petersen, “On the solution of the twodimensional boundarylayer flow equations for a nonNewtonian power law fluid,” Chemical Engineering Science, vol. 20, no. 2, pp. 101–105, 1965. View at: Publisher Site  Google Scholar
 G. D. Bixzell and J. C. Alattery, “NonNewtonian boundary layer flow,” Chemical Engineering Science, vol. 17, p. 777, 1962. View at: Publisher Site  Google Scholar
 N. Hayasi, “Similarity of twodimensional and axisymmetric boundary layer flows of nonNewtonian fluids,” Journal of Fluid Mechanics, vol. 23, pp. 293–303, 1965. View at: Publisher Site  Google Scholar
 J. N. Kapur and R. C. Srivastava, “Similar solutions of the boundary layer equations for power law fluids,” ZAMP, vol. 14, no. 4, pp. 383–389, 1963. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Y. Lee and W. F. Ames, “Similar solutions for nonNewtonian fluids,” AIChE Journal, vol. 12, p. 700, 1966. View at: Publisher Site  Google Scholar
 W. R. Schowalter, “The application of boundary layer theory to power law pseudoplastic fluids: similar solutions,” AIChE Journal, vol. 6, no. 1, pp. 24–28, 1960. View at: Publisher Site  Google Scholar
 A. G. Hansen and T. Y. Na, “Similarity solutions of laminar, incompressible boundary layer equations of nonNewtonian fluids,” in Proceedings of the ASME Winter Annual Meeting, November 1967, ASME 67WA/FEZ. View at: Google Scholar
 M. M. Denn, “Boundary layer flows for a class of elastic fluids,” Chemical Engineering Science, vol. 22, no. 3, pp. 395–405, 1967. View at: Publisher Site  Google Scholar
 R. W. Seth, “Solution of a viscoelastic boundary layer equaiions by otthog&& collocation,” Journal of Engineering Mathematics, vol. 8, no. 2, pp. 89–92, 1974. View at: Publisher Site  Google Scholar
 K. R. Rajagopal, T. Y. Na, and A. S. Gupta, “Flow of a viscoelastic fluid over a stretching sheet,” Rheologica Acta, vol. 23, no. 2, pp. 213–215, 1984. View at: Publisher Site  Google Scholar
 K. R. Rajagopal, A. S. Gupta, and T. Y. Na, “A note on the falknerskan flows of a nonNewtonian fluid,” International Journal of NonLinear Mechanics, vol. 18, no. 4, pp. 313–320, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 W. H. H. Banks and M. B. Zaturska, “Eigensolutions in boundarylayer flow adjacent to a stretching wall,” IMA Journal of Applied Mathematics, vol. 36, no. 3, pp. 263–273, 1986. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Postelnicu, “Thermal hydrodynamic instability of a Walters B viscoelastic fluid in a fluidsaturated an isotropic porous medium with fast chemical reaction,” in Proceedings of the 81st EUROTHERM Seminar on Reactive Heat Transfer in Porous Media, Ecole des Mines d'Albi, France, June 2007. View at: Google Scholar
 P. Kumar, G. J. Singh, and R. Lal, “Thermal instability of Walters B' viscoelastic fluid permeated with suspended particles in hydro magnetic in porous medium,” Thermal Science, vol. 8, no. 1, pp. 51–61, 2004. View at: Google Scholar
 V. Sharma and U. Gupta, “Stability of stratified elastoviscous Walters' (Model B') fluid in the presence of horizontal magnetic field and rotation,” Studia Geotechnica et Mechanica, vol. 32, no. 2, pp. 41–54, 2010. View at: Google Scholar
 P. Kumar and G. J. Singh, “On the stability of two stratified Walters B' viscoelastic superposed fluids,” Studia Geotechnica et Mechanica, vol. 32, no. 4, pp. 29–38, 2010. View at: Google Scholar
 A. A. Joneidi, G. Domairry, and M. Babaelahi, “Homotopy analysis method to Walter's B fluid in a vertical channel with porous wall,” Meccanica, vol. 45, no. 6, pp. 857–868, 2010. View at: Publisher Site  Google Scholar
 S. Nadeem and N. S. Akbar, “Peristaltic flow of Walter's B fluid in a uniform inclined tube,” Journal of Biorheology, vol. 24, no. 1, pp. 22–28, 2010. View at: Publisher Site  Google Scholar
 R. K. Thakur, “Comparison between Newtonian and nonNewtonian fluid,” International Journal of Educational Administration, vol. 2, no. 4, pp. 641–645, 2010. View at: Google Scholar
 S. Bariş, “Steady threedimensional flow of a Walter's B' fluid in a vertical channel,” Turkish Journal of Engineering and Environmental Sciences, vol. 26, no. 5, pp. 385–394, 2002. View at: Google Scholar
 S. G. Mohiddin, V. R. Prasad, S. V. K. Varma, and O. A. Bég, “Numerical study of unsteady free convective heta and mass transfer in a waltersB viscoelastic flow along a vertical cone,” International Journal of Applied Mathematics and Mechanics, vol. 6, no. 15, pp. 88–114, 2010. View at: Google Scholar
 N. Khabazi and K. Sadeghy, “Hydromagnetic instability of viscoelastic fluids in Blasius flow,” Nihon Reoroji Gakkaishi, vol. 37, no. 4, pp. 173–180, 2009. View at: Publisher Site  Google Scholar
 P. D. Ariel, “Generalized Gear's method for computing the flow of a viscoelastic fluid,” Computer Methods in Applied Mechanics and Engineering, vol. 142, no. 12, pp. 111–121, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Sadeghy and M. Sharifi, “Local similarity solution for the flow of a “secondgrade” viscoelastic fluid above a moving plate,” International Journal of NonLinear Mechanics, vol. 39, no. 8, pp. 1265–1273, 2004. View at: Publisher Site  Google Scholar
 H. Schlichting and K. Gersten, Boundary Layer Theory, Springer, 8th edition, 2000.
 R. W. Serth, “Solution of a viscoelastic boundary layer equation by orthogonal collocation,” Journal of Engineering Mathematics, vol. 8, no. 2, pp. 89–92, 1974. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. D. Ariel, “Hybrid method for computing the flow of viscoelastic fluids,” International Journal for Numerical Methods in Fluids, vol. 14, no. 7, pp. 757–774, 1992. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2012 Seyed Ali Madani Tonekaboni 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.