Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Performance control of distributed generation using digital estimation of signal parameters

The Queensland University of Technology (QUT) allows the presentation of a thesis for the Degree of Doctor of Philosophy in the format of published or submitted papers, where such papers have been published, accepted or submitted during the period of candidature. This thesis is composed of seven published/submitted papers, of which one has been published, three accepted for publication and the other three are under review. This project is financially supported by an Australian Research Council (ARC) Discovery Grant with the aim of proposing strategies for the performance control of Distributed Generation (DG) system with digital estimation of power system signal parameters. Distributed Generation (DG) has been recently introduced as a new concept for the generation of power and the enhancement of conventionally produced electricity. Global warming issue calls for renewable energy resources in electricity production. Distributed generation based on solar energy (photovoltaic and solar thermal), wind, biomass, mini-hydro along with use of fuel cell and micro turbine will gain substantial momentum in the near future. Technically, DG can be a viable solution for the issue of the integration of renewable or non-conventional energy resources. Basically, DG sources can be connected to local power system through power electronic devices, i.e. inverters or ac-ac converters. The interconnection of DG systems to power system as a compensator or a power source with high quality performance is the main aim of this study. Source and load unbalance, load non-linearity, interharmonic distortion, supply voltage distortion, distortion at the point of common coupling in weak source cases, source current power factor, and synchronism of generated currents or voltages are the issues of concern. The interconnection of DG sources shall be carried out by using power electronics switching devices that inject high frequency components rather than the desired current. Also, noise and harmonic distortions can impact the performance of the control strategies. To be able to mitigate the negative effect of high frequency and harmonic as well as noise distortion to achieve satisfactory performance of DG systems, new methods of signal parameter estimation have been proposed in this thesis. These methods are based on processing the digital samples of power system signals. Thus, proposing advanced techniques for the digital estimation of signal parameters and methods for the generation of DG reference currents using the estimates provided is the targeted scope of this thesis. An introduction to this research – including a description of the research problem, the literature review and an account of the research progress linking the research papers – is presented in Chapter 1. One of the main parameters of a power system signal is its frequency. Phasor Measurement (PM) technique is one of the renowned and advanced techniques used for the estimation of power system frequency. Chapter 2 focuses on an in-depth analysis conducted on the PM technique to reveal its strengths and drawbacks. The analysis will be followed by a new technique proposed to enhance the speed of the PM technique while the input signal is free of even-order harmonics. The other techniques proposed in this thesis as the novel ones will be compared with the PM technique comprehensively studied in Chapter 2. An algorithm based on the concept of Kalman filtering is proposed in Chapter 3. The algorithm is intended to estimate signal parameters like amplitude, frequency and phase angle in the online mode. The Kalman filter is modified to operate on the output signal of a Finite Impulse Response (FIR) filter designed by a plain summation. The frequency estimation unit is independent from the Kalman filter and uses the samples refined by the FIR filter. The frequency estimated is given to the Kalman filter to be used in building the transition matrices. The initial settings for the modified Kalman filter are obtained through a trial and error exercise. Another algorithm again based on the concept of Kalman filtering is proposed in Chapter 4 for the estimation of signal parameters. The Kalman filter is also modified to operate on the output signal of the same FIR filter explained above. Nevertheless, the frequency estimation unit, unlike the one proposed in Chapter 3, is not segregated and it interacts with the Kalman filter. The frequency estimated is given to the Kalman filter and other parameters such as the amplitudes and phase angles estimated by the Kalman filter is taken to the frequency estimation unit. Chapter 5 proposes another algorithm based on the concept of Kalman filtering. This time, the state parameters are obtained through matrix arrangements where the noise level is reduced on the sample vector. The purified state vector is used to obtain a new measurement vector for a basic Kalman filter applied. The Kalman filter used has similar structure to a basic Kalman filter except the initial settings are computed through an extensive math-work with regards to the matrix arrangement utilized. Chapter 6 proposes another algorithm based on the concept of Kalman filtering similar to that of Chapter 3. However, this time the initial settings required for the better performance of the modified Kalman filter are calculated instead of being guessed by trial and error exercises. The simulations results for the parameters of signal estimated are enhanced due to the correct settings applied. Moreover, an enhanced Least Error Square (LES) technique is proposed to take on the estimation when a critical transient is detected in the input signal. In fact, some large, sudden changes in the parameters of the signal at these critical transients are not very well tracked by Kalman filtering. However, the proposed LES technique is found to be much faster in tracking these changes. Therefore, an appropriate combination of the LES and modified Kalman filtering is proposed in Chapter 6. Also, this time the ability of the proposed algorithm is verified on the real data obtained from a prototype test object. Chapter 7 proposes the other algorithm based on the concept of Kalman filtering similar to those of Chapter 3 and 6. However, this time an optimal digital filter is designed instead of the simple summation FIR filter. New initial settings for the modified Kalman filter are calculated based on the coefficients of the digital filter applied. Also, the ability of the proposed algorithm is verified on the real data obtained from a prototype test object. Chapter 8 uses the estimation algorithm proposed in Chapter 7 for the interconnection scheme of a DG to power network. Robust estimates of the signal amplitudes and phase angles obtained by the estimation approach are used in the reference generation of the compensation scheme. Several simulation tests provided in this chapter show that the proposed scheme can very well handle the source and load unbalance, load non-linearity, interharmonic distortion, supply voltage distortion, and synchronism of generated currents or voltages. The purposed compensation scheme also prevents distortion in voltage at the point of common coupling in weak source cases, balances the source currents, and makes the supply side power factor a desired value.

