Modelling a biochemical reaction with computational fluid dynamics

Earlier investigations (Cartland Glover et al., 2004) into the use of computational fluid dynamics (CFD) for the modelling of gas-liquid and gas-liquid-solid flow allowed a simple biochemical reaction model to be implemented. A single plane mesh was used to represent the transport and reaction of molasses, the mould Aspergillus niger and citric acid in a bubble column with a height to diameter aspect ratio of 20:1. Two specific growth rates were used to examine the impact that biomass growth had on the local solids concentration and the effect this had on the local hydrodynamics of the bubble column.

The neural correlates of reading impairment in adults with developmental dyslexia:evidence from fMRI between group analyses and an fMRI multiple-case study

The goal of this project was to investigate the neural correlates of reading impairment in dyslexia as hypothesised by the main theories – the phonological deficit, visual magnocellular deficit and cerebellar deficit theories, with emphasis on individual differences. This research took a novel approach by: 1) contrasting the predictions in one sample of participants with dyslexia (DPs); 2) using a multiple-case study (and between-group comparisons) to investigate differences in BOLD between each DP and the controls (CPs); 3) demonstrating a possible relationship between reading impairment and its hypothesised neural correlates by using fMRI and a reading task. The multiple-case study revealed that the neural correlates of reading in dyslexia in all cases are not in agreement with the predictions of a single theory. The results show striking individual differences - even, where the neural correlates of reading in two DPs are consistent with the same theory, the areas can differ. A DP can exhibit under-engagement in an area in word, but not in pseudoword reading and vice versa, demonstrating that underactivation in that area cannot be interpreted as a ‘developmental lesion’. Additional analyses revealed complex results. Within-group analyses between behavioural measures and BOLD showed correlations in the predicted regions, areas outside ROI, and lack of correlations in some predicted areas. Comparisons of subgroups which differed on Orthography Composite supported the MDT, but only for Words. The results suggest that phonological scores are not a sufficient predictor of the under-engagement of phonological areas during reading. DPs and CPs exhibited correlations between Purdue Pegboard Composite and BOLD in cerebellar areas only for Pseudowords. Future research into reading in dyslexia should use a more holistic approach, involving genetic and environmental factors, gene by environment interaction, and comorbidity with other disorders. It is argued that multidisciplinary research, within the multiple-deficit model holds significant promise here.

The MIMIC model and formative variables:problems and solutions

The use of the multiple indicators, multiple causes model to operationalize formative variables (the formative MIMIC model) is advocated in the methodological literature. Yet, contrary to popular belief, the formative MIMIC model does not provide a valid method of integrating formative variables into empirical studies and we recommend discarding it from formative models. Our arguments rest on the following observations. First, much formative variable literature appears to conceptualize a causal structure between the formative variable and its indicators which can be tested or estimated. We demonstrate that this assumption is illogical, that a formative variable is simply a researcher-defined composite of sub-dimensions, and that such tests and estimates are unnecessary. Second, despite this, researchers often use the formative MIMIC model as a means to include formative variables in their models and to estimate the magnitude of linkages between formative variables and their indicators. However, the formative MIMIC model cannot provide this information since it is simply a model in which a common factor is predicted by some exogenous variables—the model does not integrate within it a formative variable. Empirical results from such studies need reassessing, since their interpretation may lead to inaccurate theoretical insights and the development of untested recommendations to managers. Finally, the use of the formative MIMIC model can foster fuzzy conceptualizations of variables, particularly since it can erroneously encourage the view that a single focal variable is measured with formative and reflective indicators. We explain these interlinked arguments in more detail and provide a set of recommendations for researchers to consider when dealing with formative variables.

Detection and quantification of protein oxidation in sarcopenic models:a mass spectrometry study

