Advanced multivariate analysis to assess remediation of hydrocarbons in soils

Accurate monitoring of degradation levels in soils is essential in order to understand and achieve complete degradation of petroleum hydrocarbons in contaminated soils. We aimed to develop the use of multivariate methods for the monitoring of biodegradation of diesel in soils and to determine if diesel contaminated soils could be remediated to a chemical composition similar to that of an uncontaminated soil. An incubation experiment was set up with three contrasting soil types. Each soil was exposed to diesel at varying stages of degradation and then analysed for key hydrocarbons throughout 161 days of incubation. Hydrocarbon distributions were analysed by Principal Coordinate Analysis and similar samples grouped by cluster analysis. Variation and differences between samples were determined using permutational multivariate analysis of variance. It was found that all soils followed trajectories approaching the chemical composition of the unpolluted soil. Some contaminated soils were no longer significantly different to that of uncontaminated soil after 161 days of incubation. The use of cluster analysis allows the assignment of a percentage chemical similarity of a diesel contaminated soil to an uncontaminated soil sample. This will aid in the monitoring of hydrocarbon contaminated sites and the establishment of potential endpoints for successful remediation.

Normal Modes Using Maple (Water)

The normal modes of water are obtained using Maple

Entrainment Effects in Turbulent Boundary Layers

