Nonlinear Mixed Convective Flow over a Moving Yawed Cylinder Driven by Buoyancy

Abstract

The fluid flow over a yawed cylinder is useful in understanding practical significance for undersea applications, for example, managing transference and/or separation of the boundary layer above submerged blocks and in suppressing recirculating bubbles. The present analysis examines nonlinear mixed convection flow past a moving yawed cylinder with diffusion of liquid hydrogen. The coupled nonlinear control relations and the border restrictions pertinent to the present flow problem are nondimensionalized by using nonsimilar reduction. Further, implicit finite difference schemes and Quasilinearization methods are employed to solve the nondimensional governing equations. Impact of several nondimensional parameters of the analysis on the dimensionless velocity, temperature and species concentration patterns and also on Nusselt number, Sherwood number and friction parameter defined at the cylinder shell is analyzed through numerical results presented in various graphs. Velocity profiles can be enhanced, and the coefficients of friction at the surface can be reduced, for increasing values of velocity ratio parameters along chordwise as well as spanwise directions. Species concentration profile is reduced, while the Sherwood number is enhanced, for growth of the Schmidt number and yaw angles. Furthermore, for an increasing value of yaw angle, skin-friction coefficient in chordwise direction diminishes in opposing buoyancy flow case, whereas the results exhibit the opposite trend in assisting buoyancy flow case. Moreover, very importantly, for increasing magnitude of nonlinear convection characteristic, the liquid velocity and surface friction enhance in spanwise direction. Further, for increasing magnitude of combined convection characteristics, velocity profiles and coefficient of friction at the surface enhance in both spanwise and chordwise directions. Moreover, we have observed that there is no deviation for zero yaw angle in Nusselt number and Sherwood number.

1. Introduction

Many researchers over the last few decades have given considerable attention to the investigation of heat and mass transfer characteristics of combined convection in various flow geometries. Combined convection flows appear when the temperature and species concentration variations between the wall and the external liquid are larger and, hence, become important when buoyancy forces significantly disturb the circulation, heat and species concentration patterns. When the fluid is subjected to two various density drops having various rates of diffusion, it is meant to be the double diffusion convection. The density differences may be caused by gradients in the liquid concentration, or by the changes in the temperature. The double diffusive mixed convection plays a vital role in boundary layer flow problems because of its significance in several technical and geophysical challenges including solar collectors, solar ponds, lakes, reservoirs and crystal growth [1]. Due to its vast applications, many researchers have worked on the double diffusion combined convection flow past various geometries such as sphere, exponential stretching sheet, vertical cone, slender cylinder, moving plate, etc. [2,3,4,5,6,7].

However, for viscous liquid flows with heat transfer, the impact of linear dependence of density on dimensionless temperature, that is, natural convection, seems to be highly significant in applications pertaining to industrial manufacturing processes and, therefore, cannot be ignored. The “temperature and density relation” is nonlinear for a large distinction between the border and liquid temperatures, the nonlinear density temperature differences in the buoyancy force term have a substantial effect on the circulation and energy transference features. Vajravelu and Sastri [8] have analyzed the heat transfer characteristics between the vertical borders with and without the nonlinear density temperature differences. Bhargav and Agarwal [9] have examined the fully developed free convection with temperature-dependent density within a duct. Mosta et al. [10] have examined the influence of nonlinear temperature combined with density differences in nanosuspension convection over a vertical border. In the present analysis, the “nonlinear temperature and density relation” is regarded owing to the important variation between the liquid and wall temperatures. At the same time, the liquid hydrogen diffusion [3,5] is regarded because of its cooling ability.

The boundary layer concept is the most useful and important aspect in understanding the transport processes occurring in external flows. The phenomenal results of the mixed convection flow along various geometries have been contributed by various researchers around the globe. Recently, Muthukumaran and Bathinathanan [11] have worked on mixed convection boundary layer flow over a stretching sheet, and their results have revealed that the heat transfer rate increases with a raise of the Prandtl number for both assisting and opposing flows. Halim and Noor [12] have analyzed mixed convection over a vertical stretching sheet, and their outcomes have revealed that the assisting flow has higher rates of heat and mass transfer compared to the opposing flow. Alsabery et al. [13] have investigated the mixed convection through a rotating cylinder. Khashiie et al. [14] have worked on mixed convection flow over a Riga plate. Complex flow patterns over yawed cylindrical cables, suppression of fluctuations in lift forces and the control of the drag forces are the main problems encountered in the engineering design. At the same time, the investigation of combined convection circulation about a tilted cylinder has not been analyzed so far. The examination of motion along tilted cylinder is very useful for heat exchangers design [15]. The fluid flow past yawed and unyawed cylindrical geometries extensively occurs in various engineering-related applications, like tow cables, chimney stacks, different towers, sub-sea pipelines, risers, heat exchangers and overhead cables [16,17]. King [16] has studied the vortex exited oscillations of yawed circular cylinders. This study reveals that the sustained oscillations can be observed at the yaw angles ranging between ±65°. Najafi et al. [18] have studied the undisturbed flows over the yawed cylinder, wherein it is observed that eddies near the wake of the cylinder detach due to the increase in the angle of yaw. Further, Snarski [19] has investigated the variation of wall pressure on circular cylinder due to yaw angle and has revealed that the spectra exhibit powerful narrow band energy stages related with the Strouhal vortex shedding region, for the degrees of yaw angle ranging from π/3 to π/2. In the work of Sears [20], the boundary layer motion along a yawed cylinder is analyzed, and it is noticed that the location of the boundary layer separation point is independent from the yaw angle. Moreover, Thapa et al. [21] have investigated numerically the circulation over a yawed circular cylinder near the plane border. In the work of Gupta and Sarma [22], the time-dependent circulation along a yawed infinite cylinder under the influence of cross flow has been analyzed, and this analysis reveals that the coefficient of friction along chordwise direction diminishes with larger values of nonsimilar variable (ξ). Further, Bucker and Lueptow [23] have examined the boundary layer flow on weakly yawed cylinders, and they have found that the thickness of boundary layer rises nonlinearly for small degree yaw angle. In this direction, many of the researchers have worked on the boundary layer flow past a yawed cylinder with the influence of non-uniform suction/injection [17,24,25,26,27]. Marshall [28] has studied the disturbed flow along a moving yawed cylinder and found that, as cross stream Reynolds number increases, the surface vorticity enhances. Vakil and Green [29] have performed numerical calculations pertaining to the flow over a yawed cylinder for moderate values of Reynolds number, i.e., 1 ≤ Re ≤ 40, and have validated the independence principle and also proposed an empirical relation for the lift and drag forces on the cylinder. Recently, Patil et al. [30] have scrutinized the combined convection circulation along a yawed cylinder, and the obtained data reveal that the liquid velocity and surface drag coefficient at the cylinder’s boundary in all directions rises with the heat transfer rate because of combined convection.

