Molecular dynamics simulations of sheared brush-like systems

Mit Hilfe von Molekulardynamik-Simulationen untersuchen wir bürstenartige Systeme unter guten Lösungsmittelbedingungen. Diese Systeme sind, dank ihren vielfältigen Beschaffenheiten, die von Molekularparametern und äußeren Bedingungen abhängig sind, wichtig für viele industrielle Anwendungen. Man vermutet, dass die Polymerbürsten eine entscheidende Rolle in der Natur wegen ihrer einzigartigen Gleiteigenschaften spielen. Ein vergröbertes Modell wird verwendet, um die strukturellen und dynamischen Eigenschaften zweier hochkomprimierter Polymerbürsten, die eine niedrige Reibung aufweisen, zu untersuchen. Allerdings sind die Lubrikationseigenschaften dieser Systeme, die in vielen biologischen Systemen vorhanden sind, beeinflußt. Wir untersuchen so-genannte "weiche Kolloide", die zwischen den beiden Polymerbürsten eingebettet sind, und wie diese Makroobjekte auf die Polymerbürsten wirken.rnrnNicht-Gleichgewichts-Molekulardynamik-Simulationen werden durchgeführt, in denen die hydrodynamischen Wechselwirkungen durch die Anwendung des DPD-Thermostaten mit expliziten Lösungsmittelmolekülen berücksichtigt werden. Wir zeigen, dass die Kenntnis der Gleichgewichtseigenschaften des Systems erlaubt, dynamische Nichtgleichgewichtsigenschaften der Doppelschicht vorherzusagen.rnrnWir untersuchen, wie die effektive Wechselwirkung zwischen kolloidalen Einschlüßen durch die Anwesenheit der Bürsten (in Abhängigkeit der Weichheit der Kolloide und der Pfropfdichte der Bürsten) beeinflußt wird. Als nächsten Schritt untersuchen wir die rheologische Antwort von solchen komplexen Doppelschichten auf Scherung. Wir entwickeln eine Skalen-Theorie, die die Abhängigkeit der makroskopischen Transporteigenschaften und der lateralen Ausdehnung der verankerten Ketten von der Weissenberg Zahl oberhalb des Bereichs, in dem die lineare Antwort-Theorie gilt, voraussagt. Die Vorhersagen der Theorie stimmen gut mit unseren und früheren numerischen Ergebnissen und neuen Experimenten überein. Unsere Theorie bietet die Möglichkeit, die Relaxationszeit der Doppelschicht zu berechnen. Wenn diese Zeit mit einer charakteristischen Längenskala kombiniert wird, kann auch das ''transiente'' (nicht-stationäre) Verhalten beschrieben werden.rnrnrnWir untersuchen die Antwort des Drucktensors und die Deformation der Bürsten während der Scherinvertierung für grosse Weissenberg Zahlen. Wir entwickeln eine Vorhersage für die charakteristische Zeit, nach der das System wieder den stationären Zustand erreicht.rnrnrnElektrostatik spielt eine bedeutende Rolle in vielen biologischen Prozessen. Die Lubrikationseigenschaften der Polymerbürsten werden durch die Anwesenheit langreichweitiger Wechselwirkungen stark beeinflusst. Für unterschiedliche Stärken der elektrostatischen Wechselwirkungen untersuchen wir rheologische Eigenschaften der Doppelschicht und vergleichen mit neutralen Systemen. Wir studieren den kontinuierlichen Übergang der Systemeigenschaften von neutralen zu stark geladenen Bürsten durch Variation der Bjerrumlänge und der Ladungsdichte.

Transport properties of YBa 2 Cu 3 O 7, PrBa 2 Cu 3 O 7-superlattices