Esta tesis estudia el comportamiento de la región exterior de una capa límite turbulenta sin gradientes de presiones. Se ponen a prueba dos teorías relativamente bien establecidas. La teoría de semejanza para la pared supone que en el caso de haber una pared rugosa, el fluido sólo percibe el cambio en la fricción superficial que causa, y otros efectos secundarios quedarán confinados a una zona pegada a la pared. El consenso actual es que dicha teoría es aproximadamente cierta. En el extremo exterior de la capa límite existe una región producida por la interacción entre las estructuras turbulentas y el flujo irrotacional de la corriente libre llamada interfaz turbulenta/no turbulenta. La mayoría de los resultados al respecto sugieren la presencia de fuerzas de cortadura ligeramente más intensa, lo que la hace distinta al resto del flujo turbulento. Las propiedades de esa región probablemente cambien si la velocidad de crecimiento de la capa límite aumenta, algo que puede conseguirse aumentando la fricción en la pared. La rugosidad y la ingestión de masa están entonces relacionadas, y el comportamiento local de la interfaz turbulenta/no turbulenta puede explicar el motivo por el que las capas límite sobre paredes rugosas no se comportan como en el caso de tener paredes lisas precisamente en la zona exterior. Para estudiar las capas límite a números de Reynolds lo suficientemente elevados, se ha desarrollado un nuevo código de alta resolución para la simulación numérica directa de capas límite turbulentas sin gradiente de presión. Dicho código es capaz de simular capas límite en un intervalo de números de Reynolds entre ReT = 100 — 2000 manteniendo una buena escalabilidad hasta los dos millones de hilos en superordenadores de tipo Blue Gene/Q. Se ha guardado especial atención a la generación de condiciones de contorno a la entrada correctas. Los resultados obtenidos están en concordancia con los resultados previos, tanto en el caso de simulaciones como de experimentos. La interfaz turbulenta/no turbulenta de una capa límite se ha analizado usando un valor umbral del módulo de la vorticidad. Dicho umbral se considera un parámetro para analizar cada superficie obtenida de un contorno del módulo de la vorticidad. Se han encontrado dos regímenes distintos en función del umbral escogido con propiedades opuestas, separados por una transición topológica gradual. Las características geométricas de la zona escalan con o99 cuando u^/isdgg es la unidad de vorticidad. Las propiedades del íluido relativas a la posición del contorno de vorticidad han sido analizados para una serie de umbrales utilizando el campo de distancias esféricas, que puede obtenerse con independencia de la complejidad de la superficie de referencia. Las propiedades del fluido a una distancia dada del inerfaz también dependen del umbral de vorticidad, pero tienen características parecidas con independencia del número de Reynolds. La interacción entre la turbulencia y el flujo no turbulento se restringe a una zona muy fina con un espesor del orden de la escala de Kolmogorov local. Hacia el interior del flujo turbulento las propiedades son indistinguibles del resto de la capa límite. Se ha simulado una capa límite sin gradiente de presiones con una fuerza volumétrica cerca de la pared. La el forzado ha sido diseñado para aumentar la fricción en la pared sin introducir ningún efecto geométrico obvio. La simulación consta de dos dominios, un primer dominio más pequeño y a baja resolución que se encarga de generar condiciones de contorno correctas, y un segundo dominio mayor y a alta resolución donde se aplica el forzado. El estudio de los perfiles y los coeficientes de autocorrelación sugieren que los dos casos, el liso y el forzado, no colapsan más allá de la capa logarítmica por la complejidad geométrica de la zona intermitente, y por el hecho que la distancia a la pared no es una longitud característica. Los efectos causados por la geometría de la zona intermitente pueden evitarse utilizando el interfaz como referencia, y la distancia esférica para el análisis de sus propiedades. Las propiedades condicionadas del flujo escalan con 5QQ y u/uT, las dos únicas escalas contenidas en el modelo de semejanza de pared de Townsend, consistente con estos resultados. ABSTRACT This thesis studies the characteristics of the outer region of zero-pressure-gradient turbulent boundary layers at moderate Reynolds numbers. Two relatively established theories are put to test. The wall similarity theory states that with the presence of roughness, turbulent motion is mostly affected by the additional drag caused by the roughness, and that other secondary effects are restricted to a region very close to the wall. The consensus is that this theory is valid, but only as a first approximation. At the edge of the boundary layer there is a thin layer caused by the interaction between the turbulent eddies and the irroational fluid of the free stream, called turbulent/non-turbulent interface. The bulk of results about this layer suggest the presence of some localized shear, with properties that make it distinguishable from the rest of the turbulent flow. The properties of the interface are likely to change if the rate of spread of the turbulent boundary layer is amplified, an effect that is usually achieved by increasing the drag. Roughness and entrainment are therefore linked, and the local features of the turbulent/non-turbulent interface may explain the reason why rough-wall boundary layers deviate from the wall similarity theory precisely far from the wall. To study boundary layers at a higher Reynolds number, a new high-resolution code for the direct numerical simulation of a zero pressure gradient turbulent boundary layers over a flat plate has been developed. This code is able to simulate a wide range of Reynolds numbers from ReT =100 to 2000 while showing a linear weak scaling up to around two million threads in the BG/Q architecture. Special attention has been paid to the generation of proper inflow boundary conditions. The results are in good agreement with existing numerical and experimental data sets. The turbulent/non-turbulent interface of a boundary layer is analyzed by thresholding the vorticity magnitude field. The value of the threshold is considered a parameter in the analysis of the surfaces obtained from isocontours of the vorticity magnitude. Two different regimes for the surface can be distinguished depending on the threshold, with a gradual topological transition across which its geometrical properties change significantly. The width of the transition scales well with oQg when u^/udgg is used as a unit of vorticity. The properties of the flow relative to the position of the vorticity magnitude isocontour are analyzed within the same range of thresholds, using the ball distance field, which can be obtained regardless of the size of the domain and complexity of the interface. The properties of the flow at a given distance to the interface also depend on the threshold, but they are similar regardless of the Reynolds number. The interaction between the turbulent and the non-turbulent flow occurs in a thin layer with a thickness that scales with the Kolmogorov length. Deeper into the turbulent side, the properties are undistinguishable from the rest of the turbulent flow. A zero-pressure-gradient turbulent boundary layer with a volumetric near-wall forcing has been simulated. The forcing has been designed to increase the wall friction without introducing any obvious geometrical effect. The actual simulation is split in two domains, a smaller one in charge of the generation of correct inflow boundary conditions, and a second and larger one where the forcing is applied. The study of the one-point and twopoint statistics suggest that the forced and the smooth cases do not collapse beyond the logarithmic layer may be caused by the geometrical complexity of the intermittent region, and by the fact that the scaling with the wall-normal coordinate is no longer present. The geometrical effects can be avoided using the turbulent/non-turbulent interface as a reference frame, and the minimum distance respect to it. The conditional analysis of the vorticity field with the alternative reference frame recovers the scaling with 5QQ and v¡uT already present in the logarithmic layer, the only two length-scales allowed if Townsend’s wall similarity hypothesis is valid.

Critical Velocity obtained using Simplified Models of the Railway Track: Viability and Applicability

Increased demands on the capacity of the railway network gave rise to new issues related to the dynamic response of railway tracks subjected to moving vehicles. Thus, it becomes important to evaluate the applicability of traditionally used simplified models which have a closed form solution. Regarding simplified models, transversal vibrations of a beam on a visco-elastic foundation subjected to a moving load are considered. Governing equations are obtained by Hamilton’s principle. Shear distortion, rotary inertia and effect of axial force are accounted for. The load is introduced as a time varying force moving at a constant velocity. Transversal vibrations induced by the load are solved by the normal-mode analysis. Reflected waves at the extremities of the full beam are avoided by introduction of semi-infinite elements. Firstly, the critical velocity obtained from this model is compared with results of an undamped Euler- Bernoulli formulation with zero axial force. Secondly, a finite element model in ABAQUS is examined. The new contribution lies in the introduction of semi- infinite elements and in the first step to a systematic comparison, which have not been published so fa

Crystal structure of the zinc-, cobalt-, and iron-containing adenylate kinase from Desulfovibrio gigas: a novel metal-containing adenylate kinase from Gram-negative bacteria

