#### Abstract

Thermobioconvection boundary layer flow in a suspension of water-based bionanofluid holding both nanoparticles and motile microorganisms past a wedge surface was studied. The governing nonlinear partial differential equations on reference of the Buongiorno model were transformed into a set of coupled nonlinear ordinary differential equations. Shooting technique was then used to solve the transformed nonlinear ordinary differential equations numerically. The solutions were found to be contingent on several values of the governing parameters. As highlighted, the velocity profile as well as the skin friction coefficient was affected by the pressure gradient parameter, the function of the wedge angle parameter. On the other hand, the temperature, nanoparticle concentration, and density of motile microorganism’s distributions together with its corresponding local Nusselt number, local Sherwood number, and local density of the motile microorganisms change with the thermophoresis and Brownian motion parameter and so Lewis number, Schmidt number, and bioconvection Péclet number. An experimental scheme together with sensitivity analysis on the basis of Response Surface Methodology (RSM) was applied to examine the dependency of the response parameters of interest to the input parameters’ change. Obviously, local Nusselt number was more sensitive towards the Brownian motion parameter when the Brownian motion parameter was at 0.2 and 0.3. However local Sherwood number was more sensitive towards the Lewis number for all values of Brownian motion parameter. Compatibility found by comparing results between RSM and shooting technique gave confidence for the model’s accuracy. The findings would provide initial guidelines for future device fabrication. Finally, the numerical results obtained were thoroughly inspected and verified with the existing values reported by some researchers.

#### 1. Introduction

Nanofluid is a term first innovated by Choi and Eastman [1] showing a popular thermal conductivity suspension consisting of nominal 1-100 nm size particles distributed in a conventional base fluid. Since then, the research on nanofluid focusing on its thermophysical properties in the boundary layer has not stopped but increased among researchers. The interest stems from the high thermal conductivity characteristics possessed by nanofluid for heat transfer intensification which results in energy saving. Such feature is significantly crucial and essential to a wide variety of applications and has extensive future possibilities. Recent review articles recorded the uses of nanofluid in solar industry, automobile industry, mechanical industry, electronic industry, biomedical industry, heat and cool exchanger system, geothermal system, etc. [2, 3].

On the other hand, bionanofluid refers to water-based fluid that consists of both nanoparticles and living microorganisms. The living microorganisms are found to be denser than the water-based fluid and are likely to swim in upward direction. Bioconvection is a process of heat transfer impelled from the collective upswimming motion by the motile microorganisms due to the unstable density delamination [4–6]. The mechanism of impelling upswimming motion is closely related to the species of the living microorganisms which responds to the external stimulus. For example, oxytaxis responds to the oxygen concentration gradient, phototaxis responds to light intensity, and gyrotaxis microorganisms respond to the combination of viscous torque and gravitational pull. Bioconvection is developed as a significant mechanism in bioenergy systems and biomicrosystems [7, 8]. It is worth mentioning that various features and aspects of bioconvection with suspension of nanoparticles have been conducted in the field of biomedicine and biomedical engineering. This action is aimed towards a better state of health. Therefore, the knowledge of bionanofluid thermophysical properties becomes essential when a health treatment is considered.

Wedge is a shape that is triangular in cross section. In fluid dynamic, when a surface appoints to an angle (wedge) and the free stream diverges from the surface, there would be a different fluid velocity distribution. The concept has facilitated enormous applications in the real world problem and thus gives impetus to researchers for further exploitation. Boundary layer problem of a wedge flow has large bodies of literature with Falkner and Skan [9] being the first. The extension of boundary layer flow over a shifting wedge from a static wedge has been done by Ishak et al. [10]. Michael and Boyd [11] add a slip boundary condition to the Falkner-Skan equations to determine the wedge flow. Then, there are abundant amounts of studies on the problem of boundary layer flow over a static or moving wedge devoted in nanofluid with different boundary condition considered, view [12–18]. However, the study of boundary layer flow over a wedge in a fluid containing both nanoparticles and microorganisms is not fully explored. Khan et al. [19, 20] filled the gap where nanofluid flowed across a porosity wedge in the existence of living microorganism where MHD effects is taken into account. Meanwhile, Mahdy [21] explained the consequences of the microorganisms’ behavior on nanofluid flow over an upright porous wedge.

