Evolution and breakup of viscous rotating drops

We study the evolution of a viscous fluid drop rotating about a fixed axis at constant angular velocity $Omega$ or constant angular momentum L surrounded by another viscous fluid. The problem is considered in the limit of large Ekman number and small Reynolds number. The analysis is carried out by combining asymptotic analysis and full numerical simulation by means of the boundary element method. We pay special attention to the stability/instability of equilibrium shapes and the possible formation of singularities representing a change in the topology of the fluid domain. When the evolution is at constant $Omega$, depending on its value, drops can take the form of a flat film whose thickness goes to zero in finite time or an elongated filament that extends indefinitely. When evolution takes place at constant L and axial symmetry is imposed, thin films surrounded by a toroidal rim can develop, but the film thickness does not vanish in finite time. When axial symmetry is not imposed and L is sufficiently large, drops break axial symmetry and, depending on the value of L, reach an equilibrium configuration with a 2-fold symmetry or break up into several drops with a 2- or 3-fold symmetry. The mechanism of breakup is also described

Numerical simulation of the evolution of viscous rotating drops

In this contribution we simulate numerically the evolution of a viscous fluid drop rotating about a fixed axis at constant angular velocity ? or constant angular momentum L, surrounded by another viscous fluid. The problem is considered in the limit of large Ekman number and small Reynolds number. In the lecture we will describe the numerical method we have used to solve the PDE system that describes the evolution of the drop (3D boundary element method). We will also present the results we have obtained, paying special attention to the stability/instability of the equilibrium shapes.

Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.

Bending continuous structures with SMAs: a novel robotic fish design

In this paper, we describe our research on bio-inspired locomotion systems using deformable structures and smart materials, concretely shape memory alloys (SMAs). These types of materials allow us to explore the possibility of building motor-less and gear-less robots. A swimming underwater fish-like robot has been developed whose movements are generated using SMAs. These actuators are suitable for bending the continuous backbone of the fish, which in turn causes a change in the curvature of the body. This type of structural arrangement is inspired by fish red muscles, which are mainly recruited during steady swimming for the bending of a flexible but nearly incompressible structure such as the fishbone. This paper reviews the design process of these bio-inspired structures, from the motivations and physiological inspiration to the mechatronics design, control and simulations, leading to actual experimental trials and results. The focus of this work is to present the mechanisms by which standard swimming patterns can be reproduced with the proposed design. Moreover, the performance of the SMA-based actuators’ control in terms of actuation speed and position accuracy is also addressed.

A computer program for the analysis of the dynamic bending-torsion coupling in bridges using a mini-computer

The analysis of modes and natural frequencies is of primary interest in the computation of the response of bridges. In this article the transfer matrix method is applied to this problem to provide a computer code to calculate the natural frequencies and modes of bridge-like structures. The Fortran computer code is suitable for running on small computers and results are presented for a railway bridge.

Hexagonal patterns in a model for rotating convection

We study a model equation that mimics convection under rotation in a fluid with temperature- dependent properties (non-Boussinesq (NB)), high Prandtl number and idealized boundary conditions. It is based on a model equation proposed by Segel [1965] by adding rotation terms that lead to a Kuppers-Lortz instability [Kuppers & Lortz, 1969] and can develop into oscillating hexagons. We perform a weakly nonlinear analysis to find out explicitly the coefficients in the amplitude equation as functions of the rotation rate. These equations describe hexagons and os- cillating hexagons quite well, and include the Busse?Heikes (BH) model [Busse & Heikes, 1980] as a particular case. The sideband instabilities as well as short wavelength instabilities of such hexagonal patterns are discussed and the threshold for oscillating hexagons is determined.

Defect chaos and bursts: Hexagonal rotating convection and the complex Ginzburg-Landau equation

We employ numerical computations of the full Navier-Stokes equations to investigate non-Boussinesq convection in a rotating system using water as the working fluid. We identify two regimes. For weak non- Boussinesq effects the Hopf bifurcation from steady to oscillating (whirling) hexagons is supercritical and typical states exhibit defect chaos that is systematically described by the cubic complex Ginzburg-Landau equation. For stronger non-Boussinesq effects the Hopf bifurcation becomes subcritical and the oscil- lations exhibit localized chaotic bursting, which is modeled by a quintic complex Ginzburg-Landau equation.

Frozen orbits for scientific missions using rotating tethers

We derive a semi-analytic formulation that permits to study the long-term dynamics of fast-rotating inert tethers around planetary satellites. Since space tethers are extensive bodies they generate non-keplerian gravitational forces which depend solely on their mass geometry and attitude, that can be exploited for controlling science orbits. We conclude that rotating tethers modify the geometry of frozen orbits, allowing for lower eccentricity frozen orbits for a wide range of orbital inclination, where the length of the tether becomes a new parameter that the mission analyst may use to shape frozen orbits to tighter operational constraints.

