Search

A mixed convection flow and heat transfer of pseudo-plastic power law nanofluid past a stretching vertical plate is investigated. Three types of nano-particles, such as copper (Cu), aluminum oxide (Al 2 O 3 ) and titanium oxide (TiO 2 ), are considered. The generalized Fourier law proposed by Zheng for varying thermal conductivity of nanofluids, which is dependent on the power law of velocity gradient as well as the nano-particles property, is taken into account. Dual solutions are obtained by Bvp4c with Matlab. The stability of the solution also is discussed by introducing the unsteady governing equations. Furthermore, a new interesting phenomenon is found: the local Nusselt number do not maintain the similar characteristics of Newtonian fluid near the point where the velocity ratio is equal to 0.5. [ABSTRACT FROM AUTHOR]

The modified model of magneto-hydrodynamic (MHD) flow and heat transfer of power-law fluid over a stretching permeable surface with prescribed heat flux in the presence of heat generation or absorption is investigated. The thermal conductivity is assumed to be dependent on the power law of velocity gradient. The general boundary conditions are given and Lie-group method is used to determine the appropriate boundary conditions suitable for the similarity solution. By using the achieved similarity transformations, the similarity equations are obtained and solved numerically by Bvp4c with Matlab. Graphical solutions for the dimensionless velocity and temperature as well as skin-friction coefficient and Nusselt number for various parameters are presented and discussed. [ABSTRACT FROM AUTHOR]

In this paper, Haar wavelet operational matrix(HWOM) is proposed to solve initial-boundary value problems for a class of time-space fractional partial differential equations of Caputo sense with variable coefficients in both time and space (1) ∑ i = 1 n θ i ∂ γ i u ( x , t ) ∂ t γ i = v ( x , t ) ∂ α u ( x , t ) ∂ x α + d ( x , t ) ∂ β u ( x , t ) ∂ x β + q ( x , t ) , 0 < x < 1 , 0 < t ≤ 1 , as an extension of Rehman and Khan's (2013) work. We obtain a matrix L instead of Q α in Rehman and Khan (2013). when dealing with boundary conditions. By utilizing the operational matrix of fractional integration and Hadamard product,we made an improvement of algorithm to deal with time-space coefficients and gave the error analysis of the HWOM for space-time dimensions. Some numerical results are paralleled with exact solutions to show the efficiency and precision of the presented technique. [ABSTRACT FROM AUTHOR]

In this paper, we examine the laminar film condensation of pseudo-plastic non-Newtonian fluids with variable thermal conductivity on an isothermal vertical plate. The thermal conductivity is assumed to be power-law-dependent on the velocity gradient. The dual similar solutions, which are influenced by the parameters of the power law number n and the thickness η δ of the condensation film, are obtained numerically by Runge–Kutta method coupled with shooting method. More attention is paid to discuss the first branch of the solutions with physical meaning. Especially the effects of above both parameters on the velocity and temperature distribution, condensation mass flow rate and the local Nusselt number are analyzed in detail. [ABSTRACT FROM AUTHOR]

In this paper, the flow of a micropolar fluid through a porous channel with deforming walls is investigated. The Lie group method is applied to determine symmetry reductions of partial differential equations. Compared with previous work, the similar variable for the microrotation velocity also is obtained, which is an extension of the case corresponding to the Newtonian fluids (Boutros et al., 2007). This problem we study here can be considered as an extension of the model discussed in Boutros et al. (2007). The effect of the permeation Reynolds number Re , the micropolar parameter K and the expansion ratio α are studied numerically and the results are plotted. Furthermore, the multiple solutions also can be found for some values of parameters. [ABSTRACT FROM AUTHOR]

A time-distributed order fractional continuity model is proposed to simulate the rotating electro-osmotic slip flow of generalized Maxwell fluids in an alternating electric field. The model fully exhibits the wider memory characteristics of viscoelastic fluid. Complex momentum equation is discretized by L1 and L2 algorithms and verified by constructing analytical solutions method. Due to the stronger memory of generalized Maxwell fluids, reverse flow is restrained, and the time for oscillation reaching the steady state is postponed. Moreover, with the increase of Debye–Hückle coefficient the velocity will increase rapidly. • The time-distributed fractional order model is introduced to describe the characteristics of viscoelastic fluid. • The effects of alternating current frequency and Debye–Hückle coefficient are considered. • The reverse flow is restrained due to the existence of viscoelastic characteristics. • The time to reach constant amplitude oscillation is later compared with Newtonian fluid. [ABSTRACT FROM AUTHOR]