The understanding of the coupling between superconducting YBa2Cu3O7 (YBCO) layers decoupled by non superconducting PrBa2Cu3O7 (PBCO) layers in c-axis oriented superlattices was the aim of this thesis. For this purpose two conceptually different kind of transport experiments have been performed. rnrnIn the first type of transport experiments the current is flowing parallel to the layers. Here the coupling is probed indirectly using magnetic vortex lines, which are penetrating the superlattice. Movement of the vortex segments in neighbouring YBCO layers is more or less coherent depending on the thickness of both the superconducting and non superconducting layers. This in-plane transport was measured either by sending an external current through bridges patterned in the superlattice or by an induced internal current. rnThe vortex-creep activation energy U was determined by analysis of the in-plane resistive transition in an external magnetic field B oriented along the c-axis. The activation energies for two series of superlattices were investigated. In one series the thickness of the YBCO layers was constant (nY=4 unit cells) and the number of the PBCO unit cells was varied, while in the other the number of PBCO layers was constant (nP=4) and nY varied. The correlation length of the vortex system was determined to be 80 nm along the c-axis direction. It was found that even a single PBCO unit cell in a superlattice effectively cuts the flux lines into shorter weakly coupled segments, and the coupling of the vortex systems in neighbouring layers is negligible already for a thickness of four unit cells of the PBCO layers. A characteristic variation of the activation energy for the two series of superlattices was found, where U0 is proportional to the YBCO thickness. A change in the variation of U0 with the current I in the specimen was observed, which can be explained in terms of a crossover in the vortex creep process, generated by the transport current. At low I values the dislocations mediated (plastic) vortex creep leads to thermally assisted flux-flow behaviour, whereas at high current the dc transport measurements are dominated by elastic (collective) creep.rnThe analysis of standard dc magnetization relaxation data obtained for a series superlattices revealed the occurrence of a crossover from elastic (collective) vortex creep at low temperature to plastic vortex creep at high T. The crossover is generated by the T dependent macroscopic currents induced in the sample. The existence of this creep crossover suggests that, compared with the well known Maley technique, the use of the normalized vortex creep activation energy is a better solution for the determination of vortex creep parameters.rnrnThe second type of transport experiments was to measure directly a possible Josephson coupling between superconducting CuO2 double planes in the superlattices by investigation of the transport properties perpendicular to the superconducting planes. Here three different experiments have been performed. The first one was to pattern mesa structures photolithographically as in previous works. The second used three-dimensional nanostructures cut by a focused ion beam. For the these two experiments insufficient patterning capabilities prevented an observation of the Josephson effect in the current voltage curves. rnA third experiment used a-axis and (110) oriented YBCO films, where in-plane patterning can in principle be sufficient to measure transport perpendicular to the superconducting planes. Therefore the deposition of films with this unusual growth orientation was optimized and investigated. The structural and microstructural evolution of c-axis to a-axis orientation was monitored using x-ray diffraction, scanning electron microscopy and magnetization measurements. Films with full a-axis alignment parallel to the substrate normal could be achieved on (100)SrTiO3. Due to the symmetry of the substrate the c-axis direction in-plane is twofold. Transferring the deposition conditions to films grown on (110)SrTiO3 allowed the growth of (110) oriented YBCO films with a unique in-plane c-axis orientation. While these films were of high quality by crystallographic and macroscopic visual inspection, electron microscopy revealed a coherent crack pattern on a nanoscale. Therefore the actual current path in the sample was not determined by the macroscopic patterning which prohibited investigations of the in-plane anisotropy in this case.rn

Numerical studies of colloidal systems

The interplay of hydrodynamic and electrostatic forces is of great importance for the understanding of colloidal dispersions. Theoretical descriptions are often based on the so called standard electrokinetic model. This Mean Field approach combines the Stokes equation for the hydrodynamic flow field, the Poisson equation for electrostatics and a continuity equation describing the evolution of the ion concentration fields. In the first part of this thesis a new lattice method is presented in order to efficiently solve the set of non-linear equations for a charge-stabilized colloidal dispersion in the presence of an external electric field. Within this framework, the research is mainly focused on the calculation of the electrophoretic mobility. Since this transport coefficient is independent of the electric field only for small driving, the algorithm is based upon a linearization of the governing equations. The zeroth order is the well known Poisson-Boltzmann theory and the first order is a coupled set of linear equations. Furthermore, this set of equations is divided into several subproblems. A specialized solver for each subproblem is developed, and various tests and applications are discussed for every particular method. Finally, all solvers are combined in an iterative procedure and applied to several interesting questions, for example, the effect of the screening mechanism on the electrophoretic mobility or the charge dependence of the field-induced dipole moment and ion clouds surrounding a weakly charged sphere. In the second part a quantitative data analysis method is developed for a new experimental approach, known as "Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy" (TIR-FCCS). The TIR-FCCS setup is an optical method using fluorescent colloidal particles to analyze the flow field close to a solid-fluid interface. The interpretation of the experimental results requires a theoretical model, which is usually the solution of a convection-diffusion equation. Since an analytic solution is not available due to the form of the flow field and the boundary conditions, an alternative numerical approach is presented. It is based on stochastic methods, i. e. a combination of a Brownian Dynamics algorithm and Monte Carlo techniques. Finally, experimental measurements for a hydrophilic surface are analyzed using this new numerical approach.

