Knowledge management in intelligent tutoring systems

In the last years, Intelligent Tutoring Systems have been a very successful way for improving learning experience. Many issues must be addressed until this technology can be defined mature. One of the main problems within the Intelligent Tutoring Systems is the process of contents authoring: knowledge acquisition and manipulation processes are difficult tasks because they require a specialised skills on computer programming and knowledge engineering. In this thesis we discuss a general framework for knowledge management in an Intelligent Tutoring System and propose a mechanism based on first order data mining to partially automate the process of knowledge acquisition that have to be used in the ITS during the tutoring process. Such a mechanism can be applied in Constraint Based Tutor and in the Pseudo-Cognitive Tutor. We design and implement a part of the proposed architecture, mainly the module of knowledge acquisition from examples based on first order data mining. We then show that the algorithm can be applied at least two different domains: first order algebra equation and some topics of C programming language. Finally we discuss the limitation of current approach and the possible improvements of the whole framework.

MCAO for Extremely Large Telescopes: the cases of LBT and E-ELT

Several MCAO systems are under study to improve the angular resolution of the current and of the future generation large ground-based telescopes (diameters in the 8-40 m range). The subject of this PhD Thesis is embedded in this context. Two MCAO systems, in dierent realization phases, are addressed in this Thesis: NIRVANA, the 'double' MCAO system designed for one of the interferometric instruments of LBT, is in the integration and testing phase; MAORY, the future E-ELT MCAO module, is under preliminary study. These two systems takle the sky coverage problem in two dierent ways. The layer oriented approach of NIRVANA, coupled with multi-pyramids wavefront sensors, takes advantage of the optical co-addition of the signal coming from up to 12 NGS in a annular 2' to 6' technical FoV and up to 8 in the central 2' FoV. Summing the light coming from many natural sources permits to increase the limiting magnitude of the single NGS and to improve considerably the sky coverage. One of the two Wavefront Sensors for the mid- high altitude atmosphere analysis has been integrated and tested as a stand- alone unit in the laboratory at INAF-Osservatorio Astronomico di Bologna and afterwards delivered to the MPIA laboratories in Heidelberg, where was integrated and aligned to the post-focal optical relay of one LINC-NIRVANA arm. A number of tests were performed in order to characterize and optimize the system functionalities and performance. A report about this work is presented in Chapter 2. In the MAORY case, to ensure correction uniformity and sky coverage, the LGS-based approach is the current baseline. However, since the Sodium layer is approximately 10 km thick, the articial reference source looks elongated, especially when observed from the edge of a large aperture. On a 30-40 m class telescope, for instance, the maximum elongation varies between few arcsec and 10 arcsec, depending on the actual telescope diameter, on the Sodium layer properties and on the laser launcher position. The centroiding error in a Shack-Hartmann WFS increases proportionally to the elongation (in a photon noise dominated regime), strongly limiting the performance. To compensate for this effect a straightforward solution is to increase the laser power, i.e. to increase the number of detected photons per subaperture. The scope of Chapter 3 is twofold: an analysis of the performance of three dierent algorithms (Weighted Center of Gravity, Correlation and Quad-cell) for the instantaneous LGS image position measurement in presence of elongated spots and the determination of the required number of photons to achieve a certain average wavefront error over the telescope aperture. An alternative optical solution to the spot elongation problem is proposed in Section 3.4. Starting from the considerations presented in Chapter 3, a first order analysis of the LGS WFS for MAORY (number of subapertures, number of detected photons per subaperture, RON, focal plane sampling, subaperture FoV) is the subject of Chapter 4. An LGS WFS laboratory prototype was designed to reproduce the relevant aspects of an LGS SH WFS for the E-ELT and to evaluate the performance of different centroid algorithms in presence of elongated spots as investigated numerically and analytically in Chapter 3. This prototype permits to simulate realistic Sodium proles. A full testing plan for the prototype is set in Chapter 4.

Darstellung von substituierten Phenylacetylenen und Untersuchung ihrer von Palladiumverbindungen und Phosphorliganden homogen katalysierten Dimerisierung

