Development of a Uniform Theory of Diffraction for Scattered and Surface Electromagnetic Field Analysis on a Cylinder

Con esta tesis ”Desarrollo de una Teoría Uniforme de la Difracción para el Análisis de los Campos Electromagnéticos Dispersados y Superficiales sobre un Cilindro” hemos iniciado una nueva línea de investigación que trata de responder a la siguiente pregunta: ¿cuál es la impedancia de superficie que describe una estructura de conductor eléctrico perfecto (PEC) convexa recubierta por un material no conductor? Este tipo de estudios tienen interés hoy en día porque ayudan a predecir el campo electromagnético incidente, radiado o que se propaga sobre estructuras metálicas y localmente convexas que se encuentran recubiertas de algún material dieléctrico, o sobre estructuras metálicas con pérdidas, como por ejemplo se necesita en determinadas aplicaciones aeroespaciales, marítimas o automovilísticas. Además, desde un punto de vista teórico, la caracterización de la impedancia de superficie de una estructura PEC recubierta o no por un dieléctrico es una generalización de varias soluciones que tratan ambos tipos de problemas por separado. En esta tesis se desarrolla una teoría uniforme de la difracción (UTD) para analizar el problema canónico del campo electromagnético dispersado y superficial en un cilindro circular eléctricamente grande con una condición de contorno de impedancia (IBC) para frecuencias altas. Construir una solución basada en UTD para este problema canónico es crucial en el desarrollo de un método UTD para el caso más general de una superficie arbitrariamente convexa, mediante el uso del principio de localización de los campos electromagnéticos a altas frecuencias. Esta tesis doctoral se ha llevado a cabo a través de una serie de hitos que se enumeran a continuación, enfatizando las contribuciones a las que ha dado lugar. Inicialmente se realiza una revisión en profundidad del estado del arte de los métodos asintóticos con numerosas referencias. As í, cualquier lector novel puede llegar a conocer la historia de la óptica geométrica (GO) y la teoría geométrica de la difracción (GTD), que dieron lugar al desarrollo de la UTD. Después, se investiga ampliamente la UTD y los trabajos más importantes que pueden encontrarse en la literatura. As í, este capítulo, nos coloca en la posición de afirmar que, hasta donde nosotros conocemos, nadie ha intentado antes llevar a cabo una investigación rigurosa sobre la caracterización de la impedancia de superficie de una estructura PEC recubierta por un material dieléctrico, utilizando para ello la UTD. Primero, se desarrolla una UTD para el problema canónico de la dispersión electromagnética de un cilindro circular eléctricamente grande con una IBC uniforme, cuando es iluminado por una onda plana con incidencia oblicua a frecuencias altas. La solución a este problema canónico se construye a partir de una solución exacta mediante una expansión de autofunciones de propagación radial. Entonces, ésta se convierte en una nueva expansión de autofunciones de propagación circunferencial muy apropiada para cilindros grandes, a través de la transformación de Watson. De esta forma, la expresión del campo se reduce a una integral que se evalúa asintóticamente, para altas frecuencias, de manera uniforme. El resultado se expresa según el trazado de rayos descrito en la UTD. La solución es uniforme porque tiene la importante propiedad de mantenerse continua a lo largo de la región de transición, a ambos lados de la superficie del contorno de sombra. Fuera de la región de transición la solución se reduce al campo incidente y reflejado puramente ópticos en la región iluminada del cilindro, y al campo superficial difractado en la región de sombra. Debido a la IBC el campo dispersado contiene una componente contrapolar a causa de un acoplamiento entre las ondas TEz y TMz (donde z es el eje del cilindro). Esta componente contrapolar desaparece cuando la incidencia es normal al cilindro, y también en la región iluminada cuando la incidencia es oblicua donde el campo se reduce a la solución de GO. La solución UTD presenta una muy buena exactitud cuando se compara numéricamente con una solución de referencia exacta. A continuación, se desarrolla una IBC efectiva para el cálculo del campo electromagnético dispersado en un cilindro circular PEC recubierto por un dieléctrico e iluminado por una onda plana incidiendo oblicuamente. Para ello se derivan dos impedancias de superficie en relación directa con las ondas creeping y de superficie TM y TE que se excitan en un cilindro recubierto por un material no conductor. Las impedancias de superficie TM y TE están acopladas cuando la incidencia es oblicua, y dependen de la geometría del problema y de los números de onda. Además, se ha derivado una impedancia de superficie constante, aunque con diferente valor cuando el observador se encuentra en la zona iluminada o en la zona de sombra. Después, se presenta una solución UTD para el cálculo de la dispersión de una onda plana con incidencia oblicua sobre un cilindro eléctricamente grande y convexo, mediante la generalización del problema canónico correspondiente al cilindro circular. La solución asintótica es uniforme porque se mantiene continua a lo largo de la región de transición, en las inmediaciones del contorno de sombra, y se reduce a la solución de rayos ópticos en la zona iluminada y a la contribución de las ondas de superficie dentro de la zona de sombra, lejos de la región de transición. Cuando se usa cualquier material no conductor se excita una componente contrapolar que tiende a desaparecer cuando la incidencia es normal al cilindro y en la región iluminada. Se discuten ampliamente las limitaciones de las fórmulas para la impedancia de superficie efectiva, y se compara la solución UTD con otras soluciones de referencia, donde se observa una muy buena concordancia. Y en tercer lugar, se presenta una aproximación para una impedancia de superficie efectiva para el cálculo de los campos superficiales en un cilindro circular conductor recubierto por un dieléctrico. Se discuten las principales diferencias que existen entre un cilindro PEC recubierto por un dieléctrico desde un punto de vista riguroso y un cilindro con una IBC. Mientras para un cilindro de impedancia se considera una impedancia de superficie constante o uniforme, para un cilindro conductor recubierto por un dieléctrico se derivan dos impedancias de superficie. Estas impedancias de superficie están asociadas a los modos de ondas creeping TM y TE excitadas en un cilindro, y dependen de la posición y de la orientación del observador y de la fuente. Con esto en mente, se deriva una solución UTD con IBC para los campos superficiales teniendo en cuenta las dependencias de la impedancia de superficie. La expansión asintótica se realiza, mediante la transformación de Watson, sobre la representación en serie de las funciones de Green correspondientes, evitando as í calcular las derivadas de orden superior de las integrales de tipo Fock, y dando lugar a una solución rápida y precisa. En los ejemplos numéricos realizados se observa una muy buena precisión cuando el cilindro y la separación entre el observador y la fuente son grandes. Esta solución, junto con el método de los momentos (MoM), se puede aplicar para el cálculo eficiente del acoplamiento mutuo de grandes arrays conformados de antenas de parches. Los métodos propuestos basados en UTD para el cálculo del campo electromagnético dispersado y superficial sobre un cilindro PEC recubierto de dieléctrico con una IBC efectiva suponen un primer paso hacia la generalización de una solución UTD para superficies metálicas convexas arbitrarias cubiertas por un material no conductor e iluminadas por una fuente electromagnética arbitraria. ABSTRACT With this thesis ”Development of a Uniform Theory of Diffraction for Scattered and Surface Electromagnetic Field Analysis on a Cylinder” we have initiated a line of investigation whose goal is to answer the following question: what is the surface impedance which describes a perfect electric conductor (PEC) convex structure covered by a material coating? These studies are of current and future interest for predicting the electromagnetic (EM) fields incident, radiating or propagating on locally smooth convex parts of highly metallic structures with a material coating, or by a lossy metallic surfaces, as for example in aerospace, maritime and automotive applications. Moreover, from a theoretical point of view, the surface impedance characterization of PEC surfaces with or without a material coating represents a generalization of independent solutions for both type of problems. A uniform geometrical theory of diffraction (UTD) is developed in this thesis for analyzing the canonical problem of EM scattered and surface field by an electrically large circular cylinder with an impedance boundary condition (IBC) in the high frequency regime, by means of a surface impedance characterization. The construction of a UTD solution for this canonical problem is crucial for the development of the corresponding UTD solution for the more general case of an arbitrary smooth convex surface, via the principle of the localization of high frequency EM fields. The development of the present doctoral thesis has been carried out through a series of landmarks that are enumerated as follows, emphasizing the main contributions that this work has given rise to. Initially, a profound revision is made in the state of art of asymptotic methods where numerous references are given. Thus, any reader may know the history of geometrical optics (GO) and geometrical theory of diffraction (GTD), which led to the development of UTD. Then, the UTD is deeply investigated and the main studies which are found in the literature are shown. This chapter situates us in the position to state that, as far as we know, nobody has attempted before to perform a rigorous research about the surface impedance characterization for material-coated PEC convex structures via UTD. First, a UTD solution is developed for the canonical problem of the EM scattering by an electrically large circular cylinder with a uniform IBC, when it is illuminated by an obliquely incident high frequency plane wave. A solution to this canonical problem is first constructed in terms of an exact formulation involving a radially propagating eigenfunction expansion. The latter is converted into a circumferentially propagating eigenfunction expansion suited for large cylinders, via the Watson transformation, which is expressed as an integral that is subsequently evaluated asymptotically, for high frequencies, in a uniform manner. The resulting solution is then expressed in the desired UTD ray form. This solution is uniform in the sense that it has the important property that it remains continuous across the transition region on either side of the surface shadow boundary. Outside the shadow boundary transition region it recovers the purely ray optical incident and reflected ray fields on the deep lit side of the shadow boundary and to the modal surface diffracted ray fields on the deep shadow side. The scattered field is seen to have a cross-polarized component due to the coupling between the TEz and TMz waves (where z is the cylinder axis) resulting from the IBC. Such cross-polarization vanishes for normal incidence on the cylinder, and also in the deep lit region for oblique incidence where it properly reduces to the GO or ray optical solution. This UTD solution is shown to be very accurate by a numerical comparison with an exact reference solution. Then, an effective IBC is developed for the EM scattered field on a coated PEC circular cylinder illuminated by an obliquely incident plane wave. Two surface impedances are derived in a direct relation with the TM and TE surface and creeping wave modes excited on a coated cylinder. The TM and TE surface impedances are coupled at oblique incidence, and depend on the geometry of the problem and the wave numbers. Nevertheless, a constant surface impedance is found, although with a different value when the observation point lays in the lit or in the shadow region. Then, a UTD solution for the scattering of an obliquely incident plane wave on an electrically large smooth convex coated PEC cylinder is introduced, via a generalization of the canonical circular cylinder problem. The asymptotic solution is uniform because it remains continuous across the transition region, in the vicinity of the shadow boundary, and it recovers the ray optical solution in the deep lit region and the creeping wave formulation within the deep shadow region. When a coating is present a cross-polar field term is excited, which vanishes at normal incidence and in the deep lit region. The limitations of the effective surface impedance formulas are discussed, and the UTD solution is compared with some reference solutions where a very good agreement is met. And in third place, an effective surface impedance approach is introduced for determining surface fields on an electrically large coated metallic circular cylinder. Differences in analysis of rigorouslytreated coated metallic cylinders and cylinders with an IBC are discussed. While for the impedance cylinder case a single constant or uniform surface impedance is considered, for the coated metallic cylinder case two surface impedances are derived. These are associated with the TM and TE creeping wave modes excited on a cylinder and depend on observation and source positions and orientations. With this in mind, a UTD based method with IBC is derived for the surface fields by taking into account the surface impedance variation. The asymptotic expansion is performed, via the Watson transformation, over the appropriate series representation of the Green’s functions, thus avoiding higher-order derivatives of Fock-type integrals, and yielding a fast and an accurate solution. Numerical examples reveal a very good accuracy for large cylinders when the separation between the observation and the source point is large. Thus, this solution could be efficiently applied in mutual coupling analysis, along with the method of moments (MoM), of large conformal microstrip array antennas. The proposed UTD methods for scattered and surface EM field analysis on a coated PEC cylinder with an effective IBC are considered the first steps toward the generalization of a UTD solution for large arbitrarily convex smooth metallic surfaces covered by a material coating and illuminated by an arbitrary EM source.

