Long range exhaustion-a mathematical model for the axisymmetric air flow of a local exhaust ventilation hood assisted by a turbulent radial jet

Abstract-Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art Banged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air Bow to be reduced and provide a means of identifying and assessing the various parameters that control the air Bow. The mathematical model is formulated in terms of the Stokes steam function, ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained Bow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions. | Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art flanged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air flow to be reduced and provide a means of identifying and assessing the various parameters that control the air flow. The mathematical model is formulated in terms of the Stokes steam function, Ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained flow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates. © Cambridge University Press 2013.

Structural sensitivity of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this. We investigate this spiral mode with a linear global stability analysis around the steady bubble and its wake. We obtain the linear direct and adjoint global modes of the linearized Navier-Stokes equations and overlap these to obtain the structural sensitivity of the spiral mode, which identifies the wavemaker region. We also identify regions of absolute instability with a local stability analysis. At moderate swirls, we find that the m=-1 azimuthal mode is the most unstable and that the wavemaker regions of the m=-1 mode lie around the bubble, which is absolutely unstable. The mode is most sensitive to feedback involving the radial and azimuthal components of momentum in the region just upstream of the bubble. To a lesser extent, the mode is also sensitive to feedback involving the axial component of momentum in regions of high shear around the bubble. At an intermediate swirl, in which the bubble and wake have similar absolute growth rates, other researchers have found that the wavemaker of the nonlinear global mode lies in the wake. We agree with their analysis but find that the regions around the bubble are more influential than the wake in determining the growth rate and frequency of the linear global mode. The results from this paper provide the first steps towards passive control strategies for spiral vortex breakdown. © 2013 Cambridge University Press.

The propagation of acoustic waves in a slowly varying duct with radially sheared axial mean flow

We consider the propagation of acoustic waves along a cylindrical duct carrying radially sheared axial mean flow, in which the duct radius is allowed to vary slowly along the axis. In previous work [A.J. Cooper & N. Peake, Journal of Fluid Mechanics 445 (2001) 207-234.] radially sheared axial mean flow with nonzero swirl in a slowly varying duct was considered, but in this paper we set the swirl to zero, thereby allowing simplification of the calculations of both the mean and unsteady flows. In this approach the acoustic wavenumber and corresponding eigenfunction are determined locally, while the wave amplitude is found by solving an evolution equation along the duct. Sample results are presented, including a case in which, perhaps surprisingly, the number of cut-on modes increases as the duct radius decreases. © 2013 Elsevier Ltd. All rights reserved.

Guidance of collective cell migration by substrate geometry.

Collective behavior refers to the emergence of complex migration patterns over scales larger than those of the individual elements constituting a system. It plays a pivotal role in biological systems in regulating various processes such as gastrulation, morphogenesis and tissue organization. Here, by combining experimental approaches and numerical modeling, we explore the role of cell density ('crowding'), strength of intercellular adhesion ('cohesion') and boundary conditions imposed by extracellular matrix (ECM) proteins ('constraints') in regulating the emergence of collective behavior within epithelial cell sheets. Our results show that the geometrical confinement of cells into well-defined circles induces a persistent, coordinated and synchronized rotation of cells that depends on cell density. The speed of such rotating large-scale movements slows down as the density increases. Furthermore, such collective rotation behavior depends on the size of the micropatterned circles: we observe a rotating motion of the overall cell population in the same direction for sizes of up to 200 μm. The rotating cells move as a solid body, with a uniform angular velocity. Interestingly, this upper limit leads to length scales that are similar to the natural correlation length observed for unconfined epithelial cell sheets. This behavior is strongly altered in cells that present a downregulation of adherens junctions and in cancerous cell types. We anticipate that our system provides a simple and easy approach to investigate collective cell behavior in a well-controlled and systematic manner.

The rhythm of fountains: The length and time scales of rise height fluctuations at low and high Froude numbers

The magnitude and frequency of vertical fluctuations of the top of an axisymmetric miscible Boussinesq fountain forms the focus of this work. We present measurements of these quantities for saline-aqueous fountains in uniform quiescent surroundings. Our results span source Froude numbers 0.3 ≤ Fr 0 ≤ 40 and, thereby, encompass very weak, weak, intermediate and forced classes of fountain. We identify distinct scalings, based on known quantities at the fountain source, for the frequency of fountain height fluctuations which collapse our data within bands of Fr0. Notably, our scalings reveal that the (dimensionless) frequency takes a constant value within each band. These results highlight characteristic time scales for the fluctuations which we decompose into a single, physically apparent, length scale and velocity scale within each band. Moreover, within one particular band, spanning source Froude numbers towards the lower end of the full range considered, we identify unexpectedly long-period fluctuations indicating a near balance of inertia and (opposing) buoyancy at the source. Our analysis identifies four distinct classes of fluctuation behaviour (four bands of Fr 0) and this classification matches well with existing classifications of fountains based on rise heights. As such, we show that an analysis of the behaviour of the fountain top alone, rather than the entire fountain, provides an alternative approach to classifying fountains. The similarity of classifications based on the two different methods confirms that the boundaries between classes mark tangible changes in the physics of fountains. For high Fr0 we show that the dominant fluctuations occur at the scale of the largest eddies which can be contained within the fountain near its top. Extending this, we develop a Strouhal number, Strtop, based on experimental measures of the fountain top, defined such that Strtop = 1 would suggest the dominant fluctuations are caused by a continual cycle of eddies forming and collapsing at this largest physical scale. For high- Fr 0 fountains we find Strtop ≈ 0. 9. © 2013 Cambridge University Press.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Scaling of transitionally-rough flow structures in DNSes of channels with riblets up to Reτ=550

The linear, drag-reducing effect of vanishingly small riblets breaks down once their size is in the transitionally-rough regime. We have previously reported that this breakdown is caused by the additional Reynolds stresses produced by the appearance of elongated spanwise rollers just above the riblet surface. These rollers are related with the Kelvin--Helmholtz instability of free shear layers, and to similar structures appearing over other rough and porous surfaces. However, because of the limited Reτ=180 in our previous DNSes, it could not be determined whether those structures scaled in inner or outer units. Furthermore, it is questionable if results in the transitionally-rough regime at Reτ=180 can be extrapolated to configurations of practical interest. At such small Reynolds numbers, roughness of transitional size can perturb a large portion of the boundary layer, which is not the case in most industrial and atmospheric applications. To clarify these issues we have conducted a set of DNSes at Reτ=550. Our results indicate that the spanwise rollers scale in wall units, and support the validity of the extrapolation to configurations of practical interest.

Kelvin-Helmholtz-like instability of turbulent flows over riblets

The turbulent drag reduction due to riblets is a function of their size and, for different configurations, collapses well with a length scale l+g=(A+g)1/2, based in the groove cross-section Ag. The initially linear drag reduction breaks down for l+g≈11, which agrees in our DNS with the previously reported appearance of quasi-two-dimensional spanwise rollers immediately above the riblets. They are similar to those found over porous surfaces and plant canopies, and can be traced to a Kelvin-Helmholtz-like instability associated with the relaxation of the impermeability condition for the wall-normal velocity. The extra Reynolds stress associated with them accounts quantitatively for the drag degradation. An inviscid model for the instability confirms its nature, agreeing well with the observed perturbation wavelengths and shapes. The onset of the instability is determined by a length scale L+w that, for conventional riblet geometries, is proportional to l+g. The instability onset, L+w≥4, corresponds to the empirical breakdown point l+g≈11.