Oxidised biomolecules in aged tissue could potentially be used as biomarkers for age-related diseases; however, it is still unclear whether they causatively contribute to ageing or are consequences of the ageing process. To assess the potential of using protein oxidation as markers of ageing, mass spectrometry (MS) was employed for the identification and quantification of oxidative modifications in obese (ob/ob) mice. Lean muscle mass and strength is reduced in obesity, representing a sarcopenic model in which the levels of oxidation can be evaluated for different muscular systems including calcium homeostasis, metabolism and contractility. Several oxidised residues were identified by tandem MS (MS/MS) in both muscle homogenate and isolated sarcoplasmic reticulum (SR), an organelle that regulates intracellular calcium levels in muscle. These modifications include oxidation of methionine, cysteine, tyrosine, and tryptophan in several proteins such as sarcoplasmic reticulum calcium ATPase (SERCA), glycogen phosphorylase, and myosin. Once modifications had been identified, multiple reaction monitoring MS (MRM) was used to quantify the percentage modification of oxidised residues within the samples. Preliminary data suggests proteins in ob/ob mice are more oxidised than the controls. For example SERCA, which constitutes 60-70% of the SR, had approximately a 2-fold increase in cysteine trioxidation of Cys561 in the obese model when compared to the control. Other obese muscle proteins have also shown a similar increase in oxidation for various residues. Further analysis with complex protein mixtures will determine the potential diagnostic use of MRM experiments for analysing protein oxidation in small biological samples such as muscle needle biopsies.

Protein oxidative modifications:new mass spectrometry approaches to detection in biological samples

In inflammatory diseases, release of oxidants leads to oxidative damage to proteins. The precise nature of oxidative damage to individual proteins depends on the oxidant involved. Chlorination and nitration are markers of modification by the myeloperoxidase-H2O2-Cl- system and nitric oxide-derived oxidants, respectively. Although these modifications can be detected by western blotting, currently no reliable method exists to identify the specific sites damage to individual proteins in complex mixtures such as clinical samples. We are developing novel LCMS2 and precursor ion scanning methods to address this. LC-MS2 allows separation of peptides and detection of mass changes in oxidized residues on fragmentation of the peptides. We have identified indicative fragment ions for chlorotyrosine, nitrotyrosine, hydroxytyrosine and hydroxytryptophan. A nano-LC/MS3 method involving the dissociation of immonium ions to give specific fragments for the oxidized residues has been developed to overcome the problem of false positives from ions isobaric to these immonium ions that exist in unmodified peptides. The approach has proved able to identify precise protein modifications in individual proteins and mixtures of proteins. An alternative methodology involves multiple reaction monitoring for precursors and fragment ions are specific to oxidized and chlorinated proteins, and this has been tested with human serum albumin. Our ultimate aim is to apply this methodology to the detection of oxidative post-translational modifications in clinical samples for disease diagnosis, monitoring the outcomes of therapy, and improved understanding of disease biochemistry.

The normal optic nerve variations in an optometric population:Possible ocular and systemic influences

Early detection of glaucoma relies on a detailed knowledge of how the normal optic nerve (ONH) varies within the population. The purpose of this study focused on two main areas; 1. To explore the optic nerve head appearance in the normal optometric population and compare the south Asian (principally Pakistani) with the European white population, correcting for possible ocular and non-ocular influences in a multiple regression model. The main findings were: • The optic discs of the South Asian (SA) and White European (WE) populations were not statistically different in size. The SA group possessed discs with increased cupping and thinner neuro-retinal rims (NRR) compared with the WE group. The SA group also demonstrated a more vertically oval shape than the WE population. These differences were significant at the p<0.01 level. • The upper limits of inter-eye asymmetry were: ≤0.2 for cup to disc area ratio, and 3mmHg for intra-ocular pressure (IOP) for both ethnic groups and this did not increase with age. IOP asymmetry did not vary with gender, ethnicity or a family history of glaucoma and was independent of ONH asymmetry. ONH and IOP asymmetry are therefore independent risk factors when screening for glaucoma for both ethnic groups. 2. To investigate the validity of the ISNT rule: inferior> superior> nasal> temporal NRR thickness in the optometric population. The main findings were: • As disc size increased the disc become rounder and less vertically oval in shape. Vertically oval discs had thicker superior and inferior NRRs and thinner nasal and temporal NRRs compared with rounder disc shapes due to cup shape being independent of disc shape. Vertically oval discs were therefore more likely to obey the ISNT rule than larger rounder discs. • The ISNT rule has a low adherence in our sample of normal eyes (5.7%). However, by removing the nasal sector to become the IST rule, 74.5% of normal eyes obeyed. SA eyes and female gender were more likely to obey the ISNT rule due to increased disc ovality. The IST rule is independent of disc shape and therefore more suitable for assessing discs from both ethnic backgrounds. Obeying the ISNT rule or IST rule was not related to disc or cup size.

