Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

A comparison of methods for the coupled analysis of floating structures

The dynamic analysis of a deepwater floating platform and the associated mooring/riser system should ideally be fully coupled to ensure a reliable response prediction. It is generally held that a time domain analysis is the only means of capturing the various coupling and nonlinear effects accurately. However, in recent work it has been found that for an ultra-deepwater floating system (2000m water depth), the highly efficient frequency domain approach can provide highly accurate response predictions. One reason for this is the accuracy of the drag linearization procedure over both first and second order motions, another reason is the minimal geometric nonlinearity displayed by the mooring lines in deepwater. In this paper, the aim is to develop an efficient analysis method for intermediate water depths, where both mooring/vessel coupling and geometric nonlinearity are of importance. It is found that the standard frequency domain approach is not so accurate for this case and two alternative methods are investigated. In the first, an enhanced frequency domain approach is adopted, in which line nonlinearities are linearized in a systematic way. In the second, a hybrid approach is adopted in which the low frequency motion is solved in the time domain while the high frequency motion is solved in the frequency domain; the two analyses are coupled by the fact that (i) the low frequency motion affects the mooring line geometry for the high frequency motion, and (ii) the high frequency motion affects the drag forces which damp the low frequency motion. The accuracy and efficiency of each of the methods are systematically compared. Copyright © 2007 by ASME.

Implementation of semi-active damping on a tri-axle heavy-vehicle suspension

A semi-active truck damper was developed in conjunction with a commercial shock absorber manufacturer. A linearized damper model was developed for control system design purposes. Open- and closed-loop damper force tracking control was implemented, with tests showing that an open-loop approach gave the best compromise between response speed and accuracy. A hardware-in-the-loop test facility was used to investigate performance of the damper when combined with a simulated quarter-car model. The input to the vehicle model was a set of randomly generated road profiles, each profile traversed at an appropriate speed. Modified skyhook damping tests showed a simultaneous improvement over the optimum passive case of 13 per cent in vertical body acceleration and 8 per cent in dynamic tyre forces. Full-scale vehicle tests of the damper on a heavy tri-axle trailer were carried out. Implementation of modified skyhook damping yielded a simultaneous improvement over the optimum passive case of 8 per cent in vertical body acceleration and 8 per cent in dynamic tyre forces. © IMechE 2008.

On separating propagating and non-propagating dynamics in fluid-flow equations

The ability to separate acoustically radiating and non-radiating components in fluid flow is desirable to identify the true sources of aerodynamic sound, which can be expressed in terms of the non-radiating flow dynamics. These non-radiating components are obtained by filtering the flow field. Two linear filtering strategies are investigated: one is based on a differential operator, the other employs convolution operations. Convolution filters are found to be superior at separating radiating and non-radiating components. Their ability to decompose the flow into non-radiating and radiating components is demonstrated on two different flows: one satisfying the linearized Euler and the other the Navier-Stokes equations. In the latter case, the corresponding sound sources are computed. These sources provide good insight into the sound generation process. For source localization, they are found to be superior to the commonly used sound sources computed using the steady part of the flow. Copyright © 2009 by S. Sinayoko, A. Agarwal, Z. Hu.

Reprint of: Flow filtering and the physical sources of aerodynamic sound

The physical sources of sound are expressed in terms of the non-radiating part of the flow. The non-radiating part of the flow can be obtained from convolution filtering, as we demonstrate numerically by using an axi-symmetric jet satisfying the Navier-Stokes equations. Based on the frequency spectrum of the source, we show that the sound sources exhibit more physical behaviour than sound sources based on acoustic analogies. To validate the sources of sound, one needs to let them radiate within the non-radiating flow field. However, our results suggest that the traditional Euler operator linearized about the time-averaged part of the flow should be sufficient to compute the sound field. © 2010 Published by Elsevier Ltd.

Experimental evidence of flutter and divergence instabilities induced by dry friction

Flutter and divergence instabilities have been advocated to be possible in elastic structures with Coulomb friction, but no direct experimental evidence has ever been provided. Moreover, the same types of instability can be induced by tangential follower forces, but these are commonly thought to be of extremely difficult, if not impossible, practical realization. Therefore, a clear experimental basis for flutter and divergence induced by friction or follower-loading is still lacking. This is provided for the first time in the present article, showing how a follower force of tangential type can be realized via Coulomb friction and how this, in full agreement with the theory, can induce a blowing-up vibrational motion of increasing amplitude (flutter) or an exponentially growing motion (divergence). In addition, our results show the limits of a treatment based on the linearized equations, so that nonlinearities yield the initial blowing-up vibration of flutter to reach eventually a steady state. The presented results give full evidence to potential problems in the design of mechanical systems subject to friction, open a new perspective in the realization of follower-loading systems and of innovative structures exhibiting 'unusual' dynamical behaviors. © 2011 Elsevier Ltd.

Decay of vibrations along deepwater mooring lines

Model tests for global design verification of deepwater floating structures cannot be made at reasonable scales. An overview of recent research efforts to tackle this challenge is given first, introducing the concept of line truncation techniques. In such a method the upper sections of each line are modelled in detail, capturing the wave action zone and all coupling effects with the vessel. These terminate to an approximate analytical model, that aims to simulate the remainder of the line. The rationale for this is that in deep water the transverse elastic waves of a line are likely to decay before they are reflected at the seabed. The focus of this paper is the verification of this rationale and the ongoing work, which is considering ways to produce a truncation model. Transverse dynamics of a mooring line are modelled using the equations of motion of an inextensible taut string, submerged in still water, one end fixed at the bottom the other assumed to follow the vessel response, which can be harmonic or random. Nonlinear hydrodynamic damping is included; bending and VIV effects are neglected. A dimensional analysis, supported by exact benchmark numerical solutions, has shown that it is possible to produce a universal curve for the decay of transverse vibrations along the line, which is suitable for any kind of line with any top motion. This has a significant engineering benefit, allowing for a rapid assessment of line dynamics - it is very useful in deciding whether a truncated line model is appropriate, and if so, at which point truncation might be applied. Initial efforts in developing a truncated model show that a linearized numerical solution in the frequency domain matches very closely the exact benchmark. Copyright © 2011 by ASME.

