Christoph Brücker

Entraining in trout: a behavioural and hydrodynamic analysis.

SUMMARY Rheophilic fish commonly experience unsteady flows and hydrodynamic perturbations. Instead of avoiding turbulent zones though, rheophilic fish often seek out these zones for station holding. A behaviour associated with station holding in running water is called entraining. We investigated the entraining behaviour of rainbow trout swimming in the wake of a D-shaped cylinder or sideways of a semi-infinite flat plate displaying a rounded leading edge. Entraining trout moved into specific positions close to and sideways of the submerged objects, where they often maintained their position without corrective body and/or fin motions. To identify the hydrodynamic mechanism of entraining, the flow characteristics around an artificial trout placed at the position preferred by entraining trout were analysed. Numerical simulations of the 3-D unsteady flow field were performed to obtain the unsteady pressure forces. Our results suggest that entraining trout minimise their energy expenditure during station holding by tilting their body into the mean flow direction at an angle, where the resulting lift force and wake suction force cancel out the drag. Small motions of the caudal and/or pectoral fins provide an efficient way to correct the angle, such that an equilibrium is even reached in case of unsteadiness imposed by the wake of an object.

Measuring Flow Velocity and Flow Direction by Spatial and Temporal Analysis of Flow Fluctuations

If exposed to bulk water flow, fish lateral line afferents respond only to flow fluctuations (AC) and not to the steady (DC) component of the flow. Consequently, a single lateral line afferent can encode neither bulk flow direction nor velocity. It is possible, however, for a fish to obtain bulk flow information using multiple afferents that respond only to flow fluctuations. We show by means of particle image velocimetry that, if a flow contains fluctuations, these fluctuations propagate with the flow. A cross-correlation of water motion measured at an upstream point with that at a downstream point can then provide information about flow velocity and flow direction. In this study, we recorded from pairs of primary lateral line afferents while a fish was exposed to either bulk water flow, or to the water motion caused by a moving object. We confirm that lateral line afferents responded to the flow fluctuations and not to the DC component of the flow, and that responses of many fiber pairs were highly correlated, if they were time-shifted to correct for gross flow velocity and gross flow direction. To prove that a cross-correlation mechanism can be used to retrieve the information about gross flow velocity and direction, we measured the flow-induced bending motions of two flexible micropillars separated in a downstream direction. A cross-correlation of the bending motions of these micropillars did indeed produce an accurate estimate of the velocity vector along the direction of the micropillars.

Three-dimensional nature of the glottal jet

The factors contributing to human voice production are not yet fully understood. Even normal human phonation with a symmetric glottal opening area is still the subject of extensive investigation. Among others, it has already been shown that fluid dynamics has a strong influence on the vocal process. The full characterization of the glottal jet has not been accomplished yet. Time-resolved measurement and visualization of the three-dimensional (3D) flow downstream the human vocal folds are difficult if not impossible to perform in vivo. Therefore, it is common to use mechanical and numerical models with a simplified shape and motion profile of the vocal folds. In this article, further results regarding the 3D flow structure obtained in a 3:1 up-scaled dynamic glottis model (cam model) in a water circuit are given, extending earlier work [M. Triep et al. (2005). Exp. Fluids 39, 232-245]. The model mimics the temporal variation in the 3D contour of the glottal gap while water flow reduces the characteristic frequencies by the order of 1/140. The unsteady flow processes downstream of the vocal folds are visualized in slow motion and analyzed in detail via particle imaging techniques. The visualization results show complex 3D flow behavior of lengthwise jet contraction and axis switching. In addition, the time-dependent flow rate during the phonatory oscillation cycle is measured in detail. It is shown that the pressure loss is decreased in the presence of a second constriction downstream of the glottis in form of ventricular folds and it is observed that for this case the jet is stabilized in the divergent phase of the cycle.

Mathematical Models and Numerical Schemes for the Simulation of Human Phonation

