Entropy Analysis of Carbon Nanotubes Based Nanofluid Flow Past a Vertical Cone with Thermal Radiation

Abstract

Our objective in the present study is to scrutinize the flow of aqueous based nanofluid comprising single and multi-walled carbon nanotubes (CNTs) past a vertical cone encapsulated in a permeable medium with solutal stratification. Moreover, the novelty of the problem is raised by the inclusion of the gyrotactic microorganisms effect combined with entropy generation, chemical reaction, and thermal radiation. The coupled differential equations are attained from the partial differential equations with the help of the similarity transformation technique. The set of conservation equations supported by the associated boundary conditions are solved numerically with the bvp4c MATLAB function. The influence of numerous parameters on the allied distributions is scrutinized, and the fallouts are portrayed graphically in the analysis. The physical quantities of interest including the skin friction coefficient and the rate of heat and mass transfers are evaluated versus essential parameters, and their outcomes are demonstrated in tabulated form. For both types of CNTs, it is witnessed that the velocity of the fluid is decreased for larger values of the magnetic and suction parameters. Moreover, the value of the skin friction coefficient drops versus the augmented bioconvection Rayleigh number. To corroborate the authenticity of the presented model, the obtained results (under some constraints) are compared with an already published paper, and excellent harmony is achieved in this regard.

1. Introduction

Nanofluid, characterized by copious attractive features, including outstanding chemical and mechanical steadiness, significant improvement in thermal conductivity, etc. [1], is found to serve in a number of engineering applications, for example fuel-cells [2], porous materials [3], petroleum engineering [4], and biotechnology [5,6], among others. The pioneering work was done by Choi and Eastman [7] who found that thermal conductivity of the base fluid will increase from the insertion of metallic particles. This was followed by a study by Buongiorno [8] who studied the features of Brownian motion and thermophoresis in nanofluids. Later, Makinde and Aziz [9] deliberated on the flow of Newtonian fluid past a convectively heated surface. The flow of 3D couple stress nanofluid past an exponentially stretching surface associated with zero mass flux at the surface and convective boundary conditions was deliberated by Ramzan et al. [10]. Farooq et al. [11] examined Newtonian fluid flow analytically over an exponentially stretching sheet under the influence of magneto hydrodynamics using the optimal homotopy analysis method. The nanofluid flows containing carbon nanotubes (CNTs) over a cone and an inclined permeable plate were studied numerically by Reddy et al. [12,13]. Sreedevi et al. [14] found a numerical solution for CNTs amalgamated nanofluid flow past a vertical cone with a convective boundary condition. The aqueous-silver non-Darcy Poiseuille nanofluid flow with entropy generation past a permeable media was studied by Shehzad et al. [15]. A few recent investigations highlighting nanofluid flows may be found in References [16,17,18].

CNTs are the hexagonal structure of carbon atoms that are rolled in a cylindrical shape. Carbon nanotubes possess unique features like corrosion resistance, high thermal conductivity, and exceptional strength [19]. Owing to these remarkable characteristics, CNTs are useful in numerous applications like nanotubes transistors, microwave amplifier, solar cells, chemical sensors, optics, drug delivery, prostheses, pharmacogenomics, and many other fields of engineering and material science [20,21,22]. Carbon nanotubes are labeled as multi-walled carbon nanotubes (MWCNTs) and single-walled carbon nanotubes (SWCNTs). Iijima [23] discovered carbon nanotubes in 1991. He first investigated MWCNTs utilizing the Krastschmer and Huffman method. This was followed by another exploration in 1993 by Bethune [24] who introduced the concept of SWCNTs. SWCNTs are comprised of carbon nanotubes with a diameter of 1 nm whereas MWCNTs is a collection of 2–50 carbon nanotubes with 0.34 nm spacing. Abundant studies may be found in the literature to highlight different aspects of CNTs. Ramasubramaniam et al. [25] found that single-walled CNTs are helpful in improving electrical conductivity. The idea of improved thermal conductivity using composite nanotubes was introduced by Xue [26]. Muhammad et al. [27] inspected the rotating flow of carbon nanotubes under the influence of heat generation/absorption and nonlinear thermal radiation past a linearly stretching surface. The flow of 3D viscous nanofluid containing CNTs with quartic chemical reaction and entropy generation analysis is expressed numerically by Kumar et al. [28]. The aqueous based nanofluid Darcy-Forchheimer 3D flow comprising CNTs past a permeable surface was examined analytically by Muhammad et al. [29]. The flow problem in Reference [29] is extended to homogeneous-heterogeneous reactions associated with convective boundary conditions is discussed by Alshomrani and Ullah [30]. Recent explorations studying CNTs nanofluid flow may be found in References [31,32,33] and those contained therein.

