“…The linear instability theory gives us the critical Darcy‐Rayleigh number, but it does not predict the amplitude of the convective motion, stability of steady/oscillatory bifurcations, and heat and mass transfer. Thus, the weak nonlinear stability analysis is carried out by using the well‐known perturbation method 28,29 . In the analysis, a small bifurcation parameter that indicates deviation from the critical state is introduced and the variables are expanded in powers of in the form …”