From the above literature review, it is found that the investigation of double diffusive combined convection past a moving yawed cylinder has not been attempted so far. Since, many authors such as Roy [25], Chiu and Lienhard [31], Roy and Saikrishnan [26] and Revathi et al. [17] have analyzed a boundary layer flow over a yawed cylinder. The works of these researchers have stimulated us to work on the present article, and, therefore, we have scrutinized the considered challenge. However, the present problem is formulated as an endeavor to investigate the steady combined convection boundary layer flow around a moving yawed cylinder. The novelty in the present analysis is as follows:

−

Convective flow over a moving yawed cylinder driven by buoyancy.

−

Influence of liquid hydrogen diffusion.

−

Effects of yawed angle.

−

Flow characteristics in chordwise and spanwise directions.

The considered governing equations and boundary conditions have been reduced using nonsimilar technique. Further, these equations are solved by employing the quasilinearization method and the implicit finite difference schemes [32,33,34].

2. Mathematical Simulation

Herein, we analyze the viscous, laminar, incompressible, combined convective flow over a moving yawed cylinder. The flow system is demonstrated in Figure 1, where the liquid (water) is supposed to move over a vertically yawed cylinder with radius R so that the yaw angle θ is used between 0 and π/6. Here, θ = 0 represents the vertical cylinder, while θ = π/2 represents the horizontal cylinder. To have the influence of buoyancy, the cylinder must be regarded in vertical or tilted location because the investigation is of combined convection. Therefore, the yaw angle is regarded in the range 0 ≤ θ ≤ π/2. The magnitude of angle of yaw above π/6 is not regarded because it would be nearer to stagnation point circulation and it is not the aim of the current examination of combined convection flow. For analysis, x and z are the coordinate axes in chordwise and spanwise directions, respectively, with u and w denoting the corresponding velocity components. Moreover, y is the coordinate axis drawn normal to x and z axes with v denoting the corresponding velocity component. The temperature and species concentration of the liquid at the border are denoted by T_s and C_s, while that away from the surface are represented by T_∞ and C_∞. The density changes are modeled using the Boussinesq approach [5,35].

The governing equations representing the fluid flow variations are as follows [17,24,25,26,27,30,36]

Continuity equation

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$$

(1)

Momentum equation along chordwise direction

$$\begin{array}{l}u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=u_e\frac{d{u}_e}{dx}+g\left\{{\mathsf{\beta }}_1\left(T-T_{\infty }\right)+{\mathsf{\beta }}_2{\left(T-T_{\infty }\right)}^2\right\}\sin \left(\mathsf{\theta }\right)\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}+\mathsf{\nu }\frac{{\partial }^{2}u}{\partial y^2}+g\left\{{\mathsf{\beta }}_3\left(C-C_{\infty }\right)+{\mathsf{\beta }}_4{\left(C-C_{\infty }\right)}^2\right\}\sin \left(\mathsf{\theta }\right)\end{array}$$

(2)

Momentum equation along spanwise direction

$$\begin{array}{l}u\frac{\partial w}{\partial x}+v\frac{\partial w}{\partial y}=\mathsf{\nu }\frac{{\partial }^{2}w}{\partial y^2}+g\left\{{\mathsf{\beta }}_1\left(T-T_{\infty }\right)+{\mathsf{\beta }}_2{\left(T-T_{\infty }\right)}^2\right\}\cos \left(\mathsf{\theta }\right)\\ \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}+g\left\{{\mathsf{\beta }}_3\left(C-C_{\infty }\right)+{\mathsf{\beta }}_4{\left(C-C_{\infty }\right)}^2\right\}\cos \left(\mathsf{\theta }\right)\end{array}$$

(3)

Energy equation

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\frac{k}{\mathsf{\rho }{c}_p}\frac{{\partial }^{2}T}{\partial y^2}$$

(4)

Concentration equation

$$u\frac{\partial C}{\partial x}+v\frac{\partial C}{\partial y}=D_B\frac{{\partial }^{2}C}{\partial y^2}$$

(5)

The prescribed boundary conditions are

$$\begin{array}{l}y=0:\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}u=u_s\left(x\right),\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}v=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}w=w_s\left(x\right),\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}C=C_s,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}T=T_s,\\ y\to \infty :\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}u\to u_e\left(x\right)=2u_{\infty }\sin \left(\frac{x}{R}\right),\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}w\to w_e=\hspace{0.17em}{w}_{\infty }\cos \left(\mathsf{\theta }\right),\hspace{0.17em}\hspace{0.17em}C\to C_{\infty }\hspace{0.17em},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}T\to T_{\infty }\end{array}\}$$

(6)

The nonsimilar transformations and the resulting variables are