The aforementioned literature review reveals that abundant articles are available addressing the topic of nanofluid. But this subject gets narrower once we talk about nanofluid flow over a cone with nanotubes inserted into it. Furthermore, this exploration becomes unique when the above-mentioned characteristics are supported by entropy generation and gyrotactic microorganisms (see

Table 1. The studies on nanoliquid flow comprising carbon nanotubes (CNTs).

Authors	CNTs SWCNTs/MWCNTs	Entropy Generation	Gyrotactic Microorganisms	Flow over a Cone
Reddy et al. [12]	√	×	×	√
Reddy et al. [13]	√	×	×	×
Sreedevi et al. [14]	√	×	×	√
Kumar et al. [28]	√	√	×	×
Muhammad et al. [29]	√	×	×	×
Alshomrani & Ullah [30]	√	×	×	×
Lu et al. [31]	√	×	×	×
Ramzan et al. [32]	√	√	×	×
Lu et al. [33]	√	√	×	×
Present	√	√	√	√

(√) means effect is present, and (×) means effect is absent.

). The numerical solution of the problem is acquired with requisite discussion of plotted illustrations of involved parameters versus associated distributions.

2. Mathematical Modeling

Let us assume a 2D aqueous fluid flow amalgamated with carbon nanotubes past a vertical cone in an absorbent media. The analysis is accompanied by solutal stratification, chemical reaction, and entropy generation (see Figure 1).

The flow of the fluid is along the x-axis past the cone surface. Along the y-axis, a magnetic field with strength B₀ is enforced. The fluid is an aqueous based nanofluid containing both types of CNTs, whose thermo-physical characteristics are defined in

Table 2. Values of physical features of nanoparticles and water [14].

Physical Attributes	Liquid	Nanoparticles
	H2O	SWCNTs	MWCNTs
Cp(J/kg K)	4179	425	796
$\rho$ (kg/m3)	997	2600	1600
k (W/mK)	0.613	6600	3000

Multi-walled carbon nanotubes (MWCNTs) and single-walled carbon nanotubes (SWCNTs).

. The governing system of equations representing the presented model are as follows [14]:

$$\frac{\partial (ru)}{\partial x}+\frac{\partial (rv)}{\partial y}=0,$$

(1)

$$\begin{array}{c}u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\frac{{\mu }_{nf}}{{\rho }_{nf}}\frac{{\partial }^{2}u}{\partial y^2}-\frac{{\mu }_{nf}}{{\rho }_{nf}}\frac{1}{K}u+g\left[\beta (T-T_{\infty })-{\beta }^{*}(C-C_{\infty })-{\beta }^{*}\gamma (n-n_{\infty })\mathsf{\Delta }\rho \right]\cos \gamma \\ -\frac{\sigma B_0{}^2}{{\rho }_{nf}}u,\end{array}$$

(2)

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}={\alpha }_{nf}\frac{{\partial }^{2}T}{\partial y^2}-\frac{1}{{(\rho c_p)}_{nf}}\frac{\partial q_r}{\partial y},$$

(3)

$$u\frac{\partial C}{\partial x}+v\frac{\partial C}{\partial y}=D_m\frac{{\partial }^{2}C}{\partial y^2}-K_r(C-C_{\infty }),$$

(4)

$$u\frac{\partial n}{\partial x}+v\frac{\partial n}{\partial y}+\frac{b{W}_c}{C_w-C_0}\frac{\partial }{\partial y}\left(n\frac{\partial C}{\partial y}\right)=D_n\frac{{\partial }^{2}n}{\partial y^2},$$

(5)

with the corresponding boundary conditions

