Nonlinear stability of multilayer quasi-geostrophic flow

New nonlinear stability theorems are derived for disturbances to steady basic flows in the context of the multilayer quasi-geostrophic equations. These theorems are analogues of Arnol’d's second stability theorem, the latter applying to the two-dimensional Euler equations. Explicit upper bounds are obtained on both the disturbance energy and disturbance potential enstrophy in terms of the initial disturbance fields. An important feature of the present analysis is that the disturbances are allowed to have non-zero circulation. While Arnol’d's stability method relies on the energy–Casimir invariant being sign-definite, the new criteria can be applied to cases where it is sign-indefinite because of the disturbance circulations. A version of Andrews’ theorem is established for this problem, and uniform potential vorticity flow is shown to be nonlinearly stable. The special case of two-layer flow is treated in detail, with particular attention paid to the Phillips model of baroclinic instability. It is found that the short-wave portion of the marginal stability curve found in linear theory is precisely captured by the new nonlinear stability criteria.

The stability of a two-dimensional vorticity filament under uniform strain

The quantitative effects of uniform strain and background rotation on the stability of a strip of constant vorticity (a simple shear layer) are examined. The thickness of the strip decreases in time under the strain, so it is necessary to formulate the linear stability analysis for a time-dependent basic flow. The results show that even a strain rate γ (scaled with the vorticity of the strip) as small as 0.25 suppresses the conventional Rayleigh shear instability mechanism, in the sense that the r.m.s. wave steepness cannot amplify by more than a certain factor, and must eventually decay. For γ < 0.25 the amplification factor increases as γ decreases; however, it is only 3 when γ e 0.065. Numerical simulations confirm the predictions of linear theory at small steepness and predict a threshold value necessary for the formation of coherent vortices. The results help to explain the impression from numerous simulations of two-dimensional turbulence reported in the literature that filaments of vorticity infrequently roll up into vortices. The stabilization effect may be expected to extend to two- and three-dimensional quasi-geostrophic flows.

Pulsating equilibria: stability through migration

This paper is to present a model of spatial equilibrium using a nonlinear generalization of Markov-chain type model, and to show the dynamic stability of a unique equilibrium. Even at an equilibrium, people continue to migrate among regions as well as among agent-types, and yet their overall distribution remain unchanged. The model is also adapted to suggest a theory of traffic distribution in a city.

On the role of the ocean in projected atmospheric stability changes in the Atlantic polar low region

The occurrence of destructive mesoscale ‘polar low’ cyclones in the subpolar North Atlantic is projected to decline under anthropogenic change, due to an increase in atmospheric static stability. This letter reports on the role of changes in ocean circulation in shaping the atmospheric stability. In particular, the Atlantic Meridional Overturning Circulation (AMOC) is projected to weaken in response to anthropogenic forcing, leading to a local minimum in warming in this region. The reduced warming is restricted to the lower troposphere, hence contributing to the increase in static stability. Linear correlation analysis of the CMIP3 climate model ensemble suggests that around half of the model uncertainty in the projected stability response arises from the varied response of the AMOC between models.

The importance of loop length on the stability of i-motif structures

Using UV and srCD spectroscopy it is found that loop length within the i-motif structure is important for both thermal and pH stability, but in contrast to previous statements, it is the shorter loops that exhibit the highest stability.

Structural stability of the synthetic thermoelectric ternary and nickel-substituted tetrahedrite phases

The purity and structural stability of the high thermoelectric performance Cu12Sb4S13 and Cu10.4Ni1.6Sb4S13 tetrahedrite phases, synthesized by solid–liquid–vapor reaction and Spark Plasma Sintering, were studied at high temperature by Rietveld refinement using high resolution X-ray powder diffraction data, DSC/TG measurements and high resolution transmission electron microscopy. In a complementary study, the crystal structure of Cu10.5Ni1.5Sb4S13 as a function of temperature was investigated by powder neutron diffraction. The temperature dependence of the structural stability of ternary Cu12Sb4S13 is markedly different to that of the nickel-substituted phases, providing clear evidence for the significant and beneficial role of nickel substitution on both sample purity and stability of the tetrahedrite phase. Moreover, kinetic effects on the phase stability/decomposition have been identified and discussed in order to determine the maximum operating temperature for thermoelectric applications. The thermoelectric properties of these compounds have been determined for high density samples (>98%) prepared by Spark Plasma Sintering and therefore can be used as reference values for tetrahedrite samples. The maximum ZT of 0.8 was found for Cu10.4Ni1.6Sb4S13 at 700 K.

Predicting mesh density for adaptive modelling of the global atmosphere

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.

The triggering of orographic rainbands by small-scale topography

The triggering of convective orographic rainbands by small-scale topographic features is investigated through observations of a banded precipitation event over the Oregon Coastal Range and simulations using a cloud-resolving numerical model. A quasi-idealized simulation of the observed event reproduces the bands in the radar observations, indicating the model’s ability to capture the physics of the band-formation process. Additional idealized simulations reinforce that the bands are triggered by lee waves past small-scale topographic obstacles just upstream of the nominal leading edge of the orographic cloud. Whether a topographic obstacle in this region is able to trigger a strong rainband depends on the phase of its lee wave at cloud entry. Convective growth only occurs downstream of obstacles that give rise to lee-wave-induced displacements that create positive vertical velocity anomalies w_c and nearly zero buoyancy anomalies b_c as air parcels undergo saturation. This relationship is quantified through a simple analytic condition involving w_c, b_c, and the static stability N_m^2 of the cloud mass. Once convection is triggered, horizontal buoyancy gradients in the cross-flow direction generate circulations that align the bands parallel to the flow direction.

