Kinetic modelling of growth and storage molecule production in microalgae under mixotrophic and autotrophic conditions

In order to improve algal biofuel production on a commercial-scale, an understanding of algal growth and fuel molecule accumulation is essential. A mathematical model is presented that describes biomass growth and storage molecule (TAG lipid and starch) accumulation in the freshwater microalga Chlorella vulgaris, under mixotrophic and autotrophic conditions. Biomass growth was formulated based on the Droop model, while the storage molecule production was calculated based on the carbon balance within the algal cells incorporating carbon fixation via photosynthesis, organic carbon uptake and functional biomass growth. The model was validated with experimental growth data of C. vulgaris and was found to fit the data well. Sensitivity analysis showed that the model performance was highly sensitive to variations in parameters associated with nutrient factors, photosynthesis and light intensity. The maximum productivity and biomass concentration were achieved under mixotrophic nitrogen sufficient conditions, while the maximum storage content was obtained under mixotrophic nitrogen deficient conditions. © 2014 Elsevier Ltd.

Modelling of turbulent lifted jet flames using flamelets: A priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z-c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z-c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences. © 2014 Taylor & Francis.

Temporal and spatial variability in nitrogen uptake kinetics during harmful dinoflagellate blooms in the East China Sea

Temporal and spatial variability in the kinetic parameters of uptake of nitrate (NO3-), ammonium (NH4+), urea, and glycine was measured during dinoflagellate blooms in Changjiang River estuary and East China Sea coast, 2005. Karenia mikimotoi was the dominant species in the early stage of the blooms and was succeeded by Prorocentrum donghaiense. The uptake of nitrogen (N) was determined using N-15 tracer techniques. The results of comparison kinetic parameters with ambient nutrients confirmed that different N forms were preferentially taken up during different stages of the bloom. NO3- (V-max 0.044 h(-1); K-s 60.8 mu M-N) was an important N source before it was depleted. NH4+ (V-max 0.049 h(-1); K-s 2.15 mu M-N) was generally the preferred N. Between the 2 organic N sources, urea was more preferred when K. mikimotoi dominated the bloom (V-max 0.020 h(-1); K-s 1.35 mu M-N) and glycine, considered as a dominant amino acid, was more preferred when P. donghaiense dominated the bloom (V-max 0.025 h(-1); K-s 1.76 mu M-N). The change of N uptake preference by the bloom-forming algae was also related to the variation in ambient N concentrations. Published by Elsevier B.V.

Multi-scale modelling of atomic layer deposition

High-permittivity ("high-k") dielectric materials are used in the transistor gate stack in integrated circuits. As the thickness of silicon oxide dielectric reduces below 2 nm with continued downscaling, the leakage current because of tunnelling increases, leading to high power consumption and reduced device reliability. Hence, research concentrates on finding materials with high dielectric constant that can be easily integrated into a manufacturing process and show the desired properties as a thin film. Atomic layer deposition (ALD) is used practically to deposit high-k materials like HfO2, ZrO2, and Al2O3 as gate oxides. ALD is a technique for producing conformal layers of material with nanometer-scale thickness, used commercially in non-planar electronics and increasingly in other areas of science and technology. ALD is a type of chemical vapor deposition that depends on self-limiting surface chemistry. In ALD, gaseous precursors are allowed individually into the reactor chamber in alternating pulses. Between each pulse, inert gas is admitted to prevent gas phase reactions. This thesis provides a profound understanding of the ALD of oxides such as HfO2, showing how the chemistry affects the properties of the deposited film. Using multi-scale modelling of ALD, the kinetics of reactions at the growing surface is connected to experimental data. In this thesis, we use density functional theory (DFT) method to simulate more realistic models for the growth of HfO2 from Hf(N(CH3)2)4/H2O and HfCl4/H2O and for Al2O3 from Al(CH3)3/H2O.Three major breakthroughs are discovered. First, a new reaction pathway, ’multiple proton diffusion’, is proposed for the growth of HfO2 from Hf(N(CH3)2)4/H2O.1 As a second major breakthrough, a ’cooperative’ action between adsorbed precursors is shown to play an important role in ALD. By this we mean that previously-inert fragments can become reactive once sufficient molecules adsorb in their neighbourhood during either precursor pulse. As a third breakthrough, the ALD of HfO2 from Hf(N(CH3)2)4 and H2O is implemented for the first time into 3D on-lattice kinetic Monte-Carlo (KMC).2 In this integrated approach (DFT+KMC), retaining the accuracy of the atomistic model in the higher-scale model leads to remarkable breakthroughs in our understanding. The resulting atomistic model allows direct comparison with experimental techniques such as X-ray photoelectron spectroscopy and quartz crystal microbalance.

