scholarly journals Stability Analysis of Darcy-Forchheimer Flow of Casson Type Nanofluid Over an Exponential Sheet: Investigation of Critical Points

Symmetry ◽  
2019 ◽  
Vol 11 (3) ◽  
pp. 412 ◽  
Author(s):  
Liaquat Ali Lund ◽  
Zurni Omar ◽  
Ilyas Khan ◽  
Jawad Raza ◽  
Mohsen Bakouri ◽  
...  

In this paper, steady two-dimensional laminar incompressible magnetohydrodynamic flow over an exponentially shrinking sheet with the effects of slip conditions and viscous dissipation is examined. An extended Darcy Forchheimer model was considered to observe the porous medium embedded in a non-Newtonian-Casson-type nanofluid. The governing equations were converted into nonlinear ordinary differential equations using an exponential similarity transformation. The resultant equations for the boundary values problem (BVPs) were reduced to initial values problems (IVPs) and then shooting and Fourth Order Runge-Kutta method (RK-4th method) were applied to obtain numerical solutions. The results reveal that multiple solutions occur only for the high suction case. The results of the stability analysis showed that the first (second) solution is physically reliable (unreliable) and stable (unstable).

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Golam Mortuja Sarkar ◽  
Suman Sarkar ◽  
Bikash Sahoo

Purpose This paper aims to theoretically and numerically investigate the steady two-dimensional (2D) Hiemenz flow with heat transfer of Reiner-Rivlin fluid over a linearly stretching/shrinking sheet. Design/methodology/approach The Navier–Stokes equations are transformed into self-similar equations using appropriate similarity transformations and then solved numerically by using shooting technique. A simple but effective mathematical analysis has been used to prove the existence of a solution for stretching case (λ> 0). Moreover, an attempt has been laid to carry the asymptotic solution behavior for large stretching. The obtained asymptotic solutions are compared with direct numerical solutions, and the comparison is quite remarkable. Findings It is observed that the self-similar equations exhibit dual solutions within the range [λc, −1] of shrinking parameter λ, where λc is the turning point from where the dual solutions bifurcate. Unique solution is found for all stretching case (λ > 0). It is noticed that the effects of cross-viscous parameter L and shrinking parameter λ on velocity and thermal fields show opposite character in the dual solution branches. Thus, a linear temporal stability analysis is performed to determine the basic feasible solution. The stability analysis is based on the sign of the smallest eigenvalue, where positive or negative sign leading to a stable or unstable solution. The stability analysis reveals that the first solution is stable that describes the main flow. Increase in cross-viscous parameter L resulting in a significant increment in skin friction coefficient, local Nusselt number and dual solutions domain. Originality/value This work’s originality is to examine the combined effects of cross-viscous parameter and stretching/shrinking parameter on skin friction coefficient, local Nusselt number, velocity and temperature profiles of Hiemenz flow over a stretching/shrinking sheet. Although many studies on viscous fluid and nanofluid have been investigated in this field, there are still limited discoveries on non-Newtonian fluids. The obtained results can be used as a benchmark for future studies of higher-grade non-Newtonian flows with several physical aspects. All the generated results are claimed to be novel and have not been published elsewhere.


2011 ◽  
Vol 110-116 ◽  
pp. 3184-3190
Author(s):  
Necdet Bildik ◽  
Duygu Dönmez Demir

This paper deals with the solutions of lateral heat loss equation by using collocation method with cubic B-splines finite elements. The stability analysis of this method is investigated by considering Fourier stability method. The comparison of the numerical solutions obtained by using this method with the analytic solutions is given by the tables and the figure.


2006 ◽  
Vol 18 (3) ◽  
pp. 242-248 ◽  
Author(s):  
Mizuho Shibata ◽  
◽  
Shinichi Hirai

To analyze the stability of dynamic control through asoft interface-the viscoelastic material between a manipulating finger and a manipulated object- we model dynamic control through the soft interface in continuous-discrete time. We then formulate dynamics using a modified z-transform in continuous-discrete time for feedback and feedforward control. We show that system stability depends on the viscoelasticity of the soft interface for feedback control. The relationship between material viscosity and sampling time in critical stability is not monotonous, a phenomenon we analyze by root locus. We compare stability analysis by the modified z-transform, simulations based on the Runge-Kutta method, and a regular z-transform, demonstrating that the relationship is specific to a continuous-discrete time.


2014 ◽  
Vol 31 (12) ◽  
pp. 2795-2808 ◽  
Author(s):  
Tim Rees ◽  
Adam Monahan

Abstract The stability analysis of stratified parallel shear flows is fundamental to investigations of the onset of turbulence in atmospheric and oceanic datasets. The stability analysis is performed by considering the behavior of small-amplitude waves, which is governed by the Taylor–Goldstein (TG) equation. The TG equation is a singular second-order eigenvalue problem, whose solutions, for all but the simplest background stratification and shear profiles, must be computed numerically. Accurate numerical solutions require that particular care be taken in the vicinity of critical layers resulting from the singular nature of the equation. Here a numerical method is presented for finding unstable modes of the TG equation, which calculates eigenvalues by combining numerical solutions with analytical approximations across critical layers. The accuracy of this method is assessed by comparison to the small number of stratification and shear profiles for which analytical solutions exist. New stability results from perturbations to some of these profiles are also obtained.


