Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Nucleated polymerization with secondary pathways. I. Time evolution of the principal moments.

Self-assembly processes resulting in linear structures are often observed in molecular biology, and include the formation of functional filaments such as actin and tubulin, as well as generally dysfunctional ones such as amyloid aggregates. Although the basic kinetic equations describing these phenomena are well-established, it has proved to be challenging, due to their non-linear nature, to derive solutions to these equations except for special cases. The availability of general analytical solutions provides a route for determining the rates of molecular level processes from the analysis of macroscopic experimental measurements of the growth kinetics, in addition to the phenomenological parameters, such as lag times and maximal growth rates that are already obtainable from standard fitting procedures. We describe here an analytical approach based on fixed-point analysis, which provides self-consistent solutions for the growth of filamentous structures that can, in addition to elongation, undergo internal fracturing and monomer-dependent nucleation as mechanisms for generating new free ends acting as growth sites. Our results generalise the analytical expression for sigmoidal growth kinetics from the Oosawa theory for nucleated polymerisation to the case of fragmenting filaments. We determine the corresponding growth laws in closed form and derive from first principles a number of relationships which have been empirically established for the kinetics of the self-assembly of amyloid fibrils.

Toward a human-like biped robot with compliant legs

Conventional models of bipedal walking generally assume rigid body structures, while elastic material properties seem to play an essential role in nature. On the basis of a novel theoretical model of bipedal walking, this paper investigates a model of biped robot which makes use of minimum control and elastic passive joints inspired from the structures of biological systems. The model is evaluated in simulation and a physical robotic platform by analyzing the kinematics and ground reaction force. The experimental results show that, with a proper leg design of passive dynamics and elasticity, an attractor state of human-like walking gait patterns can be achieved through extremely simple control without sensory feedback. The detailed analysis also explains how the dynamic human-like gait can contribute to adaptive biped walking. © 2007 Elsevier B.V. All rights reserved.

Toward a human-like biped robot with compliant legs

Conventional models of bipedal walking generally assume rigid body structures, while elastic material properties seem to play an essential role in nature. On the basis of a novel theoretical model of bipedal walking, this paper investigates a model of biped robot which makes use of minimum control and elastic passive joints inspired from the structures of biological systems. The model is evaluated in simulation and a physical robotic platform with respect to the kinematics and the ground reaction force. The experimental results show that the behavior of this simple locomotion model shows a considerable similarity to that of human walking. © 2006 The authors.

Generation and scattering of sound in vortical flow

We present solutions to scattering problems for unsteady disturbances to a mean swirling flow in an annular duct with a rigid 'splitter'. This situation has application to rotor-stator interaction noise in aeroengines, where the flow downstream of the fan is swirling and bifurcates into the by-pass duct and the engine core. We also consider the trailing edge extension of this problem. Inviscid mean flow in a cylindrical annulus is considered, with both axial and swirling (azimuthal) velocity components. The presence of vorticity in the mean flow couples the acoustic and vorticity modes of irrotational flow. Instead we have one combined spectrum of acoustic-vorticity waves in which the 'sonic' and 'nearly-convected' modes are fully coupled. In addition to the aeroacoustics application the results offer insight into the behaviour of these acoustic-vorticity waves, and the precise nature of the coupling between the two types of mode. Two regimes are discussed in which progress has been made, one for a specialised mean flow, uniform axial flow and rigid body swirl, and a second regime in which the frequency is assumed large, valid for any axisymmetric mean flow. The Wiener-Hopf technique is used to solve the scattering problems mathematically, and we present numerical evaluations of these solutions. Several new effects are seen to arise due to the mean vorticity, in particular the generation of sound at a trailing edge due to the scattering of a nearly convected disturbance, in contrast to the way a convected gust silently passes a trailing edge in uniform mean flow.

Modelling of long-term ground response to tunnelling under St James's Park, London