In dieser Arbeit werden Homodimerisierungen von substituierten Phenylacetylenen mit Butinsäureethylester untersucht. Als Katalysatorsysteme dienen Palladium-Phosphan-Systeme.Kinetische Untersuchungen dieser Reaktionen zeigen, dass im Allgemeinen keine einfache Reaktionsordnung vorliegt. Nur für die Dimerisierung von Phenylacetylen mit einem Katalysatorsystem aus Palladium(II)-acetat und Triphenylphosphan findet sich eine Reaktion erster Ordnung.Vergleiche von Produktzusammensetzungen nach Ablauf einiger Homodimerisierungen bei verschiedenen Temperaturen führen zu dem Schluss, dass diese Reaktionen durch Temperaturveränderungen nicht kinetisch kontrollierbar sind.Die Struktur des 2,4-Bis(4-methoxyphenyl)-but-1-en-3-ins, das bei der Dimerisierung von 4-Methoxy-phenylacetylen entsteht, konnte durch eine Röntgenstrukturanalyse bestimmt werden.

Grundlegende Untersuchungen zur Abbindereaktion von Zinkphosphatzement

Zusammenfassung Das Ziel der Arbeit bestand darin, mit der Aufklärung der Abbindereaktion von Zinkphosphatzement eine Grundlage für eine gezielte Modifikation bzw. Optimierung zu schaffen, insbesondere im Hinblick auf einen Einsatz als permanenter Füllungswerkstoff. Über den Abbindechemismus war bislang lediglich bekannt, daß es sich bei den primär gebildeten Reaktionsprodukten um röntgenamorphe Phasen handelt, die sich nach Wochen bzw. Monaten in das thermodynamisch stabile Reaktionsprodukt alpha-Hopeit (alpha-Zn3(PO4)2·4H2O) umwandeln.Im Rahmen der vorliegenden Arbeit gelang durch den Einsatz der Infrarot-Reflexionsspek-troskopie (DRIFT) die Identifikation von Dizink Cyclotetraphosphat-Octahydrat (Zn2P4O12·8H2O) als röntgenamorpher Vorläuferphase. Das kondensierte Phosphat ist hydrolyseempfindlich sowie thermodynamisch instabil. Mit Hilfe der zeitaufgelösten 1H NMR-Spektroskopie konnte gezeigt werden, daß bereits nach ca. 10 Minuten eine Phasenumwandlung (topochemische Disproportionierung) in ein zunächst röntgenamorphes Orthophosphat stattfindet. Mittels 31P Doppelquanten-NMR-Spektroskopie gelang der Nachweis, daß innerhalb des röntgenamorphen Bereiches lokal nahgeordnete (nanokristalline) Bereiche auftreten, deren Ordnung sich auf einer Längenskala von 10 bis 30 Å erstreckt. Die Nanokristallite unterliegen einem Wachstumsprozeß, der schließlich zu alpha-Hopeit-Kristallen mit Ausdehnungen im Mikrometerbereich führt. Die Ursache für die primäre Ausbildung röntgenamorpher Reaktionsprodukte kann zunächst gelösten Aluminophosphatkomplexen zugeordnet werden, die im Verlauf der Abbindereaktion zu anorganischen Polymeren aggregieren und damit als Kristallisationsinhibitoren fungieren.

Doppelresonanzmessungen an in einer Penningfalle gespeicherten Europium-Ionen zur Bestimmung des Kern-g-Faktors

