Optimisation of Kinetic Parameters for an After-treatment Catalyst

Mathematical modelling has become an essential tool in the design of modern catalytic systems. Emissions legislation is becoming increasingly stringent, and so mathematical models of aftertreatment systems must become more accurate in order to provide confidence that a catalyst will convert pollutants over the required range of conditions. 
Automotive catalytic converter models contain several sub-models that represent processes such as mass and heat transfer, and the rates at which the reactions proceed on the surface of the precious metal. Of these sub-models, the prediction of the surface reaction rates is by far the most challenging due to the complexity of the reaction system and the large number of gas species involved. The reaction rate sub-model uses global reaction kinetics to describe the surface reaction rate of the gas species and is based on the Langmuir Hinshelwood equation further developed by Voltz et al. [1] The reactions can be modelled using the pre-exponential and activation energies of the Arrhenius equations and the inhibition terms. 
The reaction kinetic parameters of aftertreatment models are found from experimental data, where a measured light-off curve is compared against a predicted curve produced by a mathematical model. The kinetic parameters are usually manually tuned to minimize the error between the measured and predicted data. This process is most commonly long, laborious and prone to misinterpretation due to the large number of parameters and the risk of multiple sets of parameters giving acceptable fits. Moreover, the number of coefficients increases greatly with the number of reactions. Therefore, with the growing number of reactions, the task of manually tuning the coefficients is becoming increasingly challenging. 
In the presented work, the authors have developed and implemented a multi-objective genetic algorithm to automatically optimize reaction parameters in AxiSuite®, [2] a commercial aftertreatment model. The genetic algorithm was developed and expanded from the code presented by Michalewicz et al. [3] and was linked to AxiSuite using the Simulink add-on for Matlab. 
The default kinetic values stored within the AxiSuite model were used to generate a series of light-off curves under rich conditions for a number of gas species, including CO, NO, C_{3H8 and C3H6. These light-off curves were used to generate an objective function.} 
This objective function was used to generate a measure of fit for the kinetic parameters. The multi-objective genetic algorithm was subsequently used to search between specified limits to attempt to match the objective function. In total the pre-exponential factors and activation energies of ten reactions were simultaneously optimized. 
The results reported here demonstrate that, given accurate experimental data, the optimization algorithm is successful and robust in defining the correct kinetic parameters of a global kinetic model describing aftertreatment processes.

Algorithms for Hidden Markov Models With Imprecisely Specified Parameters

Hidden Markov models (HMMs) are widely used models for sequential data. As with other probabilistic graphical models, they require the specification of precise probability values, which can be too restrictive for some domains, especially when data are scarce or costly to acquire. We present a generalized version of HMMs, whose quantification can be done by sets of, instead of single, probability distributions. Our models have the ability to suspend judgment when there is not enough statistical evidence, and can serve as a sensitivity analysis tool for standard non-stationary HMMs. Efficient inference algorithms are developed to address standard HMM usage such as the computation of likelihoods and most probable explanations. Experiments with real data show that the use of imprecise probabilities leads to more reliable inferences without compromising efficiency.

Absolute magnitudes and slope parameters of Pan-STARRS PS1 asteroids -preliminary results

We present the study of absolute magnitude (H) and slope parameter (G) of 170,000 asteroids observed by the Pan-STARRS1 telescope during the period of 15 months within its 3-year all-sky survey mission. The exquisite photometry with photometric errors below 0.04 mag and well-defined filter and photometric system allowed to derive H and G with statistical and systematic errors. Our new approach lies in the Monte Carlo technique simulating rotation periods, amplitudes, and colors, and deriving most-likely H, G and their systematic errors. Comparison of H_M by Muinonen's phase function (Muinonen et al., 2010) with the Minor Planet Center database revealed a negative offset of 0.22±0.29 meaning that Pan-STARRS1 asteroids are fainter. We showed that the absolute magnitude derived by Muinonen's function is systematically larger on average by 0.14±0.29 and by 0.30±0.16 when assuming fixed slope parameter (G=0.15, G_{12}=0.53) than Bowell's absolute magnitude (Bowell et al., 1989). We also derived slope parameters of asteroids of known spectral types and showed a good agreement with the previous studies within the derived uncertainties. However, our systematic errors on G and G_{12} are significantly larger than in previous work, which is caused by poor temporal and phase coverage of vast majority of the detected asteroids. This disadvantage will vanish when full survey data will be available and ongoing extended and enhanced mission will provide new data.

Pooled analysis of tiotropium Respimat® pharmacokinetics in cystic fibrosis

