No Arabic abstract
Scaled physical modeling is an important means to understand the behavior of fluids in nature. However, a common source of errors is conflicting similarity criteria. Here, we present using hypergravity to improve the scaling similarity of gravity-dominated fluid convection, e.g. natural convection and multi-phase flow. We demonstrate the validity of the approach by investigating water-brine buoyant jet experiments conducted under hypergravity created by a centrifuge. Results show that the scaling similarity increases with the gravitational acceleration. In particular, the model best represents the prototype under N3g with a spatial scale of 1/N and a time scale of 1/N2 by simultaneously satisfying the Froude and Reynolds criteria. The significance of centrifuge radius and fluid velocity in determining the accuracy of the scaled model is discussed in the light of Coriolis force and turbulence. This study demonstrates a new direction for the physical modeling of fluids subject to gravity with broad application prospects.
Multi-fluid models have recently been proposed as an approach to improving the representation of convection in weather and climate models. This is an attractive framework as it is fundamentally dynamical, removing some of the assumptions of mass-flux convection schemes which are invalid at current model resolutions. However, it is still not understood how best to close the multi-fluid equations for atmospheric convection. In this paper we develop a simple two-fluid, single-column model with one rising and one falling fluid. No further modelling of sub-filter variability is included. We then apply this model to Rayleigh-B{e}nard convection, showing that, with minimal closures, the correct scaling of the heat flux (Nu) is predicted over six orders of magnitude of buoyancy forcing (Ra). This suggests that even a very simple two-fluid model can accurately capture the dominant coherent overturning structures of convection.
When boiling occurs in a liquid flow field, the phenomenon is known as forced-convection boiling. We numerically investigate such a boiling system on a cylinder in a flow at a saturated condition. To deal with the complicated liquid-vapor phase-change phenomenon, we develop a numerical scheme based on the pseudopotential lattice Boltzmann method (LBM). The collision stage is performed in the space of central moments (CMs) to enhance numerical stability for high Reynolds numbers. The adopted forcing scheme, consistent with the CMs-based LBM, leads to a concise yet robust algorithm. Furthermore, additional terms required to ensure thermodynamic consistency are derived in a CMs framework. The effectiveness of the present scheme is successfully tested against a series of boiling processes, including nucleation, growth, and departure of a vapor bubble for Reynolds numbers varying between 30 and 30000. Our CMs-based LBM can reproduce all the boiling regimes, i.e., nucleate boiling, transition boiling, and film boiling, without any artificial input such as initial vapor phase. We find that the typical boiling curve, also known as the Nukiyama curve, appears even though the focused system is not the pool boiling but the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid-liquid contact even in the film-boiling regime. Finally, we provide quantitative comparison with the semi-empirical correlations for the forced-convection film boiling on a cylinder on the Nu-Ja diagram.
Reduced Order Modelling (ROM) has been widely used to create lower order, computationally inexpensive representations of higher-order dynamical systems. Using these representations, ROMs can efficiently model flow fields while using significantly lesser parameters. Conventional ROMs accomplish this by linearly projecting higher-order manifolds to lower-dimensional space using dimensionality reduction techniques such as Proper Orthogonal Decomposition (POD). In this work, we develop a novel deep learning framework DL-ROM (Deep Learning - Reduced Order Modelling) to create a neural network capable of non-linear projections to reduced order states. We then use the learned reduced state to efficiently predict future time steps of the simulation using 3D Autoencoder and 3D U-Net based architectures. Our model DL-ROM is able to create highly accurate reconstructions from the learned ROM and is thus able to efficiently predict future time steps by temporally traversing in the learned reduced state. All of this is achieved without ground truth supervision or needing to iteratively solve the expensive Navier-Stokes(NS) equations thereby resulting in massive computational savings. To test the effectiveness and performance of our approach, we evaluate our implementation on five different Computational Fluid Dynamics (CFD) datasets using reconstruction performance and computational runtime metrics. DL-ROM can reduce the computational runtimes of iterative solvers by nearly two orders of magnitude while maintaining an acceptable error threshold.
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.
One of the most relevant weather regimes in the mid-latitudes atmosphere is the persistent deviation from the approximately zonally symmetric jet to the emergence of blocking patterns. Such configurations are usually connected to exceptional local stability properties of the flow which come along with an improved local forecast skills during the phenomenon. It is instead extremely hard to predict onset and decay of blockings. Covariant Lyapunov Vectors (CLVs) offer a suitable characterization of the linear stability of a chaotic flow, since they represent the full tangent linear dynamics by a covariant basis which explores linear perturbations at all time scales. Therefore, we assess whether CLVs feature a signature of the blockings. As a first step, we examine the CLVs for a quasi-geostrophic beta-plane 2-layer model in a periodic channel baroclinically driven by a meridional temperature gradient $Delta T$. An orographic forcing enhances the emergence of localized blocked regimes. We detect the blocking events with a Tibaldi-Molteni scheme adapted to the periodic channel. When blocking occurs, the global growth rates of the fastest growing CLVs are significantly higher. Hence, against intuition, the circulation is globally more unstable in blocked phases. Such an increase in the finite time Lyapunov exponents with respect to the long term average is attributed to stronger barotropic and baroclinic conversion in the case of high temperature gradients, while for low values of Delta T, the effect is only due to stronger barotropic instability. In order to determine the localization of the CLVs we compare the meridionally averaged variance of the CLVs during blocked and unblocked phases. We find that on average the variance of the CLVs is clustered around the center of blocking. These results show that the blocked flow affects all time scales and processes described by the CLVs.