Motivated by the above contributions, the authors intend to investigate the behaviors and properties of a water-based bionanofluid flow over a wedge surface by employing the simplest possible boundary conditions. This investigation is associated with the possible utilization of bionanofluid in designing microfluidic devices such as biosensors for bioapplications. Similarity transformation followed by shooting method is employed to solve the governing equations numerically. Afterwards, experimental scheme is applied together with sensitivity analysis on the basis of Response Surface Methodology (RSM) to examine the dependency of the local Nusselt number and local Sherwood number on the three governing parameters, namely, Lewis number, thermophoresis parameter, and Brownian motion parameter. This action is to assess the validity of the numerical results in order to provide a guideline in device fabrication. To the best of our knowledge, no study has been performed yet for such analysis.

#### 2. Mathematical Formulation

##### 2.1. Governing Equations

A steady two-dimensional boundary layer flow of water-based bionanofluid along the wedge surface with different configuration is considered. The problem is clearly presented through the schematic diagram as in Figure 1. The Cartesian coordinate system shows - and -axes are extended parallel and normal to the wedge surface with apex as the origin respectively, associated with corresponding velocity components and . Taking and as constants where and signified to the wedge surface velocity move in the water-based bionanofluid and the free stream velocity, respectively. , , and are the variables of the uniform temperature, nanofluid volume fraction and concentration of microorganisms, while the subscripts and refer to the corresponding values near the wedge surface and distance field far away from the wedge surface, respectively. Moreover, total angle of the wedge is subjected to , where stands for the Hartree pressure gradient expressed as the function of wedge angle parameter such that .

The governing partial differential equations of continuity, momentum, thermal energy, nanoparticle concentration, and microorganisms density for the boundary layer flow of the bionanofluid are formulated as follows, on account of the Buongiorno model which includes effect of Brownian motion and thermophoresis:Here, stands for the thermal diffusivity of nanofluid, represents the ratio of effective heat capacity of the nanoparticles to the base fluid, is recognized as the Brownian diffusion coefficient, and conveys the thermophoretic diffusion coefficient while signifies the diffusivity of microorganisms. Then expresses the kinematic viscosity, , where and are exemplified as the chemotaxis constant and maximum swimming speed of microorganisms.

The equations are appointed to the following boundary conditions:

To obtain similarity solutions for the system of governing equations (1)-(5), similarity variables are indicated as below [14, 19]:where symbolizes the stream function introduced as and which satisfied the continuity equation (1).

By replacing (7) into (1)-(5), a set of ordinary differential equations are obtained as follows:The boundary conditions (6) becomeHere, , , , and represented the Prandtl number, Lewis number, bioconvection Péclet number, and Schmidt number, respectively. , , and explained the Brownian motion parameter, thermophoresis parameter, and a dimensionless constant. They are conveyed as

Skin friction coefficient , local Nusselt number , local Sherwood number , and local density of the motile microorganisms are the parameters of interest which are stipulated aswhere , , , and are, respectively, recognized as the surface shear stress, wall heat flux, wall mass flux, and wall motile microorganisms flux with the following definition:By using (7), (13), and (14), we obtainedwhere denotes as the local Reynolds number, signifies the skin friction coefficient , represents the local Nusselt number , symbolizes the local Sherwood number , and stands for the local density of the motile microorganisms .

##### 2.2. Experimental Scheme

In mathematics, an experiment is a computer simulation by means of a series of tests or runs designed by using computer toolbox software. This numerical experiment aims to investigate the output result of a computer code change due to the input factors and thus derive conclusions about the prominence of the factors. Response Surface Methodology (RSM) is the mathematical and statistical technique introduced by Box and Wilson [25] which is used to determine the model’s uncertainty and sensitivity. RSM’s main aim is to optimize the model by identifying the relationship among the output and input variables. As observed in Table 2, the local Nusselt and Sherwood numbers vary with Lewis number, thermophoresis, and Brownian motion parameter. For the purpose of identifying the correlations between the output response parameters (local Nusselt and Sherwood numbers) to the independent input parameters (Lewis number, thermophoresis, and Brownian motion parameter), mathematical experiment is designed by referring to a standard nonlinear polynomial model [26]. denotes the response local Nusselt number or local Sherwood number that rely on the three independent input parameters of the Lewis number, thermophoresis and Brownian motion parameter coded with symbols , and , respectively. One intercept term, three linear terms, three two-factor bilinear terms, and three squared terms are embedded in the following mathematical equation:

According to the RSM, twenty runs together with nineteen Degrees of Freedom are appropriate for selected three levels of parameters. These three levels are classified into low (-1), medium (0), and high (1). Table 3 shows the three input parameters with their corresponding symbols and levels.

