Sharif Digital Repository

Vibration analysis of a new type of compliant parallel mechanism with flexible intermediate links is investigated. The application of the Timoshenko beam theory to the mathematical modeling of the intermediate flexible link is described, and the equations of motion of the flexible links are obtained by using Lagrange's equation of motion. The equations of motion are obtained in the form of a set of ordinary differential equations by using assumed mode method theory. The governing differential equations of motion are solved using perturbation method. The assumed mode shapes and frequencies are to be obtained based on clamped-clamped boundary conditions. Comparing perturbation method with... 

In this work, a solution is applied to investigate the heat transfer characteristics in a pipe with turbulent decaying swirling flow by using the boundary layer integral scheme. The governing equation is solved using the forth-order Runge-Kutta scheme resulting in thermal boundary-layer thickness and dimensionless heat transfer coefficient, namely, the Nusselt number. Both forced- and free-vortex profiles are considered for the tangential velocity component. A comparison of the results obtained for the Nusselt number with available experimental data shows that this scheme has good capability in predicting the heat transfer parameters of swirling flow especially in the entrance region of a... 

In this paper, a two-link flexible manipulator is considered. For a prescribed motion, Timoshenko and Euler-Bernoulli beam models are considered. Using the Galerkin method, nonlinear equations of motion are solved. The Runge-Kutta method is employed for the time response integration method. A comparative study is made between the Euler-Bernoulli and Timoshenko beam models, with and without foreshortening effects. It is demonstrated that for two-link manipulators, both theories provide good models, and the results for both theories are very similar for all ranges of slenderness ratio. The findings suggest that for two-link manipulators with relatively high slenderness ratios, there is a... 

Background and aims of the study: Mechanical heart valves (MHV) are widely used to replace dysfunctional and failed heart valves. The bileaflet MHV design is very popular due to its superior hemodynamics. Since their introduction in 1977, the hemodynamics of bileaflet prostheses has been extensively studied. In this study the dynamic behaviour during the closing phase of a bileaflet MHV under normal physiological conditions has been investigated. Methods: Fluid analysis is based on the control volume with moving boundaries in the vicinity of the occluder. Unsteady continuity equation, unsteady momentum equation on the control volume and unsteady Bernoulli's equation have been used to... 

Theoretical and numerical analyses of rotating disks with non-uniform thickness and material properties subjected to thermo-mechanical loadings have been carried out by variable material properties (VMP), Runge-Kutta's (RK) and finite element (FE) methods. The material is assumed to be elastic-linear hardening. A power form function is used to describe the temperature gradient with the higher temperature at outer surface. Von-Mises theory has been used as failure criterion. The effects of geometry, material and thermal loading parameters as well as boundary conditions on radial, hoop and equivalent stress distributions which have not been studied in much detail in previous works have been... 

Purpose - The purpose of this paper is to study the kinetics of static recovery in cold-rolled aluminum alloy under different heating rates. Design/methodology/approach - Deformation modeling was first performed to assess the distributions of plastic strain and stress within the deformed alloy. In the next stage, thermal analysis and the rate equation of static recovery were employed to determine the progress of static recovery under non-isothermal conditions. Accordingly, a thermal finite element analysis and the Runge-Kutta method were utilized to handle the transient heat conduction and the progress of static recovery. Finally, low temperature annealing heat treatments were conducted to... 

An analytical-numerical method to determine the dynamic response of beams with various boundary conditions subjected to a moving mass under a pulsating force is explained. Governing partial differential equations of the system are changed to a convenience type of ordinary differential equations to be solved through a Runge-Kutta scheme. Pulsating force specifications influenced the dynamic response of the beam depending on the moving mass properties. Results showed the significant effect of the boundary conditions on the dynamic response of the beam, which was considered rarely in the past. Stiffening the constraints reduces the maximum stresses in the beams. Results for identical... 

In this paper, the static pull-in behavior of electrostatically actuated functionally graded (FG) micro-beams resting on an elastic medium is studied using the modified strain gradient (MSG) theory. To this end, the equilibrium equation along with classical and non-classical boundary conditions is obtained by considering the fringing field and elastic foundations effects within the principle of minimum total potential energy. Also, the elastic medium is composed of a shear layer (Pasternak foundation) and a linear normal layer (Winkler foundation). The governing differential equation is solved for cantilever and doubly fixed FG beams using an iterative numerical method. This method is a... 

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions... 

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to... 

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to... 