Economic modelling of antenatal screening and ultrasound scanning programmes for identification of fetal abnormalities

Objective Within the framework of a health technology assessment and using an economic model, to determine the most clinically and cost effective policy of scanning and screening for fetal abnormalities in early pregnancy. Design A discrete event simulation model of 50,000 singleton pregnancies. Setting Maternity services in Scotland. Population Women during the first 24 weeks of their pregnancy. Methods The mathematical model was populated with data on uptake of screening, prevalence, detection and false positive rates for eight fetal abnormalities and with costs for ultrasound scanning and serum screening. Inclusion of abnormalities was based on the relative prevalence and clinical importance of conditions and the availability of data. Six strategies for the identification of abnormalities prenatally including combinations of first and second trimester ultrasound scanning and first and second trimester screening for chromosomal abnormalities were compared. Main outcome measures The number of abnormalities detected and missed, the number of iatrogenic losses resulting from invasive tests, the total cost of strategies and the cost per abnormality detected were compared between strategies. Results First trimester screening for chromosomal abnormalities costs more than second trimester screening but results in fewer iatrogenic losses. Strategies which include a second trimester ultrasound scan result in more abnormalities being detected and have lower costs per anomaly detected. Conclusions The preferred strategy includes both first and second trimester ultrasound scans and a first trimester screening test for chromosomal abnormalities. It has been recommended that this policy is offered to all women in Scotland.

Studies of switching kinetics in ferroelectric thin films

Experimental studies are reported concerning the importance of interfacial capacitance (including electrode screening, space-charge layers, and/or chemically discrete dead layers). on domain switching behaviour in thin films of ferroelectric lead zirconate-titanate (PZT), strontium bismuth tantalate (SBT), and barium strontium titanate (BST). Emphasis is placed upon studies at applied field values very near the coercive field E, asymmetry in fatigue for positive and negative polarity coercive fields, and in the case of BST, of the coexistence of ferroelectric and paraelectric phases Studies of dielectric loss show important correlations between tan 6 and fatigue (polarization decrease) as a function of bipolar switching cycles N. This is a priori not obvious, since the former is a linear response and the latter, a nonlinear response. Modelling of enlarged interfacial,space-charge layers in PZT films and chemically distinct dead (paraelectric) layers in BST films shows contradictory tendencies of coercive-voltage changes with the growth of passive layers.

Reproductive biomass in Holcus lanatus clones that differ in their phosphate uptake kinetics and mycorrhizal colonization

