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.

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.

Microstructural evolution and mechanical properties of short-term thermally exposed 9/12Cr heat-resistant steels

The microstructural evolution during short-term (up to 3000 hours) thermal exposure of three 9/12Cr heat-resistant steels was studied, as well as the mechanical properties after exposure. The tempered martensitic lath structure, as well as the precipitation of carbide and MX type carbonitrides in the steel matrix, was stable after 3000 hours of exposure at 873 K (600 °C). A microstructure observation showed that during the short-term thermal exposure process, the change of mechanical properties was caused mainly by the formation and growth of Laves-phase precipitates in the steels. On thermal exposure, with an increase of cobalt and tungsten contents, cobalt could promote the segregation of tungsten along the martensite lath to form Laves phase, and a large size and high density of Laves-phase precipitates along the grain boundaries could lead to the brittle intergranular fracture of the steels.

Modeling of titanium alloys by using artificial neural networks

Titanium alloy exhibits an excellent combination of bio-compatibility, corrosion resistance, strength and toughness. The microstructure of an alloy influences the properties. The microstructures depend mainly on alloying elements, method of production, mechanical, and thermal treatments. The relationships between these variables and final properties of the alloy are complex, non-linear in nature, which is the biggest hurdle in developing proper correlations between them by conventional methods. So, we developed artificial neural networks (ANN) models for solving these complex phenomena in titanium alloys.

In the present work, ANN models were used for the analysis and prediction of the correlation between the process parameters, the alloying elements, microstructural features, beta transus temperature and mechanical properties in titanium alloys. Sensitivity analysis of trained neural network models were studied which resulted a better understanding of relationships between inputs and outputs. The model predictions and the analysis are well in agreement with the experimental results. The simulation results show that the average output-prediction error by models are less than 5% of the prediction range in more than 95% of the cases, which is quite acceptable for all metallurgical purposes.

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.