Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

A modular planar robotic manipulandum with end-point torque control.

Robotic manipulanda are extensively used in investigation of the motor control of human arm movements. They permit the application of translational forces to the arm based on its state and can be used to probe issues ranging from mechanisms of neural control to biomechanics. However, most current designs are optimized for studying either motor learning or stiffness. Even fewer include end-point torque control which is important for the simulation of objects and the study of tool use. Here we describe a modular, general purpose, two-dimensional planar manipulandum (vBOT) primarily optimized for dynamic learning paradigms. It employs a carbon fibre arm arranged as a parallelogram which is driven by motors via timing pulleys. The design minimizes the intrinsic dynamics of the manipulandum without active compensation. A novel variant of the design (WristBOT) can apply torques at the handle using an add-on cable drive mechanism. In a second variant (StiffBOT) a more rigid arm can be substituted and zero backlash belts can be used, making the StiffBOT more suitable for the study of stiffness. The three variants can be used with custom built display rigs, mounting, and air tables. We investigated the performance of the vBOT and its variants in terms of effective end-point mass, viscosity and stiffness. Finally we present an object manipulation task using the WristBOT. This demonstrates that subjects can perceive the orientation of the principal axis of an object based on haptic feedback arising from its rotational dynamics.

INVISCID BLADE-TO-BLADE PREDICTION OF A WAKE-GENERATED UNSTEADY FLOW.

A description is presented of a time-marking calculation of the unsteady flow generated by the interaction of upstream wakes with a moving blade row. The inviscid equations of motion are solved using a finite volume technique. Wake dissipation is modeled using an artificial viscosity. Predictions are presented for the rotor mid-span section of an axial turbine. Reasonable agreement is found between the predicted and measured unsteady blade surface static pressures and velocities. These and other results confirm that simple theories can be used to explain the phenomena of rotor-stator wake interactions.

MEASURED AND PREDICTED LOSS GENERATION IN TRANSONIC TURBINE BLADING.

For a typical transonic turbine rotor blade, designed for use with coolant ejection, the trailing edge, or base loss is three to four times the profile boundary layer loss. The base region of such a profile is dominated by viscous effects and it seems essential to attack the problem of loss prediction by solving the compressible Navier-Stokes equations. However, such an approach is inevitably compromised by both numerical accuracy and turbulence modelling constraints. This paper describes a Navier-Stokes solver written for 2D blade-blade flows and employing a simple two-layer mixing length eddy viscosity model. Then, measured and predicted losses and base pressures are presented for two transonic rotor blades and attempts are made to assess the capabilities of the Navier-Stokes solver and to outline areas for future work.

New chiral E- and Z-isomers of the 1R,4R-2-arylidene-p-menthane-3-ones in induced cholesteric and ferroelectric liquid crystals

A novel pair of the E- and Z-isomeric 1R,4R-2-(4-heptyloxyphenyl)-benzylidene-p-menthan-3-ones has been prepared and the influence of distinctions in their molecular geometry on macroscopic properties of liquid crystal systems with the induced supra-molecular helical structure has been studied. The significantly lower helical twisting power of the chiral Z-isomer in comparison with that of E- one has been confirmed in the case of induced cholesteric systems based on 4-pentyl-4-cyanobiphenyl. The phase behavior and ferroelectric characteristics have been investigated for smectic-C* compositions based on the eutectic mixture of the homological 4-hexyloxyphenyl-4'-hexyloxy- and 4-hexyloxyphenyl-4'-octyloxybenzoates containing the novel isomeric chiral dopants. The spontaneous polarisation of the opposite signs induced by the isomeric chiral components has been revealed for the compositions studied. Distinctions in phase states, absolute values of the spontaneous polarization, smectic tilt angle and rotation viscosity of the systems obtained are discussed.

Ferroelectric liquid crystal mixtures containing chiral ether and ester compounds with the 2-arylidene-p-menthan-3-one skeleton

