Open access peer-reviewed chapter

Darcy-Forchheimer Flow of Casson Nanofluid with Heat Source/Sink: A Three-Dimensional Study

By Gosikere Kenchappa Ramesh

Submitted: November 2nd 2017Reviewed: January 18th 2018Published: September 19th 2018

DOI: 10.5772/intechopen.74170

Downloaded: 642


In this chapter, three-dimensional Casson nanoliquid flow in two lateral directions past a porous space by Darcy-Forchheimer articulation is examined. The study includes the impact of uniform heat source/sink and convective boundary condition. The administering partial differential equations are shaped to utilizing comparability changes into a set of nonlinear normal differential conditions which are fathomed numerically. The self-comparative arrangements are gotten and contrasted and accessible information for uncommon cases. The conduct of parameters is displayed graphically and examined for velocity, temperature, and nanoparticle volume part. It is discovered that temperature and nanoparticle volume fraction rise for enhancement of Forchheimer and porosity parameters.


  • three-dimensional flow
  • Darcy-Forchheimer porous medium
  • Casson nanoliquid
  • uniform heat source/sink
  • convective boundary condition
  • numerical solutions

1. Introduction

In many assembling processes and for mechanical reason, the investigation of heat exchange and boundary layer flow over linearly and nonlinearly extending surface are much imperative. These procedures and applications incorporate streamlined feature forming, wire drawing, and paper generation where a specific temperature will be required for cooling the particles in the liquid. At first, the stream qualities have been analyzed by [1] overextending surfaces. The perfection of finishing up item relies upon the rate of warmth exchange at the surface of extending material. Many creators expanded crafted by [1] managed heat exchange qualities alongside the flow conduct in different physical circumstances in [2, 3, 4, 5, 6, 7, 8].

Non-Newtonian fluids can’t be portrayed because of nonexistence of single constitutive connection among stress and rate of strain. In the current year, non-Newtonian fluids have turned out to be increasingly essential because of its mechanical applications. Truth be told, the enthusiasm for boundary layer flows of non-Newtonian fluid is expanding significantly because of its extensive number of functional applications in industry producing preparing and natural fluids. Maybe, a couple of principle illustrations identified with applications are plastic polymer, boring mud, optical fibers, paper generation, hot moving, metal turning, and cooling of metallic plates in a cooling shower and numerous others. Since no single non-Newtonian model predicts every one of the properties of non-Newtonian fluid along these lines examinations proposed different non-Newtonian fluid models. These models are essentially classified into three classifications specifically differential-, rate-, and fundamental-type fluids. In non-Newtonian fluid, shear stresses and rates of strain/disfigurement are not directly related. Such fluid underthought which does not comply with Newton’s law is a straightforward non-Newtonian fluid model of respectful sort. In 1959, Casson displayed this model for the flow of viscoelastic fluids. This model has a more slow progress from Newtonian to the yield locale. This model is utilized by oil builds in the portrayal of bond slurry and is better to predict high shear-rate viscosities when just low and middle road shear-rate information are accessible. The Casson show is more exact at both high and low shear rates. Casson liquid has one of the kind attributes, which have wide application in sustenance handling, in metallurgy, in penetrating operation and bio-designing operations, and so on. The Casson show has been utilized as a part of different businesses to give more exact portrayal of high shear-rate viscosities when just low and moderate shear-rate information are accessible [9]. Toward the starting Nadeem et al. [10] introduce the idea of Casson fluid and demonstrate over an exponentially extending sheet. Numerous examinations identified with viscoelastic properties of liquid are underthought [11, 12, 13, 14, 15, 16, 17].

The nanoparticles in by and large are made of metal oxides, metallic, carbon, or some other materials [18]. Standard fluid has weaker conductivity. This weaker conductivity can be enhanced incredibly with the utilization of nanoparticles. Truly, the Brownian movement factor of nanoparticles is base fluid and is essential toward this path. An extraordinary measure of warmth is delivered in warm exchangers and microelectromechanical procedures to lessen the framework execution. Fluid thermal conductivity is enhanced by nanoparticle expansion just to cool such modern procedures. The nanoparticles have shallow significance in natural and building applications like prescription, turbine sharp-edge cooling, plasma- and laser-cutting procedure, and so on. Sizeable examinations on nanofluids have been tended to in the writing. Buongiorno [19] has investigated the mechanisms of nanofluid by means of snapshot of nanoparticles in customary base fluid. Such instruments incorporate nanoparticle measure, Magnus effect, dormancy, molecule agglomeration, Brownian movement, thermophoresis, and volume portion. Here, we introduce some imperative scientists who have been accounted for by considering the highlights of thermophoretic and Brownian movement [20, 21, 22, 23, 24, 25, 26, 27].