Donor–acceptor systems in the quest for organic semiconductors

Die letzten Jahrzehnte brachten eine Vielzahl neuer organischen Halbleiter hervor, welche erfolgreich als aktive Materialien in Bauteilen eingesetzt wurden, wie zum Beispiel Feldeffekttransistoren (FET), organische Leuchtdioden (OLED), organischen Photovoltaikzellen (OPV) und Sensoren. Einige dieser Materialien haben, obwohl sich die Technolgie noch in der „Pubertät“ befindet, die minimalen Anforderungen für eine kommerzielle Anwendung erreicht, wobei jedoch vieles noch zu entdecken, erklären und verstehen bleibt. Diese Arbeit beschreibt das Design, die Synthese und Charakterisierung neuartiger halbleitender Polymere mit speziell eingestellten optoelektronischen Eigenschaften, welche effiziente ambipolare oder n-Leitung in OFET’s und OPV’s zeigen. Das Hauptziel wurde dadurch erreicht, dass sowohl die vorteilhaften Eigenschaften des planaren, elektronenarmen heterozyklischen Bausteines Thiadiazolo[3,4-g]quinoxalin als auch von Ethinbrücken, welche den Donor (D) und den Akzeptor (A) in einem D-A-Copolymer verbinden, durch systematische Optimierung ausgenutzt wurden. Neben synthetischen Herausforderungen werden in dieser Arbeit auch detailiiete Untersuchungen der optoelektronischen Eigenschaften der hergestellten konjugierten Polymere und Modellverbindungen dargelegt. Darüber hinaus beschreibt diese Arbeit erstmals ein Beispiel für ein Polymer, welches Dreifachbindungen im Polymerrückgrat enthält, und nahezu eine ausgeglichene ambipolare Ladungsträgerleitung in OFET’s zeigt. Zusätzlich werden gemischt-valente Phenothiazine, verbrückt mittels elektronenarmen pi-Brücken wie etwa Benzo[c][2,1,3]thiadiazol, und deren Elektronentransferprozesse, im Rahmen der Marcus-Hush-Theorie, untersucht.

Charge-transport simulations in organic semiconductors

In this thesis we have extended the methods for microscopic charge-transport simulations for organic semiconductors. In these materials the weak intermolecular interactions lead to spatially localized charge carriers, and the charge transport occurs as an activated hopping process between diabatic states. In addition to weak electronic couplings between these states, different electrostatic environments in the organic material lead to a broadening of the density of states for the charge energies which limits carrier mobilities.rnThe contributions to the method development includern(i) the derivation of a bimolecular charge-transfer rate,rn(ii) the efficient evaluation of intermolecular (outer-sphere) reorganization energies,rn(iii) the investigation of effects of conformational disorder on intramolecular reorganization energies or internal site energiesrnand (iv) the inclusion of self-consistent polarization interactions for calculation of charge energies.These methods were applied to study charge transport in amorphous phases of small molecules used in the emission layer of organic light emitting diodes (OLED).rnWhen bulky substituents are attached to an aromatic core in order to adjust energy levels or prevent crystallization, a small amount of delocalization of the frontier orbital to the substituents can increase electronic couplings between neighboring molecules. This leads to improved charge-transfer rates and, hence, larger charge-mobility. We therefore suggest using the mesomeric effect (as opposed to the inductive effect) when attaching substituents to aromatic cores, which is necessary for example in deep blue OLEDs, where the energy levels of a host molecule have to be adjusted to those of the emitter.rnFurthermore, the energy landscape for charges in an amorphous phase cannot be predicted by mesoscopic models because they approximate the realistic morphology by a lattice and represent molecular charge distributions in a multipole expansion. The microscopic approach shows that a polarization-induced stabilization of a molecule in its charged and neutral states can lead to large shifts, broadening, and traps in the distribution of charge energies. These results are especially important for multi-component systems (the emission layer of an OLED or the donor-acceptor interface of an organic solar cell), if the change in polarizability upon charging (or excitation in case of energy transport) is different for the components. Thus, the polarizability change upon charging or excitation should be added to the set of molecular parameters essential for understanding charge and energy transport in organic semiconductors.rnWe also studied charge transport in self-assembled systems, where intermolecular packing motives induced by side chains can increase electronic couplings between molecules. This leads to larger charge mobility, which is essential to improve devices such as organic field effect transistors, where low carrier mobilities limit the switching frequency.rnHowever, it is not sufficient to match the average local molecular order induced by the sidernchains (such as the pitch angle between consecutive molecules in a discotic mesophase) with maxima of the electronic couplings.rnIt is also important to make the corresponding distributions as narrow as possible compared to the window determined by the closest minima of thernelectronic couplings. This is especially important in one-dimensional systems, where charge transport is limited by the smallest electronic couplings.rnThe immediate implication for compound design is that the side chains should assist the self-assemblingrnprocess not only via soft entropic interactions, but also via stronger specific interactions, such as hydrogen bonding.rnrnrnrn

