Magnetohydrodynamics Free Convection Flow of Incompressible Fluids over Corrugated Vibrating Bottom Surface with Hall Currents and Heat and Mass Transfers

Bulinda, Vincent M.; Kang’ethe, Giterere P.; Kiogora, Phineas R.

doi:https://doi.org/10.1155/2020/2589760

Journal of Applied Mathematics

On this page

Abstract Introduction Conclusion Data Availability Disclosure Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 2589760 | https://doi.org/10.1155/2020/2589760

Magnetohydrodynamics Free Convection Flow of Incompressible Fluids over Corrugated Vibrating Bottom Surface with Hall Currents and Heat and Mass Transfers

Vincent M. Bulinda,¹Giterere P. Kang’ethe,¹and Phineas R. Kiogora¹

Academic Editor: Sazzad Hossien Chowdhury

Received15 Nov 2019

Accepted13 Jan 2020

Published01 Mar 2020

Abstract

Magnetohydrodynamics free convection flow of incompressible fluids over corrugated vibrating bottom surface with Hall currents and heat and mass transfers considering heat flux is discussed. The corrugation patterns suggested are sinusoidal in nature. The governing equations are solved by the explicit finite difference numerical method of the forward-time backward-space scheme to obtain the analytical results for velocity, concentration, and temperature profiles. The unsteady resultant velocities, concentration, and temperature for various values of physical parameters are discussed in detail, and it is shown that they have significant effects on the fluid flow, and heat and mass transfers are shown graphically.

1. Introduction

More advances in modern technology calls for the need to design processing equipment that complies with the economic, environmental, and energy saving with high efficiency of mass and heat transfers.

Corrugated structures are well known for heat transfer enhancement applied in various engineering problems such as heat exchangers, air-conditioning devices, refrigeration systems, chemical reactors, and fuel rod [1]. Corrugated channels can be designed with different corrugation angles with different arrangements of channel surface [2, 3]. Corrugated ducts are designed in a way that fabricates heat [4] and improves heat transfer efficiency in an appropriate Reynolds number regime [5]. A numerical study on heat transfer of sinusoidal wavy channels at zero-degree phase shift was conducted by [6, 7]. Abo-Eldahab and Aziz [8] studied hydrodynamics and thermal behaviour of free connective flow over an isothermal, vertical flat plate and convective flow of heat-generating fluid past a vertical plate with Hall current and heat flux through a porous medium. Idowu et al. [9] studied the impact of heat and mass transfers on MHD oscillatory flow of Jeffery fluid in a porous channel with thermal conductivity, Dufour and Soret. 3D numerical predictions of fluid flow and heat transfer between cross-corrugated plates were studied in [10] and obtained for the same geometry. They found that all models used predicted practically satisfactory and acceptable within the current Reynolds number range.

MHD flows with Hall and ion-slip currents have applications in the design of MHD generators, Hall accelerators, and flight magnetohydrodynamics [5, 11, 12]. The Hall effect is important since it can be applied in determining if a given material is a semiconductor or insulator, and it can measure the magnetic field and also in sensing as they are immune to water, mud, dust, and dirt and in integrated circuits as Hall effect sensors. The heat-generated effects may alter the temperature distribution and, consequently, alter the rate of particle deposition in a nuclear reactor, electronic chips, and semiconductor wafers [13]. The magnitude of thermal diffusion effect may be considerably large, and it can be utilized as such in the separation of isotopes in a mixture of gases with very light molecular weight (hydrogen and helium) and gases with medium molecular weight (nitrogen and air). Mahmud and Mohammad [14] studied MHD fluid flow of an incompressible electrically conducting fluid along a semi-infinite vertical porous plate under strong transverse magnetic field with rotation. The momentum and energy equations were solved by the perturbation technique. Das et al. [15] investigated the effects of slip condition on an unsteady MHD flow of a viscous incompressible electrically conducting fluid past a periodically accelerated horizontal porous plate under the influence of transverse magnetic field and Hall current. The flow was generated due to the accelerated motion of the porous plate. A unified closed form analytical solution of governing equations was obtained by the Laplace transform technique. Prasad et al. [16] investigated the MHD flow and heat transfer of an electrically conducting fluid over a stretching sheet with variable thickness. The effects of external magnetic field along the sheet and the Hall currents were considered too. The governing equations were solved numerically using an implicit finite difference scheme. The results showed the effect of the Hall current which gave rise to a cross flow. Moreover, the Hall current and the magnetic field had a strong effect on the flow and heat transfer characteristics. Ram [17] employed the finite difference method to analyse the MHD Stokes problem for a vertical plate with Hall and ion-slip currents and discussed MHD Stokes problem of a convective flow of a vertical infinite plate in a dissipative rotating fluid with Hall current for the effects of various parameters.

