Modulation of a methane Bunsen flame by upstream perturbations

In this paper the effects of an upstream spatially periodic modulation acting on a turbulent Bunsen flame are investigated using direct numerical simulations of the Navier-Stokes equations coupled with the flamelet generated manifold (FGM) method to parameterise the chemistry. The premixed Bunsen flame is spatially agitated with a set of coherent large-scale structures of specific wave-number, K. The response of the premixed flame to the external modulation is characterised in terms of time-averaged properties, e.g. the average flame height ?H? and the flame surface wrinkling ?W?. Results show that the flame response is notably selective to the size of the length scales used for agitation. For example, both flame quantities ?H? and ?W? present an optimal response, in comparison with an unmodulated flame, when the modulation scale is set to relatively low wave-numbers, 4π/L ? K ? 6π/L, where L is a characteristic scale. At the agitation scales where the optimal response is observed, the average flame height, ?H?, takes a clearly defined minimal value while the surface wrinkling, ?W?, presents an increase by more than a factor of 2 in comparison with the unmodulated reference case. Combined, these two response quantities indicate that there is an optimal scale for flame agitation and intensification of combustion rates in turbulent Bunsen flames.     

Multi-scale interactions in a compressible boundary layer

The properties of spectral subranges of scales in a boundary layer at Mach=2.3 and friction Reynolds number Reτ = 570 are investigated by analysing DNS data. One major aim is to examine whether footprinting and modulation of small-scale near-wall motions by outer large structures, observed at high Reynolds numbers, also pertain to this low-Reynolds-number case, or whether the logarithmic layer simply contains a continuous hierarchy of motions without specific outer scales playing a distinctive role. To this end, the spectrum of scales is decomposed into modes by application of the “Empirical Mode Decomposition”. The properties of different scales are then investigated by means of spectra, maps of isotropy/anisotropy parameters, the premultiplied derivative of the second-order structure function, correlation coefficients and joint probability density function (PDF), the last constructed from conditionally sampled data for the small-scale motions within the large-scale footprints. A clear commonality is identified between interactions in high-Reynolds-number channel flow and the present low-Reynolds-number boundary layer.     

Dynamos at extreme magnetic Prandtl numbers: insights from shell models

We present a magnetohydrodynamic (MHD) shell model suitable for computation of various energy fluxes of MHD turbulence for very small and very large magnetic Prandtl numbers Pm; such computations are inaccessible to direct numerical simulations. For small Pm, we observe that both kinetic and magnetic energy spectra scale as k^?5/3 in the inertial range, but the dissipative magnetic energy scales as k^?11/3exp?(? k/k_η). Here the kinetic energy at large length scale feeds the large-scale magnetic field that cascades to small-scale magnetic field, which gets dissipated by Joule heating. The large-Pm dynamo has a similar behaviour except that the dissipative kinetic energy scales as k^?13/3. For this case, the large-scale velocity field transfers energy to the large-scale magnetic field, which gets transferred to small-scale velocity and magnetic fields; the energy of the small-scale magnetic field also gets transferred to the small-scale velocity field, and the energy thus accumulated is dissipated by the viscous force.     

Anisotropic Plasmonic Metasurfaces with Linear Polarization-Dependent Focusings

As a burgeoning approach to modulate electromagnetic waves, metasurfaces fascinate many researchers. However, it is a persistent problem to realize independent wavefront tailoring in different polarization states. In this work, anisotropic plasmonic metasurfaces are presented with linear polarization-dependent focusings. The designed plasmonic meta-atoms consist of cross-shaped gold, silica spacer, and gold substrate. Phase modulation in different polarization channels is independently achieved by adjusting the dimensions of cross-shaped gold in x- and y-directions. Three metasurfaces are presented to verify linear polarization-dependent focusings at 200 THz. First, a focusing metasurface is designed with focal lengths of 5.2 µm under x-polarized incidence and 6.5 µm under y-polarized incidence. Second, by applying convolution operation, the second metasurface exhibits different focusings with different deflection angles at opposite directions under orthogonally polarized incidences. Finally, a multi-focal metasurface is demonstrated, which switches between dual- and quad-focal points depending on polarization state. The work may provide a novel platform for near-infrared integrated photonics.     

The effect of the modulation of Bragg scattering in small-slope approximation

Abstract

