• Title/Summary/Keyword: Runge-Kutta technique

Search Result 60, Processing Time 0.021 seconds

Combined influence of slip parameter and Reynolds number on Casson nanofluid flowing in stretching cylinder

  • Jalil, Mudassar;Hussain, Muzamal;Khadimallah, Mohamed A.;Iqbal, Waheed;Loukil, Hassen;Mouldi, Abir;Mahmoud, S.R.;Tounsi, Abdelouahed
    • Computers and Concrete
    • /
    • v.30 no.5
    • /
    • pp.369-375
    • /
    • 2022
  • Current exertion reports the numerical analysis of boundary layer slip flow of Casson Nano fluid along a permeable cylinder that is stretching in exponential manner. The modeled PDEs are changed into nonlinear ODEs through appropriate nonlinear transformations. Numerical results are attained using a renowned numerical scheme shooting method with Runge-Kutta procedure of 6th-order. Influential role of relevant parameters like Reynolds, suction, Casson fluid and slip parameters on velocity profile is investigated. The effect of influence of slip parameter γ on temperature profile is seen through graph. To ensure the authenticity of numerical procedure, outcomes of some special cases of present work are compared with published work and strong agreement is noticed.

Optimum time history analysis of SDOF structures using free scale of Haar wavelet

  • Mahdavi, S.H.;Shojaee, S.
    • Structural Engineering and Mechanics
    • /
    • v.45 no.1
    • /
    • pp.95-110
    • /
    • 2013
  • In the recent decade, practical of wavelet technique is being utilized in various domain of science. Particularly, engineers are interested to the wavelet solution method in the time series analysis. Fundamentally, seismic responses of structures against time history loading such as an earthquake, illustrates optimum capability of systems. In this paper, a procedure using particularly discrete Haar wavelet basis functions is introduced, to solve dynamic equation of motion. In the proposed approach, a straightforward formulation in a fluent manner is derived from the approximation of the displacements. For this purpose, Haar operational matrix is derived and applied in the dynamic analysis. It's free-scaled matrix converts differential equation of motion to the algebraic equations. It is shown that accuracy of dynamic responses relies on, access of load in the first step, before piecewise analysis added to the technique of equation solver in the last step for large scale of wavelet. To demonstrate the effectiveness of this scheme, improved formulations are extended to the linear and nonlinear structural dynamic analysis. The validity and effectiveness of the developed method is verified with three examples. The results were compared with those from the numerical methods such as Duhamel integration, Runge-Kutta and Wilson-${\theta}$ method.

Free vibration analysis of cracked Timoshenko beams carrying spring-mass systems

  • Tan, Guojin;Shan, Jinghui;Wu, Chunli;Wang, Wensheng
    • Structural Engineering and Mechanics
    • /
    • v.63 no.4
    • /
    • pp.551-565
    • /
    • 2017
  • In this paper, an analytical approach is proposed for determining vibration characteristics of cracked non-uniform continuous Timoshenko beam carrying an arbitrary number of spring-mass systems. This method is based on the Timoshenko beam theory, transfer matrix method and numerical assembly method to obtain natural frequencies and mode shapes. Firstly, the beam is considered to be divided into several segments by spring-mass systems and support points, and four undetermined coefficients of vibration modal function are contained in each sub-segment. The undetermined coefficient matrices at spring-mass systems and pinned supports are obtained by using equilibrium and continuity conditions. Then, the overall matrix of undetermined coefficients for the whole vibration system is obtained by the numerical assembly technique. The natural frequencies and mode shapes of a cracked non-uniform continuous Timoshenko beam carrying an arbitrary number of spring-mass systems are obtained from the overall matrix combined with half-interval method and Runge-Kutta method. Finally, two numerical examples are used to verify the validity and reliability of this method, and the effects of cracks on the transverse vibration mode shapes and the rotational mode shapes are compared. The influences of the crack location, depth, position of spring-mass system and other parameters on natural frequencies of non-uniform continuous Timoshenko beam are discussed.

Numerical Simulation of Shock Wave Reflecting Patterns for Different Flow Conditions

  • Choi, Sung-Yoon;Oh, Se-Jong
    • International Journal of Aeronautical and Space Sciences
    • /
    • v.3 no.1
    • /
    • pp.74-85
    • /
    • 2002
  • The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

Effect of suction on flow of dusty fluid along exponentially stretching cylinder

  • Iqbal, Waheed;Jalil, Mudassar;Qazaq, Amjad;Khadimallah, Mohamed A.;Naeem, Muhammad N.;Hussain, Muzamal;Mahmoud, S.R.;Ghandourah, E.;Tounsi, Abdelouahed
    • Advances in nano research
    • /
    • v.10 no.3
    • /
    • pp.263-270
    • /
    • 2021
  • The present manuscript focuses the effects of suction on the flow of the dusty fluid along permeable exponentially stretching cylinder. Derived PDEs for this work are changed into ODEs by adopting right transformations. Numerical procedure is carried out for the obtained resultant equations by Shooting Technique in accordance with Runge-Kutta (RK-6) technique. Obtained results for the parameters namely, particle interaction parameter, suction parameter and Reynold number parameters are probed thoroughly. Some salient points are: (a) Fluid velocity decreases and the dust phase velocity rises for the higher values of particle interaction parameter; (b) more suction produces retarding velocities for both the phases; (c) high Reynold number slows down the fluid velocity while the speed of dust phase and (d) skin friction coefficient goes high for all these parameters.

