Characterization and modeling of the barrier properties in nanostructured systems

The object of the present study is the process of gas transport in nano-sized materials, i.e. systems having structural elements of the order of nanometers. The aim of this work is to advance the understanding of the gas transport mechanism in such materials, for which traditional models are not often suitable, by providing a correct interpretation of the relationship between diffusive phenomena and structural features. This result would allow the development new materials with permeation properties tailored on the specific application, especially in packaging systems. The methods used to achieve this goal were a detailed experimental characterization and different simulation methods. The experimental campaign regarded the determination of oxygen permeability and diffusivity in different sets of organic-inorganic hybrid coatings prepared via sol-gel technique. The polymeric samples coated with these hybrid layers experienced a remarkable enhancement of the barrier properties, which was explained by the strong interconnection at the nano-scale between the organic moiety and silica domains. An analogous characterization was performed on microfibrillated cellulose films, which presented remarkable barrier effect toward oxygen when it is dry, while in the presence of water the performance significantly drops. The very low value of water diffusivity at low activities is also an interesting characteristic which deals with its structural properties. Two different approaches of simulation were then considered: the diffusion of oxygen through polymer-layered silicates was modeled on a continuum scale with a CFD software, while the properties of n-alkanthiolate self assembled monolayers on gold were analyzed from a molecular point of view by means of a molecular dynamics algorithm. Modeling transport properties in layered nanocomposites, resulting from the ordered dispersion of impermeable flakes in a 2-D matrix, allowed the calculation of the enhancement of barrier effect in relation with platelets structural parameters leading to derive a new expression. On this basis, randomly distributed systems were simulated and the results were analyzed to evaluate the different contributions to the overall effect. The study of more realistic three-dimensional geometries revealed a prefect correspondence with the 2-D approximation. A completely different approach was applied to simulate the effect of temperature on the oxygen transport through self assembled monolayers; the structural information obtained from equilibrium MD simulations showed that raising the temperature, makes the monolayer less ordered and consequently less crystalline. This disorder produces a decrease in the barrier free energy and it lowers the overall resistance to oxygen diffusion, making the monolayer more permeable to small molecules.

Analisi della dinamica e composizione chimica del particolato atmosferico mediante un modello alla mesoscala

The vertical profile of aerosol in the planetary boundary layer of the Milan urban area is studied in terms of its development and chemical composition in a high-resolution modelling framework. The period of study spans a week in summer of 2007 (12-18 July), when continuous LIDAR measurements and a limited set of balloon profiles were collected in the frame of the ASI/QUITSAT project. LIDAR observations show a diurnal development of an aerosol plume that lifts early morning surface emissions to the top of the boundary layer, reaching maximum concentration around midday. Mountain breeze from Alps clean the bottom of the aerosol layer, typically leaving a residual layer at around 1500-2000 m which may survive for several days. During the last two days under analysis, a dust layer transported from Sahara reaches the upper layers of Milan area and affects the aerosol vertical distribution in the boundary layer. Simulation from the MM5/CHIMERE modelling system, carried out at 1 km horizontal resolution, qualitatively reproduced the general features of the Milan aerosol layer observed with LIDAR, including the rise and fall of the aersol plume, the residual layer in altitude and the Saharan dust event. The simulation highlighted the importance of nitrates and secondary organics in its composition. Several sensitivity tests showed that main driving factors leading to the dominance of nitrates in the plume are temperature and gas absorption process. A modelling study turn to the analysis of the vertical aerosol profiles distribution and knowledge of the characterization of the PM at a site near the city of Milan is performed using a model system composed by a meteorological model MM5 (V3-6), the mesoscale model from PSU/NCAR and a Chemical Transport Model (CTM) CHIMERE to simulate the vertical aerosol profile. LiDAR continuous observations and balloon profiles collected during two intensive campaigns in summer 2007 and in winter 2008 in the frame of the ASI/QUITSAT project have been used to perform comparisons in order to evaluate the ability of the aerosol chemistry transport model CHIMERE to simulate the aerosols dynamics and compositions in this area. The comparisons of model aerosols with measurements are carried out over a full time period between 12 July 2007 and 18 July 2007. The comparisons demonstrate the ability of the model to reproduce correctly the aerosol vertical distributions and their temporal variability. As detected by the LiDAR, the model during the period considered, predicts a diurnal development of a plume during the morning and a clearing during the afternoon, typically the plume reaches the top of the boundary layer around mid day, in this time CHIMERE produces highest concentrations in the upper levels as detected by LiDAR. The model, moreover can reproduce LiDAR observes enhancement aerosols concentrations above the boundary layer, attributing the phenomena to dust out intrusion. Another important information from the model analysis regard the composition , it predicts that a large part of the plume is composed by nitrate, in particular during 13 and 16 July 2007 , pointing to the model tendency to overestimates the nitrous component in the particular matter vertical structure . Sensitivity study carried out in this work show that there are a combination of different factor which determine the major nitrous composition of the “plume” observed and in particular humidity temperature and the absorption phenomena are the mainly candidate to explain the principal difference in composition simulated in the period object of this study , in particular , the CHIMERE model seems to be mostly sensitive to the absorption process.