Small-slope approximation (SSA) belongs to a class of ‘unifying’ scattering theories which reproduce small perturbations and semiclassic (Kirchhoff) results within appropriate limits. However, the most stringent test for such theories involves a two-scale situation when a small-scale roughness is located on a tilted plane. A ‘unifying’ theory should properly account for the effects of modulation of the scattering cross section associated with a large-scale tilt. This paper shows that SSA does properly take into account these modulation effects.     

Lees–Edwards Boundary Conditions for Lattice Boltzmann

Lees–Edwards boundary conditions (LEbc) for Molecular Dynamics simulations⁽¹⁾ are an extension of the well known periodic boundary conditions and allow the simulation of bulk systems in a simple shear flow. We show how the idea of LEbc can be implemented in isothermal lattice Boltzmann simulations and how LEbc can be used to overcome the problem of a maximum shear rate that is limited to less then 1/L _y (with L _y the transverse system size) in traditional lattice Boltzmann implementations of shear flow. The only previous Lattice Boltzmann implementation of LEbc⁽²⁾ requires a specific fourth order equilibrium distribution. In this paper we show how LEbc can be implemented with the usual quadratic equilibrium distributions.     

Observation of polarization-controlled spatial splitting of four-wave mixing in a three-level atomic system

We report experimental observations of intensity modulation and spatial splitting of four-wave mixing (FWM) signal beams which can be effectively controlled by the polarization states of the pumping laser beams. Due to the periodic change of the pumping beam’s polarization states, the intensity of the FWM beam also evolves periodically. The periodic spatial splitting phenomenon has been observed in both x- and y-directions. The cases with/without the dressing beams are compared. Such studies can be very useful in better understanding the formation of spatial solitons and for signal processing applications, such as spatial beam splitter, routing, and switching.     

Absolute cross sections for collision-induced depolarized scattering of light in krypton and xenon

Depolarized Raman spectra of binary collisional pairs of atoms in krypton and xenon are obtained at gas densities of 1–10 amagat. Absolute intensities relative to a known rotational transition of nitrogen are determined. For light of 4880 Å wavelength incident in the x-direction, polarized in the z-direction and scattered in the y-direction of a cartesian frame x, y, z, at a frequency shift of -12 cm^-1, the differential scattering cross section per unit wavenumber band times volume, is found to be 1·10 × 10^-52 cm⁶ ± 10 per cent for krypton, and 4·76 × 10^-52 cm⁶ ± 10 per cent for xenon, if the sum of both polarizations is considered. Wave-mechanical and classical computations reproduce both the shape and the intensity of the experimental spectra if the so-called point-dipole model of the anisotropy of the polarizability of collisional pairs of atoms is used. Other models of the anisotropy are seen to be overcorrected by these criteria.     

Mass Power Spectrum in a Universe Dominated by the Chaplygin Gas

The mass power spectrum for a Universe dominated by the Chaplygin gas is evaluated numerically from scales of the order of the Hubble horizon to 100 Mpc. The results are compared with a pure baryonic Universe and a cosmological constant model. In all three cases, the spectrum increases with k, the wavenumber of the perturbations. The slope of the spectrum is higher for the baryonic model and smaller for the cosmological constant model, the Chaplygin gas interpolating these two models. The results are analyzed in terms of the sound velocity of the Chaplygin gas and the moment the Universe begins to accelerate.     

Inner–outer interactions in rough-wall turbulence

Here we revisit the inner–outer interaction model (IOIM) of Marusic et al. (Science, vol. 329, 2010, pp. 193–196) that enables the prediction of statistics of the fluctuating streamwise velocity in the inner region of wall-bounded turbulent flows from a large-scale velocity signature measured in the outer region of the flow. The model is characterised by two empirically observed inner–outer interactions: superposition of energy from outer region large-scale motions; and amplitude modulation by these large-scale motions of a small-scale ‘universal’ signal (u^*), which in smooth-wall flows is Reynolds number invariant. In the present study, the inner–outer interactions in rough-wall turbulent boundary layers are examined within the framework of the IOIM. Simultaneous two-point hot-wire anemometry measurements enable quantification, via the model parameters, of the strengths of superposition and amplitude modulation effects in a rough-wall flow, and these are compared to a smooth-wall flow. It is shown that the present rough-wall significantly reduces the effects of superposition, while increasing the amplitude modulation effect. The former is true even in flows that exhibit outer region similarity. Using the model parameters obtained from the two-point measurements, predictions of inner region streamwise velocity statistics and spectra are compared to measurements over a range of friction and roughness Reynolds numbers. These results indicate that the u^* signal does depend on roughness Reynolds number (k⁺_s), but is robust to changes in friction Reynolds number (δ⁺). Additionally, the superposition strength is shown to be relatively independent of both roughness and friction Reynolds number. The implications of the present results on the suitability of the IOIM as a predictive tool in rough-wall turbulence are discussed.     