J Biol Inorg Chem (2011) 16:51–61 DOI 10.1007/s00775-010-0700-8

Ant species distribution along a topographic gradient in a "terra-firme" forest reserve in Central Amazonia

In Amazonia, topographical variations in soil and forest structure within "terra-firme" ecosystems are important factors correlated with terrestrial invertebrates' distribution. The objective of this work was to assess the effects of soil clay content and slope on ant species distribution over a 25 km² grid covering the natural topographic continuum. Using three complementary sampling methods (sardine baits, pitfall traps and litter samples extracted in Winkler sacks), 300 subsamples of each method were taken in 30 plots distributed over a wet tropical forest in the Ducke Reserve (Manaus, AM, Brazil). An amount of 26,814 individuals from 11 subfamilies, 54 genera, 85 species and 152 morphospecies was recorded (Pheidole represented 37% of all morphospecies). The genus Eurhopalothrix was registered for the first time for the reserve. Species number was not correlated with slope or clay content, except for the species sampled from litter. However, the Principal Coordinate Analysis indicated that the main pattern of species composition from pitfall and litter samples was related to clay content. Almost half of the species were found only in valleys or only on plateaus, which suggests that most of them are habitat specialists. In Central Amazonia, soil texture is usually correlated with vegetation structure and moisture content, creating different microhabitats, which probably account for the observed differences in ant community structure.

On the path integral formulation and the evaluation of quantum statistical averages