Die differentielle Hyperfeinanomalie beschreibt den Isotopieeffekt der magnetischen Hyperfeinwechselwirkung und stellt eine Testmöglichkeit für Kernmodelle dar. Nachdem in einer Reihe von Messungen die A-Faktoren der Hyperfeinwechselwirkung im Grundzustand von zwei stabilen und fünf instabilen Europium-Isotopen bestimmt wurden, sollen nun die Kern-g-Faktoren der gleichen Isotope bestimmt werden.Um die Kern-g-Faktoren zu bestimmen, wurden zunächst die Termschemata des Grundzustandesder stabilen Isotope 151,153Eu+ simuliert. Wegen der hohen Kern- und Hüllenspins I=5/2, J=4 ergeben sich 54 Zeemanzustände und ein komplexes optisches Spektrum.Etwa 10^6 Ionen der beiden stabilen Isotope 151,153Eu+ werden in einer Penningfalle beieinem Magnetfeld von ca. 1.5T gespeichert. Durch Puffergas und Anregung einer Seitenbandfrequenzder Ionenbewegung wird optisches Pumpen in metastabile D-Zustände bei gleichzeitiger Kühlung der Ionenwolke verhindert.Bei Verstimmung eines frequenzverdoppelten Ti:Sa-Lasers wurde die Ionenfluoreszenz aufgenommen. Durch die Simulation der optischen Spektren gelang es, einzelne Übergänge des Isotopes 151Eu+ zu identifizieren. Mit Laser/Mikrowellen-Doppelresonanzmessungen erhält man ein in erster Ordnung dopplerfreies Mikrowellenspektren eines Überganges zwischen zwei Zeemanzuständen im Grundzustand.Durch Messung der Übergangsfrequenzen von insgesamt fünf Delta mI=+1 Übergängen und Bestimmung des Magnetfeldes durch Messung der Zyklotronfrequenz gespeicherter Elektronen konnte der gI-Faktor des Isotopes151Eu zu 151gI = 1.37734(6) bestimmt werden.

Mixed Mode Fracture Behaviour of Piezoelectric Materials

Piezoelectrics present an interactive electromechanical behaviour that, especially in recent years, has generated much interest since it renders these materials adapt for use in a variety of electronic and industrial applications like sensors, actuators, transducers, smart structures. Both mechanical and electric loads are generally applied on these devices and can cause high concentrations of stress, particularly in proximity of defects or inhomogeneities, such as flaws, cavities or included particles. A thorough understanding of their fracture behaviour is crucial in order to improve their performances and avoid unexpected failures. Therefore, a considerable number of research works have addressed this topic in the last decades. Most of the theoretical studies on this subject find their analytical background in the complex variable formulation of plane anisotropic elasticity. This theoretical approach bases its main origins in the pioneering works of Muskelishvili and Lekhnitskii who obtained the solution of the elastic problem in terms of independent analytic functions of complex variables. In the present work, the expressions of stresses and elastic and electric displacements are obtained as functions of complex potentials through an analytical formulation which is the application to the piezoelectric static case of an approach introduced for orthotropic materials to solve elastodynamics problems. This method can be considered an alternative to other formalisms currently used, like the Stroh’s formalism. The equilibrium equations are reduced to a first order system involving a six-dimensional vector field. After that, a similarity transformation is induced to reach three independent Cauchy-Riemann systems, so justifying the introduction of the complex variable notation. Closed form expressions of near tip stress and displacement fields are therefore obtained. In the theoretical study of cracked piezoelectric bodies, the issue of assigning consistent electric boundary conditions on the crack faces is of central importance and has been addressed by many researchers. Three different boundary conditions are commonly accepted in literature: the permeable, the impermeable and the semipermeable (“exact”) crack model. This thesis takes into considerations all the three models, comparing the results obtained and analysing the effects of the boundary condition choice on the solution. The influence of load biaxiality and of the application of a remote electric field has been studied, pointing out that both can affect to a various extent the stress fields and the angle of initial crack extension, especially when non-singular terms are retained in the expressions of the electro-elastic solution. Furthermore, two different fracture criteria are applied to the piezoelectric case, and their outcomes are compared and discussed. The work is organized as follows: Chapter 1 briefly introduces the fundamental concepts of Fracture Mechanics. Chapter 2 describes plane elasticity formalisms for an anisotropic continuum (Eshelby-Read-Shockley and Stroh) and introduces for the simplified orthotropic case the alternative formalism we want to propose. Chapter 3 outlines the Linear Theory of Piezoelectricity, its basic relations and electro-elastic equations. Chapter 4 introduces the proposed method for obtaining the expressions of stresses and elastic and electric displacements, given as functions of complex potentials. The solution is obtained in close form and non-singular terms are retained as well. Chapter 5 presents several numerical applications aimed at estimating the effect of load biaxiality, electric field, considered permittivity of the crack. Through the application of fracture criteria the influence of the above listed conditions on the response of the system and in particular on the direction of crack branching is thoroughly discussed.

Development of analytical methodologies for iodine species in gaseous and particulate phases of the coastal atmosphere