$$\begin{array}{l}v=V_1,u=0,-k_{nf}\frac{\partial T}{\partial y}=h_f(T_f-T),C=C_w=C_0+dx,n=n_w,\mathrm{at}y=0,\\ u\to 0,T\to T_{\infty },C\to C_{\infty }=C_0+ex,n\to n_{\infty },\mathrm{as}y\to \infty .\end{array}$$

(6)

The hypothetical relations are characterized as follows:

$${\mu }_{nf}=\frac{{\mu }_f}{{(1-\varphi )}^{2.5}},v_{nf}=\frac{{\mu }_{nf}}{{\rho }_{nf}},$$

(7)

$${\rho }_{nf}=(1-\varphi ){\rho }_f+\varphi {\rho }_{CNT},{\alpha }_{nf}=\frac{k_{nf}}{{\rho }_{nf}{(c_p)}_{nf}},$$

(8)

$$\frac{k_{nf}}{k_f}=\frac{(1-\varphi )+2\varphi \frac{k_{CNT}}{k_{CNT}-k_f}\ln (\frac{k_{CNT}+k_f}{2k_f})}{(1-\varphi )+2\varphi \frac{k_f}{k_{CNT}-k_f}\ln (\frac{k_{CNT}+k_f}{2k_f})}.$$

(9)

Using the similarity transformations

$$\begin{array}{cc}\hfill \eta =& \frac{y}{x}R{a}_x{}^{1/4},\mathsf{\Psi }=\alpha R{a}_x{}^{1/4}f(\eta ),\theta (\eta )=\frac{T-T_{\infty }}{T_w-T_{\infty }},\hfill \\ & g(\eta )=\frac{C-C_{\infty }}{C_w-C_0},h(\eta )=\frac{n-n_{\infty }}{n_w-n_{\infty }},\hfill \end{array}$$

(10)

Equation (1) is impartially fulfilled and Equations (2) to (6) obtain

$$\begin{array}{cc}\hfill f^{\prime \prime \prime }+\frac{1}{\mathrm{Pr}}{(1-\varphi )}^{2.50}(1-\varphi +\varphi & \frac{{\rho }_{CNT}}{{\rho }_f})\left\{3f{f}^{\prime \prime }-\frac{1}{2}{f}^{{\prime }^2}\right\}-k_1{f}^{\prime }-{(1-\varphi )}^{2.5}Mf{}^{\prime }\hfill \\ & +{(1-\varphi )}^{2.50}(1-\varphi +\varphi \frac{{\rho }_{CNT}}{{\rho }_f})[\theta -N_{r}g-R_{b}h]=0,\hfill \end{array}$$

(11)

$$\frac{k_{nf}}{k_f}\left(1+R_d\right)\theta {}^{\prime \prime }+\frac{3}{4}\left[1-\varphi +\varphi \frac{{(\rho C_p)}_{CNT}}{{(\rho C_p)}_f}\right]f\theta {}^{\prime }=0,$$

(12)

$$g{}^{\prime \prime }+\frac{3}{4}{S}_{c}fg{}^{\prime }-S_{c}nf{}^{\prime }-C_{r}g=0,$$

(13)

$$h{}^{\prime \prime }+\frac{3}{4}{L}_{b}fh{}^{\prime }-P_e(h{}^{\prime }g{}^{\prime }+(h+\delta )g{}^{\prime \prime })=0,$$

(14)

and the boundary conditions (6) take the form

$$\begin{array}{l}f(0)=V_0,f^{\prime }(0)=0,\frac{k_{nf}}{k_f}\theta {}^{\prime }(0)=-B_1(1-\theta (0)),g(0)=1-n,h(0)=1,\\ f{}^{\prime }(\infty )\to 0,\theta (\infty )\to 0,g(\infty )\to 0,h(\infty )\to 0.\end{array}$$

(15)

In the aforementioned equations, the dimensionless parameters are given by:

