Muhammad Faisal

Assistant Professor of Mathematics, Department of Mathematics, Faculty of Science
The University of Azad Jammu and Kashmir


Muhammad Faisal is a Faculty Member in the Department of Mathematics, Faculty of Science, Azad Jammu & Kashmir University, Muzaffarabad, Pakistan. He received his PhD, MPhil and MSc degrees from Department of Mathematics UAJK in 2021, 2013 and 2010, respectively. His research is in the field of computational fluid dynamics, thermal engineering, engineering mathematics and physics, nanocomposites, computational methods, nanofluids, hybrid nanofluids, etc.


PhD in Mathematics
MPhil in Mathematics
MSc in Mathematics
Specialization: Applied Mathematics (Fluid Mechanics and Numerical Methods)


Computational Mathematics, Fluid Flow and Transfer Processes, Materials Science, General Physics and Astronomy


Dynamics of Nanofluids: Modeling and Simulation

To gain an understanding of the heat transfer performance  Applications of extending surfaces in nanotechnology  Stability analysis of the solution  Numerical Inspections of Newtonian and non-Newtonian nanofluid flows

Applications Invited

Dynamics of Hybrid nanofluids: Modeling and Simulation

1. Inspection of heat transport in hybrid nanofluids. 2. Discussion on thermo-physical properties of various nanoparticles. 3. Thermal conductivity of heterogeneous two-component systems. 4. Unsteady flows towards stretching devices. 5. Numerical Simulation via Keller-Box method.

Applications Invited

Dynamics of Ternary Hybrid nanofluids: Modeling and Simulation

1. To study the aspects of random motion and thermo-migration of nanoparticles. 2. Discussion on the shapes of nanoparticles. 3. Numerical Simulation via Keller-Box method. 4. Unsteady flows towards stretching devices.

Applications Invited

Scopus Publications