Interactive simulation of flexible parts

Computer simulations play an ever growing role for the development of automotive products. Assembly simulation, as well as many other processes, are used systematically even before the first physical prototype of a vehicle is built in order to check whether particular components can be assembled easily or whether another part is in the way. Usually, this kind of simulation is limited to rigid bodies. However, a vehicle contains a multitude of flexible parts of various types: cables, hoses, carpets, seat surfaces, insulations, weatherstrips... Since most of the problems using these simulations concern one-dimensional components and since an intuitive tool for cable routing is still needed, we have chosen to concentrate on this category, which includes cables, hoses and wiring harnesses. In this thesis, we present a system for simulating one dimensional flexible parts such as cables or hoses. The modeling of bending and torsion follows the Cosserat model. For this purpose we use a generalized spring-mass system and describe its configuration by a carefully chosen set of coordinates. Gravity and contact forces as well as the forces responsible for length conservation are expressed in Cartesian coordinates. But bending and torsion effects can be dealt with more effectively by using quaternions to represent the orientation of the segments joining two neighboring mass points. This augmented system allows an easy formulation of all interactions with the best appropriate coordinate type and yields a strongly banded Hessian matrix. An energy minimizing process accounts for a solution exempt from the oscillations that are typical of spring-mass systems. The use of integral forces, similar to an integral controller, allows to enforce exactly the constraints. The whole system is numerically stable and can be solved at interactive frame rates. It is integrated in the DaimlerChrysler in-house Virtual Reality Software veo for use in applications such as cable routing and assembly simulation and has been well received by users. Parts of this work have been published at the ACM Solid and Physical Modeling Conference 2006 and have been selected for the special issue of the Computer-Aided-Design Journal to the conference.

Nonlinear Control of Magnetically Actuated Spacecraft

The topic of this thesis is the feedback stabilization of the attitude of magnetically actuated spacecraft. The use of magnetic coils is an attractive solution for the generation of control torques on small satellites flying inclined low Earth orbits, since magnetic control systems are characterized by reduced weight and cost, higher reliability, and require less power with respect to other kinds of actuators. At the same time, the possibility of smooth modulation of control torques reduces coupling of the attitude control system with flexible modes, thus preserving pointing precision with respect to the case when pulse-modulated thrusters are used. The principle based on the interaction between the Earth's magnetic field and the magnetic field generated by the set of coils introduces an inherent nonlinearity, because control torques can be delivered only in a plane that is orthogonal to the direction of the geomagnetic field vector. In other words, the system is underactuated, because the rotational degrees of freedom of the spacecraft, modeled as a rigid body, exceed the number of independent control actions. The solution of the control issue for underactuated spacecraft is also interesting in the case of actuator failure, e.g. after the loss of a reaction-wheel in a three-axes stabilized spacecraft with no redundancy. The application of well known control strategies is no longer possible in this case for both regulation and tracking, so that new methods have been suggested for tackling this particular problem. The main contribution of this thesis is to propose continuous time-varying controllers that globally stabilize the attitude of a spacecraft, when magneto-torquers alone are used and when a momentum-wheel supports magnetic control in order to overcome the inherent underactuation. A kinematic maneuver planning scheme, stability analyses, and detailed simulation results are also provided, with new theoretical developments and particular attention toward application considerations.

