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.

Long range exhaustion-a mathematical model for the axisymmetric air flow of a local exhaust ventilation hood assisted by a turbulent radial jet

Abstract-Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art Banged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air Bow to be reduced and provide a means of identifying and assessing the various parameters that control the air Bow. The mathematical model is formulated in terms of the Stokes steam function, ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained Bow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions. | Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art flanged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air flow to be reduced and provide a means of identifying and assessing the various parameters that control the air flow. The mathematical model is formulated in terms of the Stokes steam function, Ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained flow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions.

A body force-based method for prediction of multiple-pure-tone noise: Validation

Emissions, fuel burn, and noise are the main drivers for innovative aircraft design. Embedded propulsion systems, such as for example used in hybrid-wing body aircraft, can offer fuel burn and noise reduction benefits but the impact of inlet flow distortion on the generation and propagation of turbomachinery noise has yet to be assessed. A novel approach is used to quantify the effects of non-uniform flow on the creation and propagation of multiple pure tone (MPT) noise. The ultimate goal is to conduct a parametric study of S-duct inlets to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the effects of distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required to capture the mechanisms at play. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the blade-to-blade flow variations that cause the MPT noise while reducing computational cost. A single, 3-D full-wheel CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted inlet flow. A new method of producing the blade-to-blade variations in the body force field for MPT noise generation has been developed and validated. The numerical dissipation inherent to the solver is quantified and used to correct for non-physical attenuation in the far-field noise spectra. Source generation, acoustic propagation and acoustic energy transfer between modes is examined in detail. The new method is validated on NASA's Source Diagnostic Test fan and inlet, showing good agreement with experimental data for aerodynamic performance, acoustic source generation, and far-field noise spectra. The next steps involve the assessment of MPT noise in serpentine inlet ducts and the development of a reduced order formulation suitable for incorporation into NASA's ANOPP framework. © 2010 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

A novel MPT noise prediction methodology for highly-integrated propulsion systems with inlet flow distortion

Embedded propulsion systems, such as for example used in advanced hybrid-wing body aircraft, can potentially offer major fuel burn and noise reduction benefits but introduce challenges in the aerodynamic and acoustic integration of the high-bypass ratio fan system. A novel approach is proposed to quantify the effects of non-uniform flow on the generation and propagation of multiple pure tone noise (MPTs). The new method is validated on a conventional inlet geometry first. The ultimate goal is to conduct a parametric study of S-duct inlets in order to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the mechanism underlying the distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the MPT noise generation mechanisms while greatly reducing computational cost. A single, 3-D full-wheel unsteady CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted mean flow. Several numerical tools were developed to enable the implementation of this new approach. Parametric studies were conducted to determine appropriate grid and time step sizes for the propagation of acoustic waves. The Ffowcs-Williams and Hawkings integral method is used to propagate the noise to far field receivers. Non-reflecting boundary conditions are implemented through the use of acoustic buffer zones. The body force modeling approach is validated and proof-of-concept studies demonstrate the generation of disturbances at both blade-passing and shaft-order frequencies using the perturbed body force method. The full methodology is currently being validated using NASA's Source Diagnostic Test (SDT) fan and inlet geometry. Copyright © 2009 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

The dynamic response of edge clamped plates loaded by spherically expanding sand shells

The dynamic deformation of both edge clamped stainless steel sandwich panels with a pyramidal truss core and equal mass monolithic plates loaded by spherically expanding shells of dry and water saturated sand has been investigated, both experimentally and via a particle based simulation methodology. The spherically expanding sand shell is generated by detonating a sphere of explosive surrounded by a shell of either dry or water saturated synthetic sand. The measurements show that the sandwich panel and plate deflections decrease with increasing stand-off between the center of the charge and the front of the test structures. Moreover, for the same charge and sand mass, the deflections of the plates are significantly higher in the water saturated sand case compared to that of dry sand. For a given stand-off, the mid-span deflection of the sandwich panel rear faces was substantially less than that of the corresponding monolithic plate for both the dry and water saturated sand cases. The experiments were simulated via a coupled discrete-particle/ finite element scheme wherein the high velocity impacting sand is modeled by interacting particles while the plate is modeled within a Lagrangian finite element setting. The simulations are in good agreement with the measurements for the dry sand impact of both the monolithic and sandwich structures. However, the simulations underestimate the effect of stand-off in the case of the water saturated sand explosion, i.e. the deflections decrease more sharply with increasing stand-off in the experiments compared to the simulations. The simulations reveal that the momentum transmitted into the sandwich and monolithic plate structures by the sand shell is approximately the same, consistent with a small fluid-structure interaction effect. The smaller deflection of the sandwich panels is therefore primarily due to the higher bending strength of sandwich structures. © 2013 The Authors. Published by Elsevier Ltd. All rights reserved.