Sequential Monte Carlo for Bayesian sequentially designed experiments for discrete data

In this paper we present a sequential Monte Carlo algorithm for Bayesian sequential experimental design applied to generalised non-linear models for discrete data. The approach is computationally convenient in that the information of newly observed data can be incorporated through a simple re-weighting step. We also consider a flexible parametric model for the stimulus-response relationship together with a newly developed hybrid design utility that can produce more robust estimates of the target stimulus in the presence of substantial model and parameter uncertainty. The algorithm is applied to hypothetical clinical trial or bioassay scenarios. In the discussion, potential generalisations of the algorithm are suggested to possibly extend its applicability to a wide variety of scenarios

Frank-starling control of a left ventricular assist device

A physiological control system was developed for a rotary left ventricular assist device (LVAD) in which the target pump flow rate (LVADQ) was set as a function of left atrial pressure (LAP), mimicking the Frank-Starling mechanism. The control strategy was implemented using linear PID control and was evaluated in a pulsatile mock circulation loop using a prototyped centrifugal pump by varying pulmonary vascular resistance to alter venous return. The control strategy automatically varied pump speed (2460 to 1740 to 2700 RPM) in response to a decrease and subsequent increase in venous return. In contrast, a fixed-speed pump caused a simulated ventricular suction event during low venous return and higher ventricular volumes during high venous return. The preload sensitivity was increased from 0.011 L/min/mmHg in fixed speed mode to 0.47L/min/mmHg, a value similar to that of the native healthy heart. The sensitivity varied automatically to maintain the LAP and LVADQ within a predefined zone. This control strategy requires the implantation of a pressure sensor in the left atrium and a flow sensor around the outflow cannula of the LVAD. However, appropriate pressure sensor technology is not yet commercially available and so an alternative measure of preload such as pulsatility of pump signals should be investigated.

A novel SNP in COMT is associated with alcohol dependence but not opiate or nicotine dependence: a case control study