Mode-coupling theory: generalizations, high dimensions and microscopic dynamics

This work contains several applications of the mode-coupling theory (MCT) and is separated into three parts. In the first part we investigate the liquid-glass transition of hard spheres for dimensions d→∞ analytically and numerically up to d=800 in the framework of MCT. We find that the critical packing fraction ϕc(d) scales as d²2^(-d), which is larger than the Kauzmann packing fraction ϕK(d) found by a small-cage expansion by Parisi and Zamponi [J. Stat. Mech.: Theory Exp. 2006, P03017 (2006)]. The scaling of the critical packing fraction is different from the relation ϕc(d)∼d2^(-d) found earlier by Kirkpatrick and Wolynes [Phys. Rev. A 35, 3072 (1987)]. This is due to the fact that the k dependence of the critical collective and self nonergodicity parameters fc(k;d) and fcs(k;d) was assumed to be Gaussian in the previous theories. We show that in MCT this is not the case. Instead fc(k;d) and fcs(k;d), which become identical in the limit d→∞, converge to a non-Gaussian master function on the scale k∼d^(3/2). We find that the numerically determined value for the exponent parameter λ and therefore also the critical exponents a and b depend on the dimension d, even at the largest evaluated dimension d=800. In the second part we compare the results of a molecular-dynamics simulation of liquid Lennard-Jones argon far away from the glass transition [D. Levesque, L. Verlet, and J. Kurkijärvi, Phys. Rev. A 7, 1690 (1973)] with MCT. We show that the agreement between theory and computer simulation can be improved by taking binary collisions into account [L. Sjögren, Phys. Rev. A 22, 2866 (1980)]. We find that an empiric prefactor of the memory function of the original MCT equations leads to similar results. In the third part we derive the equations for a mode-coupling theory for the spherical components of the stress tensor. Unfortunately it turns out that they are too complex to be solved numerically.

Closed-loop investigation of the dynamical properties of cardiomyocyte electrophysiology by means of Dynamic clamp and Mathematical modelling

The research field of my PhD concerns mathematical modeling and numerical simulation, applied to the cardiac electrophysiology analysis at a single cell level. This is possible thanks to the development of mathematical descriptions of single cellular components, ionic channels, pumps, exchangers and subcellular compartments. Due to the difficulties of vivo experiments on human cells, most of the measurements are acquired in vitro using animal models (e.g. guinea pig, dog, rabbit). Moreover, to study the cardiac action potential and all its features, it is necessary to acquire more specific knowledge about single ionic currents that contribute to the cardiac activity. Electrophysiological models of the heart have become very accurate in recent years giving rise to extremely complicated systems of differential equations. Although describing the behavior of cardiac cells quite well, the models are computationally demanding for numerical simulations and are very difficult to analyze from a mathematical (dynamical-systems) viewpoint. Simplified mathematical models that capture the underlying dynamics to a certain extent are therefore frequently used. The results presented in this thesis have confirmed that a close integration of computational modeling and experimental recordings in real myocytes, as performed by dynamic clamp, is a useful tool in enhancing our understanding of various components of normal cardiac electrophysiology, but also arrhythmogenic mechanisms in a pathological condition, especially when fully integrated with experimental data.

Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.

Synthesis and dynamic study of atropisomeric compounds containing boron-carbon bond

