Eliminating discharge imbalances in predicting shallow water flows with shocks

The shallow water equations are widely used in modelling environmental flows. Being a hyperbolic system of differential equations, they admit shocks that represent hydraulic jumps and bores. Although the water surface can be solved satisfactorily with the modern shock-capturing schemes, the predicted flow rate often suffers from imbalances where shocks occur, eg the mass conservation is violated by failing to maintain a constant discharge rate at every cross-section in a steady open channel flow. A total-variation-diminishing Lax-Wendroff scheme is developed, and used to demonstrate how to achieve an exact flux balance. The performance of the proposed methods is inspected through some test cases, which include 1- and 2-dimensional, flat and irregular bed scenarios. The proposed methods are shown to preserve the mass exactly, and can be easily extended to other shock-capturing models.

A perturbation approach to the stabilization of nonlinear cascade systems with time-delay

The effect of bounded input perturbations on the stability of nonlinear globally asymptotically stable delay differential equations is analyzed. We investigate under which conditions global stability is preserved and if not, whether semi-global stabilization is possible by controlling the size or shape of the perturbation. These results are used to study the stabilization of partially linear cascade systems with partial state feedback.

On the application of differential geometry to MDO

Multidisciplinary Design Optimization (MDO) is a methodology for optimizing large coupled systems. Over the years, a number of different MDO decomposition strategies, known as architectures, have been developed, and various pieces of analytical work have been done on MDO and its architectures. However, MDO lacks an overarching paradigm which would unify the field and promote cumulative research. In this paper, we propose a differential geometry framework as such a paradigm: Differential geometry comes with its own set of analysis tools and a long history of use in theoretical physics. We begin by outlining some of the mathematics behind differential geometry and then translate MDO into that framework. This initial work gives new tools and techniques for studying MDO and its architectures while producing a naturally arising measure of design coupling. The framework also suggests several new areas for exploration into and analysis of MDO systems. At this point, analogies with particle dynamics and systems of differential equations look particularly promising for both the wealth of extant background theory that they have and the potential predictive and evaluative power that they hold. © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

A differential lyapunov framework for contraction analysis

Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental stability is inferred from infinitesimal contraction of the Finsler metrics through integration along solutions curves. © 2013 IEEE.

Simulation of three-dimensional viscous flow in turbomachinery geometries using a solution-adaptive unstructured mesh methodology

This paper presents a numerical method for the simulation of flow in turbomachinery blade rows using a solution-adaptive mesh methodology. The fully three-dimensional, compressible, Reynolds-averaged Navier-Stokes equations with k-ε turbulence modeling (and low Reynolds number damping terms) are solved on an unstructured mesh formed from tetrahedral finite volumes. At stages in the solution, mesh refinement is carried out based on flagging cell faces with either a fractional variation of a chosen variable (like Mach number) greater than a given threshold or with a mean value of the chosen variable within a given range. Several solutions are presented, including that for the highly three-dimensional flow associated with the corner stall and secondary flow in a transonic compressor cascade, to demonstrate the potential of the new method.

Building flexible user interfaces for solving PDEs

FEniCS is a collection of software tools for the automated solution of differential equations by finite element methods. In this note, we describe how FEniCS can be used to solve a simple nonlinear model problem with varying levels of automation. At one extreme, FEniCS provides tools for the fully automated and adaptive solution of nonlinear partial differential equations. At the other extreme, FEniCS provides a range of tools that allow the computational scientist to experiment with novel solution algorithms. © 2010 American Institute of Physics.

The post-processed Galerkin method applied to non-linear shell vibrations

We present the results of a computational study of the post-processed Galerkin methods put forward by Garcia-Archilla et al. applied to the non-linear von Karman equations governing the dynamic response of a thin cylindrical panel periodically forced by a transverse point load. We spatially discretize the shell using finite differences to produce a large system of ordinary differential equations (ODEs). By analogy with spectral non-linear Galerkin methods we split this large system into a 'slowly' contracting subsystem and a 'quickly' contracting subsystem. We then compare the accuracy and efficiency of (i) ignoring the dynamics of the 'quick' system (analogous to a traditional spectral Galerkin truncation and sometimes referred to as 'subspace dynamics' in the finite element community when applied to numerical eigenvectors), (ii) slaving the dynamics of the quick system to the slow system during numerical integration (analogous to a non-linear Galerkin method), and (iii) ignoring the influence of the dynamics of the quick system on the evolution of the slow system until we require some output, when we 'lift' the variables from the slow system to the quick using the same slaving rule as in (ii). This corresponds to the post-processing of Garcia-Archilla et al. We find that method (iii) produces essentially the same accuracy as method (ii) but requires only the computational power of method (i) and is thus more efficient than either. In contrast with spectral methods, this type of finite-difference technique can be applied to irregularly shaped domains. We feel that post-processing of this form is a valuable method that can be implemented in computational schemes for a wide variety of partial differential equations (PDEs) of practical importance.

Frictional flow of damp granular material in a conical centrifuge