Structural sensitivity of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this. We investigate this spiral mode with a linear global stability analysis around the steady bubble and its wake. We obtain the linear direct and adjoint global modes of the linearized Navier-Stokes equations and overlap these to obtain the structural sensitivity of the spiral mode, which identifies the wavemaker region. We also identify regions of absolute instability with a local stability analysis. At moderate swirls, we find that the m=-1 azimuthal mode is the most unstable and that the wavemaker regions of the m=-1 mode lie around the bubble, which is absolutely unstable. The mode is most sensitive to feedback involving the radial and azimuthal components of momentum in the region just upstream of the bubble. To a lesser extent, the mode is also sensitive to feedback involving the axial component of momentum in regions of high shear around the bubble. At an intermediate swirl, in which the bubble and wake have similar absolute growth rates, other researchers have found that the wavemaker of the nonlinear global mode lies in the wake. We agree with their analysis but find that the regions around the bubble are more influential than the wake in determining the growth rate and frequency of the linear global mode. The results from this paper provide the first steps towards passive control strategies for spiral vortex breakdown. © 2013 Cambridge University Press.

Controlling the phase of an oscillator: A phase response curve approach

The paper discusses elementary control strategies to control the phase of an oscillator. Both feedforward and feedback (P and PI) control laws are designed based on the phase response curve (PRC) calculated from the linearized model. The performance is evaluated on a popular model of circadian oscillations. ©2009 IEEE.

Minimum maneuver time calculation using convex optimization

The problem of calculating the minimum lap or maneuver time of a nonlinear vehicle, which is linearized at each time step, is formulated as a convex optimization problem. The formulation provides an alternative to previously used quasi-steady-state analysis or nonlinear optimization. Key steps are: the use of model predictive control; expressing the minimum time problem as one of maximizing distance traveled along the track centerline; and linearizing the track and vehicle trajectories by expressing them as small displacements from a fixed reference. A consequence of linearizing the vehicle dynamics is that nonoptimal steering control action can be generated, but attention to the constraints and the cost function minimizes the effect. Optimal control actions and vehicle responses for a 90 deg bend are presented and compared to the nonconvex nonlinear programming solution. Copyright © 2013 by ASME.

Modal stability theory

This article contains a review of modal stability theory. It covers local stability analysis of parallel flows including temporal stability, spatial stability, phase velocity, group velocity, spatio-temporal stability, the linearized Navier-Stokes equations, the Orr-Sommerfeld equation, the Rayleigh equation, the Briggs-Bers criterion, Poiseuille flow, free shear flows, and secondary modal instability. It also covers the parabolized stability equation (PSE), temporal and spatial biglobal theory, 2D eigenvalue problems, 3D eigenvalue problems, spectral collocation methods, and other numerical solution methods. Computer codes are provided for tutorials described in the article. These tutorials cover the main topics of the article and can be adapted to form the basis of research codes. Copyright © 2014 by ASME.

Global modes with multiple saddle points

Significant progress has been made towards understanding the global stability of slowly-developing shear flows. The WKBJ theory developed by Patrick Huerre and his co-authors has proved absolutely central, with the result that both the linear and the nonlinear stability of a wide range of flows can now be understood in terms of their local absolute/convective instability properties. In many situations, the local absolute frequency possesses a single dominant saddle point in complex X-space (where X is the slow streamwise coordinate of the base flow), which then acts as a single wavemaker driving the entire global linear dynamics. In this paper we consider the more complicated case in which multiple saddles may act as the wavemaker for different values of some control parameter. We derive a frequency selection criterion in the general case, which is then validated against numerical results for the linearized third-order Ginzburg-Landau equation (which possesses two saddle points). We believe that this theory may be relevant to a number of flows, including the boundary layer on a rotating disk and the eccentric Taylor-Couette-Poiseuille flow. © 2014 Elsevier Masson SAS. All rights reserved.

Finding resonance: Adaptive frequency oscillators for dynamic legged locomotion

There is much to gain from providing walking machines with passive dynamics, e.g. by including compliant elements in the structure. These elements can offer interesting properties such as self-stabilization, energy efficiency and simplified control. However, there is still no general design strategy for such robots and their controllers. In particular, the calibration of control parameters is often complicated because of the highly nonlinear behavior of the interactions between passive components and the environment. In this article, we propose an approach in which the calibration of a key parameter of a walking controller, namely its intrinsic frequency, is done automatically. The approach uses adaptive frequency oscillators to automatically tune the intrinsic frequency of the oscillators to the resonant frequency of a compliant quadruped robot The tuning goes beyond simple synchronization and the learned frequency stays in the controller when the robot is put to halt. The controller is model free, robust and simple. Results are presented illustrating how the controller can robustly tune itself to the robot, as well as readapt when the mass of the robot is changed. We also provide an analysis of the convergence of the frequency adaptation for a linearized plant, and show how that analysis is useful for determining which type of sensory feedback must be used for stable convergence. This approach is expected to explain some aspects of developmental processes in biological and artificial adaptive systems that "develop" through the embodied system-environment interactions. © 2006 IEEE.