Global nonlinear dynamics of thin aerofoil wakes

In the present investigation of thin aerofoil wakes we compare the global nonlinear dynamics, obtained by direct numerical simulations, to the associated local instability features, derived from linear stability analyses. A given configuration depends on two control parameters: the Reynolds number Re and the adverse pressure gradient m (with m < 0) prevailing at the aerofoil trailing edge. Global instability is found to occur for large enough Re and |m|; the naturally selected frequency is determined by the local absolute frequency prevailing at the trailing edge. © 2010 Springer Science+Business Media B.V.

Linear and nonlinear model large-eddy simulations of a plane jet

The nonlinear Kosovic, and mixed Leray and α subgrid scale models are contrasted with linear Smagorinsky and Yoshizawa Large Eddy Simulations for a Re = 4000 plane jet simulation. Comparisons are made with Direct Numerical Simulation data and measurements. Global properties of the jet such as the spreading and centreline velocity decay rates are investigated. The mean-flow and turbulence parameters in the self-similar region are also studied. All models generally give encouraging agreement with the Direct Numerical Simulation data and reliable measurements. Solution differences for the models are relatively minor, none giving clear improvements for all data comparisons.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Nonlinear optics in Silicon photonic crystal cavities

Silicon is known to be a very good material for the realization of high-Q, low-volume photonic cavities, but at the same it is usually considered as a poor material for nonlinear optical functionalities like second-harmonic generation, because its second-order nonlinear susceptibility vanishes in the dipole approximation. In this work we demonstrate that nonlinear optical effects in silicon nanocavities can be strongly enhanced and even become macroscopically observable. We employ photonic crystal nanocavities in silicon membranes that are optimized simultaneously for high quality factor and efficient coupling to an incoming beam in the far field. Using a low-power, continuous-wave laser at telecommunication wavelengths as a pump beam, we demonstrate simultaneous generation of second- and third harmonics in the visible region, which can be observed with a simple camera. The results are in good agreement with a theoretical model that treats third-harmonic generation as a bulk effect in the cavity region, and second-harmonic generation as a surface effect arising from the vertical hole sidewalls. Optical bistability is also observed in the silicon nanocavities and its physical mechanisms (optical, due to two-photon generation of free carriers, as well as thermal) are investigated. © 2011 IEEE.

Modelling and control of nonlinear systems using Gaussian processes with partial model information

Gaussian processes are gaining increasing popularity among the control community, in particular for the modelling of discrete time state space systems. However, it has not been clear how to incorporate model information, in the form of known state relationships, when using a Gaussian process as a predictive model. An obvious example of known prior information is position and velocity related states. Incorporation of such information would be beneficial both computationally and for faster dynamics learning. This paper introduces a method of achieving this, yielding faster dynamics learning and a reduction in computational effort from O(Dn2) to O((D - F)n2) in the prediction stage for a system with D states, F known state relationships and n observations. The effectiveness of the method is demonstrated through its inclusion in the PILCO learning algorithm with application to the swing-up and balance of a torque-limited pendulum and the balancing of a robotic unicycle in simulation. © 2012 IEEE.

Studying adsorbent dynamics on a quartz crystal resonator using its nonlinear electrical response

The quartz crystal resonator has been traditionally employed in studying surface-confined physisorbed films and particles by measuring dissipation and frequency shifts. However, theoretical interpretation of the experimental observations is often challenged due to limited understanding of physical interaction mechanisms at the interfaces involved. Here we model a physisorbed interaction between particles and gold electrode surface of a quartz crystal and demonstrate how the nonlinear modulation of the electric response of the crystal due to the nonlinear interaction forces may be used to study the dynamics of the particles. In particular, we show that the graphs of the deviation in the third Fourier harmonic response versus oscillation amplitude provide important information about the onset, progress and nature of sliding of the particles. The graphs also present a signature of the surface-particle interaction and could be used to estimate the interaction energy profile. Interestingly, the insights gained from the model help to explain some of the experimental observations with physisorbed streptavidin-coated polystyrene microbeads on quartz resonators. © 2012 Elsevier B.V. All rights reserved.

Nonlinear phenomena in thermoacoustics: A comparison between single-mode and multi-mode methods

Nonlinear analysis of thermoacoustic instability is essential for prediction of frequencies and amplitudes of limit cycles. In frequency domain analyses, a quasi-linear transfer function between acoustic velocity and heat release rate perturbations, called the flame describing function (FDF), is obtained from a flame model or experiments. The FDF is a function of the frequency and amplitude of velocity perturbations but only contains the heat release response at the forcing frequency. While the gain and phase of the FDF provide insight into the nonlinear dynamics of the system, the accuracy of its predictions remains to be verified for different types of nonlinearity. In time domain analyses, the governing equations of the fully coupled problem are solved to find the time evolution of the system. One method is to discretize the governing equations using a suitable basis, such as the natural acoustic modes of the system. The number of modes used in the discretization alters the accuracy of the solution. In our previous work we have shown that predictions using the FDF are almost exactly the same as those obtained from the time-domain using only one mode for the discretization. We call this the single-mode method. In this paper we compare results from the single-mode and multi-mode methods, applied to a thermoacoustic system of a premixed flame in a tube. For some cases, the results differ greatly in both amplitude as well as frequency content. This study shows that the contribution from higher and subharmonics to the nonlinear dynamics can be significant and must be considered for an accurate and comprehensive analysis of thermoacoustic systems. Hence multi-mode simulations are necessary, and the single-mode method or the FDF may be insufficient to capture some of the complex nonlinear behaviour in fhermoacoustics.

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.

Wave propagation due to sinusoidal excitation in nonlinear one-dimensional soil column

The objective of the author's on-going research is to explore the feasibility of determining reliable in situ curves of shear modulus as a function of strain using the dynamic test. The purpose of this paper is limited to investigating what material stiffness is measured from a dynamic test, focusing on the harmonic excitation test. A one-dimensional discrete model with nonlinear material properties is used for this purpose. When a sinusoidal load is applied, the cross-correlation of signals from different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loops under steady-state conditions. The variables that contributed to changing the average slope of the stress-strain loop also influence the estimate of the wave velocity from cross-correlation. Copyright ASCE 2007.