Surface tension and microgravity

The behaviour of confined liquids on board an orbiting spacecraft is mainly driven by surface tension phenomena, which cause an apparently anomalous response of the liquid when compared with the behaviour that can be observed on an Earth laboratory provided that the amount of liquid is high enough. The reason is that in an orbiting spacecraft the different inertial forces acting on the bulk of the liquid are almost zero, causing thus capillary forces to be the dominant ones. Of course, since gravity forces are proportional to the liquid volume, whereas surface tension forces are proportional to the liquid surface, there are situations on Earth where capillarity can be the dominant effect, as it happens when very small volume liquid samples are considered. However, work with small size samples may require the use of sophisticated optical devices. Leaving aside the neutral buoyancy technique, a way of handling large liquid interfaces is by using drop towers, where the sample falls subjected to the action of Earth's gravity. This approach is suitable when the characteristic time of the problem under consideration is much smaller than the drop time. In this work the transformation of an out-of-use chimney into a drop tower is presented. Because of the miniaturization, hardiness and low cost of current electronic devices, a drop tower can be used as an inexpensive tool for undergraduate students to experimentally analyse a large variety of surface tension driven phenomena.

Dynamic earth pressures against a retaining wall caused by Rayleigh waves

An approximate procedure for studying harmonic soil-structure interaction problems is presented. The presence of Rayleigh waves is considered and the resulting governing equations of the dynamic soil-structure system are solved in the time domain. With this method the transient and steady states of a vibratory motion and also the nonlinear behaviour of the soil can be studied. As an example, the dynamic earth pressure against a rigid retaining wall is investigated. The loads are assumed to be harmonic Rayleigh waves with both static and dynamic surface surcharges. The dependence of the results on the excitation frequency is shown.

