Simulation studies of correlation functions and relaxation in polymeric systems

We investigate the statics and dynamics of a glassy,non-entangled, short bead-spring polymer melt with moleculardynamics simulations. Temperature ranges from slightlyabove the mode-coupling critical temperature to the liquidregime where features of a glassy liquid are absent. Ouraim is to work out the polymer specific effects on therelaxation and particle correlation. We find the intra-chain static structure unaffected bytemperature, it depends only on the distance of monomersalong the backbone. In contrast, the distinct inter-chainstructure shows pronounced site-dependence effects at thelength-scales of the chain and the nearest neighbordistance. There, we also find the strongest temperaturedependence which drives the glass transition. Both the siteaveraged coupling of the monomer and center of mass (CM) andthe CM-CM coupling are weak and presumably not responsiblefor a peak in the coherent relaxation time at the chain'slength scale. Chains rather emerge as soft, easilyinterpenetrating objects. Three particle correlations arewell reproduced by the convolution approximation with theexception of model dependent deviations. In the spatially heterogeneous dynamics of our system weidentify highly mobile monomers which tend to follow eachother in one-dimensional paths forming ``strings''. Thesestrings have an exponential length distribution and aregenerally short compared to the chain length. Thus, arelaxation mechanism in which neighboring mobile monomersmove along the backbone of the chain seems unlikely.However, the correlation of bonded neighbors is enhanced. When liquids are confined between two surfaces in relativesliding motion kinetic friction is observed. We study ageneric model setup by molecular dynamics simulations for awide range of sliding speeds, temperatures, loads, andlubricant coverings for simple and molecular fluids. Instabilities in the particle trajectories are identified asthe origin of kinetic friction. They lead to high particlevelocities of fluid atoms which are gradually dissipatedresulting in a friction force. In commensurate systemsfluid atoms follow continuous trajectories for sub-monolayercoverings and consequently, friction vanishes at low slidingspeeds. For incommensurate systems the velocity probabilitydistribution exhibits approximately exponential tails. Weconnect this velocity distribution to the kinetic frictionforce which reaches a constant value at low sliding speeds. This approach agrees well with the friction obtaineddirectly from simulations and explains Amontons' law on themicroscopic level. Molecular bonds in commensurate systemslead to incommensurate behavior, but do not change thequalitative behavior of incommensurate systems. However,crossed chains form stable load bearing asperities whichstrongly increase friction.

Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.

Size-selective investigations of spectral and emission properties of small particles and nanotubes