Acoustic data has long been harvested in fundamental voice investigations since it is easily obtained using a microphone. However, acoustic signals alone do not reveal much about the complex interplay between sound waves, structural surface waves, mechanical vibrations, and fluid flow involved in phonation. Available high speed imaging techniques have over the past ten years provided a wealth of information about the mechanical deformation of the superior surface of the larynx during phonation. Time-resolved images of the inner structure of the deformable soft tissues are not yet feasible because of low temporal resolution (MRI and ultrasound) and x-ray dose-related hazards (CT and standard x- ray). One possible approach to circumvent these challenges is to use mathematical models that reproduce observable behavior such as phonation frequency, closed quotient, onset pressure, jitter, shimmer, radiated sound pressure, and airflow. Mathematical models of phonation range in complexity from systems with relatively small degrees of freedom (multi-mass models) to models based on partial differential equations (PDEs) mostly solved by finite element (FE) methods resulting in millions of degrees-of-freedom. We will provide an overview about the current state of mathematical models for the human phonation process, since they have served as valuable tools for providing insight into the basic mechanisms of phonation and may eventually be of sufficient detail and accuracy to allow surgical planning, diagnostics, and rehabilitation evaluations on an individual basis. Furthermore, we will also critically discuss these models w.r.t. the used geometry, boundary conditions, material properties, their verification, and reproducibility.

In vitro experimental investigation of voice production.

The process of human phonation involves a complex interaction between the physical domains of structural dynamics, fluid flow, and acoustic sound production and radiation. Given the high degree of nonlinearity of these processes, even small anatomical or physiological disturbances can significantly affect the voice signal. In the worst cases, patients can lose their voice and hence the normal mode of speech communication. To improve medical therapies and surgical techniques it is very important to understand better the physics of the human phonation process. Due to the limited experimental access to the human larynx, alternative strategies, including artificial vocal folds, have been developed. The following review gives an overview of experimental investigations of artificial vocal folds within the last 30 years. The models are sorted into three groups: static models, externally driven models, and self-oscillating models. The focus is on the different models of the human vocal folds and on the ways in which they have been applied.

Single-view volumetric PIV via high-resolution scanning, isotropic voxel restructuring and 3D least-squares matching (3D-LSM)

Scanning PIV as introduced by Br?cker (1995 Exp. Fluids 19 255?63, 1996a Appl. Sci. Res. 56 157?79) has been successfully applied in the last 20?years to different flow problems where the frame rate was sufficient to ensure a ?frozen? field condition. The limited number of parallel planes however leads typically to an under-sampling in the scan direction in depth; therefore, the spatial resolution in depth is typically considerably lower than the spatial resolution in the plane of the laser sheet (depth resolution = scan shift ?z ? pixel unit in object space). In addition, a partial volume averaging effect due to the thickness of the light sheet must be taken into account. Herein, the method is further developed using a high-resolution scanning in combination with a Gaussian regression technique to achieve an isotropic representation of the tracer particles in a voxel-based volume reconstruction with cuboidal voxels. This eliminates the partial volume averaging effect due to light sheet thickness and leads to comparable spatial resolution of the particle field reconstructions in x-, y- and z-axes. In addition, advantage of voxel-based processing with estimations of translation, rotation and shear/strain is taken by using a 3D least-squares matching method, well suited for reconstruction of grey-level pattern fields. The method is discussed in this paper and used to investigate the ring vortex instability at Re = 2500 within a measurement volume of roughly 75???75???50?mm3?with a spatial resolution of 100??m/voxel (750???750???500 voxel elements). The volume has been scanned with a number of 100 light sheets and scan rates of 10?kHz. The results show the growth of the Tsai?Widnall azimuthal instabilities accompanied with a precession of the axis of the vortex ring. Prior to breakdown, secondary instabilities evolve along the core with streamwise oriented striations. The front stagnation points streamwise distance to the core starts to decrease while the rear stagnation point distance remains constant which indicates that the front part of the ring is at first losing its mass during breakdown.

