Corner effect and asymmetry in transonic channel flows

An experimental and numerical investigation into transonic shock/boundary-layer interactions in rectangular ducts has been performed. Experiments have shown that flow development in the corners of transonic shock/boundary-layer interactions in confined channels can have a significant impact on the entire flowfield. As shock strength is increased from M∞ = 1:3 to 1.5, the flowfield becomes very slightly asymmetrical. The interaction of corner flows with one another is thought to be a potential cause of this asymmetry. Thus, factors that govern the size of corner interactions (such as interaction strength) and their proximity to one another (such as tunnel aspect ratio) can affect flow symmetry. The results of the computational study show reasonable agreement with experiments, although simulations with particular turbulence models predict highly asymmetrical solutions for flows that were predominantly symmetrical in experiments. These discrepancies are attributed to the tendency of numerical schemes to overprediction corner-interaction size, and this also accounts for why computational fluid dynamics predicts the onset of asymmetry at lower shock strengths than in experiments. The findings of this study highlight the importance of making informed decisions about imposing artificial constraints on symmetry and boundary conditions for internal transonic flows. Future effort into modeling corner flows accurately is required. Copyright © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Eliminating discharge imbalances in predicting shallow water flows with shocks

The shallow water equations are widely used in modelling environmental flows. Being a hyperbolic system of differential equations, they admit shocks that represent hydraulic jumps and bores. Although the water surface can be solved satisfactorily with the modern shock-capturing schemes, the predicted flow rate often suffers from imbalances where shocks occur, eg the mass conservation is violated by failing to maintain a constant discharge rate at every cross-section in a steady open channel flow. A total-variation-diminishing Lax-Wendroff scheme is developed, and used to demonstrate how to achieve an exact flux balance. The performance of the proposed methods is inspected through some test cases, which include 1- and 2-dimensional, flat and irregular bed scenarios. The proposed methods are shown to preserve the mass exactly, and can be easily extended to other shock-capturing models.

Deposition of small particles from turbulent flows

This paper deals with particle deposition onto solid walls from turbulent flows. The aim of the study is to model particle deposition in industrial flows, such as the one in gas turbines. The numerical study has been carried out with a two fluid approach. The possible contribution to the deposition from Brownian diffusion, turbulent diffusion and shear-induced lift force are considered in the study. Three types of turbulent two-phase flows have been studied: turbulent channel flow, turbulent flow in a bent duct and turbulent flow in a turbine blade cascade. In the turbulent channel flow case, the numerical results from a two-dimensional code show good agreement with numerical and experimental results from other resources. Deposition problem in a bent duct flow is introduced to study the effect of curvature. Finally, the deposition of small particles on a cascade of turbine blades is simulated. The results show that the current two fluid models are capable of predicting particle deposition rates in complex industrial flows.

Contact and channel modelling using part and function libraries in a function-based design approach

To support the development and analysis of engineering designs at the embodiment stage, designers work iteratively with representations of those designs as they consider the function and form of their constituent parts. Detailed descriptions of "what a machine does" usually include flows of forces and active principles within the technical system, and their localization within parts and across the interfaces between them. This means that a representation should assist a designer in considering form and function at the same time and at different levels of abstraction. This paper describes a design modelling approach that enables designers to break down a system architecture into its subsystems and parts, while assigning functions and flows to parts and the interfaces between them. In turn, this may reveal further requirements to fulfil functions in order to complete the design. The approach is implemented in a software tool which provides a uniform, computable language allowing the user to describe functions and flows as they are iteratively discovered, created and embodied. A database of parts allows the user to search for existing design solutions. The approach is illustrated through an example: modelling the complex mechanisms within a humanoid robot. Copyright © 2010 by ASME.

Magnetically-actuated artificial cilia for microfluidic propulsion.

In this paper we quantitatively analyse the performance of magnetically-driven artificial cilia for lab-on-a-chip applications. The artificial cilia are fabricated using thin polymer films with embedded magnetic nano-particles and their deformation is studied under different external magnetic fields and flows. A coupled magneto-mechanical solid-fluid model that accurately captures the interaction between the magnetic field, cilia and fluid is used to simulate the cilia motion. The elastic and magnetic properties of the cilia are obtained by fitting the results of the computational model to the experimental data. The performance of the artificial cilia with a non-uniform cross-section is characterised using the numerical model for two channel configurations that are of practical importance: an open-loop and a closed-loop channel. We predict that the flow and pressure head generated by the artificial cilia can be as high as 18 microlitres per minute and 3 mm of water, respectively. We also study the effect of metachronal waves on the flow generated and show that the fluid propelled increases drastically compared to synchronously beating cilia, and is unidirectional. This increase is significant even when the phase difference between adjacent cilia is small. The obtained results provide guidelines for the optimal design of magnetically-driven artificial cilia for microfluidic propulsion.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Off-wall boundary conditions for turbulent flows obtained from buffer-layer minimal flow units

There is strong evidence that the transport processes in the buffer region of wall-bounded turbulence are common across various flow configurations, even in the embryonic turbulence in transition (Park et al., Phys. Fl. 24). We use this premise to develop off-wall boundary conditions for turbulent simulations. Boundary conditions are constructed from DNS databases using periodic minimal flow units and reduced order modeling. The DNS data was taken from a channel at Reτ=400 and a zero-pressure gradient transitional boundary layer (Sayadi et al., submitted to J. Fluid Mech.). Both types of boundary conditions were first tested on a DNS of the core of the channel flow with the aim of extending their application to LES and to spatially evolving flows.

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 velocity fluctuations in off-wall boundary conditions for turbulent flows

A model for off-wall boundary conditions for turbulent flow is investigated. The objective of such a model is to circumvent the need to resolve the buffer layer near the wall, by providing conditions in the logarithmic layer for the overlying flow. The model is based on the self-similarity of the flow at different heights in the logarithmic layer. It was first proposed by Mizuno and Jiménez (2013), imposing at the boundary plane a velocity field obtained on-the-fly from an overlying region. The key feature of the model was that the lengthscales of the field were rescaled to account for the self-similarity law. The model was successful at sustaining a turbulent logarithmic layer, but resulted in some disagreements in the flow statistics, compared to fully-resolved flows. These disagreements needed to be addressed for the model to be of practical application. In the present paper, a more refined, wavelength-dependent rescaling law is proposed, based on the wavelength-dependent dynamics in fully-resolved flows. Results for channel flow show that the new model eliminates the large artificial pressure fluctuations found in the previous one, and a better agreement is obtained in the bulk properties, the flow fluctuations, and their spectral distribution across the whole domain. © Published under licence by IOP Publishing Ltd.