Numerical simulation of von Karman swirling bioconvection nanofluid flow from a deformable rotating disk Kadir, A, Mishra, SR, Shamshuddin, M and Beg, OA Title Numerical simulation of von Karman swirling bioconvection nanofluid flow from a deformable rotating disk Authors Kadir, A, Mishra, SR, Shamshuddin, M and Beg, OA Type Conference or Workshop Item URL This version is available at: http://usir.salford.ac.uk/id/eprint/47891/ Published Date 2018 USIR is a digital collection of the research output of the University of Salford. Where copyright permits, full text material held in the repository is made freely available online and can be read, downloaded and copied for non- commercial private study or research purposes. Please check the manuscript for any further copyright restrictions. For more information, including our policy and submission procedure, please contact the Repository Team at: [email protected].
2
Embed
Numerical simulation of von Karman swirling bioconvection ...usir.salford.ac.uk/47891/1/ICHTFM 2018 poster BIOCONVECTION NA… · NUMERICAL SIMULATION OF VON KARMAN SWIRLING BIOCONVECTION
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
N u m e ric al si m ula tion of von Kar m a n s wi rling bioconvec tion
n a nofluid flow fro m a d efo r m a ble ro t a ting disk
Kadir, A, Mish r a , SR, S h a m s h u d din, M a n d Be g, OA
Tit l e N u m e ric al si m ula tion of von Kar m a n s wirling bioconvec tion n a nofluid flow fro m a d efo r m a ble ro t a tin g disk
Aut h or s Kadir, A, Mis h r a, SR, S h a m s h u d din, M a n d Beg, OA
Typ e Confe r e nc e o r Works ho p It e m
U RL This ve r sion is available a t : h t t p://usir.s alfor d. ac.uk/id/e p rin t/47 8 9 1/
P u bl i s h e d D a t e 2 0 1 8
U SIR is a digi t al collec tion of t h e r e s e a r c h ou t p u t of t h e U nive r si ty of S alford. Whe r e copyrigh t p e r mi t s, full t ex t m a t e ri al h eld in t h e r e posi to ry is m a d e fre ely availabl e online a n d c a n b e r e a d , dow nloa d e d a n d copied for no n-co m m e rcial p riva t e s t u dy o r r e s e a r c h p u r pos e s . Ple a s e c h e ck t h e m a n u sc rip t for a ny fu r t h e r copyrig h t r e s t ric tions.
For m o r e info r m a tion, including ou r policy a n d s u b mission p roc e d u r e , ple a s econ t ac t t h e Re posi to ry Tea m a t : u si r@s alford. ac.uk .
concentration, micro-organism density number and gradients of these
functions at the disk surface (radial local skin friction, local
circumferential skin friction, Local Nusselt number, Local Sherwood
number, motile microorganism mass transfer rate). Extensive
interpretation of the results is included. The work provides a useful
benchmark for further computational fluid dynamics simulations of
nano-bioconvection rotating disk reactors.
Results & Discussion
The corresponding boundary conditions:
Selected computations are shown below. For brevity we have
restricted attention to stretching (S), unsteadiness (), and Stefan
blowing/suction (fw) effects. we note that although a mass transfer
Biot number, Nd appears in the boundary conditions it is prescribed
value of 0.4. This parameter is the solutal analogy to the thermal
Biot number. Whereas thermal Biot number represents the ratio of
internal thermal conduction resistance to external thermal
convection resistance, the mass transfer Biot number symbolizes
the relative contribution of internal mass diffusion resistance to the
external mass diffusion resistance. Since Nd <1, in our simulations
the external mass diffusion resistance dominates the internal mass
diffusion resistance, and this relates to the nano-particle
concentration field, not the micro-organism species.
Fig 1
Table 1: ADM Validation with RK 45 [8]
[1] F.S. Cross, Evaluation of a rotating disc type reservoir-oxygenator, Experimental Biology and Medicine, 93, 210-214 (1956). [2] J.R. Platt, Bioconvection patterns’ in cultures of free-swimming organisms. Science, 133, 1766–1767 (1961).[3] J. Buongiorno, Convective transport in nanofluids, ASME J. Heat Transfer, 128, 240–250 (2006). [4] G. Adomian, Solving Frontier Problems in Physics: The Decomposition Method, Kluwer, Dordrecht, USA (1994). [5] Kuznetsov, A.V., Thermo-bioconvection in a suspension of oxytactic bacteria, Int. Comm. Heat Mass Transfer, 32, 991-999 (2005).[6] O. Anwar Bég, Numerical methods for multi-physical magnetohydrodynamics, Chapter 1, pp. 1-112, New Developments in Hydrodynamics Research, M. J. Ibragimov and M. A. Anisimov, Eds., Nova Science, New York, September (2012).[7] O. Anwar Bég, D. Tripathi, T. Sochi and PK Gupta, Adomian decomposition method (ADM) simulation of magneto-bio-tribological squeeze film with magnetic induction effects, J. Mechanics Medicine Biology, 15, 1550072.1-1550072.23 (2015).[8] N.A. Latiff et al, Stefan blowing effect on bioconvective flow of nanofluid over a solid rotating stretchable disk, Propulsion and Power Research, 5, 267-278 (2016).
Fig 2
Ali Kadir and Dr. O. Anwar BégUniversity of SalfordDepartment of Aeronautical and Mechanical EngineeringNewton Building, Manchester, M5 4WT, UK.Email: [email protected] & [email protected]
S.R. Mishra and M. ShamshuddinDept. Mathematics, Vaagdevi College of Engineering, Warangal, Telangana, India.
Forced bioconvection in three-dimensional time-dependent viscous
incompressible nanofluid over a rotating stretchable disk is considered
in a cylindrical polar coordinate system Velocity componets in these
coordinates are The nanofluid is assumed to be a dilute suspension
with a homogenous distribution of gyrotatcic micro-organisms. Mass
convective and no-slip boundary conditions are imposed at the disk
surface. Stefan blowing is simulated via a w-(axial) velocity condition.
The rotating disk flow model is illustrated in Fig. 1.
The nano-particle species concentration conditions close to the disksurface follow cf>cw>c, in order to simulate mass convective boundaryconditions. Intrinsic to this is the presence of a different nano-particleconcentration from the neighborhood to the disk concentration, wallconcentration and ambient concentration which yields the required masstransfer coefficient hm. With appropriate similarity variables, theconservation boundary layer equations reduce to the following system ofstrongly coupled, nonlinear, multi-order, multi-degree ordinary differentialequations for primary and secondary momentum, energy, nano-particlespecies and micro-organism density species conservation:
02
12 22
ffSgffff
02
122
ggSfggfg
022
12 fPrPrSNtNb
022
111 fS
Nb
Nt
PrLePrLe
1 1
2 0Pr 2
S f PeLb
0 , 0 1, 0 0 , 0 1,Pr0 :
0 1 0 , 0 1
: 0, 0, 0, 0 0, 0
wff g fLe
Nd
f g
Here the non-dimensional parameters are defined as follows: S= β/ is
the unsteadiness (time-dependent) parameter, Pr is Prandtl number, Nb is
Brownian motion, Nt is thermophoresis parameter, Le is ordinary Lewis
number (ratio of thermal and nano-particle species diffusivities), Lb is
bioconvection Lewis number (ratio of thermal and motile micro-organism
mass diffusivities), Pe denotes bioconvection Peclet number (ratio of the
rate of advection of micro-organisms driven by the flow to the rate of
diffusion of micro-organisms under gyrotaxis), fw is the Stefan blowing
parameter, Nd is the Biot number and is the disk stretching parameter.
It is important to note the presence of nanoparticles does not modify the
incompressibility or biological nature of the suspension considered since
the nanoparticles do not interact with the motile (i.e. self-propelled)
microorganisms. They are a separate “species". The micro-organisms
possess mean diameters which range from 1mm to 200mm and are
therefore significantly larger than nanoparticle dimensions. These
microorganisms, which are considered to be gyrotactic in the present
study, possess slightly greater densities than the water-based fluid e.g.,
several percent for algae and no more than 10% for bacteria e.g., B.
subtilis, Bacillus, Chlamydomonas, Volvox and Tetrahymena. For rotating
disk bioreactor design [5], important transport characteristics include the
radial and tangential skin friction coefficients, local Nusselt number, local
Sherwood number, local wall motile micro-organism number (wall mass
flux gradient).
ADM Numerical SolutionThe eleventh order system of coupled ordinary differential equations(ODEs) with boundary conditions does not admit exact analyticalsolutions. Recourse is therefore necessary to computational methods.ADM [6] is one of the most accurate numerical methods. It uses veryprecise polynomial expansions to achieve faster convergence than manyother procedures and has been applied by for example Bég et al. [7](magneto-rheological squeeze films). An advantage of ADM is that it canprovide analytical approximation or an approximated solution to a wideclass of nonlinear equations without linearization, perturbation closureapproximations or discretization methods. ADM deploys an infinite seriessolution for the unknown functions and utilizes recursive relations. Theunknown functions f, g, , and are expressed as infinite series’ andthen the linear and nonlinear components are decomposed by an infiniteseries of Adomian polynomials as follows:
0000
000
00
2
00
0000
00
2
0
2
00
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
T,S,fR,Q
,fP,O,N
fML,K,J
,gI,gH,gfG,gfF
fE,fD,gC,fB,ffA
u,,rNd,r,t,
,qg,g,pf,f,rf
010100010
0100010
The Adomian transformed boundary conditions:
The power-series Adomian solutions are numerically evaluated in the subroutine bvp4cin MATLAB symbolic software. To verify the accuracy of the ADM solutions, comparisonsare made with the shooting (Runge-Kutta) quadrature computations of Latiff et al. [8] forthe general model (i.e. with all parameters included) with a variation in unsteadinessparameter, S (deceleration case i.e. S < 0) for two different scenarios, namely a non-deformable disk ( = 0) and a radially stretching disk ( = 1.0). These are documented inTable 1 for radial local skin friction, local circumferential skin friction, local Nusseltnumber, local Sherwood number, motile micro-organism wall mass flux for =0..
In Fig. 2 nano-particle concentration (volume fraction of nano-particle)profiles exhibit a more complex response to unsteadiness andblowing/suction parameters. Two different types of behaviour can bedelineated in the zones can be identified, namely the near disk region andthe far field region. In the near disk zone (i.e. 0 < < 1) initially withstrong blowing the nano-particle concentration magnitudes are elevated,and are significantly higher for weak disk deceleration (S = -1) ascompared with strong disk deceleration (S = -10). With further distancefrom the disk surface, these trends are reversed i.e. blowing is observedto reduce nano-particle concentration whereas suction is observed toslightly enhance it. Again the deviation in profiles for a specific S value ismarkedly greater for the weak disk deceleration situation compared withthe strong disk deceleration scenario. The diffusion of nano-particles intothe boundary layer is therefore a function of the location from the disksurface and is not a single consistent response throughout the entireregime. Fig. 3 reveals that a much more controlled response in motilemicro-organism density numbers is produced with a change inunsteadiness (S) and Stefan blowing/suction (fw) effects. The quantity ofmicro-organisms is decreased with suction and elevated with Stefanblowing. Strong deceleration again decreases motile micro-organismdensity number compared with weak disk deceleration i.e. the motilemicro-organism species boundary layer thickness is depleted with fasterdeceleration of the disk. Inevitably this is associated with the radial andazimuthal momenta fields which are damped simultaneously..