Next Article in Journal
Analysis and Experiment for Wireless Power Transfer Systems with Two Kinds Shielding Coils in EVs
Next Article in Special Issue
CFD Analysis of Falling Film Hydrodynamics for a Lithium Bromide (LiBr) Solution over a Horizontal Tube
Previous Article in Journal
A Systematic Approach to Predict the Economic and Environmental Effects of the Cost-Optimal Energy Renovation of a Historic Building District on the District Heating System
Previous Article in Special Issue
Effect of Viscous Dissipation in Heat Transfer of MHD Flow of Micropolar Fluid Partial Slip Conditions: Dual Solutions and Stability Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heat and Mass Transfer in Hydromagnetic Second-Grade Fluid Past a Porous Inclined Cylinder under the Effects of Thermal Dissipation, Diffusion and Radiative Heat Flux

1
Department of Mathematics, Air University, PAF Complex E-9, Islamabad 44000, Pakistan
2
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam
3
Center of Excellence for Research in Engineering Materials (CEREM), King Saud University, P.O. Box 800, Al-Riyadh 11421, Saudi Arabia
4
Electrochemistry and Corrosion Laboratory, Department of Physical Chemistry, National Research Centre, El-Behoth St. 33, Dokki, Cairo 12622, Egypt
*
Author to whom correspondence should be addressed.
Energies 2020, 13(1), 278; https://doi.org/10.3390/en13010278
Submission received: 29 October 2019 / Revised: 2 December 2019 / Accepted: 30 December 2019 / Published: 6 January 2020
(This article belongs to the Special Issue Fluid Flow and Heat Transfer Ⅱ)

Abstract

:
Current disquisition is presented to excogitate heat and mass transfer features of second grade fluid flow generated by an inclined cylinder under the appliance of diffusion, radiative heat flux, convective and Joule heating effects. Mathematical modelling containing constitutive expressions by obliging fundamental conservation laws are constructed in the form of partial differential equations. Afterwards, transformations are implemented to convert the attained partial differential system into ordinary differential equations. An implicit finite difference method known as the Keller Box was chosen to extract the solution. The impact of the flow-controlling variables on velocity, temperature and concentration profiles are evaluated through graphical visualizations. Variations in skin friction, heat transfer and mass flux coefficients against primitive variables are manipulated through numerical data. It is inferred from the analysis that velocity of fluid increases for incrementing magnitude of viscoelastic parameter and curvature parameter whereas it reduces for Darcy parameter whereas skin friction coefficient decreases against curvature parameter. Assurance of present work is manifested by constructing comparison with previous published literature.

1. Literature Survey

Two types of fluids, namely Newtonian and non-Newtonian fluids, have changed the dynamics of today’s world. Initially Newtonian fluids gained pervasive focus from researchers, but due to the complexity and diversification in various fluid-dependent procedures, the essence of non-Newtonian fluids is rising. Non-Newtonian fluids have venerable significance and are required for different commercial, industrial, physiological, mechanical, physical, medicinal and technological purposes. So, for a better interpretation of such processes, narration of non-Newtonian fluids is essential. It can be done by comparing the properties of Newtonian and non-Newtonian fluids. One major difference between these two fundamental classes of fluids is that Newtonian fluid exhibits the linear relation among shear and strain rates whereas non-Newtonian fluids manifest a non-linear relationship. Due to such a relation, non-Newtonian fluids disclose intricate properties and are characterized into subdivisions like shear thinning, shear thickening, dilatant and thixotropic fluids. Rheologist’s have recognized that some fluid models which behave multiply under different conditions like second-grade fluid that explicates features of shear thinning, shear thickening and Newtonian fluid. Such dynamical features of second-grade fluid has made it notable and worthy among researchers. The thermal and dynamical analyses of second-grade fluid was done by Vejravelu and Roper [1]. Rajeswari and Rantha [2] capitalized the perturbation technique to discuss the dynamics of stagnant second-grade fluid. Garg and Rajagopal [3] followed the work conducted in Reference [2] and found total agreement with their findings. Fetecau and Fetecau [4] determined solution for second-grade fluids due to an oscillating circular cylinder by using Hankel transform method). A wealth of knowledge about the analysis of second-grade fluid over various physical configurations and conditions can be accessed through References [5,6,7,8].
Optimal production from the industrial sector requires appropriate collaboration of thermal and concentration fields. This can be done either by providing radiation, Newtonian heating, Dufour and Soret aspects, mixed convection, as well as many other ways. So, we discuss separately the abovementioned thermal aspects and their physical implications. Newtonian heating is defined as the process in which internal resistance is at a minute scale in comparison to the surface resistance. Newtonian heating is involved in thermal exchange procedures like heat exchangers, transferal of heat around fins, the oil and gas industry and many more. Merkin [9] performed premier work by considering different types of wall thermal conditions, and convective heating was one of them. Salleh et al. [10] probed changes in the thermal boundary layer and the associated heat transfer aspect by obliging Newtonian heating. Haq et al. [11] and Nadeem et al. [12] scrutinized the 2D and 3D flow of viscoelastic fluid under the existence of convective boundary sources. Innovative work in this direction is conducted and cited in [13,14,15,16,17].
In advance manufacturing and thermal procedures, the capability to generate chaotic energy emission is vital. Radiative energy flux has a conclusive impact on chaotic heat energy management in various engineering and physical processes. One of the latest investigations emphasizing the impact of radiation was conducted by Kothandapani and Parakash [18], who interpreted the impact of radiative heat transfer in peristaltic flow through tapered asymmetric arteries. The behavior of the pseudoplastic fluid model along with the suspended particles by obliging non-linearized radiative flux due to a stretchable surface was elucidated by Kumar et al. [19]. Khan et al. [20] disclosed the influence of radiative energy on a magnetized Carreau nano fluid over a convective thermal sheet. Waqas et al. [21] manipulated the modelling of a Carreau fluid by obliging the Roseland approximation. Goodarzi et al. [22] explicated the effect of radiative heat flux on laminar and turbulent mixed-convection heat transfer of a semitransparent medium in a square enclosure by implementing the finite volume method. Pordanjani et al. [23] studied the impact of radiation on the convection heat transfer rate and the nanofluid entropy generation within a diagonal rectangular chamber numerically, in the presence of a magnetic field. Aghaei et al. [24] scrutinized heat transfer aspects by placing horizontal and vertical elliptic baffles in the flow in the presence of entropy generation of multi wall carbon nano tubes MWCNTs-water nanofluid.
Thermal variations by way of radiation, viscous dissipation, convective heating across boundaries concentration and temperature gradients are generated, which yields mass/energy fluxes, respectively. An energy flux termed as diffusion-thermo (the Dufour effect) is produced due to a concentration gradient, whereas mass flux is known as thermo-diffusion (the Soret effect), and which is induced by way of temperature gradients. Such aspects are required in heat exchangers, cooling applications, steel industries and so forth. So, the investigation of flow processes under the effect of Dufour and Soret aspects are highly recommendable and relevant to researchers. In view of relevancy, a lot of information of such aspects has been gathered, which is viewed in the accompanied references. Hirshfelder et al. [25] measured the diffusion coefficient of mono atomic gases. Afify [26] investigated fluid flow generated by a stretching sheet in the presence of Dufour and Soret aspects. Bhattacharya et al. [27] investigated stagnant flow over an extendable sheet with diffusion effects. Similar to the analysis in [28], Awad et al. premediated the double diffusion (Dufour and Soret) effects over a stretching sheet by way of a computational scheme. For the power-law model, diffusion effects were elucidated by Goyal et al. [29]. Umar et al. [30] disclosed transportation of nano particles along with diffusion effects over a porous configuration. Hayat et al. [31] adumbrated the heat transfer of Newtonian fluids over an inclined stretching cylinder. Bagherzadeh et al. [32] experimentally verified the enhancement in heat transfer rate by jet injection of a nanofluid under a homogeneous magnetic field. Goshayeshi et al. [33] elaborated the effect of a ferro nanofluid on the enhancement of the thermal features of a flowing fluid within the heat pipe under the implication of a magnetic field. Goshayeshi et al. [34] experimentally justified the impact of a Fe2O3/kerosene nanofluid on fluid flow in a heated pipe to enhance the thermal performance as well as the heat transfer coefficient. Yousefzadeh et al. [35] probed laminar mixed-convection heat transfer inside an open square cavity with different heat transfer areas by finding a computational solution. Tian et al. [36] experimentally measured heat transfer in a silica DI water nano-fluid with three various surfactants on the surface of copper heaters at different concentrations at atmospheric pressure. Forced convection in a double-tube heat exchanger using nanofluids with constant and variable thermophysical properties was scrutinized by Bahman et al. [37].
The prime concern of this investigation is to anticipate the radiative flow of a second-grade fluid with Joule heating, as well as Dufour and Soret aspects over an inclined stretching cylinder with Newtonian heating. Model formulation is controlled through coupled partial differential equations and afterwards converted to ordinary differential equations by permissible transmutations. Findings are attained by implementing the Keller Box scheme. Tabular and pictorial representations regarding the impact of concerning parameters on sundry distributions are disclosed. To the best of the author′s knowledge, this work has not been done so far and it will serve as a reference study for upcoming research in this direction.

