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.

The minimal speeds (c) of the KolmogorovPetrovskyPiskunov (KPP) fronts at small diffusion ( 1) in a class of time-periodic cellular flows with chaotic streamlines is investigated in this paper. The variational principle of c reduces the computation to that of a principle eigenvalue problem on a periodic domain of a linear advectiondiffusion operator with spacetime periodic coefficients and small diffusion. To solve the advection dominated time-dependent eigenvalue problem efficiently over large time, a combination of spectral methods and finite element, as well as the associated fast solvers, are utilized to accelerate computation. In contrast to the scaling c= O ( 1/4) in steady cellular flows, a new relation c= O (1) as 1 is revealed in the time-periodic cellular flows due to the presence of chaotic streamlines. Residual propagation speed emerges from the Lagrangian chaos which is quantified as a sub

No abstract uploaded!

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.

We introduce a unified particle framework which integrates the phase-field method with multi-material simulation to allow modeling of both liquids and solids, as well as phase transitions between them. A simple elastoplastic model is used to capture the behavior of various kinds of solids, including deformable bodies, granular materials, and cohesive soils. States of matter or phases, particularly liquids and solids, are modeled using the non-conservative Allen-Cahn equation. In contrast, materials—made of different substances—are advected by the conservative Cahn-Hilliard equation. The distributions of phases and materials are represented by a phase variable and a concentration variable, respectively, allowing us to represent commonly observed fluid-solid interactions. Our multi-phase, multi-material system is governed by a unified Helmholtz free energy density. This framework provides the first method in computer graphics capable of modeling a continuous interface between phases. It is versatile and can be readily used in many scenarios that are challenging to simulate. Examples are provided to demonstrate the capabilities and effectiveness of this approach.