RANS-URANS in axial compressor, a design methodology

Turbomachinery flows are inherently unsteady. Until now during the design process, unsteadiness has been neglected, with resort merely to steady numerical simulations. Despite the assumption involved, the results obtained with steady simulations have been used with success. One of the questions arising in recent years is can unsteady simulations be used to improve the design of turbomachines? In this work the numerical simulation of a multi-stage axial compressor is carried out. Comparison of Reynolds averaged Navier-Stokes (RANS) and unsteady Reynolds averaged Navier-Stokes (URANS) calculation shows that the unsteadiness affects pressure losses and the prediction of stall limit. The unsteady inflow due to the wake passing mainly modifies the losses and whirl angle near the endwalls. The computational cost of the fully unsteady compared with a steady simulation is about four times in terms of mesh dimension and two orders of magnitude as number of iterations. A mixed RANS-URANS solution has been proposed to give the designer the possibility to simulate an unsteady stage embedded in a steady-state simulation. This method has been applied to the simulation of a four-stage axial compressor rig. The mixed RANS-URANS approach has been developed using sliding and mixing planes as interface conditions. The rotor-stator interaction has been captured physically while reducing the computational time and mesh size.

Global linear and nonlinear stability of viscous confined plane wakes with co-flow

The global stability of confined uniform density wakes is studied numerically, using two-dimensional linear global modes and nonlinear direct numerical simulations. The wake inflow velocity is varied between different amounts of co-flow (base bleed). In accordance with previous studies, we find that the frequencies of both the most unstable linear and the saturated nonlinear global mode increase with confinement. For wake Reynolds number Re = 100 we find the confinement to be stabilising, decreasing the growth rate of the linear and the saturation amplitude of the nonlinear modes. The dampening effect is connected to the streamwise development of the base flow, and decreases for more parallel flows at higher Re. The linear analysis reveals that the critical wake velocities are almost identical for unconfined and confined wakes at Re ≈ 400. Further, the results are compared with literature data for an inviscid parallel wake. The confined wake is found to be more stable than its inviscid counterpart, whereas the unconfined wake is more unstable than the inviscid wake. The main reason for both is the base flow development. A detailed comparison of the linear and nonlinear results reveals that the most unstable linear global mode gives in all cases an excellent prediction of the initial nonlinear behaviour and therefore the stability boundary. However, the nonlinear saturated state is different, mainly for higher Re. For Re = 100, the saturated frequency differs less than 5% from the linear frequency, and trends regarding confinement observed in the linear analysis are confirmed.

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.

Characterization of dynamic nonlinear absorption of carbon nanotube saturable absorber

Dynamic nonlinear absorption of composite-type single-wall carbon nanotube saturable absorbers is characterized using both femtosecond and picosecond pump pulses. Results are compared with numerical simulations based on two commonly used saturable absorber models. © 2010 Optical Society of America.

Experimental demonstration of a record high 11.25Gb/s real-time optical OFDM transceiver supporting 25km SMF end-to-end transmission in simple IMDD systems.

The fastest ever 11.25Gb/s real-time FPGA-based optical orthogonal frequency division multiplexing (OOFDM) transceivers utilizing 64-QAM encoding/decoding and significantly improved variable power loading are experimentally demonstrated, for the first time, incorporating advanced functionalities of on-line performance monitoring, live system parameter optimization and channel estimation. Real-time end-to-end transmission of an 11.25Gb/s 64-QAM-encoded OOFDM signal with a high electrical spectral efficiency of 5.625bit/s/Hz over 25km of standard and MetroCor single-mode fibres is successfully achieved with respective power penalties of 0.3dB and -0.2dB at a BER of 1.0 x 10(-3) in a directly modulated DFB laser-based intensity modulation and direct detection system without in-line optical amplification and chromatic dispersion compensation. The impacts of variable power loading as well as electrical and optical components on the transmission performance of the demonstrated transceivers are experimentally explored in detail. In addition, numerical simulations also show that variable power loading is an extremely effective means of escalating system performance to its maximum potential.

ZnO-Based FBAR resonators with carbon nanotube electrodes

Film bulk acoustic resonator (FBAR) devices with carbon nanotube (CNT) electrodes directly grown on a ZnO film by thermal chemical vapor deposition have been fabricated. CNT electrodes possess a very low density and high acoustic impedance, which reduces the intrinsic mass loading effect resulting from the electrodes' weight and better confines the longitudinal acoustic standing waves inside the resonator, in turn providing a resonator with a higher quality factor. The influence of the CNTs on the frequency response of the FBAR devices was studied by comparing two identical sets of devices; one set comprised FBARs fabricated with chromium/ gold bilayer electrodes, and the second set comprised FBARs fabricated with CNT electrodes. It was found that the CNTs had a significant effect on attenuating traveling waves at the surface of the FBARs' membranes because of their high elastic stiffness. Three-dimensional finite element analysis of the devices fabricated was carried out, and the numerical simulations were consistent with the experimental results obtained. © 2011 IEEE.

The internal structure of igniting turbulent sprays as revealed by complex chemistry DNS