Central Composite Face centered design (CCF) employed in MINITAB-18 was utilized in this study to determine the experimental result. 20 runs are required for the experimental plan by referring to the expression of , whereby is the number of factors and is the number of center points. The series of experimental plan is presented in Table 4.

The degree of model accuracy was predetermined by the ANOVA analysis, also named as the analysis of variance on the RSM model. ANOVA analysis is a statistical technique significance for the variability performance of particular variables’ dependency whereby Degree of Freedom (DOF), sum of squares, adjusted mean square, F-value, and P-value are the computing estimators. Table 5 displays the ANOVA analysis for both the local Nusselt and the Sherwood numbers.

Sensitivity is broadly defined as the derivation of the response function with respect to the model parameters. Sensitivity analysis is the study of uncertainty criteria yielded in the output of the model apportioned by the input factors which leads to the assessment of the model’s robustness. Such analysis simulates information about the key drive parameters on the model’s result which then benefits engineers in designing processes.

#### 3. Results and Discussion

Numerically computed solutions of the transformed ordinary differential equations (8)-(11) coupled with the boundary conditions (12) are accomplished by utilizing the shooting technique programmed in Maple18, the computer algebra software. The boundary layer thickness is set equal to 8. In addition, the Prandtl number is fixed at 6.2 (water at 25°C) to ensure the survivability of the living microorganisms in the medium of conductivity. The results’ veracity is verified with the results reported by previous researchers. As shown in Table 1, comparison of the values of skin friction coefficient, , for different values of in the absence of microorganisms density equation shows a favourable agreement. Table 2 presents the values of skin friction coefficient , local Nusselt number , local Sherwood number , and local density of the motile microorganisms for different values of , , , , , , and . It is observed that will only be affected by . Note that indicates a horizontal flow whereas indicates a vertical flow. increases as the value of increases, showing that the surface shear stress rises with wedge. The correlation between skin friction coefficient and the velocity profile on the pressure gradient parameter can be discovered in the transformed momentum equation (8). Then, by considering a constant , the response parameters, , , and , change with the governing parameters, , , , , , and . On the other hand, and remain constant when , , and are fixed at 1 and 0.2, respectively. In other words, variation of the local Nusselt and Sherwood numbers are intimately related on the three governing parameters. In addition, the local density of motile microorganisms increases with , , , , , and signifies the motile microorganisms performed a growing manner.

Table 5 demonstrates the ANOVA analysis to identify the correlations between the local Nusselt and Sherwood numbers to the three independent input parameters. In ANOVA analysis, F-value is the measure of data variation about the average value whereas P-value is the probability validation of the model accuracy from the statistical viewpoint. Large F-value designates a significant result while small P-value assigns strong support on the result’s significance. Consequently, F-value is always used alongside the P-value to obtain strong evidence on the result’s significance. Thereupon, impacts of the linear terms, two-factor bilinear terms, and squared terms for both response parameters (local Nusselt number, , and local Sherwood number, ) are said to be statistically significant with strong evidence of large F-value and small P-value. In general, residual error is the unexplained data point by the regression line, whereas Lack-of-fit exhibits when the model fails to describe the functional connectedness between the input parameters and the output response. At the same instant, normal probability plot of the residual for the local Nusselt and Sherwood number are displayed in Figures 2 and 3. Both plots that appear with a straight line illustrating the errors are normally distributed and so the regression model is well fitted.

The regression coefficients with its corresponding P-value for the nonlinear polynomial model in (17) for both response parameters (local Nusselt number and local Sherwood number ) are marked in Table 6. It is worth mentioning that large P-value is said to be not statistically significant which means that no associate change can be detected on the output due to the change in input. Conversely, the term with small P-value that is statistically significant elsewhere can be neglected. As a result, , , , , , , , , and are meaningful factors for while , , , , , and are important terms for . Hence, mathematical equation (17) can be rewritten asFurthermore, the values of R-squared, and adjusted R-squared, and presented in Table 6 provide information about the goodness-of-fit of the RSM model. Notice that both and are reported with high and values which are nearly close to 100% indicating a perfectly predicted correlation between the output response parameters and the independent input parameters.

As mentioned earlier, sensitivity is the partial derivative ratio of the response function with regard to the model parameters. Hence, sensitivity functions of local Nusselt number and local Sherwood number with respect to the three governing parameters, Lewis number , thermophoresis parameter , and Brownian motion parameter based on (18) and (19) are expressed asUsing (20)-(25), the sensitivity results of the responses ( and ) with relation to the input governing parameters (, , and ) are calculated as in Table 7. Noteworthy, a positive sensitivity value exemplifying increase of input governing parameter causes an increase in the response function and vice versa for the negative sensitivity value. Results of the sensitivity analysis have been plotted into vertical bar charts (Figures 4 and 5) for better comprehension.