Small-scale-velocity-induced linear instability of large-scale turbulent flows

Analysis of a simplified equation derived previously for small-scale velocity components shows that any turbulent flow of an incompressible liquid becomes unstable against infinitesimal perturbations of small-scale velocity components if the strain rate tensor for the large-scale velocity is high. Such a statement comes into conflict with the classical stability theory, which specifically asserts that the Poiseuille flow in a circular tube is linearly stable against infinitesimal perturbations.     

Visualization and characterization of small-scale spanwise vortices in turbulent channel flow

High-resolution particle-image velocimetry (PIV) measurements are made in the streamwise-wall-normal plane of turbulent channel flow at Re_τ=566, 1184 and 1759, facilitating documentation of the population trends and core diameters of small-scale spanwise vortices. Swirling strength, an unambiguous vortex-identification criterion and hence a local marker of rotation, is used to extract small-scale spanwise vortex cores from the instantaneous velocity fields. Once the small-scale vortices are properly extracted from the PIV realizations, their characteristics are studied in detail. The present results indicate that the very-near-wall region (y < 0.1h) is densely populated by spanwise vortices with clockwise (negative) rotation. This behavior supports the notion that hairpin-like vortices are generated very close to the wall and grow into the outer layer as they advect downstream. In contrast, counterclockwise (positive) spanwise vortices are scarce in the very-near-wall region, but their presence steadily increases within the logarithmic layer presumably due to a localized generation mechanism. The average core diameter of negative spanwise vortices is found to be larger than the average diameter of positive vortices, with few positive vortices having core diameters exceeding 80y.     

Modulation instability of electromagnetic excitations for a Josephson unction in a finite-thickness slab

The modulation instability of finite-amplitude uniform plane waves oscillating with a Josephson frequency and experiencing a nonlinear frequency shift in a finite-thickness slab is studied in terms of the nonlocal electrodynamics of Josephson junction. A dispersion relation for the growth rate of small amplitude perturbation is derived. The domains of modulation instability for these waves are found. Modulation instability of the waves is shown to arise when the wavevectors of long-wave amplitude perturbations fall into the finite range 0 < Q < Q _B (A, D, L). In the range Q ≥ Q _B (A, D, L), the waves are stable.     

Vorticity,backscatter and counter-gradient transport predictions using two-level simulation of turbulent flows

The two-level simulation (TLS) method evolves both the large-and the small-scale fields in a two-scale approach and has shown good predictive capabilities in both isotropic and wall-bounded high Reynolds number (Re) turbulent flows in the past. Sensitivity and ability of this modelling approach to predict fundamental features (such as backscatter, counter-gradient turbulent transport, small-scale vorticity, etc.) seen in high Re turbulent flows is assessed here by using two direct numerical simulation (DNS) datasets corresponding to a forced isotropic turbulence at Taylor’s microscale-based Reynolds number Re_λ ≈ 433 and a fully developed turbulent flow in a periodic channel at friction Reynolds number Re_τ ≈ 1000. It is shown that TLS captures the dynamics of local co-/counter-gradient transport and backscatter at the requisite scales of interest. These observations are further confirmed through a posteriori investigation of the flow in a periodic channel at Re_τ = 2000. The results reveal that the TLS method can capture both the large- and the small-scale flow physics in a consistent manner, and at a reduced overall cost when compared to the estimated DNS or wall-resolved LES cost.     

Flow-frequency effect upon Huygens front propagation

A premixed flame within a turbulent flow exhibits a decreasing enhancement of fuel consumption rate with increasing turbulence intensity, an effect known as the bending effect. Denet has shown that flow time correlations may be one cause of the bending effect. Using a Damköhler-Huygens front propagation model, we illustrate that the removal of flow components with reduced frequencies greater than unity (ω >kS _L) causes a small reduction in front area but a large reduction in the flow intensity, which is the bending effect (ω is the frequency and k is the wavenumber). To be effective in producing front area, a flow mode must have a phase velocity, ω/k, smaller than the laminar burning velocity, S _L.     

Experimental Investigation of Mixed Convection in a Channel With an Open Cavity