In normal populations of the common grass Holcus lanatus there is a polymorphism for arsenate resistance, manifested as suppressed phosphate uptake (SPU), and controlled by a major gene with dominant expression. A natural population of SPU plants had greater arbuscular-mycorrhizal colonization than wild type, nonSPU plants. It was hypothesized that, in order to survive alongside plants with a normal rate of phosphate (P) uptake, SPU plants would be more dependent on mycorrhizal associations. We performed an experiment using plants with SPU phenotypes from both arsenate mine spoils and uncontaminated soils, as well as plants with a nonSPU phenotype. They were grown with and without a mycorrhizal inoculum and added N, which altered plant P requirements. We showed that grasses with SPU phenotypes accumulated more shoot P than nonSPU plants, the opposite of the expected result. SPY plants also produced considerably more flower panicles, and had greater shoot and root biomass. The persistence of SPU phenotypes in normal populations is not necessarily related to mycorrhizal colonization as there were no differences in percentage AM colonization between the phenotypes. Being mycorrhizal reduced flower biomass production, as mycorrhizal SPU plants had lower shoot P concentrations and produced fewer flower panicles than non-mycorrhizal, nonSPU plants. We now hypothesize that the SPU phenotype is brought about by a genotype that results in increased accumulation of P in shoots, and that suppression of the rate of uptake is a consequence of this high shoot P concentration, operating by means of a homeostatic feedback mechanism. We also postulate that increased flower production is linked to a high shoot P concentration. SPU plants thus allocate more resources into seed production, leading to a higher frequency of SPU genes. Increased reproductive allocation reduces vegetative allocation and may affect competitive ability and hence survival, explaining the maintenance of the polymorphism. As mycorrhizal SPU plants behave more like nonSPU plants, AM colonization itself could play a major part in the maintenance of the SPU polymorphism.

Mathematical modelling of quantum yield enhancements of methyl orange photooxidation in aqueous TiO2 suspensions under controlled periodic UV LED illumination

Quantum yields of the photocatalytic degradation of methyl orange under controlled periodic illumination (CPI) have been modelled using existing models. A modified Langmuir-Hinshelwood (L-H) rate equation was used to predict the degradation reaction rates of methyl orange at various duty cycles and a simple photocatalytic model was applied in modelling quantum yield enhancement of the photocatalytic process due to the CPI effect. A good agreement between the modelled and experimental data was observed for quantum yield modelling. The modified L-H model, however, did not accurately predict the photocatalytic decomposition of the dye under periodic illumination.

Constraint-based modelling of mixed microbial populations: Application to polyhydroxyalkanoates production

Dissertação para obtenção do Grau de Doutor em Engenharia Química e Bioquímica

Mathematical Modelling of the Dynamics of Granular Materials in a Rotating Cylinder

The current study is aimed at the development of a theoretical simulation tool based on Discrete Element Method (DEM) to 'interpret granular dynamics of solid bed in the cross section of the horizontal rotating cylinder at the microscopic level and subsequently apply this model to establish the transition behaviour, mixing and segregation.The simulation of the granular motion developed in this work is based on solving Newton's equation of motion for each particle in the granular bed subjected to the collisional forces, external forces and boundary forces. At every instant of time, the forces are tracked and the positions velocities and accelarations of each partcle is The software code for this simulation is written in VISUAL FORTRAN 90 After checking the validity of the code with special tests, it is used to investigate the transition behaviour of granular solids motion in the cross section of a rotating cylinder for various rotational speeds and fill fraction.This work is hence directed towards a theoretical investigation based on Discrete Element Method (DEM) of the motion of granular solids in the radial direction of the horizontal cylinder to elucidate the relationship between the operating parameters of the rotating cylinder geometry and physical properties ofthe granular solid.The operating parameters of the rotating cylinder include the various rotational velocities of the cylinder and volumetric fill. The physical properties of the granular solids include particle sizes, densities, stiffness coefficients, and coefficient of friction Further the work highlights the fundamental basis for the important phenomena of the system namely; (i) the different modes of solids motion observed in a transverse crosssection of the rotating cylinder for various rotational speeds, (ii) the radial mixing of the granular solid in terms of active layer depth (iii) rate coefficient of mixing as well as the transition behaviour in terms of the bed turnover time and rotational speed and (iv) the segregation mechanisms resulting from differences in the size and density of particles.The transition behaviour involving its six different modes of motion of the granular solid bed is quantified in terms of Froude number and the results obtained are validated with experimental and theoretical results reported in the literature The transition from slumping to rolling mode is quantified using the bed turnover time and a linear relationship is established between the bed turn over time and the inverse of the rotational speed of the cylinder as predicted by Davidson et al. [2000]. The effect of the rotational speed, fill fraction and coefficient of friction on the dynamic angle of repose are presented and discussed. The variation of active layer depth with respect to fill fraction and rotational speed have been investigated. The results obtained through simulation are compared with the experimental results reported by Van Puyvelde et. at. [2000] and Ding et at. [2002].The theoretical model has been further extended, to study the rmxmg and segregation in the transverse direction for different particle sizes and their size ratios. The effect of fill fraction and rotational speed on the transverse mixing behaviour is presented in the form of a mixing index and mixing kinetics curve. The segregation pattern obtained by the simulation of the granular solid bed with respect to the rotational speed of the cylinder is presented both in graphical and numerical forms. The segregation behaviour of the granular solid bed with respect to particle size, density and volume fraction of particle size has been investigated. Several important macro parameters characterising segregation such as mixing index, percolation index and segregation index have been derived from the simulation tool based on first principles developed in this work.

