# American Institute of Mathematical Sciences

## Transitions and bifurcations of Darcy-Brinkman-Marangoni convection

 1 School of Mathematics, Southwest Jiaotong University, Chengdu, 610031, China 2 College of Mathematics and Statistics, Shenzhen University, Shenzhen, 518060, China 3 College of Mathematics, Sichuan University, Chengdu, 610065, China

* Corresponding author: Yiqiu Mao

Received  January 2021 Published  March 2021

This study examines dynamic transitions of Brinkman equation coupled with the thermal diffusion equation modeling the surface tension driven convection in porous media. First, we show that the equilibrium of the equation loses its linear stability if the Marangoni number is greater than a threshold, and the corresponding principle of exchange stability (PES) condition is then verified. Second, we establish the nonlinear transition theorems describing the possible transition types associated with the linear instability of the equilibrium by applying the center manifold theory to reduce the infinite dynamical system to a finite dimensional one together with several non-dimensional transition numbers. Finally, careful numerical computations are performed to determine the sign of these transition numbers as well as related transition types. Our result indicates that the system favors all three types of transitions. Unlike the buoyancy forces driven convections, jump and mixed type transition can occur at certain parameter regimes.

Citation: Zhigang Pan, Yiqiu Mao, Quan Wang, Yuchen Yang. Transitions and bifurcations of Darcy-Brinkman-Marangoni convection. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021106
##### References:

show all references

##### References:
Plot of critical porous Marangoni number $\text{Ma}^*$ as a function of $a\in[3, 5]$ and $b\in[3, 5]$ with $\text{Bi} = 2$ (left) and $\text{Bi} = 10$ (right) and $\text{Da} = 0.1$
Plot of $(m_c, n_c)$ as a function of $a\in[3, 4]$ and $b\in[3, 4]$ such that $Ma^* = \mathcal{M}(m_c, n_c)$, where $\text{Bi} = 2$ and $\text{Da} = 0.1$
Comparison between exact values of Marangoni number and numerical predictions, where $\text{Bi} = 2$, $\text{Da} = 0.1$, and $\lambda = 1000$
 $a, b$ Exact $Ma^*$ Numerical prediction Relative error $a=3, b=3$ 209.82420647 209.82420611 $3.6\times 10^{-7}$ $a=3.2, b=3.2$ 210.25798681 210.25798320 $3.6\times 10^{-6}$ $a=3.4, b=3.4$ 208.54692216 208.54692182 $3.4\times 10^{-7}$ $a=3.6, b=3.6$ 208.51417969 208.51417509 $4.6\times 10^{-6}$
 $a, b$ Exact $Ma^*$ Numerical prediction Relative error $a=3, b=3$ 209.82420647 209.82420611 $3.6\times 10^{-7}$ $a=3.2, b=3.2$ 210.25798681 210.25798320 $3.6\times 10^{-6}$ $a=3.4, b=3.4$ 208.54692216 208.54692182 $3.4\times 10^{-7}$ $a=3.6, b=3.6$ 208.51417969 208.51417509 $4.6\times 10^{-6}$
Numerical predictions of the sign of the transition number $r$, where $\text{Bi} = 2$, $\text{Da} = 0.1$, and $\lambda = 1000$
 $a, \; b$ $(m_c, n_c)$ sign$(r)$ $a=3.3, \; b=3.3$ $(3, 1)$ $1$ $a=3.35, \; b=3.3$ $(3, 1)$ $1$ $a=3.4, b=3.3$ $(1, 3)$ $-1$ $a=3.45, \; b=3.3$ $(1, 3)$ $-1$
 $a, \; b$ $(m_c, n_c)$ sign$(r)$ $a=3.3, \; b=3.3$ $(3, 1)$ $1$ $a=3.35, \; b=3.3$ $(3, 1)$ $1$ $a=3.4, b=3.3$ $(1, 3)$ $-1$ $a=3.45, \; b=3.3$ $(1, 3)$ $-1$
Numerical predictions of the sign of the transition number $r$, where $\text{Bi} = 2$, $\text{Da} = 0.1$, and $\lambda = 1000$
 $a, b$ $(m_c, n_c)$ sign$(r)$ $a=3.4, b=3.75$ $(0, 3)$ $-1$ $a=3.45, b=3.75$ $(0, 3)$ $-1$ $a=3.5, b=3.75$ $(0, 3)$ $-1$ $a=3.55, b=3.75$ $(0, 3)$ $-1$
 $a, b$ $(m_c, n_c)$ sign$(r)$ $a=3.4, b=3.75$ $(0, 3)$ $-1$ $a=3.45, b=3.75$ $(0, 3)$ $-1$ $a=3.5, b=3.75$ $(0, 3)$ $-1$ $a=3.55, b=3.75$ $(0, 3)$ $-1$
Numerical predictions of the sign of the transition numbers $r_1$ and $S_2$, where $\text{Bi} = 2$, $\text{Da} = 0.1$, and $\lambda = 1000$
 $a, b$ $(m_c, n_c)$ sign$(r_1)$ sign$(s_2)$ $a=1.4300024, b=2.47683688$ $(1, 1), (0, 2)$ $1$ $1$ $a=1.4300024, b= 4.95367376$ $(1, 2), (0, 4)$ $1$ $1$
 $a, b$ $(m_c, n_c)$ sign$(r_1)$ sign$(s_2)$ $a=1.4300024, b=2.47683688$ $(1, 1), (0, 2)$ $1$ $1$ $a=1.4300024, b= 4.95367376$ $(1, 2), (0, 4)$ $1$ $1$
