Plasmonic band gaps and waveguide effects in carbon nanotube arrays based metamaterials.

Highly dense periodic arrays of multiwalled carbon nanotubes behave like low-density plasma of very heavy charged particles, acting as metamaterials. These arrays with nanoscale lattice constants can be designed to display extended plasmonic band gaps within the optical regime, encompassing the crucial optical windows (850 and 1550 nm) simultaneously. We demonstrate an interesting metamaterial waveguide effect displayed by these nanotube arrays containing line defects. The nanotube arrays with lattice constants of 400 nm and radius of 50 nm were studied. Reflection experiments conducted on the nanoscale structures were in agreement with numerical calculations.

Model Refinement for the Active Control of Thermoacoustic Instability

Active control has been shown as a feasible technology for suppressing thermoacoustic instability in continuous combustion systems, and the control strategy design is substantially dependent on the reliability of the flame model. In this paper, refinement of G-equation flame model for the dynamics of lean premixed combustion is investigated. Precisely, the dynamics between the flame speed S_u and equivalence ratio phi are proposed based on numerical calculations and physical explanations. Finally, the developed model is tested on one set of experimental data.

Metamaterial filter for the near-visible spectrum

We demonstrate metamaterials operating in the near-visible regime based on two-dimensional arrays of gold-coated silicon nanopillars. The nanopillar arrays demonstrate a cutoff response at the metamaterial plasma frequency in accordance with theory and can be utilized for filtering applications. A plasma frequency in the near visible region of λ = 1 μm is calculated numerically for an array with a lattice constant of 300 nm and wire radius of 50 nm, with reflection measurements in agreement with numerical calculations. Such structures can be utilized for achieving negative-index based metamaterials for the visible spectrum. © 2012 American Institute of Physics.

Calibration of the spring constant of cantilevers of arbitrary shape using the phase signal in an atomic force microscope.

The measurement of cantilever parameters is an essential part of performing a calibrated measurement with an atomic force microscope (AFM). The thermal motion method is a widely used technique for calibrating the spring constant of an AFM cantilever, which can be applied to non-rectangular cantilevers. Given the trend towards high frequency scanning, calibration of non-rectangular cantilevers is of increasing importance. This paper presents two results relevant to cantilever calibration via the thermal motion method. We demonstrate the possibility of using the AFM's phase signal to acquire the thermal motion. This avoids the challenges associated with connecting the raw photodiode signal to a separate spectrum analyser. We also describe how numerical calculations may be used to calculate the parameters needed in a thermal motion calibration of a non-rectangular cantilever. Only accurate knowledge of the relative size of the in-plane dimensions of the cantilever is needed in this computation. We use this pair of results in the calibration of a variety of rectangular and non-rectangular cantilevers. We observe an average difference between the Sader and thermal motion values of cantilever stiffness of 10%.

Metamaterial filter for the near-visible spectrum

We demonstrate metamaterials operating in the near-visible regime based on two-dimensional arrays of gold-coated silicon nanopillars. The nanopillar arrays demonstrate a cutoff response at the metamaterial plasma frequency in accordance with theory and can be utilized for filtering applications. A plasma frequency in the near visible region of λ = 1 μm is calculated numerically for an array with a lattice constant of 300 nm and wire radius of 50 nm, with reflection measurements in agreement with numerical calculations. Such structures can be utilized for achieving negative-index based metamaterials for the visible spectrum.

Plasmonic band gaps and waveguide effects in carbon nanotube arrays based metamaterials

Highly dense periodic arrays of multiwalled carbon nanotubes behave like low-density plasma of very heavy charged particles, acting as metamaterials. These arrays with nanoscale lattice constants can be designed to display extended plasmonic band gaps within the optical regime, encompassing the crucial optical windows (850 and 1550 nm) simultaneously. We demonstrate an interesting metamaterial waveguide effect displayed by these nanotube arrays containing line defects. The nanotube arrays with lattice constants of 400 nm and radius of 50 nm were studied. Reflection experiments conducted on the nanoscale structures were in agreement with numerical calculations.

Direct comparison of the gravimetric responsivities of ZnO-based FBARs and SMRs

Film bulk acoustic resonators (FBARs) and solidly mounted resonators (SMRs) have the potential to significantly improve upon the sensitivity and minimum detection limit of traditional gravimetric sensors based on quartz crystal microbalances (QCMs) and surface acoustic wave resonators (SAWs). To date, neither FBAR nor SMR devices have been demonstrated to be superior to the other; hence the choice between them depends primarily on the users' ability to design/fabricate membranes and/or Bragg reflectors. In this work, it is shown that identically designed FBAR and SMR devices resonating at the same frequency exhibit different responsivities to mass loadings, Rm, and that the SMRs are less responsive than the FBARs. For the specific device design and resonant frequency (~2 GHz) of the resonators presented here, the FBARs' mass responsivity is ~20% greater than that of the SMRs', and although this value is not universal for all possible device designs, it clearly shows that FBAR devices should be favoured over SMRs in gravimetric sensing applications where the FBARs' fragility is not an issue. Numerical calculations based on Mason's model offer an insight into the physical mechanisms behind the greater FBARs responsivity, and it was shown that the Bragg reflector has an effect on the acoustic load at one of the facets of the piezoelectric films which is in turn responsible for the SMRs' lower responsivity to mass loadings. © 2013 Elsevier B.V.

Dynamic compression of foam supported plates impacted by high velocity soil