Transport of nanoparticles into polymersomes : a minimal model system of particles passage through biological membranes

This thesis focuses on the design and characterization of a novel, artificial minimal model membrane system with chosen physical parameters to mimic a nanoparticle uptake process driven exclusively by adhesion and softness of the bilayer. The realization is based on polymersomes composed of poly(dimethylsiloxane)-b-poly(2-methyloxazoline) (PMDS-b-PMOXA) and nanoscopic colloidal particles (polystyrene, silica), and the utilization of powerful characterization techniques. rnPDMS-b-PMOXA polymersomes with a radius, Rh ~100 nm, a size polydispersity, PD = 1.1 and a membrane thickness, h = 16 nm, were prepared using the film rehydratation method. Due to the suitable mechanical properties (Young’s modulus of ~17 MPa and a bending modulus of ~7⋅10-8 J) along with the long-term stability and the modifiability, these kind of polymersomes can be used as model membranes to study physical and physicochemical aspects of transmembrane transport of nanoparticles. A combination of photon (PCS) and fluorescence (FCS) correlation spectroscopies optimizes species selectivity, necessary for a unique internalization study encompassing two main efforts. rnFor the proof of concepts, the first effort focused on the interaction of nanoparticles (Rh NP SiO2 = 14 nm, Rh NP PS = 16 nm; cNP = 0.1 gL-1) and polymersomes (Rh P = 112 nm; cP = 0.045 gL-1) with fixed size and concentration. Identification of a modified form factor of the polymersome entities, selectively seen in the PCS experiment, enabled a precise monitor and quantitative description of the incorporation process. Combining PCS and FCS led to the estimation of the incorporated particles per polymersome (about 8 in the examined system) and the development of an appropriate methodology for the kinetics and dynamics of the internalization process. rnThe second effort aimed at the establishment of the necessary phenomenology to facilitate comparison with theories. The size and concentration of the nanoparticles were chosen as the most important system variables (Rh NP = 14 - 57 nm; cNP = 0.05 - 0.2 gL-1). It was revealed that the incorporation process could be controlled to a significant extent by changing the nanoparticles size and concentration. Average number of 7 up to 11 NPs with Rh NP = 14 nm and 3 up to 6 NPs with Rh NP = 25 nm can be internalized into the present polymersomes by changing initial nanoparticles concentration in the range 0.1- 0.2 gL-1. Rapid internalization of the particles by polymersomes is observed only above a critical threshold particles concentration, dependent on the nanoparticle size. rnWith regard possible pathways for the particle uptake, cryogenic transmission electron microscopy (cryo-TEM) has revealed two different incorporation mechanisms depending on the size of the involved nanoparticles: cooperative incorporation of nanoparticles groups or single nanoparticles incorporation. Conditions for nanoparticle uptake and controlled filling of polymersomes were presented. rnIn the framework of this thesis, the experimental observation of transmembrane transport of spherical PS and SiO2 NPs into polymersomes via an internalization process was reported and examined quantitatively for the first time. rnIn a summary the work performed in frames of this thesis might have significant impact on cell model systems’ development and thus improved understanding of transmembrane transport processes. The present experimental findings help create the missing phenomenology necessary for a detailed understanding of a phenomenon with great relevance in transmembrane transport. The fact that transmembrane transport of nanoparticles can be performed by artificial model system without any additional stimuli has a fundamental impact on the understanding, not only of the nanoparticle invagination process but also of the interaction of nanoparticles with biological as well as polymeric membranes. rn