$$\begin{array}{l}\mathrm{Pr}=\frac{{\upsilon }_f}{\alpha },k_1=\frac{x^2}{KR{a}_x{}^{1/2}},M=\frac{\sigma B_0{}^2{x}^2}{{\mu }_{f}R{a}_x{}^{1/2}},S_c=\frac{\alpha }{D_m},n=\frac{e}{d},\\ L_b=\frac{\alpha }{D_n},R_d=\frac{16T_{\infty }{}^3{\sigma }^{*}}{3k^{*}{k}_{nf}},N_r=\frac{{\beta }^{*}(C_w-C_0)}{\beta (T_f-T_{\infty })},R_b=\frac{{\beta }^{*}\gamma \mathsf{\Delta }\rho \mathsf{\Delta }{n}_w}{\beta (T_f-T_{\infty })},\\ C_r=\frac{K_r{x}^2}{D_{m}R{a}_x{}^{1/2}},B_1=\frac{h_{f}x}{R{a}_x{}^{1/4}{k}_f},P_e=\frac{b{W}_c}{D_n},\delta =\frac{n_{\infty }}{n_w-n_{\infty }}.\end{array}$$

(16)

The physically essential quantities, i.e., the skin friction, rate of heat and mass transfers, and local density of motile microorganisms, are appended below:

$$C_f=\frac{{\tau }_w}{\rho U_{\infty }{}^2},N{u}_x=\frac{x{q}_w}{k_f(T_w-T_{\infty })},S{h}_x=\frac{x{q}_m}{D_m(C_w-C_0)},N{n}_x=\frac{x{q}_n}{D_n(n_w-n_{\infty })}.$$

(17)

The aforementioned physical quantities in dimensionless form are appraised as follows:

$$\begin{array}{l}{C}_{f}R{a}_x{}^{1/4}=\frac{1}{{(1-\varphi )}^{2.5}}f\prime \prime (0),\\ N{u}_{x}R{a}_x{}^{-1/4}=-\frac{k_{nf}}{k_f}(1+R_d){\theta }^{\prime }(0),\\ S{h}_{x}R{a}_x{}^{-1/4}=-g{}^{\prime }(0),\\ N{n}_{x}R{a}_x{}^{-1/4}=-h{}^{\prime }(0).\end{array}$$

(18)

Table 3. Evaluation of the presented model with Khan et al. [34] in limiting case.

$\mathit{\varphi }$	${\mathit{f}}^{\prime \prime }(0)$				$-{\mathit{\theta }}^{\prime }(0)$
	Khan et al. [34]		Existing Results		Khan et al. [34]		Existing Results
	SWCNT	MWCNT	SWCNT	MWCNT	SWCNT	MWCNT	SWCNT	MWCNT
0.01	0.33894	0.33727	0.338910	0.337270	1.10553	1.07905	1.105710	1.079040
0.1	0.40811	0.39008	0.408120	0.390070	4.80627	4.27718	4.806290	4.277160
0.2	0.50452	0.46466	0.504530	0.464660	12.30317	10.56783	12.30352	10.56796

depicts a comparison with Khan et al. [34] for varied estimates of $\varphi$ in limiting case. An outstanding matching in both results is achieved. This reflects the corroboration of the presented outcomes.

Entropy Generation

The entropy generation of the presented model is specified as follows:

$$\begin{array}{c}{S}_{gen}{}^{{}^{\prime \prime \prime }}=\underset{HFI}{\underbrace{\frac{k_{nf}}{T_{\infty }{}^2}\left[1+\frac{16T_{\infty }{}^3{\sigma }^{*}}{3k^{*}{k}_{nf}}\right]{\left(\frac{\partial T}{\partial y}\right)}^2}}+\underset{FFI}{\underbrace{\frac{{\mu }_{nf}}{T_{\infty }}{\left(\frac{\partial u}{\partial y}\right)}^2+\frac{\sigma }{T_{\infty }}{B}_0{}^2{u}^2+\frac{{\mu }_{nf}}{T_{\infty }K}{u}^2}}\\ +\underset{Diffusiveirreversibility}{\underbrace{\frac{RD}{C_{\infty }}{\left(\frac{\partial C}{\partial y}\right)}^2+\frac{RD}{T_{\infty }}\left(\frac{\partial T}{\partial y}\right)\left(\frac{\partial C}{\partial y}\right)}}.\end{array}$$

(19)

