Viscous-slip,thermal-slip,and temperature-jump coefficients based on the linearized Boltzmann equation (and five kinetic models) with the Cercignani–Lampis boundary condition

A polynomial expansion procedure and the ADO (analytical discrete-ordinates) method are used to compute the viscous-slip coefficient, the thermal-slip coefficient, and the temperature-jump coefficient from the linearized Boltzmann equation (LBE) for rigid-sphere interactions and the Cercignani–Lampis (CL) boundary condition. These same quantities are also computed from five kinetic models, with the CL condition, and compared to the LBE result. Equivalent results for the LBE and the kinetic models, all based on the usual Maxwell boundary condition, are also reported.     

Rarefied gas flow in a channel with specular‐diffuse boundary conditions for arbitrary knudsen numbers. Onsager relations

The problem of a rarefied gas flow in a channel for arbitrary Knudsen numbers has been solved analytically for the first time in the case where the scattering of gas molecules on the channel walls can be described by speculardiffuse boundary conditions. The mean free path of gas molecules is assumed to be constant, i.e., the collision frequency is proportional to molecular velocity. The gas moves under the action of a streamwise temperature gradient. Exact relations for heat and mass fluxes and for meanmass velocity are obtained. It is shown that the Onsager relations are valid within the entire range of Knudsen numbers in the problem of heat and mass transfer in a channel. The dependence of heat and mass fluxes on the Knudsen number (channel thickness) is analyzed. A comparison with available results is performed.     

Study of laminar–turbulent flow transition under pulsatile conditions in a constricted channel

In this paper, direct numerical simulation is performed to investigate a pulsatile flow in a constricted channel to gain physical insights into laminar–turbulent–laminar flow transitions. An in-house computer code is used to conduct numerical simulations based on available high-performance shared memory parallel computing facilities. The Womersley number tested is fixed to 10.5 and the Reynolds number varies from 500 to 2000. The influences of the degree of stenosis and pulsatile conditions on flow transitions and structures are investigated. In the region upstream of the stenosis, the flow pattern is primarily laminar. Immediately after the stenosis, the flow recirculates under an adverse streamwise pressure gradient, and the flow pattern transitions from laminar to turbulent. In the region far downstream of the stenosis, the flow becomes re-laminarised. The physical characteristics of the flow field have been thoroughly analysed in terms of the mean streamwise velocity, turbulence kinetic energy, viscous wall shear stresses, wall pressure and turbulence kinetic energy spectra.     

Residual mass and flow regimes for the immiscible liquid–liquid displacement in a plane channel

The motion of two immiscible liquids in a plane channel is analyzed for the case in which the flow conditions and the interactions between the liquids and the solid surface maintain the displaced fluid attached to the wall. The Galerkin Finite Element Method is used to compute the velocity field and the configuration of the interface between the two fluids. We compare the residual mass fraction left on the wall with its two counterparts in capillary tubes, namely residual mass fraction and dimensionless layer thickness of the displaced fluid. The main result of this comparison was that although there is a qualitative similarity concerning the layer thickness between the two cases, the residual fraction of mass presented an important difference, showing that when the aspect ratio of the capillary passage is large there is an increase in the displacement efficiency. The thickness of the displaced liquid film attached to the channel walls is a function of the capillary number (Ca) and the viscosity ratio (N_μ). A map of streamlines in the Cartesian space (Ca, N_μ) with the different flow regimes of the problem is presented. We also showed that we can adapt the available analytical results obtained for gas-displacement in capillary tubes to the plane channel case, for low values of Ca.     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

Comment on “A theorem for a shear stress-free plane boundary in Stokes flow”

The results proved in a recent paper by Akhtar and Sen [Mech. Res. Commun. 38 (2011) 529–531] are not new and are already known in the literature. In fact, the main expressions – presented in the form of a theorem – in the cited article follow as a special case of the more general results for a two-fluid Stokes flow with a planar interface proved earlier by Palaniappan [Acta Mechanica 139 (2000) 1–13]. We briefly demonstrate the reduction of particular cases of the solution given by Palaniappan which also reveals the duplication of another special solution elsewhere.     

Application of solution structure theorems to Cattaneo–Vernotte heat conduction equation with non-homogeneous boundary conditions