In dieser Arbeit wurde die Elektronenemission von Nanopartikeln auf Oberflächen mittels spektroskopischen Photoelektronenmikroskopie untersucht. Speziell wurden metallische Nanocluster untersucht, als selbstorganisierte Ensembles auf Silizium oder Glassubstraten, sowie ferner ein Metall-Chalcogenid (MoS2) Nanoröhren-Prototyp auf Silizium. Der Hauptteil der Untersuchungen war auf die Wechselwirkung von fs-Laserstrahlung mit den Nanopartikeln konzentriert. Die Energie der Lichtquanten war kleiner als die Austrittsarbeit der untersuchten Proben, so dass Ein-Photonen-Photoemission ausgeschlossen werden konnte. Unsere Untersuchungen zeigten, dass ausgehend von einem kontinuierlichen Metallfilm bis hin zu Clusterfilmen ein anderer Emissionsmechanismus konkurrierend zur Multiphotonen-Photoemission auftritt und für kleine Cluster zu dominieren beginnt. Die Natur dieses neuen Mechanismus` wurde durch verschiedenartige Experimente untersucht. Der Übergang von einem kontinuierlichen zu einem Nanopartikelfilm ist begleitet von einer Zunahme des Emissionsstroms von mehr als eine Größenordnung. Die Photoemissions-Intensität wächst mit abnehmender zeitlicher Breite des Laserpulses, aber diese Abhängigkeit wird weniger steil mit sinkender Partikelgröße. Die experimentellen Resultate wurden durch verschiedene Elektronenemissions-Mechanismen erklärt, z.B. Multiphotonen-Photoemission (nPPE), thermionische Emission und thermisch unterstützte nPPE sowie optische Feldemission. Der erste Mechanismus überwiegt für kontinuierliche Filme und Partikel mit Größen oberhalb von mehreren zehn Nanometern, der zweite und dritte für Filme von Nanopartikeln von einer Größe von wenigen Nanometern. Die mikrospektroskopischen Messungen bestätigten den 2PPE-Emissionsmechanismus von dünnen Silberfilmen bei „blauer“ Laseranregung (hν=375-425nm). Das Einsetzen des Ferminiveaus ist relativ scharf und verschiebt sich um 2hν, wenn die Quantenenergie erhöht wird, wogegen es bei „roter“ Laseranregung (hν=750-850nm) deutlich verbreitert ist. Es zeigte sich, dass mit zunehmender Laserleistung die Ausbeute von niederenergetischen Elektronen schwächer zunimmt als die Ausbeute von höherenergetischen Elektronen nahe der Fermikante in einem Spektrum. Das ist ein klarer Hinweis auf eine Koexistenz verschiedener Emissionsmechanismen in einem Spektrum. Um die Größenabhängigkeit des Emissionsverhaltens theoretisch zu verstehen, wurde ein statistischer Zugang zur Lichtabsorption kleiner Metallpartikel abgeleitet und diskutiert. Die Elektronenemissionseigenschaften bei Laseranregung wurden in zusätzlichen Untersuchungen mit einer anderen Anregungsart verglichen, der Passage eines Tunnelstroms durch einen Metall-Clusterfilm nahe der Perkolationsschwelle. Die elektrischen und Emissionseigenschaften von stromtragenden Silberclusterfilmen, welche in einer schmalen Lücke (5-25 µm Breite) zwischen Silberkontakten auf einem Isolator hergestellt wurden, wurden zum ersten Mal mit einem Emissions-Elektronenmikroskop (EEM) untersucht. Die Elektronenemission beginnt im nicht-Ohmschen Bereich der Leitungsstrom-Spannungskurve des Clusterfilms. Wir untersuchten das Verhalten eines einzigen Emissionszentrums im EEM. Es zeigte sich, dass die Emissionszentren in einem stromleitenden Silberclusterfilm Punktquellen für Elektronen sind, welche hohe Emissions-Stromdichten (mehr als 100 A/cm2) tragen können. Die Breite der Energieverteilung der Elektronen von einem einzelnen Emissionszentrum wurde auf etwa 0.5-0.6 eV abgeschätzt. Als Emissionsmechanismus wird die thermionische Emission von dem „steady-state“ heißen Elektronengas in stromdurchflossenen metallischen Partikeln vorgeschlagen. Größenselektierte, einzelne auf Si-Substraten deponierte MoS2-Nanoröhren wurden mit einer Flugzeit-basierten Zweiphotonen-Photoemissions-Spektromikroskopie untersucht. Die Nanoröhren-Spektren wiesen bei fs-Laser Anregung eine erstaunlich hohe Emissionsintensität auf, deutlich höher als die SiOx Substratoberfläche. Dagegen waren die Röhren unsichtbar bei VUV-Anregung bei hν=21.2 eV. Eine ab-initio-Rechnung für einen MoS2-Slab erklärt die hohe Intensität durch eine hohe Dichte freier intermediärer Zustände beim Zweiphotonen-Übergang bei hν=3.1 eV.

Petrogenesis of dunite, wehrlite and websterite xenoliths from Kimberley, South Africa: Origin as mantle peridotites of cumulates

Dunite, wehrlite and websterite xenoliths occur amongst a large abundance of mantle xenoliths in kimberlites of the Kimberley cluster in South Africa. Up to know they have mostly been neglected. On the basis of texture, major and trace elements, oxygen isotopes as well as Re-Os isotope characteristics, they can be subdivided into two groups. A coarse-grained mantle peridotite group, comprising dunite, wehrlite and websterite xenoliths, that are similar to fertile peridotites and represent upper mantle assemblages that are differently influenced by mantle metasomatism. And a cumulate group, containing fine-grained Fe-rich dunite xenoliths that represent cumulates of flood basalt magmatism related to ~183 Ma Karoo and ~2.7 Ga Ventersdorp events in southern Africa. Dunite, wehrlite and websterite xenoliths have preserved a complex history of melt depletion and metasomatic re-enrichment events, which gives information about the different re-enrichment stages of the subcratonic lithospheric mantle and the spatial differences within the Kaapvaal craton upper mantle. Websterite xenoliths comprise orthopyroxene (40-85 Vol. %), clinopyroxene (5-42 Vol. %), garnet (4-10 Vol. %) and subordinately olivine, while dunite and wehrlite xenoliths contain predominantly olivine (65-100 Vol %) and subordinately orthopyroxene, clinopyroxene and garnet. High melt depletion and a dunitic to harzburgitic protolith composition are reflected by high forsterite (Fo90-92) and high olivine NiO contents (2800-5000 ppm) and high orthopyroxene Mg# (Mg/(Mg+Fe)) of 0.91-0.93. Re-depletion ages of predominantly 2.9 Ga reflect a minimum age of melt depletion. Melt depletion ceased in conjunction with collision of the Kimberley block with the Witwatersrand block ~2.9 Ga ago. Subduction related re-fertilisation of the previously depleted mantle xenoliths is documented by i) amoeboid textured orthopyroxene, clinopyroxene and garnet, which crystallized in schlieren along olivine grain boundaries, ii) high whole-rock SiO2, Al2O3, CaO, TiO2, FeO contents, iii) low oxygen isotope ratios in clinopyroxene and garnet of 4.8-5.4 ‰ and 4.7-5.3 ‰, respectively and iv) trace element compositions of wehrlitic clinopyroxene and garnet in equilibrium with high-pressure partial melts of eclogite. Trace element disequilibrium of orthopyroxene with clinopyroxene and garnet indicates a separate origin for orthopyroxene, on one side as primary mantle orthopyroxene in dunite and wehrlite xenoliths and on the other side as reaction product with Si-rich melts produced by partial melting of eclogite. This reaction triggered replacement of olivine by orthopyroxene in the surrounding mantle and produced the typical Si-rich composition of Kaapvaal mantle peridotites. Partial melting of eclogite at higher temperatures produced a second metasomatic melt with lower SiO2, but higher Al2O3, CaO, FeO, Ti, Zr, Hf and a low oxygen isotope ratio. This melt triggered clinopyroxene and locally garnet and rutile crystallization in percolation veins, replacing olivine and orthopyroxene in the Kaapvaal upper mantle. Additionally, websterite xenoliths have experienced late stage cryptic metasomatism by the host kimberlite melt, changing the trace element composition of clinopyroxene, orthopyroxene and garnet to different extent. Hence websterite and most fertile lherzolite xenoliths have experienced three metasomatic events: i) reaction with high-Si melt, ii) percolation of subduction related silica melt with lower SiO2 content and iii) cryptic metasomatism by kimberlite. In contrast, dunite and wehrlite xenoliths have only experienced the second metasomatic event. They represent mantle lithologies further away from metasomatising agents. The Fe-rich dunites comprise olivine neoblasts with subordinate olivine porphyroclasts and parallel-orientated needles of ilmenite, which may enclose spinel. The lower forsterite and NiO contents of olivine in Fe-rich dunites compared to mantle peridotite xenoliths (Fo87-89 vs. Fo93-95 and 1300-2800ppm vs. 2200-3900 ppm, respectively), rules out a restitic origin. Cr-rich spinels are remnants of the original cumulate mineralogy that survived a late stage metasomatic overprint related to the production of the host kimberlite, producing ilmenite and phlogopite in some samples. Olivine porphyroclasts and neoblasts have different trace element compositions, the latter having high Ti, V, Cr and Ni and low Zn, Zr and Nb contents, indicating contrasting origins for neoblasts and porphyroclasts. The dunites have high 187Os/188Os ratios (0.11-0.15) indicating young (Phanerozoic) model ages for most samples, whereas three samples show isotopic mixtures between Phanerozoic neoblasts and ancient porphyroclastic material. Most Fe-rich dunite xenoliths can be interpreted as cumulates of fractional crystallization of Karoo magmatism, whereas the porphyroclasts are interpreted to be remnants from the much earlier Archaean Ventersdorp magmatism.

Mixtures of colloidal rods and spheres in bulk and in confinement

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. Monte Carlo simulations were performed to study the phase diagram of such rod-polymer mixtures. The colloidal rods were modelled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while overlap of polymers and rods is forbidden. In this thesis the emphasis was on the depletion effects caused by the addition of spheres on the isotropic phase of rod-like particles. Although most of the present experimental studies consider systems close to or beyond the isotropic-nematic transition, the isotropic phase with depletion interactions turns out to be a not less interesting topic. First, the percolation problem was studied in canonical simulations of a system of hard rods and soft spheres, where the amount of depletant was kept low to prevent phase separation of the mixture. The lowering of the percolation threshold seen in experiment is confirmed to be due to the depletion interactions. The local changes in the structure of the fluid of rods, which were measured in the simulations, indicated that the depletion forces enhance local alignment and aggregation of the rods. Then, the phase diagram of isotropic-isotropic demixing of short spherocylinders was calculated using grand canonical ensemble simulations with successive umbrella sampling. Finite size scaling analysis allowed to estimate the location of the critical point. Also, estimates for the interfacial tension between the coexisting isotropic phases and analyses of its power-law behaviour on approach of the critical point are presented. The obtained phase diagram was compared to the predictions of the free volume theory. After an analysis of the bulk, the phase behaviour in confinement was studied. The critical point of gas-liquid demixing is shifted to higher concentrations of rods and smaller concentrations of spheres due to the formation of an orientationally ordered surface film. If the separation between the walls becomes very small, the critical point is shifted back to smaller concentrations of rods because the surface film breaks up. A method to calculate the contact angle of the liquid-gas interface with the wall is introduced and the wetting behaviour on the approach to the critical point is analysed.

Experimental and Numerical Analysis of Gas Forced Convection through Microtubes and Micro Heat Exchangers

The last decade has witnessed very fast development in microfabrication technologies. The increasing industrial applications of microfluidic systems call for more intensive and systematic knowledge on this newly emerging field. Especially for gaseous flow and heat transfer at microscale, the applicability of conventional theories developed at macro scale is not yet completely validated; this is mainly due to scarce experimental data available in literature for gas flows. The objective of this thesis is to investigate these unclear elements by analyzing forced convection for gaseous flows through microtubes and micro heat exchangers. Experimental tests have been performed with microtubes having various inner diameters, namely 750 m, 510 m and 170 m, over a wide range of Reynolds number covering the laminar region, the transitional zone and also the onset region of the turbulent regime. The results show that conventional theory is able to predict the flow friction factor when flow compressibility does not appear and the effect of fluid temperature-dependent properties is insignificant. A double-layered microchannel heat exchanger has been designed in order to study experimentally the efficiency of a gas-to-gas micro heat exchanger. This microdevice contains 133 parallel microchannels machined into polished PEEK plates for both the hot side and the cold side. The microchannels are 200 µm high, 200 µm wide and 39.8 mm long. The design of the micro device has been made in order to be able to test different materials as partition foil with flexible thickness. Experimental tests have been carried out for five different partition foils, with various mass flow rates and flow configurations. The experimental results indicate that the thermal performance of the countercurrent and cross flow micro heat exchanger can be strongly influenced by axial conduction in the partition foil separating the hot gas flow and cold gas flow.

Design and Characterization of Curved and Spherical Flexure Hinges for Planar and Spatial Compliant Mechanisms

A flexure hinge is a flexible connector that can provide a limited rotational motion between two rigid parts by means of material deformation. These connectors can be used to substitute traditional kinematic pairs (like bearing couplings) in rigid-body mechanisms. When compared to their rigid-body counterpart, flexure hinges are characterized by reduced weight, absence of backlash and friction, part-count reduction, but restricted range of motion. There are several types of flexure hinges in the literature that have been studied and characterized for different applications. In our study, we have introduced new types of flexures with curved structures i.e. circularly curved-beam flexures and spherical flexures. These flexures have been utilized for both planar applications (e.g. articulated robotic fingers) and spatial applications (e.g. spherical compliant mechanisms). We have derived closed-form compliance equations for both circularly curved-beam flexures and spherical flexures. Each element of the spatial compliance matrix is analytically computed as a function of hinge dimensions and employed material. The theoretical model is then validated by comparing analytical data with the results obtained through Finite Element Analysis. A case study is also presented for each class of flexures, concerning the potential applications in the optimal design of planar and spatial compliant mechanisms. Each case study is followed by comparing the performance of these novel flexures with the performance of commonly used geometries in terms of principle compliance factors, parasitic motions and maximum stress demands. Furthermore, we have extended our study to the design and analysis of serial and parallel compliant mechanisms, where the proposed flexures have been employed to achieve spatial motions e.g. compliant spherical joints.

Volatile organic compounds at the air-sea interface : gas exchange rates, oceanic emissions and the effect of ocean acidification

Oceans are key sources and sinks in the global budgets of significant atmospheric trace gases, termed Volatile Organic Compounds (VOCs). Despite their low concentrations, these species have an important role in the atmosphere, influencing ozone photochemistry and aerosol physics. Surprisingly, little work has been done on assessing their emissions or transport mechanisms and rates between ocean and atmosphere, all of which are important when modelling the atmosphere accurately.rnA new Needle Trap Device (NTD) - GC-MS method was developed for the effective sampling and analysis of VOCs in seawater. Good repeatability (RSDs <16 %), linearity (R2 = 0.96 - 0.99) and limits of detection in the range of pM were obtained for DMS, isoprene, benzene, toluene, p-xylene, (+)-α-pinene and (-)-α-pinene. Laboratory evaluation and subsequent field application indicated that the proposed method can be used successfully in place of the more usually applied extraction techniques (P&T, SPME) to extend the suite of species typically measured in the ocean and improve detection limits. rnDuring a mesocosm CO2 enrichment study, DMS, isoprene and α-pinene were identified and quantified in seawater samples, using the above mentioned method. Based on correlations with available biological datasets, the effects of ocean acidification as well as possible ocean biological sources were investigated for all examined compounds. Future ocean's acidity was shown to decrease oceanic DMS production, possibly impact isoprene emissions but not affect the production of α-pinene. rnIn a separate activity, ocean - atmosphere interactions were simulated in a large scale wind-wave canal facility, in order to investigate the gas exchange process and its controlling mechanisms. Air-water exchange rates of 14 chemical species (of which 11 VOCs) spanning a wide range of solubility (dimensionless solubility, α = 0:4 to 5470) and diffusivity (Schmidt number in water, Scw = 594 to 1194) were obtained under various turbulent (wind speed at ten meters height, u10 = 0:8 to 15ms-1) and surfactant modulated (two different sized Triton X-100 layers) surface conditions. Reliable and reproducible total gas transfer velocities were obtained and the derived values and trends were comparable to previous investigations. Through this study, a much better and more comprehensive understanding of the gas exchange process was accomplished. The role of friction velocity, uw* and mean square slope, σs2 in defining phenomena such as waves and wave breaking, near surface turbulence, bubbles and surface films was recognized as very significant. uw* was determined as the ideal turbulent parameter while σs2 described best the related surface conditions. A combination of both uw* and σs2 variables, was found to reproduce faithfully the air-water gas exchange process. rnA Total Transfer Velocity (TTV) model provided by a compilation of 14 tracers and a combination of both uw* and σs2 parameters, is proposed for the first time. Through the proposed TTV parameterization, a new physical perspective is presented which provides an accurate TTV for any tracer within the examined solubility range. rnThe development of such a comprehensive air-sea gas exchange parameterization represents a highly useful tool for regional and global models, providing accurate total transfer velocity estimations for any tracer and any sea-surface status, simplifying the calculation process and eliminating inevitable calculation uncertainty connected with the selection or combination of different parameterizations.rnrn

Prediction and measurement of forced response of a composite blade undergoing nonlinear vibration

The main objective of this project is to experimentally demonstrate geometrical nonlinear phenomena due to large displacements during resonant vibration of composite materials and to explain the problem associated with fatigue prediction at resonant conditions. Three different composite blades to be tested were designed and manufactured, being their difference in the composite layup (i.e. unidirectional, cross-ply, and angle-ply layups). Manual envelope bagging technique is explained as applied to the actual manufacturing of the components; problems encountered and their solutions are detailed. Forced response tests of the first flexural, first torsional, and second flexural modes were performed by means of a uniquely contactless excitation system which induced vibration by using a pulsed airflow. Vibration intensity was acquired by means of Polytec LDV system. The first flexural mode is found to be completely linear irrespective of the vibration amplitude. The first torsional mode exhibits a general nonlinear softening behaviour which is interestingly coupled with a hardening behaviour for the unidirectional layup. The second flexural mode has a hardening nonlinear behaviour for either the unidirectional and angle-ply blade, whereas it is slightly softening for the cross-ply layup. By using the same equipment as that used for forced response analyses, free decay tests were performed at different airflow intensities. Discrete Fourier Trasform over the entire decay and Sliding DFT were computed so as to visualise the presence of nonlinear superharmonics in the decay signal and when they were damped out from the vibration over the decay time. Linear modes exhibit an exponential decay, while nonlinearities are associated with a dry-friction damping phenomenon which tends to increase with increasing amplitude. Damping ratio is derived from logarithmic decrement for the exponential branch of the decay.

Experimental characterization and modelling of a servo-pneumatic system for a knee loading apparatus

The new knee test rig developed in University of Bologna used pneumatic cylinder as actuator system. Specific characterization and modelling about the pneumatic cylinder and the related devices are needed in better controlling the test rig. In this thesis, an experimental environment for the related device is set up with data acquisition system using Real-time Windows Target, Simulink, MatLab. Based on the experimental data, a fitted model for the pneumatic cylinder friction is found.

Equations for the Drag Force and Aerodynamic Roughness Length of Urban Areas with Random Building Heights

We use a conceptual model to investigate how randomly varying building heights within a city affect the atmospheric drag forces and the aerodynamic roughness length of the city. The model is based on the assumptions regarding wake spreading and mutual sheltering effects proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992). It is applied both to canopies having uniform building heights and to those having the same building density and mean height, but with variability about the mean. For each simulated urban area, a correction is determined, due to height variability, to the shear stress predicted for the uniform building height case. It is found that u (*)/u (*R) , where u (*) is the friction velocity and u (*R) is the friction velocity from the uniform building height case, is expressed well as an algebraic function of lambda and sigma (h) /h (m) , where lambda is the frontal area index, sigma (h) is the standard deviation of the building height, and h (m) is the mean building height. The simulations also resulted in a simple algebraic relation for z (0)/z (0R) as a function of lambda and sigma (h) /h (m) , where z (0) is the aerodynamic roughness length and z (0R) is z (0) found from the original Raupach formulation for a uniform canopy. Model results are in keeping with those of several previous studies.

Modeling And Predicting The Performance Of Infiltrated Solid Oxide Fuel Cell Anodes

Fuel cells are a topic of high interest in the scientific community right now because of their ability to efficiently convert chemical energy into electrical energy. This thesis is focused on solid oxide fuel cells (SOFCs) because of their fuel flexibility, and is specifically concerned with the anode properties of SOFCs. The anodes are composed of a ceramic material (yttrium stabilized zirconia, or YSZ), and conducting material. Recent research has shown that an infiltrated anode may offer better performance at a lower cost. This thesis focuses on the creation of a model of an infiltrated anode that mimics the underlying physics of the production process. Using the model, several key parameters for anode performance are considered. These are the initial volume fraction of YSZ in the slurry before sintering, the final porosity of the composite anode after sintering, and the size of the YSZ and conducting particles in the composite. The performance measures of the anode, namely percolation threshold and effective conductivity, are analyzed as a function of these important input parameters. Simple two and three-dimensional percolation models are used to determine the conditions at which the full infiltrated anode would be investigated. These more simple models showed that the aspect ratio of the anode has no effect on the threshold or effective conductivity, and that cell sizes of 303 are needed to obtain accurate conductivity values. The full model of the infiltrated anode is able to predict the performance of the SOFC anodes and it can be seen that increasing the size of the YSZ decreases the percolation threshold and increases the effective conductivity at low conductor loadings. Similar trends are seen for a decrease in final porosity and a decrease in the initial volume fraction of YSZ.

The Effect of Various Electrode Microstructure Parameters on the Effective Conductivity and Three-Phase Boundary Density of Infiltrated SOFC Electrodes

Solid oxide fuel cell (SOFC) technology has the potential to be a significant player in our future energy technology repertoire based on its ability to convert chemical energy into electrical energy. Infiltrated SOFCs, in particular, have demonstrated improved performance and at lower cost than traditional SOFCs. An infiltrated electrode comprises porous ceramic scaffolding (typically constructed from the oxygen ion conducting material) that is infiltrated with electron conducting and catalytic particles. Two important SOFC electrode properties are effective conductivity and three phase boundary density (TPB). Researchers study these electrode properties separately, and fail to recognize them as competing properties. This thesis aims to (1) develop a method to model the TPB density and use it to determine the effect of porosity, scaffolding particle size, and pore former size on TPB density as well as to (2) compare the effect of porosity, scaffolding particle size, and pore former size on TPB density and effective conductivity to determine a desired set of parameters for infiltrated SOFC electrode performance. A computational model was used to study the effect of microstructure parameters on the effective conductivity and TPB density of the infiltrated SOFC electrode. From this study, effective conductivity and TPB density are determined to be competing properties of SOFC electrodes. Increased porosity, scaffolding particle size, and pore former particle size increase the effective conductivity for a given infiltrate loading above percolation threshold. Increased scaffolding particle size and pore former size ratio, however, decreases the TPB density. The maximum TPB density is achievable between porosities of 45% and 60%. The effect of microstructure parameters are more prominent at low loading with scaffolding particle size being the most significant factor and pore former size ratio being the least significant factor.

Analyzing the effects of vertical acceleration and deceleration in multi-phase, debris-flow sediment-column experiments / by Christine Marie Williams.

Water-saturated debris flows are among some of the most destructive mass movements. Their complex nature presents a challenge for quantitative description and modeling. In order to improve understanding of the dynamics of these flows, it is important to seek a simplified dynamic system underlying their behavior. Models currently in use to describe the motion of debris flows employ depth-averaged equations of motion, typically assuming negligible effects from vertical acceleration. However, in many cases debris flows experience significant vertical acceleration as they move across irregular surfaces, and it has been proposed that friction associated with vertical forces and liquefaction merit inclusion in any comprehensive mechanical model. The intent of this work is to determine the effect of vertical acceleration through a series of laboratory experiments designed to simulate debris flows, testing a recent model for debris flows experimentally. In the experiments, a mass of water-saturated sediment is released suddenly from a holding container, and parameters including rate of collapse, pore-fluid pressure, and bed load are monitored. Experiments are simplified to axial geometry so that variables act solely in the vertical dimension. Steady state equations to infer motion of the moving sediment mass are not sufficient to model accurately the independent solid and fluid constituents in these experiments. The model developed in this work more accurately predicts the bed-normal stress of a saturated sediment mass in motion and illustrates the importance of acceleration and deceleration.

Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.