An analysis is given of velocity and pressure-dependent sliding flow of a thin layer of damp granular material in a spinning cone. Integral momentum equations for steady state, axisymmetric flow are derived using a boundary layer approximation. These reduce to two coupled first-order differential equations for the radial and circumferential sliding velocities. The influence of viscosity and friction coefficients and inlet boundary conditions is explored by presentation of a range of numerical results. In the absence of any interfacial shear traction the flow would, with increasing radial and circumferential slip, follow a trajectory from inlet according to conservation of angular momentum and kinetic energy. Increasing viscosity or friction reduces circumferential slip and, in general, increases the residence time of a particle in the cone. The residence time is practically insensitive to the inlet velocity. However, if the cone angle is very close to the friction angle then the residence time is extremely sensitive to the relative magnitude of these angles. © 2011 Authors.

Development of realistic liquid equations of state for multi-phase CFD modelling

Several equations of state (EOS) have been incorporated into a novel algorithm to solve a system of multi-phase equations in which all phases are assumed to be compressible to varying degrees. The EOSs are used to both supply functional relationships to couple the conservative variables to the primitive variables and to calculate accurately thermodynamic quantities of interest, such as the speed of sound. Each EOS has a defined balance of accuracy, robustness and computational speed; selection of an appropriate EOS is generally problem-dependent. This work employs an AUSM+-up method for accurate discretisation of the convective flux terms with modified low-Mach number dissipation for added robustness of the solver. In this paper we show a newly-developed time-marching formulation for temporal discretisation of the governing equations with incorporated time-dependent source terms, as well as considering the system of eigenvalues that render the governing equations hyperbolic.

Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Reconstruction of arbitrary biochemical reaction networks: A compressive sensing approach

Reconstruction of biochemical reaction networks (BRN) and genetic regulatory networks (GRN) in particular is a central topic in systems biology which raises crucial theoretical challenges in system identification. Nonlinear Ordinary Differential Equations (ODEs) that involve polynomial and rational functions are typically used to model biochemical reaction networks. Such nonlinear models make the problem of determining the connectivity of biochemical networks from time-series experimental data quite difficult. In this paper, we present a network reconstruction algorithm that can deal with ODE model descriptions containing polynomial and rational functions. Rather than identifying the parameters of linear or nonlinear ODEs characterised by pre-defined equation structures, our methodology allows us to determine the nonlinear ODEs structure together with their associated parameters. To solve the network reconstruction problem, we cast it as a compressive sensing (CS) problem and use sparse Bayesian learning (SBL) algorithms as a computationally efficient and robust way to obtain its solution. © 2012 IEEE.

Delayed control of a Moore-Greitzer axial compressor model

Several feedback control laws have appeared in the literature concerning the stabilization of the nonlinear Moore-Greitzer axial compression model. Motivated by magnitude and rate limitations imposed by the physical implementation of the control law, Larsen et al. studied a dynamic implementation of the S-controller suggested by Sepulchre and Kokotović. They showed the potential benefit of implementing the S-controller through a first-order lag: while the location of the closed-loop equilibrium achieved with the static control law was sensitive to poorly known parameters, the dynamic implementation resulted in a small limit cycle at a very desirable location, insensitive to parameter variations. In this paper, we investigate the more general case when the control is applied with a time delay. This can be seen as an extension of the model with a first-order lag. The delay can either be a result of system constraints or be deliberately implemented to achieve better system behavior. The resulting closed-loop system is a set of parameter-dependent delay differential equations. Numerical bifurcation analysis is used to study this model and investigate whether the positive results obtained for the first-order model persist, even for larger values of the delay.

Asymptotic stability for time-variant systems and observability: Uniform and nonuniform criteria

This paper presents some new criteria for uniform and nonuniform asymptotic stability of equilibria for time-variant differential equations and this within a Lyapunov approach. The stability criteria are formulated in terms of certain observability conditions with the output derived from the Lyapunov function. For some classes of systems, this system theoretic interpretation proves to be fruitful since - after establishing the invariance of observability under output injection - this enables us to check the stability criteria on a simpler system. This procedure is illustrated for some classical examples.

Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band-gap edge: An investigation into the role of the density of photon states near resonance

In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e., at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called "mirrorless"cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in two-dimensional photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first-order differential equations, solved using the Green's-function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data. © 2013 American Physical Society.

Spontaneous emission from radiative chiral nematic liquid crystals at the photonic band-gap edge: an investigation into the role of the density of photon states near resonance.

In this article, we investigate the spontaneous emission properties of radiating molecules embedded in a chiral nematic liquid crystal, under the assumption that the electronic transition frequency is close to the photonic edge mode of the structure, i.e., at resonance. We take into account the transition broadening and the decay of electromagnetic field modes supported by the so-called "mirrorless"cavity. We employ the Jaynes-Cummings Hamiltonian to describe the electron interaction with the electromagnetic field, focusing on the mode with the diffracting polarization in the chiral nematic layer. As known in these structures, the density of photon states, calculated via the Wigner method, has distinct peaks on either side of the photonic band gap, which manifests itself as a considerable modification of the emission spectrum. We demonstrate that, near resonance, there are notable differences between the behavior of the density of states and the spontaneous emission profile of these structures. In addition, we examine in some detail the case of the logarithmic peak exhibited in the density of states in two-dimensional photonic structures and obtain analytic relations for the Lamb shift and the broadening of the atomic transition in the emission spectrum. The dynamical behavior of the atom-field system is described by a system of two first-order differential equations, solved using the Green's-function method and the Fourier transform. The emission spectra are then calculated and compared with experimental data.