Free convection in channel flow has many important applications in designing ventilating and heating of buildings, cooling of electronic components of a nuclear reactor, bed thermal storage, and heat sink in the turbine blades. Convective flows driven by temperature difference of the bounding walls of channels are important in industrial applications [18, 19]. Soundalgekar [3] studied magnetohydrodynamic free convection flow past an infinite vertical plate oscillating in its own plane. An isothermal plate was considered, and it was observed that a phase angle of brings about cooling of the plate making the flow unstable. Bhaskar et al. [20] studied the effects of Hall currents and radiation on MHD free convection of a viscous incompressible fluid confined between two vertical walls in a rotating medium. The governing equations were solved analytically using the Laplace transform technique. Mahendar and Rao [21] investigated the effect of Hall current on an unsteady MHD free convective flow along a vertical porous plate with thermal radiation. The dimensionless governing coupled, nonlinear boundary layer partial differential equations were solved by an efficient finite element method for velocity, temperature, concentration field, and also the expression for shearing stress, Nusselt number, and Sherwood number. MHD free convective heat and mass transfers of a heat generating fluid past an impulsively started infinite vertical porous plate with Hall current and radiation absorption was studied in [22]. The effects of Hall currents on MHD free convective flow past a semi-infinite vertical plate with mass transfer for different parameters was studied in [23].

Although different authors studied heat and mass transfers on different rotating geometric structures with or without viscous dissipation and porosity applying different numerical techniques, this study concentrates on corrugated vibrating surface with hall currents and heat and mass transfers that is sinusoidal in nature and nonporous bottom by applying the FDM explicit scheme because of stability.

2. Mathematical Formulation

Consider an unsteady MHD flow of an electrically conducting incompressible fluid on a vibrating bottom surface that is sinusoidal as shown in Figure 1. The -axis is taken along the infinite surface and the -axis normal to it. Initially the corrugated surface and the fluid are at the same temperature, and fluid is allowed to flow in the upward direction. At time , the plate starts vibrating with a frequency of oscillation and reference velocity forming boundary velocity , where is the reference velocity having a linear relation with the vibration amplitude and the vibration frequency. A magnetic field of uniform strength is applied perpendicular to the plate along the positive -axis. The gravitational acceleration is acting downwards.

The magnetic Reynolds number is assumed to be small since the electric intensity is zero at the plate; therefore, it is assumed to be zero everywhere within the flow [24]. Considering the equation of conservation of electric charge, [25] given as , the components of electric current density along the , , and -axes are , , and , respectively. Since the bottom vibrating surface is electrically nonconducting, and is zero everywhere in the flow [26]. Therefore, considering Hall current in the problem, the generalized Ohm’s law [25] is given as

For the electron pressure gradient for weakly ionized fluid in equation (1), ion-slip and thermoelectric effect are neglected. Taking these assumptions into consideration and in the absence of electric field, it simplifies towhere is the Hall parameter. By solving equations (2) and (3) for and , it gives

Applying Lorentz force yields

If there is an appreciable temperature difference between the vibrating surface and the ambient fluid, the temperature-dependent heat source or sink is considered which may exert a strong influence on the heat transfer characteristics within the fluid [27]. Moreover, if is the temperature-dependent volumetric heat generation parameter, then represents a heat source, and represents a heat sink. Incorporating equations (2) to (4) into the equations of momentum, the following respective equations of momentum, energy, and concentration are obtained:

Equations of momentum:

Energy equation:

Equation of concentration:

The initial and boundary conditions for the current problem are specified as follows:

3. Nondimensionalization

The nondimensional equations of momentum, energy, and concentration are transforming the equations governing the flow to their respective nondimensionalised form given aswhere , the nondimensional heat, is the source or sink parameter [28].