Visual field severity indices demonstrate dose-dependent visual loss from Vigabatrin therapy

Purpose: The aims of this study were to develop an algorithm to accurately quantify Vigabatrin (VGB)-induced central visual field loss and to investigate the relationship between visual field loss and maximum daily dose, cumulative dose and duration of dose. Methods: The sample comprised 31 patients (mean age 37.9 years; SD 14.4 years) diagnosed with epilepsy and exposed to VGB. Each participant underwent standard automated static visual field examination of the central visual field. Central visual field loss was determined using continuous scales quantifying severity in terms of area and depth of defect and additionally by symmetry of defect between the two eyes. A simultaneous multiple regression model was used to explore the relationship between these visual field parameters and the drug predictor variables. Results: The regression model indicated that maximum VGB dose was the only factor to be significantly correlated with individual eye severity (right eye: p = 0.020; left eye: p = 0.012) and symmetry of visual field defect (p = 0.024). Conclusions: Maximum daily dose was the single most reliable indicator of those patients likely to exhibit visual field defects due to VGB. These findings suggest that high maximum dose is more likely to result in visual field defects than high cumulative doses or those of long duration.

Kinetic study of the pyrolysis of miscanthus and its acid hydrolysis residue by thermogravimetric analysis

The kinetic parameters of the pyrolysis of miscanthus and its acid hydrolysis residue (AHR) were determined using thermogravimetric analysis (TGA). The AHR was produced at the University of Limerick by treating miscanthus with 5 wt.% sulphuric acid at 175 °C as representative of a lignocellulosic acid hydrolysis product. For the TGA experiments, 3 to 6 g of sample, milled and sieved to a particle size below 250 μm, were placed in the TGA ceramic crucible. The experiments were carried out under non-isothermal conditions heating the samples from 50 to 900 °C at heating rates of 2.5, 5, 10, 17 and 25 °C/min. The activation energy (EA) of the decomposition process was determined from the TGA data by differential analysis (Friedman) and three isoconversional methods of integral analysis (Kissinger–Akahira–Sunose, Ozawa–Flynn–Wall, Vyazovkin). The activation energy ranged from 129 to 156 kJ/mol for miscanthus and from 200 to 376 kJ/mol for AHR increasing with increasing conversion. The reaction model was selected using the non-linear least squares method and the pre-exponential factor was calculated from the Arrhenius approximation. The results showed that the best fitting reaction model was the third order reaction for both feedstocks. The pre-exponential factor was in the range of 5.6 × 1010 to 3.9 × 10+ 13 min− 1 for miscanthus and 2.1 × 1016 to 7.7 × 1025 min− 1 for AHR.

Thermal processing of miscanthus, sugarcane bagasse, sugarcane trash and their acid hydrolysis residues