$$\begin{array}{l}\xi ={\displaystyle \underset{0}{\overset{x}{\int }}\frac{u_e}{u_{\infty }} d\left(\frac{x}{R}\right)}\hspace{0.17em},\mathsf{\eta }={\left(\frac{Re}{2\xi }\right)}^{\frac{1}{2}}\left(\frac{u_e}{u_{\infty }}\right)\frac{y}{R},u_{\infty }=w_{\infty }\sin \left(\mathsf{\theta }\right),\\ w=w_e\cdot S=w_{\infty }\cos \left(\mathsf{\theta }\right)\cdot S,G=\frac{T-T_{\infty }}{T_s-T_{\infty }},H=\frac{C-C_{\infty }}{C_s-C_{\infty }},\\ \mathsf{\psi }\left(x,y\right)=u_{\infty }R{\left(\frac{2\mathsf{\xi }}{Re}\right)}^{\frac{1}{2}}f\left(\mathsf{\xi },\mathsf{\eta }\right),u=\frac{\partial \mathsf{\psi }}{\partial y},v=-\frac{\partial \mathsf{\psi }}{\partial x},Re=\frac{\mathsf{\rho }{u}_{\infty }R}{\mathsf{\mu }}\end{array}\}$$

(7)

In view of (7), we have

$$\begin{array}{l}u=2u_{\infty }\sin \left(\frac{x}{R}\right)F\\ v=u_{\infty }\left\{-\frac{y}{R}\left(2\cos \left(\frac{x}{R}\right)-\frac{{\sin }^2\left(\frac{x}{R}\right)}{1-\cos \left(\frac{x}{R}\right)}\right)F-{\left(\frac{2}{\mathrm{Re}\cdot \mathsf{\xi }}\right)}^{\frac{1}{2}}f\cdot \sin \left(\frac{x}{R}\right)-2{\left(\frac{2\mathsf{\xi }}{\mathrm{Re}}\right)}^{\frac{1}{2}}\cdot \sin \left(\frac{x}{R}\right)f_{\mathsf{\xi }}\right\}\end{array}$$

Now, utilizing the transformations (7), Equations (2)–(5) can be written as

$$\begin{array}{l}{F}_{\mathsf{\eta }\mathsf{\eta }}+f{F}_{\mathsf{\eta }}+\mathsf{\beta }\left(\mathsf{\xi }\right)\left(1-F^2\right)+2\mathsf{\xi }\left(F_{\mathsf{\eta }}{f}_{\mathsf{\xi }}-F{F}_{\mathsf{\xi }}\right)+\\ +Ri\hspace{0.17em}s\left(\mathsf{\xi }\right)\left\{\left(1+{\mathsf{\beta }}_{T}G\right)G+\left(1+{\mathsf{\beta }}_{c}H\right)NcH\right\}\sin \left(\mathsf{\theta }\right)=0\end{array}$$

(8)

$$S_{\mathsf{\eta }\mathsf{\eta }}+f\cdot S_{\mathsf{\eta }}+2\mathsf{\xi }\left(S_{\mathsf{\eta }}{f}_{\mathsf{\xi }}-F{S}_{\mathsf{\xi }}\right)+Rip\left(\mathsf{\xi }\right)\left\{\left(1+{\mathsf{\beta }}_{T}G\right)G+\left(1+{\mathsf{\beta }}_{c}H\right)NcH\right\}\sin \left(\mathsf{\theta }\right)=0$$

(9)

$$G_{\mathsf{\eta }\mathsf{\eta }}+Pr\cdot f\cdot G_{\mathsf{\eta }}+2\mathsf{\xi }\left(G_{\mathsf{\eta }}{f}_{\mathsf{\xi }}-F{G}_{\mathsf{\xi }}\right)Pr=0$$

(10)

$$H_{\mathsf{\eta }\mathsf{\eta }}+Sc\cdot f\cdot H_{\mathsf{\eta }}+2\mathsf{\xi }\left(H_{\mathsf{\eta }}{f}_{\mathsf{\xi }}-F{H}_{\mathsf{\xi }}\right)Sc=0$$

(11)

The employed boundary conditions are given below

$$\begin{array}{l}\mathrm{at}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }=0:\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}F={\mathsf{\epsilon }}_1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}S={\mathsf{\epsilon }}_2,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}G=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}H=1,\\ \mathrm{as}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }\to \infty ,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}F=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}S=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}G=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}H=0\end{array}\}$$

(12)

The nondimensional characteristics arising in this investigation are

$$\begin{array}{l}\mathsf{\beta }\left(\mathsf{\xi }\right)=\frac{\mathsf{\xi }\cos \left(\overline{x}\right)}{{\sin }^2\left(\overline{x}\right)}, Ri=\frac{g{\mathsf{\beta }}_1\left(T_s-T_{\infty }\right)R}{u_{\infty }^2}, s\left(\mathsf{\xi }\right)=\frac{\mathsf{\xi }}{4{\sin }^3\left(\overline{x}\right)}, Nc=\frac{{\mathsf{\beta }}_3\left(C_s-C_{\infty }\right)}{{\mathsf{\beta }}_1\left(T_s-T_{\infty }\right)},\\ p\left(\mathsf{\xi }\right)=\frac{\mathsf{\xi }}{2{\sin }^2\left(\overline{x}\right)}, {\mathsf{\beta }}_T=\frac{{\mathsf{\beta }}_2\left(T_s-T_{\infty }\right)}{{\mathsf{\beta }}_1}, {\mathsf{\beta }}_C=\frac{{\mathsf{\beta }}_4\left(C_s-C_{\infty }\right)}{{\mathsf{\beta }}_3}, \mathrm{where} \overline{x}=\frac{x}{R}\end{array}$$

(13)