The nondimensional form of the initial and boundary conditions for this problem is given as

The equations that govern the flow on a vibrating corrugated surface in the presence of heat source or sink are coupled and highly nonlinear. They are represented by equations (10)–(12) and are to be solved under the initial and boundary conditions (13).

4. Numerical Technique

The numerical method used to solve the equations governing the flow together with the initial and boundary conditions applies the forward-time backward-space (FTBS) difference method which is stable and convergent. A necessary condition for time stability is the CFL condition which depends on time and space discretization.

The finite difference form of the equations of momentum (14) and (15), concentration (16), and energy) (17) are given as

A mesh size is fixed at , , and , where a sinusoidal shape is formed by the -axis and -axis considered appropriate for calculation of velocities, temperature, and concentration. The FTBS scheme is unconditionally stable, and thus compatibility and stability ensures convergence. Hence, the present employed scheme is always convergent (Figure 2).

5. Discussion of the Results

In order to ascertain the accuracy of the numerical results, the present study is compared with the available solution in literature. The velocity, temperature, and concentration fields have been discussed by assigning numerical values to various nondimensional parameters, that is , , , , , , , , , , and . The values of the Schmidt number used are 0.22, 0.62, and 0.78 and represent hydrogen, water vapour, and ammonia, respectively, at and one atmosphere pressure. It is clear from Figures 1 and 3 that, in the presence of a heat source, the increase in Hall parameter leads to an increase in the magnitude of primary and secondary velocity profiles. The Hall parameter increases the transverse velocity profile up to a value ; however, for the values , the transverse velocity profiles decreases since the large values of turn very small, and hence the resistive effect of magnetic field diminishes. This is due to the fact that an increase in decreases the effective conductivity and hence the magnetic damping.

An increase in the magnetic parameter leads to a decrease in the magnitude of primary velocity profiles and secondary velocity profiles, respectively. A large value of causes a reversal in the direction of the secondary velocity profiles. The effect of results to a decreasing velocity distribution across the boundary layer due to the effect of transverse magnetic field which gives rise to Lorentz force that has a tendency to flow down the fluid motion. A higher increases the transverse velocity field and reaches a peak value near the vicinity of the boundary layer and approaches zero, and the less magnetic field has less peak value as it reaches zero far away from the boundary.

It is clear from Figure 4 that, for lower values of , it decreases the concentration distribution in the boundary layer thickness, which is associated with the reduction in the concentration profiles. The particles have greater diffusivity conditioned by viscosity of the medium for them to cross from laminar to turbulent-characterized fluid flows in which there are simultaneous momentum and mass diffusion convection processes. This means that an increase in means a decrease of molecular diffusion; hence, the concentration of the species is higher for large values of and lower for small values of . It is noticed that the effect of increasing values of results in a decreasing concentration profiles across the boundary layer. Physically, the rise in the value of Schmidt number means reduction of molecular diffusion. Hence, the concentration of the species is advanced for smaller values of and lesser for larger values of Sc.

It is observed from Figure 5 that an increase in the Prandtl number decreases the velocity of the flow field. Since Prandtl number is the ratio of kinematic viscosity to thermal diffusivity, this implies that an increase in the kinematic viscosity of the fluid dominates the thermal diffusivity of the fluid which leads to decreasing the velocity of the flow field.

Figure 6 presents the effect of the Prandtl number on the temperature profiles as a result of the heat source (). An increase in the value of has the tendency to decrease the fluid temperature in the boundary layer as well as the thermal boundary layer thickness. This causes the wall slope of the temperature to decrease as increases. This means that a mixture of flue gas () and hydrogen () diffuses faster than does a mixture of air and hydrogen; and a mixture of air and hydrogen diffuses faster than a mixture of water and hydrogen. Figure 5 shows that, in the presence of a heat source, an increase in the Schmidt number causes a decrease in the velocity profiles.

A similar condition is experienced as shown in Figure 7, where an increase over a given time causes a decrease in concentration profile.

Figures 8 and 9 show that an increase in the magnitude of the heat sink parameter () leads to a decrease in the magnitude of temperature and velocity profiles, respectively. A contrary effect is shown in Figure 8, where an increase in leads to an increase in concentration profiles since the presence of a heat sink produces a cooling effect that decreases velocity of the convection currents that move upwards next to the surface of the vibrating surface, leading to higher concentration profiles.

