Non-linear dispersive waves with a small dissipation

The general theory of Whitham for slowly-varying non-linear wavetrains is extended to the case where some of the defining partial differential equations cannot be put into conservation form. Typical examples are considered in plasma dynamics and water waves in which the lack of a conservation form is due to dissipation; an additional non-conservative element, the presence of an external force, is treated for the plasma dynamics example. Certain numerical solutions of the water waves problem (the Korteweg-de Vries equation with dissipation) are considered and compared with perturbation expansions about the linearized solution; it is found that the first correction term in the perturbation expansion is an excellent qualitative indicator of the deviation of the dissipative decay rate from linearity.

A method for deriving necessary and sufficient conditions for the existence of a general uniform wavetrain solution is presented and illustrated in the plasma dynamics problem. Peaking of the plasma wave is demonstrated, and it is shown that the necessary and sufficient existence conditions are essentially equivalent to the statement that no wave may have an amplitude larger than the peaked wave.

A new type of fully non-linear stability criterion is developed for the plasma uniform wavetrain. It is shown explicitly that this wavetrain is stable in the near-linear limit. The nature of this new type of stability is discussed.

Steady shock solutions are also considered. By a quite general method, it is demonstrated that the plasma equations studied here have no steady shock solutions whatsoever. A special type of steady shock is proposed, in which a uniform wavetrain joins across a jump discontinuity to a constant state. Such shocks may indeed exist for the Korteweg-de Vries equation, but are barred from the plasma problem because entropy would decrease across the shock front.

Finally, a way of including the Landau damping mechanism in the plasma equations is given. It involves putting in a dissipation term of convolution integral form, and parallels a similar approach of Whitham in water wave theory. An important application of this would be towards resolving long-standing difficulties about the "collisionless" shock.

A perturbation procedure for nonlinear oscillations (The dynamics of two oscillators with weak nonlinear coupling)

This thesis considers in detail the dynamics of two oscillators with weak nonlinear coupling. There are three classes of such problems: non-resonant, where the Poincaré procedure is valid to the order considered; weakly resonant, where the Poincaré procedure breaks down because small divisors appear (but do not affect the O(1) term) and strongly resonant, where small divisors appear and lead to O(1) corrections. A perturbation method based on Cole's two-timing procedure is introduced. It avoids the small divisor problem in a straightforward manner, gives accurate answers which are valid for long times, and appears capable of handling all three types of problems with no change in the basic approach.

One example of each type is studied with the aid of this procedure: for the nonresonant case the answer is equivalent to the Poincaré result; for the weakly resonant case the analytic form of the answer is found to depend (smoothly) on the difference between the initial energies of the two oscillators; for the strongly resonant case we find that the amplitudes of the two oscillators vary slowly with time as elliptic functions of ϵ t, where ϵ is the (small) coupling parameter.

Our results suggest that, as one might expect, the dynamical behavior of such systems varies smoothly with changes in the ratio of the fundamental frequencies of the two oscillators. Thus the pathological behavior of Whittaker's adelphic integrals as the frequency ratio is varied appears to be due to the fact that Whittaker ignored the small divisor problem. The energy sharing properties of these systems appear to depend strongly on the initial conditions, so that the systems not ergodic.

The perturbation procedure appears to be applicable to a wide variety of other problems in addition to those considered here.

A singular perturbation method of calculating the behavior of supercavitating hydrofoils with rounded noses

A simple, direct and accurate method to predict the pressure distribution on supercavitating hydrofoils with rounded noses is presented. The thickness of body and cavity is assumed to be small. The method adopted in the present work is that of singular perturbation theory. Far from the leading edge linearized free streamline theory is applied. Near the leading edge, however, where singularities of the linearized theory occur, a non-linear local solution is employed. The two unknown parameters which characterize this local solution are determined by a matching procedure. A uniformly valid solution is then constructed with the aid of the singular perturbation approach.

The present work is divided into two parts. In Part I isolated supercavitating hydrofoils of arbitrary profile shape with parabolic noses are investigated by the present method and its results are compared with the new computational results made with Wu and Wang's exact "functional iterative" method. The agreement is very good. In Part II this method is applied to a linear cascade of such hydrofoils with elliptic noses. A number of cases are worked out over a range of cascade parameters from which a good idea of the behavior of this type of important flow configuration is obtained.

Some of the computational aspects of Wu and Wang's functional iterative method heretofore not successfully applied to this type of problem are described in an appendix.

Loss comparison between an SPWM and harmonic elimination excited small, (<1kW) induction motor drive using Pspice simulation and calorimetry