It has been demonstrated that iodine does have an important influence on atmospheric chemistry, especially the formation of new particles and the enrichment of iodine in marine aerosols. It was pointed out that the most probable chemical species involved in the production or growth of these particles are iodine oxides, produced photochemically from biogenic halocarbon emissions and/or iodine emission from the sea surface. However, the iodine chemistry from gaseous to particulate phase in the coastal atmosphere and the chemical nature of the condensing iodine species are still not understood. A Tenax / Carbotrap adsorption sampling technique and a thermo-desorption / cryo-trap / GC-MS system has been further developed and improved for the volatile organic iodine species in the gas phase. Several iodo-hydrocarbons such as CH3I, C2H5I, CH2ICl, CH2IBr and CH2I2 etc., have been measured in samples from a calibration test gas source (standards), real air samples and samples from seaweeds / macro-algae emission experiments. A denuder sampling technique has been developed to characterise potential precursor compounds of coastal particle formation processes, such as molecular iodine in the gas phase. Starch, TMAH (TetraMethylAmmonium Hydroxide) and TBAH (TetraButylAmmonium Hydroxide) coated denuders were tested for their efficiencies to collect I2 at the inner surface, followed by a TMAH extraction and ICP/MS determination, adding tellurium as an internal standard. The developed method has been proved to be an effective, accurate and suitable process for I2 measurement in the field, with the estimated detection limit of ~0.10 ng∙L-1 for a sampling volume of 15 L. An H2O/TMAH-Extraction-ICP/MS method has been developed for the accurate and sensitive determination of iodine species in tropospheric aerosol particles. The particle samples were collected on cellulose-nitrate filters using conventional filter holders or on cellulose nitrate/tedlar-foils using a 5-stage Berner impactor for size-segregated particle analysis. The water soluble species as IO3- and I- were separated by anion exchanging process after water extraction. Non-water soluble species including iodine oxide and organic iodine were digested and extracted by TMAH. Afterwards the triple samples were analysed by ICP/MS. The detection limit for particulate iodine was determined to be 0.10~0.20 ng•m-3 for sampling volumes of 40~100 m3. The developed methods have been used in two field measurements in May 2002 and September 2003, at and around the Mace Head Atmospheric Research Station (MHARS) located at the west coast of Ireland. Elemental iodine as a precursor of the iodine chemistry in the coastal atmosphere, was determined in the gas phase at a seaweed hot-spot around the MHARS, showing I2 concentrations were in the range of 0~1.6 ng∙L-1 and indicating a positive correlation with the ozone concentration. A seaweed-chamber experiment performed at the field measurement station showed that the I2 emission rate from macro-algae was in the range of 0.019~0.022 ng•min-1•kg-1. During these experiments, nanometer-particle concentrations were obtained from the Scanning Mobility Particle Sizer (SMPS) measurements. Particle number concentrations were found to have a linear correlation with elemental iodine in the gas phase of the seaweeds chamber, showing that gaseous I2 is one of the important precursors of the new particle formation in the coastal atmosphere. Iodine contents in the particle phase were measured in both field campaigns at and around the field measurement station. Total iodine concentrations were found to be in the range of 1.0 ~ 21.0 ng∙m-3 in the PM2.5 samples. A significant correlation between the total iodine concentrations and the nanometer-particle number concentrations was observed. The particulate iodine species analysis indicated that iodide contents are usually higher than those of iodate in all samples, with ratios in the range of 2~5:1. It is possible that those water soluble iodine species are transferred through the sea-air interface into the particle phase. The ratio of water soluble (iodate + iodide) and non-water soluble species (probably iodine oxide and organic iodine compounds) was observed to be in the range of 1:1 to 1:2. It appears that higher concentrated non-water soluble species, as the products of the photolysis from the gas phase into the particle phase, can be obtained in those samples while the nucleation events occur. That supports the idea that iodine chemistry in the coastal boundary layer is linked with new particle formation events. Furthermore, artificial aerosol particles were formed from gaseous iodine sources (e.g. CH2I2) using a laboratory reaction-chamber experiment, in which the reaction constant of the CH2I2 photolysis was calculated to be based upon the first order reaction kinetic. The end products of iodine chemistry in the particle phase were identified and quantified.

Characterization of igneous terranes by zircon dating: implications for the UHP relicts occurrences and suture identification in the Central Rhodope, Northern Greece

