Prediction of tread block forces for a free-rolling tyre in contact with a rough road

This paper describes a new approach to model the forces on a tread block for a free-rolling tyre in contact with a rough road. A theoretical analysis based on realistic tread mechanical properties and road roughness is presented, indicating partial contact between a tread block and a rough road. Hence an asperity-scale indentation model is developed using a semi-empirical formulation, taking into account both the rubber viscoelasticity and the tread block geometry. The model aims to capture the essential details of the contact at the simplest level, to make it suitable as part of a time-domain dynamic analysis of the coupled tyre-road system. The indentation model is found to have a good correlation with the finite element (FE) predictions and is validated against experimental results using a rolling contact rig. When coupled to a deformed tyre belt profile, the indentation model predicts normal and tangential force histories inside the tyre contact patch that show good agreement with FE predictions. © 2012 Elsevier B.V..

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

Modelling and Fragility Analysis of Non-Ductile Reinforced Concrete Buildings in Low-to-Moderate Seismic Zones

Reinforced concrete buildings in low-to-moderate seismic zones are often designed only for gravity loads in accordance with the non-seismic detailing provisions. Deficient detailing of columns and beam-column joints can lead to unpredictable brittle failures even under moderate earthquakes. Therefore, a reliable estimate of structural response is required for the seismic evaluation of these structures. For this purpose, analytical models for both interior and exterior slab-beam-column subassemblages and for a 1/3 scale model frame were implemented into the nonlinear finite element platform OpenSees. Comparison between the analytical results and experimental data available in the literature is carried out using nonlinear pushover analyses and nonlinear time history analysis for the subassemblages and the model frame, respectively. Furthermore, the seismic fragility assessment of reinforced concrete buildings is performed on a set of non-ductile frames using nonlinear time history analyses. The fragility curves, which are developed for various damage states for the maximum interstory drift ratio are characterized in terms of peak ground acceleration and spectral acceleration using a suite of ground motions representative of the seismic hazard in the region.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

When will a flame describing function approach to thermoacoustics work well?

In any thermoacoustic analysis, it is important not only to predict linear frequencies and growth rates, but also the amplitude and frequencies of any limit cycles. The Flame Describing Function (FDF) approach is a quasi-linear analysis which allows the prediction of both the linear and nonlinear behaviour of a thermoacoustic system. This means that one can predict linear growth rates and frequencies, and also the amplitudes and frequencies of any limit cycles. The FDF achieves this by assuming that the acoustics are linear and that the flame, which is the only nonlinear element in the thermoacoustic system, can be adequately described by considering only its response at the frequency at which it is forced. Therefore any harmonics generated by the flame's nonlinear response are not considered. This implies that these nonlinear harmonics are small or that they are sufficiently filtered out by the linear dynamics of the system (the low-pass filter assumption). In this paper, a flame model with a simple saturation nonlinearity is coupled to simple duct acoustics, and the success of the FDF in predicting limit cycles is studied over a range of flame positions and acoustic damping parameters. Although these two parameters affect only the linear acoustics and not the nonlinear flame dynamics, they determine the validity of the low-pass filter assumption made in applying the flame describing function approach. Their importance is highlighted by studying the level of success of an FDF-based analysis as they are varied. This is achieved by comparing the FDF's prediction of limit-cycle amplitudes to the amplitudes seen in time domain simulations.

Nonlinear drillstring dynamics analysis

This paper studies the dynamical response of a rotary drilling system with a drag bit, using a lumped parameter model that takes into consideration the axial and torsional vibration modes of the bit. These vibrations are coupled through a bit-rock interaction law. At the bit-rock interface, the cutting process introduces a state-dependent delay, while the frictional process is responsible for discontinuous right-hand sides in the equations governing the motion of the bit. This complex system is characterized by a fast axial dynamics compared to the slow torsional dynamics. A dimensionless formulation exhibits a large parameter in the axial equation, enabling a two-time-scales analysis that uses a combination of averaging methods and a singular perturbation approach. An approximate model of the decoupled axial dynamics permits us to derive a pseudoanalytical expression of the solution of the axial equation. Its averaged behavior influences the slow torsional dynamics by generating an apparent velocity weakening friction law that has been proposed empirically in earlier work. The analytical expression of the solution of the axial dynamics is used to derive an approximate analytical expression of the velocity weakening friction law related to the physical parameters of the system. This expression can be used to provide recommendations on the operating parameters and the drillstring or the bit design in order to reduce the amplitude of the torsional vibrations. Moreover, it is an appropriate candidate model to replace empirical friction laws encountered in torsional models used for control. © 2009 Society for Industrial and Applied Mathematics.

