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.

A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

A critical analysis of loudness calculation methods

The physical meaning and methods of determining loudness were reviewed Loudness is a psychoacoustic metric which closely corresponds to the perceived intensity of a sound stimulus. It can be determined by graphical procedures, numerical methods, or by commercial software. These methods typically require the consideration of the 1/3 octave band spectrum of the sound of interest. The sounds considered in this paper are a 1 kHz tone and pink noise. The loudness of these sounds was calculated in eight ways using different combinations of input data and calculation methods. All the methods considered are based on Zwicker loudness. It was determined that, of the combinations considered, only the commercial software dBSonic and the loudness calculation procedure detailed in DIN 45631 using 1/3 octave band levels filtered using ANSI S1.11-1986 gave the correct values of loudness for a 1 kHz tone. Comparing the results between the sources also demonstrated the difference between sound pressure level and loudness. It was apparent that the calculation and filtering methods must be considered together, as a given calculation will produce different results for different 1/3 octave band input. In the literature reviewed, no reference provided a guide to the selection of the type of filtering that should be used in conjunction with the loudness computation method.

Phase noise analysis of a mechanical autonomous impact oscillator with a MEMS resonator

In this paper, phase noise analysis of a mechanical autonomous impact oscillator with a MEMS resonator is performed. Since the circuit considered belongs to the class of hybrid systems, methods based on the variational model for the evaluation of either phase noise or steady state solutions cannot be directly applied. As a matter of fact, the monodromy matrix is not defined at impact events in these systems. By introducing saltation matrices, this limit is overcome and the aforementioned methods are extended. In particular, the unified theory developed by Demir is used to analyze the phase noise after evaluating the asymptotically stable periodic solution of the system by resorting to the shooting method. Numerical results are presented to show how noise sources affect the phase noise performances. © 2011 IEEE.

Effects of initial soil fabric and mode of shearing on quasi-steady state line for monotonic undrained behavior

The effects of initial soil fabric and mode of shearing on quasi-steady state line in void ratiostress space are studied by employing the Distinct Element Method numerical analysis. The results show that the initial soil fabric and the mode of shearing have a profound effect on the location of the quasi-steady state line. The evolution of the soil fabric during the course of undrained shearing shows that the specimens with different initial soil fabrics reach quasi-steady state at various soil fabric conditions. At quasi-steady state, the soil fabric has a significant adjustment to change its behavior from contractive to dilative. As the stress state approaches the steady state, the soil fabrics of different initial conditions become similar. The numerical analysis results are compared qualitatively with the published experimental data and the effects of specimen reconstitution methods and mode of shearing found in the experimental studies canbe systematically explained by the numerical analysis. © 2009 Taylor & Francis Group.

LES of heat transfer in electronics

Numerical methods based on the Reynolds Averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) equations are applied to the thermal prediction of flows representative of those found in and around electronics systems and components. Low Reynolds number flows through a heated ribbed channel, around a heated cube and within a complex electronics system case are investigated using linear and nonlinear LES models, hybrid RANS-LES and RANS-Numerical-LES (RANS-NLES) methods. Flow and heat transfer predictions using these techniques are in good agreement with each other and experimental data for a range of grid resolutions. Using second order central differences, the RANS-NLES method performs well for all simulations. © 2011 Elsevier Inc.

Review of computational aeroacoustics for application in electronics cooler noise

A review of computational aeroacoustics (CCA) was made for application in electronics cooler noise. Computational aeroacoustics encompasses all numerical methods where the purposes is to predict the noise emissions from a simulated flow. Numerical simulation of the flow inside and around heat sinks and fans can lead to a prediction of the emitted noise while they are still in the design phase. Direct CCA is theoretically the best way to predict flow-based acoustic phenomena numerically. It is typically used only for low-frequency sound prediction. The boundary element method offers low computational cost and does not use a computational grid, but instead use vortex-surface calculations to determine tonal noise. Axial fans are commonly used to increase the airflow and thus the heat transfer over the heat sinks within the computer cases. Very detailed source simulations in the fan and heat sink region coupled with the use of analogy methods could result in excellent simulation results with a reasonable computational effort.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

Autonomous extraction of optimal flame fronts in OH planar laser-induced fluorescence images.

The location of a flame front is often taken as the point of maximum OH gradient. Planar laser-induced fluorescence of OH can be used to obtain the flame front by extracting the points of maximum gradient. This operation is typically performed using an edge detection algorithm. The choice of operating parameters a priori poses significant problems of robustness when handling images with a range of signal-to-noise ratios. A statistical method of parameter selection originating in the image processing literature is detailed, and its merit for this application is demonstrated. A reduced search space method is proposed to decrease computational cost and render the technique viable for large data sets. This gives nearly identical output to the full method. These methods demonstrate substantial decreases in data rejection compared to the use of a priori parameters. These methods are viable for any application where maximum gradient contours must be accurately extracted from images of species or temperature, even at very low signal-to-noise ratios.

Progress towards the design and numerical analysis of a 3D microchannel biochip separator

This paper reports the design and numerical analysis of a three-dimensional biochip plasma blood separator using computational fluid dynamics techniques. Based on the initial configuration of a two-dimensional (2D) separator, five three-dimensional (3D) microchannel biochip designs are categorically developed through axial and plenary symmetrical expansions. These include the geometric variations of three types of the branch side channels (circular, rectangular, disc) and two types of the main channel (solid and concentric). Ignoring the initial transient behaviour and assuming that steady-state flow has been established, the behaviour of the blood fluid in the devices is algebraically analysed and numerically modelled. The roles of the relevant microchannel mechanisms, i.e. bifurcation, constriction and bending channel, on promoting the separation process are analysed based on modelling results. The differences among the different 3D implementations are compared and discussed. The advantages of 3D over 2D separator in increasing separation volume and effectively depleting cell-free layer fluid from the whole cross section circumference are addressed and illustrated. © 2011 John Wiley & Sons, Ltd.