In this study, a non-Fourier heat conduction problem formulated using the Cattaneo–Vernotte (C–V) model with non-homogeneous boundary conditions is solved with the superposition principle in conjunction with solution structure theorems. It is well known that the aforementioned analytical method is not suitable for such a class of thermal problems. However, by performing a functional transformation, the original non-homogeneous partial differential equation governing the physical problem can be cast into a new form such that it consists of a homogeneous part and an additional auxiliary function. As a result, the modified homogeneous governing equation can then be solved with solution structure theorems for temperatures inside a finite planar medium. The methodology provides a convenient, accurate, and efficient solution to the C–V heat conduction equation with non-homogeneous boundary conditions.     

Direct simulation of liquid–gas–solid flow with a free surface lattice Boltzmann method

Direct numerical simulation of liquid–gas–solid flows is uncommon due to the considerable computational cost. As the grid spacing is determined by the smallest involved length scale, large grid sizes become necessary – in particular, if the bubble–particle aspect ratio is on the order of 10 or larger. Hence, it arises the question of both feasibility and reasonability. In this paper, we present a fully parallel, scalable method for direct numerical simulation of bubble–particle interaction at a size ratio of 1–2 orders of magnitude that makes simulations feasible on currently available super-computing resources. With the presented approach, simulations of bubbles in suspension columns consisting of more than 100,000 fully resolved particles become possible. Furthermore, we demonstrate the significance of particle-resolved simulations by comparison to previous unresolved solutions. The results indicate that fully resolved direct numerical simulation is indeed necessary to predict the flow structure of bubble–particle interaction problems correctly.     

Nusselt number and friction factor in thermally stratified turbulent channel flow under Non-Oberbeck–Boussinesq conditions

In stably stratified turbulence, computations under Oberbeck–Boussinesq (OB) hypothesis of temperature-independent fluid properties may lead to inaccurate representation of the flow field and to wrong estimates of momentum/heat transfer coefficients. This is clearly assessed here comparing direct numerical simulations of stratified turbulence under OB conditions to simulations under NOB (Non-Oberbeck–Boussinesq) conditions of temperature-dependent fluid viscosity and thermal expansion coefficient. Compared to the OB case, NOB conditions may induce local flow relaminarization with significant variations (up to 30%) of heat and momentum transfer coefficients. Together with DNS results, we propose a phenomenological model (based on turbulent bursts) for heat transfer prediction in stratified turbulence under OB and NOB conditions. Implications of NOB assumptions on mixing efficiency (i.e. flux Richardson number Ri_f) and turbulent Prandtl number (Pr_t) are also discussed. These results are of specific importance in RANS modelling, where the condition Pr_t = 1 is usually assumed (Reynolds analogy). Although this assumption is valid in some situations (i.e. boundary layer, pipe flow) there is uncertainty about its validity for stably-stratified turbulence. We demonstrate that this assumption is inaccurate when NOB effects become significant.     

The PELskin project—part I: fluid–structure interaction for a row of flexible flaps: a reference study in oscillating channel flow

Previous studies of flexible flaps attached to the aft part of a cylinder have demonstrated a favourable effect on the drag and lift force fluctuation. This observation is thought to be linked to the excitation of travelling waves along the flaps and as a consequence of that, periodic shedding of the von Kármán vortices is altered in phase. A more general case of such interaction is studied herein for a limited row of flaps in an oscillating flow; representative of the cylinder case since the transversal flow in the wake-region shows oscillating character. This reference case is chosen to qualify recently developed numerical methods for the simulation of fluid–structure interaction in the context of the EU funded ‘PELskin’ project. The simulation of the two-way coupled dynamics of the flexible elements is achieved via a structure model for the flap motion, which was implemented and coupled to two different fluid solvers via the immersed boundary method. The results show the waving behaviour observed at the tips of the flexible elements in interaction with the fluid flow and the formation of vortices in the gaps between the flaps. In addition, formation of vortices upstream of the leading and downstream of the trailing flap is seen, which interact with the formation of the shear-layer on top of the row. This leads to a phase shift in the wave-type motion along the row that resembles the observation in the cylinder case.     

Thermo-fluid-dynamic analysis of the flow in a rotating channel with a sharp “U” turn

Infrared thermography has been employed to carry out a detailed convective heat transfer measurements at Re?=?20,000 in a two-pass square channel both for the static case (absence of channel rotation) and for the rotating case (Ro?=?0.3). At the same time, the main and secondary flow fields have been measured by means of particle image velocimetry with the aim to investigate how the flow behavior affects the local distributions of the convective heat transfer coefficient for the two cases. The normal-to-wall velocity component (w) and the turbulent kinetic energy, both measured close to the heat exchanging wall, have been used to formulate an empirical heat transfer correlation within an attempt to identify the role performed by these two quantities on the convective heat transfer coefficient distributions. The latter ones have been reported in terms of normalized Nusselt number (Nu/Nu*) maps, where Nu* is the Nusselt number evaluated with the classical Dittus-B?lter correlation.     