Hydrodynamic forces on heave plates for offshore systems oscillating close to the seabed or the free surface

La energía eólica marina es uno de los recursos energéticos con mayor proyección pudiendo contribuir a reducir el consumo de combustibles fósiles y a cubrir la demanda de energía en todo el mundo. El concepto de aerogenerador marino está basado en estructuras fijas como jackets o en plataformas flotantes, ya sea una semisumergible o una TLP. Se espera que la energía eólica offshore juegue un papel importante en el perfil de producción energética de los próximos años; por tanto, las turbinas eólicas deben hacerse más fables y rentables para ser competitivas frente a otras fuentes de energía. Las estructuras flotantes pueden experimentar movimientos resonantes en estados de la mar con largos períodos de oleaje. Estos movimientos disminuyen su operatividad y pueden causar daños en los componentes eléctricos de las turbinas y en las palas, también en los risers y moorings. La respuesta de la componente vertical del movimiento puede reducirse mediante diferentes actuaciones: (1) aumentando la amortiguación del sistema, (2) manteniendo el período del movimiento vertical fuera del rango de la energía de la ola, y (3) reduciendo las fuerzas de excitación verticales. Un ejemplo típico para llevar a cabo esta reducción son las "Heave Plates". Las heave plates son placas que se utilizan en la industria offshore debido a sus características hidrodinámicas, ya que aumentan la masa añadida y la amortiguación del sistema. En un análisis hidrodinámico convencional, se considera una estructura sometida a un oleaje con determinadas características y se evalúan las cargas lineales usando la teoría potencial. El amortiguamiento viscoso, que juega un papel crucial en la respuesta en resonancia del sistema, es un dato de entrada para el análisis. La tesis se centra principalmente en la predicción del amortiguamiento viscoso y de la masa añadida de las heave plates usadas en las turbinas eólicas flotantes. En los cálculos, las fuerzas hidrodinámicas se han obtenido con el f n de estudiar cómo los coeficientes hidrodinámicos de masa añadida5 y amortiguamiento varían con el número de KC, que caracteriza la amplitud del movimiento respecto al diámetro del disco. Por otra parte, se ha investigado la influencia de la distancia media de la ‘heave plate’ a la superficie libre o al fondo del mar, sobre los coeficientes hidrodinámicos. En este proceso, un nuevo modelo que describe el trabajo realizado por la amortiguación en función de la enstrofía, es descrito en el presente documento. Este nuevo enfoque es capaz de proporcionar una correlación directa entre el desprendimiento local de vorticidad y la fuerza de amortiguación global. El análisis también incluye el estudio de los efectos de la geometría de la heave plate, y examina la sensibilidad de los coeficientes hidrodinámicos al incluir porosidad en ésta. Un diseño novedoso de una heave plate, basado en la teoría fractal, también fue analizado experimentalmente y comparado con datos experimentales obtenidos por otros autores. Para la resolución de las ecuaciones de Navier Stokes se ha usado un solver basado en el método de volúmenes finitos. El solver usa las librerías de OpenFOAM (Open source Field Operation And Manipulation), para resolver un problema multifásico e incompresible, usando la técnica VOF (volume of fluid) que permite capturar el movimiento de la superficie libre. Los resultados numéricos han sido comparados con resultados experimentales llevados a cabo en el Canal del Ensayos Hidrodinámicos (CEHINAV) de la Universidad Politécnica de Madrid y en el Canal de Experiencias Hidrodinámicas (CEHIPAR) en Madrid, al igual que con otros experimentos realizados en la Escuela de Ingeniería Mecánica de la Universidad de Western Australia. Los principales resultados se presentan a continuación: 1. Para pequeños valores de KC, los coeficientes hidrodinámicos de masa añadida y amortiguamiento incrementan su valor a medida que el disco se aproxima al fondo marino. Para los casos cuando el disco oscila cerca de la superficie libre, la dependencia de los coeficientes hidrodinámicos es más fuerte por la influencia del movimiento de la superficie libre. 2. Los casos analizados muestran la existencia de un valor crítico de KC, donde la tendencia de los coeficientes hidrodinámicos se ve alterada. Dicho valor crítico depende de la distancia al fondo marino o a la superficie libre. 3. El comportamiento físico del flujo, para valores de KC cercanos a su valor crítico ha sido estudiado mediante el análisis del campo de vorticidad. 4. Introducir porosidad al disco, reduce la masa añadida para los valores de KC estudiados, pero se ha encontrado que la porosidad incrementa el valor del coeficiente de amortiguamiento cuando se incrementa la amplitud del movimiento, logrando un máximo de damping para un disco con 10% de porosidad. 5. Los resultados numéricos y experimentales para los discos con faldón, muestran que usar este tipo de geometrías incrementa la masa añadida cuando se compara con el disco sólido, pero reduce considerablemente el coeficiente de amortiguamiento. 6. Un diseño novedoso de heave plate basado en la teoría fractal ha sido experimentalmente estudiado a diferentes calados y comparado con datos experimentales obtenidos por otro autores. Los resultados muestran un comportamiento incierto de los coeficientes y por tanto este diseño debería ser estudiado más a fondo. ABSTRACT Offshore wind energy is one of the promising resources which can reduce the fossil fuel energy consumption and cover worldwide energy demands. Offshore wind turbine concepts are based on either a fixed structure as a jacket or a floating offshore platform like a semisubmersible, spar or tension leg platform. Floating offshore wind turbines have the potential to be an important part of the energy production profile in the coming years. In order to accomplish this wind integration, these wind turbines need to be made more reliable and cost efficient to be competitive with other sources of energy. Floating offshore artifacts, such oil rings and wind turbines, may experience resonant heave motions in sea states with long peak periods. These heave resonances may increase the system downtime and cause damage on the system components and as well as on risers and mooring systems. The heave resonant response may be reduced by different means: (1) increasing the damping of the system, (2) keeping the natural heave period outside the range of the wave energy, and (3) reducing the heave excitation forces. A typical example to accomplish this reduction are “Heave Plates”. Heave plates are used in the offshore industry due to their hydrodynamic characteristics, i.e., increased added mass and damping. Conventional offshore hydrodynamic analysis considers a structure in waves, and evaluates the linear and nonlinear loads using potential theory. Viscous damping, which is expected to play a crucial role in the resonant response, is an empirical input to the analysis, and is not explicitly calculated. The present research has been mainly focused on the prediction of viscous damping and added mass of floating offshore wind turbine heave plates. In the calculations, the hydrodynamic forces have been measured in order to compute how the hydrodynamic coefficients of added mass1 and damping vary with the KC number, which characterises the amplitude of heave motion relative to the diameter of the disc. In addition, the influence on the hydrodynamic coefficients when the heave plate is oscillating close to the free surface or the seabed has been investigated. In this process, a new model describing the work done by damping in terms of the flow enstrophy, is described herein. This new approach is able to provide a direct correlation between the local vortex shedding processes and the global damping force. The analysis also includes the study of different edges geometry, and examines the sensitivity of the damping and added mass coefficients to the porosity of the plate. A novel porous heave plate based on fractal theory has also been proposed, tested experimentally and compared with experimental data obtained by other authors for plates with similar porosity. A numerical solver of Navier Stokes equations, based on the finite volume technique has been applied. It uses the open-source libraries of OpenFOAM (Open source Field Operation And Manipulation), to solve 2 incompressible, isothermal immiscible fluids using a VOF (volume of fluid) phase-fraction based interface capturing approach, with optional mesh motion and mesh topology changes including adaptive re-meshing. Numerical results have been compared with experiments conducted at Technical University of Madrid (CEHINAV) and CEHIPAR model basins in Madrid and with others performed at School of Mechanical Engineering in The University of Western Australia. A brief summary of main results are presented below: 1. At low KC numbers, a systematic increase in added mass and damping, corresponding to an increase in the seabed proximity, is observed. Specifically, for the cases when the heave plate is oscillating closer to the free surface, the dependence of the hydrodynamic coefficients is strongly influenced by the free surface. 2. As seen in experiments, a critical KC, where the linear trend of the hydrodynamic coefficients with KC is disrupted and that depends on the seabed or free surface distance, has been found. 3. The physical behavior of the flow around the critical KC has been explained through an analysis of the flow vorticity field. 4. The porosity of the heave plates reduces the added mass for the studied porosity at all KC numbers, but the porous heave plates are found to increase the damping coefficient with increasing amplitude of oscillation, achieving a maximum damping coefficient for the heave plate with 10% porosity in the entire KC range. 5. Another concept taken into account in this work has been the heave plates with flaps. Numerical and experimental results show that using discs with flaps will increase added mass when compared to the plain plate but may also significantly reduce damping. 6. A novel heave plate design based on fractal theory has tested experimentally for different submergences and compared with experimental data obtained by other authors for porous plates. Results show an unclear behavior in the coefficients and should be studied further. Future work is necessary in order to address a series of open questions focusing on 3D effects, optimization of the heave plates shapes, etc.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

