Journal of Applied and Computational MechanicsJournal of Applied and Computational Mechanics
http://iranjournals.nlai.ir/1348/
Tue, 25 Feb 2020 03:18:40 +0100FeedCreatorJournal of Applied and Computational Mechanics
http://iranjournals.nlai.ir/1348/
Feed provided by Journal of Applied and Computational Mechanics. Click to visit.An Efficient Implementation of Phase Field Method with Explicit Time Integration
http://iranjournals.nlai.ir/1348/article_661965_49855.html
The phase field method integrates the Griffith theory and damage mechanics approach to predict crack initiation, propagation, and branching within one framework. No crack tracking topology is needed, and complex crack shapes can be captures without user intervention. In this paper, a detailed description of how the phase field method is implemented with explicit dynamics into LS-DYNA is provided. The displacement field and the damage field are solved in a staggered approach and the phase field equation is solved every Nth time step (N is refered to as calculation cycle) to save computational time. An N value smaller than 1/400 of the total time step numbers is suggested. Several simulations are presented to demonstrate the feasibility of this solving scheme.Tue, 30 Jun 2020 19:30:00 +0100Finite Element Analysis of Low Velocity Impact on Carbon Fibers/Carbon Nanotubes Reinforced ...
http://iranjournals.nlai.ir/1348/article_661966_49855.html
An effort is made to gain insight on the effect of carbon nanotubes (CNTs) on the impact response of carbon fiber reinforced composites (CFRs) under low velocity impact. Certain amount of CNTs could lead improvements in mechanical properties of composites. In the present investigation, ABAQUS/Explicit finite element code (FEM) is employed to investigate various damages modes of nano composites including matrix cracking, fiber damage and delamination by employing Hashin’s criterion and cohesive zone modeling. The obtained results for 0, 0.5, 1, 2 and 4% CNTs demonstrate that by including CNTs in composite plates, damage could be reduced. However, adding further CNTs causes sudden reduction of impact tolerance capability of the composite plates, particularly, damage due to delamination.Tue, 30 Jun 2020 19:30:00 +0100Simple Two Variable Refined Theory for Shear Deformable Isotropic Rectangular Beams
http://iranjournals.nlai.ir/1348/article_661967_49855.html
In this paper, a displacement-based, variationally consistent, two variable refined theory for shear deformable beams is presented. The beam is assumed to be of linearly elastic, homogeneous, isotropic material and has a uniform rectangular cross-section. In this theory, the beam axial displacement and beam transverse displacement consist of bending components and shearing components. The assumed displacement field of this theory is such that, bending components do not take part in the cross-sectional shearing force, and shearing components do not take part in the cross-sectional bending moment. This theory utilizes linear strain-displacement relations. The displacement functions give rise to the beam transverse shear strain (and hence to the beam transverse shear stress) which varies quadratically through the beam thickness and maintains transverse shear stress-free beam surface conditions. Hence the shear correction factor is not required. Hamilton’s principle is utilized to derive governing differential equations and variationally consistent boundary conditions. This theory involves only two governing differential equations of fourth-order. These governing equations are only inertially coupled for the case of dynamics and are decoupled for the case of statics. This theory is simple and has a strong resemblance with the Bernoulli-Euler beam theory. To demonstrate the efficacy of the present theory, illustrative examples pertain to the static bending and free vibrations of shear deformable isotropic rectangular beams are presented.Tue, 30 Jun 2020 19:30:00 +0100Vibration Analysis of Different Types of Porous FG Conical Sandwich Shells in Various Thermal ...
http://iranjournals.nlai.ir/1348/article_661968_49855.html
Vibration behavior of different types of porous functionally graded (FG) conical sandwich shells are investigated based on a modified high order sandwich shells theory for the first time. Sandwich shell includes FG face sheets covering a homogeneous core and the second one includes homogeneous face sheets and a FG core. Power law rule modified by considering two types of porosity distributions is used to model the functionally graded materials. All materials are temperature dependent and uniform, linear and nonlinear temperature distributions are used to model the effect of the temperature variation in the sandwiches. Governing equations are obtained by the Hamilton's energy principle and solved with Galerkin method. To verify the results, they are compared with ones achieved by finite element method obtained by Abaqus software for special cases with the results in literatures.Tue, 30 Jun 2020 19:30:00 +0100Heat Transfer Analysis of Nanofluid Flow with Porous Medium through Jeffery Hamel ...
http://iranjournals.nlai.ir/1348/article_661969_49855.html
In this paper, flow and heat transfer of nanofluid through a converging or diverging channel with porous medium is investigated. The fluid constantly flows under the effect of magnetic field through the channel. The diverging/converging fluid motion is modeled using the momentum and energy equations. The influence of some parameters such as opening channel angle, Reynolds number and Darcy’s number when the nanofluid flows through the non-parallel plates are studied. It is seen that high Reynolds number enhances the fluid viscosity while decreases velocity. Similarly, heat transfer reduces at high Darcy’s number owing to decreased flow consequently internal friction reduces. The obtained results in comparison with the similar studies in the literatures show satisfactory agreement.Tue, 30 Jun 2020 19:30:00 +0100On Green and Naghdi Thermoelasticity Model without Energy Dissipation with Higher Order Time ...
http://iranjournals.nlai.ir/1348/article_661970_49855.html
In the present work, a modified model of heat conduction including higher order of time derivative is derived by extending Green and Naghdi theory without energy dissipation. We introduce two phase lag times to include the thermal displacement gradient and the heat flux in the heat conduction and depict microscopic responses more precisely. The constructed model is applied to study thermoelastic waves in a homogeneous and isotropic perfect conducting unbounded solid body containing a spherical cavity. We use the Laplace transform method to analyze the problem. The solutions for the field functions are obtained numerically using the numerical Laplace inversion technique. The results are analyzed in different tables and graphs and compared with those obtained earlier in the contexts of some other theories of thermoelasticity.Tue, 30 Jun 2020 19:30:00 +0100Thermal Buckling Analysis of Functionally Graded Euler-Bernoulli Beams with ...
http://iranjournals.nlai.ir/1348/article_661971_49855.html
Thermal buckling behavior of functionally graded Euler-Bernoulli beams in thermal conditions is investigated analytically. The beam with material and thermal properties dependent on the temperature and position is considered. Based on the transformed-section method, the functionally graded beam is considered as an equivalent homogeneous Euler-Bernoulli beam with an effective bending rigidity under an eccentric thermal load. Then, the thermal elastic buckling equation associated with the bending deflection about the neutral axis is established. The easily usable closed-form solutions for the critical thermal buckling temperature of functionally graded beams under uniform and non-linear temperature rise are obtained and used to calculate the thermal buckling temperature. Some results are evaluated and compared with those by other investigators to validate the accuracy of the presented method. The effects of material compositions, temperature-dependent material properties, slenderness ratios and restraint conditions on thermal buckling behaviors are discussed. It is believed that the proposed model provides engineers and designers an easy and useful method to investigate the effects of various parameters affecting the thermal buckling characteristics of functionally graded beams.Tue, 30 Jun 2020 19:30:00 +0100Stability Assessment of the Flexible System using Redundancy
http://iranjournals.nlai.ir/1348/article_661972_49855.html
In this study, dynamic behavior of a mooring line in a floating system is analyzed by probability approaches. In dynamics, most researches have shown the system model and environments by mathematical expression. We called this process as the forward dynamics. However, there is a limit to define the exact environments because of uncertainty. To consider uncertainty, we introduce the redundancy in flexible system, mooring line. For verifying the effectiveness and stability of the mooring line, criterion of axial breaking load of the mooring line is applied to joint reaction forces according to the various path of the mooring line. To cover the limits for defining the non-linearity of the environments, various responses of the mooring line along the redundancy that is used in Robotics, are derived by probability distribution. By using the Newton-Euler formulation, the inverse kinematics and the linear acceleration theorem to get joint displacements, velocities and accelerations, the joint reaction forces and moments are calculated and probability distribution of the mooring about stability and compatibility is investigated. Lastly, we simulate the flexible systems in various null motions, calculated each joint torque and force, and evaluated failure probabilities using the Monte-Carlo method.Tue, 30 Jun 2020 19:30:00 +0100Emotional Learning Based Intelligent Controller for MIMO Peripheral Milling Process
http://iranjournals.nlai.ir/1348/article_661973_49855.html
During the milling process, one of the most important factors in reducing tool life expectancy and quality of workpiece is the chattering phenomenon due to self-excitation. The milling process is considered as a MIMO strongly coupled nonlinear plant with time delay terms in cutting forces. We stabilize the plant using two independent Emotional Learning-based Intelligent Controller (ELIC) in parallel. Control inputs are considered as forces Ux and Uy in two directions x and y, which are applied by the piezoelectrics. The ELIC consists of three elements; Critic, TSK controller and the learning element. The results of the ELIC have been compared with a Sliding Mode Controller (SMC). The simulation for the nominal plant shows better performance of the ELIC in IAE and ITSE values at least 86% in the x-direction and 79% in the y-direction. Similar simulation for an uncertain plant also shows an improvement of at least 89% in the x-direction and 97% in the y-direction.Tue, 30 Jun 2020 19:30:00 +0100A FEM Multiscale Homogenization Procedure using Nanoindentation for High Performance Concrete
http://iranjournals.nlai.ir/1348/article_661974_49855.html
This paper aims to develop a numerical multiscale homogenization method for prediction of elasto-viscoplastic properties of a high performance concrete (HPC). The homogenization procedure is separated into two-levels according to the microstructure of the HPC: the mortar or matrix level and the concrete level. The elasto-viscoplastic behavior of individual microstructural phases of the matrix are identified from nanoindentation data using an inverse identification method. The micromechanical results are then used as input parameters for numerical elasto-viscoplastic homogenization at microscale. The mortar level is analyzed with numerical homogenization by using the finite element simulation to predict the overall elasto-viscoplastic properties of HPC. The results are compared with macroscopic experimental and analytical results from the literature showing a good agreement.Tue, 30 Jun 2020 19:30:00 +0100Using the Finite Element Analysis Method to Study the 3-point Bending Test for the ...
http://iranjournals.nlai.ir/1348/article_661975_49855.html
An elastic finite element analysis was conducted to evaluate the stress distribution in the initiation zone of the adhesive rupture during the 3-point bending test. This test is used to measure the adherence between a polyepoxy adhesive and aluminum alloy with different surface treatments. The purpose is to compare, in the high stress concentration areas, the stress fields calculated using finite element method with the experimental data obtained in different configurations. Focusing on the load level at crack initiation, on the localization and the size of adhesive failure initiation, a local criterion for adhesive fracture is proposed based on the value of the stress normal to the interface.Tue, 30 Jun 2020 19:30:00 +0100Modelling of Love Waves in Fluid Saturated Porous Viscoelastic Medium resting over an ...
http://iranjournals.nlai.ir/1348/article_661976_49855.html
The present article is devoted to a theoretical study on Love wave vibration in a pre-stressed fluid-saturated anisotropic porous viscoelastic medium embedded over an inhomogeneous isotropic half-space influenced by gravity. The expression of dispersion has been achieved with the help of mathematical tools such as variable separable method and Whittaker’s function’s expansion under certain boundary conditions. After that, the obtained result has been coincided with the pre-established classical equation of Love wave, as shown in the section of particular case and validation. The substantial influence of various affecting factors like gravity, initial stress, porosity, viscosity and inhomogeneity on dispersion curves of Love wave has been investigated extensively by means of graphical depictions and discussions accomplished by numerical results.Tue, 30 Jun 2020 19:30:00 +0100Numerical Scrutinization of Three Dimensional Casson-Carreau Nano Fluid Flow
http://iranjournals.nlai.ir/1348/article_661977_49855.html
This study presents the computational analysis of three dimensional Casson and Carreau nanofluid flow concerning the convective conditions. To do so, the flow equations are modified to nonlinear system of ODEs after using appropriate self-similarity functions. The solution for the modified system is evaluated by numerical techniques. The results show the impacts of involving variables on flow characteristics and the outcomes of the friction factors are evaluated as well. In this study, the outcomes to local Nusselt number and Sherwood numbers are evaluated. Favourable comparison is performed with previously available outcomes. The achieved results are similar to solutions obtained by other researchers. The results are presented for flow characteristics in the case of Casson and Carreau fluids. Velocities are reduced for the growing values of permeability and velocity slip parameters in case of Casson and Carreau nanofluids. Temperature field enhances with the hike in the estimations of thermophoresis parameter and the thermal Biot number in case of Casson and Carreau nanofluids. Enhancing values of velocity slip parameter results in decrease in the skin friction coefficients and the rate of heat transfer, and rise in the rate of mass transfer in case of Casson and Carreau nanofluids.Tue, 30 Jun 2020 19:30:00 +0100Sealing Performance of the End Fitting of a Marine Unbonded Flexible Pipe under Pressure Penetration
http://iranjournals.nlai.ir/1348/article_661978_49855.html
The sealing performance of end fittings is very important for offshore oil and gas pipelines. To investigate the sealing behavior of a ring-shaped wedge seal, global and local numerical models of the ring–pipe interaction have been developed based on the finite-element method. First, the sealing process of the ring under pressure is simulated. Second, a criterion for the penetration of fluid pressure is applied in these models to assess how the sealing capacity changes along the contact surface. Finally, the contact magnitude of interference and the shape of the sawtooth heaves on the sealing ring are predicted and compared. The results show an interesting concentration of von Mises stress in the sealing ring and also that the peak contact pressure appears in the sealing zone. However, the penetration of fluid pressure has obvious effects on the distributions of von Mises stress and contact pressure. The best sealing performance is when the axial displacement of the sealing ring is 1.4 mm and the contact magnitude of interference is 0.3 mm. Given the sawtooth heave of sealing ring, semicircular heave gives the better sealing capacity compared with trapezoidal heave.Tue, 30 Jun 2020 19:30:00 +0100Effect of Tool Shoulder and Pin Cone Angles in Friction Stir Welding using Non-circular Tool Pin
http://iranjournals.nlai.ir/1348/article_661979_49855.html
In friction stir welding frictional heat is generated by the rotating tool, sliding over the stationary plate along the weld centre. Tool being the only source of heat producing member, its geometrical design influences the heat generation rate. In this present work, effects of variation in tool shoulder and tool pin taper angles on thermal history during joining are analysed. Tools with triangular and hexagonal tool pins are used to understand the influence of tool pin shape on process temperature. An analytical heat input model is developed for tools with non-circular tool pins and a comparative study is carried out between the hexagonal and triangular tool pins on temperature distribution using a three dimensional Matlab model. Proposed model is validated through experimental analysis. Apart from this, regression model based comparative study is carried out on the variation in temperature response to the change in tool pin shape, tool shoulder and tool pin taper angle.Tue, 30 Jun 2020 19:30:00 +0100Impact of Blood Vessel Wall Flexibility on the Temperature and Concentration Dispersion
http://iranjournals.nlai.ir/1348/article_661980_49855.html
The analysis of solute and thermal dispersion in pulsatile flow through the stenotic tapered blood vessel is presented. The present problem is an extension of the work done by Ramana et al. who considered the time-invariant arterial wall. In the present model, the flexible nature of the arterial wall through the obstruction (called stenosis) is considered and it is achieved with the help of period trigonometric function. In the present study, the impact of the time-dependent arterial wall on the blood flow dynamics is discussed in details. The rheology of the blood is modeled as a couple stress fluid. The proposed fluid model is the isothermal inclusion of temperature-sensitive drug coated Titanium dioxide Nano-particles in the couple stress fluid for examining the concentration and temperature dispersion. The effects of the catheter and permeability of the stenosis are considered in the model. Care has been taken to model the thermo-physical properties of the fluid with the immersed nanoparticle, e.g., TiO2, Ag and Cu. The modeled non-linear and coupled equations are solved by using the Homotopy Perturbation Method. The temperature and concentration dispersion effects are in the flexible stenotic arterial vessel under the pulsatile physiological pressure gradient are studied and reported in details. The alterations in the axial velocity, resistance to the flow, and wall shear stress are studied and found out that the high intense vortex regions are identified in the stenotic region. The model has direct applications in the pharmaceutical industry in design and developing the drug to treat stenotic conditions.Tue, 30 Jun 2020 19:30:00 +0100Prediction of Entrance Length for Magnetohydrodynamics Channels Flow using Numerical simulation ...
http://iranjournals.nlai.ir/1348/article_661981_49855.html
This paper focuses on using the numerical finite volume method (FVM) and artificial neural network (ANN) in order to propose a correlation for computing the entrance length of laminar magnetohydrodynamics (MHD) channels flow. In the first step, for different values of the Reynolds (Re) and Hartmann (Ha) numbers (600ReL increases.Tue, 30 Jun 2020 19:30:00 +0100A Hybrid Particle Swarm Optimization and Genetic Algorithm for Truss Structures with Discrete ...
http://iranjournals.nlai.ir/1348/article_661982_49855.html
A new hybrid algorithm of Particle Swarm Optimization and Genetic Algorithm (PSOGA) is presented to get the optimum design of truss structures with discrete design variables. The objective function chosen in this paper is the total weight of the truss structure, which depends on upper and lower bounds in the form of stress and displacement limits. The Particle Swarm Optimization basically modeled the social behavior of birds on the basis of the fact that Individual birds exchange information about their position, velocity, fitness, and on the basis that the behavior of the flock is then influenced to increase the probability of migration to other regions with high fitness. One of the problems of PSO is that it is easily trapped at the local point due to its non-uniform movement. The present study uses the mutation, random selection, and reproduction to reach the best genetic algorithm with the operators of natural genetics. Therefore, only identical chromosomes or particles can be converged. In other words, PSO and GA algorithm goes from one point in the search space to another point, interacting with each other. In this way, this helps them to find the optimum design by means of deterministic and probabilistic rules. The present study merged the two algorithms together in order to design several benchmark truss structures, and then the results of the new algorithm compared to those of other evolutionary optimization methods.Tue, 30 Jun 2020 19:30:00 +0100Multivariate Jeffrey Fluid Flow past a Vertical Plate through Porous Medium
http://iranjournals.nlai.ir/1348/article_661983_49855.html
An analysis is suggested to study the impact of Hall currents in Jeffrey fluid which is chemically reactive through a porous medium limited by a semi-infinite vertical permeable plate within the existence of heat generation. An evenly distributed magnetic field turns vertically on the porous surface which absorbs the Jeffrey fluid with a changed suction velocity with time. The analytical expressions are solved by means of three terms harmonic and non-harmonic functions. Statistical calculations are carried out for the point of resultant outcomes which are shown graphically and the impacts of the parameters velocity, temperature and concentration are listed. In addition, the results of skin-friction coefficient (τ), Nusselt number (Nu) and Sherwood number (Sh) are taken in to consideration. It is revealed that the impact of the Hall parameter on the channel velocities and skin friction coefficient is subjected to the estimation of the wall suction parameter.Tue, 30 Jun 2020 19:30:00 +0100On the Six Node Hexagon Elements for Continuum Topology Optimization of Plates Carrying in ...
http://iranjournals.nlai.ir/1348/article_661984_49855.html
The need of polygonal elements to represent the domain is gaining interest among structural engineers. The objective is to perform static analysis and topology optimization of a given continuum domain using the rational fraction type shape functions of six node hexagonal elements. In this paper, the main focus is to perform the topology optimization of two-dimensional plate structures using Evolutionary Swarm Intelligence Firefly Algorithms (ESIFA) and three-dimensional shell structures using optimality criteria. The optimization of plates carrying in plane loading is performed with minimum weight as objective. Two different types of shell structures are optimized using maximum strain energy as criteria. The optimal distribution of the material in the design domain obtained using six node hexagon elements is compared with the optimal distribution of material obtained using quadrilateral elements. A few problems from the literature have been solved and this study has proved that hexagon element gives better results over traditional quadrilateral elements.Tue, 30 Jun 2020 19:30:00 +0100Effects of Non-uniform Suction, Heat Generation/Absorption and Chemical Reaction with ...
http://iranjournals.nlai.ir/1348/article_661985_49855.html
In the present investigation, the magnetohydrodynamic Falkner-Skan flow of tangent hyperbolic nanofluids over a stretching/shrinking wedge with variable suction, internal heat generation/absorption and chemical reaction with activation energy have been scrutinized. Nanofluid model is composed of “Brownian motion’’ and “Thermophoresis’’. Transformed non-dimensional coupled non-linear equations are solved by adopting the fourth-order R-K method along with the shooting technique. A comprehensive analysis of nanofluid velocity, the relative temperature, and its concentration profiles has been addressed. The major outcomes of the current study include that augmentation in the Weissenberg parameter, Hartmann number along with suction impede fluid flow and the shrinkage of the related boundary layer while internal heating develops an ascending thermal boundary layer for static and moving (stretching/shrinking) wedge. An increment in reaction rate undermines the nanoparticle concentration while that of activation energy exhibits a reverse trend.Tue, 30 Jun 2020 19:30:00 +0100Effect of Chemical Reaction on Bioconvective Flow in Oxytactic Microorganisms Suspended Porous ...
http://iranjournals.nlai.ir/1348/article_661986_49855.html
In this paper, the bioconvective flow in a porous square cavity containing oxytactic microorganism in the presence of chemical reaction is investigated. The bioconvection flow and heat transfer in porous media are formulated based on the Darcy model of Boussinesq approximation. The governing partial differential equations are solved using the Galerkin finite element method. The computational numerical results are exhibited by the streamlines, isotherms, isoconcentrations of oxygen, isoconcentrations of microorganisms, average Nusselt number, average Sherwood numbers of oxygen concentration and microorganisms. The effects of key parameters such as bioconvection Rayleigh number (Rb), chemical reaction parameter (Kr) and thermal Rayleigh number (Ra) are presented and analyzed. It can be deduced that the chemical reaction reduces the strength of isoconcentrations of both oxygen and microorganisms. It has been revealed that the chemical reaction has a greater effect on the swimming of the microorganisms, average Nusselt number, and average density number.Tue, 30 Jun 2020 19:30:00 +0100Nonlinear Buckling and Post-buckling of Shape Memory Alloy Shallow Arches
http://iranjournals.nlai.ir/1348/article_661987_49855.html
In this work, the nonlinear buckling and post-buckling behavior of shallow arches made of Shape Memory Alloy (SMA) is investigated. Arches are susceptible to large deflections, due to their slenderness, especially when the external load exceeds the serviceability limit point. Beyond this, loss of stability may occur, the famous snap-through buckling. For this reason, curved beams can be used in passive vibration control devices for seismic response mitigation, and the geometrically nonlinear analysis is needed for the accurate prediction of their response. Thus, in this research effort, the assumptions of the Euler-Bernoulli beam theory are considered, and the Von Karman strain field is employed to account for large deflections. The formulation of the problem is displacement-based regarding the axial (tangential) and transverse (normal) displacements, while the two governing equations are coupled and nonlinear. In order to introduce the SMA constitutive law, the stress-strain experimental curves described in the literature are employed together with a fiber approach at specific control cross-sections along the beam. The numerical solution of the longitudinal problem is achieved using the Analog Equation Method (AEM), a Boundary Element Method (BEM) based technique, and the iterative procedure is based on a Newton-Raphson scheme by using a displacement control algorithm to trace the fully nonlinear equilibrium path and overcome the limit points. Several representative examples are studied, not only to validate the proposed model but also to investigate the nonlinear buckling and post-buckling of SMA shallow arches.Tue, 30 Jun 2020 19:30:00 +0100Laplace Variational Iteration Method for Modified Fractional Derivatives with Non-singular Kernel
http://iranjournals.nlai.ir/1348/article_661988_49855.html
A universal approach by Laplace transform to the variational iteration method for fractional derivatives with the nonsingular kernel is presented; in particular, the Caputo-Fabrizio fractional derivative and the Atangana-Baleanu fractional derivative with the non-singular kernel is considered. The analysis elaborated for both non-singular kernel derivatives is shown the necessity of considering the modified Caputo-Fabrizio fractional derivative and the analogous modifications for the Atangana-Baleanu fractional derivative with non-singular Mittag-Leffler kernel in order to satisfy the initial conditions for some fractional differential equations.Tue, 30 Jun 2020 19:30:00 +0100The Density-Driven Nanofluid Convection in an Anisotropic Porous Medium Layer with Rotation and ...
http://iranjournals.nlai.ir/1348/article_661989_49855.html
In this study, a numerical examination of the significance of rotation and changeable gravitational field on the start of nanofluid convective movement in an anisotropic porous medium layer is shown. A model that accounts for the impact of Brownian diffusion and thermophoresis is used for nanofluid, while Darcy’s law is taken for the porous medium. The porous layer is subjected to uniform rotation and changeable downward gravitational field which fluctuates with the height from the layer by linearly or parabolic. The higher-order Galerkin technique is applied to obtain the numerical solutions. The outcomes demonstrate that the rotation parameter TD, the thermal anisotropy parameterh and the gravity variation parameter λ slow the beginning of convective motion, whereas the mechanical anisotropy parameter ξ, the nanoparticle Rayleigh-Darcy number Rnp, the modified diffusivity ratio NAnf and the modified nanofluid Lewis number Lenf quick the start of convective motion. For instance, by rising the gravity variation parameterfrom zero to 1.4, the critical nanofluid thermal Rayleigh-Darcy number Rnf,c and the critical wave numberboost maximum around 133% and 7%, respectively for linear variation of the gravity field, while it were 47% and 2.8% for parabolic variation of the gravity field. It is also observed that the system is more unstable for the parabolic variation of the gravity field.Tue, 30 Jun 2020 19:30:00 +0100