Converting the nanodomains of a diblock-copolymer thin film from spheres to cylinders with an external electric field

We investigate the ability of an applied electric field to convert the morphology of a diblock-copolymer thin film from a monolayer of spherical domains embedded in the matrix to cylindrical domains that penetrate through the matrix. As expected, the applied field increases the relative stability of cylindrical domains, while simultaneously reducing the energy barrier that impedes the transition to cylinders. The effectiveness of the field is enhanced by a large dielectric contrast between the two block-copolymer components, particularly when the low-dielectric contrast component forms the matrix. Furthermore, the energy barrier is minimized by selecting sphere-forming diblock copolymers that are as compositionally symmetric as possible. Our calculations, which are the most quantitatively reliable to date, are performed using a numerically precise spectral algorithm based on self-consistent-field theory supplemented with an exact treatment for linear dielectric materials.

Consistent approximate models of the global atmosphere: shallow, deep, hydrostatic, quasi-hydrostatic and non-hydrostatic

We study global atmosphere models that are at least as accurate as the hydrostatic primitive equations (HPEs), reviewing known results and reporting some new ones. The HPEs make spherical geopotential and shallow atmosphere approximations in addition to the hydrostatic approximation. As is well known, a consistent application of the shallow atmosphere approximation requires omission of those Coriolis terms that vary as the cosine of latitude and of certain other terms in the components of the momentum equation. An approximate model is here regarded as consistent if it formally preserves conservation principles for axial angular momentum, energy and potential vorticity, and (following R. Müller) if its momentum component equations have Lagrange's form. Within these criteria, four consistent approximate global models, including the HPEs themselves, are identified in a height-coordinate framework. The four models, each of which includes the spherical geopotential approximation, correspond to whether the shallow atmosphere and hydrostatic (or quasi-hydrostatic) approximations are individually made or not made. Restrictions on representing the spatial variation of apparent gravity occur. Solution methods and the situation in a pressure-coordinate framework are discussed. © Crown copyright 2005.

A comparison of two methods for developing the linearization of a shallow-water model

We develop the linearization of a semi-implicit semi-Lagrangian model of the one-dimensional shallow-water equations using two different methods. The usual tangent linear model, formed by linearizing the discrete nonlinear model, is compared with a model formed by first linearizing the continuous nonlinear equations and then discretizing. Both models are shown to perform equally well for finite perturbations. However, the asymptotic behaviour of the two models differs as the perturbation size is reduced. This leads to difficulties in showing that the models are correctly coded using the standard tests. To overcome this difficulty we propose a new method for testing linear models, which we demonstrate both theoretically and numerically. © Crown copyright, 2003. Royal Meteorological Society

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.

Ionic calcium and pH as predictors of stability of milk to UHT processing

This paper describes the use of pH and calcium ion electrodes for investigating factors affecting the heat stability of UHT milk with added calcium chloride. Calcium chloride was added to raw milk to manipulate ionic calcium and pH to within the range that may be typically encountered in raw milk of different compositions and microbial quality. Addition of only 5 mM calcium chloride was sufficient to induce considerable changes in pH, ionic calcium and ethanol stability and alter its stability to UHT treatment. There was a strong relationship between pH decrease and increase in ionic calcium when pH was reduced, whether by addition of calcium chloride or by acidification. Calcium chloride addition was found to increase sediment formation in UHT treated milk. However, sediment could be reduced by addition of stabilizers. Those most effective were ones which decreased ionic calcium and increased pH, such as trisodium citrate and disodium hydrogen phosphate. Sediment formation following UHT treatment was only slight for milk samples whose ethanol stability was greater than 80%.

Resolving the energy paradox of chaperone/usher-mediated fibre assembly

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone-subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone-subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone-fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone-subunit preassembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation.

Effects of copper on the antioxidant activity of olive polyphenols in bulk oil and oil-in-water emulsions

The antioxidant activity and interactions with copper of four olive oil phenolic compounds, namely oleuropein, hydroxytyrosol, 3,4- dihydroxyphenylethanol- elenolic acid ( 1), and 3,4- dihydroxyphenyl-ethanolelenolic acid dialdehyde ( 2), in olive oil and oil- in- water emulsions stored at 60 degrees C were studied. All four phenolic compounds significantly extended the induction time of lipid oxidation in olive oil with the order of activity being hydroxytyrosol > compound 1 > compound 2 > oleuropein > alpha- tocopherol; but in the presence of Cu( II), the stability of oil samples containing phenolic compounds decreased by at least 90%, and the antioxidant activity of hydroxytyrosol and compounds 1 and 2 became similar. In oil- in- water emulsions prepared from olive oil stripped of tocopherols, hydroxytyrosol enhanced the prooxidant effect of copper at pH 5.5 but not at pH 7.4. The stability of samples containing copper at pH 5.5 was not significantly different if oleuropein was present from that of the control. Oleuropein at pH 7.4, and compounds 1 and 2 at both pH values tested, reduced the prooxidant effect of copper. The lower stability and the higher reducing capacity of all compounds at pH 7.4 could not explain the higher stability of emulsions containing phenolic compounds at this pH value. However, mixtures containing hydroxytyrosol or oleuropein with copper showed higher 1,1-diphenyl- 2- picrylhydrazyl radical scavenging activity at pH 7.4 than at pH 5.5. Moreover, the compound 2- copper complex showed higher radical scavenging activity then the uncomplexed compound at pH 5.5. It can be concluded that the formation of a copper complex with radical scavenging activity is a key step in the antioxidant action of the olive oil phenolic compounds in an emulsion containing copper ions.