On wave action and phase in the non-canonical Hamiltonian formulation

The long time–evolution of disturbances to slowly–varying solutions of partial differential equations is subject to the adiabatic invariance of the wave action. Generally, this approximate conservation law is obtained under the assumption that the partial differential equations are derived from a variational principle or have a canonical Hamiltonian structure. Here, the wave action conservation is examined for equations that possess a non–canonical (Poisson) Hamiltonian structure. The linear evolution of disturbances in the form of slowly varying wavetrains is studied using a WKB expansion. The properties of the original Hamiltonian system strongly constrain the linear equations that are derived, and this is shown to lead to the adiabatic invariance of a wave action. The connection between this (approximate) invariance and the (exact) conservation laws of pseudo–energy and pseudomomentum that exist when the basic solution is exactly time and space independent is discussed. An evolution equation for the slowly varying phase of the wavetrain is also derived and related to Berry's phase.

A search for invariant relative satellite motion

This report presents the canonical Hamiltonian formulation of relative satellite motion. The unperturbed Hamiltonian model is shown to be equivalent to the well known Hill-Clohessy-Wilshire (HCW) linear formulation. The in°uence of perturbations of the nonlinear Gravitational potential and the oblateness of the Earth; J2 perturbations are also modelled within the Hamiltonian formulation. The modelling incorporates eccentricity of the reference orbit. The corresponding Hamiltonian vector ¯elds are computed and implemented in Simulink. A numerical method is presented aimed at locating periodic or quasi-periodic relative satellite motion. The numerical method outlined in this paper is applied to the Hamiltonian system. Although the orbits considered here are weakly unstable at best, in the case of eccentricity only, the method ¯nds exact periodic orbits. When other perturbations such as nonlinear gravitational terms are added, drift is signicantly reduced and in the case of the J2 perturbation with and without the nonlinear gravitational potential term, bounded quasi-periodic solutions are found. Advantages of using Newton's method to search for periodic or quasi-periodic relative satellite motion include simplicity of implementation, repeatability of solutions due to its non-random nature, and fast convergence. Given that the use of bounded or drifting trajectories as control references carries practical di±culties over long-term missions, Principal Component Analysis (PCA) is applied to the quasi-periodic or slowly drifting trajectories to help provide a closed reference trajectory for the implementation of closed loop control. In order to evaluate the e®ect of the quality of the model used to generate the periodic reference trajectory, a study involving closed loop control of a simulated master/follower formation was performed. 2 The results of the closed loop control study indicate that the quality of the model employed for generating the reference trajectory used for control purposes has an important in°uence on the resulting amount of fuel required to track the reference trajectory. The model used to generate LQR controller gains also has an e®ect on the e±ciency of the controller.

The physical pendulum in an advanced undergraduate course in mechanics

The physical pendulum treated with a Hamiltonian formulation is a natural topic for study in a course in advanced classical mechanics. For the past three years, we have been offering a series of problem sets studying this system numerically in our third-year undergraduate courses in mechanics. The problem sets investigate the physics of the pendulum in ways not easily accessible without computer technology and explore various algorithms for solving mechanics problems. Our computational physics is based on Mathematica with some C communicating with Mathematica, although nothing in this paper is dependent on that choice. We have nonetheless found this system, and particularly its graphics, to be a good one for use with undergraduates.

A well-posed integral equation formulation for three-dimensional rough surface scattering

We consider the problem of scattering of time-harmonic acoustic waves by an unbounded sound-soft rough surface. Recently, a Brakhage Werner type integral equation formulation of this problem has been proposed, based on an ansatz as a combined single- and double-layer potential, but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half-space Green's function. Moreover, it has been shown in the three-dimensional case that this integral equation is uniquely solvable in the space L-2 (Gamma) when the scattering surface G does not differ too much from a plane. In this paper, we show that this integral equation is uniquely solvable with no restriction on the surface elevation or slope. Moreover, we construct explicit bounds on the inverse of the associated boundary integral operator, as a function of the wave number, the parameter coupling the single- and double-layer potentials, and the maximum surface slope. These bounds show that the norm of the inverse operator is bounded uniformly in the wave number, kappa, for kappa > 0, if the coupling parameter h is chosen proportional to the wave number. In the case when G is a plane, we show that the choice eta = kappa/2 is nearly optimal in terms of minimizing the condition number.

Evaluating a new LISFLOOD-FP formulation with data from the summer 2007 floods in Tewkesbury, UK

