On the distortion of turbulence by a progressive surface wave

A rapid-distortion model is developed to investigate the interaction of weak turbulence with a monochromatic irrotational surface water wave. The model is applicable when the orbital velocity of the wave is larger than the turbulence intensity, and when the slope of the wave is sufficiently high that the straining of the turbulence by the wave dominates over the straining of the turbulence by itself. The turbulence suffers two distortions. Firstly, vorticity in the turbulence is modulated by the wave orbital motions, which leads to the streamwise Reynolds stress attaining maxima at the wave crests and minima at the wave troughs; the Reynolds stress normal to the free surface develops minima at the wave crests and maxima at the troughs. Secondly, over several wave cycles the Stokes drift associated with the wave tilts vertical vorticity into the horizontal direction, subsequently stretching it into elongated streamwise vortices, which come to dominate the flow. These results are shown to be strikingly different from turbulence distorted by a mean shear flow, when `streaky structures' of high and low streamwise velocity fluctuations develop. It is shown that, in the case of distortion by a mean shear flow, the tendency for the mean shear to produce streamwise vortices by distortion of the turbulent vorticity is largely cancelled by a distortion of the mean vorticity by the turbulent fluctuations. This latter process is absent in distortion by Stokes drift, since there is then no mean vorticity. The components of the Reynolds stress and the integral length scales computed from turbulence distorted by Stokes drift show the same behaviour as in the simulations of Langmuir turbulence reported by McWilliams, Sullivan & Moeng (1997). Hence we suggest that turbulent vorticity in the upper ocean, such as produced by breaking waves, may help to provide the initial seeds for Langmuir circulations, thereby complementing the shear-flow instability mechanism developed by Craik & Leibovich (1976). The tilting of the vertical vorticity into the horizontal by the Stokes drift tends also to produce a shear stress that does work against the mean straining associated with the wave orbital motions. The turbulent kinetic energy then increases at the expense of energy in the wave. Hence the wave decays. An expression for the wave attenuation rate is obtained by scaling the equation for the wave energy, and is found to be broadly consistent with available laboratory data.

Referencing of street-level flows measured during the DAPPLE 2004 campaign

Street-level mean flow and turbulence govern the dispersion of gases away from their sources in urban areas. A suitable reference measurement in the driving flow above the urban canopy is needed to both understand and model complex street-level flow for pollutant dispersion or emergency response purposes. In vegetation canopies, a reference at mean canopy height is often used, but it is unclear whether this is suitable for urban canopies. This paper presents an evaluation of the quality of reference measurements at both roof-top (height = H) and at height z = 9H = 190 m, and their ability to explain mean and turbulent variations of street-level flow. Fast response wind data were measured at street canyon and reference sites during the six-week long DAPPLE project field campaign in spring 2004, in central London, UK, and an averaging time of 10 min was used to distinguish recirculation-type mean flow patterns from turbulence. Flow distortion at each reference site was assessed by considering turbulence intensity and streamline deflection. Then each reference was used as the dependent variable in the model of Dobre et al. (2005) which decomposes street-level flow into channelling and recirculating components. The high reference explained more of the variability of the mean flow. Coupling of turbulent kinetic energy was also stronger between street-level and the high reference flow rather than the roof-top. This coupling was weaker when overnight flow was stratified, and turbulence was suppressed at the high reference site. However, such events were rare (<1% of data) over the six-week long period. The potential usefulness of a centralised, high reference site in London was thus demonstrated with application to emergency response and air quality modelling.

Scalability of efficient parallel K-Means

Clustering is defined as the grouping of similar items in a set, and is an important process within the field of data mining. As the amount of data for various applications continues to increase, in terms of its size and dimensionality, it is necessary to have efficient clustering methods. A popular clustering algorithm is K-Means, which adopts a greedy approach to produce a set of K-clusters with associated centres of mass, and uses a squared error distortion measure to determine convergence. Methods for improving the efficiency of K-Means have been largely explored in two main directions. The amount of computation can be significantly reduced by adopting a more efficient data structure, notably a multi-dimensional binary search tree (KD-Tree) to store either centroids or data points. A second direction is parallel processing, where data and computation loads are distributed over many processing nodes. However, little work has been done to provide a parallel formulation of the efficient sequential techniques based on KD-Trees. Such approaches are expected to have an irregular distribution of computation load and can suffer from load imbalance. This issue has so far limited the adoption of these efficient K-Means techniques in parallel computational environments. In this work, we provide a parallel formulation for the KD-Tree based K-Means algorithm and address its load balancing issues.

On the initiation of surface waves by turbulent shear flow

An analytical model is developed for the initial stage of surface wave generation at an air-water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417-445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate F, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer. For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as proportional to 1/Gamma), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.

The calculation of accurate normal coordinates I: general theory, and application to methane

The calculation of accurate and reliable vibrational potential functions and normal co-ordinates is discussed, for such simple polyatomic molecules as it may be possible. Such calculations should be corrected for the effects of anharmonicity and of resonance interactions between the vibrational states, and should be fitted to all the available information on all isotopic species: particularly the vibrational frequencies, Coriolis zeta constants and centrifugal distortion constants. The difficulties of making these corrections, and of making use of the observed data are reviewed. A programme for the Ferranti Mercury Computer is described by means of which harmonic vibration frequencies and normal co-ordinate vectors, zeta factors and centrifugal distortion constants can be calculated, from a given force field and from given G-matrix elements, etc. The programme has been used on up to 5 × 5 secular equations for which a single calculation and output of results takes approximately l min; it can readily be extended to larger determinants. The best methods of using such a programme and the possibility of reversing the direction of calculation are discussed. The methods are applied to calculating the best possible vibrational potential function for the methane molecule, making use of all the observed data.

Microwave spectrum of CHD2CN and CHD2NC, and the harmonic force field and rz structure of methyl cyanide and isocyanide

The microwave spectra of CHD2CN and CHD2NC have been measured from 18 to 40 GHz; about 20 type A and 30 type C transitions have been observed for each molecule. These have been fitted to a Hamiltonian using 3 rotational constants, and 5 quartic and 4 sextic distortion constants, in the IrS reduction of Watson [in “Vibrational spectra and structure” Vol. 6 (1977)]; the standard error of the fit is 26 kHz. For methyl cyanide the 5 quartic distortion constants have been used to further refine the recent harmonic force field of Duncan et al. [J. Mol. Spectrosc. 69, 123 (1978)], but the changes are small. Finally, for both molecules, the harmonic force field has been used to determine zero point average moments of inertia Iz from the ground state rotational constants for many isotopic species, and these have been used to determine an rz structure. The results are compared with rs structure calculations.

Fourier transform infrared spectra of H2CS and D2CS

Infrared spectra of thoformaldehyde, H2CS and D2CS, were observed in the gas phase at a resolution of better than 0.1 cm−1 from 4000 to 400 cm−1 using a Nicolet FTIR system. Vibrational band origins and rotational constants were determined for ν2, ν3, ν4, and ν6 of H2CS and for ν1, ν2, ν3, ν4, and ν6 of D2CS. The ν3, ν4, and ν6 bands of H2CS were analyzed as a set of three Coriolis interacting bands, and three Coriolis constants were determined; similarly the ν4 and ν6 bands of D2CS were analyzed as a pair of interacting bands and one Coriolis constant was determined. A general harmonic force field was determined, without constraints, to fit the vibrational wavenumbers, Coriolis constants, and centrifugal distortion constants. A zero-point (rz) structure was determined from the ground-state rotational constants, and the equilibrium (re) bond lengths were estimated.

Excited vibrational state microwave spectra, structure, and force-field of trifluorosilane

The J = 2−1 microwave spectrum of six isotopic species of HSiF3 has been observed and assigned in excited states of five of the six fundamental vibrations. The assignment is based on relative intensities, double resonance experiments, and trial anharmonic force constant calculations. Analysis of the spectra leads to experimental values for five of the constants, all three l-doubling constants qt, one Fermi resonance constant φ233, and one zeta constant. The harmonic force field has been refined to all the available data on vibration wavenumbers, centrifugal distortion constants, and zeta constants. The cubic anharmonic force field has been refined to the data on and qt constants, using two models: a valence force model with two cubic force constants for SiH and SiF stretching, and a more sophisticated model. With the help of these calculations, the following equilibrium structure has been determined: re(SiH) = 1.4468(±5) Å, re(SiF) = 1.5624(±1) Å, HSiF = 110.64(±3)°,

Vibration-rotation spectra and the harmonic force field of HOCl

Vibration-rotation spectra of HOCl have been measured at a resolution of 0.05 cm−1 to determine vibration rotation constants, and 35–37 Cl isotope shifts in the vibration frequencies. The spectrum of DOCl has also been recorded, and a preliminary analysis for the band origins has been made. The vibrational frequency data and centrifugal distortion constants have been used to determine the harmonic force field in a least-squares refinement; the force field obtained also gives a good fit to data on the vibrational contributions to the inertial defect. The equilibrium rotational constants of HOCl have been obtained, and an equilibrium structure has been estimated.

The harmonic force field and rz structure of HNCO

The presently available microwave, millimeter wave, and far-infrared data of five isotopic species of isocyanic acid, namely, HNCO, H15NCO, HN13CO, HNC18O, and DNCO, have been used to obtain improved values of the ground-state rotational constants, the five quartic distortion constants, and some higher-order distortion constants in the Ir S reduced Hamiltonian of Watson. The appropriate planarity relation among the quartic centrifugal distortion constants has been imposed in the fitting procedure. The general harmonic force field of isocyanic acid has been determined using all existing data, and assuming a trans bent equilibrium geometry of the molecule with an NCO angle of 170°. Finally an rz structure has been obtained using the Az, Bz, and Cz rotational constants of five isotopic species. The bending of the NCO chain is found to be 8° in the trans configuration.

The microwave spectrum of thietan-2,2,4,4-d4

The microwave spectrum for thietan-2,2,4,4-d4 is analysed in six of its lowest puckering states and up to J = 25. The close lying pairs of states of vp = 0/1, 2/3 and 4/5 are treated with a vibration-rotation hamiltonian which includes an off-diagonal coupling term in vp. Additional corrections to this coupling term in higher powers of the angular momentum operator are derived and their importance for improving the fit of calculated to observed data is tested. The variation of the centrifugal distortion constants with vp follows the model of Creswell and Mills (1974, J. molec. Spectrosc., 52, 392). A value is determined for the derivative with respect to the puckering coordinate of the ac-component of the inverse moment of inertia tensor.

The microwave spectrum of 1-pyrazoline

The microwave spectrum of 1-pyrazoline has been observed from 18 to 40 GHz in the six lowest states of the ring-puckering vibration. It is an a-type spectrum of a near oblate asymmetric top. Each vibrational state has been fitted to a separate effective Hamiltonian, and the vibrational dependence of both the rotational constants and the quartic centrifugal distortion constants has been observed and analyzed. The v = 0 and 1 states have also been analyzed using a coupled Hamiltonian; this gives consistent results, with an improved fit to the high J data. The preferred choice of Durig et al. [J. Chem. Phys. 52, 6096 (1970)] for the ring-puckering potential is confirmed as essentially correct, but the A and B inertial axes are shown to be interchanged from those assumed by Durig et al. in their analysis of the mid-infrared spectrum.

Microwave spectrum of deuterated silyl isocyanate, SiD3NCO

The microwave spectrum of SiD3NCO has been observed and analyzed for 18 different vibrational states in the ν10 manifold. Some accidental resonances have been observed and analyzed. The vibrational dependence of the rotational and l-doubling constant and centrifugal distortion constant DJK has been successfully interpreted in terms of the two-dimensional anharmonic oscillator model.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.