Three-dimensional direct numerical simulations of autoignition in turbulent non-premixed flows with simple and complex chemistry

3D Direct Numerical Simulations (DNS) of autoignition in turbulent non-premixed flows between fuel and hotter air have been carried out using both 1-step and complex chemistry consisting of a 22 species n-heptane mechanism to investigate spontaneous ignition timing and location. The simple chemistry results showed that the previous findings from 2D DNS that ignition occurred at the most reactive mixture fraction (ξMR) and at small values of the conditional scalar dissipation rate (N|ξMR) are valid also for 3D turbulent mixing fields. Performing the same simulation many times with different realizations of the initial velocity field resulted in a very narrow statistical distribution of ignition delay time, consistent with a previous conjecture that the first appearance of ignition is correlated with the low-N content of the conditional probability density function of N. The simulations with complex chemistry for conditions outside the Negative Temperature Coefficient (NTC) regime show behaviour similar to the single-step chemistry simulations. However, in the NTC regime, the most reactive mixture fraction is very rich and ignition seems to occur at high values of scalar dissipation. Copyright © 2006 by ASME.

Optimisation of the chirp performance of electroabsorption modulators using a numerical system model

A model to simulate an electroabsorption modulator in a dispersive communications system is described and confirmed experimentally for a 5Gbit/s 100km system. Optimisation of the device shows that transmission of 10Gbit/s over 100km is possible.

Numerical simulation of the growth of ZnO nanostructures in a tube furnace by physical vapour deposition

Zinc oxide is a versatile II-VI naturally n-type semiconductor that exhibits piezoelectric properties. By controlling the growth kinetics during a simple carbothermal reduction process a wide range of 1D nanostructures such as nanowires, nanobelts, and nanotetrapods have been synthesized. The driving force: for the nanostructure growth is the Zn vapour supersaturation and supply rate which, if known, can be used to predict and explain the type of crystal structure that results. A model which attempts to determine the Zn vapour concentration as a function of position in the growth furnace is described. A numerical simulation package, COMSOL, was used to simultaneously model the effects of fluid flow, diffusion and heat transfer in a tube furnace made specifically for ZnO nanostructure growth. Parameters such as the temperature, pressure, and flow rate are used as inputs to the model to show the effect that each one has on the Zn concentration profile. An experimental parametric study of ZnO nanostructure growth was also conducted and compared to the model predictions for the Zn concentration in the tube. © 2008 Materials Research Society.

On a classical problem of acoustic wave scattering by a free vortex: Numerical modelling

The scattering of sound from a point source by a Rankine vortex is investigated numerically by solving the Euler equations with the novel high-resolution CABARET method. For several Mach numbers of the vortex, the time-average amplitudes of the scattered field obtained from the numerical modeling are compared with the theoretical scaling laws' predictions. Copyright © 2009 by Sergey Karabasov.

An experimental and numerical study of an oscillating transonic shock wave in a duct

A combined experimental and numerical study of a transonic shock wave in a parallel walled duct subject to downstream pressure perturbations has been conducted. Experiments and simulations have been carried out with a shock strength of M∞ = 1.4 for pressure perturbation frequencies in the range 16-90 Hz. The dynamics of unsteady shock motion and the interaction structure between the unsteady transonic shock wave and the turbulent tunnel floor boundary layer have been investigated. It is found that the (experimentally measured) dynamics of shock motion are generally well predicted by the computational scheme, especially at relatively low (≈ 40 Hz) frequencies. However, at higher frequencies (≈ 90 Hz), some subtle differences between the shock dynamics measured in experiments and those predicted by Computational Fluid Dynamics (CFD) exist. There is evidence from experiments that variations in shock / boundary layer interaction (SBLI) structure caused by shock motion are responsible for a change in the nature of shock dynamics between low and high frequency. In contrast, numerical results at low and high frequencies do not differ significantly and this suggests that the numerical method is not fully capturing the physics of the unsteady flow. Possible reasons for this are considered and a number of areas where CFD is unable to replicate experimental observations are identified. Significantly, CFD predicts changes in SBLI structure due to shock motion that are much too large and this may explain why none of the subtle effects on shock dynamics seen in experiments occur in CFD. Further work developing numerical methods that demonstrate a more realistic sensitivity of SBLI structure to unsteady shock motion is required. Copyright © 2010 by P.J.K. Bruce.

Numerical analysis of the current and voltage sharing issues for resistive fault current limiter using YBCO coated conductors

YBaCuO-coated conductors offer great potential in terms of performance and cost-saving for superconducting fault current limiter (SFCL). A resistive SFCL based on coated conductors can be made from several tapes connected in parallel or in series. Ideally, the current and voltage are shared uniformly by the tapes when quench occurs. However, due to the non-uniformity of property of the tapes and the relative positions of the tapes, the currents and the voltages of the tapes are different. In this paper, a numerical model is developed to investigate the current and voltage sharing problem for the resistive SFCL. This model is able to simulate the dynamic response of YBCO tapes in normal and quench conditions. Firstly, four tapes with different Jc 's and n values in E-J power law are connected in parallel to carry the fault current. The model demonstrates how the currents are distributed among the four tapes. These four tapes are then connected in series to withstand the line voltage. In this case, the model investigates the voltage sharing between the tapes. Several factors that would affect the process of quenches are discussed including the field dependency of Jc, the magnetic coupling between the tapes and the relative positions of the tapes. © 2010 IEEE.

Numerical back-analysis of energy pile test at Lambeth College, London

Energy Piles present an efficient solution for long-term carbon emission reduction and sustainable construction. However, they have received only partial acceptance by the industry, because of concerns regarding the impact of cyclic thermal changes on the serviceability of energy pile foundations. This paper investigates the applicability of the hybrid load transfer approach to load-settlement analysis of single piles behavior during thermal energy exchange processes. Back-analysis results in terms of the thermal and mechanical response of energy piles show good agreement with field test results from Lambeth College in London. © ASCE 2011.

Mandrel peeling of a flexible laminate with a pressure-sensitive adhesive: Comparison of experiment with numerical analysis

Mandrel peel tests with mandrels or rollers of varying diameters have been carried out using Mylar backing of several thicknesses and a commercial synthetic acrylic adhesive. The results are critically compared with the numerical predictions of the peeling software package ICPeel. In addition, a finite element model of the mandrel peeling process has been completed which gives good agreement with experiment provided appropriate mechanical properties of adherend and adhesive are used which must include the effects of adherent constraint. The influence of the thickness of the backing is also considered and both experiment and analysis confirm that there is a backing thickness at which the peel force for a laminate of this sort will show a maximum. © 2010 Blackwell Publishing Ltd.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.