In displaying the flow in permeable media, Darcy’s law is a standout among the most prominent models. Particularly, flow in permeable media is exceptionally valuable in grain stockpiling, development of water in repositories, frameworks of groundwater, fermentation process, raw petroleum generation, and oil assets. In any case, it is by and large perceived that Darcy’s model may over anticipate the convective streams when the inertial drag and vorticity dissemination coproductive are considered. The expansion of established Darcy demonstrates incorporates inertial drag and vorticity dispersion. To think about the inertial drag and vorticity dispersion, Forchheimer [28] joined the root mean square. Forchheimer’s term was named by Muskat [29] and inferred that the consideration of Forchheimer’s term is substantial for high Reynolds number. Pal and Mondal [30] examined the hydromagnetic Darcy-Forchheimer flow for variable liquid property. Utilizing HAM strategy, Hayat et al. [31] got the systematic answer for Darcy-Forchheimer stream of Maxwell liquid by considering the Cattaneo-Christov hypothesis. Vishnu Ganesh et al. [32] analyzed the thick and Ohmic dispersals, and the second-order slip consequences for Darcy-Forchheimer flow of nanoliquid past an extending/contracting surface. Scientific model for Darcy-Forchheimer stream of Maxwell fluid with attractive field and convective boundary condition is given by Sadiq and Hayat [33]. Utilizing Keller’s box strategy, Ishak et al. [34] numerically examined the magnetohydrodynamic flow and heat exchange exhibitions over an extending cylinder. Mixed convective flow and the related warmth and mass exchange qualities over a vertical sheet immersed in a permeable medium have been explored by Pal and Mondal [35] by considering different impacts, for example, Soret, Dufour, and thermal radiation.

The principle objective here is to uncover the Darcy-Forchheimer connection on a three-dimensional Casson nanofluid flow past a stretching sheet. Heat transfer is handled with regular heat generation/absorption and convective-type boundary condition.


2. Mathematical formulation

Three-dimensional flow of Casson nanofluid filling permeable space by Darcy-Forchheimer connection is considered. Flow is bidirectional extending surface. Nanofluid model depicts the properties of Brownian dispersion and thermophoresis. Concurrent states of heat convective and heat source/sink are executed. We receive the Cartesian coordinate in such a way to the point that and pivot are picked along xand yordinary to the stretchable surface. Let Uwx=axand Vwy=bybe the extending velocity along xand ydirections separately. The surface temperature is directed by a convective heating procedure which is depicted by heat exchange coefficient hfand temperature of hot liquid Tfunder the surface (see Figure 1). The boundary layer equations for flow underthought are


and boundary conditions of the problem is


Figure 1.

Geometry of the problem.

Here u,vand wrepresent as components of velocity in x,yand zdirections, respectively; ν=μρfstands for kinematic viscosity; μfor dynamic viscosity; ρffor density of base liquid; Kfor permeability of porous medium; F=CbxK12for nonuniform inertia coefficient of porous medium; Cbfor drag coefficient; Tfor temperature; αm=Kρcffor thermal diffusivity; q0volumetric rate of heat generation and absorption; βis the Casson parameter; kfor thermal conductivity; ρcpfor effective heat capacity of nanoparticles; ρcffor heat capacity of fluid; Cfor concentration; DBfor Brownian diffusion coefficient; DTfor thermophoretic diffusion coefficient; Tfor ambient fluid temperature; Cfor ambient fluid concentration; and aand bfor positive constants.

Selecting similarity transformations are

u=axf 'ς,v=ayg'ς,w=aν12f+g,

Applying Eq. (7) in (1) is verified. and Eqs. (2)(5) are


Boundary conditions of Eq. (6) become

f0=g0=0,f '0=1,g'0=α,θ'0=Bi1θ0,ϕ0=1,
f '0,g'0,θ0,φ0E12

In the above expressions, λstands for porosity parameter, Frfor Forchheimer parameter, αfor ratio parameter, Prfor Prandtl number, Lefor Schmidt number, Bifor Biot number, Ntfor thermophoresis parameter, Nbfor Brownian motion parameter, and Qheat source/sink parameter.

These dimensionless variables are given by


Dimensionless relations of skin friction coefficient, local Nusselt number, and local Sherwood number are as follows:


where Rex=Uwxνand Rey=Vwyνdepict the local Reynolds numbers.


3. Results and discussion

