Optimal control with stochastic PDE constraints and uncertain controls

The optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed. The Lagrangian that defines the problem is postulated in terms of stochastic functions, with the control function possibly decomposed into an unknown deterministic component and a known zero-mean stochastic component. The extra freedom provided by the stochastic dimension in defining cost functionals is explored, demonstrating the scope for controlling statistical aspects of the system response. One-shot stochastic finite element methods are used to find approximate solutions to control problems. It is shown that applying the stochastic collocation finite element method to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional, thereby forgoing the primary advantage of the collocation method over the stochastic Galerkin method for the considered problem. The application of the presented methods is demonstrated through a number of numerical examples. The presented framework is sufficiently general to also consider a class of inverse problems, and numerical examples of this type are also presented. © 2011 Elsevier B.V.

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.

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.

Physical consequences of a nonparametric uncertainty model in structural dynamics

One of the main claims of the nonparametric model of random uncertainty introduced by Soize (2000) [3] is its ability to account for model uncertainty. The present paper investigates this claim by examining the statistics of natural frequencies, total energy and underlying dispersion equation yielded by the nonparametric approach for two simple systems: a thin plate in bending and a one-dimensional finite periodic massspring chain. Results for the plate show that the average modal density and the underlying dispersion equation of the structure are gradually and systematically altered with increasing uncertainty. The findings for the massspring chain corroborate the findings for the plate and show that the remote coupling of nonadjacent degrees of freedom induced by the approach suppresses the phenomenon of mode localization. This remote coupling also leads to an instantaneous response of all points in the chain when one mass is excited. In the light of these results, it is argued that the nonparametric approach can deal with a certain type of model uncertainty, in this case the presence of unknown terms of higher or lower order in the governing differential equation, but that certain expectations about the system such as the average modal density may conflict with these results. © 2012 Elsevier Ltd.

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.

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.

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.

Optimization algorithms and ODE's in MDO

There is a need for a stronger theoretical understanding of Multidisciplinary Design Optimization (MDO) within the field. Having developed a differential geometry framework in response to this need, we consider how standard optimization algorithms can be modeled using systems of ordinary differential equations (ODEs) while also reviewing optimization algorithms which have been derived from ODE solution methods. We then use some of the framework's tools to show how our resultant systems of ODEs can be analyzed and their behaviour quantitatively evaluated. In doing so, we demonstrate the power and scope of our differential geometry framework, we provide new tools for analyzing MDO systems and their behaviour, and we suggest hitherto neglected optimization methods which may prove particularly useful within the MDO context. Copyright © 2013 by ASME.

instability of standing wave, global existence and blowup for the klein-gordon-zakharov system with different-degree nonlinearities

This paper discusses the Klein–Gordon–Zakharov system with different-degree nonlinearities in two and three space dimensions. Firstly, we prove the existence of standing wave with ground state by applying an intricate variational argument. Next, by introducing an auxiliary functional and an equivalent minimization problem, we obtain two invariant manifolds under the solution flow generated by the Cauchy problem to the aforementioned Klein–Gordon–Zakharov system. Furthermore, by constructing a type of constrained variational problem, utilizing the above two invariant manifolds as well as applying potential well argument and concavity method, we derive a sharp threshold for global existence and blowup. Then, combining the above results, we obtain two conclusions of how small the initial data are for the solution to exist globally by using dilation transformation. Finally, we prove a modified instability of standing wave to the system under study.

Hybrid numerical methods for exponential models of growth

A series of novel numerical methods for the exponential models of growth are proposed. Based on these methods, hybrid predictor-corrector methods are constructed. The hybrid numerical methods can increase the accuracy and the computing speed obviously, as well as enlarge the stability domain greatly. (c) 2005 Published by Elsevier Inc.

A transfer matrix approach to conductance in quantum waveguides

A transfer matrix approach is presented for the study of electron conduction in an arbitrarily shaped cavity structure embedded in a quantum wire. Using the boundary conditions for wave functions, the transfer matrix at an interface with a discontinuous potential boundary is obtained for the first time. The total transfer matrix is calculated by multiplication of the transfer matrix for each segment of the structure as well as numerical integration of coupled second-order differential equations. The proposed method is applied to the evaluation of the conductance and the electron probability density in several typical cavity structures. The effect of the geometrical features on the electron transmission is discussed in detail. In the numerical calculations, the method is found to be more efficient than most of the other methods in the literature and the results are found to be in excellent agreement with those obtained by the recursive Green's function method.