Holistic finite differences accurately model the dynamics of the Kuramoto-Sivashinsky equation

We analyse the nonlinear Kuramoto-Sivashinsky equation to develop an accurate finite difference approximation to its dynamics. The analysis is based upon centre manifold theory so we are assured that the finite difference model accurately models the dynamics and may be constructed systematically. The theory is applied after dividing the physical domain into small elements by introducing insulating internal boundaries which are later removed. The Kuramoto-Sivashinsky equation is used as an example to show how holistic finite differences may be applied to fourth order, nonlinear, spatio-temporal dynamical systems. This novel centre manifold approach is holistic in the sense that it treats the dynamical equations as a whole, not just as the sum of separate terms

Accurate macroscale modelling of spatial dynamics in multiple dimensions

Developments in dynamical systems theory provides new support for the macroscale modelling of pdes and other microscale systems such as Lattice Boltzmann, Monte Carlo or Molecular Dynamics simulators. By systematically resolving subgrid microscale dynamics the dynamical systems approach constructs accurate closures of macroscale discretisations of the microscale system. Here we specifically explore reaction-diffusion problems in two spatial dimensions as a prototype of generic systems in multiple dimensions. Our approach unifies into one the modelling of systems by a type of finite elements, and the `equation free' macroscale modelling of microscale simulators efficiently executing only on small patches of the spatial domain. Centre manifold theory ensures that a closed model exist on the macroscale grid, is emergent, and is systematically approximated. Dividing space either into overlapping finite elements or into spatially separated small patches, the specially crafted inter-element/patch coupling also ensures that the constructed discretisations are consistent with the microscale system/PDE to as high an order as desired. Computer algebra handles the considerable algebraic details as seen in the specific application to the Ginzburg--Landau PDE. However, higher order models in multiple dimensions require a mixed numerical and algebraic approach that is also developed. The modelling here may be straightforwardly adapted to a wide class of reaction-diffusion PDEs and lattice equations in multiple space dimensions. When applied to patches of microscopic simulations our coupling conditions promise efficient macroscale simulation.Comment: some figures with 3D interaction when viewed in Acrobat Reader. arXiv admin note: substantial text overlap with arXiv:0904.085

Accurately model the Kuramoto--Sivashinsky dynamics with holistic discretisation

We analyse the nonlinear Kuramoto--Sivashinsky equation to develop accurate discretisations modeling its dynamics on coarse grids. The analysis is based upon centre manifold theory so we are assured that the discretisation accurately models the dynamics and may be constructed systematically. The theory is applied after dividing the physical domain into small elements by introducing isolating internal boundaries which are later removed. Comprehensive numerical solutions and simulations show that the holistic discretisations excellently reproduce the steady states and the dynamics of the Kuramoto--Sivashinsky equation. The Kuramoto--Sivashinsky equation is used as an example to show how holistic discretisation may be successfully applied to fourth order, nonlinear, spatio-temporal dynamical systems. This novel centre manifold approach is holistic in the sense that it treats the dynamical equations as a whole, not just as the sum of separate terms.Comment: Without figures. See http://www.sci.usq.edu.au/staff/aroberts/ksdoc.pdf to download a version with the figure

Metamodelling of multivariable engine models for real-time flight simulation.

Sophisticated real-time distributed flight simulation environments may be constructed from a wide range of modelling and simulation tools. In this way accuracy, detail and model flexibility may be incorporated into the simulator. Distributed components may be constructed by a wide range of methods, from high level environments such as Matlab, through coded environments such as C or Fortran to hardware-in-the- loop. In this paper the Response Surface Methodology is combined with a hyper-heuristic (evolutionary algorithm) and applied to the representation of computationally intensive non-linear multivariable engine modelling. The paper investigates the potential for metamodelling (models of models) dynamic models which were previously too slow to be included in multi-component, high resolution real-time simulation environments. A multi-dimensional gas turbine model with five primary control inputs, six environmental inputs and eleven outputs is considered. An investigation has been conducted to ascertain to what extent these systems can be approximated by response surfaces with experiments which have been designed by hyper-heuristics as a first step towards automatic modelling methodology

Neighbour coverage: a dynamic probabilistic route discovery for mobile ad hoc networks

Blind flooding is extensively use in ad hoc routing protocols for on-demand route discovery, where a mobile node blindly rebroadcasts received route request (RREQ) packets until a route to a particular destination is established. This can potentially lead to high channel contention, causing redundant retransmissions and thus excessive packet collisions in the network. Such a phenomenon induces what is known as broadcast storm problem, which has been shown to greatly increase the network communication overhead and end-to-end delay. In this paper, we show that the deleterious impact of such a problem can be reduced if measures are taken during the dissemination of RREQ packets. We propose a generic probabilistic method for route discovery, that is simple to implement and can significantly reduce the overhead associated with the dissemination of RREQs. Our analysis reveals that equipping AODV with probabilistic route discovery can result in significant reduction of routing control overhead while achieving good throughput

Lattice QCD Production on Commodity Clusters at Fermilab

We describe the construction and results to date of Fermilab's three Myrinet-networked lattice QCD production clusters (an 80-node dual Pentium III cluster, a 48-node dual Xeon cluster, and a 128-node dual Xeon cluster). We examine a number of aspects of performance of the MILC lattice QCD code running on these clusters.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 6 pages, LaTeX, 8 eps figures. PSN TUIT00

Photoemission Quasi-Particle Spectra of Sr2_2RuO4_4

Multi-band quasi-particle calculations based on perturbation theory and dynamical mean field methods show that the creation of a photoemission hole state in Sr$_2$RuO$_4$ is associated with a highly anisotropic self-energy. Since the narrow Ru-derived $d_{xz,yz}$ bands are more strongly distorted by Coulomb correlations than the wide $d_{xy}$ band, charge is partially transferred from $d_{xz,yz}$ to $d_{xy}$, thereby shifting the $d_{xy}$ van Hove singularity close to the Fermi level.Comment: 4 pages, to be published in PRB Rapid Com

Development of a contra-rotating tidal current turbine and analysis of performance

A contra-rotating marine current turbine has a number of attractive features: nearzero reactive torque on the support structure, near-zero swirl in the wake, and high relative inter-rotor rotational speeds. Modified blade element modelling theory has been used to design and predict the characteristics of such a turbine, and a model turbine and test rig have been constructed. Tests in a towing tank demonstrated the feasibility of the concept. Power coefficients were high for such a small model and in excellent agreement with predictions, confirming the accuracy of the computational modelling procedures. Highfrequency blade loading data were obtained in the course of the experiments. These show the anticipated dynamic components for a contra-rotating machine. Flow visualization of the wake verified the lack of swirl behind the turbine. A larger machine is presently under construction for sea trials

Validity of the Adiabatic Approximation

We analyze the validity of the adiabatic approximation, and in particular the reliability of what has been called the "standard criterion" for validity of this approximation. Recently, this criterion has been found to be insufficient. We will argue that the criterion is sufficient only when it agrees with the intuitive notion of slowness of evolution of the Hamiltonian. However, it can be insufficient in cases where the Hamiltonian varies rapidly but only by a small amount. We also emphasize the distinction between the adiabatic {\em theorem} and the adiabatic {\em approximation}, two quite different although closely related ideas.Comment: 4 pages, 1 figur