One of the key for the understanding of an orogenic belt is the characterization of the terranes involved and the identification of the suture(s) separating crustal blocks: these are essential information for large-scale paleo-reconstructions. In addition, the structural relationships between the terranes involved in the collisional processes and the eventual UHP relicts may provide first order inputs to exhumation models of subducted rocks. The structure of the Rhodope Massif (northern Greece and southern Bulgaria) results from the stacking of high-grade nappes during a continental collision, which age is comprised between Latest-Jurassic and Early-Cenozoic. UHP and HP relicts, associated with oceanic and ultramafic material, suggest the presence of a dismembered suture zone within the massif. The location of this suture remains unclear; furthermore, up to now, the UHP and eclogitic localities represent isolated spots and no synthesis on their structural position within the massif has been proposed. The first aim of this work is to define the relationships between HP-UHP relicts, crustal blocks, shear zones and amphibolitic material. To achieve this objective, we characterized the accreted blocks in terms of protoliths ages of the orthogneisses mainly along two cross sections on the Greek part of the belt. Geochemical affinities of meta-igneous rocks served as a complementary tool for terrane characterization and geodynamic interpretation. Single-zircon Pb-Pb evaporation and zircon U-Pb SHRIMP dating of orthogneiss protoliths define two groups of intrusion-ages: Permo-Carboniferous and Late Jurassic-Early Cretaceous. Structurally, these two groups correspond to distinct units: the Late Jurassic gneissic complex overthrusts the one bearing the Permo-Carboniferous orthogneisses. Mylonites, eclogites, amphibolites of oceanic affinities, and UHP micaschists, mark a “melange” zone, intensively sheared towards the SW, which separates the two units. Thus, we interpret them as two distinct terranes, the Rhodope and Thracia terranes, separated by the Nestos suture. The correlation of our findings in northern Greece to the Bulgarian part of the Massif suggests a northern rooting of the Nestos Suture. This configuration results of the closure of a marginal oceanic basin of the Tethys system by a north-directed subduction. This interpretation is supported by the geochemical affinities of the orthogneisses: the Late-Jurassic igneous rocks formed by subduction-related magmatism, pprobably the same north-directed subduction that gave rise to the UHP metamorphism of the metasediments of the “melange” zone. It is noteworthy that the UHP-HP relicts seem to be restricted to the contact between the two terranes suggesting that the UHP relicts are exhumed only within the suture zone. Furthermore, the singularity of the suture suggests that the Late-Jurassic subduction explains the occurrence of UHP and eclogite relicts in the Central Rhodope despite the large age range previously attributed the UHP and/or HP stage.

Numerical and analytical approaches to strongly correlated electron systems

In this thesis we consider three different models for strongly correlated electrons, namely a multi-band Hubbard model as well as the spinless Falicov-Kimball model, both with a semi-elliptical density of states in the limit of infinite dimensions d, and the attractive Hubbard model on a square lattice in d=2. In the first part, we study a two-band Hubbard model with unequal bandwidths and anisotropic Hund's rule coupling (J_z-model) in the limit of infinite dimensions within the dynamical mean-field theory (DMFT). Here, the DMFT impurity problem is solved with the use of quantum Monte Carlo (QMC) simulations. Our main result is that the J_z-model describes the occurrence of an orbital-selective Mott transition (OSMT), in contrast to earlier findings. We investigate the model with a high-precision DMFT algorithm, which was developed as part of this thesis and which supplements QMC with a high-frequency expansion of the self-energy. The main advantage of this scheme is the extraordinary accuracy of the numerical solutions, which can be obtained already with moderate computational effort, so that studies of multi-orbital systems within the DMFT+QMC are strongly improved. We also found that a suitably defined Falicov-Kimball (FK) model exhibits an OSMT, revealing the close connection of the Falicov-Kimball physics to the J_z-model in the OSM phase. In the second part of this thesis we study the attractive Hubbard model in two spatial dimensions within second-order self-consistent perturbation theory. This model is considered on a square lattice at finite doping and at low temperatures. Our main result is that the predictions of first-order perturbation theory (Hartree-Fock approximation) are renormalized by a factor of the order of unity even at arbitrarily weak interaction (U->0). The renormalization factor q can be evaluated as a function of the filling n for 00, the q-factor vanishes, signaling the divergence of self-consistent perturbation theory in this limit. Thus we present the first asymptotically exact results at weak-coupling for the negative-U Hubbard model in d=2 at finite doping.