This study numerically simulates the two-dimensional flow of Oldroyd-B fluid around an isosceles right-angled triangular cylinder with five orientations. The log-conformation reformulation is employed to stabilize the numerical simulations. By adjusting the triangular orientation angle (θ), three types of fluids development process can be observed: from steady to vortex shedding at θ = 0 and π , keeping the vortex shedding at θ = π 4 and 3 π 4 , and from vortex shedding to steady state at θ = π 2 . When the triangular cylinder faces the incoming stream with the inclined plane, the elastic effect acting on the cylinder is strong, otherwise it is weak. For θ = π 2 , the effects of the viscosity ratio (β), the Reynolds number (R e), and the Weissenberg number (W i) are further investigated. When the elasticity is reduced by changing the viscosity ratio (β) that ranged from 0 to 0.9, the final flow state will transition from stable to vortex shedding state, which indicates the restraining effect of elasticity on wake instability. In the high elastic Oldroyd-B fluid, the critical Reynolds number for vortex shedding is about 110 for W i = 1. Besides, the Weissenberg numbers (W i) ranged from 0.25 to 8 are discussed at R e = 100. With the increase of W i , four different flow states of the wake are observed: periodic vortex shedding at low Weissenberg number W i = 0. 25 , stabilizing for W i ranged from 0.5 to 1, semi-periodic strong vortex shedding for W i is about 2, and chaos when W i ≥ 4. The results indicate that excessively strong elastic effects may also lead to unstable flows. Finally, the flow states corresponding to each Reynolds number and Weissenberg number in a certain range (90 ≤ R e ≤ 120 and 0. 125 ≤ W i ≤ 4) are given in this study. • The strong elasticity mainly acts on the hypotenuse of the leading edge of the triangular cylinder. • The vortex shedding can be inhibited by the elasticity. • Inertial elastic instability is observed at high Weissenberg numbers. [ABSTRACT FROM AUTHOR]

The finite difference schemes for equations depicting the fluid flow, heat and mass transfer of viscoelastic fluid in porous media based on fractional constitutive model are analyzed, when Soret and Dufour effects are taken into account. Different from the conventional fractional partial differential equations, the governing momentum equation of this system not only contains convection terms, but also a term: D t β 0 (∂ 2 u / ∂ y 2) , which are great challenges when we do the stability analysis. In addition, the temperature equation and concentration equation are coupled with each other, which increases the difficulty of numerical analysis. By using the energy method, the stability and convergence of the finite difference schemes employed for the governing equations are obtained. The accuracy of the scheme for momentum equation is O (τ + h x 2 + h y 2) , while the one for temperature and concentration equations is O (τ 2 + h x 2 + h y). Finally, numerical examples are presented to verify methods and demonstrate the analysis effectiveness. [ABSTRACT FROM AUTHOR]

This work investigates unsteady electro-osmotic flow of a fractional second-order viscoelastic fluid in a micro-rectangle channel with considering depletion effects. Due to the existence of depletion, there are depletion layer near the plate and bulk flow adjacent to depletion layer for the flow in the whole domain. In order to express the boundary conditions at the interface, we introduce Maxwell stress to describe the interaction between Newtonian fluid and viscoelastic fluid with fractional derivative. To ensure the validity of the semi-analytical solution obtained by Laplace transform, the numerical solutions are acquired by finite difference method and Grünwald-Letnikov approximation and both of them have a better agreement. The increasing viscosity of viscoelastic fluid leads to the increasing driving difficulty of fluid in non-depleted layer and the decreasing velocity amplitude at steady state. Moreover, the increasing interfacial zeta potential difference can promote the fluid flow, and the change of interface charge density makes the velocity at the interface change dramatically. [ABSTRACT FROM AUTHOR]

