No Arabic abstract
The effects of volume exclusion and long-range intermolecular attraction are investigated by the simplified kinetic model for surface-confined inhomogeneous fluids. Gas dynamics of the ideal gas, the hard-sphere fluid and the real gas are simulated by the Boltzmann equation, the Enskog equation and the simple kinetic equation, respectively. Only the Knudsen minimum appears for the ideal gas, while both the Knudsen minimum and the Knudsen maximum occur for the hard-sphere fluid and the real gas under certain confinements, beyond which the maximum and minimum may disappear. The Boltzmann equation and the Enskog equation overestimates and underestimates the mass flow rate of the real gas dynamics under confinement, respectively, where the volume exclusion and the long-range intermolecular attractive potential among molecules are not ignorable. With the increase of the channel width, gas dynamics of the hard-sphere fluid and the real gas tends to the Boltzmann prediction gradually. The density inhomogeneity, which hinders the flow under confinement, is more obvious when the solid fraction is larger. The anomalous slip occurs for real gas under constant confinement. The flow at a smaller Knudsen number (larger solid fraction or channel width) contributes more practical amount of mass transfer, although the rarefaction effects is more prominent at larger Knudsen numbers. The temperature has no effect on density and velocity profiles of the ideal gas and the hard-sphere fluid, but the energy parameter among the real gas molecules decreases with the increasing temperature and the real gas dynamics tends to the hard-sphere ones consequently.
We study, using direct numerical simulations, the effect of geometrical confinement on heat transport and flow structure in Rayleigh-Benard convection in fluids with different Prandtl numbers. Our simulations span over two decades of Prandtl number $Pr$, $0.1 leq Pr leq 40$, with the Rayleigh number $Ra$ fixed at $10^8$. The width-to-height aspect ratio $Gamma$ spans between $0.025$ and $0.25$ while the length-to-height aspect ratio is fixed at one. We first find that for $Pr geq 0.5$, geometrical confinement can lead to a significant enhancement in heat transport as characterized by the Nusselt number $Nu$. For those cases, $Nu$ is maximal at a certain $Gamma = Gamma_{opt}$. It is found that $Gamma_{opt}$ exhibits a power-law relation with $Pr$ as $Gamma_{opt}=0.11Pr^{-0.06}$, and the maximal relative enhancement generally increases with $Pr$ over the explored parameter range. As opposed to the situation of $Pr geq 0.5$, confinement-induced enhancement in $Nu$ is not realized for smaller values of $Pr$, such as $0.1$ and $0.2$. The $Pr$ dependence of the heat transport enhancement can be understood in its relation to the coverage area of the thermal plumes over the thermal boundary layer (BL) where larger coverage is observed for larger $Pr$ due to a smaller thermal diffusivity. We further show that $Gamma_{opt}$ is closely related to the crossing of thermal and momentum BLs, and find that $Nu$ declines sharply when the thickness ratio of the thermal and momentum BLs exceeds a certain value of about one. In addition, through examining the temporally averaged flow fields and 2D mode decomposition, it is found that for smaller $Pr$ the large-scale circulation is robust against the geometrical confinement of the convection cell.
We numerically investigate the effect of non-condensable gas inside a vapor bubble on bubble dynamics, collapse pressure and pressure impact of spherical and aspherical bubble collapses. Free gas inside a vapor bubble has a damping effect that can weaken the pressure wave and enhance the bubble rebound. To estimate this effect numerically, we derive and validate a multi-component model for vapor bubbles containing gas. For the cavitating liquid and the non-condensable gas, we employ a homogeneous mixture model with a coupled equation of state for all components. The cavitation model for the cavitating liquid is a barotropic thermodynamic equilibrium model. Compressibility of all phases is considered in order to capture the shock wave of the bubble collapse. After validating the model with an analytical energy partitioning model, simulations of collapsing wall-attached bubbles with different stand-off distances are performed. The effect of the non-condensable gas on rebound and damping of the emitted shock wave is well captured.
Numerical work on shockwave/boundary-layer interactions (SBLIs) to date has largely focused on span-periodic quasi-2D configurations that neglect the influence lateral confinement has on the core flow. The present study is concerned with the effect of flow confinement on Mach 2 laminar SBLI in rectangular ducts. An oblique shock generated by a 2 degree wedge forms a conical swept SBLI with sidewall boundary layers before reflecting from the bottom wall of the domain. Multiple large regions of flow-reversal are observed on the sidewalls, bottom wall and at the corner intersection. The main interaction is found to be strongly three-dimensional and highly dependent on the geometry of the duct. Comparison to quasi-2D span-periodic simulations showed sidewalls strengthen the interaction by 31% for the baseline configuration with an aspect ratio of one. The length of the shock generator and subsequent trailing edge expansion fan position was shown to be a critical parameter in determining the central separation length. By shortening the length of the shock generator, control of the interaction and suppression of the central interaction is demonstrated. Parametric studies of shock strength and duct aspect ratio were performed to find limiting behaviours. For the largest aspect ratio of four, three-dimensionality was visible across 30% of the span width away from the wall. Topological features of the three-dimensional separation are identified and shown to be consistent with `owl-like separations of the first kind. The reflection of the initial conical swept SBLIs is found to be the most significant factor determining the flow structures downstream of the main interaction.
Truncated Taylor expansions of smooth flow maps are used in Hamiltons principle to derive a multiscale Lagrangian particle representation of ideal fluid dynamics. Numerical simulations for scattering of solutions at one level of truncation are found to produce solutions at higher levels. These scattering events to higher levels in the Taylor expansion are interpreted as modeling a cascade to smaller scales.
The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. (2014b) show that the liquid-air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier-Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square root scaling relationship. We also predict the effect of interfacial tension at the liquid-gas interface on the drop impact, showing that it mediates the lift-off behavior.