In this thesis we studied the stereodynamic behavior of 1,2-azaborines variously substituted on boron (7a, 7b, 13). Depending on the hindrance of the asymmetric aryl substituent the resulting conformations could be stereolabile or configurationally stable. Through dynamic NMR and lineshape simulation, the energy rotational barriers of the different conformers are obtained. When the barrier is higher than 22-23 kcal/mol stable atropisomers that are fisically separable could be obtained (case of compound 13) and the free activation energy barrier is determinable by kinetic analysis. Absolute configuration of two atropisomers were assigned by comparison between computational calculations and experimental ECD. Isosteric compound 21 is then synthesized in order to compare the rotational barrier around B-Caryl with the one around Cnaphth-Caryl bond.

Reservoir characterization of the Cardium Formation in the Ferrier Oilfield, west-central Alberta, Canada.

This thesis has the main aim of defining the lithostratigraphy, depositional architecture, post-depositional modifications and reservoir characteristics of the Cardium Formation in the Ferrier Oilfield, and how these characteristics can have great impact over production rates, GOR and produced fluid discrimination. In the Ferrier area, the Cardium Formation is composed by a NE prograding clastic sequence made up of offshore to shoreface deposits sealed by marine shales. The main reservoir is composed by sandstones and conglomerates interpreted to have deposited in a shoreface depositional environment. Lithofacies and net reservoir thickness mapping led to more detailed understanding of the 3D reservoir architecture, and cross-sections shed light on the Cardium depositional architecture and post-deposition sediment erosion in the Ferrier area. Detailed core logging, thin section, SEM and CL analyses were used to study the mineralogy, texture and pore characterization of the Cardium reservoir, and three main compartments have been identified based on production data and reservoir characteristics. Finally, two situations showing odd production behaviour of the Cardium were resolved. This shed light on the effect of structural features and reservoir quality and thickness over hydrocarbon migration pathways. The Ferrier example offers a unique case of fluid discrimination in clastic reservoirs due both to depositional and post-depositional factors, and could be used as analogue for similar situations in the Western Canadian Sedimentary Basin.

Teeth restored using fiber-reinforced posts: in vitro fracture tests and finite element analysis

In dentistry the restoration of decayed teeth is challenging and makes great demands on both the dentist and the materials. Hence, fiber-reinforced posts have been introduced. The effects of different variables on the ultimate load on teeth restored using fiber-reinforced posts is controversial, maybe because the results are mostly based on non-standardized in vitro tests and, therefore, give inhomogeneous results. This study combines the advantages of in vitro tests and finite element analysis (FEA) to clarify the effects of ferrule height, post length and cementation technique used for restoration. Sixty-four single rooted premolars were decoronated (ferrule height 1 or 2 mm), endodontically treated and restored using fiber posts (length 2 or 7 mm), composite fillings and metal crowns (resin bonded or cemented). After thermocycling and chewing simulation the samples were loaded until fracture, recording first damage events. Using UNIANOVA to analyze recorded fracture loads, ferrule height and cementation technique were found to be significant, i.e. increased ferrule height and resin bonding of the crown resulted in higher fracture loads. Post length had no significant effect. All conventionally cemented crowns with a 1-mm ferrule height failed during artificial ageing, in contrast to resin-bonded crowns (75% survival rate). FEA confirmed these results and provided information about stress and force distribution within the restoration. Based on the findings of in vitro tests and computations we concluded that crowns, especially those with a small ferrule height, should be resin bonded. Finally, centrally positioned fiber-reinforced posts did not contribute to load transfer as long as the bond between the tooth and composite core was intact.

An insulin infusion advisory system based on autotuning nonlinear model-predictive control

This paper aims at the development and evaluation of a personalized insulin infusion advisory system (IIAS), able to provide real-time estimations of the appropriate insulin infusion rate for type 1 diabetes mellitus (T1DM) patients using continuous glucose monitors and insulin pumps. The system is based on a nonlinear model-predictive controller (NMPC) that uses a personalized glucose-insulin metabolism model, consisting of two compartmental models and a recurrent neural network. The model takes as input patient's information regarding meal intake, glucose measurements, and insulin infusion rates, and provides glucose predictions. The predictions are fed to the NMPC, in order for the latter to estimate the optimum insulin infusion rates. An algorithm based on fuzzy logic has been developed for the on-line adaptation of the NMPC control parameters. The IIAS has been in silico evaluated using an appropriate simulation environment (UVa T1DM simulator). The IIAS was able to handle various meal profiles, fasting conditions, interpatient variability, intraday variation in physiological parameters, and errors in meal amount estimations.