Nonlinear phenomena in thermoacoustic systems with premixed flames

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies, amplitudes and stability of limit cycles. Limit cycles in thermoacoustic systems are reached when the energy input from driving processes and energy losses from damping processes balance each other over a cycle of the oscillation. In this paper an integral relation for the rate of change of energy of a thermoacoustic system is derived. This relation is analogous to the well-known Rayleigh criterion in thermoacoustics, but can be used to calculate the amplitudes of limit cycles, as well as their stability. The relation is applied to a thermoacoustic system of a ducted slot-stabilized 2-D premixed flame. The flame is modelled using a nonlinear kinematic model based on the G-equation, while the acoustics of planar waves in the tube are governed by linearised momentum and energy equations. Using open-loop forced simulations, the flame describing function (FDF) is calculated. The gain and phase information from the FDF is used with the integral relation to construct a cyclic integral rate of change of energy (CIRCE) diagram that indicates the amplitude and stability of limit cycles. This diagram is also used to identify the types of bifurcation the system exhibits and to find the minimum amplitude of excitation needed to reach a stable limit cycle from another linearly stable state, for single- mode thermoacoustic systems. Furthermore, this diagram shows precisely how the choice of velocity model and the amplitudedependence of the gain and the phase of the FDF influence the nonlinear dynamics of the system. Time domain simulations of the coupled thermoacoustic system are performed with a Galerkin discretization for acoustic pressure and velocity. Limit cycle calculations using a single mode, as well as twenty modes, are compared against predictions from the CIRCE diagram. For the single mode system, the time domain calculations agree well with the frequency domain predictions. The heat release rate is highly nonlinear but, because there is only a single acoustic mode, this does not affect the limit cycle amplitude. For the twenty-mode system, however, the higher harmonics of the heat release rate and acoustic velocity interact resulting in a larger limit cycle amplitude. Multimode simulations show that in some situations the contribution from higher harmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Copyright © 2012 by ASME.

Fibre optic monitoring of a deep circular excavation

This paper describes part of the monitoring undertaken at Abbey Mills shaft F, one of the main shafts of Thames Water's Lee tunnel project in London, UK. This shaft, with an external diameter of 30 m and 73 m deep, is one of the largest ever constructed in the UK and consequently penetrates layered and challenging ground conditions (Terrace Gravel, London Clay, Lambeth Group, Thanet Sand Formation, Chalk Formation). Three out of the twenty 1-2 m thick and 84 m deep diaphragm wall panels were equipped with fibre optic instrumentation. Bending and circumferential hoop strains were measured using Brillouin optical time-domain reflectometry and analysis technologies. These measurements showed that the overall radial movement of the wall was very small. Prior to excavation during a dewatering trial, the shaft may have experienced three-dimensional deformation due to differential water pressures. During excavation, the measured hoop and bending strains of the wall in the chalk exceeded the predictions. This appears to be related to the verticality tolerances of the diaphragm wall and lower circumferential hoop stiffness of the diaphragm walls at deep depths. The findings from this case study provide valuable information for future deep shafts in London. © ICE Publishing: All rights reserved.

Ultrafast carrier dynamics in undoped and p-doped InAs/GaAs quantum dots characterized by pump-probe reflection measurements

We investigate the dependence of the differential reflection on the structure parameters of quantum dot (QD) heterostructures in pump-probe reflection measurements by both numerical simulations based on the finite-difference time-domain technique and theoretical calculations based on the theory of dielectric films. It is revealed that the value and sign of the differential reflection strongly depend on the thickness of the cap layer and the QD layer. In addition, a comparison between the carrier dynamics in undoped and p-doped InAs/GaAs QDs is carried out by pump-probe reflection measurements. The carrier capture time from the GaAs barrier into the InAs wetting layer and that from the InAs wetting layer into the InAs QDs are extracted by appropriately fitting differential reflection spectra. Moreover, the dependence of the carrier dynamics on the injected carrier density is identified. A detailed analysis of the carrier dynamics in the undoped and p-doped QDs based on the differential reflection spectra is presented, and its difference with that derived from the time-resolved photoluminescence is discussed. (C) 2008 American Institute of Physics.