Vortex dynamics in the wake of a mechanical fish

This study focuses on the three-dimensional flow around a mechanical fish model, which reproduces the typical undulatory body and fin motion of a carangiform swimmer. The mechanical model consists of a flexible skeleton embedded in a soft transparent silicone body, which is connected with two cams to a flapping and bending hinge generating a traveling wave motion with increasing amplitude from anterior to posterior, extending to a combined heaving and pitching motion at the fin. The model is submerged in a water tank and towed at the characteristic swimming speed for the neutral swimming mode at U/V = 1. The method of Scanning Particle Image Velocimetry was used to analyze the three-dimensional time-dependent flow field in the axial and saggital planes. The results confirm the earlier observations that the wake develops into a chain of vortex rings which travel sidewards perpendicular to the swimming direction. However, instead of one single vortex shed at each tail beat half-cycle we observed a pair of two vortex rings being shed. Each pair consists of a larger main vortex ring corresponding to the tail beat start-stop vortex, while the second vortex ring is due to the body bending motion. The existence of the second vortex reflects the role of the body in undulatory swimming. A simplified model of the fish body comparing it to a plate with a hinged flap demonstrates the link between the sequence of kinematics and vortex shedding.

Asymmetric glottal jet deflection: Differences of two- and three-dimensional models

Flow is studied through a channel with an oscillating orifice mimicking the motion of the glottal-gap during phonation. Simulations with prescribed flow and wall-motion are carried out for different orifice geometries, a 2D slit-like and a 3D lens-like one. Although the jet emerges from a symmetric orifice a significant deflection occurs in case of the slit-like geometry, contrary to the 3D lens-like one. The results demonstrate the dependency of jet entrainment and vortex dynamics on the orifice geometry and the interpretation of asymmetric jet deflection with regard to the relevance of the Coanda effect in the process of human phonation.

Diving-Flight Aerodynamics of a Peregrine Falcon (Falco peregrinus)

This study investigates the aerodynamics of the falcon Falco peregrinus while diving. During a dive peregrines can reach velocities of more than 320 km h−1. Unfortunately, in freely roaming falcons, these high velocities prohibit a precise determination of flight parameters such as velocity and acceleration as well as body shape and wing contour. Therefore, individual F. peregrinus were trained to dive in front of a vertical dam with a height of 60 m. The presence of a well-defined background allowed us to reconstruct the flight path and the body shape of the falcon during certain flight phases. Flight trajectories were obtained with a stereo high-speed camera system. In addition, body images of the falcon were taken from two perspectives with a high-resolution digital camera. The dam allowed us to match the high-resolution images obtained from the digital camera with the corresponding images taken with the high-speed cameras. Using these data we built a life-size model of F. peregrinus and used it to measure the drag and lift forces in a wind-tunnel. We compared these forces acting on the model with the data obtained from the 3-D flight path trajectory of the diving F. peregrinus. Visualizations of the flow in the wind-tunnel uncovered details of the flow structure around the falcon’s body, which suggests local regions with separation of flow. High-resolution pictures of the diving peregrine indicate that feathers pop-up in the equivalent regions, where flow separation in the model falcon occurred.

Multiple-plane particle image velocimetry using a light-field camera.

Planar velocity fields in flows are determined simultaneously on parallel measurement planes by means of an in-house manufactured light-field camera. The planes are defined by illuminating light sheets with constant spacing. Particle positions are reconstructed from a single 2D recording taken by a CMOS-camera equipped with a high-quality doublet lens array. The fast refocusing algorithm is based on synthetic-aperture particle image velocimetry (SAPIV). The reconstruction quality is tested via ray-tracing of synthetically generated particle fields. The introduced single-camera SAPIV is applied to a convective flow within a measurement volume of 30 x 30 x 50 mm³.