Influence of curing rate on softening in ethanol, degree of conversion, and wear of resin composite

PURPOSE: To investigate the effect of curing rate on softening in ethanol, degree of conversion, and wear of resin composites. METHOD: With a given energy density and for each of two different light-curing units (QTH or LED), the curing rate was reduced by modulating the curing mode. Thus, the irradiation of resin composite specimens (Filtek Z250, Tetric Ceram, Esthet-X) was performed in a continuous curing mode and in a pulse-delay curing mode. Wallace hardness was used to determine the softening of resin composite after storage in ethanol. Degree of conversion was determined by infrared spectroscopy (FTIR). Wear was assessed by a three-body test. Data were submitted to Levene's test, one and three-way ANOVA, and Tukey HSD test (alpha = 0.05). Results: Immersion in ethanol, curing mode, and material all had significant effects on Wallace hardness. After ethanol storage, resin composites exposed to the pulse-delay curing mode were softer than resin composites exposed to continuous cure (P< 0.0001). Tetric Ceram was the softest material followed by Esthet-X and Filtek Z250 (P< 0.001). Only the restorative material had a significant effect on degree of conversion (P< 0.001): Esthet-X had the lowest degree of conversion followed by Filtek Z250 and Tetric Ceram. Curing mode (P= 0.007) and material (P< 0.001) had significant effect on wear. Higher wear resulted from the pulse-delay curing mode when compared to continuous curing, and Filtek Z250 showed the lowest wear followed by Esthet-X and Tetric Ceram.

A simple, economical, and effective portable paediatric mock circulatory system

Ventricular assist devices (VADs) and total artificial hearts have been in development for the last 50 years. Since their inception, simulators of the circulation with different degrees of complexity have been produced to test these devices in vitro. Currently, a new path has been taken with the extensive efforts to develop paediatric VADs, which require totally different design constraints. This paper presents the manufacturing details of an economical simulator of the systemic paediatric circulation. This simulator allows the insertion of a paediatric VAD, includes a pumping ventricle, and is adjustable within the paediatric range. Rather than focusing on complexity and physiological simulation, this simulator is designed to be simple and practical for rapid device testing. The simulator was instrumented with medical sensors and data were acquired under different conditions with and without the new PediaFlowTM paediatric VAD. The VAD was run at different impeller speeds while simulator settings such as vascular resistance and stroke volume were varied. The hydraulic performance of the VAD under pulsatile conditions could be characterized and the magnetic suspension could be tested via manipulations such as cannula clamping. This compact mock loop has proven to be valuable throughout the PediaFlow development process and has the advantage that it is uncomplicated and can be manufactured cheaply. It can be produced by several research groups and the results of different VADs can then be compared easily.

Cell therapies for heart function recovery: focus on myocardial tissue engineering and nanotechnologies

Cell therapies have gained increasing interest and developed in several approaches related to the treatment of damaged myocardium. The results of multiple clinical trials have already been reported, almost exclusively involving the direct injection of stem cells. It has, however, been postulated that the efficiency of injected cells could possibly be hindered by the mechanical trauma due to the injection and their low survival in the hostile environment. It has indeed been demonstrated that cell mortality due to the injection approaches 90%. Major issues still need to be resolved and bed-to-bench followup is paramount to foster clinical implementations. The tissue engineering approach thus constitutes an attractive alternative since it provides the opportunity to deliver a large number of cells that are already organized in an extracellular matrix. Recent laboratory reports confirmed the interest of this approach and already encouraged a few groups to investigate it in clinical studies. We discuss current knowledge regarding engineered tissue for myocardial repair or replacement and in particular the recent implementation of nanotechnological approaches.