Scopus Publications

  • On convective heat and mass transport of radiative double diffusive Williamson hybrid nanofluid by a Riga surface
    Muhammad Faisal, F. Mabood, K. K. Asogwa, and I. A. Badruddin

    World Scientific Pub Co Pte Ltd
    Convective heat and mass transport of radiative Williamson hybrid [Formula: see text] nanofluid (NF) by a Riga surface with the novel features of Cattaneo–Christov double-diffusion has been investigated. Thermal contributions of internal heat mechanism and Arrhenius energy in Darcy–Forchheimer medium have also been incorporated in the modeling. Mathematical modeling has been completed by using suitable mathematical expressions for thermophysical features of hybrid nanofluid (HNF). Transport partial differential equations (PDEs) have been transformed into ordinary differential equations (ODEs) by means of similarity variables. Numerical approximation of the transformed system has been obtained by using shooting-based Runge–Kutta–Fehlberg approach. Results have been presented through various graphs and discussed physically in detail. Solution is validated for limited cases. Concentration of the hybrid mixture is reduced for progressive concentration-relaxation parameter. Temperature is alleviated for developing thermal-relaxation parameter. Nusselt number is observed to be higher for Williamson HNF than simple ordinary NF.

  • Analysis of entropy generation and activation energy in unsteady double diffusive flow over a stretching cylinder
    Muhammad Faisal, Iftikhar Ahmad, and Muhammad Awais Awan

    AbstractEntropy generation with activation energy has garnered worldwide attention from researchers due to its extensive applications in thermodynamic design, chemical engineering, and optimization processes. One of the primary reasons for focusing on entropy is the global energy crisis, driven by massive consumption against limited resources. This study elucidates the double diffusion process in the dynamics of a viscous fluid over a stretching cylinder, incorporating considerations of entropy generation and activation energy. The non‐Fourier heat flux model is employed to describe the thermodynamics of the thermal system, while non‐Fick's law is used to elucidate the mass diffusion process. The unsteady phases of axisymmetric flow of the viscous fluid are thoroughly discussed. Transport equations in cylindrical configuration are transformed into ordinary differential equations, and an efficient mathematical method (i.e., homotopy analysis method) is applied to solve the transformed system. The graphical examination reveals the impact of flow parameters on velocity configuration, entropy configuration, temperature configurations, and concentration configurations. Additionally, a comparative benchmark is established with the results from previously published work for authentication and validation purposes. It is noted that concentration and temperature configurations are directly related to the unsteady stretching parameter. Entropy generation exhibits an inverse relationship with the unsteady stretching parameter when close to the stretching cylinder, whereas it shows a direct relationship when the fluid is flowing away from the surface.

  • Darcy–Forchheimer dynamics of hybrid nanofluid due to a porous Riga surface capitalizing Cattaneo–Christov theory
    Muhammad Faisal, Kanayo Kenneth Asogwa, Fazle Mabood, and Irfan Anjum Badruddin

    Informa UK Limited

  • Numerical analysis of heat and mass transfer with viscous dissipation, Joule dissipation, and activation energy
    G. Dharmaiah, B. Shankar Goud, Nehad Ali Shah, and Muhammad Faisal

    Informa UK Limited

  • Entropy optimization in squeezed nanofluidic dissipative transport of radiative water conveying aluminum alloys
    Muhammad Faisal, Kanayo Kenneth Asogwa, F. Mabood, and I. A. Badruddin

    World Scientific Pub Co Pte Ltd
    In this paper, squeezing transport of radiative water conveying aluminum alloys (i.e., AA7072 and AA7075) mobilized by entropy generation and dissipative energy is analyzed. Problem is formulated in a rotating frame with the consideration of magnetohydrodynamics and Joule heating aspects. Maxwell model for nanofluid has been used to incorporate the thermophysical properties of nanoelements. Formulated governing expressions have been transformed into system of ODEs by introducing similarity variables. The transformed system of ODEs is then numerically solved by Runge–Kutta–Fehlberg (RKF) method based on shooting background. The physical quantities (i.e., skin-friction coefficient, Nusselt and Bejan numbers) of scientific interest are formulated and illustrated via various plots. Graphical representations of squeezing function, temperature profile and velocity profile have been made to examine the effects of involved parameters. Streamlines and isotherms patterns have been formed and discussed. To authenticate the validity of model, skin-friction values have been compared with published literature for limited version of the model. Entropy and temperature of the system are improved with the involvement of aluminum alloys in water. Symmetrical behavior of streamlines is observed for positive approach of squeezing parameter.

  • Morphological nanolayer impact on hybrid nanofluids flow due to dispersion of polymer/CNT matrix nanocomposite material
    M Zubair Akbar Qureshi, M Faisal, Qadeer Raza, Bagh Ali, Thongchai Botmart, and Nehad Ali Shah

    American Institute of Mathematical Sciences (AIMS)
    <abstract> <p>The objective of this study is to explore the heat transfer properties and flow features of an MHD hybrid nanofluid due to the dispersion of polymer/CNT matrix nanocomposite material through orthogonal permeable disks with the impact of morphological nanolayer. Matrix nanocomposites (MNC) are high-performance materials with unique properties and design opportunities. These MNC materials are beneficial in a variety of applications, spanning from packaging to biomedical applications, due to their exceptional thermophysical properties. The present innovative study is the dispersion of polymeric/ceramic matrix nanocomposite material on magnetized hybrid nanofluids flow through the orthogonal porous coaxial disks is deliberated. Further, we also examined the numerically prominence of the permeability ($ {\\mathrm{A}}_{\\mathrm{*}} $) function consisting of the Permeable Reynold number associated with the expansion/contraction ratio. The morphological significant effects of these nanomaterials on flow and heat transfer characteristics are explored. The mathematical structure, as well as empirical relations for nanocomposite materials, are formulated as partial differential equations, which are then translated into ordinary differential expressions using appropriate variables. The Runge–Kutta and shooting methods are utilized to find the accurate numerical solution. Variations in skin friction coefficient and Nusselt number at the lower and upper walls of disks, as well as heat transfer rate measurements, are computed using important engineering physical factors. A comparison table and graph of effective nanolayer thermal conductivity (ENTC) and non-effective nanolayer thermal conductivity are presented. It is observed that the increment in nanolayer thickness (0.4−1.6), enhanced the ENTC and thermal phenomena. By the enhancement in hybrid nanoparticles volume fraction (2% to 6%), significant enhancement in Nusselt number is noticed. This novel work may be beneficial for nanotechnology and relevant nanocomponents.</p> </abstract>

  • Dynamics of MHD tangent hyperbolic nanofluid with prescribed thermal conditions, random motion and thermo-migration of nanoparticles
    Muhammad Faisal, Iftikhar Ahmad, and Tariq Javed

    Informa UK Limited
    Investigation of expanding sheet flow in the existence of tiny particles is an interesting field of research and this type of mathematical model is usually expressed in the form of partial differen...

  • Nonlinear Radiative Nanofluidic Hydrothermal Unsteady Bidirectional Transport with Thermal/Mass Convection Aspects
    Muhammad Faisal, Kanayo Kenneth Asogwa, Nazek Alessa, and Karuppusamy Loganathan

    The collective effect of thermal and mass convection along with the significance of thermal radiation, heat source/sink, and magneto-nanofluid are considered. A bi-directional stretching device is used to generate the symmetry of the flowing structure. Nonlinear behavior of thermal radiation is considered here. The magnetic field is considered non-uniform and vertically upward. Significances of pedesis motion and Ludwig–Soret are also revealed in an innovative way with heat source/sink effects. The concept of symmetry is used to transmute the transport equations from PDE type to nonlinear ODE type. We solved the transformed setup numerically by adopting Keller-box method criteria with the targeted accuracy rate. Graphical interpretations are explored with code verification. It is important to conclude that friction coefficients decline for incremental values of stretching parameter (0.1≤α≤0.9), magnetic field (0.3≤M≤0.9), and unsteady parameter (0.2≤Λ≤0.9) along with the bidirectional velocity components, and the rate of heat transmission rises with temperature ratio (1.3≤Γ≤1.7) and temperature Biot number (0.3≤BiT≤0.9) amplification. Moreso, the rate of mass transfer is enhanced with growing values of pedesis motion (0.2≤Nb≤0.6), unsteady parameter and concentration Biot number (0.3≤BiC≤0.9) with opposite effect when the Ludwig–Soret parameter (0.3≤Nt≤0.6) is boosted.

  • Activation Energy Performance through Magnetized Hybrid Fe<inf>3</inf>O<inf>4</inf>–PP Nanofluids Flow with Impact of the Cluster Interfacial Nanolayer
    M. Zubair Akbar Qureshi, Qadeer Raza, Aroosa Ramzan, M. Faisal, Bagh Ali, Nehad Ali Shah, and Wajaree Weera

    The current work investigated the mass and heat transfer of the MHD hybrid nanofluid flow subject to the impact of activation energy and cluster interfacial nanolayer. The heat transport processes related to the interfacial nanolayer between nanoparticles and base fluids enhanced the base fluid’s thermal conductivity. The tiny particles of Fe3O4 and PPy were considered due to the extraordinary thermal conductivity which is of remarkable significance in nanotechnology, electronic devices, and modern shaped heat exchangers. Using the similarity approach, the governing higher-order nonlinear coupled partial differential equation was reduced to a system of ordinary differential equations (ODEs). Fe3O4–PPy hybrid nanoparticles have a considerable influence on thermal performance, and when compared to non-interfacial nanolayer thermal conductivity, the interfacial nanolayer thermal conductivity model produced substantial findings. The increase in nanolayer thickness from level 1 to level 5 had a significant influence on thermal performance improvement. Further, the heat and mass transfer rate was enhanced with higher input values of interfacial nanolayer thickness.

  • Insight into the motion of water conveying three kinds of nanoparticles shapes on a horizontal surface: Significance of thermo-migration and Brownian motion
    S. Saleem, I.L. Animasaun, Se-Jin Yook, Qasem M. Al-Mdallal, Nehad Ali Shah, and Muhammad Faisal

    Elsevier BV

  • Dynamics of Nanoplatelets in Mixed Convective Radiative Flow of Hybridized Nanofluid Mobilized by Variable Thermal Conditions
    Iftikhar Ahmad, Muhammad Faisal, Qazi Zan-Ul-Abadin, K. Loganathan, Tariq Javed, and Balachandra Pattanaik

    Hindawi Limited
    The present mathematical model discloses the effects of Boussinesq and Rosseland approximations on unsteady 3D dynamics of water-driven hybridized nanomaterial with the movements of nanoplatelets (molybdenum disulfide, MoS 2 and graphene oxide, GO ). Variable thermal conditions, namely, VST (variable surface temperature) and VHF (variable heat flux), are opted to provide temperature to the surface. MHD effects have also been used additionally to make the study more versatile. In order to transmute the transportation equations into nondimensionlized forms, similarity transformations have been adopted. The Keller-Box technique has been applied to obtain a numerical simulation of the modeled problem. The convergence of the solution for both VST and VHF cases is presented via the grid independence tactic. Thermal setup against escalating choices of power indices and nonlinear thermal radiation parameter is discussed via graphical illustrations. The rate of heat transaction has been discussed with the growing choices of mixed convection, thermal radiation, and unsteady parameters through various tabular arrangements. It is observed through the present analysis that mixed convection parameter, radiation parameter, temperature maintaining indices r , s , and unsteady parameter magnify the rate of heat transference under the control of platelet-shaped nanoparticles.

  • Convective Heat Transport in Bidirectional Water Driven Hybrid Nanofluid Using Blade Shaped Cadmium Telluride and Graphite Nanoparticles under Electromagnetohydrodynamics Process
    Iftikhar Ahmad, Qazi Zan-Ul-Abadin, Muhammad Faisal, K. Loganathan, Tariq Javed, and Ngawang Namgyel

    Hindawi Limited
    The recent progress in nanotechnology provides the concept of hybrid-class nano-fluids having advanced thermal features comparing to regular nanofluids. The idea of a hybrid nanofluid has motivated many researchers because of its convincible performance in thermal systems. The novel theme of the present effort is to scrutinize the consequences of convective heat transfer in bidirectional water driven hybridclass nanofluid involving blade shaped cadmium telluride CdTe and graphite C nanoparticles with electromagnetohydrodynamics (EMHD) process. The transport equations representing the aforementioned topic are firstly nondimensionalized by using scaling-group transformation and then tackled by the Keller-box method, numerically. The significant results for pertinent parameters have been simulated and presented graphically as well as in tabular forms. Coefficients of drag forces are diminished with the more loadings of cadmium telluride and graphite nanoparticles and opposite results are noticed in the case of the Nusselt number. Heat transport has been improved significantly with more loadings of nanoparticles from 1 wt% to 10 wt%. A comparison benchmark for a limited version of the investigation is made with the previously published data.

  • Prescribed Thermal Activity in the Radiative Bidirectional Flow of Magnetized Hybrid Nanofluid: Keller-Box Approach
    I. Ahmad, Qazi Zan-Ul-Abadin, M. Faisal, K. Loganathan, T. Javed and D. Chaudhary

    In this exploration, we decided to investigate the significance of prescribed thermal conditions on unsteady 3D dynamics of water-based radiative hybrid nanofluid with the impact of cylindrical-shaped nanosized particles (alumina ( Al 2 O 3 ) and titania ( Ti O 2 )). For physical relevancy, the impact of the Lorentz force is also included. The combination of suitable variables has been used to transform the transport equations into the system of ordinary differential equations and then numerically solved via the Keller-Box approach. Graphical illustrations have been used to predict the impact of the involved parameters on the thermal setup. Convergence analysis is presented via the grid independence approach. Skin frictions and local Nusselt numbers against various choices of involved parameters are plotted and arranged in tabular forms. It is observed through the present investigation that temperature distribution is increased with the higher choices of radiation parameter 0.0 ≤ R d ≤ 2.0 and decreased with the improvement in the choices of temperature maintaining indices (i.e., − 2.0 ≤ r , s ≤ 2.0 ). Moreover, the thermophysical properties except specific heat for hybrid nanofluid are improved with the involvement of cylindrical-shaped nanoparticles. The temperature of the hybrid nanofluid is observed to be higher for variable thermal conditions as compared to uniform thermal conditions. Outfalls for a limited version of the report have been compared with a previous published paper.

  • Entropy Analysis in Bidirectional Hybrid Nanofluid Containing Nanospheres with Variable Thermal Activity
    I. Ahmad, Qazi Zan-Ul-Abadin, M. Faisal, K. Loganathan, T. Javed and Sonam Gyeltshen

    An investigation of the thermal performance of water-conveying nanospheres (magnetite ( Fe 3 O 4 ) and silver ( Ag )) subject to variable thermal controls, namely, variable surface temperature and variable surface normal heat flux, has been made. A bidirectionally elongating surface is used to generate an unsteady flow mechanism with the action of the Lorentz force. Derived equations of basic laws are firstly nondimensionalized and then numerically solved by applying the Keller-Box method. The local Nusselt number for both the thermal cases is calculated and discussed. Percent-wise enhancement in the rate of heat transport has also been included in the analysis. It was concluded through the present exploration that at lower volume fractions of magnetite and silver, the rate of heat transport is observed to be dominant. The rate of heat transference has attained identical values for both the provided thermal conditions at the surface. Moreover, intensities of velocity and thermal profiles diminish with the appreciation of the choice of unsteadiness. The temperature-controlling indices also affect the thermal profile, and it is reduced with the intensification in the considerations of these indices. The values of thermal conductivity, density, and electrical conductivity have been improved with the inclusion of nanospheres (magnetite ( Fe 3 O 4 ) and silver ( Ag )), whereas the value of specific heat is reduced with the mixture of these nanospheres. The Nusselt number is increased up to 5% with the involvement of magnetite nanospheres, and it is enhanced up to 4% with the involvement of silver nanospheres.

  • Nonlinear Mixed Convective Bidirectional Dynamics of Double Stratified Radiative Oldroyd-B Nanofluid Flow with Heat Source/Sink and Higher-Order Chemical Reaction
    Iftikhar Ahmad, Muhammad Faisal, K. Loganathan, Muhammad Zaheer Kiyani, and Ngawang Namgyel

    Hindawi Limited
    Little is known in the literature about the concept of nonuniform heat source/sink and higher-order chemical reaction for the dynamics of Oldroyd-B nanoparticles. Therefore, the present article addresses the nonuniform heat source/sink and higher-order chemical reaction features in nonlinear mixed convection bidirectional MHD dynamics of Oldroyd-B nanoparticles with thermal radiation aspects through porous space. Stratification effects for both the temperature and concentration setups are also used in the mathematical model with the significance of random movement and thermodiffusion of nanoparticles. Shape-preserving transformations have been employed to convert the transport equations into solvable forms. An innovative analytical tactic, namely, homotopy analysis method, has been adopted to find the solution of the modeled problem. Behaviors of pertinent parameters on thermal and concentration profiles have been discussed through various graphs. Inspection of heat/mass transport against appropriate varieties of pertinent parameters has been made and explained physically. Thermal profile is augmented with the higher estimations of space and temperature-dependent heat source/sink links. Concentration profile is diminished with the augmentation of higher-order chemical reaction parameter. Sherwood number is improved with the estimation of 0 ≤ β t ≤ 100 and is reduced with the growth of 0 ≤ β c ≤ 100 . Nusselt number is declined with the upgraded amounts of 0 ≤ N b ≤ 3 and 0 ≤ N t ≤ 5 .

  • Entropy and Heat Transfer Analysis for MHD Flow of Cu/Ag-Water-Based Nanofluid on a Heated 3D Plate with Nonlinear Radiation
    S. Eswaramoorthi, S. Divya, Muhammad Faisal, and Ngawang Namgyel

    Hindawi Limited
    This paper scrutinizes the consequences of radiation and heat consumption of MHD convective flow of nanofluid on a heated stretchy plate with injection/suction and convective heating/cooling conditions. The nanofluid encompasses with C u and A g nanoparticles. We enforce the suited transformation to remodel the governing mathematical models to ODE models. The HAM (homotopy analysis method) idea is applied to derive the series solutions. The divergence of fluid velocity, temperature, skin friction coefficient, local Nusselt number, entropy generation, and Bejan number on disparate governing parameters is exhibited via graphs and tables. It is seen that the fluid velocity in both directions is subsided when elevating the magnetic field and Forchheimer number. Also, the C u nanoparticles possess hefty speed compared to A g nanoparticles because the density of A g nanoparticles is high compared to that of C u nanoparticles. The fluid temperature upturns when enlarging the heat generation and radiation parameters. The skin friction coefficients and local Nusselt number are high in A g nanoparticles than in C u nanoparticles.

  • Unsteady 3D heat transport in hybrid nanofluid containing brick shaped ceria and zinc-oxide nanocomposites with heat source/sink
    Iftikhar Ahmad, Muhammad Faisal, Qazi Zan-Ul-Abadin, Tariq Javed, and K. Loganathan

    Informa UK Limited

  • Numerical investigation for mixed convective 3D radiative flow of chemically reactive Williamson nanofluid with power law heat/mass fluxes
    Iftikhar Ahmad, Muhammad Faisal, Tariq Javed, Ayesha Mustafa, and Muhammad Zaheer Kiyani

    Elsevier BV

  • Insight into the relationship between unsteady Cattaneo-Christov double diffusion, random motion and thermo-migration of tiny particles
    Iftikhar Ahmad, Muhammad Faisal, Tariq Javed, and Isaac Lare Animasaun

    Elsevier BV
    Abstract Significant impact of unsteady Cattaneo-Christov double diffusion, random motion and thermo-migration of tiny particles on thermal storage and in heat exchangers devices has attracted us to investigate the role of these thermal engineering terms in the existence of convectively heating and zero-mass flux at the bidirectionally stretchable wall. Combination of similarity expressions is adopted to transform the prevailing partial differential equations into ordinary differential equations and then solved, numerically, by using Keller-Box simulation technique. Temperature and concentration phases are discussed under the influences of involved parameters, graphically. It is detected that thermal relaxation dependent on heat flux maximizes the temperature gradient, whereas thermal relaxation dependent on mass flux minimizes the concentration gradient. It is also scrutinized that escalating choice of random motion of tiny particles provides constant temperature gradient because of Nield’s expression. Moreover, augmenting choice of thermo-migration of tiny particles enhances the concentration gradient. Furthermore, intensifying choice of unsteady factor abridges the temperature plus concentration distributions and also it condenses all the three boundary layer thicknesses.

  • Unsteady rotating flow of nanofluid with prescribed thermal aspects
    Iftikhar Ahmad, Muhammad Faisal, and Tariq Javed

    World Scientific Pub Co Pte Lt
    Motivated by the significant role of nanofluid in pollution cleaning and energy recovery, we decided to explore the unsteady three-dimensional rotating flow of nanofluid driven by the movement of a flat surface with the potencies of prescribed heat distributions. The modeling of the physical model is completed with the help of Buongiorno nanofluid model. Suitable arrangement of similarity variables is implemented to transform the model equations into strongly nonlinear ordinary differential equations. Numerical inspection of the model is made by employing Keller–Box algorithm. Influences of involved parameters on the distributions of heat and mass are discussed graphically, while the potencies of influential parameters on reduced Nusselt and reduced Sherwood numbers are physically discussed through tabular arrangements. It is deduced that increasing the values of Prandtl factor and heat controlling indices diminishes the temperature and concentration distributions, whereas intensification in the amount of rotation factor enhances the temperature as well as concentration distribution. Moreover, negative trends in the amounts of reduced Nusselt and Sherwood numbers are achieved with the escalations in the values of rotation and thermophoresis factors, whereas opposite trend is achieved with the intensification in the choice of Prandtl factor.