In Equation (19), entropy is comprised of three terms, namely (i) HFI (heat transfer irreversibility), (ii) FFI (fluid friction irreversibility), and (iii) diffusion irreversibility. The entropy generation $N_G$ is defined as:

$$N_G=\frac{S^{{}^{\prime \prime \prime }}{}_{gen}}{S_0{}^{{}^{\prime \prime \prime }}},$$

(20)

where $S^{{}^{\prime \prime \prime }}{}_{gen}$ and $S_0{}^{{}^{\prime \prime \prime }}$ characterize the entropy generation rate and characteristic entropy generation rate, respectively, such that

$$\begin{array}{c}{N}_G=\frac{k_{nf}}{k_f}\left(1+R\right)R{a}_x\theta {{}^{\prime }}^2+\frac{1}{{(1-\varphi )}^{2.5}}\frac{BrR{a}_x}{\alpha }(f{{}^{\prime \prime }}^2+k_{1}f{{}^{\prime }}^2)+\frac{R{a}_{x}BrM}{\alpha }f{{}^{\prime }}^2\\ +\lambda {\left(\frac{\zeta }{\alpha }\right)}^{2}R{a}_{x}g{{}^{\prime }}^2+\frac{\zeta }{\alpha }R{a}_x\lambda \theta {}^{\prime }g{}^{\prime }.\end{array}$$

(21)

Parameter used in above equation are define as,

$$\alpha =\frac{\mathsf{\Delta }T}{T_{\infty }},Br=\frac{{\mu }_f{u}_w}{k_f\mathsf{\Delta }T},\zeta =\frac{\mathsf{\Delta }C}{C_{\infty }},\lambda =\frac{RD{C}_{\infty }}{k_f}.$$

(22)

3. Results and Discussion

This section is devoted to comprehend the discussion of the graphical illustrations. The impressions of the miscellaneous parameters on entangled profiles are given in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. The numerical values of the parameters used are taken to be fixed as: $\varphi =0.01,N_r=P_e=k_1=0.5=L_b,S_c=B_1=1.0=M=V_0,R_b=n=R_d=0.1=C_r=\delta \mathrm{and}\mathrm{Pr}=6.2$ unless otherwise stated. The ranges of parameters defined in the figures are $0.4\le M\le 1.0,$ $0.1\le k_1\le 0.7,$ $0.2\le N_r\le 0.4,$ $0.1\le R_b\le 0.3,$ $0.01\le \varphi \le 0.03,$ $0.5\le B_1\le 1.5,$ $0.1\le R_d\le 0.7,$ $0.5\le S_c\le 1.5,$ $0.1\le n\le 0.5,$ $0.5\le P_e\le 0.9,$ $0.5\le L_b\le 0.7,$ $0.1\le \alpha \le 0.3,$ and $0.1\le \lambda \le 0.5.$

3.1. Velocity Profile

The trend of axial velocity versus different parameters’ effects is described in Figure 2, Figure 3, Figure 4 and Figure 5. Figure 2 depicts the impact of the magnetic parameter $M$ on the velocity field. The velocity of the fluid diminishes for increasing values of $M$. This is because of the strong Lorentz force that presents resistance to the fluid’s movement that eventually lowers the fluid’s movement. In Figure 3, the consequence of porous parameter $k_1$ on the velocity profile is sketched. It is understood that the velocity is a decreasing function of $k_1$. Physically, more resistance against the fluid’s movement is witnessed due to the augmented thickness of the permeable medium that results in feeble fluid velocity. The impact of the buoyancy ratio parameter $N_r$ and the bioconvection Rayleigh number $R_b$ on the velocity profile for both CNTs is depicted in Figure 4 and Figure 5, respectively. It is witnessed that the velocity profile declines with increasing values of $N_r$ and $R_b$. Higher values of the buoyancy ratio parameter mean an increase in the number of CNTs immersed into the aqueous solution, which increases the viscosity of the fluid and results in a decrease in fluid’s velocity. Similarly, the velocity of the fluid is affected by the growth in the bioconvection Rayleigh number. This is due to the inertia force of the fluid motion being surpassed by the bioconvection. There is a decrease in speed of roughly 15.48% with an approximate increase in $R_b$ of 400% [35].

3.2. Temperature Profile