Field theory in superfluid 3He: What are the lessons for particle physics, gravity, and high-temperature superconductivity?

There are several classes of homogeneous Fermi systems that are characterized by the topology of the energy spectrum of fermionic quasiparticles: (i) gapless systems with a Fermi surface, (ii) systems with a gap in their spectrum, (iii) gapless systems with topologically stable point nodes (Fermi points), and (iv) gapless systems with topologically unstable lines of nodes (Fermi lines). Superfluid 3He-A and electroweak vacuum belong to the universality class 3. The fermionic quasiparticles (particles) in this class are chiral: they are left-handed or right-handed. The collective bosonic modes of systems of class 3 are the effective gauge and gravitational fields. The great advantage of superfluid 3He-A is that we can perform experiments by using this condensed matter and thereby simulate many phenomena in high energy physics, including axial anomaly, baryoproduction, and magnetogenesis. 3He-A textures induce a nontrivial effective metrics of the space, where the free quasiparticles move along geodesics. With 3He-A one can simulate event horizons, Hawking radiation, rotating vacuum, etc. High-temperature superconductors are believed to belong to class 4. They have gapless fermionic quasiparticles with a “relativistic” spectrum close to gap nodes, which allows application of ideas developed for superfluid 3He-A.

Imaging of spin waves in atomically designed nanomagnets