Design and optimization of photonic crystal coupling layer for bi-color quantum well infrared photodetectors

Finite difference time domain (FDTD) method is used for the simulation and analysis of electromagnetic field in the top coupling layer of GaAs/AlGaAs quantum well infrared photodetector (QWIP). Simulation results demonstrated the coupling efficiencies and distributions of electromagnetic (EM) field in a variety of 2D photonic crystal coupling layer structures. A photonic crystal structure for bi-color-QWIP is demonstrated with high coupling efficiency for two wavelengths.

Numerical study of sub-wavelength plasmonic waveguide

The authors present an analysis of a plasmonic waveguide, simulated using a two-dimensional finite-difference time-domain technique. With the surface structures located on the surface of the metal, the device is able to confine and guide light waves in a sub-wavelength scale. And two waveguides can be placed within 150 nm (similar to 6% of the incident wavelength) that will helpful for the optoelectronic integration. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly two to four times the peak intensity of the incident light, and the propagation length can reach approximately 40 Pm at the wavelength of 2.44 mu m. (c) 2007 Elsevier B.V. All rights reserved.

Design of quantum cascade microcavity lasers based on Q factor versus etching depth

The choice of the etching depth for semiconductor microcavities is a compromise between a high Q factor and a difficult technique in a practical fabricating process. In this paper, the influences of the etching depth on mode Q factors for mid-infrared quantum cascade microcylinder and microsquare lasers around 4.8 and 7.8 mu m are simulated by three-dimensional (3D) finite-difference time-domain (FDTD) techniques. For the microcylinder and the microsquare resonators, the mode Q factors of the whispering-gallery modes (WGMs) increase exponentially and linearly with the increase in the etching depth, respectively Furthermore, the mode Q factors of some higher order transverse WGMs may be larger than that of the fundamental transverse WGM in 3D microsquares. Based on the field distribution of the vertical multilayer slab waveguide and the mode Q factors versus the etching depth, the necessary etching depth is chosen at the position where the field amplitude is 1% of the peak value of the slab waveguide. In addition, the influences of sidewall roughness on the mode Q factors are simulated for microsquare resonators by 2D FDTD simulation. (C) 2009 Optical Society of America

Resonant Enhanced Wave Filter and Waveguide via Surface Plasmons

The authors present an analysis of plasmonic wave filter and curved waveguide, simulated using a 2-D finite-difference time-domain technique. With different dielectric materials or surface structures located on the interface of the metal/dielectric, the resonant enhanced wave filter can divide light waves of different wavelengths and guide them with low losses. And the straight or curved waveguide can confine and guide light waves in a subwavelength scale. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly four times the peak intensity of the incident light.

Design, fabrication, and characterization of an ultracompact low-loss photonic crystal corner mirror

An ultracompact, low-loss, and broad-band corner mirror, based on photonic crystals, is investigated in this paper. Based on the theoretical analysis of the loss mechanism, the boundary layers of the photonic crystal region are revised to improve the extra losses, and the transmission characteristics are evaluated by using the 3-D finite-difference time-domain method. The device with optimized structure was fabricated on silicon-on-insulator substrate by using electron-beam lithography and inductively coupled plasma etching. The measured extra losses are about 1.1 +/- 0.4 dB per corner mirror for transverse-electronic polarization for the scanning wavelength range of 1510-1630 nm. Dimensions of the achieved PC corner mirror are less than ;7 x 7 mu m(2), which are only about one tenth of conventional wave-guide corner mirrors.

Prediction and suppression of strong dispersive coupling in microracetrack channel drop filters

For a four-port microracetrack channel drop filter, unexpected transmission characteristics due to strong dispersive coupling are demonstrated by the light tunneling between the input-output waveguides and the resonator, where a large dropping transmission at off-resonance wavelengths is observed by finite-difference time-domain simulation. It causes a severe decline of the extinction ratio and finesse. An appropriate decrease of the coupling strength is found to suppress the dispersive coupling and greately increase the extinction ratio and finesse, a decreased coupling strength can be realized by the application of an asymmetrical coupling waveguide structure. In addition, the profile of the coupling dispersion in the transmission spectra can be predicted based on a coupled mode theory analysis of an equivalent system consisting of two coupling straight waveguides. The effects of structure parameters on the transmission spectra obtained by this method agree well with the numerical results. It is useful to avoid the strong dispersive coupling region in the filter design. (c) 2007 Optical Society of America.