**(a) Sensitivity of Nu_x, for A=0 and B=-1**

**(b) Sensitivity of Nu_x, for A=0 and B=0**

**(c) Sensitivity of Nu_x, for A=0 and B=1**

**(a) Sensitivity of Sh_x, for A=0 and B=-1**

**(b) Sensitivity of Sh_x, for A=0 and B=0**

**(c) Sensitivity of Sh_x, for A=0 and B=1**

As observed in Figures 4(a)–4(c), an overall trend shows sensitivity of the local Nusselt number that falls with increasing governing parameters for all values of the Brownian motion parameter. Yet, sensitivity of the local Nusselt number raises with increasing Brownian motion parameter values from 0.2 to 0.4 ( to 1). Also, sensitivity of the local Nusselt number towards the Brownian motion parameter is higher than the Lewis number and thermophoresis parameter for and . For the case of high Brownian motion parameter , local Nusselt number seems to have a high sensitivity towards the Lewis number instead of Brownian motion and thermophoresis parameter.

Sensitivity of local Sherwood number towards different governing parameters at different values of Brownian motion parameter is shown in Figures 5(a)–5(c). In general terms, sensitivity of the local Sherwood number towards Lewis number and thermophoresis parameter increases with increasing Brownian motion parameter but the sensitivity towards Brownian motion parameter is kept constant despite the increased Brownian motion parameter. Besides, sensitivity of the local Sherwood number towards Lewis number and thermophoresis parameter shows a decreasing manner for values of from -1 to 1. To discuss sensitivity of the local Sherwood number towards the Brownian motion parameter in details, Figure 5(a) shows similar positive sensitivity at low thermophoresis parameter for all values of Brownian motion parameter. However, zero sensitivity of local Sherwood number towards Brownian motion parameter is found at thermophoresis parameter equal to 0.3 (Figure 5(b)). Then, Figure 5(c) displays that, at high thermophoresis parameter, similar negative sensitivity of the local Sherwood number towards the Brownian motion parameter is found for all values of Brownian motion parameter.

Accuracy of RSM model is benchmarked by comparing and results between RSM and shooting method for and . Figures 6(a)-6(b) show a good compatibility between RSM results and shooting results for the local Nusselt number and local Sherwood number. The maximum errors found between RSM and shooting results for and are 1.52% and 1.47%, respectively.

**(a) Nu_x, for A=1 and C=0**

**(b) Sh_x, for A=1 and C=0**

#### 4. Conclusion

A numerical study on bioconvection boundary layer flow in a bionanofluid across a wedge surface has been carried out. Bionanofluid is believed to have a positive result on the heat and mass transfer intensification in term of stability which is crucial in designing microfluidic devices. The governing nonlinear partial differential equations are transformed to ordinary differential equations by using similarity transformation followed by the shooting technique. The numerical results then undergo a sensitivity analysis to extrapolate the dependency of the parameters of interest on the different governing parameters. The findings provided initial guidelines for future device fabrication. The findings are highlighted as follows:(i)Wedge surface with different value of Hartree pressure gradient, , influences the skin friction coefficient and velocity profile.(ii)The rate of heat transfer directly proportional to the local Nusselt number decreases with Lewis number, thermophoresis, and Brownian motion parameter.(iii)The nanoparticle concentration directly proportional to the local Sherwood number increases with Lewis number, thermophoresis, and Brownian motion parameter.(iv)Sensitivity of the local Nusselt number decreases with increasing of the three governing parameters (Lewis number, thermophoresis parameter, and Brownian motion parameter).(v)The local Nusselt number is more sensitive towards the Brownian motion parameter by increasing the values of Brownian motion parameter from 0.2 to 0.3. Yet, the local Nusselt number is more sensitive towards Lewis number when Brownian motion parameter is high which is equal to 0.4.(vi)Sensitivity of the local Sherwood number increases with the increasing governing parameters, namely, Lewis number and thermophoresis parameter.(vii)The local Sherwood number is more sensitive towards the Lewis number for three different values of thermophoresis parameter and all values of Brownian motion parameter.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

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

#### Acknowledgments

Authors would like to express gratitude to the Research Management Centre (RMC) Universiti Tun Hussein Onn Malaysia (Project Code: GPPS-U812) and the Ministry of Higher Education, Malaysia (Project Code: FRGS-1644), for the financial support.