Route to chaos for moderate Prandtl number convection in a porous layer heated from below

The route to chaos for moderate Prandtl number gravity driven convection in porous media is analysed by using Adomian's decomposition method which provides an accurate analytical solution in terms of infinite power series. The practical need to evaluate numerical values from the infinite power series, the consequent series truncation, and the practical procedure to accomplish this task, transform the otherwise analytical results into a computational solution achieved up to a desired but finite accuracy. The solution shows a transition to chaos via a period doubling sequence of bifurcations at a Rayleigh number value far beyond the critical value associated with the loss of stability of the convection steady solution. This result is extremely distinct from the sequence of events leading to chaos in low Prandtl number convection in porous media, where a sudden transition from steady convection to chaos associated with an homoclinic explosion occurs in the neighbourhood of the critical Rayleigh number (unless mentioned otherwise by 'the critical Rayleigh number' we mean the value associated with the loss of stability of the convection steady solution). In the present case of moderate Prandtl number convection the homoclinic explosion leads to a transition from steady convection to a period-2 periodic solution in the neighbourhood of the critical Rayleigh number. This occurs at a slightly sub-critical value of Rayleigh number via a transition associated with a period-1 limit cycle which seem to belong to the sub-critical Hopf bifurcation around the point where the convection steady solution looses its stability. The different regimes are analysed and periodic windows within the chaotic regime are identified. The significance of including a time derivative term in Darcy's equation when wave phenomena are being investigated becomes evident from the results.