Some 1R,4R-2-(4-phenylbenzylidene)-p-menthane-3-one derivatives containing the ether or ester linking group between benzene rings of the arylidene fragment have been studied as chiral dopants in ferroelectric liquid crystal systems based on the eutectic mixture (1:1) of two phenylbenzoate derivatives CmH2m+1OC6H4COOC6 H4OCnH2n+1 (n = 6; m = 8, 10). The ferroelectric properties of these compositions (spontaneous polarization, rotation viscosity, smectic tilt angle as well as quantitative characteristics of their concentration dependences) were compared with those for systems including chiral dopants containing no linking group. Ferroelectric parameters of the induced ferroelectric compositions studied have been shown to depend essentially on the presence of the linking group between benzene rings and its nature as well as on the number of the benzene rings in the rigid molecular core of the chiral dopants used. For all ferroelectric liquid crystal systems studied, the influence of the chiral dopants on the thermal stability of N*, SmA and SmC* mesophases has been quantified. The influence of the linking group nature in the dopant molecules on the characteristics of the systems studied is discussed taking into account results of the conformational analysis carried out by the semi-empirical AM1 and PM3 methods.

Analysis of the piston ring/liner oil film development during warm-up for an SI-engine

A one-dimensional ring-pack lubrication model developed at MIT is applied to simulate the oil film behavior during the warm-up period of a Kohler spark ignition engine [1]. This is done by making assumptions for the evolution of the oil temperatures during warm-up and that the oil control ring during downstrokes is fully flooded. The ring-pack lubrication model includes features such as three different lubrication regimes, i.e. pure hydrodynamic lubrication, boundary lubrication and pure asperity contact, non-steady wetting of both inlet and outlet of the piston ring, capability to use all ring face profiles that can be approximated by piece-wise polynomials and, finally, the ability to model the rheology of multi-grade oils. Not surprisingly, the simulations show that by far the most important parameter is the temperature dependence of the oil viscosity. This dependence is subsequently examined further by choosing different oils. The baseline oil is SAE 10W30 and results are compared to those using the SAE 30 and the SAE 10W50 oils.

Physical properties of some novel nematic bimesogens for the flexoelectric effect

In a chiral nematic liquid crystal, the flexoelectric effect consists of a fast and linear coupling with an applied electric field. One difficulty to overcome is the unwinding of the helix that occurs at higher fields due to dielectric coupling. The use of bimesogens, which possess very low molecular dielectric anisotropy can improve flexoelectric characteristics. New bimesogen compounds have recently been synthesised that exhibit switching angles of 45° for applied fields of about 9 V.μm-1. In this paper, results from dielectric, electro-optic and dynamic light scattering measurements are reported for the new bimesogenic mixture. The dielectric anisotropy Δε changes sign with temperature and its values range between -0.2 and 0.3 for the temperature range studied. For Δε weakly positive, no electric field Freedericksz transition could be induced but Williams domains are observed instead. The large decrease in the bend elastic constant to viscosity coefficient ratio is attributed to a large increase in the bend viscosity coefficient. © 2001 OPA (Overseas Publishers Association) N.V. Published by license under the Gordon and Breach Science Publishers imprint, a member of the Taylor & Francis Group.

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number W e = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20 % of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

Stabilizing effect of surrounding gas flow on a plane liquid sheet

The stability of a plane liquid sheet is studied experimentally and theoretically, with an emphasis on the effect of the surrounding gas. Co-blowing with a gas velocity of the same order of magnitude as the liquid velocity is studied, in order to quantify its effect on the stability of the sheet. Experimental results are obtained for a water sheet in air at Reynolds number Rel = 3000 and Weber number We = 300, based on the half-thickness of the sheet at the inlet, water mean velocity at the inlet, the surface tension between water and air and water density and viscosity. The sheet is excited with different frequencies at the inlet and the growth of the waves in the streamwise direction is measured. The growth rate curves of the disturbances for all air flow velocities under study are found to be within 20% of the values obtained from a local spatial stability analysis, where water and air viscosities are taken into account, while previous results from literature assuming inviscid air overpredict the most unstable wavelength with a factor 3 and the growth rate with a factor 2. The effect of the air flow on the stability of the sheet is scrutinized numerically and it is concluded that the predicted disturbance growth scales with (i) the absolute velocity difference between water and air (inviscid effect) and (ii) the square root of the shear from air on the water surface (viscous effect).