A parametric study of spark ignition in a uniform monodisperse turbulent spray is performed with complex chemistry three-dimensional Direct Numerical Simulations in order to improve the understanding of the structure of the ignition kernel. The heat produced by the kernel increases with the amount of fuel evaporated inside the spark volume. Moreover, the heat sink by evaporation is initially higher than the heat release and can have a negative effect on ignition. With the sprays investigated, heat release occurs over a large range of mixture fractions, being high within the nominal flammability limits and finite but low below the lean flammability limit. The burning of very lean regions is attributed to the diffusion of heat and species from regions of high heat release, and from the spark, to lean regions. Two modes of spray ignition are reported. With a relatively dilute spray, nominally flammable material exists only near the droplets. Reaction zones are created locally near the droplets and have a non-premixed character. They spread from droplet to droplet through a very lean interdroplet spacing. With a dense spray, the hot spark region is rich due to substantial evaporation but the cold region remains lean. In between, a large surface of flammable material is generated by evaporation. Ignition occurs there and a large reaction zone propagates from the rich burned region to the cold lean region. This flame is wrinkled due to the stratified mixture fraction field and evaporative cooling. In the dilute spray, the reaction front curvature pdf contains high values associated with single droplet combustion, while in the dense spray, the curvature is lower and closer to the curvature associated with gaseous fuel ignition kernels. © 2011 The Combustion Institute.

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.

Breakdown voltage for superjunction power devices with charge imbalance: An analytical model valid for both punch through and non punch through devices

An analytical model for the electric field and the breakdown voltage (BV) of an unbalanced superjunction (SJ) device is presented in this paper. The analytical technique uses a superposition approach treating the asymmetric charge in the pillars as an excess charge component superimposed on a balanced charge component. The proposed double-exponentialmodel is able to accurately predict the electric field and the BV for unbalanced SJ devices in both punch through and non punch through conditions. The model is also reasonably accurate at extremely high levels of charge imbalance when the devices behave similarly to a PiN diode or to a high-conductance layer. The analytical model is compared against numerical simulations of charge unbalanced SJ devices and against experimental results. © 2009 IEEE.

Jetting behaviour of polymer solutions in drop-on-demand inkjet printing

The jetting of dilute polymer solutions in drop-on-demand printing is investigated. A quantitative model is presented which predicts three different regimes of behaviour depending upon the jet Weissenberg number Wi and extensibility of the polymer molecules. In regime I (Wi < ½) the polymer chains are relaxed and the fluid behaves in a Newtonian manner. In regime II (½ < Wi < L) where L is the extensibility of the polymer chain the fluid is viscoelastic, but the polymer do not reach their extensibility limit. In regime III (Wi > L) the chains remain fully extended in the thinning ligament. The maximum polymer concentration at which a jet of a certain speed can be formed scales with molecular weight to the power of (1-3ν), (1-6ν) and -2ν in the three regimes respectively, where ν is the solvent quality coefficient. Experimental data obtained with solutions of mono-disperse polystyrene in diethyl phthalate with molecular weights between 24 - 488 kDa, previous numerical simulations of this system, and previously published data for this and another linear polymer in a variety of “good” solvents, all show good agreement with the scaling predictions of the model.

The analysis of impact forces in randomly vibrating elastic systems

This work is concerned with the characteristics of the impact force produced when two randomly vibrating elastic bodies collide with each other, or when a single randomly vibrating elastic body collides with a stop. The impact condition includes a non-linear spring, which may represent, for example, a Hertzian contact, and in the case of a single body, closed form approximate expressions are derived for the duration and magnitude of the impact force and for the maximum deceleration at the impact point. For the case of two impacting bodies, a set of algebraic equations are derived which can be solved numerically to yield the quantities of interest. The approach is applied to a beam impacting a stop, a plate impacting a stop, and to two impacting beams, and in each case a comparison is made with detailed numerical simulations. Aspects of the statistics of impact velocity are also considered, including the probability that the impact velocity will exceed a specified value within a certain time. © 2012 Elsevier Ltd. All rights reserved.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Global nonlinear dynamics of thin aerofoil wakes

In the present investigation of thin aerofoil wakes we compare the global nonlinear dynamics, obtained by direct numerical simulations, to the associated local instability features, derived from linear stability analyses. A given configuration depends on two control parameters: the Reynolds number Re and the adverse pressure gradient m (with m < 0) prevailing at the aerofoil trailing edge. Global instability is found to occur for large enough Re and |m|; the naturally selected frequency is determined by the local absolute frequency prevailing at the trailing edge. © 2010 Springer Science+Business Media B.V.

Compact three-dimensional silicon termination solutions for high voltage SOI SuperJunction

This paper demonstrates and discusses novel "three dimensional" silicon based junction isolation/termination solutions suitable for high density ultra-low-resistance Lateral Super-Junction structures. The proposed designs are both compact and effective in safely distributing the electrostatic potential away from the active device area. The designs are based on the utilization of existing layers in the device fabrication line, hence resulting in no extra complexity or cost increase. The study/demonstration is done through extensive experimental measurements and numerical simulations. © 2012 IEEE.

A coordination-based approach to elasticity of floppy and stiff random networks

We study the role of connectivity on the linear and nonlinear elastic behavior of amorphous systems using a two-dimensional random network of harmonic springs as a model system. A natural characterization of these systems arises in terms of the network coordination relative to that of an isostatic network $\delta z$; a floppy network has $\delta z<0$, while a stiff network has $\delta z>0$. Under the influence of an externally applied load we observe that the response of both floppy and rigid network are controlled by the same critical point, corresponding to the onset of rigidity. We use numerical simulations to compute the exponents which characterize the shear modulus, the amplitude of non-affine displacements, and the network stiffening as a function of $\delta z$, derive these theoretically and make predictions for the mechanical response of glasses and fibrous networks.