Maclaurin-series Assisted Linear and Nonlinear Stability Analyses of Rayleigh-Bénard Convection with Variable Gravity Fields

Document Type : Research Paper

Authors
S. Noor Arshika 1 , 2
P.G. Siddheshwar 2
S. Tarannum 3
D. Laroze 4
1 Department of Mathematics, Dayananda Sagar College of Engineering, Kumaraswamy Layout, Bengaluru, 560078, India
2 Centre for Mathematical Needs, Department of Mathematics, Christ University, Hosur road, Bengaluru, 560029, India
3 Department of Professional Studies, Christ University, Hosur road, Bengaluru, 560029, India
4 Instituto de Alta Investigaciόn, Universidad de Tarapacá, Casilla 7 D, Chile
10.22055/jacm.2025.48418.5220
Abstract
The effect of variable gravity fields that vary through the height of a viscous fluid layer in a Rayleigh-Bénard convection (RBC) is investigated in the paper. A minimal Fourier-series expansion leads to a boundary eigenvalue problem with variable coefficients. Using the Maclaurin-series approach, the recurrence relations for six cases of gravity fields are generated. The eigenvalue of the problem is located using the Newton-Raphson method with an error tolerance of 10-8. The main novelty of the present work is studying the influence of variable gravity fields on the nonlinear dynamics of the problem. A weakly nonlinear stability analysis is performed by first identifying the convective mode and by further arriving at the scaled Lorenz model. Comparison is made with the results of three boundary conditions viz., free-free, rigid-free and rigid-rigid. In the absence of the gravity variation parameter, the results of RBC with constant gravity are recovered. The influence of varying gravity fields on the dynamics is studied using indicators: rH-plots, bifurcation-diagram, periodicity-diagram and trapping-region of the trajectories. It is found that the effect of increasing the strength of gravity is to delay the appearance of chaos for all cases of gravity fields except for the positively linear case. Furthermore, by varying the gravity fields, one can witness shifts in the chaotic and periodic regimes. Moreover, the first largest periodic burst almost overlaps for cubic and biquadratic cases of gravity fields. The shrinking in the size of the ellipsoid with change in the gravity fields is highlighted in the paper.
Keywords
Variable gravity
Convection
Boundary conditions
Chaos
Maclaurin-series
Newton-Raphson
Runge-Kutta
Subjects

Publisher’s Note Shahid Chamran University of Ahvaz remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

[1] Chen, C.F., Chen, F., Onset of salt finger convection in a gravity gradient, Physics of Fluids, 4, 1992, 451-452.
[2] Hirt, C., Claessens, C.S., Fecher, T., Kuhn, M., Pail, R., Rexer, M., New ultrahigh-resolution picture of earth’s gravity field, Geophysical Research Letters, 40, 2013, 4279-4283.
[3] Tapley, B.D., Bettadpur, S., Ries, J.C., Thompson, P.F., Watkins, M.M., GRACE measurements of mass variability in the earth system, Science, 305, 2004, 503-505.
[4] Pradhan, G.K., Samal, P.C., Thermal stability of a fluid layer under variable body forces, Journal of Mathematical Analysis and Applications, 122, 1987, 487-495.
[5] Straughan, B., Convection in a variable gravity field, Journal of Mathematical Analysis and Applications, 140, 1989, 467-475.
[6] Rionero, S., Straughan, B., Convection in a porous medium with internal heat source and variable gravity effects, International Journal of Engineering Science, 28, 1990, 497-503.
[7] Kaloni, P.N., Qiao, Z., Non-linear convection in a porous medium with inclined temperature gradient and variable gravity effects, International Journal of Heat and Mass Transfer, 44, 2001, 1585-1591.
[8] Straughan, B., The Energy Method, Stability, and Nonlinear Convection, Springer, 2004.
[9] Arnone, G., Capone, F., Iovanna, F., Massa G., Variable gravity effects on penetrative porous convection, International Journal of Non-Linear Mechanics, 158, 2024, 104579.
[10] Gheorghiu, C.I., Dragomirescu, F.I., Spectral methods in linear stability. Applications to thermal convection with variable gravity field, Applied Numerical Mathematics, 59, 2009, 1290-1302.
[11] Mahajan, A., Sharma, M.K., The onset of convection in a magnetic nanofluid layer with variable gravity effects, Applied Mathematics and Computation, 339, 2018, 622-635.
[12] Yadav, D., Numerical investigation of the combined impact of variable gravity field and throughflow on the onset of convective motion in a porous medium layer, International Communications in Heat and Mass Transfer, 108, 2019, 104274.
[13] Mahabaleshwar, U.S., Basavaraja, D., Wang, S., Lorenzini, G., Lorenzini, E., Convection in a porous medium with variable internal heat source and variable gravity, International Journal of Heat and Mass Transfer, 111, 2017, 651-656.
[14] Saltzman, B., Finite amplitude free convection as an initial value problem-I, Journal of the Atmospheric Sciences, 19, 1962, 329-341.
[15] Lorenz, E.N., Deterministic non periodic flow, Journal of the Atmospheric Sciences, 20, 1963, 130-141.
[16] Laroze, D., Siddheshwar, P.G., Pleiner, H., Chaotic convection in a ferrofluid, Communications in Nonlinear Science and Numerical Simulation, 18, 2013, 2436-2447.
[17] Siddheshwar, P.G., Laroze, D., Kanchana, C., Two-dimensional Rayleigh-Bénard convection of viscoelastic liquids in Cartesian and cylindrical coordinates: regular and chaotic regimes, The European Physical Journal Plus, 137, 2022, 922.
[18] Barman, D., Srinivasacharya, D., Effect of variable gravity field on the onset of convection in a Brinkman porous medium under convective boundary conditions, International Communications in Heat and Mass Transfer, 144, 2023, 106777.
[19] Yadav, D., Chu, Y.M., Li, Z., Examination of the nanofluid convective instability of vertical constant throughflow in a porous medium layer with variable gravity, Applied Nanoscience, 13, 2023, 353-366.
[20] Yadav, D., The Density-Driven Nanofluid Convection in an Anisotropic Porous Medium Layer with Rotation and Variable Gravity Field: A Numerical Investigation, Journal of Applied and Computational Mechanics, 6(3), 2020, 699-712.
[21] Kanchana, C., Siddheshwar, P.G., Zhao, Y., The effect of boundary conditions on the onset of chaos in Rayleigh-Bénard convection using energy conserving Lorenz models, Applied Mathematical Modelling, 88, 2020, 349-366.
[22] Kanchana, C., Siddheshwar, P.G., Pérez, L.M., Laroze, D., Comparison of the effect of suction-injection-combination on Rayleigh–Bénard convection in the case of asymmetric boundaries with those of symmetric ones, Physics of Fluids, 35(5), 2023, 053615.
[23] Hart, J.E., Toomre, J., Deane, A.E., Hurlburt, N.E., Glatzmaier, G.A., Fichtl, G.H., Leslie, F., Fowlis, W.W., Gilman, P.A., Laboratory experiments on planetary and stellar convection performed on Spacelab 3, NASA Technical Reports Server, 1986.
[24] Rodot, H., Regel, L.L., Turtchaninov, A.M., Crystal growth of IV–VI semiconductors in a centrifuge, Journal of Crystal Growth, 104, 1990, 280-284.
Journal of Applied and Computational Mechanics
Volume 12, Issue 2
Spring 2026
Pages 502-514