Furthermore, $f\left(\mathsf{\xi },\mathsf{\eta }\right)={\displaystyle \underset{0}{\overset{\mathsf{\eta }}{\int }}Fd\mathsf{\eta }+f_s;\hspace{0.17em}\hspace{0.17em}\mathrm{where}\hspace{0.17em}\hspace{0.17em}{f}_s=0}$ represents an impermeable yawed cylinder.

The velocity distribution is

$$u_e\left(\overline{x}\right)=2u_{\infty }\sin \left(\overline{x}\right),w_e\left(\overline{x}\right)=w_{\infty }\cos \left(\mathsf{\theta }\right)=\mathrm{constant}, \mathrm{where} \overline{x}=\frac{x}{R}$$

(14)

Therefore, $\mathsf{\xi },\hspace{0.17em}\hspace{0.17em}\mathsf{\beta }\left(\mathsf{\xi }\right),\hspace{0.17em}s\left(\mathsf{\xi }\right)$ and $p\left(\mathsf{\xi }\right)$ can be defined in terms of $\overline{x}$

$$\begin{array}{l}\mathsf{\xi }=2\left[1-\cos \left(\overline{x}\right)\right], \mathsf{\beta }\left(\overline{x}\right)=\frac{2\cos \left(\overline{x}\right)}{1+\cos \left(\overline{x}\right)},\\ s\left(\overline{x}\right)=\frac{1}{2\sin \left(\overline{x}\right)\left[1+\cos \left(\overline{x}\right)\right]}, p\left(\overline{x}\right)=\frac{1}{1+\cos \left(\overline{x}\right)}\end{array}$$

(15)

Using the following relation, equations are expressed in $\overline{x}$ instead of ξ

$$\mathsf{\xi }\frac{\partial }{\partial \mathsf{\xi }}=B\left(\overline{x}\right)\frac{\partial }{\partial \overline{x}}$$

(16)

$$\mathrm{where} B\left(\overline{x}\right)=\tan \left(\frac{\overline{x}}{2}\right)$$

(17)

In view of Equations (16) and (17), Equations (8)–(11) become as follows

$$\begin{array}{l}{F}_{\mathsf{\eta }\mathsf{\eta }}+f{F}_{\mathsf{\eta }}+\mathsf{\beta }\left(\overline{x}\right)\left(1-F^2\right)+2B\left(\overline{x}\right)\left(F_{\mathsf{\eta }}{f}_{\overline{x}}-F{F}_{\overline{x}}\right)+\\ +Ri\cdot \hspace{0.17em}s\left(\overline{x}\right)\left\{\left(1+{\mathsf{\beta }}_{T}G\right)G+\left(1+{\mathsf{\beta }}_{c}H\right)NcH\right\}\sin \left(\mathsf{\theta }\right)=0\end{array}$$

(18)

$$S_{\mathsf{\eta }\mathsf{\eta }}+f\cdot S_{\mathsf{\eta }}+2B\left(\overline{x}\right)\left(S_{\mathsf{\eta }}{f}_{\overline{x}}-F{S}_{\overline{x}}\right)+Ri\cdot p\left(\overline{x}\right)\left\{\left(1+{\mathsf{\beta }}_{T}G\right)G+\left(1+{\mathsf{\beta }}_{c}H\right)NcH\right\}\sin \left(\mathsf{\theta }\right)=0$$

(19)

$$G_{\mathsf{\eta }\mathsf{\eta }}+Pr\cdot f\cdot G_{\mathsf{\eta }}+2B\left(\overline{x}\right)\left(G_{\mathsf{\eta }}{f}_{\overline{x}}-F{G}_{\overline{x}}\right)Pr=0$$

(20)

$$H_{\mathsf{\eta }\mathsf{\eta }}+Sc\cdot f\cdot H_{\mathsf{\eta }}+2B\left(\overline{x}\right)\left(H_{\mathsf{\eta }}{f}_{\overline{x}}-F{H}_{\overline{x}}\right)Sc=0$$

(21)

The relevant boundary conditions are

$$\begin{array}{l}\mathrm{at}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }=0\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}F={\mathsf{\epsilon }}_1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}S={\mathsf{\epsilon }}_2,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}G=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}H=1,\\ \mathrm{as}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }\to \infty \hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}F=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}S=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}G=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}H=0\end{array}\}$$

(22)

Skin friction coefficient along chordwise direction can be defined as follows $C_f=\frac{2{\left(\mathsf{\mu }\frac{\partial u}{\partial y}\right)}_w}{\mathsf{\rho }{w}_{\infty }^2}.$ Taking into account Equation (7) one can define this skin friction parameter as $C_f=\frac{2\mathsf{\mu }\left(4\frac{u_{\infty }}{R}{\sin }^2\left(\overline{x}\right){\left(\frac{Re}{2\mathsf{\xi }}\right)}^{0.5}{F}_{\mathsf{\eta }}\left(\overline{x},0\right)\right){\sin }^2\left(\mathsf{\theta }\right)}{\mathsf{\rho }{u}_{\infty }^2},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\left(u_{\infty }=w_{\infty }\sin \left(\mathsf{\theta }\right)\right)$ and as a result

$$\begin{array}{l}{C}_f=\frac{4\left(1+\cos \left(\overline{x}\right)\right){\left(1-\cos \left(\overline{x}\right)\right)}^{0.5}{\sin }^2\left(\mathsf{\theta }\right)F_{\mathsf{\eta }}\left(\overline{x},0\right)}{R{e}^{0.5}},\\ \mathrm{i}.\mathrm{e}., R{e}^{0.5}{C}_f=4\left(1+\cos \left(\overline{x}\right)\right){\left(1-\cos \left(\overline{x}\right)\right)}^{0.5}{\sin }^2\left(\mathsf{\theta }\right)F_{\mathsf{\eta }}\left(\overline{x},0\right)\end{array}$$