The response of back-supported buffer plates comprising a solid face sheet and foam core backing impacted by a column of high velocity particles (sand slug) is investigated via a lumped parameter model and coupled discrete/continuum simulations. The buffer plate is either resting on (unattached) or attached to a rigid stationary foundation. The lumped parameter model is used to construct maps of the regimes of behaviour with axes of the ratio of the height of the sand slug to core thickness and the normalised core strength. Four regimes of behaviour are identified based on whether the core compression ends prior to the densification of the sand slug or vice versa. Coupled discrete/continuum simulations are also reported and compared with the lumped parameter model. While the model predicted regimes of behaviour are in excellent agreement with numerical simulations, the lumped parameter model is unable to predict the momentum transmitted to the supports as it neglects the role of elasticity in both the buffer plate and the sand slug. The numerical calculations show that the momentum transfer is minimised for intermediate values of the core strength when the so-called "soft-catch" mechanism is in play. In this regime the bounce-back of the sand slug is minimised which reduces the momentum transfer. However, in this regime, the impulse reduction is small (less than 10% of that transferred to a rigid structure). For high values of the core strength, the response of the buffer plate resembles a rigid plate with nearly no impulse mitigation while at low values of core strength, a slap event occurs when the face sheet impinges against the foundation due to full densification of the foam core. This slap event results in a significant enhancement of the momentum transfer to the foundation. The results demonstrate that appropriately designed buffer plates have potential as impulse mitigators in landmine loading situations. © 2013 Elsevier Ltd. All rights reserved.

Numerical stability of the predictor-corrector method in Monte Carlo burnup calculations of critical reactors

Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the simplest methods, such as the beginning-of-step and middle-of-step constant flux approximations, are numerically unstable in fuel cycle calculations of critical reactors. Here we show that even the predictor-corrector methods that are implemented in established Monte Carlo burnup codes can be numerically unstable in cycle calculations of large systems. © 2013 Elsevier Ltd. All rights reserved.

Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

A matrix-calculation-based algorithm for numerical change propagation analysis

Engineering changes (ECs) are raised throughout the lifecycle of engineering products. A single change to one component produces knock-on effects on others necessitating additional changes. This change propagation significantly affects the development time and cost and determines the product's success. Predicting and managing such ECs is, thus, essential to companies. Some prediction tools model change propagation by algorithms, whereof a subgroup is numerical. Current numerical change propagation algorithms either do not account for the exclusion of cyclic propagation paths or are based on exhaustive searching methods. This paper presents a new matrix-calculation-based algorithm which can be applied directly to a numerical product model to analyze change propagation and support change prediction. The algorithm applies matrix multiplications on mutations of a given design structure matrix accounting for the exclusion of self-dependences and cyclic propagation paths and delivers the same results as the exhaustive search-based Trail Counting algorithm. Despite its factorial time complexity, the algorithm proves advantageous because of its straightforward matrix-based calculations which avoid exhaustive searching. Thereby, the algorithm can be implemented in established numerical programs such as Microsoft Excel which promise a wider application of the tools within and across companies along with better familiarity, usability, practicality, security, and robustness. © 1988-2012 IEEE.

An investigation into numerical analysis alternatives for predicting reingestion in turbine disc rim cavities

Reliable means of predicting ingestion in cavities adjacent to the main gas path are increasingly being sought by engineers involved in the design of gas turbines. In this paper, analysis is to be presented that results from an extended research programme, MAGPI, sponsored by the EU and several leading gas turbine manufactures and universities. Extensive use is made of CFD modelling techniques to understand the aerodynamic behaviour of a turbine stator well cavity, focusing on the interaction of cooling air supply with the main annulus gas. The objective of the study has been to benchmark a number of CFD codes and numerical techniques covering RANS and URANS calculations with different turbulence models in order to assess the suitability of the standard settings used in the industry for calculating the mechanics of the flow travelling between cavities in a turbine through the main gas path. The modelling methods employed have been compared making use of experimental data gathered from a dedicated two-stage turbine rig, running at engine representative conditions. Extensive measurements are available for a range of flow conditions and alternative cooling arrangements. The limitations of the numerical methods in calculating the interaction of the cooling flow egress and the main stream gas, and subsequent ingestion into downstream cavities in the engine (i.e. re-ingestion), have been exposed. This has been done without losing sight of the validation of the CFD for its use for predicting heat transfer, which was the main objective of the partners of the MAGPI Work- Package 1 consortium. Copyright © 2012 by ASME.

Monitoring and preventing numerical oscillations in 3D simulations with coupled Monte Carlo codes

Previous studies have reported that different schemes for coupling Monte Carlo (MC) neutron transport with burnup and thermal hydraulic feedbacks may potentially be numerically unstable. This issue can be resolved by application of implicit methods, such as the stochastic implicit mid-point (SIMP) methods. In order to assure numerical stability, the new methods do require additional computational effort. The instability issue however, is problem-dependent and does not necessarily occur in all cases. Therefore, blind application of the unconditionally stable coupling schemes, and thus incurring extra computational costs, may not always be necessary. In this paper, we attempt to develop an intelligent diagnostic mechanism, which will monitor numerical stability of the calculations and, if necessary, switch from simple and fast coupling scheme to more computationally expensive but unconditionally stable one. To illustrate this diagnostic mechanism, we performed a coupled burnup and TH analysis of a single BWR fuel assembly. The results indicate that the developed algorithm can be easily implemented in any MC based code for monitoring of numerical instabilities. The proposed monitoring method has negligible impact on the calculation time even for realistic 3D multi-region full core calculations. © 2014 Elsevier Ltd. All rights reserved.