This paper mainly analyzes the heat transfer of Carreau fluid over an infinite rotating disk, and two models with variable thermal conductivity are considered. Firstly, the governing partial differential equations (PDEs) are transformed into ordinary differential equations (ODEs) by similarity transformation. Then, the boundary value problem is solved by improved bvp4c method, which reduced sensitivity to the initial values. For shear-thinning and shear-thickening fluids, the effects of Carreau fluid index n and Prandtl number Pr on velocity and temperature fields are shown and analyzed when 0.5 ≤ n ≤ 1.5 , 1 ≤ P r ≤ 10. Furthermore, the thermal conductivity is computed under two cases. [ABSTRACT FROM AUTHOR]

This paper is concerned with multiple solutions of a singular nonlinear boundary value problem (BVP) on the interval [ 0 , 1 ] , which arises in a study of the laminar flow in a porous pipe with an expanding or contracting wall. For the singular nonlinear BVP, the correct boundary conditions are derived to guarantee that its linearization has a unique smooth solution. Then a numerical technique is proposed to find all possible multiple solutions. For the suction driven pipe flow with the expanding wall (e.g. α = 2 ), we find a new solution numerically and classify it as a type VI solution. The computed results agree well with what can be obtained by the bifurcation package AUTO. In addition, we also construct asymptotic solutions for a few cases of parameters, which agree well with numerical solutions. These serve as validations of our numerical results. Thus we believe that the numerical technique designed in the paper is reliable, and may be further applied to solve a variety of nonlinear equations that arise from other flow problems. [ABSTRACT FROM AUTHOR]

• A Caputo fractional-order COVID-19 model with generalized incidence rate is established, and some dynamic behaviors are analyzed. • Identify the parameters of the model according to relevant statistical data, which verifies the effectiveness of the model for depicting the COVID-19 pandemic. • Based on the relevant statistical data, the corresponding fractional optimal control problem is proposed and analyzed. • The numerical format of Riemann–Liouville adjoint system with transversality condition is given, and the fractional optimal control problem of the COVID-19 model is solved numerically. In view of the spread of corona virus disease 2019 (COVID-19), this paper proposes a fractional-order generalized SEIR model. The non-negativity of the solution of the model is discussed. Based on the established threshold R 0 , the existence of the disease-free equilibrium and endemic equilibrium is analyzed. Then, sufficient conditions are established to ensure the local asymptotic stability of the equilibria. The parameters of the model are identified based on the statistical data of COVID-19 cases. Furthermore, the validity of the model for describing the COVID-19 outbreak is verified. Meanwhile, the accuracy of the relevant theoretical results are also verified. Considering the relevant strategies of COVID-19 prevention and control, the fractional optimal control problem (FOCP) is proposed. Numerical schemes for Riemann–Liouville (R–L) fractional-order adjoint system with transversal conditions is presented. Based on the relevant statistical data, the corresponding FOCP is numerically solved, and the control effect of the COVID-19 outbreak under the optimal control strategy is discussed. [ABSTRACT FROM AUTHOR]

This paper is concerned with asymptotic solutions of a nonlinear boundary value problem (BVP), which arises in a study of laminar flow in a uniformly porous channel with retractable walls and an applied transverse magnetic field. For different ranges of the control parameters (i.e. α , Re and M ) arising in the BVP, four cases are considered using different singular perturbation methods. For the first case, unlike those in the existing literature, we make use of the Lighthill method and successfully construct an asymptotic solution with high-order derivatives at the center of the channel. For the second case, under large suction we consider M 2 = O (1) and M 2 = O ( Re ), respectively, which will further extend the applying range of asymptotic solutions. In other cases, asymptotic solutions with a boundary layer are successfully constructed. In addition, numerical solutions presented for each case agree well with asymptotic solutions, which illustrates that the asymptotic solutions constructed in this paper are more reliable. Finally, the influences of some parameters on flow field are discussed to develop a better understanding of the flow problem. [ABSTRACT FROM AUTHOR]

In order to study the mass separation or mixture of EOF for biological or chemical fluids with viscoelastic characteristics, the EOF of two layers consisting of Newtonian fluid and fractional second-order fluid in a rotating frame in parallel microchannel is investigated. The numerical solutions are obtained by finite difference method. Two interesting phenomena can be found. One of them is that reverse flow can be found near the interface as the rotation velocity is big enough. Another is that the streamwise velocity will increase firstly and decrease subsequently before the steady flow reaches on account of the existence of angular velocity. [ABSTRACT FROM AUTHOR]