(23)

Skin friction coefficient along spanwise direction can be defined as follows ${\overline{C}}_f=\frac{2{\left(\mathsf{\mu }\frac{\partial w}{\partial y}\right)}_w}{\mathsf{\rho }{w}_{\infty }^2}.$ Taking into account Equation (7) one can define this skin friction parameter as ${\overline{C}}_f=\frac{2\mathsf{\mu }\left(2w_e\frac{\sin \left(\overline{x}\right)}{R}{\left(\frac{Re}{2\mathsf{\xi }}\right)}^{1/2}{S}_{\mathsf{\eta }}\left(\overline{x},0\right)\right){\sin }^2\left(\mathsf{\theta }\right)}{\mathsf{\rho }{u}_{\infty }^2},\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\left(w_e=w_{\infty }\cos \left(\mathsf{\theta }\right)\right)$ and as a result

$${\overline{C}}_f=\frac{2^{1.5}\cos \left(\frac{\overline{x}}{2}\right)\cos \left(\mathsf{\theta }\right)\sin \left(\mathsf{\theta }\right)S_{\mathsf{\eta }}\left(\overline{x},0\right)}{R{e}^{0.5}}, \mathrm{i}.\mathrm{e}., R{e}^{0.5}{\overline{C}}_f=2^{1.5}\cos \left(\frac{\overline{x}}{2}\right)\cos \left(\mathsf{\theta }\right)\sin \left(\mathsf{\theta }\right)S_{\mathsf{\eta }}\left(\overline{x},0\right)$$

(24)

Nusselt number reflecting the heat transfer rate can be defined as $Nu=-\frac{R{\left(\frac{\partial T}{\partial y}\right)}_w}{\left(T_w-T_{\infty }\right)}.$ Taking into account Equation (7), one can find $Nu=-\frac{R\left(\left(T_w-T_{\infty }\right)\frac{2\sin \left(\overline{x}\right)}{R}{\left(\frac{Re}{2\mathsf{\xi }}\right)}^{0.5}{G}_{\mathsf{\eta }}\left(\overline{x},0\right)\right)}{\left(T_w-T_{\infty }\right)}$ and $Nu=-2\sin \left(\overline{x}\right){\left(\frac{Re}{2\mathsf{\xi }}\right)}^{0.5}{G}_{\mathsf{\eta }}\left(\overline{x},0\right)$. Using Equation (15) for ξ one can find ${Re}^{-0.5}Nu=-\frac{\sin \left(\overline{x}\right)}{{\left(1-\cos \left(\overline{x}\right)\right)}^{1/2}}{G}_{\mathsf{\eta }}\left(\overline{x},0\right)$ and ${Re}^{-0.5}Nu=-\left(\frac{2\sin \left(\frac{\overline{x}}{2}\right)\cos \left(\frac{\overline{x}}{2}\right)}{\sqrt{2}\cdot \sin \left(\frac{\overline{x}}{2}\right)}{G}_{\mathsf{\eta }}\left(\overline{x},0\right)\right)$, where $\left(1-\cos \left(\overline{x}\right)\right)=2\cdot {\sin }^2\left(\frac{\overline{x}}{2}\right)$. As a result

$$Nu=\frac{-\sqrt{2}\cos \left(\frac{\overline{x}}{2}\right)G_{\mathsf{\eta }}\left(\overline{x},0\right)}{R{e}^{-0.5}}, \mathrm{i}.\mathrm{e}., R{e}^{-0.5}Nu=-\sqrt{2}\cos \left(\frac{\overline{x}}{2}\right)G_{\mathsf{\eta }}\left(\overline{x},0\right)$$

(25)

Sherwood number reflecting the mass transfer rate can be defined as $Sh=-\frac{R{\left(\frac{\partial C}{\partial y}\right)}_w}{\left(C_w-C_{\infty }\right)}$. Taking into account Equation (7) one can find $Sh=-\frac{R\left(\left(C_w-C_{\infty }\right)\frac{2\sin \left(\overline{x}\right)}{R}{\left(\frac{Re}{2\mathsf{\xi }}\right)}^{0.5}{H}_{\mathsf{\eta }}\left(\overline{x},0\right)\right)}{\left(C_w-C_{\infty }\right)}$ and $Sh=-2\sin \left(\overline{x}\right){\left(\frac{Re}{2\mathsf{\xi }}\right)}^{0.5}{H}_{\mathsf{\eta }}\left(\overline{x},0\right)$. Using Equation (15) for ξ one can find ${Re}^{-0.5}Sh=-\frac{\sin \left(\overline{x}\right)}{{\left(1-\cos \left(\overline{x}\right)\right)}^{0.5}}{H}_{\mathsf{\eta }}\left(\overline{x},0\right)$ and ${Re}^{-0.5}Sh=-\left(\frac{2\sin \left(\frac{\overline{x}}{2}\right)\cos \left(\frac{\overline{x}}{2}\right)}{\sqrt{2}\cdot \sin \left(\frac{\overline{x}}{2}\right)}{H}_{\mathsf{\eta }}\left(\overline{x},0\right)\right)$, where $\left(1-\cos \left(\overline{x}\right)\right)=2\cdot {\sin }^2\left(\frac{\overline{x}}{2}\right)$. As a result

$$Sh=\frac{-\sqrt{2}\cos \left(\frac{\overline{x}}{2}\right)H_{\mathsf{\eta }}\left(\overline{x},0\right)}{R{e}^{-0.5}}, \mathrm{i}.\mathrm{e}., R{e}^{-0.5}Sh=-\sqrt{2}\cos \left(\frac{\overline{x}}{2}\right)H_{\mathsf{\eta }}\left(\overline{x},0\right)$$