Design of a new harrow type wool transport mechanism to reduce fibre entanglement

The wool is entangled at several stages of its processing. In the conventional scouring machines, the prongs or the rakes agitate the wool and lead the fiber entanglement. Several scouring systems have been commercialized in order to reduce the fiber entanglement. In spite of the existing technologies, the conventional scouring machines are widely used in wool processing. In this thesis, a new approach for the harrow type wool transport mechanism has been introduced. The proposed mechanism has been designed based on the motion of the conventional harrow type wool transport mechanism by exploiting new synthesis concepts. The developed structure has been synthesized based on the Hrones and Nelson's "Atlas of four bar linkages". The four bar linkage has been applied for the desired trajectory of the developed wool transport mechanism. The prongs of the developed mechanism immerse the wool into the scouring liquor and gently propel forward toward the end of the machine with approximately straight line motion in a certain length instead of circular or elliptical motion typical of the conventional machines.

Charge transport in organic photovoltaic cells

Die vorliegende Dissertation dient dazu, das Verständnis des Ladungstransportes in organischen Solarzellen zu vertiefen. Mit Hilfe von Computersimulationen wird die Bewegung von Ladungsträgern in organischen Materialien rekonstruiert, und zwar ausgehend von den quantenmechanischen Prozessen auf mikroskopischer Ebene bis hin zur makroskopischen Skala, wo Ladungsträgermobilitäten quantifizierbar werden. Auf Grundlage dieses skalenübergreifenden Ansatzes werden Beziehungen zwischen der chemischen Struktur organischer Moleküle und der makroskopischen Mobilität hergestellt (Struktur-Eigenschafts-Beziehungen), die zu der Optimierung photovoltaischer Wirkungsgrade beitragen. Das Simulationsmodell beinhaltet folgende drei Schlüsselkomponenten. Erstens eine Morphologie, d. h. ein atomistisch aufgelöstes Modell der molekularen Anordnung in dem untersuchten Material. Zweitens ein Hüpfmodell des Ladungstransportes, das Ladungswanderung als eine Abfolge von Ladungstransferreaktionen zwischen einzelnen Molekülen beschreibt. Drittens ein nichtadiabatisches Modell des Ladungstransfers, das Übergangsraten durch drei Parameter ausdrückt: Reorganisationsenergien, Lageenergien und Transferintegrale. Die Ladungstransport-Simulationen richten sich auf die Materialklasse der dicyanovinyl-substituierten Oligothiophene und umfassen Morphologien von Einkristallen, Dünnschichten sowie amorphen/smektischen Mesophasen. Ein allgemeiner Befund ist, dass die molekulare Architektur, bestehend aus einer Akzeptor-Donor-Akzeptor-Sequenz und einem flexiblen Oligomergerüst, eine erhebliche Variation molekularer Dipolmomente und damit der Lageenergien bewirkt. Diese energetische Unordnung ist ungewöhnlich hoch in den Kristallen und umso höher in den Mesophasen. Für die Einkristalle wird beobachtet, dass Kristallstrukturen mit ausgeprägter π-Stapelung und entsprechend großer Transferintegrale zu verhältnismäßig niedrigen Mobilitäten führen. Dieses Verhalten wird zurückgeführt auf die Ausbildung bevorzugter Transportrichtungen, die anfällig für energetische Störungen sind. Für die Dünnschichten bestätigt sich diese Argumentation und liefert ein mikroskopisches Verständnis für experimentelle Mobilitäten. In der Tat korrelieren die Simulationsergebnisse sowohl mit gemessenen Mobilitäten als auch mit photovoltaischen Wirkungsgraden. Für die amorphen/smektischen Systeme steigt die energetische Unordnung mit der Oligomerlänge, sie führt aber auch zu einer unerwarteten Mobilitätsabnahme in dem stärker geordneten smektischen Zustand. Als Ursache dafür erweist sich, dass die smektische Schichtung der räumlichen Korrelation der energetischen Unordnung entgegensteht.

