Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTwo-fluid modeling of heat transfer in flows of dense suspensions
90
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We develop a two-fluid model (TFM) for simulation of thermal transport coupled to particle migration in flows of non-Brownian suspensions. Specifically, we propose a closure relation for the inter-phase heat transfer coefficient of the TFM as a function of the particle volume fraction, particle diameter, magnitude of the particle phases shear-rate tensor, and the thermal diffusivity of the particles. The effect of shear-induced migration in the particulate phase is captured through the use of state-of-the-art rheological closures. We validate the proposed interphase heat transfer coupling by calibrating it against previous experiments in a Couette cell. We find that, when the shear rate is controlled by the rotation of the inner cylinder, the shear and thermal gradients aid each other to increase the particle migration when temperature difference between the inner and outer walls, $Delta T = T_mathrm{in} - T_mathrm{out} < 0$. Meanwhile, for $Delta T > 0$, the shear and thermal gradients oppose each other, resulting in diminished particle migration, and a more uniform distribution of the particulate phase across the gap. Within the TFM framework, we identify the origin and functional form of a thermo-rheological migration force that rationalizes our observations. We also investigate the interplay of shear and thermal gradients in the presence of recirculating regions in an eccentric Couette cell (with offset axis and rotating inner cylinder). Simulations reveal that the system Nusselt number increases with the eccentricity $E$ for $Delta T > 0$, but a maximum occurs for $Delta T < 0$ at $E = 0.4$. This observation is explained by showing that, for $E>0.4$ and $Delta T < 0$, significant flow recirculation enhances particle inhomogeneity, which in turn reduces heat transfer in the system (compared to $Delta T > 0$).
rate research
Read More
We present results of interface-resolved simulations of heat transfer in suspensions of finite-size neutrally-buoyant spherical particles for solid volume fractions up to 35% and bulk Reynolds numbers from 500 to 5600. An Immersed Boundary-Volume of Fluid method is used to solve the energy equation in the fluid and solid phase. We relate the heat transfer to the regimes of particle motion previously identified, i.e. a viscous regime at low volume fractions and low Reynolds number, particle-laden turbulence at high Reynolds and moderate volume fraction and particulate regime at high volume fractions. We show that in the viscous dominated regime, the heat transfer is mainly due to thermal diffusion with enhancement due to the particle-induced fluctuations. In the turbulent-like regime, we observe the largest enhancement of the global heat transfer, dominated by the turbulent heat flux. In the particulate shear-thickening regime, however, the heat transfer enhancement decreases as mixing is quenched by the particle migration towards the channel core. As a result, a compact loosely-packed core region forms and the contribution of thermal diffusion to the total heat transfer becomes significant once again. The global heat transfer becomes, in these flows at volume fractions larger than 25%, lower than in single-phase turbulence.
We numerically investigate turbulent Rayleigh-Benard convection within two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $Nu$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $Ra=10^8$, Prandtl number $Pr=4.38$, and Weber number $We=5$. We vary the relative thickness of the upper layer between $0.01 le alpha le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 le lambda_k le 10$. Two flow regimes are observed: In the first regime at $0.04lealphale0.96$, convective flows appear in both layers and $Nu$ is not sensitive to $alpha$. In the second regime at $alphale0.02$ or $alphage0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $Nu$ is sensitive to $alpha$. To predict $Nu$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann-Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $Nu$ and for the temperature at the interface well agree with our numerical results and previous experimental data.
A method for density-based topology optimization of heat exchangers with two fluids is proposed. The goal of the optimization process is to maximize the heat transfer from one fluid to the other, under maximum pressure drop constraints for each of the fluid flows. A single design variable is used to describe the physical fields. The solid interface and the fluid domains are generated using an erosion-dilation based identification technique, which guarantees well-separated fluids, as well as a minimum wall thickness between them. Under the assumption of laminar steady flow, the two fluids are modelled separately, but in the entire computational domain using the Brinkman penalization technique for ensuring negligible velocities outside of the respective fluid subdomains. The heat transfer is modelled using the convection-diffusion equation, where the convection is driven by both fluid flows. A stabilized finite element discretization is used to solve the governing equations. Results are presented for two different problems: a two-dimensional example illustrating and verifying the methodology; and a three-dimensional example inspired by shell-and-tube heat exchangers. The optimized designs for both cases show an improved heat transfer compared to the baseline designs. For the shell-and-tube case, the full freedom topology optimization approach is shown to yield performance improvements of up to 113% under the same pressure drop.
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 two-phase, low-Mach-number flow solver is proposed for variable-density liquid and gas with phase change. The interface is captured using a split Volume-of-Fluid method, which solves the advection of the reference phase, generalized for the case where the liquid velocity is not divergence-free and both phases exchange mass. A sharp interface is identified by using PLIC. Mass conservation is achieved in the limit of incompressible liquid, but not with the liquid compressibility and mass exchange. This is a relevant modeling choice for two-phase mixtures at near-critical and supercritical pressure conditions for the liquid but away from the mixture critical temperature. Under this thermodynamic environment, the dissolution of lighter gas species into the liquid phase is enhanced and vaporization or condensation can occur simultaneously at different interface locations. The numerical challenge of solving two-phase, supercritical-pressure flows is greater than simpler two-phase solvers because: a) local phase equilibrium is imposed at each interface cell to determine temperature, composition, or surface tension coefficient; b) a real-fluid thermodynamic model is used to obtain fluid properties; and c) necessary phase-wise values for certain variables are obtained via extrapolation techniques. To alleviate the increased numerical cost, the pressure Poisson equation (PPE) used to solve the low-Mach-number flow is split into a constant-coefficient implicit part and a variable-coefficient explicit part. Thus, a Fast Fourier Transform method can be used for the PPE. Various verification tests are performed to show the accuracy and viability of the present approach. The growth of surface instabilities in a binary system composed of liquid n-decane and gaseous oxygen at supercritical pressures for n-decane is analyzed. Other features of supercritical liquid injection are also shown.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