(26)

3. Solution Technique

Equations (18)–(21) are linearized using quasilinearization technique as below

$$F_{\mathsf{\eta }\mathsf{\eta }}^{i+1}+A_1^i{F}_{\mathsf{\eta }}^{i+1}+A_2^i{F}_{\overline{x}}^{i+1}+A_3^i{F}^{i+1}+A_4^i{G}^{i+1}+A_5^i{H}^{i+1}=A_6^i$$

(27)

$$S_{\mathsf{\eta }\mathsf{\eta }}^{i+1}+B_1^i{S}_{\mathsf{\eta }}^{i+1}+B_2^i{S}_{\overline{x}}^{i+1}+B_3^i{F}^{i+1}+B_4^i{G}^{i+1}+B_5^i{H}^{i+1}=B_6^i$$

(28)

$$G_{\mathsf{\eta }\mathsf{\eta }}^{i+1}+C_1^i{G}_{\mathsf{\eta }}^{i+1}+C_2^i{G}_{\overline{x}}^{i+1}+C_3^i{F}^{i+1}=C_4^i$$

(29)

$$H_{\mathsf{\eta }\mathsf{\eta }}^{i+1}+D_1^i{H}_{\mathsf{\eta }}^{i+1}+D_2^i{H}_{\overline{x}}^{i+1}+D_3^i{F}^{i+1}=D_4^i$$

(30)

where the coefficients at the (i + 1)th iteration are defined employing the known ith iterative parameters. The dimensionless boundary conditions are

$$\begin{array}{l}{F}^{i+1}={\epsilon }_1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{S}^{i+1}={\epsilon }_2,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{G}^{i+1}=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{H}^{i+1}=1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathrm{at}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }=0,\\ F^{i+1}=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{S}^{i+1}=1,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{G}^{i+1}=0,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}{H}^{i+1}=0\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\mathrm{at}\hspace{0.17em}\hspace{0.17em}\mathsf{\eta }={\mathsf{\eta }}_{\infty }\end{array}\}$$

(31)

where the boundary layer edge is denoted by η_∞.

The coefficients in Equations (27)–(30) are given below

$$\begin{array}{l}{A}_1^i=f+2B\left(\overline{x}\right)f_{\overline{x}},A_3^i=-2B\left(\overline{x}\right)F,A_3^i=-2\mathsf{\beta }\left(\overline{x}\right)F-2B\left(\overline{x}\right)F_{\overline{x}},\\ A_4^i=Ri\cdot s\left(\overline{x}\right)\left(1+2{\mathsf{\beta }}_{T}G\right)\sin \left(\mathsf{\theta }\right), A_5^i=Ri\cdot s\left(\overline{x}\right)\left(1+2{\mathsf{\beta }}_{C}H\right)Nc\sin \left(\mathsf{\theta }\right),\\ A_6^i=-2B\left(\overline{x}\right)F{F}_{\overline{x}}-\mathsf{\beta }\left(\overline{x}\right)\left(1+F^2\right)+Ri\cdot s\left(\overline{x}\right)\left({\mathsf{\beta }}_T{G}^2+Nc{\mathsf{\beta }}_C{H}^2\right)\sin \left(\mathsf{\theta }\right),\\ B_1^i=f+2B\left(\overline{x}\right)f_{\overline{x}},B_2^i=-2B\left(\overline{x}\right)F,B_3^i=-2B\left(\overline{x}\right)S_{\overline{x}},\\ B_4^i=Ri\cdot p\left(\overline{x}\right)\left(1+2{\mathsf{\beta }}_{T}G\right)\sin \left(\mathsf{\theta }\right), B_5^i=Ri\cdot p\left(\overline{x}\right)\left(1+2{\mathsf{\beta }}_{C}H\right)Nc\sin \left(\mathsf{\theta }\right),\\ B_6^i=-2B\left(\overline{x}\right)F{S}_{\overline{x}}+Ri\cdot p\left(\overline{x}\right)\left({\mathsf{\beta }}_T{G}^2+Nc{\mathsf{\beta }}_C{H}^2\right)\sin \left(\mathsf{\theta }\right),\\ C_1^i=Pr\left(f+2B\left(\overline{x}\right)f_{\overline{x}}\right),C_2^i=-2PrB\left(\overline{x}\right)F,C_3^i=-2PrB\left(\overline{x}\right)G_{\overline{x}},C_4^i=-2PrB\left(\overline{x}\right)F{G}_{\overline{x}},\\ D_1^i=Sc\left(f+2B\left(\overline{x}\right)f_{\overline{x}}\right),D_2^i=-2ScB\left(\overline{x}\right)F,D_3^i=-2ScB\left(\overline{x}\right)H_{\overline{x}},D_4^i=-2ScB\left(\overline{x}\right)F{H}_{\overline{x}}\end{array}$$