Two-dimensional flood inundation modelling is a widely used tool to aid flood risk management. In urban areas, where asset value and population density are greatest, the model spatial resolution required to represent flows through a typical street network (i.e. < 10m) often results in impractical computational cost at the whole city scale. Explicit diffusive storage cell models become very inefficient at such high resolutions, relative to shallow water models, because the stable time step in such schemes scales as a quadratic of resolution. This paper presents the calibration and evaluation of a recently developed new formulation of the LISFLOOD-FP model, where stability is controlled by the Courant–Freidrichs–Levy condition for the shallow water equations, such that, the stable time step instead scales linearly with resolution. The case study used is based on observations during the summer 2007 floods in Tewkesbury, UK. Aerial photography is available for model evaluation on three separate days from the 24th to the 31st of July. The model covered a 3.6 km by 2 km domain and was calibrated using gauge data from high flows during the previous month. The new formulation was benchmarked against the original version of the model at 20 m and 40 m resolutions, demonstrating equally accurate performance given the available validation data but at 67x faster computation time. The July event was then simulated at the 2 m resolution of the available airborne LiDAR DEM. This resulted in a significantly more accurate simulation of the drying dynamics compared to that simulated by the coarse resolution models, although estimates of peak inundation depth were similar.

Microwave spectra and centrifugal distortion constants of oxetane

The microwave spectra of oxetane (trimethylene oxide) and its three symmetrically deuterated isotopic species have been observed on a Hewlett-Packard microwave spectrometer from 26.5 to 40 GHz. For the parent species, the β-d2 and the αα′-d4 species, about 300 lines have been assigned for each molecule, and for the d6 species more than 600 lines have been assigned. The assignments range from v = 0 to v = 5 in the puckering vibration; although they are mostly Q transitions, either 3 or 4 R transitions have been observed for each vibrational state. The spectra have been interpreted using an effective rotational hamiltonian for each vibrational state, including five quartic distortion constants according to Watson's formulation, and a variable number of sextic distortion constants; in general, the lines are fitted to about ± 10 kHz. The distortion constants show an anomalous zig-zag dependence on the puckering vibrational quantum number, similar to that first observed for the rotational constants by Gwinn and coworkers. This is interpreted according to a simple modification of the standard theory of centrifugal distortion, involving the double minimum potential function in the puckering coordinate.

Formulation of a live bacterial vaccine for stable room temperature storage results in loss of acid, bile and bile salt resistance

Live bacterial vaccines have great promise both as vaccines against enteric pathogens and as heterologous antigen vectors against diverse diseases. Ideally, room temperature stable dry formulations of live bacterial vaccines will allow oral vaccination without cold-chain storage or injections. Attenuated Salmonella can cross the intestinal wall and deliver replicating antigen plus innate immune activation signals directly to the intestinal immune tissues, however the ingested bacteria must survive firstly gastric acid and secondly the antimicrobial defences of the small intestine. We found that the way in which cells are grown prior to formulation markedly affects sensitivity to acid and bile. Using a previously published stable storage formulation that maintained over 10% viability after 56 days storage at room temperature, we found dried samples of an attenuated S. typhimurium vaccine lost acid and bile resistance compared to the same bacteria taken from fresh culture. The stable formulation utilised osmotic preconditioning in defined medium plus elevated salt concentration to induce intracellular trehalose accumulation before drying. Dried bacteria grown in rich media without osmotic preconditioning showed more resistance to bile, but less stability during storage, suggesting a trade-off between bile resistance and stability. Further optimization is needed to produce the ultimate room-temperature stable oral live bacterial vaccine formulation.

Formulation and delivery of the bacterial antagonist Bacillus subtilis for management of lentil vascular wilt caused by Fusarium oxysporum f. sp lentis

Different formulations of Bacillus subtilis were prepared using standard laboratory protocols. Bacillus subtilis survived in glucose and talc powders at 8.6 and 7.8 log(10) CFU/g, respectively, for 1 year of storage at room temperature compared with 3.5 log(10) CFU/g on a peat formulation. Glasshouse experiments using soil and seed treatments were conducted to test the efficacy of B. subtilis for protecting lentil against the wilt disease caused by Fusariumoxysporum f. sp. lentis. Seed treatments with formulations of B. subtilis on glucose, talc and peat significantly enhanced its biocontrol activity against Fusarium compared with a treatment in which spores were applied directly to seed. The formulations decreased disease severity by reducing colonization of plants by the pathogen, promoting their growth and increased the dry weight of lentil plants. Of these treatments the glucose and talc-based powder formulations were more effective than the peat formulation and the spore application without a carrier. It was shown that the B. subtilis spores applied with glucose were viable for longer than those applied with other carriers. Seed treatment with these formulated spores is an effective delivery system that can provide a conducive environment for B. subtilis to suppress vascular wilt disease on lentil and has the potential for utilization in commercial field application.

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently. Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J = 0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N - 7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N - 7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J > 0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion. (C) 2003 Elsevier B.V. All rights reserved.