Adopting square wave excitation to drive induction motors (IMs) can substantially reduce inverter switching losses. However, the low-order time harmonics inherent in the output voltage generates parasitic torques that degrade motor performance and reduce efficiency. In this paper, a novel harmonic elimination modulation technique with full voltage control is studied as an interesting and alternative means of operating small (<1kW) IM drives efficiently. A fully verified harmonic elimination scheme, which removes the 5th, 7th, 11th, 13th and 17 th time harmonics, was implemented and applied to an IGBT driven IM. The power losses incurred in the inverter and the IM as a result of the switching scheme have been determined. © 2008 Crown copyright.

Electrostatic mode envelope excitations in e-p-i plasmas - application in warm pair ion plasmas with a small fraction of stationary ions

The nonlinear propagation of amplitude-modulated electrostatic wavepackets in an electron-positron-ion (e-p-i) plasma is considered, by employing a two-fluid plasma model. Considering propagation parallel to the external magnetic field, two distinct electrostatic modes are obtained, namely a quasi-thermal acoustic-like lower mode and a Langmuir-like optic-type upper one. These results equally apply in warm pair ion ( e. g. fullerene) plasmas contaminated by a small fraction of stationary ions ( or dust), in agreement with experimental observations and theoretical predictions in pair plasmas. Considering small yet weakly nonlinear deviations from equilibrium, and adopting a multiple-scales perturbation technique, the basic set of model equations is reduced to a nonlinear Schrodinger (NLS) equation for the slowly varying electric field perturbation amplitude. The analysis reveals that the lower ( acoustic) mode is mostly stable for large wavelengths, and may propagate in the form of a dark-type envelope soliton ( a void) modulating a carrier wavepacket, while the upper linear mode is intrinsically unstable, and thus favours the formation of bright-type envelope soliton ( pulse) modulated wavepackets. The stability ( instability) range for the acoustic ( Langmuir-like optic) mode shifts to larger wavenumbers as the positive-to-negative ion temperature ( density) ratio increases. These results may be of relevance in astrophysical contexts, where e-p-i plasmas are encountered, and may also serve as prediction of the behaviour of doped ( or dust-contaminated) fullerene plasmas, in the laboratory.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.

Evaluating the Small Population Paradigm for Rare Large-Bodied Woodpeckers, with Implications for the Ivory-billed Woodpecker

Six large-bodied, ≥ 120 g, woodpecker species are listed as near-threatened to critically endangered by the International Union for Conservation of Nature (IUCN). The small population paradigm assumes that these populations are likely to become extinct without an increase in numbers, but the combined influences of initial population size and demographic rates, i.e., annual adult survival and fecundity, may drive population persistence for these species. We applied a stochastic, stage-based single-population model to available demographic rates for Dryocopus and Campephilus woodpeckers. In particular, we determined the change in predicted extinction rate, i.e., proportion of simulated populations that went extinct within 100 yr, to concomitant changes in six input parameters. To our knowledge, this is the first study to evaluate the combined importance of initial population size and demographic rates for the persistence of large-bodied woodpeckers. Under a worse-case scenario, the median time to extinction was 7 yr (range: 1–32). Across the combinations of other input values, increasing initial population size by one female induced, on average, 0.4%–3.2% (range: 0%–28%) reduction in extinction rate. Increasing initial population size from 5–30 resulted in extinction rates < 0.05 under limited conditions: (1) all input values were intermediate, or (2) Allee effect present and annual adult survival ≥ 0.8. Based on our model, these species can persist as rare, as few as five females, and thus difficult-to-detect, populations provided they maintain ≥ 1.1 recruited females annually per adult female and an annual adult survival rate ≥ 0.8. Athough a demographic-based population viability analysis (PVA) is useful to predict how extinction rate changes across scenarios for life-history attributes, the next step for modeling these populations should incorporate more easily acquired data on changes in patch occupancy to make predictions about patch colonization and extinction rates.

Perturbation growth at the convective scale for CSIP IOP18

The Met Office Unified Model is run for a case observed during Intensive Observation Period 18 (IOP18) of the Convective Storms Initiation Project (CSIP). The aims are to identify the physical processes that lead to perturbation growth at the convective scale in response to model-state perturbations and to determine their sensitivity to the character of the perturbations. The case is strongly upper-level forced but with detailed mesoscale/convective-scale evolution that is dependent on smaller-scale processes. Potential temperature is perturbed within the boundary layer. The effects on perturbation growth of both the amplitude and typical scalelength of the perturbations are investigated and perturbations are applied either sequentially (every 30 min throughout the simulation) or at specific times. The direct effects (within one timestep) of the perturbations are to generate propagating Lamb and acoustic waves and produce generally small changes in cloud parameters and convective instability. In exceptional cases a perturbation at a specific gridpoint leads to switching of the diagnosed boundary-layer type or discontinuous changes in convective instability, through the generation or removal of a lid. The indirect effects (during the entire simulation) are changes in the intensity and location of precipitation and in the cloud size distribution. Qualitatively different behaviour is found for strong (1K amplitude) and weak (0.01K amplitude) perturbations, with faster growth after sunrise found only for the weaker perturbations. However, the overall perturbation growth (as measured by the root-mean-square error of accumulated precipitation) reaches similar values at saturation, regardless of the perturbation characterisation.