Symmetry ◽  
2019 ◽  
Vol 11 (9) ◽  
pp. 1078 ◽  
Author(s):  
Mohamad Mustaqim Junoh ◽  
Fadzilah Md Ali ◽  
Ioan Pop

The revised Buongiorno’s nanofluid model with the effect of induced magnetic field on steady magnetohydrodynamics (MHD) stagnation-point flow of nanofluid over a stretching or shrinking sheet is investigated. The effects of zero mass flux and suction are taken into account. A similarity transformation with symmetry variables are introduced in order to alter from the governing nonlinear partial differential equations into a nonlinear ordinary differential equations. These governing equations are numerically solved using the bvp4c function in Matlab solver, a very adequate finite difference method. The influences of considered parameters ( P r , M, χ , L e , N b , N t , S, and λ ) on velocity, induced magnetic, temperature, and concentration profiles together with the reduced skin friction and heat transfer rate are discussed. Results from these criterion exposed the existence of dual solutions when magnetic field and suction are applied for a specific range of λ . The stability of the solutions obtained is carried out by performing a stability analysis.


Symmetry ◽  
2020 ◽  
Vol 12 (10) ◽  
pp. 1637
Author(s):  
Liaquat Ali Lund ◽  
Zurni Omar ◽  
Sumera Dero ◽  
Dumitru Baleanu ◽  
Ilyas Khan

This article aims to study numerically the rotating, steady, and three-dimensional (3D) flow of a hybrid nanofluid over an exponentially shrinking sheet with the suction effect. We considered water as base fluid and alumina (Al2O3), and copper (Cu) as solid nanoparticles. The system of governing partial differential equations (PDEs) was transformed by an exponential similarity variable into the equivalent system of ordinary differential equations (ODEs). By applying a three-stage Labatto III-A method that is available in bvp4c solver in the Matlab software, the resultant system of ODEs was solved numerically. In the case of the hybrid nanofluid, the heat transfer rate improves relative to the viscous fluid and regular nanofluid. Two branches were obtained in certain ranges of the involved parameters. The results of the stability analysis revealed that the upper branch is stable. Moreover, the results also indicated that the equations of the hybrid nanofluid have a symmetrical solution for different values of the rotation parameter (Ω).


2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
Puntani Pongsumpun ◽  
I-Ming Tang

The respiratory disease caused by the Influenza A Virus is occurring worldwide. The transmission for new strain of the H1N1 Influenza A virus is studied by formulating a SEIQR (susceptible, exposed, infected, quarantine, and recovered) model to describe its spread. In the present model, we have assumed that a fraction of the infected population will die from the disease. This changes the mathematical equations governing the transmission. The effect of repetitive contact is also included in the model. Analysis of the model by using standard dynamical modeling method is given. Conditions for the stability of equilibrium state are given. Numerical solutions are presented for different values of parameters. It is found that increasing the amount of repetitive contacts leads to a decrease in the peak numbers of exposed and infectious humans. A stability analysis shows that the solutions are robust.


2008 ◽  
Vol 05 (03) ◽  
pp. 613-642 ◽  
Author(s):  
ALY FATHY ◽  
CHENG WANG ◽  
JOSHUA WILSON ◽  
SONGNAN YANG

The Maxwell equations are solved by a long-stencil fourth order finite difference method over a Yee grid, in which different physical variables are located at staggered mesh points. A careful treatment of the numerical values near the boundary is introduced, which in turn leads to a "symmetric image" formula at the "ghost" grid points. Such a symmetric formula assures the stability of the boundary extrapolation. In turn, the fourth order discrete curl operator for the electric and magnetic vectors gives a complete set of eigenvalues in the purely imaginary axis. To advance the dynamic equations, the four-stage Runge–Kutta method is utilized, which results in a full fourth order accuracy in both time and space. A stability constraint for the time step is formulated at both the theoretical and numerical levels, using an argument of stability domain. An accuracy check is presented to verify the fourth order precision, using a comparison between exact solution and numerical solutions at a fixed final time. In addition, some numerical simulations of a loss-less rectangular cavity are also carried out and the frequency is measured precisely.


Author(s):  
Nur Auni Baharum ◽  
Zanariah Abdul Majid ◽  
Norazak Senu

The performance of the numerical computation based on the diagonally implicit multistep block method for solving Volterra integrodifferential equations (VIDE) of the second kind has been analyzed. The numerical solutions of VIDE will be computed at two points concurrently using the proposed numerical method and executed in the predictor-corrector (PECE) mode. The strategy to obtain the numerical solution of an integral part is discussed and the stability analysis of the diagonally implicit multistep block method was investigated. Numerical results showed the competence of diagonally implicit multistep block method when solving Volterra integrodifferential equations compared to the existing methods.


Sign in / Sign up

Export Citation Format

Share Document