Background It is well established that COMT is a strong candidate gene for substance use disorder and schizophrenia. Recently we identified two SNPs in COMT (rs4680 and rs165774) that are associated with schizophrenia in an Australian cohort. Individuals with schizophrenia were more than twice as likely to carry the GG genotype compared to the AA genotype for both the rs165774 and rs4680 SNPs. Association of both rs4680 and rs165774 with substance dependence, a common comorbidity of schizophrenia has not been investigated. Methods To determine whether COMT is important in substance dependence, rs165774 and rs4680 were genotyped and haplotyped in patients with nicotine, alcohol and opiate dependence. Results The rs165774 SNP was associated with alcohol dependence. However, it was not associated with nicotine or opiate dependence. Individuals with alcohol dependence were more than twice as likely to carry the GG or AG genotypes compared to the AA genotype, indicating a dominant mode of inheritance. The rs4680 SNP showed a weak association with alcohol dependence at the allele level that did not reach significance at the genotype level but it was not associated with nicotine or opiate dependence. Analysis of rs165774/rs4680 haplotypes also revealed association with alcohol dependence with the G/G haplotype being almost 1.5 times more common in alcohol-dependent cases. Conclusions Our study provides further support for the importance of the COMT in alcohol dependence in addition to schizophrenia. It is possible that the rs165774 SNP, in combination with rs4680, results in a common molecular variant of COMT that contributes to schizophrenia and alcohol dependence susceptibility. This is potentially important for future studies of comorbidity. As our participant numbers are limited our observations should be viewed with caution until they are independently replicated.

Design and fabrication of tubular scaffolds via direct writing in a melt electrospinning mode

Flexible tubular structures fabricated from solution electrospun fibers are finding increasing use in tissue engineering applications. However it is difficult to control the deposition of fibers due to the chaotic nature of the solution electrospinning jet. By using non-conductive polymer melts instead of polymer solutions the path and collection of the fiber becomes predictable. In this work we demonstrate the melt electrospinning of polycaprolactone in a direct writing mode onto a rotating cylinder. This allows the design and fabrication of tubes using 20 μm diameter fibers with controllable micropatterns and mechanical properties. A key design parameter is the fiber winding angle, where it allows control over scaffold pore morphology (e.g. size, shape, number and porosity). Furthermore, the establishment of a finite element model as a predictive design tool is validated against mechanical testing results of melt electrospun tubes to show that a lesser winding angle provides improved mechanical response to uniaxial tension and compression. In addition, we show that melt electrospun tubes support the growth of three different cell types in vitro and are therefore promising scaffolds for tissue engineering applications.

Practical stability of approximating discrete-time filters with respect to model mismatch

This paper establishes practical stability results for an important range of approximate discrete-time filtering problems involving mismatch between the true system and the approximating filter model. Practical stability is established in the sense of an asymptotic bound on the amount of bias introduced by the model approximation. Our analysis applies to a wide range of estimation problems and justifies the common practice of approximating intractable infinite dimensional nonlinear filters by simpler computationally tractable filters.

Surface effects on the dual-mode vibration of <110> silver nanowires with different cross-sections

Dual-mode vibration of nanowires has been reported experimentally through actuation of the nanowire at its resonance frequency, which is expected to open up a variety of new modalities for the NEMS that could operate in the nonlinear regime. In the present work, we utilize large scale molecular dynamics simulations to investigate the dual-mode vibration of <110> Ag nanowires with triangular, rhombic and truncated rhombic cross-sections. By incorporating the generalized Young-Laplace equation into Euler-Bernoulli beam theory, the influence of surface effects on the dual-mode vibration is studied. Due to the different lattice spacing in principal axes of inertia of the {110} atomic layers, the NW is also modeled as a discrete system to reveal the influence from such specific atomic arrangement. It is found that the <110> Ag NW will under a dual-mode vibration if the actuation direction is deviated from the two principal axes of inertia. The predictions of the two first mode natural frequencies by the classical beam model appear underestimated comparing with the MD results, which are found to be enhanced by the discrete model. Particularly, the predictions by the beam theory with the contribution of surface effects are uniformly larger than the classical beam model, which exhibit better agreement with MD results for larger cross-sectional size. However, for ultrathin NWs, current consideration of surface effects is still experiencing certain inaccuracy. In all, for all different cross-sections, the inclusion of surface effects is found to reduce the difference between the two first mode natural frequencies. This trend is observed consistent with MD results. This study provides a first comprehensive investigation on the dual-mode vibration of <110> oriented Ag NWs, which is supposed to benefit the applications of NWs that acting as a resonating beam.

