Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation-this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL. (C) 2015 AIP Publishing LLC.

Detection and estimation of liquid flow through a pipe in a tissue-like object with ultrasound-assisted diffuse correlation spectroscopy

Using coherent light interrogating a turbid object perturbed by a focused ultrasound (US) beam, we demonstrate localized measurement of dynamics in the focal region, termed the region-of-interest (ROI), from the decay of the modulation in intensity autocorrelation of light. When the ROI contains a pipe flow, the decay is shown to be sensitive to the average flow velocity from which the mean-squared displacement (MSD) of the scattering centers in the flow can be estimated. While the MSD estimated is seen to be an order of magnitude higher than that obtainable through the usual diffusing wave spectroscopy (DWS) without the US, it is seen to be more accurate as verified by the volume flow estimated from it. It is further observed that, whereas the MSD from the localized measurement grows with time as tau(alpha) with alpha approximate to 1.65, without using the US, a is seen to be much less. Moreover, with the local measurement, this super-diffusive nature of the pipe flow is seen to persist longer, i.e., over a wider range of initial tau, than with the unassisted DWS. The reason for the super-diffusivity of flow, i.e., alpha < 2, in the ROI is the presence of a fluctuating (thermodynamically nonequilibrium) component in the dynamics induced by the US forcing. Beyond this initial range, both methods measure MSDs that rise linearly with time, indicating that ballistic and near-ballistic photons hardly capture anything beyond the background Brownian motion. (C) 2015 Optical Society of America

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

A characteristics study on the performance of a 2-stage light gas gun

In order to obtain an overall and systematic understanding of the performance of a two-stage light gas gun (TLGG), a numerical code to simulate the process occurring in a gun shot is advanced based on the quasi-one-dimensional unsteady equations of motion with the real gas effect,;friction and heat transfer taken into account in a characteristic formulation for both driver and propellant gas. Comparisons of projectile velocities and projectile pressures along the barrel with experimental results from JET (Joint European Tons) and with computational data got by the Lagrangian method indicate that this code can provide results with good accuracy over a wide range of gun geometry and loading conditions.

The motion of solid particles in a liquid with vertical temperature-gradient

A ground-experiment study on the motions of solid particles in liquid media with vertical temperature gradient is performed in this paper. The movement of solid spheres toward the heating end of a close cell is observed. The behavior and features of the motions examined are quite similar to thermocapillary migration of bubbles and drops in a liquid. The motion velocities of particles measured are about 10(-3) to 10(-4) mm\s. The velocity is compared with the velocity of particles floated in two liquid media. The physical mechanism of motion is explored.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

A Comparison of the Range and Movements of Acmaea digitalis Eschscholtz, Acmaea scutum Eschscholtz, Acmaea limatula Carpenter and Lottia gigantea Sowerby

Ten limpets (4 Acmaea digitalis , 4 Acmaea scutum, 1 Acmaea limatula, and 1 Lottia gigantea) were marked and their movements observed over a thirteen day period. Recordings of positions were made on a map, and the path of each was drawn on the map from day to day. Acmaea digitalis showed the greatest range, mostly in a vertical direction, and moved usually at night during high tide. Acmaea scutum showed a more limited range in a horizontal direction, and moved both day and night during high tide. Acmaea limatula had a horizontal range similar to A. scutum,, but exhibited no movement during the day time. Lottia gigantea had the most restricted range of any limpet studied, and moved only at night during high tide. This is a student paper done for a University of California Berkeley Zoology class. Since UCB didn't have its own marine lab at the time, it rented space at Hopkins Marine Station where this work was done. Gene Haderlie went on to earn his Ph.D. from Berkeley and later became a Professor at the Naval Post Graduate School in Monterey. (PDF contains 23 pages)

Distribution of skipjack in the Pacific Ocean, based on records of incidental catches of the Japanese longline tuna fishery