A comparison of two time-domain analysis procedures in the determination of VO2 kinetics by pseudorandom binary sequence exercise testing

The purpose of this study was to apply and compare two time-domain analysis procedures in the determination of oxygen uptake (VO2) kinetics in response to a pseudorandom binary sequence (PRBS) exercise test. PRBS exercise tests have typically been analysed in the frequency domain. However, the complex interpretation of frequency responses may have limited the application of this procedure in both sporting and clinical contexts, where a single time measurement would facilitate subject comparison. The relative potential of both a mean response time (MRT) and a peak cross-correlation time (PCCT) was investigated. This study was divided into two parts: a test-retest reliability study (part A), in which 10 healthy male subjects completed two identical PRBS exercise tests, and a comparison of the VO2 kinetics of 12 elite endurance runners (ER) and 12 elite sprinters (SR; part B). In part A, 95% limits of agreement were calculated for comparison between MRT and PCCT. The results of part A showed no significant difference between test and retest as assessed by MRT [mean (SD) 42.2 (4.2) s and 43.8 (6.9) s] or by PCCT [21.8 (3.7) s and 22.7 (4.5) s]. Measurement error (%) was lower for MRT in comparison with PCCT (16% and 25%, respectively). In part B of the study, the VO2 kinetics of ER were significantly faster than those of SR, as assessed by MRT [33.4 (3.4) s and 39.9 (7.1) s, respectively; P<0.01] and PCCT [20.9 (3.8) s and 24.8 (4.5) s; P < 0.05]. It is possible that either analysis procedure could provide a single test measurement Of VO2 kinetics; however, the greater reliability of the MRT data suggests that this method has more potential for development in the assessment Of VO2 kinetics by PRBS exercise testing.

Respirometric techniques for assessment of biological kinetics in constructed wetland

The aim of this research was the development of a procedure to measure biological kinetics of organic matter oxidation and nitrification in constructed wetland, by using respirometric techniques. Columns simulating cores of vertical subsurface flow systems were investigated. The oxygen uptake rate (OUR) of the columns was calculated on the basis of the difference of DO concentrations measured continuously at the top and at the bottom of the column. From the respirogram, the following kinetic parameters have been evaluated: maximum rate of oxidation of readily biodegradable COD, maximum rate of nitrification, endogenous respiration of the biomass grown inside the bed. In order to improve the interpretation of the respirograms, additional respirometric tests were carried out on the wetland columns by using pure substrates, such as acetate (carbon source) and ammonium (substrate for nitrification). The kinetic parameters obtained from respirograms can be useful for control and design of constructed wetlands or for improving nutrient and carbon mass balances.