The outcome of solid volume fraction $\varphi$ on the temperature field is evident in Figure 6. The temperature profile is enhanced with increasing estimates of the solid volume fraction of the nanoparticles. It is also understood that the thermal boundary layer thickness is enhanced by augmenting the estimation of the solid volume fraction $\varphi$ for both nanotubes. This is all because of the enhancement in thermal conductivity of CNTs with solid volume fraction that becomes the main cause for augmented temperature. The effect of Biot number $B_1$ is studied in Figure 7. It is perceived that with an upsurge in $B_1$, temperature distribution escalates for both SWCNT and MWCNT. Physically, larger estimates of $B_1$ means more thermal resistance inside the cone in comparison to the boundary layer; consequently, a higher temperature of the fluid in the boundary layer area is witnessed. Figure 8 shows the influence of radiation parameter $R_d$ on temperature profile. It was determined that larger values of $R_d$ result in more energy being produced, which eventually raises the temperature of the fluid.

3.3. Concentration Profile

The impact of numerous parameters on the concentration field is presented in Figure 9 and Figure 10. Figure 9 depicts the influence of Schmidt number $S_c$ on concentration distribution for both nanotubes. It is understood from the figure that the concentration field is a diminishing function of $S_c$. Since the Schmidt number is the proportion of kinematic viscosity and the molecular diffusion coefficient, higher values of $S_c$ leads to a reduced molecular diffusion that ultimately lowers the concentration of the fluid. In Figure 10, the graph of concentration profile versus the solutal stratification $n$ is depicted. It is clear that for improved values of $n$, the concentration of the fluid is diminished for both SWCNT and MWCNT nanotubes. In actuality, the lowering of concentration field is due of the concentration differences between the ambient fluid and the surface of the cone.

3.4. Density of Motile Microorganism Profile

Figure 11 and Figure 12 demonstrate the impacts of the bioconvection Péclet number and the bioconvection Lewis number on the density of motile microorganisms, respectively. It is observed that motile microorganisms decrease for both bioconvection Péclet number and bioconvection Lewis numbers. Indeed, higher estimates of the bioconvection Péclet and bioconvection Lewis numbers result in a decline in the microorganism diffusion, which ultimately results in the decay of the density and boundary layer thickness of motile microorganisms.

3.5. Entropy Generation

From Figure 13, it is seen that increasing the temperature difference parameter $\alpha$ decreases the entropy generation number $N_G$ for both nanoparticles. Similarly, the local entropy generation increases for growing estimates of the diffusive constant parameter $\lambda$ for both SWCNT and MWCNT, which is displayed in Figure 14.

Table 4. Numerical value of $\frac{1}{{(1-\varphi )}^{2.5}}f{}^{\prime \prime }(0)$.

$\mathit{\varphi }$	${\mathit{k}}_1$	${\mathit{V}}_0$	${\mathit{R}}_{\mathit{b}}$	$\mathit{M}$	$\frac{1}{{(1-\mathit{\varphi })}^{2.5}}{\mathit{f}}^{\prime \prime }(0)$
					SWCNTs	MWCNTs
0.1	0.5	1.0	0.1	1.0	1.11420	0.57810
0.2					1.23160	1.00280
0.3					1.47580	1.11930
	0.2				1.29720	1.26970
	0.3				1.22650	1.16470
	0.4				1.16610	1.07700
		0.5			1.04840	0.87042
		0.6			1.07060	0.89693
		0.7			1.09350	0.92345
			0.2		1.10830	0.94259
			0.3		1.05020	0.88198
			0.4		0.99171	0.82097
				0.5	1.46240	1.25400
				0.6	1.38720	1.19320
				0.7	1.32120	1.13840

shows that the skin friction coefficient is enhanced with increase in the values of the solid volume fraction of nanoparticles and suction parameter, while it declines for higher values of the porous medium, magnetic parameter, and bioconvection Rayleigh number.

Table 5. Numerical value of $-\frac{k_{nf}}{k_f}(1+R_d){\theta }^{\prime }(0)$.

