On the diffraction of acoustic waves by a quarter-plane

This paper follows the work of A.V. Shanin on diffraction by an ideal quarter-plane. Shanin's theory, based on embedding formulae, the acoustic uniqueness theorem and spherical edge Green's functions, leads to three modified Smyshlyaev formulae, which partially solve the far-field problem of scattering of an incident plane wave by a quarter-plane in the Dirichlet case. In this paper, we present similar formulae in the Neumann case, and describe a numerical method allowing a fast computation of the diffraction coefficient using Shanin's third modified Smyshlyaev formula. The method requires knowledge of the eigenvalues of the Laplace-Beltrami operator on the unit sphere with a cut, and we also describe a way of computing these eigenvalues. Numerical results are given for different directions of incident plane wave in the Dirichlet and the Neumann cases, emphasising the superiority of the third modified Smyshlyaev formula over the other two. © 2011 Elsevier B.V.

Adjoint algorithms for the Navier-Stokes equations in the low Mach number limit

This paper describes a derivation of the adjoint low Mach number equations and their implementation and validation within a global mode solver. The advantage of using the low Mach number equations and their adjoints is that they are appropriate for flows with variable density, such as flames, but do not require resolution of acoustic waves. Two versions of the adjoint are implemented and assessed: a discrete-adjoint and a continuous-adjoint. The most unstable global mode calculated with the discrete-adjoint has exactly the same eigenvalue as the corresponding direct global mode but contains numerical artifacts near the inlet. The most unstable global mode calculated with the continuous-adjoint has no numerical artifacts but a slightly different eigenvalue. The eigenvalues converge, however, as the timestep reduces. Apart from the numerical artifacts, the mode shapes are very similar, which supports the expectation that they are otherwise equivalent. The continuous-adjoint requires less resolution and usually converges more quickly than the discrete-adjoint but is more challenging to implement. Finally, the direct and adjoint global modes are combined in order to calculate the wavemaker region of a low density jet. © 2011 Elsevier Inc.

Oscillatory sensitivity patterns for global modes in wakes

Globally unstable wakes with co-flow at intermediate Reynolds numbers are studied, to quantify important spatial regions for the development and control of the global instability. One region of high structural sensitivity is found close to the inlet for all wakes, in agreement with previous findings for cylinder wakes. A second, elongated region of high structural sensitivity is seen downstream of the first one for unconfined wakes at Re = 400. When base flow modifications are considered, a spatially oscillating sensitivity pattern is found inside the downstream high structural sensitivity region. This implies that the same change in the base flow can either destabilize or stabilize the flow, depending on the exact position where it is applied. It is shown that the sensitivity pattern remains unchanged for different choices of streamwise boundary conditions and numerical resolution. Actual base flow modifications are applied in selected configurations, and the linear global modes recomputed. It is confirmed that the linear global eigenvalues move according to the predicted sensitivity pattern for small amplitude base flow modifications, for which the theory applies. We also look at the implications of a small control cylinder on the flow. Only the upstream high sensitivity region proves to be robust in terms of control, but one should be careful not to disturb the flow in the downstream high sensitivity region, in order to achieve control. The findings can have direct implications on the numerical resolution requirements for wakes at higher Reynolds numbers. Furthermore, they provide one more possible explanation to why confined wakes have a more narrow frequency spectrum than unconfined wakes.

Sensor graphs for guaranteed cooperative localization performance

A group of mobile robots can localize cooperatively, using relative position and absolute orientation measurements, fused through an extended Kalman filter (ekf). The topology of the graph of relative measurements is known to affect the steady-state value of the position error covariance matrix. Classes of sensor graphs are identified, for which tight bounds for the trace of the covariance matrix can be obtained based on the algebraic properties of the underlying relative measurement graph. The string and the star graph topologies are considered, and the explicit form of the eigenvalues of error covariance matrix is given. More general sensor graph topologies are considered as combinations of the string and star topologies, when additional edges are added. It is demonstrated how the addition of edges increases the trace of the steady-state value of the position error covariance matrix, and the theoretical predictions are verified through simulation analysis.

Development of realistic liquid equations of state for multi-phase CFD modelling

Several equations of state (EOS) have been incorporated into a novel algorithm to solve a system of multi-phase equations in which all phases are assumed to be compressible to varying degrees. The EOSs are used to both supply functional relationships to couple the conservative variables to the primitive variables and to calculate accurately thermodynamic quantities of interest, such as the speed of sound. Each EOS has a defined balance of accuracy, robustness and computational speed; selection of an appropriate EOS is generally problem-dependent. This work employs an AUSM+-up method for accurate discretisation of the convective flux terms with modified low-Mach number dissipation for added robustness of the solver. In this paper we show a newly-developed time-marching formulation for temporal discretisation of the governing equations with incorporated time-dependent source terms, as well as considering the system of eigenvalues that render the governing equations hyperbolic.

Swirling mean flow effect on fan-trailing edge broadband noise in a lined annular duct

The present study aims at investigating the effect of a swirling mean flow and a lined annular duct on rotor trailing-edge noise. The objectives are to investigate these effects on the eigenvalues and a tailored Green's function on one hand and on the realistic case of the fan trailing-edge noise on the other hand. Indeed, the mean flow in between the rotor and the stator of the fan is highly swirling. Moreover, interstage liners are used to reduce the noise produced by the fan stage. The extension of Ffowcs-Williams & Hawkings' acoustic analogy in a medium at rest with moving surfaces, of Goldstein's acoustic analogy in a hardwall circular duct with uniform mean flow and of Rienstra & Tester's Green's function in an annular lined duct with uniform mean flow to a swirling mean flow in an annular duct with liner is introduced. First, the eigenvalues and the Green's function are investigated showing a strong effect of the swirl and of the liner. Second, a rotor trailing-edge noise model accounting for both the effects of the annular duct with lined walls and the swirling mean flow is developed and applied to a realistic fan rotor with different swirling mean flows (and as a result different associated blade stagger angles). The benchmark cases are built from the Boeing 18-inch Fan Rig Broadband Noise Test. In all cases the swirling mean flow has a strong effect on the absolute noise level. The overall liner insertion loss is little changed by the swirl in the studied cases.

Shock amplification, curve veering and the role of damping

The circumstances are investigated under which high peak acceleration can occur in the internal parts of a system when subjected to impulsive driving on the outside. Previous work using a coupled beam model has highlighted the importance of veering pairs of modes. Such a veering pair can be approximated by a lumped system with two degrees of freedom. The worst case of acceleration amplification is shown to occur when the two oscillators are tuned to the same frequency, and for this case closed-form expressions are derived to show the parameter dependence of the acceleration ratio on the mass ratio and coupling strength. Sensitivity analysis of the eigenvalues and eigenvectors indicates that mass ratio is the most sensitive parameter for altering the veering behaviour in an undamped system. Non-proportional damping is also shown to have a strong influence on the veering behaviour. The study gives design guidelines to allow permissible acceleration levels to be achieved by the choice of the effective mass and damping of the indirectly driven subsystem relative to the directly driven subsystem. © 2013 Elsevier Ltd.