Mixed convection in an open cavity with a heated wall bounded by a horizontally unheated plate is investigated experimentally. The cavity has the heated wall on the inflow side. Mixed convection fluid flow and heat transfer within the cavity is governed by the buoyancy parameter, Richardson number (Ri), and Reynolds number (Re). The results are reported in terms of wall temperature profiles of the heated wall and flow visualization for Re = 100 and 1000, Ri in the range 30–110 (for Re = 1000) and 2800–8700 (for Re = 100), the ratio of the length to the height of cavity (L/D) is in the range 0.5–1.5, and the ratio of the channel height to cavity height (H/D) is in the range of 0.5 and 1.0. The present results show that the maximum dimensional temperature rise values decrease as the Reynolds and the Richardson numbers decrease. The flow visualization points out that for Re = 1000 there are two nearly distinct fluid motions: a parallel forced flow in the channel and a recirculation flow inside the cavity. For Re = 100 the effect of a stronger buoyancy determines a penetration of thermal plume from the heated plate wall into the upper channel. Nusselt numbers increase when L/D increase in the considered range of Richardson numbers.     

Wavelet multi-resolution analysis of initial condition effects on near-wake turbulent structures

The effects of initial conditions on turbulence structures of various scales in a near wake have been investigated for two wake generators with the same characteristic dimension, i.e., a circular cylinder and a screen of 50% solidity, based on the wavelet multi-resolution analysis. The experimental investigation used two orthogonal arrays of sixteen X-wires, eight in the (x, y)-plane, and eight in the (x, z)-plane. Measurements were made atx/h (x is the streamwise distance downstream of the cylinder andh is the height of the wake generator) = 20. The wavelet multi-resolution technique was applied to decomposing the velocity data, obtained in the wakes generated by the two generators, into a number of wavelet components based on the central frequencies. The instantaneous sectional streamlines and vorticity field were thus ‘visualized’ for each wavelet component or central frequency. It was found that the behavior of large- and intermediate-scale structures depend on the initial conditions and the small-scale structures are independent of the initial conditions. The contributions from the wavelet components to the time-averaged Reynolds stresses and vorticity were estimated. Both the large-scale and intermediate longitudinal structures make the most significant contributions to Reynolds stresses in the circular cylinder wake, but the contribution from the large-scale structures appears dominating in the screen wake. The relatively small scale structures of the circular cylinder wake contribute most to the total rms spanwise vorticity.     

采用密集波分复用技术的光纤水听器阵列研究

利用密集波分复用和时分多路复用技术相结合的大规模阵列结构,以Mach-Zehnder干涉型光纤水听器为例,分析了采用相位产生载波技术的频分多路复用,提出了密集波分复用技术在干涉型光纤水听器阵列应用的新方法,给出了复用体系结构,并分析了其在工程上可行性.     

Experimental study of initial condition dependence on turbulent mixing in shock-accelerated Richtmyer–Meshkov fluid layers

Experimental evidence is needed to verify the hypothesis that the memory of initial conditions is retained at late times in variable density flows. If true, this presents an opportunity to “design” and “control” late-time turbulence, with an improved understanding in the prediction of inertial confinement fusion and other general fluid mixing processes. In this communication, an experimental and theoretical study on the effects of initial condition parameters, namely, the amplitude δ₀ and wavenumber κ₀ , where λ₀ is the initial wavelength) of perturbations, on late-time turbulence and mixing in shock-driven Richtmyer–Meshkov (R-M) unstable fluid layers in a 2D plane is presented. Single and multi-mode membrane-free initial conditions in the form of a gas curtain having a light-heavy-light configuration (air-SF₆-air) with an Atwood number of A= 0.57 were used in our experiments. A planar shock wave with a shock Mach number M= 1.21 drives the R-M instability, and the evolution of this instability after incident shock is captured using high resolution simultaneous planar laser induced fluorescence (PLIF) and particle image velocimetry (PIV) diagnostics. Time evolution of statistics such as amplitude of the mixing layer, 2D turbulent kinetic energy, Reynolds number, rms of velocity fluctuations, probability density functions, and density-specific volume correlation were observed to quantify the amount of mixing and understand the nature of turbulence in this flow. Based on these results, it was found that the R-M mixing layer is asymmetric and non-Boussinesq. There is a correlation between initial condition parameters and large-scale, and small-scale mixing at late times, indicating an initial condition dependence on R-M mixing.