Muscle damage slows oxygen uptake kinetics during moderate-intensity exercise performed at high pedal rate

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Oxygen uptake kinetics and time to exhaustion in cycling and running: A comparison between trained and untrained subjects

The objective of the present study was to compare pulmonary gas exchange kinetics (VO 2 kinetics) and time to exhaustion (Tlim) between trained and untrained individuals during severe exercise performed on a cycle ergometer and treadmill. Eleven untrained males in running (UR) and cycling (UC), nine endurance cyclists (EC), and seven endurance runners (ER) were submitted to the following tests on separate days: (i) incremental test for determination of maximal oxygen uptake (VO 2max) and the intensity associated with the achievement of VO 2max (IVO 2max) on a mechanical braked cycle ergometer (EC and UC) and on a treadmill (ER and UR); (ii) all-out exercise bout performed at IVO 2max to determine the time to exhaustion at IVO 2max (Tlim) and the time constant of oxygen uptake kinetics (τ). The τ was significantly faster in trained group, both in cycling (EC = 28.2 ± 4.7 s; UC = 63.8 ± 25.0 s) and in running (ER = 28.5 ± 8.5 s; UR = 59.3 ± 12.0 s). Tlim of untrained was significantly lower in cycling (EC = 384.4 ± 66.6 s vs. UC; 311.1 ± 105.7 s) and higher in running (ER = 309.2 ± 176.6 s vs. UR = 439.8 ± 104.2 s). We conclude that the VO 2 kinetic response at the onset of severe exercise, carried out at the same relative intensity is sensitive to endurance training, irrespective of the exercise type. The endurance training seems to differently influence Tlim during exercise at IVO 2max in running and cycling. © 2003 Taylor & Francis Ltd.

Pharmacokinetic and pharmacodynamic modelling of intravenous, intramuscular and subcutaneous buprenorphine in conscious cats

Objective To describe simultaneous pharmacokinetics (PK) and thermal antinociception after intravenous (IV), intramuscular (IM) and subcutaneous (SC) buprenorphine in cats. Study design Randomized, prospective, blinded, three period crossover experiment. Animals Six healthy adult cats weighing 4.1±0.5kg. Methods Buprenorphine (0.02mgkg-1) was administered IV, IM or SC. Thermal threshold (TT) testing and blood collection were conducted simultaneously at baseline and at predetermined time points up to 24hours after administration. Buprenorphine plasma concentrations were determined by liquid chromatography tandem mass spectrometry. TT was analyzed using anova (p<0.05). A pharmacokinetic-pharmacodynamic (PK-PD) model of the IV data was described using a model combining biophase equilibration and receptor association-dissociation kinetics. Results TT increased above baseline from 15 to 480minutes and at 30 and 60minutes after IV and IM administration, respectively (p<0.05). Maximum increase in TT (mean±SD) was 9.3±4.9°C at 60minutes (IV), 4.6±2.8°C at 45minutes (IM) and 1.9±1.9°C at 60minutes (SC). TT was significantly higher at 15, 60, 120 and 180minutes, and at 15, 30, 45, 60 and 120minutes after IV administration compared to IM and SC, respectively. IV and IM buprenorphine concentration-time data decreased curvilinearly. SC PK could not be modeled due to erratic absorption and disposition. IV buprenorphine disposition was similar to published data. The PK-PD model showed an onset delay mainly attributable to slow biophase equilibration (t1/2ke0=47.4minutes) and receptor binding (kon=0.011mL ng-1minute-1). Persistence of thermal antinociception was due to slow receptor dissociation (t1/2koff=18.2minutes). Conclusions and clinical relevance IV and IM data followed classical disposition and elimination in most cats. Plasma concentrations after IV administration were associated with antinociceptive effect in a PK-PD model including negative hysteresis. At the doses administered, the IV route should be preferred over the IM and SC routes when buprenorphine is administered to cats. © 2012 Association of Veterinary Anaesthetists and the American College of Veterinary Anesthesiologists.