Kinematics of Bulkley–Herschel fluid flow with a free surface during the filling of a channel

The non-Newtonian fluid flow with a free surface occurring during the filling of a plane channel in the gravity field is modeled. The mathematical formulation of the problem using the rheological Bulkley–Herschel model is presented. A numerical finite-difference algorithm for solving this problem is developed. A parametric investigation of the main characteristics of the process as functions of the control parameters is performed. The effect of the rheological parameters of the fluid on the distribution of the quasisolid motion zones is demonstrated.     

Four-way coupled Eulerian–Lagrangian Direct Numerical Simulations in a vertical laminar channel flow

Direct Numerical Simulations of a laminar two-phase flow into a vertical channel are investigated. An Eulerian–Lagrangian approach allows tracking each bubble position with a four-way coupling strategy, i.e. taking into account bubble-fluid and bubble-bubble interactions. The flow configuration has been chosen to highlight the buoyancy effects due to significant values of void fraction (high numbers of bubbles); hence the bubbles collisions and wall effects are the critical parameters to ensure the dispersion of the bubble plume. The DNS approach is self-consistent and does not rely on closure relations or empirical correlations for describing the collective bubble dynamics. It is found that the DNS predicts well the behavior of the bubble plume and its back effect on the liquid phase when compared with a mixture model and experimental data. The elastic nature of collisions, the sensitivity of the mean and RMS values of velocities and void fraction to the mesh quality are explored.     

Thermal boundary layer flow over a stretching sheet in a micropolar fluid with radiation effect

In the present paper, we study the effects of radiation on the thermal boundary layer flow induced by a linearly stretching sheet immersed in an incompressible micropolar fluid with constant surface temperature. Similarity transformation is employed to transform the governing partial differential equations into ordinary ones, which are then solved numerically using the Runge-Kutta-Fehlberg method. Results for the local Nusselt number as well as the temperature profiles are presented for different values of the governing parameters. It is found that the heat transfer rate at the surface decreases in the presence of radiation. Comparison with known results for certain particular cases is excellent.     

General soliton solutions to a nonlocal long-wave–short-wave resonance interaction equation with nonzero boundary condition

Nonlinear Dynamics - Under investigation in this work is a newly proposed nonlocal long-wave–short-wave resonance interaction (LSRI) equation with the self-induced parity-time (PT) symmetric...     

The flow of an Oldroyd-B fluid past a cylinder in a channel: adaptive viscosity vorticity (DAVSS-ω) formulation

A parallel unstructured finite volume method (FVM) is developed and implemented under a distributed computing environment through the parallel virtual machine (PVM) libraries, and is used to simulate the channel flow of the Oldroyd-B fluid past a circular cylinder. Differing from our previous work 11, 12, a discrete elastic viscous split stress (DEVSS) formulation together with an independent interpolation of the vorticity (DEVSS-ω) is proposed in this paper. This method has almost the same stability behavior as the elastic viscous split stress (EVSS) formulation, and is suitable for complex constitutive models. To further improve the stability at high Deborah numbers, we combine the idea of the discrete adaptive elastic viscous split stress (DAVSS) formulation [7] with the independent interpolation of the vorticity to arrive at the DAVSS-ω method. The numerical implementation is based on the unstructured FVM method and the semi-implicit method for pressure-linked equations revised (SIMPLER) algorithm. The parallelization of the program is implemented by a domain decomposition strategy and using PVM software libraries. The results are compared with those by the EVSS, DEVSS, and the plain Oldroyd-B formulation (without splitting the stress). It is found that the drag coefficient first decreases and then increases with the De number, for a channel half width to cylinder radius ratio of h/R = 2. It is also confirmed that the drag enhancement at high Deborah number is due to the increasing extension effect in the regions near the front and the rear stagnation points.     

The flow pattern map of a two-phase non-Newtonian liquid–gas flow in the vertical pipe

A map for the determination of flow pattern for two-phase flow of gas and non-Newtonian liquid in the vertical pipe has been presented. Our own experimental data confirm applicability of such a map.