Four problems of physical interest have been solved in this thesis using the path integral formalism. Using the trigonometric expansion method of Burton and de Borde (1955), we found the kernel for two interacting one dimensional oscillators• The result is the same as one would obtain using a normal coordinate transformation, We next introduced the method of Papadopolous (1969), which is a systematic perturbation type method specifically geared to finding the partition function Z, or equivalently, the Helmholtz free energy F, of a system of interacting oscillators. We applied this method to the next three problems considered• First, by summing the perturbation expansion, we found F for a system of N interacting Einstein oscillators^ The result obtained is the same as the usual result obtained by Shukla and Muller (1972) • Next, we found F to 0(Xi)f where A is the usual Tan Hove ordering parameter* The results obtained are the same as those of Shukla and Oowley (1971), who have used a diagrammatic procedure, and did the necessary sums in Fourier space* We performed the work in temperature space• Finally, slightly modifying the method of Papadopolous, we found the finite temperature expressions for the Debyecaller factor in Bravais lattices, to 0(AZ) and u(/K/ j,where K is the scattering vector* The high temperature limit of the expressions obtained here, are in complete agreement with the classical results of Maradudin and Flinn (1963) .

Caractérisation de l’ADNc et de l’ADN génomique codant l’intégrine plaquettaire αIIb‐β3 chez un cheval atteint de Thrombasthénie de Glanzmann

La thrombasthénie de Glanzmann (TG) est une maladie caractérisée par un défaut d’agrégation plaquettaire. C’est une maladie génétique autosomale récessive causée par une anomalie du récepteur plaquettaire pour le fibrinogène. Ce récepteur est une intégrine localisée à la surface plasmatique qui est formée par un complexe composé des sous‐unités αIIb et β3. Nous avons identifié un cheval démontrant les caractéristiques clinicopathologiques de la TG. Des études par cytométrie de flux ont révélé une déficience au niveau de la portion αIIb du récepteur. Ces résultats suggèrent une ou plusieurs mutations au niveau du gène codant pour cette portion αIIb du récepteur. L’objectif de notre étude était de caractériser l’ADNc et l’ADN génomique codant pour les gènes ITGA2B et ITGB3 codant respectivement pour les deux sous‐unités αIIb et β3 chez un cheval atteint de la TG. L’ADNc a été synthétisé par RT‐PCR en utilisant l’ARN total récolté à partir des plaquettes. L’ADN génomique a été extrait à partir des globules blancs. Des amorces spécifiques ont été utilisées pour l’amplification par PCR d’ITGA2B et d’ITGB3. Les séquences d’ADNc et d’ADN génomique de notre patient ont été caractérisées par séquençage et comparées par l’analyse BLAST (GenBank). Une substitution d’une guanine par une cytosine a été mise en évidence au niveau de l’exon 2 d’ITGA2B amenant à la substitution d’une arginine (Arg72) par une proline (Pro72). Ce changement d’acide aminé pourrait résulter en une conformation structurelle anormale qui amènerait à une sous‐unité αIIb inactive. L’analyse de l’ADN génomique a démontré que ce cheval était homozygote pour cette mutation. Le séquençage de l’ADN génomique des parents et de la grand‐mère du patient a démontré que ces individus étaient hétérozygotes pour cette mutation. Le séquençage d’ITGB3 n’a démontré aucune anomalie.

Harmonic force fields from scaled SCF calculations: Program ASYM40

We report an extended version of our normal coordinate program ASYM40, which may be used to transform Cartesian force constants from ab initio calculations to a force field in nonredundant internal (symmetry) coordinates. When experimental data are available, scale factors for the theoretical force field may then be optimized by least-squares refinement. The alternative of refining an empirical force field to fit a wide variety of data, as with the previous version ASYM20, has been retained. We compare the results of least-squares refinement of the full harmonic force field with least-squares refinement of only the scale factors for an SCF calculated force field and conclude that the latter approach may be useful for large molecules where more sophisticated calculations are impractical. The refinement of scale factors for a theoretical force field is also useful when there are only limited spectroscopic data. The program will accept ab initio calculated force fields from any program that presents Cartesian force constants as output. The program is available through Quantum Chemistry Program Exchange.

General derivative relations for anharmonic force fields

The brace notation, introduced by Allen and Csaszar (1993, J. chem. Phys., 98, 2983), provides a simple and compact way to deal with derivatives of arbitrary non-tensorial quantities. One of its main advantages is that it builds the permutational symmetry of the derivatives directly into the formalism. The brace notation is applied to formulate the general nth-order Cartesian derivatives of internal coordinates, and to provide closed forms for general, nth-order transformation equations of anharmonic force fields, expressed as Taylor series, from internal to Cartesian or normal coordinate spaces.

High resolution photoacoustic overtone spectroscopy of monofluoroacetylene, HCCF, with a titanium:sapphire ring laser

Intracavity photoacoustic overtone spectrum of monofluoroacetylene, HCCF, has been recorded in the wave number region 10 750–14 500 cm−1 with a titanium:sapphire ring laser. The spectrum contains many dense vibration–rotation band systems which can be resolved with Doppler limited resolution. Altogether 58 individual overtone bands have been analyzed rotationally. Many of the observed bands show perturbations of which some have been attributed to anharmonic resonance interactions. A Fermi resonance model based on conventional rectilinear normal coordinate theory has been used to assign vibrationally bands from this work and from earlier studies. Many of the observed vibrational term values and rotational constants can be reproduced well with this model. The results show the importance of the Fermi resonance interactions at the high overtone excitations.

Quartic anharmonic resonances in acetylene molecules

Formulas are derived for the quartic anharmonic resonance coefficients observed to be important between C–H stretching and the combination of one quantum of C≡C stretching and two quanta of H–C≡C bending in a number of acetylene molecules. Examples of this resonance are ν3 with ν2+ν4+ν5 in 12C2H2, ν1 with ν2+2ν5 in 13C2H2, and ν1 with ν2+2ν4 in monofluoroacetylene and monochloroacetylene. The coefficients characterizing the resonances in these examples, which we denote K3,245, K1,255, and K1,244, arise from cubic and quartic terms in the anharmonic force field, in the normal coordinate representation, through second order and first order perturbation treatments respectively, where the second order resonances are calculated by a Van Vleck resonance formalism. The experimentally determined values of these coefficients are compared with values calculated from model anharmonic force fields.

Infrared intensities of furan, pyrrole and thiophene: beyond the double harmonic approximation

Infrared intensities of the fundamental, overtone and combination transitions in furan, pyrrole and thiophene have been calculated using the variational normal coordinate code MULTIMODE. We use pure vibrational wavefunctions, and quartic force fields and cubic dipole moment vector surfaces, generated by density functional theory. The results are compared graphically with second-order perturbation calculations and with relative intensities from experiment for furan and pyrrole.

The vibrations and tunnelling of malonaldehyde on a Moller-Plesset surface

The vibrations and tunnelling motion of malonaldehyde have been studied in their full dimensionality using an internal coordinate path Hamiltonian. In this representation there is one large amplitude internal coordinate s and 3N - 7 (=20) normal coordinates Q which are orthogonal to the large amplitude motion at all points. It is crucial that a high accuracy potential energy surface is used in order to obtain a good representation for the tunneling motion; we use a Moller-Plesset (MP2) surface. Our methodology is variational, that is we diagonalize a sufficiently large matrix in order to obtain the required vibrational levels, so an exact representation for the kinetic energy operator is used. In a harmonic valley representation (s, Q) complete convergence of the normal coordinate motions and the internal coordinate motions has been obtained; for the anharmonic valley in which we use two- and three-body terms in the surface (s, Q(1), Q(2)), we also obtain complete convergence. Our final computed stretching fundamentals are deficient because our potential energy surface is truncated at quartic terms in the normal coordinates, but our lower fundamentals are good.

Baroclinic instability on hot extrasolar planets

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value cal- culations are carried out – focusing on the freely-evolving, adiabatic situation. Using a high- resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot ex- trasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qual- itative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks – but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.