The research presented in this thesis was developed as part of DIBANET, an EC funded project aiming to develop an energetically self-sustainable process for the production of diesel miscible biofuels (i.e. ethyl levulinate) via acid hydrolysis of selected biomass feedstocks. Three thermal conversion technologies, pyrolysis, gasification and combustion, were evaluated in the present work with the aim of recovering the energy stored in the acid hydrolysis solid residue (AHR). Mainly consisting of lignin and humins, the AHR can contain up to 80% of the energy in the original feedstock. Pyrolysis of AHR proved unsatisfactory, so attention focussed on gasification and combustion with the aim of producing heat and/or power to supply the energy demanded by the ethyl levulinate production process. A thermal processing rig consisting on a Laminar Entrained Flow Reactor (LEFR) equipped with solid and liquid collection and online gas analysis systems was designed and built to explore pyrolysis, gasification and air-blown combustion of AHR. Maximum liquid yield for pyrolysis of AHR was 30wt% with volatile conversion of 80%. Gas yield for AHR gasification was 78wt%, with 8wt% tar yields and conversion of volatiles close to 100%. 90wt% of the AHR was transformed into gas by combustion, with volatile conversions above 90%. 5volO2%-95vol%N2 gasification resulted in a nitrogen diluted, low heating value gas (2MJ/m3). Steam and oxygen-blown gasification of AHR were additionally investigated in a batch gasifier at KTH in Sweden. Steam promoted the formation of hydrogen (25vol%) and methane (14vol%) improving the gas heating value to 10MJ/m3, below the typical for steam gasification due to equipment limitations. Arrhenius kinetic parameters were calculated using data collected with the LEFR to provide reaction rate information for process design and optimisation. Activation energy (EA) and pre-exponential factor (ko in s-1) for pyrolysis (EA=80kJ/mol, lnko=14), gasification (EA=69kJ/mol, lnko=13) and combustion (EA=42kJ/mol, lnko=8) were calculated after linearly fitting the data using the random pore model. Kinetic parameters for pyrolysis and combustion were also determined by dynamic thermogravimetric analysis (TGA), including studies of the original biomass feedstocks for comparison. Results obtained by differential and integral isoconversional methods for activation energy determination were compared. Activation energy calculated by the Vyazovkin method was 103-204kJ/mol for pyrolysis of untreated feedstocks and 185-387kJ/mol for AHRs. Combustion activation energy was 138-163kJ/mol for biomass and 119-158 for AHRs. The non-linear least squares method was used to determine reaction model and pre-exponential factor. Pyrolysis and combustion of biomass were best modelled by a combination of third order reaction and 3 dimensional diffusion models, while AHR decomposed following the third order reaction for pyrolysis and the 3 dimensional diffusion for combustion.

Probabilistic multiple model neural network based leak detection system:experimental study

This paper presents an effective decision making system for leak detection based on multiple generalized linear models and clustering techniques. The training data for the proposed decision system is obtained by setting up an experimental pipeline fully operational distribution system. The system is also equipped with data logging for three variables; namely, inlet pressure, outlet pressure, and outlet flow. The experimental setup is designed such that multi-operational conditions of the distribution system, including multi pressure and multi flow can be obtained. We then statistically tested and showed that pressure and flow variables can be used as signature of leak under the designed multi-operational conditions. It is then shown that the detection of leakages based on the training and testing of the proposed multi model decision system with pre data clustering, under multi operational conditions produces better recognition rates in comparison to the training based on the single model approach. This decision system is then equipped with the estimation of confidence limits and a method is proposed for using these confidence limits for obtaining more robust leakage recognition results.

Parametric pattern selection in a reaction-diffusion model

We compare spot patterns generated by Turing mechanisms with those generated by replication cascades, in a model one-dimensional reaction-diffusion system. We determine the stability region of spot solutions in parameter space as a function of a natural control parameter (feed-rate) where degenerate patterns with different numbers of spots coexist for a fixed feed-rate. While it is possible to generate identical patterns via both mechanisms, we show that replication cascades lead to a wider choice of pattern profiles that can be selected through a tuning of the feed-rate, exploiting hysteresis and directionality effects of the different pattern pathways.

Bifurcation and stability in a low-dimensional model for multiple-frequency mode-locked lasers

Recent theoretical investigations have demonstrated that the stability of mode-locked solutions of multiple frequency channels depends on the degree of inhomogeneity in gain saturation. In this article, these results are generalized to determine conditions on each of the system parameters necessary for both the stability and the existence of mode-locked pulse solutions for an arbitrary number of frequency channels. In particular, we find that the parameters governing saturable intensity discrimination and gain inhomogeneity in the laser cavity also determine the position of bifurcations of solution types. These bifurcations are completely characterized in terms of these parameters. In addition to influencing the stability of mode-locked solutions, we determine a balance between cubic gain and quintic loss, which is necessary for the existence of solutions as well. Furthermore, we determine the critical degree of inhomogeneous gain broadening required to support pulses in multiple-frequency channels. © 2010 The American Physical Society.