Weak constraints in four-dimensional variational data assimilation

The formulation of four-dimensional variational data assimilation allows the incorporation of constraints into the cost function which need only be weakly satisfied. In this paper we investigate the value of imposing conservation properties as weak constraints. Using the example of the two-body problem of celestial mechanics we compare weak constraints based on conservation laws with a constraint on the background state.We show how the imposition of conservation-based weak constraints changes the nature of the gradient equation. Assimilation experiments demonstrate how this can add extra information to the assimilation process, even when the underlying numerical model is conserving.

How to save a bad element with weak boundary conditions

The P-1-P-1 finite element pair is known to allow the existence of spurious pressure (surface elevation) modes for the shallow water equations and to be unstable for mixed formulations. We show that this behavior is strongly influenced by the strong or the weak enforcement of the impermeability boundary conditions. A numerical analysis of the Stommel model is performed for both P-1-P-1 and P-1(NC)-P-1 mixed formulations. Steady and transient test cases are considered. We observe that the P-1-P-1 element exhibits stable discrete solutions with weak boundary conditions or with fully unstructured meshes. (c) 2005 Elsevier Ltd. All rights reserved.

A weak formulation of Roe's approximate Riemann Solver applied to 'Barotropic' flows

A weak formulation of Roe's approximate Riemann solver is applied to the equations of ‘barotropic’ flow, including the shallow water equations, and it is shown that this leads to an approximate Riemann solver recently presented for such flows.

Equivalence of weak formulations of the steady water waves equations

We prove the equivalence of three weak formulations of the steady water waves equations, namely: the velocity formulation, the stream function formulation and the Dubreil-Jacotin formulation, under weak Hölder regularity assumptions on their solutions.

Cloud-resolving model simulations with one- and two-way couplings via the weak temperature gradient approximation

A cloud-resolving model is modified to implement the weak temperature gradient approximation in order to simulate the interactions between tropical convection and the large-scale tropical circulation. The instantaneous domain-mean potential temperature is relaxed toward a reference profile obtained from a radiative–convective equilibrium simulation of the cloud-resolving model. For homogeneous surface conditions, the model state at equilibrium is a large-scale circulation with its descending branch in the simulated column. This is similar to the equilibrium state found in some other studies, but not all. For this model, the development of such a circulation is insensitive to the relaxation profile and the initial conditions. Two columns of the cloud-resolving model are fully coupled by relaxing the instantaneous domain-mean potential temperature in both columns toward each other. This configuration is energetically closed in contrast to the reference-column configuration. No mean large-scale circulation develops over homogeneous surface conditions, regardless of the relative area of the two columns. The sensitivity to nonuniform surface conditions is similar to that obtained in the reference-column configuration if the two simulated columns have very different areas, but it is markedly weaker for columns of comparable area. The weaker sensitivity can be understood as being a consequence of a formulation for which the energy budget is closed. The reference-column configuration has been used to study the convection in a local region under the influence of a large-scale circulation. The extension to a two-column configuration is proposed as a methodology for studying the influence on local convection of changes in remote convection.

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.

The effect of copper(II), iron(II) sulphate and vitamin C combinations on the weak antimicrobial activity of (+)-catechin against Staphylococcus aureus and other microbes

Few attempts have been made to improve the activity of plant compounds with low antimicrobial efficacy. (+)-Catechin, a weak antimicrobial tea flavanol, was combined with putative adjuncts and tested against different species of bacteria. Copper(II) sulphate enhanced (+)-catechin activity against Pseudomonas aeruginosa but not Staphylococcus aureus, Proteus mirabilis or Escherichia coli. Attempts to raise the activity of (+)-catechin against two unresponsive species, S. aureus and E. coli, with iron(II) sulphate, iron(III) chloride, and vitamin C, showed that iron(II) enhanced (+)-catechin against S. aureus, but not E. coli; neither iron(III) nor combined iron(II) and copper(II), enhanced (+)-catechin activity against either species. Vitamin C enhanced copper(II) containing combinations against both species in the absence of iron(II). Catalase or EDTA added to active samples removed viability effects suggesting that active mixtures had produced H2O2via the action of added metal(II) ions. H2O2 generation by (+)-catechin plus copper(II) mixtures and copper(II) alone could account for the principal effect of bacterial growth inhibition following 30 minute exposures as well as the antimicrobial effect of (+)-catechin–iron(II) against S. aureus. These novel findings about a weak antimicrobial flavanol contrast with previous knowledge of more active flavanols with transition metal combinations. Weak antimicrobial compounds like (+)-catechin within enhancement mixtures may therefore be used as efficacious agents. (+)-Catechin may provide a means of lowering copper(II) or iron(II) contents in certain crop protection and other products.