$\mathit{\varphi }$	${\mathit{R}}_{\mathit{d}}$	${\mathit{B}}_1$	$\mathit{M}$	$-\frac{{\mathit{k}}_{\mathit{n}\mathit{f}}}{{\mathit{k}}_{\mathit{f}}}(1+{\mathit{R}}_{\mathit{d}}){\mathit{\theta }}^{\prime }(0)$
				SWCNTs	MWCNTs
0.01	0.1	1.0	1.0	0.45621	0.45560
0.02				0.46268	0.46097
0.03				0.47205	0.46855
	0.2			0.47736	0.47659
	0.3			0.49760	0.49666
	0.4			0.51704	0.51594
		0.5		0.31751	0.31731
		0.7		0.38387	0.38351
		1.0		0.45621	0.45560
			1.0	0.45621	0.45560
			2.0	0.45279	0.45238
			3.0	0.45094	0.45063

demonstrates the numerical values of the Nusselt number for different varying parameters. It is found that the Nusselt number rises for larger estimates of solid volume fraction, radiation parameter, and Biot number, and decreases when the values of the magnetic parameter are increased.

Table 6. Numerical values of $-g^{\prime }(0)$.

${\mathit{S}}_{\mathit{c}}$	${\mathit{C}}_{\mathit{r}}$	$\mathit{n}$	${\mathit{N}}_{\mathit{r}}$	$-{\mathit{g}}^{\prime }(0)$
				SWCNTs	MWCNTs
0.1	0.1	0.1	0.5	0.31891	0.31882
0.5				0.50221	0.50155
0.9				0.74207	0.74087
	0.1			0.80642	0.80511
	0.2			0.88714	0.88613
	0.3			0.95695	0.95612
		0.2		0.73573	0.73379
		0.3		0.66795	0.66532
		0.4		0.60326	0.59988
			0.6	0.79903	0.79771
			0.7	0.79130	0.78997
			0.8	0.78319	0.78185

displays the numerical value of the Sherwood number for varied parameters. The Sherwood number boosted with an increase in the values of the chemical reaction parameter and the Schmidt number, while it decreases with an increase in the values of the concentration stratification and buoyancy ratio parameter.

Table 7. Numerical values of $N{n}_{x}R{a}_x{}^{-1/4}$.

${\mathit{L}}_{\mathit{b}}$	${\mathit{P}}_{\mathit{e}}$	${\mathit{R}}_{\mathit{b}}$	$\mathit{\delta }$	$-{\mathit{h}}^{\prime }(0)$
				SWCNTs	MWCNTs
0.5	0.5	0.1	0.1	0.83681	0.83535
0.6				0.91358	0.91207
0.7				0.99024	0.98870
	0.1			0.48494	0.48380
	0.2			0.57267	0.57146
	0.3			0.66056	0.65927
		0.2		0.82280	0.82132
		0.3		0.80759	0.80609
		0.4		0.79092	0.78938
			0.2	0.87679	0.87530
			0.3	0.91679	0.91528
			0.4	0.95681	0.95528

depicts the numerical value of motile density number versus different parameters. The motile density number increases for larger estimates of the Péclet number and microorganism concentration difference parameter and decreases for increasing values of the Rayleigh number.

4. Final Remarks

The flow of water based carbon nanotubes (SWCNT and MWCNT) fluid past a cone erected vertically is discussed numerically. The analysis is performed in the presence of motile organisms with solutal stratification in spongy media. Furthermore, the attributes of thermal radiation and chemical species are explored in the presence of entropy generation. The main outcomes of the analysis are:

➢

The velocity of the fluid diminishes with increasing values of the magnetic and suction parameters in the case of both nanotubes.

➢

The fluid’s concentration is diminished for both SWCNT and MWCNT nanotubes versus higher values of solutal stratification.

➢

When increasing the temperature difference parameter, the entropy generation number decreases for both nanoparticles.

➢

The Sherwood number increases with increasing values of the chemical reaction parameter and the Schmidt number, while it decreases with increasing estimates of solutal stratification.

➢

The motile density number decreases with increasing values of the Péclet number.

➢

The skin friction coefficient increases for the suction parameter while decreasing for the bioconvection Rayleigh number.

➢

It is found that the Nusselt number increases with an increase in the estimates of solid volume fraction, radiation parameter, and Biot number, whereas it decreases with increasing values of the magnetic parameter.