The correct arrangements do not appear to be achievable for an entire arrangement of Eqs. (8)(11) with proper limit conditions given in Eq. (12) in light of the nonlinear shape. This reality compels one to get the arrangement of the issue numerically. Suitable likeness change is received to change the overseeing incomplete differential conditions into an arrangement of non-straight customary differential conditions. The resultant limit esteem issue is understood numerically utilizing an efficient fourth-order Runge-Kutta method alongside shooting method (see Ramesh and Gireesha [27]). Facilitate the union examination is available in Table 1. For the verification of precision of the connected numerical plan, an examination of the present outcomes compared to the 1+1βf0and 1+1βg0(nonappearance of Forchheimer parameter, porosity parameter) with the accessible distributed consequences of Ahmad and Nazar [16] and Nadeem et al. [15] is made and exhibited in Table 2, demonstrates a great understanding in this manner give confidence that the numerical outcomes got are precise.


Table 1.

Convergence analysis of the present work.

1+1βf01+1βf ″01+1βg01+1βg0
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]
[15, 16]

Table 2.

Current numerical values and validation for friction factor 1+1βf0.

This section is fundamentally arranged to depict the effect of different correlated physical parameters on velocity profiles f ′ς,gς, temperature profile θς, nanoparticle part ϕς, skin friction coefficient, and the local Nusselt and local Sherwood number through Figures 216. Give the first focus on the impact of extending parameter αon velocity profile as shown in Figure 2. It is noticed that for expanding benefits of extending parameter α, it decreases the speed f ′ς, while gςfluctuates for different benefits of extending parameter α. It see that for α=0, exhibit wonders lessen the instance of two dimensional linear stretching, while for ​α = 1​, sheet will extended along the two bearings with a similar proportion (axisymmetric case), and third and last case identify with extending proportion parameter αother than 0 or 1; at that point, the flow conduct along both directions will be extraordinary.

Figure 2.

Influence ofαon velocity field.

Attributes of Forchheimer parameter Fron f ′ς,gςis plotted in Figure 3. Obviously, f ′ς,gςis a diminishing capacity of Fr. Figure 4 features the impact of Fron θς,ϕς. Of course θς,ϕςand related thickness of boundary layer are higher when Frincrements. Figure 5 illustrates the varieties in f ′ς,gςfor particular estimations of λ. Both f ′ς,gςand related layer thickness decay for larger λ. Physically, nearness of permeable media is to upgrade the protection from smooth movement which makes decay in fluid speed and thickness of energy boundary layer. Figure 6 depicts the impact of λon θς,ϕς. It is discovered that bigger λcauses an addition in θς,ϕς.

Figure 3.

Influence ofFron velocity fields.

Figure 4.

Influence ofFron temperature and concentration fields.

Figure 5.

Influence ofλon velocity field.

Figure 6.

Influence ofλon temperature and concentration fields.

Figure 7 showed the impacts of non-Newtonian Casson fluid parameter βon velocity profiles f ′ς,gς. It is seen that with the influence of βinfers an abatement in the yield worry of the Casson fluid. This successfully encourages flow of the fluid, i.e., quickens the boundary layer flow near the extending surface, as appeared in Figure 7. In addition, it is found that with substantial estimations of β, the fluid is nearer to the Newtonian fluid. Truth be told, expanding estimations of the Casson fluid parameter βupgrade both temperature and nanoparticle division which is shown in Figure 8.

Figure 7.

Influence ofβon velocity fields.

Figure 8.

Influence ofβon temperature and concentration fields.

The variety in dimensionless temperature profile θςwith expanding estimations of generalized Prandtl number Pris appeared through Figure 9. The temperature profile diminishes with an expansion in the estimations of Prandtl number Pr, as Prandtl number is the proportion of energy diffusivity to thermal diffusivity. So, an expanding estimation of Prandtl number Prinfers a moderate rate of thermal dissemination which thus lessens the thermal boundary layer thickness. It can be directly seen that Prandtl number has more noticeable impact on Newtonian liquid when contrasted with non-Newtonian liquid.

Figure 9.

Influence ofPron temperature fields.

Figure 10 displays the temperature θςand the nanoparticle division ϕςfor variable estimations of Brownian movement Nb. The fluid velocity is found to increment with expanding Nb, while in nanoparticle fraction decreases as Nbexpansion which consequently improves the nanoparticle’s concentration at the sheet. This might be because of the way the Brownian movement parameter diminishes the mass exchange of a nanofluid. The diagram of thermophoresis parameter Nton the temperature θςand the nanoparticle part ϕςprofiles is portrayed in Figure 11. From these plots, it is seen that the impact of expanding estimations of Ntis to build the temperature and nanoparticle fraction.