Budget-aware role based access control

The suitability of Role Based Access Control (RBAC) is being challenged in dynamic environments like healthcare. In an RBAC system, a user's legitimate access may be denied if their need has not been anticipated by the security administrator at the time of policy specification. Alternatively, even when the policy is correctly specified an authorised user may accidentally or intentionally misuse the granted permission. The heart of the challenge is the intrinsic unpredictability of users' operational needs as well as their incentives to misuse permissions. In this paper we propose a novel Budget-aware Role Based Access Control (B-RBAC) model that extends RBAC with the explicit notion of budget and cost, where users are assigned a limited budget through which they pay for the cost of permissions they need. We propose a model where the value of resources are explicitly defined and an RBAC policy is used as a reference point to discriminate the price of access permissions, as opposed to representing hard and fast rules for making access decisions. This approach has several desirable properties. It enables users to acquire unassigned permissions if they deem them necessary. However, users misuse capability is always bounded by their allocated budget and is further adjustable through the discrimination of permission prices. Finally, it provides a uniform mechanism for the detection and prevention of misuses.

Advanced battery storage control for an autonomous microgrid

A new control method for battery storage to maintain acceptable voltage profile in autonomous microgrids is proposed in this article. The proposed battery control ensures that the bus voltages in the microgrid are maintained during disturbances such as load change, loss of micro-sources, or distributed generations hitting power limit. Unlike the conventional storage control based on local measurements, the proposed method is based on an advanced control technique, where the reference power is determined based on the voltage drop profile at the battery bus. An artificial neural network based controller is used to determine the reference power needed for the battery to hold the microgrid voltage within regulation limits. The pattern of drop in the local bus voltage during power imbalance is used to train the controller off-line. During normal operation, the battery floats with the local bus voltage without any power injection. The battery is charged or discharged during the transients with a high gain feedback loop. Depending on the rate of voltage fall, it is switched to power control mode to inject the reference power determined by the proposed controller. After a defined time period, the battery power injection is reduced to zero using slow reverse-droop characteristics, ensuring a slow rate of increase in power demand from the other distributed generations. The proposed control method is simulated for various operating conditions in a microgrid with both inertial and converter interfaced sources. The proposed battery control provides a quick load pick up and smooth load sharing with the other micro-sources in a disturbance. With various disturbances, maximum voltage drop over 8% with conventional energy storage is reduced within 2.5% with the proposed control method.

Adaptive stepsize based on control theory for stochastic differential equations

The numerical solution of stochastic differential equations (SDEs) has been focused recently on the development of numerical methods with good stability and order properties. These numerical implementations have been made with fixed stepsize, but there are many situations when a fixed stepsize is not appropriate. In the numerical solution of ordinary differential equations, much work has been carried out on developing robust implementation techniques using variable stepsize. It has been necessary, in the deterministic case, to consider the "best" choice for an initial stepsize, as well as developing effective strategies for stepsize control-the same, of course, must be carried out in the stochastic case. In this paper, proportional integral (PI) control is applied to a variable stepsize implementation of an embedded pair of stochastic Runge-Kutta methods used to obtain numerical solutions of nonstiff SDEs. For stiff SDEs, the embedded pair of the balanced Milstein and balanced implicit method is implemented in variable stepsize mode using a predictive controller for the stepsize change. The extension of these stepsize controllers from a digital filter theory point of view via PI with derivative (PID) control will also be implemented. The implementations show the improvement in efficiency that can be attained when using these control theory approaches compared with the regular stepsize change strategy.

A class B switch-mode assisted linear amplifier

A switch-mode assisted linear amplifier (SMALA) combining a linear (Class B) and a switch-mode (Class D) amplifier is presented. The usual single hysteretic controlled half-bridge current dumping stage is replaced by two parallel buck converter stages, in a parallel voltage controlled topology. These operate independently: one buck converter sources current to assist the upper Class B output device, and a complementary converter sinks current to assist the lower device. This topology lends itself to a novel control approach of a dead-band at low power levels where neither class D amplifier assists, allowing the class B amplifier to supply the load without interference, ensuring high fidelity. A 20 W implementation demonstrates 85% efficiency, with distortion below 0.08% measured across the full audio bandwidth at 15 W. The class D amplifier begins assisting at 2 W, and below this value, the distortion was below 0.03%. Complete circuitry is given, showing the simplicity of the additional class D amplifier and its corresponding control circuitry.