The spin dynamics of all ferromagnetic materials are governed by two types of collective phenomenon: spin waves and domain walls. The fundamental processes underlying these collective modes, such as exchange interactions and magnetic anisotropy, all originate at the atomic scale. However, conventional probing techniques based on neutron1 and photon scattering2 provide high resolution in reciprocal space, and thereby poor spatial resolution. Here we present direct imaging of standing spin waves in individual chains of ferromagnetically coupled S = 2 Fe atoms, assembled one by one on a Cu2N surface using a scanning tunnelling microscope. We are able to map the spin dynamics of these designer nanomagnets with atomic resolution in two complementary ways. First, atom-to-atom variations of the amplitude of the quantized spin-wave excitations are probed using inelastic electron tunnelling spectroscopy. Second, we observe slow stochastic switching between two opposite magnetization states3, 4, whose rate varies strongly depending on the location of the tip along the chain. Our observations, combined with model calculations, reveal that switches of the chain are initiated by a spin-wave excited state that has its antinodes at the edges of the chain, followed by a domain wall shifting through the chain from one end to the other. This approach opens the way towards atomic-scale imaging of other types of spin excitation, such as spinon pairs and fractional end-states5, 6, in engineered spin chains.

Magnetic charachterization of surface sediments in the South Atlantic