ENGLISH: All available longline data on skipjack captured in the Pacific Ocean by Japanese research vessels (1949-1965) and from incidental skipjack catches by Japanese commercial vessels (1956-1964) were analyzed. As skipjack are not specifically sought by longline vessels, the data are limited. Considering this it was found that: longline gear captures skipjack of wider size-range and is more selective for larger skipjack than conventional fishing methods, i.e. pole-and-line and purse-seine; skipjack are widely and almost continuously distributed across the Pacific; throughout the year average hook-rates are greater in the southeastern Pacific than in the northwestern Pacific; areas of high hook-rate shift south during the second and third quarters and north during the first and fourth quarters; in the western Pacific the north-south range of the catch distribution was greatest in the first and fourth quarters; skipjack hook-rates are relatively high in the northwestern Pacific east of Japan only during the first and fourth quarters; the highest hook-rates were recorded in extensive areas along the equator (from lO°N to 20°8 between approximately 155°W-100°W); generally more skipjack were captured by research longline gear in water temperature ranges approaching both the upper and lower temperature limits of skipjack distribution (18-21C and 26-28C), than is the case in surface skipjack fisheries; tentative comparisons of longline skipjack catch distributions with Pacific current systems, suggests low skipjack abundance in both North Pacific Central and North Pacific Equatorial water; the sex ratio was 95 males : 63 females in a small sample of skipjack examined; longlines capture skipjack of three, and possibly more, age groups; in skipjack size-composition samples studied, the smaller modal group (65 cm) observed in January-March in the northwestern Pacific (1600E-180oE and 20oN-45°N) corresponds in size to the larger modal group appearing in the late-summer surface fishery off the Izu-Bonin Islands southeast of Japan, and also compares in modal size to the skipjack taken in the Hawaiian fishery in spring time; the analysis of skipjack catches by hook position on the longline and by death-rate studies, indicates that part of the catch is made while the gear is in motion near the surface, and a lesser part of the catch is made when the gear is stabilized at a depth of 70 to 140 m. A brief discussion is given, in the light of new information presented, on several hypotheses by other authors concerning the population structure and migration of skipjack in the Pacific Ocean. SPANISH: Se analizaron todos los datos disponibles de la pesca con palangre de barriletes capturados en el Océano Pacífico por barcos japoneses de investigación (1949-1965) y por las capturas incidentales de los barcos comerciales japoneses (1956-1964). Como los barcos palangreros específicamente, no persiguen al barrilete, los datos son limitados. Considerando ésto, se encontró: que el arte palangrero obtiene barriletes con una distribución más amplia de tallas, y es más selectivo en cuanto a los barriletes de mayor talla, que los métodos convencionales de pesca, Le. cañas de pescar y redes de cerco; el barrilete se encuentra amplia y casi continuamente distribuido a través del Pacífico; en todo el año, las tasas promedio de captura por anzuelo son superiores en el Pacífico sudoriental que las del Pacífico noroeste; las áreas con una tasa alta de captura por anzuelo, se cambian hacia el sur durante los trimestres segundo y tercero, y durante los trimestres primero y cuarto hacia el norte; en el Pacífico occidental la amplitud de la distribución de captura norte-sur, fue superior en los trimestres primero y cuarto; las tasas de captura por anzuelo de barrilete, son relativamente altas en el Pacífico noroeste al este del Japón, únicamente durante los trimestres primero y cuarto; las tasas de captura por anzuelo más altas fueron registradas en extensas áreas a lo largo del ecuador (desde los 10°N hasta los 20°S, aproximadamente entre los 155°W-100°W) ; generalmente las artes palangreras de investigación capturaron más barrilete en aguas en las que la temperatura se aproximaba a los límites más altos o bajos de la temperatura en la distribución del barrilete (18-21 C y 26-28 C), que en el caso de la pesca superficial de barrilete; las comparaciones tentativas de la captura de barrilete con palangre, con el sistema de las corrientes del Pacífico, sugieren una abundancia inferior de barrilete tanto en las aguas del Pacífico central del norte como en las del Pacífico ecuatorial del norte; la proporcíon sexual examinada en una pequeña muestra de barriletes, fue de 95 machos y 63 hembras; los palangreros capturan barriletes de tres grupos de edad y posiblemente de más; en las muestras estudiadas de la composición de las tallas de barrilete, el grupo modal más pequeño (65 cm), observado en enero-marzo en el Pacífico noroeste (160 0E-180° y 20 oN-45°N), corresponde en talla al grupo modal más grande que aparece en la pesca de superficie a fines del verano frente a las Islas Izu-Bonín al sudeste del Japón, y se compara también con la talla modal del barrilete obtenido en la pesca hawaiana en la época de primavera; el análisis de las capturas de barrilete por medio del estudio de la posición de los anzuelos en el palangre y por la tasa de mortalidad, indica que parte de la captura se efectúa cuando el equipo está en movimiento cerca a la superficie y una parte inferior de la captura se realiza, cuando las artes se estabilizan a una profundidad de 70 a 140 m. Se ofrece una breve discusión sobre varias hipótesis de otros autores, en vista de la nueva información presentada referente a la estructura poblacional y a la migración del barrilete en el Océano Pacífico. (PDF contains 100 pages.)

Effects of volumetric allocation on heave response of semisubmersible in deep sea