Measurement of Joint Resistance of $(Bi,Pb)_2Sr_2Ca_2Cu_3O_x$/Ag Superconducting Tape by Field decay Technique (자장감쇠법을 이용한 $(Bi,Pb)_2Sr_2Ca_2Cu_3O_x$/Ag 초전도선재의 접합저항 측정)

  • Kim, Jung-Ho;Lee, Seung-Muk;Joo, Jin-Ho
    • Progress in Superconductivity
    • /
    • v.14 no.1
    • /
    • pp.1-10
    • /
    • 2012
  • We fabricated a closed coils by using resistive-joint method and the joint resistance of the coils were estimated by field decay technique in liquid nitrogen. We used the Runge-kutta method for the numerical analysis to calculate the decay properties. The closed coil was wound by $(Bi,Pb)_2Sr_2Ca_2Cu_3O_x$/Ag tape. Both ends the tape were overlapped and soldered to each other. The current was induced in a closed coils by external magnetic flux density. Its decay characteristic was observed by means of measuring the magnetic flux density generated by induced current at the center of the closed coil with hall sensor. The joint resistance was calculated as the ratio of the inductance of the loop to the time constants. The joint resistances were evaluated as a function of critical current of loop, contact length, sweep time, and external magnetic flux density in a contact length of 7 cm. It was observed that joint resistance was dependent on contact length of a closed coil, but independent of critical current, sweep time, and external magnetic flux density. The joint resistance was measured to be higher for a standard four-probe method, compared with that for the field decay technique. This implies that noise of measurement in a standard four-probe method is larger than that of field decay technique. It was estimated that joint resistance was $8.0{\times}10^{-9}{\Omega}$ to $11.4{\times}10^{-9}{\Omega}$ for coils of contact length for 7 cm. It was found that 40Pb/60Sn solder are unsuitable for persistent mode.

Quantification of nonlinear seismic response of rectangular liquid tank

  • Nayak, Santosh Kumar;Biswal, Kishore Chandra
    • Structural Engineering and Mechanics
    • /
    • v.47 no.5
    • /
    • pp.599-622
    • /
    • 2013
  • Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

Effect of dynamic absorber on the nonlinear vibration of SFG cylindrical shell

  • Foroutan, Kamran;Ahmadi, Habib
    • Advances in aircraft and spacecraft science
    • /
    • v.7 no.4
    • /
    • pp.291-308
    • /
    • 2020
  • In this paper, a numerical method is utilized to study the effect of a new vibration absorber on vibration response of the stiffened functionally graded (SFG) cylindrical shell under a couple of axial and transverse compressions. The material composition of the stiffeners and shell is continuously changed through the thickness. The vibration absorber consists of a mass-spring-damper system which is connected to the ground utilizing a linear local damper. To simplify, the spring element of the vibration absorber is called global potential. The von Kármán strain-displacement kinematic nonlinearity is employed in the constitutive laws of the shell and stiffeners. To consider the stiffeners in the model, the smeared stiffener technique is used. After obtaining the governing equations, the Galerkin method is applied to discretize the nonlinear dynamic equation of system. In order to find the nonlinear vibration responses, the fourth order Runge-Kutta method is utilized. The influence of the stiffeners, the dynamic absorber parameters on the vibration behavior of the SFG cylindrical shell is investigated. Also, the influences of material parameters of the system on the vibration response are examined.

IRK vs Structural Integrators for Real-Time Applications in MBS

  • Dopico D.;Lugris U.;Gonzalez M.;Cuadrado J.
    • Journal of Mechanical Science and Technology
    • /
    • v.19 no.spc1
    • /
    • pp.388-394
    • /
    • 2005
  • Recently, the authors have developed a method for real-time dynamics of multibody systems, which combines a semi-recursive formulation to derive the equations of motion in dependent relative coordinates, along with an augmented Lagrangian technique to impose the loop closure conditions. The following numerical integration procedures, which can be grouped into the so-called structural integrators, were tested : trapezoidal rule, Newmark dissipative schemes, HHT rule, and the Generalized-${\alpha}$ family. It was shown that, for large multi body systems, Newmark dissipative was the best election since, provided that the adequate parameters were chosen, excellent behavior was achieved in terms of efficiency and robustness with acceptable levels of accuracy. In the present paper, the performance of the described method in combination with another group of integrators, the Implicit Runge-Kutta family (IRK), is analyzed. The purpose is to clarify which kind of IRK algorithms can be more suitable for real-time applications, and to see whether they can be competitive with the already tested structural family of integrators. The final objective of the work is to provide some practical criteria for those interested in achieving real-time performance for large and complex multibody systems.

Marangoni convection radiative flow of dusty nanoliquid with exponential space dependent heat source

  • Mahanthesh, Basavarajappa;Gireesha, Bijjanal Jayanna;PrasannaKumara, Ballajja Chandra;Shashikumar, Nagavangala Shankarappa
    • Nuclear Engineering and Technology
    • /
    • v.49 no.8
    • /
    • pp.1660-1668
    • /
    • 2017
  • The flow of liquids submerged with nanoparticles is of significance to industrial applications, specifically in nuclear reactors and the cooling of nuclear systems to improve energy efficiency. The application of nanofluids in water-cooled nuclear systems can result in a significant improvement of their economic performance and/or safety margins. Therefore, in this paper, Marangoni thermal convective boundary layer dusty nanoliquid flow across a flat surface in the presence of solar radiation is studied. A two phase dusty liquid model is considered. Unlike classical temperature-dependent heat source effects, an exponential space-dependent heat source aspect is considered. Stretching variables are utilized to transform the prevailing partial differential system into a nonlinear ordinary differential system, which is then solved numerically via the Runge-Kutta-Fehlberg approach coupled with a shooting technique. The roles of physical parameters are focused in momentum and heat transport distributions. Graphical illustrations are also used to consider local and average Nusselt numbers. We examined the results under both linear and quadratic variation of the surface temperature. Our simulations established that the impact of Marangoni flow is useful for an enhancement of the heat transfer rate.