Single polymers at surfaces: adsorption and detachment

This thesis is concerned with the adsorption and detachment of polymers at planar, rigid surfaces. We have carried out a systematic investigation of adsorption of polymers using analytical techniques as well as Monte Carlo simulations with a coarse grained off-lattice bead spring model. The investigation was carried out in three stages. In the first stage the adsorption of a single multiblock AB copolymer on a solid surface was investigated by means of simulations and scaling analysis. It was shown that the problem could be mapped onto an effective homopolymer problem. Our main result was the phase diagram of regular multiblock copolymers which shows an increase in the critical adsorption potential of the substrate with decreasing size of blocks. We also considered the adsorption of random copolymers which was found to be well described within the annealed disorder approximation. In the next phase, we studied the adsorption kinetics of a single polymer on a flat, structureless surface in the regime of strong physisorption. The idea of a ’stem-flower’ polymer conformation and the mechanism of ’zipping’ during the adsorption process were used to derive a Fokker-Planck equation with reflecting boundary conditions for the time dependent probability distribution function (PDF) of the number of adsorbed monomers. The numerical solution of the time-dependent PDF obtained from a discrete set of coupled differential equations were shown to be in perfect agreement with Monte Carlo simulation results. Finally we studied force induced desorption of a polymer chain adsorbed on an attractive surface. We approached the problem within the framework of two different statistical ensembles; (i) by keeping the pulling force fixed while measuring the position of the polymer chain end, and (ii) by measuring the force necessary to keep the chain end at fixed distance above the adsorbing plane. In the first case we treated the problem within the framework of the Grand Canonical Ensemble approach and derived analytic expressions for the various conformational building blocks, characterizing the structure of an adsorbed linear polymer chain, subject to pulling force of fixed strength. The main result was the phase diagram of a polymer chain under pulling. We demonstrated a novel first order phase transformation which is dichotomic i.e. phase coexistence is not possible. In the second case, we carried out our study in the “fixed height” statistical ensemble where one measures the fluctuating force, exerted by the chain on the last monomer when a chain end is kept fixed at height h over the solid plane at different adsorption strength ε. The phase diagram in the h − ε plane was calculated both analytically and by Monte Carlo simulations. We demonstrated that in the vicinity of the polymer desorption transition a number of properties like fluctuations and probability distribution of various quantities behave differently, if h rather than the force, f, is used as an independent control parameter.

Higher-order perturbative relativistic corrections to energies and properties

Relativistic effects need to be considered in quantum-chemical calculations on systems including heavy elements or when aiming at high accuracy for molecules containing only lighter elements. In the latter case, consideration of relativistic effects via perturbation theory is an attractive option. Among the available techniques, Direct Perturbation Theory (DPT) in its lowest order (DPT2) has become a standard tool for the calculation of relativistic corrections to energies and properties.In this work, the DPT treatment is extended to the next order (DPT4). It is demonstrated that the DPT4 correction can be obtained as a second derivative of the energy with respect to the relativistic perturbation parameter. Accordingly, differentiation of a suitable Lagrangian, thereby taking into account all constraints on the wave function, provides analytic expressions for the fourth-order energy corrections. The latter have been implemented at the Hartree-Fock level and within second-order Møller-Plesset perturbaton theory using standard analytic second-derivative techniques into the CFOUR program package. For closed-shell systems, the DPT4 corrections consist of higher-order scalar-relativistic effects as well as spin-orbit corrections with the latter appearing here for the first time in the DPT series.Relativistic corrections are reported for energies as well as for first-order electrical properties and compared to results from rigorous four-component benchmark calculations in order to judge the accuracy and convergence of the DPT expansion for both the scalar-relativistic as well as the spin-orbit contributions. Additionally, the importance of relativistic effects to the bromine and iodine quadrupole-coupling tensors is investigated in a joint experimental and theoretical study concerning the rotational spectra of CH2BrF, CHBrF2, and CH2FI.