In this paper, distributions of stress and strain components of rotating disks with non-uniform thickness and material properties subjected to thermo-elasto-plastic loading are obtained by semi-exact method of Liao's homotopy analysis method (HAM) and finite element method (FEM). The materials are assumed to be elastic-linear strain hardening and isotropic. The analysis of rotating disk is based on Von Mises' yield criterion. A two dimensional plane stress analysis is used. The distribution of temperature is assumed to have power forms with the hotter point located at the outer surface of the disk. A mathematical technique of transformation has been proposed to solve the homotopy equations... 

In this paper, distributions of stress and strain components of rotating disks with non-uniform thickness and material properties subjected to thermo-elastic loading under different boundary conditions are obtained by semi-exact methods of Liao's homotopy analysis method (HAM), Adomian's decomposition method and He's variational iteration method (VIM). The materials are assumed to be perfectly elastic and isotropic. A two dimensional plane stress analysis is used. The distribution of temperature over the disk radius is assumed to have power forms with the higher temperature at the outer surface. The results of the three methods are compared with those obtained by Runge-Kutta's numerical... 

In this study, the numerical simulation of the fluid flow and acoustic field of a supersonic jet is performed by using high-order discretization and the vorticity confinement (VC) method on coarse grids. The three-dimensional Navier-Stokes equations are considered in the generalized curvilinear coordinate system and the high-order compact finite-difference scheme is applied for the space discretization, and the time integration is performed by the fourth-order Runge-Kutta scheme. A low-pass high-order filter is applied to stabilize the numerical solution. The non-reflecting boundary conditions are adopted for all the free boundaries, and the Kirchhoff surface integration is utilized to... 

A new perspective suitable for understanding the details of nonlinear pumping (formation of traveling shocks) inside a pressurized cavity is constructed in this paper. Full compressible axisymmetric three-dimensional Navier-Stokes equations are used as the starting point to cover all complexities of the problem that exceedingly increase for particular ranges of Mach, Reynolds and Prandtl numbers. Then a very high-order numerical method is introduced to preserve the user-defined order of accuracy for practical simulations. For removal of spurious waves, higher-order compact filters are derived. All equations are marched in time using the classical Runge-Kutta algorithm which is appropriate... 

The main goal of this study is to assess the application of high-order compact finite-difference schemes for the solution of the Euler equations in conjunction with the compressible vorticity confinement method on both uniform Cartesian and curvilinear grids. Here, the spatial discretization of the governing equations is performed by the fourth-order compact finite-difference scheme and the temporal term is discretized by the fourth-order Runge-Kutta method. To stabilize the numerical solution, appropriate dissipation terms are applied and a detail assessment is performed to study the effects of the values of confinement and dissipation coefficients on the solution to reasonably preserve the... 

In the present study, the shock-disturbances interaction in hypersonic inviscid flows with real gas effects is studied by applying a high-order accurate numerical method with the shock capturing technique. To consider real gas effects, the equilibrium air model is utilized here. The strong conservative form of the unsteady compressible Euler equations in the 2D generalized curvilinear coordinates is formulated and the resulting system of equations for the equilibrium air model is discretized by using the fifth-order finite-difference WENO scheme in space and the explicit third-order TVD Runge–Kutta scheme in time to provide a highly accurate and robust equilibrium airflow solver. The... 

In this paper, the compressible gas flow through a pipe subjected to wall heat flux in unsteady condition in the entrance region is investigated numerically. The coupled conservation equations governing turbulent compressible viscous flow in the developing region of a pipe are solved numerically under different thermal boundary conditions. The numerical procedure is a finite-volume-based finite-element method applied to unstructured grids. The convection terms are discretized by the well-defined Roe method, whereas the diffusion terms are discretized by a Galerkin finite-element formulation. The temporal terms are evaluated based on an explicit fourth-order Runge-Kutta scheme. The effect of... 

A high-order compact finite-difference scheme is applied and assessed for the numerical simulation of structural dynamics. The two-dimensional elastic stress-strain equations are considered in the generalized curvilinear coordinates and the spatial derivatives in the resulting equations are discretized by a fourth-order compact finite-difference scheme. For the time integration, an implicit second-order dual time-stepping method is utilized in which a fourth-order Runge-Kutta scheme is used to integrate in the pseudo-time level. The accuracy and robustness of the solution procedure proposed are investigated through simulating different two-dimensional benchmark test cases in structural... 

Nonlinear vibration of a three-dimensional moving gantry crane carrying a trolley hoisting a swinging object is studied in this paper.A finite element method is used to solve nonlinear coupled governing equations of the structure. A combinational technique (Newmark-Runge-Kutta) is employed for direct integration procedure. To develop a comprehensive parametric study and sensitivity analysis of the coupled nonlinear system, sequence of numerical simulations are carried out. Parametric study is directed to find out how different parameters like speed and acceleration of the trolley and gantry crane as well as the mass of the moving trolley and swinging object may affect the linear and...