Conceptual design of the beryllium rotating target for the ESS-Bilbao facility

The ESS-Bilbao facility, hosted by the University of the Basque Country (UPV/EHU), envisages the operation of a high-current proton accelerator delivering beams with energies up to 50 MeV. The time-averaged proton current will be 2.25 mA, delivered by 1.5 ms proton pulses with a repetition rate of 20 Hz. This beam will feed a neutron source based upon the Be (p,n) reaction, which will enable the provision of relevant neutron experimentation capabilities. The neutron source baseline concept consists in a rotating beryllium target cooled by water. The target structure will comprise a rotatable disk made of 6061-T6 aluminium alloy holding 20 beryllium plates. Heat dissipation from the target relies upon a distribution of coolant-flow channels. The practical implementation of such a concept is here described with emphasis put on the beryllium plates thermo-mechanical optimization, the chosen coolant distribution system as well as the mechanical behavior of the assembly.

Discussion of "Design Footing Area with Biaxial Bending" by Ramon Jarquio and Victor Jarquio (October, 1983)

The authors present a very interesting criterion for choosing a rectangular foundation.The writers should like to point out that the obtention of minimum area can be reduced to the problem of finding the minimum of x*+y*, subjected to the condition x*.y*=k2 whose solution is evidently x*=y*=k

Influence of the Second Bending Mode on the Response of High-Speed Bridges at Resonance

This paper deals with the assessment of the contribution of the second bending mode to the dynamic behavior of simply supported railway bridges. Traditionally the contributions of modes higher than the fundamental have been considered of little importance for the computation of the magnitudes of interest to structural engineers (vertical deflections, bending moments, etc.). Starting from the dimensionless equations of motion of a simply supported beam subjected to moving loads, the key parameters governing the dynamic behavior are identified. Then, a parametric study over realistic ranges of values of those parameters is conducted, and the influence of the second mode examined in detail. The main purpose is to decide whether the second mode should be taken into account for the determination of the maximum displacement and acceleration in high-speed bridges. In addition, the reasons that cause the contribution of the second bending mode to be relevant in some situations are highlighted, particularly with regard to the computation of the maximum acceleration.

A comparative analysis between FEM and Boundary Methods in bending of Plates

Two different methods of analysis of plate bending, FEM and BM are discussed in this paper. The plate behaviour is assumed to be represented by using the linear thin plate theory where the Poisson-Kirchoff assumption holds. The BM based in a weighted mean square error technique produced good results for the problem of plate bending. The computational effort demanded in the BM is smaller than the one needed in a FEM analysis for the same level of accuracy. The general application of the FEM cannot be matched by the BM. Particularly, different types of geometry (plates of arbitrary geometry) need a similar but not identical treatment in the BM. However, this loss of generality is counterbalanced by the computational efficiency gained in the BM in the solution achievement

Analysis of plate bending by means of high order finite hyperelements

The possibilities and limitations of high order hyperelements in plate bending analysis are discussed. Explicit shape functions for some types of triangular elements are given. These elements are applied to simple cases to assess their computational efficiency.

A C1 finite element for Kirchhoff plate bending

After a short introduction the possibilities and limitations of polynomial simple elements with C1 continuity are discussed with reference to plate bending analysis. A family of this kind of elements is presented.. These elements are applied to simple cases in order to assess their computational efficiency. Finally some conclusions are shown, and future research is also proposed.

Bending reinforcement of timber beams with composite carbon fiber and basalt fiber materials

En el presente trabajo se lleva a cabo un estudio basado en datos obtenidos experimentalmente mediante el ensayo a flexión de vigas de madera de pino silvestre reforzadas con materiales compuestos. Las fibras que componen los tejidos utilizados para la ejecución de los refuerzos son de basalto y de carbono. En el caso de los compuestos de fibra de basalto se aplican en distintos gramajes, y los de carbono en tejido unidireccional y bidireccional. El material compuesto se realizó in situ, simultáneamente a la ejecución del refuerzo. Se aplicaron en una y en dos capas, según el caso, y la forma de colocación fue en ?U?, adhiriéndose al canto inferior y a las caras laterales de la viga mediante resina o mortero epoxi. Se analiza el comportamiento de las vigas según las variables de refuerzo aplicadas y se comparan con los resultados de vigas ensayadas sin reforzar. Con este trabajo queda demostrado el buen funcionamiento del FRP de fibra de basalto aplicado en el refuerzo de vigas de madera y de los tejidos de carbono bidireccionales con respecto a los unidireccionales.