Following a tunnel excavation in low-permeability soil, it is commonly observed that the ground surface continues to settle and ground loading on the tunnel lining changes, as the pore pressures in the ground approach a new equilibrium condition. The monitored ground response following the tunnelling under St James's Park, London, shows that the mechanism of subsurface deformation is composed of three different zones: swelling, consolidation and rigid body movement. The swelling took place in a confined zone above the tunnel crown, extending vertically to approximately 5 m above it. On the sides of the tunnel, the consolidation of the soil occurred in the zone primarily within the tunnel horizon, from the shoulder to just beneath the invert, and extending laterally to a large offset from the tunnel centreline. Above these swelling and consolidation zones the soil moved downward as a rigid body. In this study, soil-fluid coupled three-dimensional finite element analyses were performed to simulate the mechanism of long-term ground response monitored at St James's Park. An advanced critical state soil model, which can simulate the behaviour of London Clay in both drained and undrained conditions, was adopted for the analyses. The analysis results are discussed and compared with the field monitoring data. It is found that the observed mechanism of long-term subsurface ground and tunnel lining response at St James's Park can be simulated accurately only when stiffness anisotropy, the variation of permeability between different units within the London Clay and non-uniform drainage conditions for the tunnel lining are considered. This has important implications for future prediction of the long-term behaviour of tunnels in clays.

The hemodynamic effects of in-tandem carotid artery stenosis: implications for carotid endarterectomy.

OBJECTIVES: It remains controversial whether patients with severe disease of the internal carotid artery and a coexisting stenotic lesion downstream would benefit from a carotid endarterectomy (CEA) of the proximal lesion. The aim of this study was to simulate the hemodynamic and wall shear effects of in-tandem internal carotid artery stenosis using a computational fluid dynamic (CFD) idealized model to give insight into the possible consequences of CEA on these lesions. METHODS: A CFD model of steady viscous flow in a rigid tube with two asymmetric stenoses was introduced to simulate blood flow in arteries with multiple constrictions. The effect of varying the distance between the two stenoses, and the severity of the upstream stenosis on the pressure and wall shear stress (WSS) distributions on the second plaque, was investigated. The influence of the relative positions of the two stenoses was also assessed. RESULTS: The distance between the plaques was found to have minimal influence on the overall hemodynamic effect except for the presence of a zone of low WSS (range -20 to 30 dyne/cm2) adjacent to both lesions when the two stenoses were sufficiently close (<4 times the arterial diameter). The upstream stenosis was protective if it was larger than the downstream stenosis. The relative positions of the stenoses were found to influence the WSS but not the pressure distribution. CONCLUSIONS: The geometry and positions of the lesions need to be considered when considering the hemodynamic effects of an in-tandem stenosis. Low WSS is thought to cause endothelial dysfunction and initiate atheroma formation. The fact that there was a flow recirculation zone with low WSS in between the two stenoses may demonstrate how two closely positioned plaques may merge into one larger lesion. Decision making for CEA may need to take into account the hemodynamic situation when an in-tandem stenosis is found. CFD may aid in the risk stratification of patients with this problem.

MEAN FLOW EFFECTS OF THE LOW-WAVENUMBER PRESSURE SPECTRUM ON A FLEXIBLE SURFACE.

The Lighthill theory is extended so that it may be used to determine the flow noise induced by a turbulent boundary layer over a plane homogeneous flexible surface. The influence of the surface properties and the mean flow on the sound generation is brought out explicitly through the use of a Green function. The form of the low-wavenumber wall-pressure spectrum on a rigid surface with an arbitrary mean flow profile is determined. The effect of a coating layer is investigated.

The modified Kanerva model for automatic speech recognition

A parallel processing network derived from Kanerva's associative memory theory Kanerva 1984 is shown to be able to train rapidly on connected speech data and recognize further speech data with a label error rate of 0·68%. This modified Kanerva model can be trained substantially faster than other networks with comparable pattern discrimination properties. Kanerva presented his theory of a self-propagating search in 1984, and showed theoretically that large-scale versions of his model would have powerful pattern matching properties. This paper describes how the design for the modified Kanerva model is derived from Kanerva's original theory. Several designs are tested to discover which form may be implemented fastest while still maintaining versatile recognition performance. A method is developed to deal with the time varying nature of the speech signal by recognizing static patterns together with a fixed quantity of contextual information. In order to recognize speech features in different contexts it is necessary for a network to be able to model disjoint pattern classes. This type of modelling cannot be performed by a single layer of links. Network research was once held back by the inability of single-layer networks to solve this sort of problem, and the lack of a training algorithm for multi-layer networks. Rumelhart, Hinton & Williams 1985 provided one solution by demonstrating the "back propagation" training algorithm for multi-layer networks. A second alternative is used in the modified Kanerva model. A non-linear fixed transformation maps the pattern space into a space of higher dimensionality in which the speech features are linearly separable. A single-layer network may then be used to perform the recognition. The advantage of this solution over the other using multi-layer networks lies in the greater power and speed of the single-layer network training algorithm. © 1989.