The effects of environmental perturbation and measurement error on estimates of the shape parameter in the theta-logistic model of population regulation

The theta-logistic is a widely used generalisation of the logistic model of regulated biological processes which is used in particular to model population regulation. Then the parameter theta gives the shape of the relationship between per-capita population growth rate and population size. Estimation of theta from population counts is however subject to bias, particularly when there are measurement errors. Here we identify factors disposing towards accurate estimation of theta by simulation of populations regulated according to the theta-logistic model. Factors investigated were measurement error, environmental perturbation and length of time series. Large measurement errors bias estimates of theta towards zero. Where estimated theta is close to zero, the estimated annual return rate may help resolve whether this is due to bias. Environmental perturbations help yield unbiased estimates of theta. Where environmental perturbations are large, estimates of theta are likely to be reliable even when measurement errors are also large. By contrast where the environment is relatively constant, unbiased estimates of theta can only be obtained if populations are counted precisely Our results have practical conclusions for the design of long-term population surveys. Estimation of the precision of population counts would be valuable, and could be achieved in practice by repeating counts in at least some years. Increasing the length of time series beyond ten or 20 years yields only small benefits. if populations are measured with appropriate accuracy, given the level of environmental perturbation, unbiased estimates can be obtained from relatively short censuses. These conclusions are optimistic for estimation of theta. (C) 2008 Elsevier B.V All rights reserved.

On the numerical integration of a class of singular perturbation problems

A three-point difference scheme recently proposed in Ref. 1 for the numerical solution of a class of linear, singularly perturbed, two-point boundary-value problems is investigated. The scheme is derived from a first-order approximation to the original problem with a small deviating argument. It is shown here that, in the limit, as the deviating argument tends to zero, the difference scheme converges to a one-sided approximation to the original singularly perturbed equation in conservation form. The limiting scheme is shown to be stable on any uniform grid. Therefore, no advantage arises from using the deviating argument, and the most accurate and efficient results are obtained with the deviation at its zero limit.

The effect of a single supernova explosion on the cuspy density profile of a small-mass dark matter halo

Some observations of galaxies, and in particular dwarf galaxies, indicate a presence of cored density profiles in apparent contradiction with cusp profiles predicted by dark matter N-body simulations. We constructed an analytical model, using particle distribution functions (DFs), to show how a supernova (SN) explosion can transform a cusp density profile in a small-mass dark matter halo into a cored one. Considering the fact that an SN efficiently removes matter from the centre of the first haloes, we study the effect of mass removal through an SN perturbation in the DFs. We find that the transformation from a cusp into a cored profile occurs even for changes as small as 0.5 per cent of the total energy of the halo, which can be produced by the expulsion of matter caused by a single SN explosion.

Stability of the curvature perturbation in dark sectors` mutual interacting models

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations: however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density. (C) 2008 Elsevier B.V. All rights reserved.

Small-signal stability assessment of active distribution networks with dynamic loads

This paper investigates small-signal stability of a distribution system with distributed generator and induction motor load, as a dynamic element. The analysis is carried out over a distribution test system with different types of induction motor loads. The system is linearised by the perturbation method. Eigenvalues and participation factors are calculated to see the modal interaction of the system. The study indicates that load voltage dynamics significantly influence the damping of a newly identified voltage mode. This mode has frequency of oscillation between the electromechanical and subsynchronous oscillation of power systems. To justify the validity of the modal analysis time domain simulation is also carried out. Finally, significant parameters of the system that affect the damping and frequency of the oscillation are identified.

Time-periodic perturbation of a Lienard equation with an unbounded homoclinic loop

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

LARGE-ORDER PERTURBATION EXPANSIONS FOR THE CHARGED HARMONIC-OSCILLATOR

The charged oscillator, defined by the Hamiltonian H = -d2/dr2+ r2 + lambda/r in the domain [0, infinity], is a particular case of the family of spiked oscillators, which does not behave as a supersingular Hamiltonian. This problem is analysed around the three regions lambda --> infinity, lambda --> 0 and lambda --> -infinity by using Rayleigh-Ritz large-order perturbative expansions. A path is found to connect the large lambda regions with the small lambda region by means of the renormalization of the series expansions in lambda. Finally, the Riccati-Pade method is used to construct an implicit expansion around lambda --> 0 which extends to very large values of Absolute value of lambda.