Friction factor analysis for a nanofluid circulating in a microchannel filled with a homogeneous porous medium

This work presents the numerical solution for different velocity profiles and friction factors on a rectangular porous microchannel fully saturated by the flow of a nanofluid introducing different viscosity models, including one nanofluid density model. The Darcy-Brinkman-Forchheimer equation was used to solve the momentum equation in the porous medium. The results show that the relative density of the fluid, the nanoparticle diameters and their volumetric concentration have a direct influence on the velocity profiles only when the inertial effects caused by the presence of the porous matrix are important. Finally, it was found that only viscosity models that depend on temperature and nanoparticle diameter reduce the friction factor by seventy percent compared to a base fluid without nanoparticles; furthermore, these models show a velocity reduction of even ten percent along the symmetry axis of the microchannel.

Large Eddy simulation of tracer gas dispersion in a cavity

This paper assesses the prediction of inert tracer gas dispersion within a cavity of height (H) 1.0 m, and unity aspect ratio, using large Eddy simulation (LES). The flow Reynolds number was 67 000, based on the freestream velocity and cavity height. The flow upstream of the cavity was laminar, producing a cavity shear layer which underwent a transition to turbulence over the cavity. Three distinct meshes are used, with grid spacings of H / 100 (coarse), H / 200 (intermediate), and H / 400 (fine) respectively. The Smagorinsky, WALE, and Germano-Lilly subgrid-scale models are used on each grid to quantify the effects of subgrid-scale modelling on the simulated flow. Coarsening the grid led to small changes in the predicted velocity field, and to substantial over-prediction of the tracer gas concentration statistics. Quantitative metric analysis of the tracer gas statistics showed that the coarse grid simulations yielded results outside of acceptable tolerances, while the intermediate and fine grids produced acceptable output. Interrogation of the fluid dynamics present in each simulation showed that the evolution of the cavity shear layer is heavily influenced by the grid and subgrid scale model. On the coarse and intermediate grids the development of the shear layer is delayed, inhibiting the entrainment and mixing of the tracer gas into the shear layer, reducing the removal of the tracer gas from the cavity. On the fine grid, the shear layer developed more rapidly, resulting in enhanced removal of the tracer gas from the cavity. Concentration probability density functions showed that the fine grid simulations accurately predicted the range, and the most probable value, of the tracer gas concentration towards both walls of the cavity. The results presented in this paper show that the WALE and Germano-Lilly models may be advantageous over the standard Smagorinsky model for simulations of pollutant dispersion in the urban environment.

The effects of surface roughness on the flow in multiple connected fractures

We present a novel computationally efficient approach for investigating the effect of surface roughness on the fluid flow in discrete fracture networks at low Reynolds number. The effect of parallel and series fracture arrangements on the flow rate and hydraulic resistance was studied numerically by patching Hele-Shaw (HS) cells to represent the network. In this analysis, the impact of surface roughness was studied in different arrangements of the network. For this aim, four models with different sequences of fracture connections were studied. The validity of the models was assessed by comparing the results with solutions of the full Navier-Stokes equations (NSE). The approximate hydraulic resistance and flow rate calculated by the HS method were found to be in good agreement with the NSE (less than 7% deviation). Results suggest a quadratic relationship between the network hydraulic resistance and the joint roughness coefficient (JRC). Notably, an increase in surface roughness caused a growth in hydraulic resistance and a fall in flow rate. Further insight was provided by drawing an analogy between resistors in electrical circuits and fractures in networks.

The Okubo-Weiss criterion in hydrodynamic flows:Geometric aspects and further extension

The Okubo [5]-Weiss [6] criterion has been extensively used as a diagnostic tool to divide a two-dimensional (2D) hydrodynamical flow field into hyperbolic and elliptic regions and to serve as a useful qualitative guide to the complex quantitative criteria. The Okubo-Weiss criterion is frequently validated on empirical grounds by the results ensuing its application. So, we will explore topological implications into the Okubo-Weiss criterion and show the Okubo-Weiss parameter is, to within a positive multiplicative factor, the negative of the Gaussian curvature of the underlying vorticity manifold. The Okubo-Weiss criterion is reformulated in polar coordinates, and is validated via several examples including the Lamb- Oseen vortex, and the Burgers vortex. These developments are then extended to 2D quasi- geostrophic (QG) flows. The Okubo-Weiss parameter is shown to remain robust under the -plane approximation to the Coriolis parameter. The Okubo-Weiss criterion is shown to be able to separate the 2D flow-field into coherent elliptic structures and hyperbolic flow configurations very well via numerical simulations of quasi-stationary vortices in QG flows. An Okubo-Weiss type criterion is formulated for 3D axisymmetric flows, and is validated via application to the round Landau-Squire Laminar jet flow.

On the derivation of the Landau-Lifshitz frame in relativistic kinetic theory