Figure 10.

Influence ofNbon temperature and concentration fields.

Figure 11.

Influence ofNton temperature and concentration fields.

Figure 12 shows the impact of Lewis number Leon temperature θςand the nanoparticle portion ϕςprofiles. It is take note of that the temperature of the liquid increments however nanoparticle portion of the fluid diminishes with increment in Le. Physically truth that the bigger estimations of Lewis number makes the mass diffusivity littler, subsequently it diminishes the fixation field. The impacts of heat source/sink parameter Qcan be found in Figure 13. For Q>0(heat source), it can be seen that the thermal boundary layer produces the vitality, and this causes the temperature in the thermal boundary layer increments with increment in Q. Though Q<0(heat sink) prompts diminish in the thermal boundary layer. Q=0speaks to the nonattendance of heat source/sink.

Figure 12.

Influence ofLeon temperature and concentration fields.

Figure 13.

Influence ofQon temperature field.

Impacts of the Biot number Bion temperature are shown in Figure 14. Physically, the Biot number is communicated as the convection at the surface of the body to the conduction inside the surface of the body. At the point when thermal angle is connected at first glance, the proportion representing the temperature inside a body fluctuates significantly, while the body heats or cools over a period. Regularly, for uniform temperature field inside the surface, we consider Bi<<1. In any case, Bi>>1portrays that temperature field inside the surface is not uniform. In Figure 14, we have examined the impacts of Biot number Bion the temperature and nanoparticle portion profiles in two ways. The first one is the situation.

Figure 14.

Influence ofBion temperature and concentration fields.

when Bi<1. It is seen from Figure 14 that for the littlest estimations of the Biot number Bi<1, the variety of temperature inside the body is slight and can sensibly be approximated as being uniform. While in the second case, Bi>1delineates that the temperature inside the body is not performing a uniform conduct (see Figure 14).

The impact of physical parameter on nearby Nusselt θ0and Sherwood number ϕ0is displayed in Figure 15. We can see through Figure 15 that non-Darcy-Forchheimer connection creates the low heat and mass at the divider when contrasted with the Darcy-Forchheimer connection. Thus, it is seen with an expansion of the two reasons for speeding up in the λandFr. From Figure 16, the expanding estimations of the heat source/sink parameter Qimprove the local Nusselt number θ0and decrease the local Sherwood number ϕ0with Bi. A similar conduct is likewise found for the variety in NtandNbwhich is portrayed in Figure 17.

Figure 15.

Variation ofθ0andϕ0withβfor different values ofFr,λ.

Figure 16.

Variation ofθ0andϕ0withBifor different values ofQ.

Figure 17.

Variation ofθ0andϕ0withNtfor different values ofNb.


4. Conclusions

Three-dimensional flow of Casson nanoliquid within the sight of Darcy-Forchheimer connection, uniform warmth source/sink, and convective type boundary condition is considered. Numerical plan prompts the arrangements of physical marvel. From this investigation, we analyzed that the expanding Casson parameter compares to bring down velocity and higher temperature fields. The nearness of Frand λcaused a lessening in velocity and increasing speed on temperature and nanoparticle portion. The bigger Biot number improved the temperature and nanoparticle division. Additionally, for vast estimations of Biot number, there are no noteworthy changes in θ0and ϕ0, which are available in Table 3. In heat exchange issues, heat sink parameter controls the relative thickening of the force and the thermal boundary layers.


Table 3.

Computational values of local Nusselt number and local Sherwood number for several values of Bi.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Gosikere Kenchappa Ramesh (September 19th 2018). Darcy-Forchheimer Flow of Casson Nanofluid with Heat Source/Sink: A Three-Dimensional Study, Heat and Mass Transfer - Advances in Modelling and Experimental Study for Industrial Applications, Yong Ren, IntechOpen, DOI: 10.5772/intechopen.74170. Available from:

chapter statistics

642total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Modelling of Bubbly Flow in Bubble Column Reactors with an Improved Breakup Kernel Accounting for Bubble Shape Variations

By Weibin Shi, Jie Yang, Guang Li, Yuan Zong and Xiaogang Yang

Related Book

Advances in Microfluidic Technologies for Energy and Environmental Applications

Edited by Yong Ren

First chapter

Droplet Microfluidic Device for Rapid and Efficient Metals Separation Using Host-Guest Chemistry

By Yehezkiel Steven Kurniawan, Ramachandra Rao Sathuluri and Keisuke Ohto

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us