Transport of Hydraulic Fracturing Water and Wastes in the Susquehanna River Basin, Pennsylvania

The development of the Marcellus Shale gas play in Pennsylvania and the northeastern United States has resulted in significant amounts of water and wastes transported by truck over roadways. This study used geographic information systems (GIS) to quantify truck travel distances via both the preferred routes (minimum distance while also favoring higher-order roads) as well as, where available, the likely actual distances for freshwater and waste transport between pertinent locations (e. g., gas wells, treatment facilities, freshwater sources). Results show that truck travel distances in the Susquehanna River Basin are greater than those used in prior life-cycle assessments of tight shale gas. When compared to likely actual transport distances, if policies were instituted to constrain truck travel to the closest destination and higher-order roads, transport mileage reductions of 40-80% could be realized. Using reasonable assumptions of current practices, greenhouse gas (GHG) emissions associated with water and waste hauling were calculated to be 70-157 MT CO2 eq per gas well. Furthermore, empty so-called backhaul trips, such as to freshwater withdrawal sites or returning from deep well injection sites, were found to increase emissions by an additional 30%, underscoring the importance of including return trips in the analysis. The results should inform future life-cycle assessments of tight shale gases in managed watersheds and help local and regional governments plan for impacts of transportation on local infrastructure. (C) 2013 American Society of Civil Engineers.

Diffusion-driven transport in clayrock formations

Clay mineral-rich sedimentary formations are currently under investigation to evaluate their potential use as host formations for installation of deep underground disposal facilities for radioactive waste (e.g. Boom Clay (BE), Opalinus Clay (CH), Callovo-Oxfordian argillite (FR)). The ultimate safety of the corresponding repository concepts depends largely on the capacity of the host formation to limit the flux towards the biosphere of radionuclides (RN) contained in the waste to acceptably low levels. Data for diffusion-driven transfer in these formations shows extreme differences in the measured or modelled behaviour for various radionuclides, e. g. between halogen RN (Cl-36, I-129) and actinides (U-238,U-235, Np-237, Th-232, etc.), which result from major differences between RN of the effects on transport of two phenomena: diffusion and sorption. This paper describes recent research aimed at improving understanding of these two phenomena, focusing on the results of studies carried out during the EC Funmig IP on clayrocks from the above three formations and from the Boda formation (HU). Project results regarding phenomena governing water, cation and anion distribution and mobility in the pore volumes influenced by the negatively-charged surfaces of clay minerals show a convergence of the modelling results for behaviour at the molecular scale and descriptions based on electrical double layer models. Transport models exist which couple ion distribution relative to the clay-solution interface and differentiated diffusive characteristics. These codes are able to reproduce the main trends in behaviour observed experimentally, e.g. D-e(anion) < D-e(HTO) < D-e(cation) and D-e(anion) variations as a function of ionic strength and material density. These trends are also well-explained by models of transport through ideal porous matrices made up of a charged surface material. Experimental validation of these models is good as regards monovalent alkaline cations, in progress for divalent electrostatically-interacting cations (e.g. Sr2+) and still relatively poor for 'strongly sorbing', high K-d cations. Funmig results have clarified understanding of how clayrock mineral composition, and the corresponding organisation of mineral grain assemblages and their associated porosity, can affect mobile solute (anions, HTO) diffusion at different scales (mm to geological formation). In particular, advances made in the capacity to map clayrock mineral grain-porosity organisation at high resolution provide additional elements for understanding diffusion anisotropy and for relating diffusion characteristics measured at different scales. On the other hand, the results of studies focusing on evaluating the potential effects of heterogeneity on mobile species diffusion at the formation scale tend to show that there is a minimal effect when compared to a homogeneous property model. Finally, the results of a natural tracer-based study carried out on the Opalinus Clay formation increase confidence in the use of diffusion parameters measured on laboratory scale samples for predicting diffusion over geological time-space scales. Much effort was placed on improving understanding of coupled sorption-diffusion phenomena for sorbing cations in clayrocks. Results regarding sorption equilibrium in dispersed and compacted materials for weakly to moderately sorbing cations (Sr2+, Cs+, Co2+) tend to show that the same sorption model probably holds in both systems. It was not possible to demonstrate this for highly sorbing elements such as Eu(III) because of the extremely long times needed to reach equilibrium conditions, but there does not seem to be any clear reason why such elements should not have similar behaviour. Diffusion experiments carried out with Sr2+, Cs+ and Eu(III) on all of the clayrocks gave mixed results and tend to show that coupled diffusion-sorption migration is much more complex than expected, leading generally to greater mobility than that predicted by coupling a batch-determined K-d and Ficks law based on the diffusion behaviour of HTO. If the K-d measured on equivalent dispersed systems holds as was shown to be the case for Sr, Cs (and probably Co) for Opalinus Clay, these results indicate that these cations have a D-e value higher than HTO (up to a factor of 10 for Cs+). Results are as yet very limited for very moderate to strongly sorbing species (e.g. Co(II), Eu(III), Cu(II)) because of their very slow transfer characteristics. (C) 2011 Elsevier Ltd. All rights reserved.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.