Seismic response of stone masonry spires: Analytical modeling

Stone masonry spires are vulnerable to seismic loading. Computational methods are often used to predict the dynamic linear elastic response of masonry towers and spires, but this approach is only applicable until the first masonry joint begins to open, limiting the ability to predict collapse. In this paper, analytical modeling is used to investigate the uplift, rocking and collapse of stone spires. General equations for static equilibrium of the spire under lateral acceleration are first presented, and provide a reasonable lower bound for predicting collapse. The dynamic response is then considered through elastic modal analysis and rigid body rocking. Together, these methods are used to provide uplift curves and single impulse overturning collapse curves for a complete range of possible spire geometries. Results are used to evaluate the historic collapse of two specific stone spires. © 2012 Elsevier Ltd.

Density stratified environments: The double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the 'double-tank' method proposed by Oster (Sci Am 213:70-76, 1965). We refer to Oster's method as the 'forced-drain' approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the 'free-drain' approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the 'free-drain' extension of Oster's method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established. © 2008 Springer-Verlag.

An experimental study on engine dynamics model based in-cylinder pressure estimation

The information provided by the in-cylinder pressure signal is of great importance for modern engine management systems. The obtained information is implemented to improve the control and diagnostics of the combustion process in order to meet the stringent emission regulations and to improve vehicle reliability and drivability. The work presented in this paper covers the experimental study and proposes a comprehensive and practical solution for the estimation of the in-cylinder pressure from the crankshaft speed fluctuation. Also, the paper emphasizes the feasibility and practicality aspects of the estimation techniques, for the real-time online application. In this study an engine dynamics model based estimation method is proposed. A discrete-time transformed form of a rigid-body crankshaft dynamics model is constructed based on the kinetic energy theorem, as the basis expression for total torque estimation. The major difficulties, including load torque estimation and separation of pressure profile from adjacent-firing cylinders, are addressed in this work and solutions to each problem are given respectively. The experimental results conducted on a multi-cylinder diesel engine have shown that the proposed method successfully estimate a more accurate cylinder pressure over a wider range of crankshaft angles. Copyright © 2012 SAE International.

Guidance of collective cell migration by substrate geometry.

Collective behavior refers to the emergence of complex migration patterns over scales larger than those of the individual elements constituting a system. It plays a pivotal role in biological systems in regulating various processes such as gastrulation, morphogenesis and tissue organization. Here, by combining experimental approaches and numerical modeling, we explore the role of cell density ('crowding'), strength of intercellular adhesion ('cohesion') and boundary conditions imposed by extracellular matrix (ECM) proteins ('constraints') in regulating the emergence of collective behavior within epithelial cell sheets. Our results show that the geometrical confinement of cells into well-defined circles induces a persistent, coordinated and synchronized rotation of cells that depends on cell density. The speed of such rotating large-scale movements slows down as the density increases. Furthermore, such collective rotation behavior depends on the size of the micropatterned circles: we observe a rotating motion of the overall cell population in the same direction for sizes of up to 200 μm. The rotating cells move as a solid body, with a uniform angular velocity. Interestingly, this upper limit leads to length scales that are similar to the natural correlation length observed for unconfined epithelial cell sheets. This behavior is strongly altered in cells that present a downregulation of adherens junctions and in cancerous cell types. We anticipate that our system provides a simple and easy approach to investigate collective cell behavior in a well-controlled and systematic manner.

On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - o realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 AIAA.