No Arabic abstract
In a recent paper (El Omari and Le Guer, IJHMT, 53, 2010) we have investigated mixing and heat transfer enhancement in a mixer composed of two circular rods maintained vertically in a cylindrical tank. The rods and tank can rotate around their revolution axes while their surfaces were maintained at a constant temperature. In the present study we investigate the differences in the thermal mixing process arising from the utilization of a constant heat flux as a boundary condition. The study concerns a highly viscous fluid with a high Prandtl number $Pr = 10,000$ for which this chaotic mixer is suitable. Chaotic flows are obtained by imposing temporal modulations of the rotational velocities of the walls. By solving numerically the flow and energy equations, we studied the effects of different stirring protocols and flow configurations on the efficiency of mixing and heat transfer. For this purpose, we used different statistical indicators as tools to characterize the evolution of the fluid temperature and its homogenization. Fundamental differences have been reported between these two modes of heating or cooling: while the mixing with an imposed temperature results in a homogeneous temperature field, with a fixed heat flux we observe a constant difference between the maximal and minimal temperatures that establish in the fluid; the extent of this difference is governed by the efficiency of the mixing protocol.
A homogenization approach is proposed for the treatment of porous wall boundary conditions in the computation of compressible viscous flows. Like any other homogenization approach, it eliminates the need for pore-resolved fluid meshes and therefore enables practical flow simulations in computational fluid domains with porous wall boundaries. Unlike alternative approaches however, it does not require prescribing a mass flow rate and does not introduce in the computational model a heuristic discharge coefficient. Instead, it models the inviscid flux through a porous wall surrounded by the flow as a weighted average of the inviscid flux at an impermeable surface and that through pores. It also introduces a body force term in the governing equations to account for friction loss along the pore boundaries. The source term depends on the thickness of the porous wall and the concept of an equivalent single pore. The feasibility of the latter concept is demonstrated using low-speed permeability test data for the fabric of the Mars Science Laboratory parachute canopy. The overall homogenization approach is illustrated with a series of supersonic flow computations through the same fabric and verified using supersonic, pore-resolved numerical simulations.
Understanding the generation mechanism of the heating flux is essential for the design of hypersonic vehicles. We proposed a novel formula to decompose the heat flux coefficient into the contributions of different terms by integrating the conservative equation of the total energy. The reliability of the formula is well demonstrated by the direct numerical simulation results of a hypersonic transitional boundary layer. Through this formula, the exact process of the energy transport in the boundary layer can be explained and the dominant contributors to the heat flux can be explored, which are beneficial for the prediction of the heat and design of the thermal protection devices
We report an experimental study aiming to clarify the role of boundary conditions (BC) in high Rayleigh number $10^8 < {rm{Ra}} < 3 times 10^{12}$ turbulent thermal convection of cryogenic helium gas. We switch between BC closer to constant heat flux (CF) and constant temperature (CT) applied to the highly conducting bottom plate of the aspect ratio one cylindrical cell 30 cm in size, leading to dramatic changes in the temperature probability density function and in power spectral density of the temperature fluctuations measured at the bottom plate, while the dynamic thermal behaviour of the top plate and bulk convective flow remain unaffected. Within our experimental accuracy, we find no appreciable changes in Reynolds number Re(Ra) scaling, in the dimensionless heat transfer efficiency expressed via Nusselt number Nu(Ra) scaling, nor in the rate of direction reversals of large scale circulation.
We present an effective thermal open boundary condition for convective heat transfer problems on domains involving outflow/open boundaries. This boundary condition is energy-stable, and it ensures that the contribution of the open boundary will not cause an ``energy-like temperature functional to increase over time, irrespective of the state of flow on the open boundary. It is effective in coping with thermal open boundaries even in flow regimes where strong vortices or backflows are prevalent on such boundaries, and it is straightforward to implement. Extensive numerical simulations are presented to demonstrate the stability and effectiveness of our method for heat transfer problems with strong vortices and backflows occurring on the open boundaries. Simulation results are compared with previous works to demonstrate the accuracy of the presented method.
A series of direct numerical simulations of Rayleigh-Benard convection, the flow in a fluid layer heated from below and cooled from above, were conducted to investigate the effect of mixed insulating and conducting boundary conditions on convective flows. Rayleigh numbers between $text{Ra}=10^7$ and $text{Ra}=10^9$ were considered, for Prandtl numbers $text{Pr}=1$ and $text{Pr}=10$. The bottom plate was divided into patterns of conducting and insulating stripes. The size ratio between these stripes was fixed to unity and the total number of stripes was varied. Global quantities such as the heat transport and average bulk temperature and local quantities such as the temperature just below the insulating boundary wall were investigated. For the case with the top boundary divided into two halves, one conducting and one insulating, the heat transfer was found to be approximately two thirds of the fully conducting case. Increasing the pattern frequency increased the heat transfer which asymptotically approached the fully conducting case, even if only half of the surface is conducting. Fourier analysis of the temperature field revealed that the imprinted pattern of the plates is diffused in the thermal boundary layers, and cannot be detected in the bulk. With conducting-insulating patterns on both plates, the trends previously described were similar, however, the half-and-half division led to a heat transfer of about a half of the fully conducting case instead of two-thirds. The effect of the ratio of conducting and insulating areas was also analyzed, and it was found that even for systems with a top plate with only $25%$ conducting surface, heat-transport of $60%$ of the fully conducting case can be seen. Changing the 1D stripe pattern to 2D checkerboard tessellations does not result in a significantly different response of the system.