Vibration effects on weak turbulent natural convection in a porous layer heated from below

Solutions for the amplitudes of stream function and temperature, as a function of time presented in this paper are evaluated via Fortran programming, by using the subroutine DVERK by Hull et al. [12] and applying a fifth- and a sixth-order Runge-Kutta-Verner numerical method of solution. The results are presented in graphs showing the amplitudes as a function of time as well as on phase diagrams with the aim of identifying the different parameter regimes and possibly identifying the regimes that enhance the heat transfer. In order to establish the stability type of the results and the type of attracting trajectory, the results were used to evaluate their corresponding Lyapunov exponents.