Numerical predictions of the sign of the transition numbers $r_1$ and $S_2$, where $\text{Bi} = 10$, $\text{Da} = 0.1$, and $\lambda = 1000$
 $a, b$ $(m_c, n_c)$ sign$(r_1)$ sign$(s_2)$ $a=1.21836059, b= 2.11026245$ $(1, 1), (0, 2)$ $1$ $-1$ $a=1.21836059, b= 4.2205249$ $(1, 2), (0, 4)$ $1$ $-1$
 $a, b$ $(m_c, n_c)$ sign$(r_1)$ sign$(s_2)$ $a=1.21836059, b= 2.11026245$ $(1, 1), (0, 2)$ $1$ $-1$ $a=1.21836059, b= 4.2205249$ $(1, 2), (0, 4)$ $1$ $-1$
 [1] Julian Tugaut. Captivity of the solution to the granular media equation. Kinetic & Related Models, 2021, 14 (2) : 199-209. doi: 10.3934/krm.2021002 [2] Wenjing Liu, Rong Yang, Xin-Guang Yang. Dynamics of a 3D Brinkman-Forchheimer equation with infinite delay. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021052 [3] Simone Cacace, Maurizio Falcone. A dynamic domain decomposition for the eikonal-diffusion equation. Discrete & Continuous Dynamical Systems - S, 2016, 9 (1) : 109-123. doi: 10.3934/dcdss.2016.9.109 [4] Yohei Yamazaki. Center stable manifolds around line solitary waves of the Zakharov–Kuznetsov equation with critical speed. Discrete & Continuous Dynamical Systems, 2021, 41 (8) : 3579-3614. doi: 10.3934/dcds.2021008 [5] Yves Capdeboscq, Shaun Chen Yang Ong. Quantitative jacobian determinant bounds for the conductivity equation in high contrast composite media. Discrete & Continuous Dynamical Systems - B, 2020, 25 (10) : 3857-3887. doi: 10.3934/dcdsb.2020228 [6] Shuting Chen, Zengji Du, Jiang Liu, Ke Wang. The dynamic properties of a generalized Kawahara equation with Kuramoto-Sivashinsky perturbation. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021098 [7] Takeshi Saito, Kazuyuki Yagasaki. Chebyshev spectral methods for computing center manifolds. Journal of Computational Dynamics, 2021  doi: 10.3934/jcd.2021008 [8] M. Phani Sudheer, Ravi S. Nanjundiah, A. S. Vasudeva Murthy. Revisiting the slow manifold of the Lorenz-Krishnamurthy quintet. Discrete & Continuous Dynamical Systems - B, 2006, 6 (6) : 1403-1416. doi: 10.3934/dcdsb.2006.6.1403 [9] Matthias Ebenbeck, Harald Garcke, Robert Nürnberg. Cahn–Hilliard–Brinkman systems for tumour growth. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021034 [10] Chris Guiver, Nathan Poppelreiter, Richard Rebarber, Brigitte Tenhumberg, Stuart Townley. Dynamic observers for unknown populations. Discrete & Continuous Dynamical Systems - B, 2021, 26 (6) : 3279-3302. doi: 10.3934/dcdsb.2020232 [11] Matthias Erbar, Jan Maas. Gradient flow structures for discrete porous medium equations. Discrete & Continuous Dynamical Systems, 2014, 34 (4) : 1355-1374. doi: 10.3934/dcds.2014.34.1355 [12] Charles Amorim, Miguel Loayza, Marko A. Rojas-Medar. The nonstationary flows of micropolar fluids with thermal convection: An iterative approach. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2509-2535. doi: 10.3934/dcdsb.2020193 [13] Zhigang Pan, Chanh Kieu, Quan Wang. Hopf bifurcations and transitions of two-dimensional Quasi-Geostrophic flows. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021025 [14] Alexandr Mikhaylov, Victor Mikhaylov. Dynamic inverse problem for Jacobi matrices. Inverse Problems & Imaging, 2019, 13 (3) : 431-447. doi: 10.3934/ipi.2019021 [15] Manoel J. Dos Santos, Baowei Feng, Dilberto S. Almeida Júnior, Mauro L. Santos. Global and exponential attractors for a nonlinear porous elastic system with delay term. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2805-2828. doi: 10.3934/dcdsb.2020206 [16] Guillaume Bal, Wenjia Jing. Homogenization and corrector theory for linear transport in random media. Discrete & Continuous Dynamical Systems, 2010, 28 (4) : 1311-1343. doi: 10.3934/dcds.2010.28.1311 [17] Huy Dinh, Harbir Antil, Yanlai Chen, Elena Cherkaev, Akil Narayan. Model reduction for fractional elliptic problems using Kato's formula. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021004 [18] Xiaofei Liu, Yong Wang. Weakening convergence conditions of a potential reduction method for tensor complementarity problems. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021080 [19] Palash Sarkar, Subhadip Singha. Classical reduction of gap SVP to LWE: A concrete security analysis. Advances in Mathematics of Communications, 2021  doi: 10.3934/amc.2021004 [20] Marian Gidea, Rafael de la Llave, Tere M. Seara. A general mechanism of instability in Hamiltonian systems: Skipping along a normally hyperbolic invariant manifold. Discrete & Continuous Dynamical Systems, 2020, 40 (12) : 6795-6813. doi: 10.3934/dcds.2020166

2019 Impact Factor: 1.27