STOCHASTIC CHARGE TRANSPORT IN MULTI-ISLAND SINGLE-ELECTRON TUNNELING DEVICES

The physics of the operation of singe-electron tunneling devices (SEDs) and singe-electron tunneling transistors (SETs), especially of those with multiple nanometer-sized islands, has remained poorly understood in spite of some intensive experimental and theoretical research. This computational study examines the current-voltage (IV) characteristics of multi-island single-electron devices using a newly developed multi-island transport simulator (MITS) that is based on semi-classical tunneling theory and kinetic Monte Carlo simulation. The dependence of device characteristics on physical device parameters is explored, and the physical mechanisms that lead to the Coulomb blockade (CB) and Coulomb staircase (CS) characteristics are proposed. Simulations using MITS demonstrate that the overall IV characteristics in a device with a random distribution of islands are a result of a complex interplay among those factors that affect the tunneling rates that are fixed a priori (e.g. island sizes, island separations, temperature, gate bias, etc.), and the evolving charge state of the system, which changes as the source-drain bias (VSD) is changed. With increasing VSD, a multi-island device has to overcome multiple discrete energy barriers (up-steps) before it reaches the threshold voltage (Vth). Beyond Vth, current flow is rate-limited by slow junctions, which leads to the CS structures in the IV characteristic. Each step in the CS is characterized by a unique distribution of island charges with an associated distribution of tunneling probabilities. MITS simulation studies done on one-dimensional (1D) disordered chains show that longer chains are better suited for switching applications as Vth increases with increasing chain length. They are also able to retain CS structures at higher temperatures better than shorter chains. In sufficiently disordered 2D systems, we demonstrate that there may exist a dominant conducting path (DCP) for conduction, which makes the 2D device behave as a quasi-1D device. The existence of a DCP is sensitive to the device structure, but is robust with respect to changes in temperature, gate bias, and VSD. A side gate in 1D and 2D systems can effectively control Vth. We argue that devices with smaller island sizes and narrower junctions may be better suited for practical applications, especially at room temperature.

The effect of demand, tank parameters, and pumping stations on energy use in municipal drinking water distribution systems

