We study wetting front (traveling wave) solutions to the Richards equation that describe the vertical infiltration of water through one-dimensional periodically layered unsaturated soils. We prove the existence, uniqueness, and large time asymptotic stability of the traveling wave solutions under prescribed flux boundary conditions and certain constitutive conditions. The traveling waves are connections between two steady state solutions that form near the ground surface and towards the underground water table. We found a closed form expression of the wave speed. The speed of a traveling wave is equal to the ratio of the flux difference and the difference of the spatial averages of the two steady states. We give both analytical and numerical examples showing that the wave speeds in the periodic soils can be larger or smaller than those in the homogeneous soils which have the same mean diffusivity and conductivity

We analyze the power spectra and structure functions (SFs) of the temperature and radial velocity fields, calculated in the radial and azimuthal directions, in annular centrifugal Rayleigh–Bénard convection (ACRBC) for Rayleigh number Ra ∈[108,1011], Prandtl number Pr = 10.7, and inverse Rossby number Ro−1=16 using the spatial data obtained by quasi-two-dimensional direct numerical simulation. Bolgiano and Obukhov-like (BO59-like) scalings for the energy spectrum in both the azimuthal and radial directions and thermal spectrum in the azimuthal direction are observed. The range of BO59-like scaling becomes wider as Ra increases. At Ra=1011, it is found that BO59-like scaling Eu(kr)∼kr−11/5 spans nearly two decades for the energy spectrum calculated in the radial direction. Power-law fittings in the range larger than the Bolgiano scales, the scaling exponents of transverse and longitudinal velocity SFs vs the order coincide with the theoretical prediction of BO59 scaling ζ u p =3p/5 basically. The second-order temperature SFs exhibit a gradual transition from the Obukhov–Corrsin behavior at scales smaller than the Bolgiano scales to the BO59 behavior at scales larger than the Bolgiano scales. The slopes from the third to sixth-order temperature SFs are similar, which is similar to classical Rayleigh–Bénard convection and Rayleigh–Taylor turbulence. The probability density functions (p.d.f.) of temperature fluctuations δT/σT reveal the cold plumes are strong and the p.d.f. in different regions at high Ra are similar. The stronger turbulent-mixing and larger centrifugal buoyancy in ACRBC may result in the BO59-like scaling.

We perform a numerical study of the heat transfer and flow structure of Rayleigh-Benard (RB) convection in (in most cases regular) porous media, which are comprised of circular, solid obstacles located on a square lattice. This study is focused on the role of porosity in the flow properties during the transition process from the traditional RB convection with (so no obstacles included) to Darcy-type porous-media convection with approaching 0. Simulations are carried out in a cell with unity aspect ratio, for Rayleigh number from to and varying porosities , at a fixed Prandtl number , and we restrict ourselves to the two-dimensional case. For fixed , the Nusselt number is found to vary non-monotonically as a function of ; namely, with decreasing , it first increases, before it decreases for approaching 0. The non-monotonic behaviour of originates from two competing effects of the porous structure on the heat transfer. On the one hand, the flow coherence is enhanced in the porous media, which is beneficial for the heat transfer. On the other hand, the convection is slowed down by the enhanced resistance due to the porous structure, leading to heat transfer reduction. For fixed , depending on , two different heat transfer regimes are identified, with different effective power-law behaviours of versus , namely a steep one for low when viscosity dominates, and the standard classical one for large . The scaling crossover occurs when the thermal boundary layer thickness and the pore scale are comparable. The influences of the porous structure on the temperature and velocity fluctuations, convective heat flux and energy dissipation rates are analysed, further demonstrating the competing effects of the porous structure to enhance or reduce the heat transfer.

In the long-time scale, we consider the fluid dynamical limits for the kinetic equations when the fluctuation is decomposed into even and odd parts with respect to the microscopic velocity with different scalings. It is shown that when the background state is an absolute Maxwellian, the limit fluid dynamical equations are the incompressible Navier-Stokes equations with viscous heating. This is different from the case when the even and odd parts of the fluctuation have the same scaling where the standard incompressible Navier-Stokes equations without viscous heating are obtained. On the other hand, when the background is a local Maxwellian, it is shown that the above even-odd decomposition leads to a non-classical fluid dynamical system without viscous heating which has been used to describe the ghost effect in the kinetic theory. In addition, the above even-odd decomposition is justified rigorously for the Boltzmann equation for the former case when the background is an absolute Maxwellian.

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.