Periodicities in the Hodgkin-Huxley model and versions of this model with stochastic input

A field of computational neuroscience develops mathematical models to describe neuronal systems. The aim is to better understand the nervous system. Historically, the integrate-and-fire model, developed by Lapique in 1907, was the first model describing a neuron. In 1952 Hodgkin and Huxley [8] described the so called Hodgkin-Huxley model in the article “A Quantitative Description of Membrane Current and Its Application to Conduction and Excitation in Nerve”. The Hodgkin-Huxley model is one of the most successful and widely-used biological neuron models. Based on experimental data from the squid giant axon, Hodgkin and Huxley developed their mathematical model as a four-dimensional system of first-order ordinary differential equations. One of these equations characterizes the membrane potential as a process in time, whereas the other three equations depict the opening and closing state of sodium and potassium ion channels. The membrane potential is proportional to the sum of ionic current flowing across the membrane and an externally applied current. For various types of external input the membrane potential behaves differently. This thesis considers the following three types of input: (i) Rinzel and Miller [15] calculated an interval of amplitudes for a constant applied current, where the membrane potential is repetitively spiking; (ii) Aihara, Matsumoto and Ikegaya [1] said that dependent on the amplitude and the frequency of a periodic applied current the membrane potential responds periodically; (iii) Izhikevich [12] stated that brief pulses of positive and negative current with different amplitudes and frequencies can lead to a periodic response of the membrane potential. In chapter 1 the Hodgkin-Huxley model is introduced according to Izhikevich [12]. Besides the definition of the model, several biological and physiological notes are made, and further concepts are described by examples. Moreover, the numerical methods to solve the equations of the Hodgkin-Huxley model are presented which were used for the computer simulations in chapter 2 and chapter 3. In chapter 2 the statements for the three different inputs (i), (ii) and (iii) will be verified, and periodic behavior for the inputs (ii) and (iii) will be investigated. In chapter 3 the inputs are embedded in an Ornstein-Uhlenbeck process to see the influence of noise on the results of chapter 2.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Development of Dynamic Constraint Models for a Model Based Transient Calibration Process

Model based calibration has gained popularity in recent years as a method to optimize increasingly complex engine systems. However virtually all model based techniques are applied to steady state calibration. Transient calibration is by and large an emerging technology. An important piece of any transient calibration process is the ability to constrain the optimizer to treat the problem as a dynamic one and not as a quasi-static process. The optimized air-handling parameters corresponding to any instant of time must be achievable in a transient sense; this in turn depends on the trajectory of the same parameters over previous time instances. In this work dynamic constraint models have been proposed to translate commanded to actually achieved air-handling parameters. These models enable the optimization to be realistic in a transient sense. The air handling system has been treated as a linear second order system with PD control. Parameters for this second order system have been extracted from real transient data. The model has been shown to be the best choice relative to a list of appropriate candidates such as neural networks and first order models. The selected second order model was used in conjunction with transient emission models to predict emissions over the FTP cycle. It has been shown that emission predictions based on air-handing parameters predicted by the dynamic constraint model do not differ significantly from corresponding emissions based on measured air-handling parameters.

Using the Natural Period of a Structure as an Indicator of the Significance of Second-Order Effects

This study investigates the feasibility of predicting the momentamplification in beam-column elements of steel moment-resisting frames using the structure's natural period. Unlike previous methods, which perform moment-amplification on a story-by-story basis, this study develops and tests two models that aim to predict a global amplification factor indicative of the largest relevant instance of local moment amplification in the structure. To thisend, a variety of two-dimensional frames is investigated using first and secondorder finite element analysis. The observed moment amplification is then compared with the predicted amplification based on the structure's natural period, which is calculated by first-order finite element analysis. As a benchmark, design moment amplification factors are calculated for each story using the story stiffness approach, and serve to demonstrate the relativeconservatism and accuracy of the proposed models with respect to current practice in design. The study finds that the observed moment amplification factors may vastly exceed expectations when internal member stresses are initially very small. Where the internal stresses are small relative to the member capacities, thesecases are inconsequential for design. To qualify the significance of the observed amplification factors, two parameters are used: the second-order moment normalized to the plastic moment capacity, and the combined flexural and axial stress interaction equations developed by AISC