Heat source is an indication that generation of heat from the surface (this is due to ), which increases the temperature in the flow field. This physically implies that when the heat source parameter increased, the temperature increases steeply and exponentially from the surface. The influence of heat source parameter on velocity and temperature profiles is very much significantly related to the heat sink parameter (Figure 10).

As changes, separation of flows in the weak region of the corrugated vortices is created. At lower , no separation occurs, and an increase in unstable separations are seen as a result of vortex shedding at certain frequencies, and it applies that the viscous force is predominant, thus imposing drug in the fluid and thus decreasing the fluid flow as shown in Figure 11. The effect of the flow velocity on the rod vibration can be represented by the Reynolds number, which expresses the extent of the flow turbulence.

The Nusselt number physically indicates the rate of heat transfer and the local Sherwood number as well as the rate of mass transfer. The shearing stress on the surface of the vibrating surface is defined. The local skin friction coefficients due to the primary and the secondary velocity profiles are given as

The length units along a vibrating surface is expressed in Sherwood number as

Furthermore, the local Nusselt number Nu may be expressed as

Heat and mass transfers are considered. This is done by computing velocity, concentration, and temperature in equations (18) to (21) in MATLAB code to generate the values of skin friction, Nusselt number, and Sherwood number as shown in Table 1.

It is clear from the table that(i)As time increases, a decrease in Nusselt number, Sherwood number, and Skin friction occurs as it physically implies that shear stresses decrease with an increase in time.(ii)An increase in the heat source parameter leads to an increase in both skin frictions, meaning that an increase in enhances convection currents on the flow and a reverse effect in Sherwood number and Nusselt number since a thicker thermal boundary layer leading to lower temperature gradients.(iii)A high value of radiation parameter leads to an increase in magnitudes of skin frictions and Nusselt number as a result of an increase in the rate of species concentration. It is clear that the effect of radiation is to decrease the rate of energy transport to the fluid, thereby decreasing the temperature of the fluid but decreasing the Sherwood number.

6. Conclusion

A numerical study has been conducted on free convective heat and mass transfers of an incompressible electrically conducting fluid on a vibrating bottom surface with Hall currents and heat and mass transfers. It is clear that(i)The velocity profiles decreases with an increase in the Prandtl number , Schmidt number , Magnetic parameter , and Soret number .(ii)Concentration profiles increase with an increase in the Reynolds number and decreases with an increase in the Prandtl number and Schmidt number .(iii)Both hydrodynamic and thermal behaviours are relevant to the Reynolds number, and the patterns are shown by oscillations.(iv)Fluid damping is generated as the bottom surface vibrates normal to the flow, and this is proportional to the surface velocity. Therefore, an increase in damping causes an increase in the velocity.(v)It is shown that an increase in the parameters under consideration causes reduction in the skin friction.

Nomenclature

:	Magnetic field strength, Wbm⁻¹
:	Magnetic flux density, Wbm⁻²
:	Schmidt number
:	Soret number
:	Dufour number
:	Molecular diffusion coefficient
:	Magnetic parameter
:	Dimensional temperature of the fluid, K
R:	Joules heating parameter
N:	Radiative parameter
:	Acceleration due to gravity (m·s⁻²)
:	Prandtl number
:	The Eckert number
:	Electric flux density, (C·m⁻²)
:	Magnetic field intensity vector, (A·m⁻¹)
L:	Characteristic length (m)
:	Phase angle
:	Local temperature Grashof number
:	Local mass Grashof number
:	Frequency of oscillation
:	Fluid temperature far from the surface of the object) (K)
:	Wall surface temperature (K)
:	Concentration in the fluid away from the surface
:	Hydromagnetic Reynold’s number
:	Electric field strength, V·m⁻¹
:	Magnetic parameter
:	Dimensional temperature of the fluid, K
:	Acceleration due to gravity (m·s⁻²)
, :	Dimensionless velocity components
, :	Cartesian coordinates
: :	Nondimensional velocity components in and directions.
:	Dimensionless heat source parameter
:	Volumetric heat generating parameter (W·m⁻³)
:	Concentration susceptibility parameter (kmol·m⁻³)
:	Specific heat at constant pressure (J·Kg⁻¹·K⁻¹)
:	Hall parameter
:	Charge density, (C·m⁻³)
:	The local skin friction coefficient due the primary velocity profiles.
:	The local skin friction coefficient due the secondary velocity profiles.
:	Viscous stress due to primary velocity profiles, (N·m⁻²)
:	Viscous stress due to secondary velocity profiles, (N·m⁻²)
:	Gradient operator