The Landau-Lifshitz frame has been widely used to represent the macroscopic quantities of relativistic hydrodynamics in the presence of the dissipative process. In this paper, we derive the Landau-Lifshitz frame in the near-equilibrium regime under self-contained assumptions that do not require comparison with the Eckart frame. And then we revisit the relativistic BGK model proposed by Anderson and Witting to provide an application example of the Landau-Lifshitz frame.

Energy-conserving model of hexagonal pattern in Rayleigh-Bénard convection

We have constructed an energy-conserving sixteen mode dynamical system to model hexagonal pattern in Rayleigh-Bénard convection of Boussinesq fluids with symmetric stress-free thermally conducting boundaries. The model shows stable roll pattern at the onset of convection. Hexagon is found to appear in the system via sausage and (or) stationary rhombus patterns. Both up and down hexagons arise periodically or chaotically with roll, sausage and rhombus patterns. Hexagonal patterns exist for all values of the Prandtl number, 1 ≤ Pr ≤ 5 explored here. However the pattern is more prominent for small Pr and k < kc , where k denotes the wave number. The plot of Nusselt number matches with previous theoretical result. In dissipationless limit, the total energy and the unavailable energy are constants though the kinetic energy, the potential energy and the available energy vary with time. The derived model does not diverge for large values of Rayleigh number Ra.

Shape and scaling of the mean-velocity profile in thermally-stratified plane-Couette flows

It has long been known from measurements that buoyant motions cause the mean-velocity profile (MVP) in thermally-stratified, wall-bounded turbulent flows to significantly deviate from its constant-density counterpart. Theoretical analysis has restricted attention to an “intermediate layer” of the MVP, akin to the celebrated “log layer” in the constant-density case. Here, for thermally-stratified plane-Couette flows, we study the shape and scaling of the whole MVP. We elucidate the mechanisms that dictate the shape of the MVP by using the framework of the spectral link (Gioia et al.; 2010), and obtain scaling laws for the whole MVP by generalizing the Monin-Obukhov similarity theory.

Determination of hydraulic resistance of channels using spectral geometry methods

In this paper, we propose a new method for calculation of hydraulic resistance of channels with constant cross-section. The method is based on the obtained estimates for the average energy dissipation rate in a turbulent flow. The first part of the paper is devoted to theoretical justification of the method. The second part is devoted to calculation of hydraulic resistance of various channels using the abovementioned method and comparison of these values with the known results. The proposed method allows for calculation of hydraulic resistance of various channels with sufficiently high accuracy and is based only on the information about the channel geometry.

Algebraic stability modes in rotational shear flow

We investigate the two-dimensional (2D) stability of rotational shear flows in an unbounded domain. The eigenvalue problem is formulated by using a novel algebraic mode decomposition distinct from the normal modes with temporal evolution $\exp(\omega t)$. Based on the work of \citeasnoun{NoldOberlack2013}, we show how these new modes can be constructed from the symmetries of the linearized stability equation. For the azimuthal base flow velocity $V(r)=r^{-1}$ an additional symmetry exists, such that a mode with algebraic temporal evolution $t^s$ is found. $s$ refers to an eigenvalue for the algebraic growth or decay of the kinetic energy of the perturbations. An eigenvalue problem for the viscous and inviscid stability using algebraic modes is formulated on an infinite domain with $r \to \infty$. An asymptotic analysis of the eigenfunctions shows that the flow is linearly stable under 2D perturbations. We find stable modes with the algebraic mode ansatz, which can not be obtained by a normal mode analysis. The stability results are in line with Rayleigh's inflection point theorem.

New single vortex generation method for flame/vortex interaction using flow behind a rotating cylinder

This paper investigates a new experimental method to generate a single two-dimensional translated vortex for flame/vortex interaction studies. A rotating cylinder is immersed in a uniform flow and, its rotating speed is impulsively reduced. This sudden action triggers the generation of a single vortex when both the initial and the final rotation speeds are in the range of a steady-state regime. Flow visualization allows confirming the applicability of this method, while a complementary two-dimensional numerical simulation is conducted to understand the vortex formation process. A vorticity layer is detached from the cylinder, initiating a feeding process and gradual growth of a single leading vortex. The feeding process is saturated at a specific distance from the cylinder and, vortex separation from the vorticity layer is observed. At the final stage of the formation process, the generated vortex is advected away and, a steady-state regime is again established behind the cylinder. The vortex characteristics appear to be related to the normalized reduction in the rotation rate ∆α, defined as the initial and final rotation rates difference normalized by the initial rotation rate. Several combinations of initial and final rotation rates corresponding to different normalized reductions are investigated experimentally and numerically. The results allow understanding the effect of this parameter; a higher normalized reduction generates a stronger, more rapidly growing vortex. However, its trajectory is related to the wake deviation corresponding to the final rotation rate.