Two of the indicators of the UN Millennium Development Goals ensuring environmental sustainability are energy use and per capita carbon dioxide emissions. The increasing urbanization and increasing world population may require increased energy use in order to transport enough safe drinking water to communities. In addition, the increase in water use would result in increased energy consumption, thereby resulting in increased green-house gas emissions that promote global climate change. The study of multiple Municipal Drinking Water Distribution Systems (MDWDSs) that relates various MDWDS aspects--system components and properties--to energy use is strongly desirable. The understanding of the relationship between system aspects and energy use aids in energy-efficient design. In this study, components of a MDWDS, and/or the characteristics associated with the component are termed as MDWDS aspects (hereafter--system aspects). There are many aspects of MDWDSs that affect the energy usage. Three system aspects (1) system-wide water demand, (2) storage tank parameters, and (3) pumping stations were analyzed in this study. The study involved seven MDWDSs to understand the relationship between the above-mentioned system aspects in relation with energy use. A MDWDSs model, EPANET 2.0, was utilized to analyze the seven systems. Six of the systems were real and one was a hypothetical system. The study presented here is unique in its statistical approach using seven municipal water distribution systems. The first system aspect studied was system-wide water demand. The analysis involved analyzing seven systems for the variation of water demand and its impact on energy use. To quantify the effects of water use reduction on energy use in a municipal water distribution system, the seven systems were modeled and the energy usage quantified for various amounts of water conservation. It was found that the effect of water conservation on energy use was linear for all seven systems and that all the average values of all the systems' energy use plotted on the same line with a high R 2 value. From this relationship, it can be ascertained that a 20% reduction in water demand results in approximately a 13% savings in energy use for all seven systems analyzed. This figure might hold true for many similar systems that are dominated by pumping and not gravity driven. The second system aspect analyzed was storage tank(s) parameters. Various tank parameters: (1) tank maximum water levels, (2) tank elevation, and (3) tank diameter were considered in this part of the study. MDWDSs use a significant amount of electrical energy for the pumping of water from low elevations (usually a source) to higher ones (usually storage tanks). The use of electrical energy has an effect on pollution emissions and, therefore, potential global climate change as well. Various values of these tank parameters were modeled on seven MDWDSs of various sizes using a network solver and the energy usage recorded. It was found that when averaged over all seven analyzed systems (1) the reduction of maximum tank water level by 50% results in a 2% energy reduction, (2) energy use for a change in tank elevation is system specific, and (2) a reduction of tank diameter of 50% results in approximately a 7% energy savings. The third system aspect analyzed in this study was pumping station parameters. A pumping station consists of one or more pumps. The seven systems were analyzed to understand the effect of the variation of pump horsepower and the number of booster stations on energy use. It was found that adding booster stations could save energy depending upon the system characteristics. For systems with flat topography, a single main pumping station was found to use less energy. In systems with a higher-elevation neighborhood, however, one or more booster pumps with a reduced main pumping station capacity used less energy. The energy savings for the seven systems was dependent on the number of boosters and ranged from 5% to 66% for the analyzed five systems with higher elevation neighborhoods (S3, S4, S5, S6, and S7). No energy savings was realized for the remaining two flat topography systems, S1, and S2. The present study analyzed and established the relationship between various system aspects and energy use in seven MDWDSs. This aids in estimating the amount of energy savings in MDWDSs. This energy savings would ultimately help reduce Greenhouse gases (GHGs) emissions including per capita CO 2 emissions thereby potentially lowering the global climate change effect. This will in turn contribute to meeting the MDG of ensuring environmental sustainability.

On Transport Service Selection in Intermodal Rail/Road Distribution Networks

Intermodal rail/road freight transport constitutes an alternative to long-haul road transport for the distribution of large volumes of goods. The paper introduces the intermodal transportation problem for the tactical planning of mode and service selection. In rail mode, shippers either book train capacity on a per-unit basis or charter block trains completely. Road mode is used for short-distance haulage to intermodal terminals and for direct shipments to customers. We analyze the competition of road and intermodal transportation with regard to freight consolidation and service cost on a model basis. The approach is applied to a distribution system of an industrial company serving customers in eastern Europe. The case study investigates the impact of transport cost and consolidation on the optimal modal split.

Performance Availability and Anticipatory Change Planning of Intralogistics Systems: A Simulation-Based Approach

In recent years, the ability to respond to real time changes in operations and reconfigurability in equipment are likely to become essential characteristics for next generation intralogistics systems as well as the level of automation, cost effectiveness and maximum throughput. In order to cope with turbulences and the increasing level of dynamic conditions, future intralogistics systems have to feature short reaction times, high flexibility in processes and the ability to adapt to frequent changes. The increasing autonomy and complexity in processes of today’s intralogistics systems requires new and innovative management approaches, which allow a fast response to (un)anticipated events and adaptation to changing environment in order to reduce the negative consequences of these events. The ability of a system to respond effectively a disruption depends more on the decisions taken before the event than those taken during or after. In this context, anticipatory change planning can be a usable approach for managers to make contingency plans for intralogistics systems to deal with the rapidly changing marketplace. This paper proposes a simulation-based decision making framework for the anticipatory change planning of intralogistics systems. This approach includes the quantitative assessments based on the simulation in defined scenarios as well as the analysis of performance availability that combines the flexibility corridors of different performance dimensions. The implementation of the approach is illustrated on a new intralogistics technology called the Cellular Transport System.