Hybrid microstrip and carbon nanotubes based patch antenna for wireless applications

In this we have looked at the concept of introducing carbon nanotubes on the surfaces of the microstrip patch antennas. We examined the performance improvements in a patch antenna through finite difference time domain simulations to increase the efficiency of the antenna. The results suggest that carbon nanotubes lead to a higher gain due to their electrical properties. A high gain antenna with low power requirements resulted in achieving a higher overall bandwidth. The designed antenna's gain, bandwidth and directivity are analyzed before and after introducing carbon nanotubes. © 2013 IEEE.

Control of material transport through pulse shape manipulation-a development toward designer pulses

The variety of laser systems available to industrial laser users is growing and the choice of the correct laser for a material target application is often based on an empirical assessment. Industrial master oscillator power amplifier systems with tuneable temporal pulse shapes have now entered the market, providing enormous pulse parameter flexibility in an already crowded parameter space. In this paper, an approach is developed to design interaction parameters based on observations of material responses. Energy and material transport mechanisms are studied using pulsed digital holography, post process analysis techniques and finite-difference modelling to understand the key response mechanisms for a variety of temporal pulse envelopes incident on a silicon (1/1/1) substrate. The temporal envelope is shown to be the primary control parameter of the source term that determines the subsequent material response and the resulting surface morphology. A double peak energy-bridged temporal pulse shape designed through direct application of holographic imaging data is shown to substantially improve surface quality. © 2014 IEEE.

Mode Simulation for Midinfrared Microsquare Resonators With Sloped Sidewalls and Confined Metals

EQUILATERAL-TRIANGLE; MU-M; LASERS; MICROLASERS; MICRODISK Abstract: Mode characteristics for midinfrared microsquare resonators with sloped sidewalls and confined metal layers are investigated by finite-difference time-domain (FDTD) techniques. For a microsquare with a side length of 10 mu m, the mode quality (Q)-factors of 8329, 4772, and 2053 are obtained for TM5,7 mode at wavelength 7.1 mu m by three-dimensional FDTD simulations, as the tilting angles of the side walls are 90 degrees, 88 degrees, and 86 degrees, respectively. Furthermore, microsquare resonators laterally surrounded by SiO2 and metal layers are investigated by the two-dimensional FDTD technique for the metal layers of Au, Ti-Au, Ag-Au, and Ti-Ag-Au, respectively.

Multimode resonances in metallically confined square-resonator microlasers

Directional emission InP/AlGaInAs square-resonator microlasers with a side length of 20 mu m are fabricated by standard photolithography and inductively coupled-plasma etching technique. Multimode resonances with about seven distinct mode peaks in a free-spectral range are observed from 1460 to 1560 nm with the free-spectral range of 12.1 nm near the wavelength of 1510 nm, and the mode refractive index versus the photon energy E (eV) as 3.07152+0.18304E are obtained by fitting the laser spectra with an analytical mode wavelength formula derived by light ray method. In addition, mode field pattern is simulated for cold cavity by two dimensional finite-difference time-domain technique.

Mode Characteristics for Square Resonators With a Metal Confinement Layer

Microsquare resonators laterally confined by SiO2/Au/air multilayer structure are investigated by light ray method with reflection phase-shift of the multiple layers and two-dimensional (2-D) finite-difference time-domain (FDTD) technique. The reflectivity and phase shift of the mode light ray on the sides of the square resonator with the semiconductor/SiO2/Au/air multilayer structure are calculated for TE and TM modes by transfer matrix method. Based on the reflection phase shift and the reflectivity, the mode wavelength and factor are calculated by the resonant condition and the mirror loss, which are in agreement well with that obtained by the FDTD simulation. We find that the mode factor increases greatly with the increase of the SiO2 layer thickness, especially as d < 0.3 mu m. For the square resonator with side length 2 mu m and refractive index 3.2, anticrossing mode couplings are found for confined TE modes at wavelength about 1.6 mu m at d = 0.11 mu m, and confined TM modes at d = 0.71 mu m, respectively.

