Wave-activity conservation laws for the three-dimensional anelastic and Boussinesq equations with a horizontally homogeneous background flow

Wave-activity conservation laws are key to understanding wave propagation in inhomogeneous environments. Their most general formulation follows from the Hamiltonian structure of geophysical fluid dynamics. For large-scale atmospheric dynamics, the Eliassen–Palm wave activity is a well-known example and is central to theoretical analysis. On the mesoscale, while such conservation laws have been worked out in two dimensions, their application to a horizontally homogeneous background flow in three dimensions fails because of a degeneracy created by the absence of a background potential vorticity gradient. Earlier three-dimensional results based on linear WKB theory considered only Doppler-shifted gravity waves, not waves in a stratified shear flow. Consideration of a background flow depending only on altitude is motivated by the parameterization of subgrid-scales in climate models where there is an imposed separation of horizontal length and time scales, but vertical coupling within each column. Here we show how this degeneracy can be overcome and wave-activity conservation laws derived for three-dimensional disturbances to a horizontally homogeneous background flow. Explicit expressions for pseudoenergy and pseudomomentum in the anelastic and Boussinesq models are derived, and it is shown how the previously derived relations for the two-dimensional problem can be treated as a limiting case of the three-dimensional problem. The results also generalize earlier three-dimensional results in that there is no slowly varying WKB-type requirement on the background flow, and the results are extendable to finite amplitude. The relationship A E =cA P between pseudoenergy A E and pseudomomentum A P, where c is the horizontal phase speed in the direction of symmetry associated with A P, has important applications to gravity-wave parameterization and provides a generalized statement of the first Eliassen–Palm theorem.

A Hamiltonian weak-wave model for shallow-water flow

A reduced dynamical model is derived which describes the interaction of weak inertia–gravity waves with nonlinear vortical motion in the context of rotating shallow–water flow. The formal scaling assumptions are (i) that there is a separation in timescales between the vortical motion and the inertia–gravity waves, and (ii) that the divergence is weak compared to the vorticity. The model is Hamiltonian, and possesses conservation laws analogous to those in the shallow–water equations. Unlike the shallow–water equations, the energy invariant is quadratic. Nonlinear stability theorems are derived for this system, and its linear eigenvalue properties are investigated in the context of some simple basic flows.

On the generalized eigenvalue problem for the Rossby wave vertical velocity in the presence of mean flow and topography

In a series of papers, Killworth and Blundell have proposed to study the effects of a background mean flow and topography on Rossby wave propagation by means of a generalized eigenvalue problem formulated in terms of the vertical velocity, obtained from a linearization of the primitive equations of motion. However, it has been known for a number of years that this eigenvalue problem contains an error, which Killworth was prevented from correcting himself by his unfortunate passing and whose correction is therefore taken up in this note. Here, the author shows in the context of quasigeostrophic (QG) theory that the error can ulti- mately be traced to the fact that the eigenvalue problem for the vertical velocity is fundamentally a non- linear one (the eigenvalue appears both in the numerator and denominator), unlike that for the pressure. The reason that this nonlinear term is lacking in the Killworth and Blundell theory comes from neglecting the depth dependence of a depth-dependent term. This nonlinear term is shown on idealized examples to alter significantly the Rossby wave dispersion relation in the high-wavenumber regime but is otherwise irrelevant in the long-wave limit, in which case the eigenvalue problems for the vertical velocity and pressure are both linear. In the general dispersive case, however, one should first solve the generalized eigenvalue problem for the pressure vertical structure and, if needed, diagnose the vertical velocity vertical structure from the latter.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

The physics of orographic gravity wave drag

The drag and momentum fluxes produced by gravity waves generated in flow over orography are reviewed, focusing on adiabatic conditions without phase transitions or radiation effects, and steady mean incoming flow. The orographic gravity wave drag is first introduced in its simplest possible form, for inviscid, linearized, non-rotating flow with the Boussinesq and hydrostatic approximations, and constant wind and static stability. Subsequently, the contributions made by previous authors (primarily using theory and numerical simulations) to elucidate how the drag is affected by additional physical processes are surveyed. These include the effect of orography anisotropy, vertical wind shear, total and partial critical levels, vertical wave reflection and resonance, non-hydrostatic effects and trapped lee waves, rotation and nonlinearity. Frictional and boundary layer effects are also briefly mentioned. A better understanding of all of these aspects is important for guiding the improvement of drag parametrization schemes.

Impact of non-hydrostatic effects and trapped lee waves on mountain wave drag in directionally sheared flow

The orographic gravity wave drag produced in flow over an axisymmetric mountain when both vertical wind shear and non-hydrostatic effects are important was calculated using a semi-analytical two-layer linear model, including unidirectional or directional constant wind shear in a layer near the surface, above which the wind is constant. The drag behaviour is determined by partial wave reflection at the shear discontinuity, wave absorption at critical levels (both of which exist in hydrostatic flow), and total wave reflection at levels where the waves become evanescent (an intrinsically non-hydrostatic effect), which produces resonant trapped lee wave modes. As a result of constructive or destructive wave interference, the drag oscillates with the thickness of the constant-shear layer and the Richardson number within it (Ri), generally decreasing at low Ri and when the flow is strongly non-hydrostatic. Critical level absorption, which increases with the angle spanned by the wind velocity in the constant-shear layer, shields the surface from reflected waves, keeping the drag closer to its hydrostatic limit. While, for the parameter range considered here, the drag seldom exceeds this limit, a substantial drag fraction may be produced by trapped lee waves, particularly when the flow is strongly non-hydrostatic, the lower layer is thick and Ri is relatively high. In directionally sheared flows with Ri = O(1), the drag may be misaligned with the surface wind in a direction opposite to the shear, a behaviour which is totally due to non-trapped waves. The trapped lee wave drag, whose reaction force on the atmosphere is felt at low levels, may therefore have a distinctly different direction from the drag associated with vertically propagating waves, which acts on the atmosphere at higher levels.

Flow in the high latitude ionosphere: measurements at 15s resolution made using the EISCAT ‘Polar’ experiment

We present a first overview of flows in the high latitude ionosphere observed at 15 s resolution using the U.K.-Polar EISCAT experiment. Data are described from experiments conducted on two days, 27 October 1984 and 29 August 1985, which together span the local times between about 0200 and 2130MLT and cover five different regions of ionospheric flow. With increasing local time, these are: the dawn auroral zone flow cell, the dayside region of low background flows equatorward of the flow cells, the dusk auroral zone flow cell, the boundary region between the dusk auroral zone and the polar cap, and the evening polar cap. Flows in both the equatorward and poleward portions of the auroral zone cells appear to be relatively smooth, while in the central region of high speed flow considerable variations are generally present. These have the form of irregular fluctuations on a wide range of time scales in the early morning dawn cell, and impulsive wave-like variations with periods of a few minutes in the afternoon dusk cell. In the dayside region between the flow cells, the ionosphere is often essentially stagnant for long intervals, but low amplitude ULF waves with a period of about 5 min can also occur and persist for many cycles. These conditions are punctuated at one to two hour intervals by sudden ‘flow burst’ events with impulsively generated damped wave trains. Initial burst flows are generally directed poleward and can peak at line-of-sight speeds in excess of 1 km s^{−1} after perhaps 45 s. Flows in the polar cap are reasonably smooth on time scales of a few minutes and show no evidence for the presence of ULF waves. Under most, but not all, of the above conditions, the beam-swinging algorithm used to determine background vector flows should produce meaningful results. Comparison of these flow data with simultaneous plasma and magnetic field measurements in the solar wind, made by the AMPTE IRM and UKS spacecraft, emphasizes the strong control exerted on high latitude flows by the north-south component of the IMF.

Analysis of acoustic frequency and wave amplitude generated by the Oscillatory Thoracic Thixotropic Device (Diottix(r)) in human chest

Diottix(r) was calibrated at 25 Hz to achieve the frequency indicated in literature as being effective to mobilize the airways secretions. However, the amplitude and frequency of the waves generated by the equipment in different regions of the chest still need to be investigated. The objective of this study was to analyze the frequency and amplitude of waves generated by Diottix(r) in chests of healthy subjects. Diottix(r) was used in the anterior and posterior regions of the chest. The mechanical waves were captured using stethoscopes connected to electret microphones, which were connected to a digital oscilloscope. Frequency and amplitude data were recorded by the stethoscope, positioned in six points in the anterior region and six in the posterior region of the chest, following the positions commonly used in pulmonary auscultation. Signals were recorded and transferred to a computer with software for their analysis. The frequency of waves did not present a significant change (from 24.9 to 26.4 Hz). The wave amplitude in the anterior versus the posterior region in each area of the lung, the upper, middle and lower, had differences. Diottix(r) produces frequencies in the chest according to the calibrated; thus, it can be a complementary resource to bronchial hygiene maneuvers. The amplitudes of waves seem to be affected by other structures like bone parts and heart.

On the equivalence between the wave-packet phase space representation (WPPSR) and the phase space generated by the squeezed states

We show that the parametrized Wave-Packet Phase Space representation, which has been studied earlier by one of the authors, is equivalent to a Squeezed States Phase Space Representation of quantum mechanics. © 1988.

Added mass and wave radiation damping for flow-induced rotational vibrations of skinplates of hydraulic gates

One observed vibration mode for Tainter gate skinplates involves the bending of the skinplate about a horizontal nodal line. This vibration mode can be approximated as a streamwise rotational vibration about the horizontal nodal line. Such a streamwise rotational vibration of a Tainter gate skinplate must push away water from the portion of the skinplate rotating into the reservoir and draw water toward the gate over that portion of the skinplate receding from the reservoir. The induced pressure is termed the push-and-draw pressure. In the present paper, this push-and-draw pressure is analyzed using the potential theory developed for dissipative wave radiation problems. In the initial analysis, the usual circular-arc skinplate is replaced by a vertical, flat, rigid weir plate so that theoretical calculations can be undertaken. The theoretical push-and-draw pressure is used in the derivation of the non-dimensional equation of motion of the flow-induced rotational vibrations. Non-dimensionalization of the equation of motion permits the identification of the dimensionless equivalent added mass and the wave radiation damping coefficients. Free vibration tests of a vertical, flat, rigid weir plate model, both in air and in water, were performed to measure the equivalent added mass and the wave radiation damping coefficients. Experimental results compared favorably with the theoretical predictions, thus validating the theoretical analysis of the equivalent added mass and wave radiation damping coefficients as a prediction tool for flow-induced vibrations. Subsequently, the equation of motion of an inclined circular-arc skinplate was developed by incorporating a pressure correction coefficient, which permits empirical adaptation of the results from the hydrodynamic pressure analysis of the vertical, flat, rigid weir plate. Results from in-water free vibration tests on a 1/31-scale skinplate model of the Folsom Dam Tainter gate are used to demonstrate the utility of the equivalent added mass coefficient.

From the chiral magnetic wave to the charge dependence of elliptic flow

The quark-gluon plasma formed in heavy ion collisions contains charged chiral fermions evolving in an external magnetic field. At finite density of electric charge or baryon number (resulting either from nuclear stopping or from fluctuations), the triangle anomaly induces in the plasma the Chiral Magnetic Wave (CMW). The CMW first induces a separation of the right and left chiral charges along the magnetic field; the resulting dipolar axial charge density in turn induces the oppositely directed vector charge currents leading to an electric quadrupole moment of the quark-gluon plasma. Boosted by the strong collective flow, the electric quadrupole moment translates into the charge dependence of the elliptic flow coefficients, so that $v_2(\pi^+) < v_2(\pi^-)$ (at positive net charge). Using the latest quantitative simulations of the produced magnetic field and solving the CMW equation, we make further quantitative estimates of the produced $v_2$ splitting and its centrality dependence. We compare the results with the available experimental data.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.

Effects of Thoratec pulsatile ventricular assist device (PVAD) timing on the abdominal aortic wave intensity pattern

Arterial waves are seen as possible independent mediators of cardiovascular risks, and the wave intensity analysis (WIA) has therefore been proposed as a method for patient selection for ventricular assist device (VAD) implantation. Interpreting measured wave intensity (WI) is challenging and complexity is increased by the implantation of a VAD. The waves generated by the VAD interact with the waves generated by the native heart, and this interaction varies with changing VAD settings. Eight sheep were implanted with a pulsatile VAD (PVAD) through ventriculo-aortic cannulation. The start of PVAD ejection was synchronized to the native R-wave and delayed between 0 % - 90 % of the cardiac cycle in 10 % steps or phase shifts (PS). Pressure and velocity signals were registered, using a combined Doppler and pressure wire positioned in the abdominal aorta, and used to calculate the WI. Depending on the PS, different wave interference phenomena occurred. Maximum unloading of the left ventricle (LV) coincided with constructive interference and maximum blood flow pulsatility, and maximum loading of the LV coincided with destructive interference and minimum blood flow pulsatility. We believe, that non-invasive WIA could potentially be used clinically to assess the mechanical load of the LV, and to monitor the peripheral hemodynamics such as blood flow pulsatility and risk of intestinal bleeding.

Bed forms generated in the laboratory under an oscillatory flow : analytical and experimental study /

"June 1969."