The configuration of semisubmersibles consisting of pontoons and columns and their corresponding heave motion response in incident progressive waves are examined. The purpose of the present study is to provide a theoretical approach to estimating the effects of volumetric allocation on natural period and response amplitude operator (RAO) in heave motion. We conclude that the amplitude of heave motion response can be considerably suppressed by appropriately adjusting volumetric allocation so that the natural heave period keeps away from the range of wave energy. The theoretical formulae are found in good agreement with the corresponding computational results by WAMIT.

The hydro-kinetic theory of liquid helium II

In Part I the kinetic theory of excitations in flowing liquid He II is developed to a higher order than that carried out previously, by Landau and Khalatnikov, in order to demonstrate the existence of non-equilibrium terms of a new nature in the hydrodynamic equations. It is then shown that these terms can lead to spontaneous destabilization in counter currents when the relative velocity of the normal and super fluids exceeds a critical value that depends on the temperature, but not on geometry. There are no adjustable parameters in the theory. The critical velocities are estimated to be in the 14-20 m/sec range for T ≤ 2.0° K, but tend to zero as T → T_λ. The possibility that these critical velocities may be related to the experimentally observed "intrinsic" critical velocities is discussed.

Part II consists of a semi-classical investigation of rotonquantized vortex line interactions. An essentially classical model is used for the collision and the behavior of the roton in the vortex field is investigated in detail. From this model it is possible to derive the HVBK mutual friction terms that appear in the phenomenalogical equations of motion for rotating liquid He II. Estimates of the Hall and Vinen B and B' coefficients are in good agreement with experiments. The claim is made that the theory does not contain any arbitrary adjustable parameters.

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

Design for all: designing for the motion-impaired user

Designing for all requires the adaptation and modification of current design best practices to encompass a broader range of user capabilities. This is particularly the case in the design of the human-product interface. Product interfaces exist everywhere and when designing them, there is a very strong temptation to jump to prescribing a solution with only a cursory attempt to understand the nature of the problem. This is particularly the case when attempting to adapt existing designs, optimised for able-bodied users, for use by disabled users. However, such approaches have led to numerous products that are neither usable nor commercially successful. In order to develop a successful design approach it is necessary consider the fundamental structure of the design process being applied. A three stage design process development strategy which includes problem definition, solution development and solution evaluation, should be adopted. This paper describes the development of a new design approach based on the application of usability heuristics to the design of interfaces. This is illustrated by reference to a particular case study of the re-design of a computer interface for controlling an assistive device.

Three-dimensional nonlinear inelastic analysis of steel moment-frame buildings damaged by earthquake excitations

The Northridge earthquake of January 17, 1994, highlighted the two previously known problems of premature fracturing of connections and the damaging capabilities of near-source ground motion pulses. Large ground motions had not been experienced in a city with tall steel moment-frame buildings before. Some steel buildings exhibited fracture of welded connections or other types of structural degradation.

A sophisticated three-dimensional nonlinear inelastic program is developed that can accurately model many nonlinear properties commonly ignored or approximated in other programs. The program can assess and predict severely inelastic response of steel buildings due to strong ground motions, including collapse.

Three-dimensional fiber and segment discretization of elements is presented in this work. This element and its two-dimensional counterpart are capable of modeling various geometric and material nonlinearities such as moment amplification, spread of plasticity and connection fracture. In addition to introducing a three-dimensional element discretization, this work presents three-dimensional constraints that limit the number of equations required to solve various three-dimensional problems consisting of intersecting planar frames.

Two buildings damaged in the Northridge earthquake are investigated to verify the ability of the program to match the level of response and the extent and location of damage measured. The program is used to predict response of larger near-source ground motions using the properties determined from the matched response.

A third building is studied to assess three-dimensional effects on a realistic irregular building in the inelastic range of response considering earthquake directivity. Damage levels are observed to be significantly affected by directivity and torsional response.

Several strong recorded ground motions clearly exceed code-based levels. Properly designed buildings can have drifts exceeding code specified levels due to these ground motions. The strongest ground motions caused collapse if fracture was included in the model. Near-source ground displacement pulses can cause columns to yield prior to weaker-designed beams. Damage in tall buildings correlates better with peak-to-peak displacements than with peak-to-peak accelerations.

Dynamic response of tall buildings shows that higher mode response can cause more damage than first mode response. Leaking of energy between modes in conjunction with damage can cause torsional behavior that is not anticipated.

Various response parameters are used for all three buildings to determine what correlations can be made for inelastic building response. Damage levels can be dramatically different based on the inelastic model used. Damage does not correlate well with several common response parameters.

Realistic modeling of material properties and structural behavior is of great value for understanding the performance of tall buildings due to earthquake excitations.