Greek Symbols

:	Thermal expansion coefficient (K⁻¹)
:	Electrical conductivity, ⁻¹ m⁻¹
:	Density (kg·m⁻³)
:	Electron frequency (s⁻¹)
:	Electron collision time (t)
:	Stefan–Boltzmann constant
:	Magnetic permeability, (H·m⁻¹)
:	Thermal expansion coefficient
:	Mass expansion coefficient
:	Dynamic viscosity (m²·s⁻¹)
:	Kinematic viscosity (m²·s⁻¹)
:	Thermal diffusivity (m²·s⁻¹).

Data Availability

The data was generated by MATLAB version 8.5.0.197613 (R2015a) computer program.

Disclosure

This research work is an academic paper as a requirement for PhD completion, therefore self-sponsored by the student.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

Y. Perets, B. Rivin, I. Bronstein, R. Harari, and E. Sher, “Flow induced vibration of corrugated rod subjected to axial flow,” WIT Transactions on the Built Environment, vol. 105, pp. 129–133, 2009.
View at: Publisher Site | Google Scholar
M. Amielh, F. Anselmet, Y. Jiang et al., “Aeroacoustic source analysis in a corrugated flow pipe using low-frequency mitigation,” Journal of Turbulence, vol. 15, no. 10, pp. 650–676, 2014.
View at: Publisher Site | Google Scholar
O. Veli, T. Nehir, and S. Besir, “Flow characteristics and heat transfer enhancement in 2DCorrugated channels,” International Journal of Mechanical, Aerospace, Industrial and Mechatronics Engineering, vol. 7, 2013.
View at: Google Scholar
V. M. Soundalgekar, “Free convection effects on the flow past an infinite vertical oscillating plate,” Astrophysics and Space Science, vol. 64, no. 1, pp. 165–171, 1979.
View at: Publisher Site | Google Scholar
R. B. Bird and W. E. S. Lightfoot, Transport Phenomena, Wiley, New York, NY, USA, 1966.
V. Garg and P. Maji, “Flow and heat-transfer in a sinusoidally curved channel,” International Journal of Engineering Fluid Mechanics, pp. 293–319, 1988.
View at: Google Scholar
J. A. Gbadeyan, T. L. Oyekunle, P. F. Fasogbon, and J. U. Abubakar, “Soret and Dufour effects on heat and mass transfer in chemically reacting MHD flow through a wavy channel,” Journal of Taibah University of Science, vol. 12, pp. 630–651, 2018.
View at: Publisher Site | Google Scholar
E. M. Abo-Eldahab and E. Aziz, “Hydrodymagnetic three-dimensional flow over a stretching surface with heat and mass transfer,” Heat and Mass Transfer, vol. 41, no. 8, pp. 734–743, 2005.
View at: Publisher Site | Google Scholar
A. S. Idowu, A. Jimoh, and L. O. Ahmed, “Impact of heat and mass transfer on MHD oscillatory flow of Jeffery fluid in a porous channel with thermal conductivity, Dufour and Soret,” Journal of Applied Sciences Environment Management, vol. 19, no. 4, pp. 819–830, 2015.
View at: Google Scholar
L. Zhang and D. che, “Turbulence models for fluid flow and heat transfer between cross- corrugated plates,” Numerical Heat Transfer, Part A: Applications, vol. 60, no. 5, pp. 410–440, 2011.
View at: Publisher Site | Google Scholar
S. K. Tania and M. A. Samad, “Effects of radiation, heat generation and viscous dissipation on MHD free convection flow along a stretching sheet,” Research Journal of Applied Sciences, Engineering and Technology, vol. 2, pp. 368–377, 2010.
View at: Google Scholar
J. O’Brien and E. Sparrow, “Corrugated-duct heat transfer, pressure drop, and flow visualization,” ASME Journal of Heat Transfer, vol. 104, no. 1, pp. 410–416, 1982.
View at: Publisher Site | Google Scholar
L. Goldstein and E. M. Sparrow, “Heat/mass transfer characteristics for flow in a corrugated wall channel,” Journal of Heat Transfer, vol. 99, no. 2, pp. 187–195, 1977.
View at: Publisher Site | Google Scholar
A. Mahmud and R. Mohammad, “MHD fluid flow of an incompressible electrically conducting fluid along a semi-infinite vertical porous plate under strong transverse magnetic field with rotation,” in Proceedings of the 8th BSME International Conference on Thermal Engineering, Dhaka, Bangladesh, December 2019.
View at: Publisher Site | Google Scholar
S. Das, B. Tarafdar, and R. N. Jana, “Hall effects on unsteady MHD rotating flow past a periodically accelerated porous plate with slippage,” European Journal of Mechanics—B/Fluids, vol. 72, pp. 135–143, 2018.
View at: Publisher Site | Google Scholar
K. V. Prasad, K. Vajravelu, and H. Vaidya, “Hall effect on MHD flow and heat transfer over a stretching sheet with variable thickness,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 17, no. 4, pp. 288–297, 2016.
View at: Publisher Site | Google Scholar
P. C. Ram, “Hall effects on hydromagnetic convective flow in a rotating fluid through porous medium,” Journal of Engineering Science, vol. 5, pp. 43–53, 1991.
View at: Google Scholar
P. Dulal and T. Babulal, “Influence of Hall current and thermal radiation on MHD convective heat and mass transfer in a rotating porous channel with chemical reaction,” International Journal of Engineering Mathematics, Article ID 367064, pp. 2013–13, 2013.
View at: Publisher Site | Google Scholar
M. Gradeck, B. Hoareau, and M. Lebouché, “Local analysis of heat transfer inside corrugated channel,” International Journal of Heat and Mass Transfer, vol. 48, no. 10, pp. 1909–1915, 2005.
View at: Publisher Site | Google Scholar
C. Bhaskar, D. Sanatan, and R. J. Rabindra, “The Effects of Hall currents and radiation on MHD free convection of a viscous incompressible fluid confined between two vertical walls in a rotating,” Advances in Applied Science Research, vol. 3, pp. 3766–3787, 2012.
View at: Google Scholar
D. Mahendar and P. S. Rao, “The effect of hall current on an unsteady MHD free convective flow along a vertical porous plate with thermal radiation,” IOSR Journal of Mathematics, vol. 11, pp. 122–141, 2015.
View at: Google Scholar
M. Kinyanjui, J. K. Kwanza, and S. M. Uppal, “MHD free convection heat and mass transfer of a heat generating fluid past an impulsively started infinite vertical porous plate with Hall current and radiation absorption,” Energy Conversion and Management, vol. 42, no. 8, pp. 17–931, 2001.
View at: Publisher Site | Google Scholar
E. M. AboEldahab, “Convective heat transfer, by the presence of radiation, in an electrically conducting fluid at a stretching surface,” Canadian Journal of Physics, vol. 79, no. 6, pp. 929–937, 2001.
View at: Publisher Site | Google Scholar
K. Giterere, M. Kinyanjui, and S. M. Uppal, “MHD flow in porous media over a stretching surface in rotating system with heat and mass transfer,” International Electronic Journal of Pure and Applied Mathematics, vol. 4, no. 1, pp. 9–32, 2012.
View at: Google Scholar
R. C. Meyer, “On reducing aerodynamic heat-transfer rates by magnetohydrodynamic techniques,” Journal of the Aerospace Sciences, vol. 25, no. 9, pp. 561–566, 1958.
View at: Publisher Site | Google Scholar
N. Ida, “Faraday’s law and induction,” in Engineering Electromagnetics, Springer, New York, NY, USA, 2000.
View at: Google Scholar
M. Q. Brewster, Thermal Radiative Transfer and Properties, John Wiley & Sons, New York, NY, USA, 1992.
W. M. Rohsenow, J. P. Hartnett, and Y. I. Cho, Handbook of Heat Transfer, McGraw-Hill, New York, NY, USA, 3rd edition, 1998.

Copyright

Copyright © 2020 Vincent M. Bulinda 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.

PDF Download Citation

Download other formats

Order printed copies