Surface sediment samples representative for the tropical and subtropical South Atlantic (15°N to 40°S) were investigated by isothermal magnetic methods to delineate magnetic mineral distribution patterns and to identify their predominant Holocene climatic and oceanographic controls. Individual parameters reveal distinct, yet frequently overlapping, regional sedimentation characteristics. A probabilistic ('fuzzy c-means') cluster analysis was applied to five concentration independent magnetic properties assessing magnetite to hematite ratios and diagnostic of bulk and fine-particle magnetite grain size and coercivity spectra. The resultant 10 cluster structures establish an oceanwide magnetic sediment classification scheme tracing the major terrigenous eolian and fluvial fluxes, authigenic biogenic magnetite accumulation in high-productivity areas, transport by ocean current systems, and effects of bottom water velocity on depositional regimes. Distinct dissimilarities in magnetic mineral inventories between the eastern and western basins of the South Atlantic reflect prominent contrasts of both oceanic and continental influences.

(Table 2) Planktic foraminifera and their diversity indices of sediment surface samples from the Atlantic Ocean

Species distribution patterns in planktonic foraminiferal assemblages are fundamental to the understanding of the determinants of their ecology. Until now, data used to identify such distribution patterns was mainly acquired using the standard >150 µm sieve size. However, given that assemblage shell size-range in planktonic foraminifera is not constant, this data acquisition practice could introduce artefacts in the distributional data. Here, we investigated the link between assemblage shell size-range and diversity in Recent planktonic foraminifera by analysing multiple sieve-size fractions in 12 samples spanning all bioprovinces of the Atlantic Ocean. Using five diversity indices covering various aspects of community structure, we found that counts from the >63 µm fraction in polar oceans and the >125 µm elsewhere sufficiently approximate maximum diversity in all Recent assemblages. Diversity values based on counts from the >150 µm fraction significantly underestimate maximum diversity in the polar and surprisingly also in the tropical provinces. Although the new methodology changes the shape of the diversity/sea-surface temperature (SST) relationship, its strength appears unaffected. Our analysis reveals that increasing diversity in planktonic foraminiferal assemblages is coupled with a progressive addition of larger species that have distinct, offset shell-size distributions. Thus, the previously documented increase in overall assemblage shell size-range towards lower latitudes is linked to an expanding shell-size disparity between species from the same locality. This observation supports the idea that diversity and shell size-range disparity in foraminiferal assemblages are the result of niche separation. Increasing SST leads to enhanced surface water stratification and results in vertical niche separation, which permits ecological specialisation. Specific deviations from the overall diversity and shell-size disparity latitudinal pattern are seen in regions of surface-water instability, indicating that coupled shell-size and diversity measurements could be used to reconstruct water column structures of past oceans.