The nonlinear coupled partial differential Equations (18)–(21) under the boundary conditions (22) have been solved numerically using an implicit finite difference scheme in combination with the quasilinearization technique. The quasilinearization technique can be viewed as a generalization of the Newton–Raphson approximation method in functional space. An iterative sequence of linear equations is carefully constructed to approximate the nonlinear Equations (18)–(21) under the boundary conditions (22) achieving quadratic convergence and monotonicity. Applying the quasilinearization technique, the nonlinear coupled partial differential Equations (18)–(21) with boundary conditions (22) are replaced by the sequence of linear ordinary differential equations. Because the method is presented for ordinary differential equations by Inouye and Tate [37] and for partial differential equations in a recent study by Singh and Roy [38], its detailed description is not provided here. At each iteration step, the sequence of linear partial differential Equations (27)–(30) is expressed in difference form using the central difference scheme in the $\overline{x}$-direction and the backward difference scheme in η direction. Thus, in each step, the resulting equations have been then reduced to a system of linear algebraic equations with a block tri-diagonal matrix, which is solved by Varga’s algorithm [39]. To ensure the convergence of the numerical solution to the exact solution, step sizes $\Delta \overline{x}$ and Δη are taken as 0.01 and 0.01. A convergence criterion based on the relative difference between the current and previous iteration values is employed. When the difference reaches 0.0001, the solution is assumed to have converged and the iteration process is terminated, i.e., $\max \left\{\left|{\left(F_{\eta }\right)}_s^{i+1}-{\left(F_{\eta }\right)}_s^i\right|,\left|{\left(G_{\eta }\right)}_s^{i+1}-{\left(G_{\eta }\right)}_s^i\right|,\left|{\left(H_{\eta }\right)}_s^{i+1}-{\left(H_{\eta }\right)}_s^i\right|,\left|{\left(S_{\eta }\right)}_s^{i+1}-{\left(S_{\eta }\right)}_s^i\right|\right\}\le {10}^{-4}.$

4. Results and Discussion

In order to examine the qualitative behavior of fluid flow characteristics over a moving yawed cylinder, comprehensive numerical computation has been carried out for numerous values of nondimensional parameters that characterize the flow, heat and mass transfer. The simulated outcomes are presented graphically. Considering water as the working liquid, the value of Prandtl number Pr = 7 is selected. The realistic Schmidt numbers are Sc = 160, 240 and 340 which correspond to liquid hydrogen, liquid nitrogen and liquid oxygen, respectively. The values of yaw angle (θ), combined convection parameter (Ri) and the ratio of buoyancy force parameter (Nc) are varied in the ranges 0 ≤ θ ≤ π/3, –3 ≤ Ri ≤ 10 and 0 ≤ Nc ≤ 1, respectively. It should be noted that considered values of governing parameters allow one to show the features of liquid behavior and heat and mass transfer.

Further, θ = 0 corresponds to the perfect geometry of a vertical cylinder. Moreover, $\overline{x}=0$ and $\overline{x}\ne 0$ represent the similarity and nonsimilarity cases, respectively. Additionally, ${\mathsf{\epsilon }}_1=u_s/u_e$ and ${\mathsf{\epsilon }}_2=w_s/w_e$ are two velocity ratio parameters along chordwise and spanwise directions, respectively.

4.1. Impacts of Combined Convection Parameter and Velocity Ratio Parameter

Figure 2, Figure 3 and Figure 4 display the influence of combined convection parameter (Ri) and velocity ratio parameters $\left({\mathsf{\epsilon }}_1\hspace{0.17em}\mathrm{and}\hspace{0.17em}{\mathsf{\epsilon }}_2\right)$ on velocity patterns $\left(F\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)\hspace{0.17em}\mathrm{and}\hspace{0.17em}S\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)\right)$ and the coefficients of friction at the surface $\left(R{e}^{1/2}{C}_f\hspace{0.17em}\mathrm{and}\hspace{0.17em}R{e}^{1/2}{\overline{C}}_f\right)$ along chordwise and spanwise directions. Increasing magnitude of combined convection characteristics enhances the velocity pattern and the surface drag coefficient in chordwise and spanwise directions. The positive magnitudes of Ri designate the essential impacts of buoyancy force over inertia force.

Thus, the buoyancy for liquid enhances the liquid motion and simultaneously enhances the liquid’s velocity and respective frictions at the surface. Further, for enhancing values of velocity ratio characteristics (ε₁), the velocity profile enhances, while the surface drag coefficient in the chordwise direction diminishes. Moreover, the same behavior is observed for larger values of velocity ratio parameters (ε₂) in the spanwise direction. The physical reason is that, for ε₁ < 1 and ε₂ < 1, the free stream velocity dominates over surface velocity and causes such variations. Moreover, the combined impacts of assisting buoyancy force because of temperature and concentration gradients along with a rise of velocity ratio parameter acts like a favorable pressure gradient which enhances the liquid flow. For $\overline{x}=0.5$, Pr = 7.0, β_T = 0.1, θ = π/6, Nc = 0.1, ε₁ = ε₂ = 1.5, β_C = 0.1 and Sc = 160, the coefficients of friction along chordwise and spanwise directions at the wall are increased by about 22% and 62%, respectively, as Ri varies from opposing buoyancy flow case (Ri = −1) to assisting buoyancy flow case (Ri = 3).

4.2. Impacts of Nonlinear Convection Parameter and Yaw Angle

Figure 5, Figure 6, Figure 7 and Figure 8 display the influence of velocity ratio characteristics (ε₁) and angle of yaw (θ) on velocity profile $F\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)$ and the coefficients of friction at the surface $\left(R{e}^{0.5}{C}_f\right)$ along the chordwise direction. Moreover, the influence of nonlinear convection parameter (β_T) and angle of yaw (θ) on velocity profile $S\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)$ and the coefficient of friction at the surface $\left({Re}^{0.5}{\overline{C}}_f\right)$ along the spanwise direction has been observed. For enhancing values of velocity ratio parameter (ε₁), the velocity pattern increases, while surface drag coefficient in chordwise direction diminishes. The physical reason is that, for ε₁ < 1, the free stream velocity dominates over surface velocity and causes such variations. Further, for the increasing magnitudes of β_T, the velocity patterns and coefficient of friction at the surface increase along the spanwise direction. The larger magnitudes of β_T characterize the greater variation between the plate and the environmental temperature. Therefore, for larger magnitudes of β_T, the larger temperature variation reasons stronger convection and as a result enhances the liquid velocity and the friction between the cylinder and the liquid.