The molecular structure of the borate mineral inderite Mg(H4B3O7)(OH)⋅5H2O – A vibrational spectroscopic study

We have undertaken a study of the mineral inderite Mg(H4B3O7)(OH)⋅5H2O a hydrated hydroxy borate mineral of magnesium using scanning electron microscopy, thermogravimetry and vibrational spectroscopic techniques. The structure consists of [B3O3(OH)5]2-[B3O3(OH)5]2- soroborate groups and Mg(OH)2(H2O)4 octahedra interconnected into discrete molecules by the sharing of two OH groups. Thermogravimetry shows a mass loss of 47.2% at 137.5 °C, proving the mineral is thermally unstable. Raman bands at 954, 1047 and 1116 cm−1 are assigned to the trigonal symmetric stretching mode. The two bands at 880 and 916 cm−1 are attributed to the symmetric stretching mode of the tetrahedral boron. Both the Raman and infrared spectra of inderite show complexity. Raman bands are observed at 3052, 3233, 3330, 3392 attributed to water stretching vibrations and 3459 cm−1 with sharper bands at 3459, 3530 and 3562 cm−1 assigned to OH stretching vibrations. Vibrational spectroscopy is used to assess the molecular structure of inderite.

Residential dissonance and mode choice

Residential dissonance refers to the mismatch in land-use patterns between individuals’ preferred residential neighbourhood type and the type of neighbourhood in which they currently reside. Current knowledge regarding the impact of residential dissonance is limited to short-term travel behaviours in urban vs. suburban, and rural vs. urban areas. Although the prevailing view is that dissonants adjust their orientation and lifestyle around their surrounding land use over time, empirical evidence is lacking to support this proposition. This research identifies both short-term mode choice behaviour and medium-term mode shift behaviour of dissonants in transit oriented development (TODs) vs. non-TOD areas in Brisbane, Australia. Natural groupings of neighbourhood profiles (e.g. residential density, land use diversity, intersection density, cul-de-sac density, and public transport accessibility levels) of 3957 individuals were identified as living either in a TOD (510 individuals) or non-TOD (3447 individuals) areas in Brisbane using the TwoStep cluster analysis technique. Levels of dissonance were measured based on a factor analysis of 16 items representing the travel attitudes/preferences of individuals. Two multinomial logistic (MNL) regression models were estimated to understand mode choice behaviour of (1) TOD dissonants, and (2) non-TOD dissonants in 2009, controlling for socio-demographics and environmental characteristics. Two additional MNL regression models were estimated to investigate mode shift behaviour of (3) TOD dissonants, and (4) non-TOD dissonants between 2009 and 2011, also controlling for socio-demographic, changes in socio-demographic, and built environmental factors. The findings suggest that travel preference is relatively more influential in transport mode choice decisions compared with built environment features. Little behavioural evidence was found to support the adjustment of a dissonant orientation toward a particular land use feature and mode accessibility they represent (e.g. a modal shift to greater use of the car for non-TOD dissonants). TOD policies should focus on reducing the level of dissonance in TODs in order to enhance transit ridership.

Improved catalyst, particularly for emission control

A catalyst comprising a catalytic material supported on a support, characterized in that the support comprises particles predominantly having a max. particle size of less than 1000 nm and an aspect ratio of greater than, and the. catalytic material is mainly present in the form of discrete islands of catalytic material supported on the support, with a substantial proportion of the islands of catalytic material being sep. and isolated from other islands of catalytic material. The islands of catalytic material are sep. and isolated from other islands of catalytic material such that diffusion and growth of the islands of catalytic material at elevated temp. is minimized or avoided. The disclosure and examples pertain to emission control catalysts. [on SciFinder(R)]