Sea-surface temperature reconstruction for sediment core GeoB1023-5

Sediment core GeoB 1023-5 from the eastern South Atlantic was investigated at high temporal resolution for variations of sea-surface temperature (SST) during the past 22 kyr, using the alkenone (UK'37) method. SSTs increased by 3.5°C from about 18°C during the Last Ice Age (21±2 cal kyr BP) to about 21.5°C at 14.5 cal kyr BP. This warming trend associated with the deglaciation phase was followed by a cooling event with lowest SSTs near 20°C, persisting for about 1000 years between 13 and 12 cal kyr BP. The SSTs then continued to increase to about 22.5°C at the Holocene climatic optimum at 7 cal kyr BP, and decreased again during the Late Holocene to a core-top value of 19.8°C that is comparable to modern annual mean SST values. When compared with alkenone SST records from the eastern North Atlantic, our SST record indicates continuous warming throughout the deglaciation phase in the Benguela Current, while its northern counterpart, the Canary Current, experienced prominent cooling during 'Heinrich Event 1' (H1). On the other hand, for the time period corresponding to the 'Younger Dryas' (YD) cooling event, the Benguela SST record exhibits a cold-temperature interval that corresponds to that observed in the eastern North Atlantic SST records. This observation suggests that interhemispheric climate response in Atlantic eastern boundary current systems was different with respect to the two abrupt climate events associated with Termination I. For the H1, the eastern South Atlantic SST record strongly supports the hypothesis that an 'anti-phase' thermal behavior in South Atlantic surface waters was forced by the slowdown of the North Atlantic Deep Water formation during cold spells in the North Atlantic. In contrast, the abrupt cooling in the eastern South Atlantic coincident with the YD period was probably induced by more vigorous global atmospheric circulation, enhancing the upwelling intensity in both eastern boundary current systems. This atmospheric control may have overridden any effect caused by changes in thermohaline circulation on the South Atlantic SSTs during the YD, which leads to the assumption that the thermohaline circulation was already much closer to its interglacial mode during the YD than during the H1.

Sea-surface temperature reconstruction in the Southwest Pacific

Uniquely in the Southern Hemisphere the New Zealand micro-continent spans the interface between a subtropical gyre and the Subantarctic Circumpolar Current. Its 20° latitudinal extent includes a complex of submerged plateaux, ridges, saddles and basins which, in the present interglacial, are partial barriers to circulation and steer the Subtropical (STF) and Subantarctic (SAF) fronts. This configuration offers a singular opportunity to assess the influence of bottom topography on oceanic circulation through Pleistocene glacial - interglacial (G/I) cycles, its effect on the location and strength of the fronts, and its ability to generate significant differences in mixed layer thermal history over short distances. For this study we use new planktic foraminiferal based sea-surface temperature (SST) estimates spanning the past 1 million years from a latitudinal transect of four deep ocean drilling sites. We conclude that: 1. the effect of the New Zealand landmass was to deflect the water masses south around the bathymetric impediments; 2. the effect of a shallow submerged ridge on the down-current side (Chatham Rise), was to dynamically trap the STF along its crest, in stark contrast to the usual glacial-interglacial (G-I) meridional migration that occurs in the open ocean; 3. the effect of more deeply submerged, downstream plateaux (Campbell, Bounty) was to dynamically trap the SAF along its steep southeastern margin; 4. the effects of saddles across the submarine plateaux was to facilitate the development of jets of subtropical and subantarctic surface water through the fronts, forming localized downstream gyres or eddies during different phases in the G-I climate cycles; 5. the deep Pukaki Saddle across the Campbell-Bounty Plateaux guided a branch of the SAF to flow northwards during each glacial, to form a strong gyre of circumpolar surface water in the Bounty Trough, especially during the mid-Pleistocene Climate Transition (MIS 22-16) when exceptionally high SST gradients existed across the STF; 6. the shallower Mernoo Saddle, at the western end of the Chatham Rise, provided a conduit for subtropical water to jet southwards across the STF in the warmest interglacial peaks (MIS 11, 5.5) and for subantarctic water to flow northwards during glacials; 7. although subtropical or subantarctic drivers can prevail at a particular phase of a G-I cycles, it appears that the Antarctic Circumpolar Current is the main influence on the regional hydrography. Thus complex submarine topography can affect distinct differences in the climate records over short distances with implications for using such records in interpreting global or regional trends. Conversely, the local topography can amplify the paleoclimate record in different ways in different places, thus enhancing its value for the study of more minor paleoceanographic influences that elsewhere are more difficult to detect. Such sites include DSDP 594, which like some other Southern Ocean sites, has the typical late Pleistocene asymmetrical saw-tooth G-I climate pattern transformed to a gap-tooth pattern of quasi-symmetrical interglacial spikes that interrupt extended periods of minimum glacial temperatures.

Sea surface temperatures derived from alkenones from Late Quaternary sediments of the East Equatorial South Atlantic

Angola Basin and Walvis Ridge records of past sea surface temperatures (SST) derived from the alkenone Uk 37 index are used to reconstruct the surface circulation in the east equatorial South Atlantic for the last 200,000 years. Comparison of SST estimates from surface sediments between 5° and 20°S with modern SST data suggests that the alkenone temperatures represent annual mean values of the surface mixed layer. Alkenone-derived temperatures for the warm climatic maxima of the Holocene and the penultimate interglacial are 1 to 4°C higher than latest Holocene values. All records show glacial to interglacial differences of about 3.5°C in annual mean SST, which is about 1.5°C greater than the difference estimated by CLIMAP (1981) for the eastern Angola Basin. At the Walvis Ridge, significant SST variance is observed at all of the Earth's orbital periodicities. SST records from the Angola Basin vary predominantly at 23- and 100-kyr periodicities. For the precessional cycle, SST changes at the Walvis Ridge correspond to variations of boreal summer insolation over Africa and lead ice volume changes, suggesting that the east equatorial South Atlantic is sensitive to African monsoon intensity via trade-wind zonality. Angola Basin SST records lag those from the Walvis Ridge and the equatorial Atlantic by about 3 kyr. The comparison of Angola Basin and Walvis Ridge SST records implies that the Angola-Benguela Front (ABF) (currently at about 14-16°S) has remained fairly stationary between 12° and 20°S (the limits of our cores) during the last two glacial-interglacial cycles. The temperature contrast associated with the ABF exhibits a periodic 23-kyr variability which is coherent with changes in boreal summer insolation over Africa. These observations suggest that surface waters north of the present ABF have not directly responded to monsoon-modulated changes in the trade-wind vector, that the central field of zonally directed trades in the southern hemisphere was not shifted or extended northward by several degrees of latitude during glacials, and that a cyclonic gyre circulation has existed in the east equatorial South Atlantic over the last 200,000 years. This scenario contradicts former assumptions of glacial intensification of the Benguela Current into the eastern Angola Basin and increased coastal upwelling off Angola.

Nd isotope ratios of surface sediments from the Pacific Ocean (Table 2)

The neodymium isotopic composition of the silicate fraction of Holocene pelagic sediments from the North Pacific define two provinces: a central North Pacific province characterized by unradiogenic and remarkably homogeneous end (-10.2 +/- 0.5) and a narrow circum-Pacific marginal province characterized by more radiogenic and variable end (-4.2 +/- 3.8). The silicate fraction in the central North Pacific is exclusively eolian; based on prevailing wind patterns, meteorological data, and neodymium isotopic data, the only significant sediment source is Chinese loess. Leaching experiments on Chinese loess confirm that leachable Nd is isotopically indistinguishable from bulk and residual silicate Nd. Silicates in the circum-North Pacific marginal province comprise eolian loess, volcanic ash, and hemipelagic sediments derived from volcanic arcs. A compilation of Pacific seawater and Mn nodule epsilon-Nd data shows no clear spatial variation except for a general decrease from surface to deep waters from -3 to -4 and slightly lower epsilon-Nd in bottom waters along the western North Pacific due to the incursion of Antarctic Bottom Water. The relative homogeneity of bottom water epsilon-Nd, which contrasts sharply with the distinctive variation in sediment epsilon-Nd, plus the large difference between the average end of bottom waters and the central North Pacific eolian silicates (-4 vs. -10), suggests that any contribution of REE to seawater from eolian materials is insignificant. Furthermore, leaching of REE from eolian particles as they sink though the water column must be insignificant because Nd in shallow waters is more radiogenic than Nd in deeper waters. That there is no contrast in the Nd isotopic composition of bottom waters that overlie the central and marginal sediment provinces suggests that the ash and hemipelagic sediments derived from Pacific rim volcanic arcs also contribute minimal REE to seawater. The elimination of eolian, ash, and hemipelagic sediments leaves only near-shore riverine particulates as a possibly significant particulate source of REE to seawater.