Tiotropium is the first bronchodilator to be studied systematically in cystic fibrosis (CF). We investigated whether any intrinsic or extrinsic factors affected pharmacokinetic (PK) parameters of inhaled tiotropium delivered by Respimat® in adults and children with CF. Tiotropium PK in patients with CF was compared with that of healthy volunteers and patients with chronic obstructive pulmonary disease (COPD). This pooled analysis summarizes the PK parameters of inhaled tiotropium Respimat® across 9 early- and late-phase trials involving 27 healthy volunteers (1 trial), 409 patients with CF (3 trials), and 281 patients with COPD (5 trials). Patients with CF aged 5 to 11, 12 to 17, and ≥18 years had similar tiotropium plasma concentrations (geometric mean C_{0.083,ss,norm}: 2.22pg/mL/μg; not determined for patients aged <5 years). The fraction excreted unchanged in the urine was 3.4-fold lower for patients aged 0.4 to <5 years than for those aged 5 to 11 years (fe_0-4,ss: 1.19% vs 4.09%). Tiotropium PK parameters were similar between CF patients and COPD patients.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Absolute magnitudes and slope parameters for 250,000 asteroids observed by Pan-STARRS PS1 - preliminary results

We present the results of a Monte Carlo technique to calculate the absolute magnitudes (H) and slope parameters (G) of about 240,000 asteroids observed by the Pan-STARRS1 telescope during the first 15 months of its 3-year all-sky survey mission. The system's exquisite photometry with photometric errors asteroids rotation period, amplitude and color to derive the most-likely H and G, but its major advantage is in estimating realistic statistical+systematic uncertainties and errors on each parameter. The method was confirmed by comparison with the well-established and accurate results for about 500 asteroids provided by Pravec et al. (2012) and then applied to determining H and G for the Pan-STARRS1 asteroids using both the Muinonen et al. (2010) and Bowell et al. (1989) phase functions. Our results confirm the bias in MPC photometry discovered by ( Jurić et al., 2002).

An improved TLBO with elite strategy for parameters identification of PEM fuel cell and solar cell models

Clean and renewable energy generation and supply has drawn much attention worldwide in recent years, the proton exchange membrane (PEM) fuel cells and solar cells are among the most popular technologies. Accurately modeling the PEM fuel cells as well as solar cells is critical in their applications, and this involves the identification and optimization of model parameters. This is however challenging due to the highly nonlinear and complex nature of the models. In particular for PEM fuel cells, the model has to be optimized under different operation conditions, thus making the solution space extremely complex. In this paper, an improved and simplified teaching-learning based optimization algorithm (STLBO) is proposed to identify and optimize parameters for these two types of cell models. This is achieved by introducing an elite strategy to improve the quality of population and a local search is employed to further enhance the performance of the global best solution. To improve the diversity of the local search a chaotic map is also introduced. Compared with the basic TLBO, the structure of the proposed algorithm is much simplified and the searching ability is significantly enhanced. The performance of the proposed STLBO is firstly tested and verified on two low dimension decomposable problems and twelve large scale benchmark functions, then on the parameter identification of PEM fuel cell as well as solar cell models. Intensive experimental simulations show that the proposed STLBO exhibits excellent performance in terms of the accuracy and speed, in comparison with those reported in the literature.

Optimising drug solubilisation in amorphous polymer dispersions: rational selection of hot-melt extrusion processing parameters

The aim of this article was to construct a T–ϕ phase diagram for a model drug (FD) and amorphous polymer (Eudragit® EPO) and to use this information to understand the impact of how temperature–composition coordinates influenced the final properties of the extrudate. Defining process boundaries and understanding drug solubility in polymeric carriers is of utmost importance and will help in the successful manufacture of new delivery platforms for BCS class II drugs. Physically mixed felodipine (FD)–Eudragit® EPO (EPO) binary mixtures with pre-determined weight fractions were analysed using DSC to measure the endset of melting and glass transition temperature. Extrudates of 10 wt% FD–EPO were processed using temperatures (110°C, 126°C, 140°C and 150°C) selected from the temperature–composition (T–ϕ) phase diagrams and processing screw speed of 20, 100 and 200rpm. Extrudates were characterised using powder X-ray diffraction (PXRD), optical, polarised light and Raman microscopy. To ensure formation of a binary amorphous drug dispersion (ADD) at a specific composition, HME processing temperatures should at least be equal to, or exceed, the corresponding temperature value on the liquid–solid curve in a F–H T–ϕ phase diagram. If extruded between the spinodal and liquid–solid curve, the lack of thermodynamic forces to attain complete drug amorphisation may be compensated for through the use of an increased screw speed. Constructing F–H T–ϕ phase diagrams are valuable not only in the understanding drug–polymer miscibility behaviour but also in rationalising the selection of important processing parameters for HME to ensure miscibility of drug and polymer.