No Arabic abstract
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
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.
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.
We present a numerical study of quasistatic magnetoconvection in a cubic Rayleigh-Benard (RB) convection cell subjected to a vertical external magnetic field. For moderate values of the Hartmann number Ha, we find an enhancement of heat transport. Furthermore, a maximum heat transport enhancement is observed at certain optimal $Ha_{opt}$. The enhanced heat transport may be understood as a result of the increased coherency of the thermal plumes, which are elementary heat carriers of the system. To our knowledge this is the first time that a heat transfer enhancement by the stabilising Lorentz force in quasistatic magnetoconvection has been observed. We further found that the optimal enhancement may be understood in terms of the crossing between the thermal and the momentum boundary layers (BL) and the fact that temperature fluctuations are maximum near the position where the BLs cross. These findings demonstrate that the heat transport enhancement phenomenon in the quasistatic magnetoconvection system belongs to the same universality class of stabilising$-$destabilising ($S$-$D$) turbulent flows as the systems of confined Rayleigh-Benard (CRB), rotating Rayleigh-Benard (RRB) and double-diffusive convection (DDC). This is further supported by the findings that the heat transport, boundary layer ratio and the temperature fluctuations in magnetoconvection at the boundary layer crossing point are similar to the other three cases.
A practical application of universal wall scalings is near-wall turbulence modeling. In this paper, we exploit temperatures semi-local scaling [Patel, Boersma, and Pecnik, {Scalar statistics in variable property turbulent channel flows}, Phys. Rev. Fluids, 2017, 2(8), 084604] and derive an eddy conductivity closure for wall-modeled large-eddy simulation of high-speed flows. We show that while the semi-local scaling does not collapse high-speed direct numerical simulation (DNS) data, the resulting eddy conductivity and the wall model work fairly well. The paper attempts to answer the following outstanding question: why the semi-local scaling fails but the resulting eddy conductivity works well. We conduct DNSs of Couette flows at Mach numbers from $M=1.4$ to 6. We add a source term in the energy equation to get a cold, a close-to-adiabatic wall, and a hot wall. Detailed analysis of the flows energy budgets shows that aerodynamic heating is the answer to our question: aerodynamic heating is not accounted for in Patel et al.s semi-local scaling but is modeled in the equilibrium wall model. We incorporate aerodynamic heating in semi-local scaling and show that the new scaling successfully collapses the high-speed DNS data. We also show that incorporating aerodynamic heating or not, the semi-local scaling gives rise to the exact same eddy conductivity, thereby answering the outstanding question.
A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement framework of AMReX. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds $-$ acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries $-$ an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating the conservative nature of the moving boundary algorithm. Viscous test cases involve that of a horizontally moving cylinder at $Re=40$, an inline oscillating cylinder at $Re=100$, and a transversely oscillating cylinder at $Re=185$. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.