Flow reversals in two-dimensional Rayleigh-Benard convection led by non-Oberbeck-Boussinesq (NOB) effects due to large temperature differences arc studied by direct numerical simulation. Perfect gas is chosen as the working fluid and the Prandtl number is 0.71 for the reference state. 11 NOB effects are included, the flow pattern P-11 with only one dominant roll often becomes unstable by the growth of the cold corner roll, which sometimes results in cession-led flow reversals. By exploiting the vorticity transport equation, it is found that the asymmetries of buoyancy and viscous forces are responsible for the growth of the cold corner roll because both such asymmetries cause an imbalance between the corner rolls and the large-scale circulation (I.SC). The buoyancy force near the cold wall increases and decreases near the hot wall originating from the temperature-dependent isobaric thermal expansion coefficient alpha = 1/T if NOB effects are included; Moreover, the decreased dissipation due to lower viscosity is favourable for the growth of the cold corner roll, while the increased viscosity further suppresses the growth of the hot corner roll. Finally, it is found that the boundary layer near the cold wall plays an important role in the mass transport from LSC to corner rolls subject to mass conservation.

We report on a three-dimensional direct numerical simulation study of flow structure and heat transport in the annular centrifugal Rayleigh-Bdnard convection (ACRBC) system, with cold inner and hot outer cylinders corotating axially, for the Rayleigh number range Ra is an element of [10(6), 10(8)] and radius ratio range eta = R-i/R-o is an element of [0.3, 0.9] (R-i and R-o are the radius of the inner and outer cylinders, respectively). This study focuses on the dependence of flow dynamics, heat transport and asymmetric mean temperature fields on the radius ratio eta. For the inverse Rossby number Ro(-1) = 1, as the Coriolis force balances inertial force, the flow is in the inertial regime. The mechanisms of zonal flow revolving in the prograde direction in this regime are attributed to the asymmetric movements of plumes and the different curvatures of the cylinders. The number of roll pairs is smaller than the circular roll hypothesis as the convection rolls are probably elongated by zonal flow. The physical mechanism of zonal flow is verified by the dependence of the drift frequency of the large-scale circulation (LSC) rolls and the space- and time-averaged azimuthal velocity on eta. The larger eta is, the weaker the zonal flow becomes. We show that the heat transport efficiency increases with eta. It is also found that the bulk temperature deviates from the arithmetic mean temperature and the deviation increases as eta decreases. This effect can be explained by a simple model that accounts for the curvature effects and the radially dependent centrifugal force in ACRBC.

A key difficulty in the analysis and numerical approximation of the shallow water equations is the nonconservative product of measures due to the gravitational force acting on a sloped bottom. Solutions may be nonunique, and numerical schemes are not only consistent discretizations of the shallow water equations, but they also determine how to model the physics. Our derivation is based on a subcell reconstruction using infinitesimal singular layers at the cell boundaries, as inspired by S. Noelle, Y. Xing, and C.-W. Shu [J. Comput. Phys., 226 (2007), pp. 29-58]. One key step is to separate the singular measures. Another aspect is the reconstruction of the solution variables in the singular layers. We study three reconstructions. The first leads to the well-known scheme of Audusse et al., [SIAM J. Sci. Comput., 25 (2004), pp. 2050-2065], which introduces the hydrostatic reconstruction. The second is a modification proposed in [T. Morales de Luna, M. J. Castro D´ıaz, and C. Par´es, Appl. Math. Comput., 219 (2013), pp. 9012-9032], which analyzes whether a wave has enough energy to overcome a step. The third is our new scheme, which borrows its structure from the wet-dry front. For a number of cases discussed in recent years, where water runs down a hill, Audusse’s scheme converges slowly or fails. Morales’ scheme gives a visible improvement. Both schemes are clearly outperformed by our new scheme.

We report a numerical study of Rayleigh-Benard convection through random porous media using pore-scale modelling, focusing on the Lagrangian dynamics of fluid particles and heat transfer for varied porosities . Due to the interaction between the porous medium and the coherent flow structures, the flow is found to be highly heterogeneous, consisting of convection channels with strong flow strength and stagnant regions with low velocities. The modifications of flow field due to porous structure have a significant influence on the dynamics of fluid particles. Evaluation of the particle displacement along the trajectory reveals the emergence of anomalous transport for long times as is decreased, which is associated with the long-time correlation of Lagrangian velocity of the fluid. As porosity is decreased, the cross-correlation between the vertical velocity and temperature fluctuation is enhanced, which reveals a mechanism to enhance the heat transfer in porous-media convection.

The shape stability of the reaction interface for reactive flow in a layered porous medium is studied. This is done using a complete linearized stability analysis in the setting of a free boundary model of this phenomenon. The spectrum of the linearized problem is obtained in the general case, and it is compared with that obtained for homogeneous media in two typical cases.