Further, velocity profiles in chordwise and spanwise directions and the skin friction coefficient in chordwise and spanwise directions both increase with the increasing values of yaw angle. Higher values of yaw angle (i.e., the cylinder inclines more) cause pressure increase in the fluid flow and increase fluid velocity, and the inclination of the cylinder causes the enhancement in the surface friction along chordwise and spanwise directions. Moreover, in Figure 7, at θ = 0, the lines obtained for different values of β_T merge together. For $\overline{x}=0.5$, Pr = 7.0, Ri = 10.0, β_T = 0.5, Nc = 0.1, ε₁ = 0.8, ε₂ = 0.5, β_C = 0.1 and Sc = 160, the coefficient of friction along spanwise directions at the wall increases by about 50%, as θ varies from π/6 to π/3.

4.3. Impacts of Yaw Angle and Schmidt Number

Figure 9 and Figure 10 display the impact of Schmidt number (Sc) as well as angle of yaw (θ) on concentration profile $H\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)$ and the mass transfer rate $\left({Re}^{-0.5}Sh\right)$, respectively. It is observed that the concentration profile reduces while the Sherwood number enhances, for the larger Schmidt numbers and angles of yaw. Mass diffusivity reduces for the higher values of Sc. As Sc enhances, the thickness of the concentration boundary layer reduces, and as a result the species concentration profile diminishes. Accordingly, the Sherwood number increases. Further, we observe that there is no deviation for θ = 0 which represents the vertical cylinder, as compared to other yaw angles. For $\overline{x}=0.5$, Ri = 10.0, β_T = 0.1, Nc = 0.1, ε₁ = 0.8, ε₂ = 0.5, β_C = 0.1 and θ = π/6, the Sherwood number enhances by about 27%, when Sc value is reduced from 340 to 160.

4.4. Impact of Yaw Angle and Combined Convection Characteristics

Figure 11 and Figure 12 depict the influence of combined convection parameter (Ri) and angle of yaw (θ) on dimensionless temperature pattern $G\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)$ and heat transfer rate $\left({Re}^{-0.5}Nu\right)$. Temperature pattern and heat transfer rate are reduced for an enhancing magnitudes of combined convection parameter. The magnitudes Ri > 0 illustrate the essential influence of buoyancy force compared to the inertia force. Thus, a rise of the buoyancy impacts illustrates more essential fluid flow that characterizes a reduction of the liquid temperature and an appearance of the cool liquid close to the cylinder’s border. This in turn enhances the energy transport strength from the cylinder’s border to the liquid. Further, as yaw angle enhances temperature profile and heat transfer rate, for enhancing magnitudes of combined convection parameter, we observe the curves obtained by varying the yaw angle; the curves overlap into a single curve with irrespective values of combined convection parameter.

Further, we observe that there is no deviation for θ = 0 which represents the vertical cylinder, as compared to other yaw angles. Moreover, at $\overline{x}=0.5$ and Ri = 10.0, as yaw angle enhances from θ = π/6 to θ = π/3, Nusselt number increases approximately about 7%.

Velocity profile $F\left(\overline{x},\hspace{0.17em}\mathsf{\eta }\right)$ is compared with the already existing results of Eswara and Nath [40] for a particular case by assuming A = 0, Ec = 0, θ = 0. The results are found in an excellent agreement and the comparisons are shown in Figure 13.

Table 1. Surface drag coefficient along chordwise direction $\left({Re}^{0.5}{C}_f\right)$, heat transfer rate $\left({Re}^{-0.5}Nu\right)$ and mass transfer rate $\left({Re}^{-0.5}Sh\right)$ values, for different values of yaw angle (θ = 0 to θ = π/3) when Ri = 10.0, β_T = 0.1, Nc = 0.1, ε₁ = 0.8, ε₂ = 0.5, β_C = 0.1 and Sc = 160.

Yaw Angle (θ)	$\left|{\mathit{R}\mathit{e}}^{0.5}{\mathit{C}}_{\mathit{f}}\right|$	$\left|{\mathit{R}\mathit{e}}^{-0.5}\mathit{N}\mathit{u}\right|$	$\left|{\mathit{R}\mathit{e}}^{-0.5}\mathit{S}\mathit{h}\right|$
0	0	0.94430	4.41621
π/12	0.12140	0.94281	4.42190
π/6	0.68785	0.95536	4.43383
π/3	3.10383	0.96857	4.45684

demonstrates the variations of friction parameter at the yawed cylinder along chordwise direction, Nusselt and Sherwood numbers, for various magnitudes of yaw angle. From the Table 1, it is found that the surface drag coefficient, Nusselt and Sherwood numbers are raised for increasing yaw angles. The friction parameter, rates of heat transfer and mass transfer are enhanced by about 31%, 2% and 1%, respectively, at $\overline{x}=0.5$, as yaw angle varies from π/12 to π/6.

5. Conclusions

This research considers the double-diffusive combined convection around a moved yawed cylinder. In that, the influence of varying strength of yaw angle, mixed convection and nonlinear convection characteristics on the velocity patterns and skin friction coefficients in chordwise and spanwise directions, nondimensional temperature, concentration profile, mass and heat transfer rates are analyzed employing different graphs. Taking into account this detailed analysis, the obtained outcomes are detailed as follows:

−

Velocity profiles can be enhanced, while the coefficients of friction at the surface diminish, for increasing values of velocity ratio parameters in spanwise and chordwise directions.

−

For enhancing magnitudes of nonlinear convection coefficient, the velocity profile and the skin friction parameter in spanwise direction are increased.

−

Concentration profile diminishes, while the Sherwood number enhances, for increasing values of Schmidt number and yaw angle.

−

Velocity profiles in spanwise and chordwise directions and skin friction coefficient at the border in chordwise and spanwise directions are enhanced with growing values of yaw angle.

−

Increasing magnitude of combined convection characteristics, enhancing the velocity profiles and surface drag coefficient in spanwise and chordwise directions.