2. Mathematical Modelling

Let us consider the steady, axisymmetric and incompressible flow of a second-grade fluid flowing over an inclined stretchable cylinder under the appliance of a transverse magnetic field.
The tensor describing the flow analysis for concerning fluid is presented as follows:
T = p I + μ A 1 + α 1 A 2 + α 2 A 1 2 .
In the above expression α 1 and α 2 are the material moduli and A 2 is given by
A 2 = d A 1 d t + A 1 ( g r a d   V ) + ( g r a d   V ) T A 1 ,
where we can write
A 2 = [ 4 ( u r ) 2 + 2 v r ( u x + w r ) 0 2 2 u r x + 2 2 w r x + w x ( u x + w r ) + u r ( u x + w r ) 0 4 u 2 r 2 0 2 2 u r x + 2 2 w r x + w x ( u x + w r ) + u r ( u x + w r ) 0 4 ( w x ) 2 + 2 u x ( u x + w r ) ] .
Thus, the equations for motion, energy and concentration are expressed as
u u x + w u r   = ν ( 2 u r 2 + 1 r   u r ) + g β ( T T ) C o s ( α ) + g β c ( C C ) C o s ( α ) ν K u + α 1 ρ ( w 3 u r 3 + u 3 u x 2 r + u x 2 u r 2 u r 2 w r 2 + 1 r ( w 2 u r 2 + u 2 u r x + u r u x u r w r ) ) σ B 0 2 ρ u ,
u T x + w T r   =   α ( 1 + 4 3 R ) 1 r r ( r T r ) + μ ρ c p ( u r ) 2 + D m k T c p c s 1 r r ( r C r ) + α 1 ρ c p ( w u r 2 u r 2 +     u u r 2 u x r ) + σ B 0 2 ρ c p u 2 ,
u C x + w C r   =   D m 1 r r ( r C r ) +   D m k T T m   1 r r ( r T r ) .
The present flow situation is subjected to the following boundary conditions,
{ u   =   U o x l , w = 0 ,   k T r   =   h ( T f T ) , D m C r   =   k m ( C f C ) ,     a t   r   =   a   u 0 ,     T T ,     C C ,       a s     r       .
Invoking the similarity variables,
u   =   U o x l f ( η ) , w = ( ν U o l ) f ( η ) , θ ( η ) = T T T f T , ϕ ( η ) = C C C f C , η = r 2 a 2 2 a U o ν l .
Employing the transformations, Equations (4)–(6) taken as the form, we get
( 1 + 2 γ η ) f + 2 γ f + f f f 2 D a f + 4 γ β ( f f f f ) + β ( 1 + 2 γ η ) ( 2 f f + f 2 f f ) H a 2 f + G r θ C o s ( α ) + G c ϕ C o s ( α ) = 0 ,
( 1 + 4 3 R ) ( ( 1 + 2 γ η ) θ + 2 γ θ ) + P r f θ + D u ( ( 1 + 2 γ η ) ϕ + 2 γ ϕ ) P r E c β γ f f 2 + P r E c ( 1 + 2 γ η ) ( f 2 β f f f + β f f 2 ) + P r E c H a 2 f 2 = 0 ,
( 1 + 2 γ η ) ϕ + 2 γ ϕ + S c f ϕ + S r ( ( 1 + 2 γ η ) θ + 2 γ θ ) = 0 .
Here the derivatives of the function are indicated by a prime with respect to η , and β   = α 1 U o ρ ν l and H a = σ B 0 2 l ρ U o are the viscoelastic parameter and Hartmann number, respectively.
{ f = 0 ,     f = 1 ,   θ ( 0 ) = B i 1 ( 1 θ ( 0 ) ) ,   ϕ ( 0 ) = B i 2 ( 1 ϕ ( 0 ) ) ,                 a t           η = 0 f = 0 ,           θ = 0 ,           ϕ = 0 ,           a t   η = .  
in which
γ = ν l a 2 U o   , G r = g β ( T f T ) l 2 x U o 2 , G c = g β c ( C f C ) l 2 x U o 2 , D a = ν l K U o , P r = ν α , D u = D m k T c p c s ρ c p k C f C T f T , S r = k T ( T f T ) T m ( C f C ) , S c = ν D m , E c = U o 2 ( x l ) 2 c p ( T f T ) , B i 1 = h k ν l U o   , B i 2 = k m D m ν l U o , R = 4 σ T 3 k k ,
where γ is the curvature parameter, D a is the Darcy number, D u is the Dufour number, S r is the Soret number, G r is the thermal Grashof number, G c is the solutal Grashof number, R is the radiation parameter, and B i 1 and B i 2 are the thermal and concentration Biot numbers, respectively.
Nusselt and Sherwood mass flux coefficients are expressed as N u x and S h x , and with the skin-friction co-efficient C f x , are the physical quantities of interest expressed as
N u x = x q w k ( T f T ) ,     q w = k T r 16 σ T 3 3 k T r | r = a ,
S h x = J w x D m ( C f C ) ,     J w = D m C r | r = a ,
C f x = 2 τ w ρ u w 2 ,     u w = U 0 ( x l )   a n d   τ w = μ u r | r = a .
In dimensionless form
1 2 C f x R e x 1 / 2 = f ( 0 ) ,
N u x / R e x 1 2 = ( 1 + 4 3 R ) θ ( 0 ) ,
S h x / R e x 1 2 = ϕ ( 0 ) .

3. Solution Methodology

The problem comprising Equations (9)–(11) is handled numerically by way of the Keller Box scheme. Firstly, Equations (9)–(11) are reduced to the initial value problem by utilizing a new variable given by
f = u ,   u = v ,   θ = p ,   ϕ = g ,
then, Equations (9)–(11) become
( 1 + 2 γ η ) v + 2 γ v + f v u 2 D a   u + G r θ C o s ( α ) + G c ϕ C o s ( α ) = 0 ,
( 1 + 4 3 R ) ( ( 1 + 2 γ η ) p + 2 γ p ) + P r f p + D u ( ( 1 + 2 γ η ) g + 2 γ g ) + P r E c ( 1 + 2 γ η ) v 2 = 0 ,
( 1 + 2 γ η ) g + 2 γ g + S c f g + S r ( ( 1 + 2 γ η ) p + 2 γ p ) = 0 ,
along with the B. C′s
{ f ( 0 ) = 0 ,   u ( 0 ) = 1 ,   p ( 0 ) = B i 1 ( 1 θ ( 0 ) ) ,   g ( 0 ) = B i 2 ( 1 ϕ ( 0 ) ) , u ( ) = 0 ,   θ ( ) = 0 ,   ϕ ( ) = 0 .
By using the central difference gradients and the average at the mid points of the net derivatives, the following relation is approximated:
x 0 = 0 ,   x n = x n 1 + k n ,     n = 1 , 2 , , J ,
η 0 = 0 ,   η j = η j 1 + h j ,     j = 1 , 2 , , J   n J = n .
Equations (20)–(23) then become
f j f j 1 h j   = u j 1 2   , u j u j 1 h j   = v j 1 2   , θ j θ j 1 h j   = p j 1 2   , ϕ j ϕ j 1 h j   = g j 1 2 ,
( 1 + 2 γ η ) ( v j v j 1 h j ) + 2 γ ( v j 1 2 ) + ( f j 1 2 ) ( v j 1 2 ) ( u j 1 2 ) 2 D a ( u j 1 2 ) + G r C o s ( α ) ( θ j 1 2 ) + G c C o s ( α ) ( ϕ j 1 2 ) = 0 ,
( 1 + 4 3 R ) ( ( 1 + 2 γ η ) ( p j p j 1 h j ) + 2 γ ( p j 1 2 ) ) + P r ( f j 1 2 ) ( p j 1 2 ) + D u ( ( 1 + 2 γ η ) ( g j g j 1 h j ) + 2 γ ( g j 1 2 ) ) + P r E c ( 1 + 2 γ η ) ( v j 1 2 ) 2 = 0 ,
( 1 + 2 γ η ) ( g j g j 1 h j ) + 2 γ ( g j 1 2 ) + S c ( f j 1 2 ) ( g j 1 2 ) + S r ( ( 1 + 2 γ η ) ( p j p j 1 h j ) + 2 γ ( p j 1 2 ) ) = 0 ,
where f j 1 2 = f j + f j 1 2 , etc.
Equations (26)–(28) are in algebraic form and therefore have to be linear equations before the implementation of factorization. We write the Newton iterations in the following way:
f j ( i + 1 ) = f j ( i ) + δ f j ( i ) ,   u j ( i + 1 ) = u j ( i ) + δ u j ( i ) ,   v j ( i + 1 ) = v j ( i ) + δ v j ( i ) ,
θ j ( i + 1 ) = θ j ( i ) + δ θ j ( i ) ,   ϕ j ( i + 1 ) = ϕ j ( i ) + δ ϕ j ( i ) ,   p j ( i + 1 ) = p j ( i ) + δ p j ( i ) ,     g j ( i + 1 ) = g j ( i ) + δ g j ( i ) .
By substituting these expressions into Equations (26)–(28) and neglecting the second-order and higher-order terms in δ , a linear tri-diagonal system of equations will be found as follows:
δ f j δ f j 1 h j 2 ( δ u j + δ u j 1 ) = ( r 1 ) j 1 2 ,
δ u j δ u j 1 h j 2 ( δ v j + δ v j 1 ) = ( r 5 ) j 1 2 ,
δ θ j δ f θ j 1 h j 2 ( δ p j + δ p j 1 ) = ( r 6 ) j 1 2 ,
δ ϕ j δ ϕ j 1 h j 2 ( δ g j + δ g j 1 ) = ( r 7 ) j 1 2 ,
( ξ 1 ) j δ f j + ( ξ 2 ) j δ f j 1 + ( ξ 3 ) j δ u j + ( ξ 4 ) j δ u j 1 + ( ξ 5 ) j δ v j + ( ξ 6 ) j δ v j 1 + ( ξ 7 ) j δ θ j + ( ξ 8 ) j δ θ j 1 + ( ξ 9 ) j δ ϕ j + ( ξ 10 ) j δ ϕ j 1 = ( r 2 ) j 1 2 ,
( ψ 1 ) j δ f j + ( ψ 2 ) j δ f j 1 + ( ψ 3 ) j δ v j + ( ψ 4 ) j δ v j 1 + ( ψ 5 ) j δ p j + ( ψ 6 ) j δ p j 1 + ( ψ 7 ) j δ g j + ( ψ 8 ) j δ g j 1 = ( r 3 ) j 1 2 ,
( λ 1 ) j δ f j + ( λ 2 ) j δ f j 1 + ( λ 3 ) j δ p j + ( λ 4 ) j δ p j 1 + ( λ 5 ) j δ g j + ( λ 6 ) j δ g j 1 = ( r 4 ) j 1 2 .
Subject to boundary conditions
{ δ f 0 = 0 ,   δ u 0 = 0 , δ p 0 B i 1 δ θ 0 = p 0 + B i 1 ( θ 0 + 1 ) , δ g 0 B i 2 δ ϕ 0 = g 0 + B i 2 ( ϕ 0 + 1 ) , δ u J = 0 ,   δ θ J = 0 ,   δ ϕ J = 0 ,
where
( ξ 1 ) j = ( ξ 2 ) j = h 4 ( v j + v j 1 ) ,
( ξ 3 ) j = ( ξ 4 ) j = h j 2 ( D a + ( u j + u j 1 ) ) ,
( ξ 5 ) j = 1 + γ ( η j + η j 1 ) + h j γ + h j 4 ( f j + f j 1 ) ,
( ξ 6 ) j = 1 γ ( η j + η j 1 ) + h j γ + h j 4 ( f j + f j 1 ) , ( ξ 7 ) j = ( ξ 8 ) j = h j 2 G r C o s ( α ) ,
( ξ 9 ) j = ( ξ 10 ) j = h j 2 G c C o s ( α ) ,
( ψ 1 ) j = ( ψ 2 ) j = P r   h j 4 ( p j + p j 1 ) ,
( ψ 3 ) j = ( ψ 4 ) j = P r E c h j 4 ( 1 + γ ( η j + η j 1 ) ) ( v j + v j 1 ) ,
( ψ 5 ) j = Pr h j 4 ( f j + f j 1 ) + ( 1 + 4 3 R ) ( 1 + γ ( η j + η j 1 ) + h j γ ) ,
( ψ 6 ) j = Pr h j 4 ( f j + f j 1 ) + ( 1 + 4 3 R ) ( 1 γ ( η j + η j 1 ) + h j γ ) ,
( ψ 7 ) j = D u ( 1 + γ ( η j + η j 1 ) + h j γ ) ,
( ψ 8 ) j = D u ( 1 + γ ( η j + η j 1 ) + h j γ ) ,
( λ 1 ) j = ( λ 2 ) j = S c h 4 ( g j + g j 1 ) ,
( λ 3 ) j = S r ( 1 + γ ( η j + η j 1 ) + h j γ ) ,
( λ 4 ) j = S r ( 1 γ ( η j + η j 1 ) + h j γ ) ,
( λ 5 ) j = 1 + γ ( η j + η j 1 ) + h j γ + S c h 4 ( f j + f j 1 ) ,
( λ 6 ) j = 1 γ ( η j + η j 1 ) + h j γ + S c h 4 ( f j + f j 1 ) ,
( r 2 ) j 1 2 = ( 1 + γ ( η j + η j 1 ) ) ( v j v j 1 ) 2 γ h j ( v j 1 2 ) h j ( f j 1 2 ) ( v j 1 2 ) + h j ( u j 1 2 ) 2 + D a   h j ( u j 1 2 ) G r C o s ( α ) h j ( θ j 1 2 ) G c C o s ( α ) h j ( ϕ j 1 2 ) ,
( r 3 ) j 1 2 = ( 1 + 4 3 R ) ( ( 1 + γ ( η j + η j 1 ) ) ( p j p j 1 ) + 2 γ h j ( p j 1 2 ) ) P r h j ( f j 1 2 ) ( p j 1 2 ) D u ( 1 + γ ( η j + η j 1 ) ) ( g j g j 1 ) 2 D u γ h j ( g j 1 2 ) P r E c ( 1 + 2 γ η ) h j ( v j 1 2 ) 2 ,
( r 4 ) j 1 2 = ( 1 + γ ( η j + η j 1 ) ) ( g j g j 1 ) 2 γ h j ( g j 1 / 2 ) S c h j ( f j 1 / 2 ) ( g j 1 / 2 ) S r ( ( 1 + γ ( η j + η j 1 ) ) ( p j p j 1 ) + 2 γ h j ( p j 1 2 ) ) .
In the matrix vector form, we can write
A δ = r ,
in which
A = [ [ A 1 ] [ C 1 ] [ B 2 ] [ A 2 ] [ C 2 ] [ B J 1 ] [ A J 1 ] [ C J 1 ] [ B J ] [ A J ] ]
δ = [ [ δ 1 ] [ δ 2 ] [ δ J 1 ] [ δ J ] ] ,   r = [ [ r 1 ] [ r 2 ] [ r J 1 ] [ r J ] ] .
where in Equation (37) the components are defined by
[ A 1 ] = [ 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 B i 1 1 0 0 0 0 0 0 0 B i 2 1 0 1 C j 0 0 0 0 0 0 0 1 C j 0 0 0 0 0 0 0 1 C j ] , c j = 1 2 h j ,
[ A j ] = [ 1 C j 0 0 0 0 0 ( ξ 1 ) j ( ξ 3 ) j ( ξ 5 ) j ( ξ 7 ) j 0 ( ξ 9 ) j ( ψ 1 ) j 0 ( ψ 3 ) j 0 ( ψ 5 ) j 0 ( ψ 7 ) j ( λ 1 ) j 0 0 0 ( λ 3 ) j 0 ( λ 5 ) j 0 1 C i 0 0 0 0 0 0 0 1 C j 0 0 0 0 0 0 0 1 C j ] , 2 < j < J 1 ,
[ A J ] = [ 1 C J 0 0 0 0 0 ( ξ 1 ) J ( ξ 3 ) J ( ξ 5 ) J ( ξ 7 ) J 0 ( ξ 9 ) J 0 ( ψ 1 ) J 0 ( ψ 3 ) J 0 ( ψ 5 ) J 0 ( ψ 7 ) J ( λ 1 ) J 0 0 0 ( λ 3 ) J 0 ( λ 5 ) J 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 ] ,
[ B j ] = [ 1 C j 0 0 0 0 0 ( ξ 2 ) j ( ξ 4 ) j ( ξ 6 ) j ( ξ 8 ) j 0 ( ξ 10 ) j 0 ( ψ 2 ) j 0 ( ψ 4 ) j 0 ( ψ 6 ) j 0 ( ψ 8 ) j ( λ 2 ) j 0 0 0 ( λ 4 ) j 0 ( λ 6 ) j 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] , 2 < j < J ,
[ C j ] = [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 c j 0 0 0 0 0 0 0 1 c j 0 0 0 0 0 0 0 1 c j ] , 1 < j < J 1 ,
[ δ j ] = [ δ f j δ u j δ v j δ θ j δ p j δ ϕ j δ g j ]   ,   1 < j < J ,
[ r j ] = [ ( r 1 ) j 1 / 2 ( r 2 ) j 1 / 2 ( r 3 ) j 1 / 2 ( r 4 ) j 1 / 2 ( r 5 ) j 1 / 2 ( r 6 ) j 1 / 2 ( r 7 ) j 1 / 2 ]   ,   1 < j < J ,
Now, we let
A = L U
where
L = [ [ a 1 ] [ B 2 ] [ a 2 ] [ B J ] [ a J ] ]   ,
U = [ [ I ] [ Γ 1 ] [ I ] [ Γ 1 ] [ I ] [ Γ J 1 ] [ I ] ]   ,
where [ I ] is the unit matrix while [ a i ] and [ Γ i ] are 7 × 7 matrices whose elements are found by the following expressions:
[ a i ] = [ A 1 ] ,
[ A 1 ] [ Γ 1 ] = [ C 1 ] ,  
[ a j ] = [ A j ] [ B j ] [ Γ j 1 ] ,   j = 2 , 3 , . . , J   ,  
[ a j ] [ Γ j ] = [ C j ] ,   j = 2 , 3 , . . , J 1 .
Equation (38) can be substituted into Equation (37) to get
L U δ = r .
Once the element of δ is found, Equation (38) can be used to find the ( i + 1 ) th iteration. The procedure described above is implemented in Mathematica and results are displayed in the next section.

4. Results Interpretation

To scrutinize the flow features of the second-grade fluid along with the effect of the magnetic field, radiative heat flux, viscous dissipation, and a double diffusion (Dufour and Soret) solution of Equations (9)–(11) is essential. Before getting to the solution of the coupled system, it is highly important to look at the nature and complexity of it. After sighting the system, it is clear that the momentum equation mentioned in Equation (9) is a third order equation, whereas Equations (10)–(11), i.e., the temperature and concentration equations, are second-order equations. So, there are various approaches (finite difference, finite element and finite volume) to get to the solution of these equations. Among them the best is the finite difference approach for such type of problems, so the Keller Box differencing scheme is implemented. Figure 1, Figure 2 and Figure 3 are displayed to portray the behavior of the velocity distribution against the relative parameters. To see the impacts of viscoelastic parameter β on second grade velocity field SVF, Figure 1 is plotted. Variation in velocity against ( β ) is measured for β = 0.0 ,   0.2 ,   0.4 ,   0.6 . An increasing trend in velocity is sketched against an increasing ( β ) . This behavior is justified by the mathematical representation of ( β   = α 1 U o ρ ν l ) that by increasing the magnitude of ( β ) , viscosity decreases as a result of the velocity of the fluid mounts and chaotic behavior uplifts. Here, it is productive to mention that for β = 0.0 , the present problem reduces to the Newtonian case. From a boundary layer point of view, the thickness of the fluid increases with an increase in ( β ) . Figure 2 indicates the incrementing behavior of velocity against inciting the magnitude of the curvature parameter ( γ = 0.0 ,   0.2 ,   0.4 ,   0.6 ,   ). The reason is that by the increasing curvature parameter, bending of the surface as well as the radius decreases. Hence, less friction will be offered to fluid molecules by the surface and velocity uplifts. A declining trend in velocity is depicted against the Darcy parameter (Da) in Figure 3. Since the Darcy parameter (Da) represents the presence of a porous medium, suggests porosity creates a high resistance to fluid molecules and as an outcome velocity declines. Variation in the thermal profile against influencing parameters are disclosed in Figure 4, Figure 5 and Figure 6. Figure 4 is manifested to excogitate the impression of radiation parameter (R) on the thermal field. An increasing pattern in temperature is observed by choosing R = 0.5, 1.0, 1.5, 2.0 and by fixing γ = 0.1 , Pr = 1 ,   D u = 0.5 ,   E c = 0.1 . It is because of the fact that energy flux increases and consequently fluid temperature increases. In Figure 5, the effect of the thermal Biot number ( B i 1   ) on θ(η) is presented. By definition, the thermal Biot number is directly related to the heat transfer coefficient generated by the hot fluid. Thus, as the thermal Biot number increases, the convection due to the hot fluid raises and the temperature mounts. Impacts of the Dufour (Du) and Soret (Sr) aspects on the thermal field is disclosed in Figure 6. In the present graph, the values of (Du) and (Sr) are selected in such a way that their product will give a constant magnitude. By growing the Dufour number (i.e., by declining the Soret number) the thermal change between the hot and ambient fluid increases, which enhances the temperature. Figure 7, Figure 8 and Figure 9 are plotted to predict the changes in the concentration profile with respect to the involved parameters like γ ,     B i 2   ,   D u   a n d   S r . Figure 7 presents the impression of the curvature parameter γ on the concentration profile. Duality in concentration features are interpreted against the curvature parameter. For small values of η, i.e., (η < 2), the concentration decreases whereas it increases when η > 3. The upshot of the concentration Biot number ( B i 2 ) on ϕ(η) is captured in Figure 8. From the drawn curves it is portrayed that by increasing the ( B i 2 ) concentration field the associated boundary layer thickness increases. The joint conspiration of (Sr) and (Du) on the concentration field is divulged in Figure 9. It is interpreted that by decreasing (Sr) and increasing (Du) the strength of the intermolecular forces weakens and as a consequence the concentration field declines. Table 1, Table 2 and Table 3 are enumerated to record the variation in the Nusselt number ( N u ) , Sherwood number ( S h ) and coefficient of skin-friction ( C f ) for numerous values of involved parameters. Excellent correlation is noticed between our results and findings made by Hayat et al. [31]. This comparison between attained figures assures the credibility of the current work.

5. Outcomes of Analysis

The current communication is devoted to explicate the flow features of magnetohydrodynamic boundary layer flow of a second-grade fluid over an inclined porous cylinder under the impacts of viscos dissipation, radiation, convective heating, and Dufour and Soret effects. The mathematical modelling is attained in the form of partial differential equations and converted into ordinary differential equations by employing compatible transformations. Afterwards the solution is accomplished by using an implicit finite difference method known as the Keller Box scheme. The impact of the involved variables on velocity, temperature and concentration profiles is excogitated through graphical visualizations. Assurance of the present findings is achieved by constructing comparative analysis with previously published literature. The key findings are summarized as follows:
(i)
Increasing aptitude of velocity within the boundary layer region is depicted against the curvature parameter;
(ii)
It is found that the velocity profile upsurges against the viscoelastic parameter whereas it declines by uplifting the Darcy parameter;
(iii)
The thermal Biot number raises the temperature profile and also enriches the magnitude of the Nusselt number;
(iv)
Concentration of the Biot number causes growth in the concentration profile and augments the mass flux coefficient;
(v)
Dufour and Soret effects enhance the temperature field and depreciate the concentration profile.

Author Contributions

Formulation done by S.B. and A.H.M.; problem solved by R.M. and S.B.; results computed by I.K. and A.H.S.; results discussed by E.-S.M.S.; all authors contributed equally in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia.

Acknowledgments

Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u , v , w Dimensional velocity components f Dimensionless velocity component
TDimensional temperatureθDimensionless temperature
CDimensional concentration ϕ Dimensionless concentration
ρ Density D m Mass diffusivity
γ Curvature parameter D a Darcy parameter
D u Dufour number S r Soret number
G r Thermal Grashof numberGcSolutal Grashof number
B i 1 Thermal Biot number R Radiation parameter
η Dimensionless similarity variable B i 2 Concentration Biot number
N u x Nusselt number C f x Skin friction coefficient
S h x Sherwood number K Dimensional curvature parameter
β Viscoelastic parameter T f Fluid temperature
T Ambient fluid temperature C Ambient fluid concentration

References

  1. Vejravelu, K.; Roper, T. Flow and heat transfer in a second-grade fluid over a stretching sheet. Int. J. Non Linear Mech. 1999, 34, 1031–1036. [Google Scholar] [CrossRef]
  2. Rajeswari, G.K.; Rathna, S.L. Flow of a Particular Class of non-Newtonian Visco-elastic and Visco-elastic Fluids near a Stagnation Point. Z. Angew. Math. Phys. 1962, 13, 43–57. [Google Scholar] [CrossRef]
  3. Garg, V.K.; Rajagopal, K.R. Flow of a non-Newtonian Fluid Past a Wedge. Acta Mech. 1991, 88, 113–123. [Google Scholar] [CrossRef]
  4. Fetecau, C. Starting Solutions for the Motion of a Second Grade Fluid Due to Longitudnal and Torsional Oscillations of a Circular Cylinder. Int. J. Eng. Sci. 2006, 44, 788–796. [Google Scholar] [CrossRef]
  5. Bilal, S.; Mustafa, Z.; Rehman, K.U.; Malik, M.M. MHD Second Grade NanoFluid Flow Induced by a Rotatory Cone. J. Nanofluids 2019, 8, 876–884. [Google Scholar] [CrossRef]
  6. Abel, M.S.; Mahesha, N.; Malipatil, S.B. Heat Transfer Due to MHD Slip Flow of a Second-Grade Liquid Over a Stretching Sheet Through a Porous Medium with Non Uniform Heat Source/Sink. Chem. Eng. Commun. 2011, 198, 191–213. [Google Scholar] [CrossRef]
  7. Mushtaq, M.; Asghar, S.; Hossain, M.A. Mixed Convection Flow of Second Grade Fluid along a Vertical Stretching Flat Surface with Variable Surface Temperature. Heat Mass Transf. 2007, 43, 1049–1061. [Google Scholar] [CrossRef]
  8. Cortell, R. MHD Flow and Mass Transfer of an Electrically Conducting Fluid of Second Grade in a Porous Medium over a Stretching Sheet with Chemically Reactive Species. Chem. Eng. Process. 2007, 46, 721–728. [Google Scholar] [CrossRef]
  9. Merkin, J.H. Natural Convection Boundary-Layer Flow on a Vertical Surface with Newtonian Heating. Int. J. Heat Fluid Flow 1994, 15, 392–398. [Google Scholar] [CrossRef]
  10. Salleh, M.Z.; Nazar, R.; Pop, I. Boundary layer flow and heat transfer over a stretching sheet with Newtonian heating. J. Taiwan Inst. Chem. Eng. 2010, 41, 651–655. [Google Scholar] [CrossRef]
  11. Haq, R.U.; Nadeem, S.; Khan, Z.; Okedayo, T. Convective heat transfer and MHD effects on Casson nanofluid flow over a shrinking sheet. Cent. Eur. J. Phys. 2014, 129, 862–871. [Google Scholar] [CrossRef]
  12. Nadeem, S.; Haq, R.U.; Akbar, N.S. MHD three-dimensional boundary layer flow of Casson nanofluid past a linearly stretching sheet with convective boundary condition. IEEE Trans. Nanotechnol. 2014, 13, 1326–1332. [Google Scholar] [CrossRef]
  13. Kameswaran, P.K.; Sibanda, P.; Murti, A.S.N. Nanofluid flow over a permeable surface with convective boundary conditions and radiative heat transfer. Math. Probl. Eng. 2013, 201219. [Google Scholar] [CrossRef]
  14. Bakar, N.A.A.; Zaimi, W.M.K.A.W.; Hamid, R.A.; Bidin, B.; Ishak, A. Boundary layer flow over a stretching sheet with a convective boundary condition and slip effect. World Appl. Sci. J. 2012, 17, 49–53. [Google Scholar]
  15. Hayat, T.; Shehzad, S.A.; Qasim, M.; Obaidat, S. Flow of a second-grade fluid with convective boundary conditions. Int. J. Therm. Sci. 2011, 15, S253–S261. [Google Scholar] [CrossRef]
  16. Makinde, O.D.; Aziz, A. Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition. Int. J. Therm. Sci. 2011, 50, 1326–1332. [Google Scholar] [CrossRef]
  17. Merkin, J.H.; Pop, I. The forced convection flow of a uniform stream over a flat surface with a convective surface boundary condition, Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 3602–3609. [Google Scholar] [CrossRef]
  18. Kothandapani, M.; Prakash, J. Effects of thermal radiation parameter and magnetic field on the peristaltic motion of Williamson nanofluid in a tapered asymmetric channel. Int. J. Heat Mass Tran. 2015, 81, 234–245. [Google Scholar] [CrossRef]
  19. Kumar, K.G.; Rudraswamy, N.G.; Gireesha, B.J.; Manjunatha, S. Nonlinear thermal radiation effect on Williamson fluid with particle-liquid suspension over a stretching surface. Results Phys. 2017, 7, 3196–3202. [Google Scholar] [CrossRef]
  20. Khan, M.; Hussain, M.; Azam, M. Magnetohydrodynamic flow of Carreau fluid over a convectively heated surface in the presence of thermal radiation. J. Magn. Magn. Mater. 2016, 412, 63–68. [Google Scholar] [CrossRef]
  21. Waqas, M.; Khan, M.I.; Hayat, T.; Alsaedi, A. Numerical simulation for magneto Carreau nanofluid model with thermal radiation: A revised model. Comput. Method. Appl. M. 2017, 324, 640–653. [Google Scholar] [CrossRef]
  22. Goodarzi, M.; Safaei, M.R.; Oztop, H.K.; Karimipour, A.; Sadighinezhad, E.; Dahari, M.; Kazi, S.N.; Jomhari, N. Numerical study of entropy generation due to coupled laminar and turbulent mixed convection and radiation in an enclosure filled with a semitransparent medium. Sci. World J. 2014. [Google Scholar] [CrossRef] [PubMed]
  23. Podarjani, H.; Aghakhani, S.; Karimpour, A.; Afrand, M.; Goodarzi, M. Investigation of free convection heat transfer and entropy generation of nanofluid flow inside a cavity affected by magnetic field by magnetic field and radiation. J. Therm. Anal. Calorim. 2019, 37, 997–1019. [Google Scholar]
  24. Aghaei, A.; Sheikhzadeh, G.A.; Goodarzi, M.; Hasani, H.; Damirchi, H.; Afrand, M. Effect of horizontal and vertical elliptic baffles inside an enclosure on the mixed convection of a MWCNTs-water nanofluid and its entropy generation. Eur. Phys. J. Plus 2018, 133, 486–497. [Google Scholar] [CrossRef]
  25. Hirshfelder, J.O.; Bird, R.C. Molecular Theory of Gases and Liquids; John Wiley: New York, NY, USA, 1954. [Google Scholar]
  26. Afify, A.A. Similarity solution in MHD effects of thermal diffusion and diffusionthermo on free convective heat and mass transfer over a stretching surfaceconsidering suction or injection. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 2202–2214. [Google Scholar] [CrossRef]
  27. Bhaatacharyya, K.; Layek, G.C.; Seth, G.S. Soret and Dufour effects on convective heat and mass transfer in stagnation point flow towards a shrinking sheet. Phys. Scr. 2014, 89, 095203. [Google Scholar] [CrossRef]
  28. Awad, F.G.; Sibanda, P.; Khidir, A.A. Thermo diffusion effects on magneto-nano fluid flow over a stretching sheet. Bound. Value Prob. 2013. [Google Scholar] [CrossRef] [Green Version]
  29. Goyal, M.; Bhrgava, R. Numerical study of thermodiffusion effects on boundary layer flow of nanofluids over a power law stretching sheet. Inst. Mech. Eng. 2014, 17, 591–604. [Google Scholar] [CrossRef]
  30. Khan, U.; Mohyuddin, S.T.; Mohsin, B.B. Convective heat transfer and thermo-diffusion effects on flow of nanofluid towards a permeable stretching sheetsaturated by a porous medium. Aerosp. Sci. Technol. 2015, 50, 196–203. [Google Scholar] [CrossRef]
  31. Hayat, T.; Saeed, Y.; Asad, S.; Alsaedi, A. Convective heat and mass transfer in flow by an inclined stretching cylinder. J. Mol. Liq. 2016, 220, 573–580. [Google Scholar] [CrossRef]
  32. Bagherzadeh, S.A.; Jalali, E.; Sarafraz, M.M.; Akbari, O.A.; Karimipour, A.; Goodarzi, M.; Bach, Q.V. Effects of magnetic field on micro cross jet injection of dispersed nanoparticles in a microchannel. Int. J. Numer. Meth. Fluid Flow 2019. [Google Scholar] [CrossRef]
  33. Goshayeshi, H.M.; Goodarzi, M.; Safaei, M.R.; Dahari, M. Experimental study on the effect of inclination angle on heat transfer enhancement of a ferrofluid in a closed loop oscillating heat pipe under magnetic field. Exper. Therm. Fluid Sci. 2016, 74, 265–270. [Google Scholar] [CrossRef]
  34. Goshayeshi, H.R.; Goodarzi, M.; Dahari, M. Effect of magnetic field on the heat transfer rate of kerosene/Fe2O3 nanofluid in a copper oscillating heat pipe. Exper. Therm. Fluid Sci. 2015, 68, 663–668. [Google Scholar] [CrossRef]
  35. Yousefzadeh, S.; Rajabi, H.; Ghajari, N.; Sarafraz, M.M.; Akbari, O.A.; Goodarz, M. Numerical investigation of mixed convection heat transfer behavior of nanofluid in a cavity with different heat transfer areas. J. Therm. Analy. Calorimetry 2019, 1–25. [Google Scholar] [CrossRef]
  36. Tian, Z.; Etedali, S.; Afrand, M.; Abdollahi, A.; Goodarzi, M. Experimental study of the effect of various surfactants on surface sediment and pool boiling heat transfer coefficient of silica/DI water nano-fluid. Powder Technol. 2019, 356, 391–402. [Google Scholar] [CrossRef]
  37. Bahmani, M.H.; Akbari, O.A.; Zarringhalam, M.; Shabani, G.A.S.; Goodarzi, M. Forced convection in a double tube heat exchanger using nanofluids with constant and variable thermophysical properties. Int. J. Numer. Methd. Fluid Flow 2019. [Google Scholar] [CrossRef]
Figure 1. f ( η ) for different β .
Figure 1. f ( η ) for different β .
Energies 13 00278 g001
Figure 2. f ( η ) for different γ .
Figure 2. f ( η ) for different γ .
Energies 13 00278 g002
Figure 3. f ( η ) for different D a .
Figure 3. f ( η ) for different D a .
Energies 13 00278 g003
Figure 4. θ ( η ) for different R .
Figure 4. θ ( η ) for different R .
Energies 13 00278 g004
Figure 5. θ ( η ) for different β i 1 .
Figure 5. θ ( η ) for different β i 1 .
Energies 13 00278 g005
Figure 6. θ ( η ) for different S r and D u .
Figure 6. θ ( η ) for different S r and D u .
Energies 13 00278 g006
Figure 7. ϕ ( η ) for different γ .
Figure 7. ϕ ( η ) for different γ .
Energies 13 00278 g007
Figure 8. ϕ ( η ) for different β i 2 .
Figure 8. ϕ ( η ) for different β i 2 .
Energies 13 00278 g008
Figure 9. ϕ ( η ) for different S r and D u .
Figure 9. ϕ ( η ) for different S r and D u .
Energies 13 00278 g009
Table 1. Variation in the local Nusselt number N u x / R e x 1 / 2 for various parameters.
Table 1. Variation in the local Nusselt number N u x / R e x 1 / 2 for various parameters.
γ G r G c α D a P r Hayat et al. [31]Present Results
0.10.20.2 π / 4 0.210.14270.142683
0.0 0.13970.139746
0.1 0.14270.142683
0.12 0.14330.143252
0.1 0.14170.141627
0.3 0.14360.143621
0.5 0.14520.145236
0.2 0.14270.142683
0.4 0.14390.143921
0.6 0.14500.145011
π / 6 0.14340.143397
π / 4 0.14270.142683
π / 3 0.14170.141653
0.1 0.14450.144551
0.3 0.14090.140885
0.5 0.13760.137455
0.90.13940.139384
1.00.14270.142683
1.20.14800.148027
Table 2. Variation in the local Sherwood number S h x / R e x 1 / 2 for various parameters.
Table 2. Variation in the local Sherwood number S h x / R e x 1 / 2 for various parameters.
R B i 1 B i 2 D u E c S c Hayat et al. [31]Present Results
0.10.20.20.50.10.90.14270.142683
0.3 0.17120.171117
0.5 0.19830.198211
0.7 0.20490.224187
0.3 0.19000.189984
0.4 0.22780.227806
0.5 0.25880.258763
0.3 0.13910.139068
0.4 0.13620.136198
0.5 0.13390.133862
0.3 0.14750.147474
0.7 0.13790.137879
1.1 0.12820.128200
0.2 0.13100.130941
0.3 0.11940.119422
0.4 0.10810.108115
0.70.14420.144269
1.20.14110.141060
1.60.13970.139623
Table 3. Variation in the skin-friction coefficient C f x / R e x 1 / 2 for various parameters.
Table 3. Variation in the skin-friction coefficient C f x / R e x 1 / 2 for various parameters.
γ P r R B i 1 B i 2 D u E c S c S r Hayat et al. [31]Present Results
0.11.00.10.20.20.50.10.90.20.14450.144411
0.0 0.14210.142061
0.1 0.14450.144411
0.12 0.14490.144864
0.9 0.14470.144664
1.0 0.14450.144411
1.2 0.14410.144013
0.3 0.14480.144730
0.5 0.14510.145036
0.7 0.14540.145323
0.3 0.14320.143138
0.4 0.14220.142127
0.5 0.14130.141304
0.3 0.19240.192309
0.4 0.23070.230657
0.5 0.26210.262084
0.3 0.14410.144029
0.6 0.14460.144601
0.9 0.14520.145170
0.3 0.14650.146412
0.4 0.14750.147384
0.5 0.14840.148338
0.9 0.14450.144411
1.2 0.15200.151954
1.6 0.15860.158605
0.40.14110.141043
0.50.13950.139376
0.60.13780.137719

Share and Cite

MDPI and ACS Style

Bilal, S.; Majeed, A.H.; Mahmood, R.; Khan, I.; Seikh, A.H.; Sherif, E.-S.M. Heat and Mass Transfer in Hydromagnetic Second-Grade Fluid Past a Porous Inclined Cylinder under the Effects of Thermal Dissipation, Diffusion and Radiative Heat Flux. Energies 2020, 13, 278. https://doi.org/10.3390/en13010278

AMA Style

Bilal S, Majeed AH, Mahmood R, Khan I, Seikh AH, Sherif E-SM. Heat and Mass Transfer in Hydromagnetic Second-Grade Fluid Past a Porous Inclined Cylinder under the Effects of Thermal Dissipation, Diffusion and Radiative Heat Flux. Energies. 2020; 13(1):278. https://doi.org/10.3390/en13010278

Chicago/Turabian Style

Bilal, Sardar, Afraz Hussain Majeed, Rashid Mahmood, Ilyas Khan, Asiful H. Seikh, and El-Sayed M. Sherif. 2020. "Heat and Mass Transfer in Hydromagnetic Second-Grade Fluid Past a Porous Inclined Cylinder under the Effects of Thermal Dissipation, Diffusion and Radiative Heat Flux" Energies 13, no. 1: 278. https://doi.org/10.3390/en13010278

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop