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 unified theory of balance in the extratropics

Many physical systems exhibit dynamics with vastly different time scales. Often the different motions interact only weakly and the slow dynamics is naturally constrained to a subspace of phase space, in the vicinity of a slow manifold. In geophysical fluid dynamics this reduction in phase space is called balance. Classically, balance is understood by way of the Rossby number R or the Froude number F; either R ≪ 1 or F ≪ 1. We examined the shallow-water equations and Boussinesq equations on an f -plane and determined a dimensionless parameter _, small values of which imply a time-scale separation. In terms of R and F, ∈= RF/√(R^2+R^2 ) We then developed a unified theory of (extratropical) balance based on _ that includes all cases of small R and/or small F. The leading-order systems are ensured to be Hamiltonian and turn out to be governed by the quasi-geostrophic potential-vorticity equation. However, the height field is not necessarily in geostrophic balance, so the leading-order dynamics are more general than in quasi-geostrophy. Thus the quasi-geostrophic potential-vorticity equation (as distinct from the quasi-geostrophic dynamics) is valid more generally than its traditional derivation would suggest. In the case of the Boussinesq equations, we have found that balanced dynamics generally implies hydrostatic balance without any assumption on the aspect ratio; only when the Froude number is not small and it is the Rossby number that guarantees a timescale separation must we impose the requirement of a small aspect ratio to ensure hydrostatic balance.

On the interpretation of Andrews’ theorem

Andrews (1984) has shown that any flow satisfying Arnol'd's (1965, 1966) sufficient conditions for stability must be zonally-symmetric if the boundary conditions on the flow are zonally-symmetric. This result appears to place very strong restrictions on the kinds of flows that can be proved to be stable by Arnol'd's theorems. In this paper, Andrews’ theorem is re-examined, paying special attention to the case of an unbounded domain. It is shown that, in that case, Andrews’ theorem generally fails to apply, and Arnol'd-stable flows do exist that are not zonally-symmetric. The example of a circular vortex with a monotonic vorticity profile is a case in point. A proof of the finite-amplitude version of the Rayleigh stability theorem for circular vortices is also established; despite its similarity to the Arnol'd theorems it seems not to have been put on record before.

On Rossby waves modified by weak sinusoidal shear

The spatial structure of beta-plane Rossby waves in a sinusoidal basic zonal flow U 0cos(γ,y) is determined analytically in the (stable) asymptotic limit of weak shear, U 0γ2 0/β≈1. The propagating neutral normal modes are found to take their greatest amplitude in the region of maximum westerly flow, while their most rapid phase variation is achieved in the region of maximum easterly flow. These results are shown to be consistent with what is obtained by ray-tracing methods in the limit of small meridional disturbance wavelength.

Mean motions induced by baroclinic instability in a jet

A study is made of the zonal-mean motions induced by a growing baroclinic wave in several contexts, under the framework of three different analysis schemes: the conventional Eulerian mean (EM), the transformed Eulerian mean (TEM), and the generalized Lagrangian mean (GLM). The effect of meridional shear in the initial jet on these induced mean motions is considered by treating the instability problem in the context of the two-layer model. The conceptual simplicity of the TEM formulation is shown to be useful in diagnosing the dynamics of instability, much as it has been found helpful in many problems of wave, mean-flow interaction. In addition, it is found that the TEM vertical velocity is a very good indicator of the GLM vertical velocity. However, the GLM meridional velocity is always convergent towards the centre of instability activity, and is not at all well represented by the nondivergent TEM meridional velocity. In comparing the results with Uryu's (1979) calculation of the GLM circulation induced by a growing Eady wave, it is found that the inclusion of meridional jet shear in the present work leads to some strikingly different effects in the GLM zonal wind acceleration. In the case of pure baroclinic instability treated by Uryu, the Eulerian and Stokes accelerations nearly cancel each other in the centre of the channel, leaving a weak Lagrangian acceleration opposed to the Eulerian one. In the more general case of mixed baroclinic-barotropic instability, however, the Eulerian and Stokes accelerations can reinforce one another, leading to a very strong Lagrangian zonal wind

Comments on 'Balance and the slow quasimanifold: some explicit results'

The concept of slow vortical dynamics and its role in theoretical understanding is central to geophysical fluid dynamics. It leads, for example, to “potential vorticity thinking” (Hoskins et al. 1985). Mathematically, one imagines an invariant manifold within the phase space of solutions, called the slow manifold (Leith 1980; Lorenz 1980), to which the dynamics are constrained. Whether this slow manifold truly exists has been a major subject of inquiry over the past 20 years. It has become clear that an exact slow manifold is an exceptional case, restricted to steady or perhaps temporally periodic flows (Warn 1997). Thus the concept of a “fuzzy slow manifold” (Warn and Ménard 1986) has been suggested. The idea is that nearly slow dynamics will occur in a stochastic layer about the putative slow manifold. The natural question then is, how thick is this layer? In a recent paper, Ford et al. (2000) argue that Lighthill emission—the spontaneous emission of freely propagating acoustic waves by unsteady vortical flows—is applicable to the problem of balance, with the Mach number Ma replaced by the Froude number F, and that it is a fundamental mechanism for this fuzziness. They consider the rotating shallow-water equations and find emission of inertia–gravity waves at O(F2). This is rather surprising at first sight, because several studies of balanced dynamics with the rotating shallow-water equations have gone beyond second order in F, and found only an exponentially small unbalanced component (Warn and Ménard 1986; Lorenz and Krishnamurthy 1987; Bokhove and Shepherd 1996; Wirosoetisno and Shepherd 2000). We have no technical objection to the analysis of Ford et al. (2000), but wish to point out that it depends crucially on R 1, where R is the Rossby number. This condition requires the ratio of the characteristic length scale of the flow L to the Rossby deformation radius LR to go to zero in the limit F → 0. This is the low Froude number scaling of Charney (1963), which, while originally designed for the Tropics, has been argued to be also relevant to mesoscale dynamics (Riley et al. 1981). If L/LR is fixed, however, then F → 0 implies R → 0, which is the standard quasigeostrophic scaling of Charney (1948; see, e.g., Pedlosky 1987). In this limit there is reason to expect the fuzziness of the slow manifold to be “exponentially thin,” and balance to be much more accurate than is consistent with (algebraic) Lighthill emission.

Evidence of component merging equatorward of the cusp

The Polar spacecraft passed through a region near the dayside magnetopause on May 29, 1996, at a geocentric distance of similar to 8 R-E and high, northern magnetic latitudes. The interplanetary magnetic field (IMF) was northward during the pass. Data from the Thermal Ion Dynamics Experiment revealed the existence of low-speed (similar to 50 km s(-1)) ion D-shaped distributions mixed with cold ions (similar to 2 eV) over a period of 2.5 hours. These ions were traveling parallel to the magnetic field toward the Northern Hemisphere ionosphere and were convecting primarily eastward. The D-shaped distributions are distinct from a convecting Maxwellian and, along with the magnetic field direction, are taken as evidence that the spacecraft was inside the magnetosphere and not in the magnetosheath. Furthermore, the absence of ions in the antiparallel direction is taken as evidence that low-shear merging was occurring at a location southward of the spacecraft and equatorward of the Southern Hemisphere cusp. The cold ions were of ionospheric origin, with initially slow field-aligned speeds, which were accelerated upon reflection from the magnetopause. These observations provide significant new evidence consistent with component magnetic merging sites equatorward of the cusp for northward IMF.

The non-conservation of potential vorticity by a dynamical core compared with the effects of parametrized physical processes

Numerical models of the atmosphere combine a dynamical core, which approximates solutions to the adiabatic, frictionless governing equations for fluid dynamics, with tendencies arising from the parametrization of other physical processes. Since potential vorticity (PV) is conserved following fluid flow in adiabatic, frictionless circumstances, it is possible to isolate the effects of non-conservative processes by accumulating PV changes in an air-mass relative framework. This “PV tracer technique” is used to accumulate separately the effects on PV of each of the different non-conservative processes represented in a numerical model of the atmosphere. Dynamical cores are not exactly conservative because they introduce, explicitly or implicitly, some level of dissipation and adjustment of prognostic model variables which acts to modify PV. Here, the PV tracers technique is extended to diagnose the cumulative effect of the non-conservation of PV by a dynamical core and its characteristics relative to the PV modification by parametrized physical processes. Quantification using the Met Office Unified Model reveals that the magnitude of the non-conservation of PV by the dynamical core is comparable to those from physical processes. Moreover, the residual of the PV budget, when tracing the effects of the dynamical core and physical processes, is at least an order of magnitude smaller than the PV tracers associated with the most active physical processes. The implication of this work is that the non-conservation of PV by a dynamical core can be assessed in case studies with a full suite of physics parametrizations and directly compared with the PV modification by parametrized physical processes. The nonconservation of PV by the dynamical core is shown to move the position of the extratropical tropopause while the parametrized physical processes have a lesser effect at the tropopause level.

Molecular detection of enterotoxins in environmental strains of Aeromonas hydrophila and Aeromonas jandaei

Aeromonas species are widely distributed in aquatic environments and recent studies include the genus in the emergent pathogens group because of its frequent association with local and systemic infections in immunocompetent humans. Aiming to search for virulence genes in environmental strains of Aeromonas hydrophila and Aeromonas jandaei, we designed specific primers to detect act/hly A/aer complex and alt genes. Primers described elsewhere were used to detect ast. Eighty-seven strains previously identified using phenotypic and genotypic tests as A. hydrophila (41) and A. jandaei (46) were analysed for the presence of the virulence genes using PCR. DNA fragments of expected size were purified and directly sequenced. Among the 41 strains of A. hydrophila 70.7% (29), 97.6% (40) and 26.8% (11) possessed act/hly A/aer complex, ast and alt genes, respectively. Among the 46 strains of A. jandaei, 4.4% (2), 0% (0) and 32.6% (15) were positive for act/hly A/aer complex, ast and alt genes, respectively. Sequencing allowed for the confirmation of amplified products using BLAST. The present work proposes a specific and rapid diagnostic method to detect the main virulence determinants of Aeromonas, a genus potentially pathogenic to humans.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

Front tracking with moving-least-squares surfaces

The representation of interfaces by means of the algebraic moving-least-squares (AMLS) technique is addressed. This technique, in which the interface is represented by an unconnected set of points, is interesting for evolving fluid interfaces since there is]to surface connectivity. The position of the surface points can thus be updated without concerns about the quality of any surface triangulation. We introduce a novel AMLS technique especially designed for evolving-interfaces applications that we denote RAMLS (for Robust AMLS). The main advantages with respect to previous AMLS techniques are: increased robustness, computational efficiency, and being free of user-tuned parameters. Further, we propose a new front-tracking method based on the Lagrangian advection of the unconnected point set that defines the RAMLS surface. We assume that a background Eulerian grid is defined with some grid spacing h. The advection of the point set makes the surface evolve in time. The point cloud can be regenerated at any time (in particular, we regenerate it each time step) by intersecting the gridlines with the evolved surface, which guarantees that the density of points on the surface is always well balanced. The intersection algorithm is essentially a ray-tracing algorithm, well-studied in computer graphics, in which a line (ray) is traced so as to detect all intersections with a surface. Also, the tracing of each gridline is independent and can thus be performed in parallel. Several tests are reported assessing first the accuracy of the proposed RAMLS technique, and then of the front-tracking method based on it. Comparison with previous Eulerian, Lagrangian and hybrid techniques encourage further development of the proposed method for fluid mechanics applications. (C) 2008 Elsevier Inc. All rights reserved.

PLANT LEAF IDENTIFICATION BASED ON VOLUMETRIC FRACTAL DIMENSION

Texture is an important visual attribute used to describe the pixel organization in an image. As well as it being easily identified by humans, its analysis process demands a high level of sophistication and computer complexity. This paper presents a novel approach for texture analysis, based on analyzing the complexity of the surface generated from a texture, in order to describe and characterize it. The proposed method produces a texture signature which is able to efficiently characterize different texture classes. The paper also illustrates a novel method performance on an experiment using texture images of leaves. Leaf identification is a difficult and complex task due to the nature of plants, which presents a huge pattern variation. The high classification rate yielded shows the potential of the method, improving on traditional texture techniques, such as Gabor filters and Fourier analysis.

Performance evaluation of ventilation radiators

A ventilation radiator is a combined ventilation and heat emission unit currently of interest due to its potential for increasing energy efficiency in exhaust ventilated buildings with warm water heating. This paper presents results of performance tests of several ventilation radiator models conducted under controlled laboratory conditions.   The purpose of the study was to validate results achieved by Computational Fluid Dynamics (CFD) in an earlier study and indentify possible improvements in the performance of such systems. The main focus was on heat transfer from internal convection fins, but comfort and health aspects related to ventilation rates and air temperatures were also considered.   The general results from the CFD simulations were confirmed; the heat output of ventilation radiators may be improved by at least 20 % without sacrificing ventilation efficiency or thermal comfort.   Improved thermal efficiency of ventilation radiators allows a lower supply water temperature and energy savings both for heating up and distribution of warm water in heat pumps or district heating systems. A secondary benefit is that a high ventilation rate can be maintained all year around without risk for cold draught.

Radiation properties of coil-coated steel in building envelope surfaces and the influence on building thermal performance

Recent studies have shown that the optical properties of building exterior surfaces are important in terms of energy use and thermal comfort. While the majority of the studies are related to exterior surfaces, the radiation properties of interior surfaces are less thoroughly investigated. Development in the coil-coating industries has now made it possible to allocate different optical properties for both exterior and interior surfaces of steel-clad buildings. The aim of this thesis is to investigate the influence of surface radiation properties with the focus on the thermal emittance of the interior surfaces, the modeling approaches and their consequences in the context of the building energy performance and indoor thermal environment. The study consists of both numerical and experimental investigations. The experimental investigations include parallel field measurements on three similar test cabins with different interior and exterior surface radiation properties in Borlänge, Sweden, and two ice rink arenas with normal and low emissive ceiling in Luleå, Sweden. The numerical methods include comparative simulations by the use of dynamic heat flux models, Building Energy Simulation (BES), Computational Fluid Dynamics (CFD) and a coupled model for BES and CFD. Several parametric studies and thermal performance analyses were carried out in combination with the different numerical methods. The parallel field measurements on the test cabins include the air, surface and radiation temperatures and energy use during passive and active (heating and cooling) measurements. Both measurement and comparative simulation results indicate an improvement in the indoor thermal environment when the interior surfaces have low emittance. In the ice rink arenas, surface and radiation temperature measurements indicate a considerable reduction in the ceiling-to-ice radiation by the use of low emittance surfaces, in agreement with a ceiling-toice radiation model using schematic dynamic heat flux calculations. The measurements in the test cabins indicate that the use of low emittance surfaces can increase the vertical indoor air temperature gradients depending on the time of day and outdoor conditions. This is in agreement with the transient CFD simulations having the boundary condition assigned on the exterior surfaces. The sensitivity analyses have been performed under different outdoor conditions and surface thermal radiation properties. The spatially resolved simulations indicate an increase in the air and surface temperature gradients by the use of low emittance coatings. This can allow for lower air temperature at the occupied zone during the summer. The combined effect of interior and exterior reflective coatings in terms of energy use has been investigated by the use of building energy simulation for different climates and internal heat loads. The results indicate possible energy savings by the smart choice of optical properties on interior and exterior surfaces of the building. Overall, it is concluded that the interior reflective coatings can contribute to building energy savings and improvement of the indoor thermal environment. This can be numerically investigated by the choice of appropriate models with respect to the level of detail and computational load. This thesis includes comparative simulations at different levels of detail.