Design of a compact polarizing beam splitter based on a photonic crystal ring resonator with a triangular lattice

We propose and simulate a new kind of compact polarizing beam splitter (PBS) based on a photonic crystal ring resonator (PCRR) with complete photonic bandgaps. The two polarized states are separated far enough by resonant and nonresonant coupling between the waveguide modes and the microring modes. Some defect holes are utilized to control the beam propagation. The simulated results obtained by the finite-difference time-domain method show that high transmission (over 95%) is obtained and the polarization separation is realized with a length as short as 3.1 mu m. The design of the proposed PBS can be flexible, thanks to the advantages of PCRRs.

Investigation of mode radiation loss for microdisk resonators with pedestals by FDTD technique

Mode radiation loss for microdisk resonators with pedestals is investigated by three-dimensional (3D) finite-difference time-domain (FDTD) technique. For the microdisk with a radius of 1 mu m, a thickness of 0.2 mu m, and a refractive index of 3.4, on a pedestal with a refractive index of 3.17, the mode quality (Q) factor of the whispering-gallery mode (WGM) quasi-TE7,1 first increases with the increase of the radius of the pedestal, and then quickly decreases as the radius is larger than 0.75 mu m. The mode radiation loss is mainly the vertical radiation loss induced by the mode coupling between the WGM and vertical radiation mode in the pedestal, instead of the scattering loss around the perimeter of the round pedestal. The WG M can keep the high Q factor when the mode coupling is forbidden.

Enhancement of quality factor for TE whispering-gallery modes in microcylinder resonators

The enhancement of quality factor for TE whispering-gallery modes is analyzed for three-dimensional microcylinder resonators based on the destructive interference between vertical leakage modes. In the microcylinder resonator, the TE whispering-gallery modes can couple with vertical propagation modes, which results in vertical radiation loss and low quality factors. However, the vertical loss can be canceled by choosing appropriate thickness of the upper cladding layer or radius of the microcylinder. A mode quality factor increase by three orders of magnitude is predicted by finite-difference time-domain simulation. Furthermore, the condition of vertical leakage cancellation is analyzed.

Mode characteristics of metallically coated square microcavity connected with an output waveguide

Mode characteristics of a square microcavity with an output waveguide on the middle of one side, laterally confined by an insulating layer SiO2 and a p-electrode metal Au, are investigated by two-dimensional finite-difference time-domain technique. The mode quality (Q) factors versus the width of the output waveguide are calculated for Fabry-Peacuterot type and whispering-gallery type modes in the square cavity. Mode coupling between the confined modes in the square cavity and the guided modes in the output waveguide determines the mode Q factors, which is greatly influenced by the symmetry behaviors of the modes. Fabry-Peacuterot type modes can also have high Q factors due to the high reflectivity of the Au layer for the vertical incident mode light rays. For the square cavity with side length 4 mu m and refractive index 3.2, the mode Q factors of the Fabry-Peacuterot type modes can reach 10(4) at the mode wavelength of 1.5 mu m as the output waveguide width is 0.4 mu m.

Analysis of coupled microcircular resonators coupled to a bus waveguide with high output efficiency

Coupled microcircular resonators tangentially coupled to a bus waveguide, which is between the resonators, are numerically investigated by the finite-difference time-domain technique. For symmetrically coupled microcircular resonators with refractive index of 3.2, radius of 2 mu m, and width of the bus waveguide of 0.4 mu m, a mode Q factor of the order of 105 is obtained for a mode at the frequency of 243 THz. An output coupling efficiency of as high as 0.99 is calculated for a mode with a Q factor ranging from 10(3) to 10(4). The mode Q factor is 2 orders larger than that of the modes confined in a single circular resonator tangentially coupled to the same bus waveguide. Furthermore, the high Q traveling modes in the coupled microcircular resonators are suitable for optical single processing.