Improved models for drop-on-demand ink-jet shortening

A previous 1-D model for the shortening of an unbroken drop-on-demand ink-jet ligament has been extended to the case of an arbitrary attached tail mass, and can also include extensional viscosity (which has ∼ 2% effect) as well as linear elasticity in the fluid. Predictions from the improved model are shown to be very similar to results from 2-D axisymmetric numerical simulations of DoD ink-jet ligaments and also to the results of recent experiments on Newtonian fluids jetted without satellite formation.

Frictional flow of damp granular material in a conical centrifuge

An analysis is given of velocity and pressure-dependent sliding flow of a thin layer of damp granular material in a spinning cone. Integral momentum equations for steady state, axisymmetric flow are derived using a boundary layer approximation. These reduce to two coupled first-order differential equations for the radial and circumferential sliding velocities. The influence of viscosity and friction coefficients and inlet boundary conditions is explored by presentation of a range of numerical results. In the absence of any interfacial shear traction the flow would, with increasing radial and circumferential slip, follow a trajectory from inlet according to conservation of angular momentum and kinetic energy. Increasing viscosity or friction reduces circumferential slip and, in general, increases the residence time of a particle in the cone. The residence time is practically insensitive to the inlet velocity. However, if the cone angle is very close to the friction angle then the residence time is extremely sensitive to the relative magnitude of these angles. © 2011 Authors.

The non-local nature of structure functions

Kolmogorov's two-thirds, ((Δv) 2) ∼ e 2/ 3r 2/ 3, and five-thirds, E ∼ e 2/ 3k -5/ 3, laws are formally equivalent in the limit of vanishing viscosity, v → 0. However, for most Reynolds numbers encountered in laboratory scale experiments, or numerical simulations, it is invariably easier to observe the five-thirds law. By creating artificial fields of isotropic turbulence composed of a random sea of Gaussian eddies whose size and energy distribution can be controlled, we show why this is the case. The energy of eddies of scale, s, is shown to vary as s 2/ 3, in accordance with Kolmogorov's 1941 law, and we vary the range of scales, γ = s max/s min, in any one realisation from γ = 25 to γ = 800. This is equivalent to varying the Reynolds number in an experiment from R λ = 60 to R λ = 600. While there is some evidence of a five-thirds law for g > 50 (R λ > 100), the two-thirds law only starts to become apparent when g approaches 200 (R λ ∼ 240). The reason for this discrepancy is that the second-order structure function is a poor filter, mixing information about energy and enstrophy, and from scales larger and smaller than r. In particular, in the inertial range, ((Δv) 2) takes the form of a mixed power-law, a 1+a 2r 2+a 3r 2/ 3, where a 2r 2 tracks the variation in enstrophy and a 3r 2/ 3 the variation in energy. These findings are shown to be consistent with experimental data where the polution of the r 2/ 3 law by the enstrophy contribution, a 2r 2, is clearly evident. We show that higherorder structure functions (of even order) suffer from a similar deficiency.

Improving model quality through computer controlled saturation

In dynamic centrifuge modelling, fluids with enhanced viscosity are often used to correct for the discrepancy in time scaling. However, great care must be taken when using a high viscosity fluid that hydraulic gradients during saturation do not become high enough to cause excessive model disturbance. This paper introduces the CAM-Sat system which aims to improve the saturation process by continually controlling the fluid flow into the model, limiting it to rates low enough to avoid model disturbance. A new method for measuring the fluid flow rate is then described, and its implementation & improvement to the system is discussed. © 2010 Taylor & Francis Group, London.