Euler–Euler large eddy simulations for dispersed turbulent bubbly flows

In this paper we present detailed Euler–Euler Large Eddy Simulations (LES) of dispersed bubbly flow in a rectangular bubble column. The motivation of this study is to investigate the potential of this approach for the prediction of bubbly flows, in terms of mean quantities. The physical models describing the momentum exchange between the phases including drag, lift and wall force were chosen according to previous experiences of the authors. Experimental data, Euler–Lagrange LES and unsteady Euler–Euler Reynolds-Averaged Navier–Stokes results are used for comparison. It is found that the present model combination provides good agreement with experimental data for the mean flow and liquid velocity fluctuations. The energy spectrum obtained from the resolved velocity of the Euler–Euler LES is presented as well.     

On a model for the pedestrians-induced lateral vibrations of footbridges

We study the synchronization of the pedestrians motion with the lateral motion of slender footbridges, which attracted the attention of researchers mainly after the Millennium Bridge well known event. We consider a bridge-pedestrians model developed by Strogatz et al. To better understand some aspects of the underlying mechanical phenomena and to increase the agreement of the model results to the effective behaviour of the walkers (as observed during the Millennium Bridge’s opening day and during Arup’s tests on site) we have introduced some modifications to the model and we have performed extensive parametric investigations, supported by many numerical simulations. Our approach is therefore computational, by means of a self-made code. This permits to highlight the parameters which mainly affect the trigger and the development of the phenomenon of synchronous lateral excitation, thus allowing a good understanding of the physical event and an evaluation of the engineering reliability of the Strogatz et al. model.     

On the Existence and Uniqueness of a Solution to the Interior Transmission Problem for Piecewise-Homogeneous Solids

The interior transmission problem (ITP), which plays a fundamental role in inverse scattering theories involving penetrable defects, is investigated within the framework of mechanical waves scattered by piecewise-homogeneous, elastic or viscoelastic obstacles in a likewise heterogeneous background solid. For generality, the obstacle is allowed to be multiply connected, having both penetrable components (inclusions) and impenetrable parts (cavities). A variational formulation is employed to establish sufficient conditions for the existence and uniqueness of a solution to the ITP, provided that the excitation frequency does not belong to (at most) countable spectrum of transmission eigenvalues. The featured sufficient conditions, expressed in terms of the mass density and elasticity parameters of the problem, represent an advancement over earlier works on the subject in that (i) they pose a precise, previously unavailable provision for the well-posedness of the ITP in situations when both the obstacle and the background solid are heterogeneous, and (ii) they are dimensionally consistent, i.e., invariant under the choice of physical units. For the case of a viscoelastic scatterer in an elastic solid it is further shown, consistent with earlier studies in acoustics, electromagnetism, and elasticity that the uniqueness of a solution to the ITP is maintained irrespective of the vibration frequency. When applied to the situation where both the scatterer and the background medium are viscoelastic, i.e., dissipative, on the other hand, the same type of analysis shows that the analogous claim of uniqueness does not hold. Physically, such anomalous behavior of the “viscoelastic-viscoelastic” case (that has eluded previous studies) has its origins in a lesser known fact that the homogeneous ITP is not mechanically insulated from its surroundings—a feature that is particularly cloaked in situations when either the background medium or the scatterer are dissipative. A set of numerical results, computed for ITP configurations that meet the sufficient conditions for the existence of a solution, is included to illustrate the problem. Consistent with the preceding analysis, the results indicate that the set of transmission values is indeed empty in the “elastic-viscoelastic” case, and countable for “elastic-elastic” and “viscoelastic-viscoelastic” configurations.     

On the Asymptotic Behaviour of the Solutions of the Equations of Linear Elastodynamics in Unbounded Domains

We prove, among other things, that if the acoustic tensor satisfies a suitable growth condition at infinity (the hyperbolicity condition) and the total initial energy is summable with a suitable weight, then the solution to the initial boundary value problem of linear elastodynamics in unbounded domains decays at infinity, at every instant, with a rate depending on the weight. Moreover, we show that the hyperbolicity condition is necessary and sufficient for the equipartition in mean of the total energy.     

Lattice–Boltzmann simulations for analysing the detachment of micron-sized spherical particles from surfaces with large-scale roughness structures

Fully resolved numerical simulations of a micron-sized spherical particle residing on a surface with large-scale roughness are performed by using the Lattice–Boltzmann method. The aim is to investigate the influence of surface roughness on the detachment of fine drug particles from larger carrier particles for transporting fine drug particles in a DPI (dry powder inhaler). Often the carrier surface is modified by mechanical treatments for modifying the surface roughness in order to reduce the adhesion force of drug particles. Therefore, drug particle removal from the carrier surface is equivalent to the detachment of a sphere from a rough plane surface. Here a sphere with a diameter of 5 μm at a particle Reynolds number of 1.0, 3.5 and 10 are considered. The surface roughness is described as regularly spaced semi-cylindrical asperities (with the axes oriented normal to the flow direction) on a smooth surface. The influence of asperity distance and size ratio (i.e. the radius of the semi-cylinder to the particle radius, R_c/R_d) on particle adhesion and detachment are studied. The asperity distance is varied in the range 1.2 < L/R_d < 2 and the semi-cylinder radius between 0.5 < R_c/R_d < 0.75. The required particle resolution and domain size are appropriately selected based on numerical studies, and a parametric analysis is performed to investigate the relationship between the contact distance (i.e. half the distance between the particle contact points on two neighbouring semi-cylinders), the asperity distance, the size ratio, and the height of the particle centroid from the plane wall. The drag, lift and torque acting on the spherical particle are measured for different particle Reynolds numbers, asperity distances and sizes or diameters. The detachment of particles from rough surfaces can occur through lift-off, sliding and rolling, and the corresponding detachment models are constructed for the case of rough surfaces. These studies will be the basis for developing Lagrangian detachment models that eventually should allow the optimisation of dry powder inhaler performance through computational fluid dynamics.     

Phenomenological continuous contact–impact modelling for multibody simulations of pedestrian–vehicle contact interactions based on experimental data

Multibody modelling of pedestrian collisions requires the definition of contact–impact between the pedestrian and the vehicle. An examination of relevant impact test data reveals large rate-dependent components of the reaction force, permanent indentation, and concomitant energy loss. Contact–impact models previously used in simulations of pedestrian impacts typically have not adequately modelled one, two or all three of these phenomena. This paper presents a phenomenological contact–impact model based on the Hunt–Crossley model of impact, which includes rate-dependent damping, and is extended to include permanent indentation. The proposed model suitably characterises impact test data in a form that can also be implemented in the multibody simulation code MADYMO (TASS-Safe, Netherlands). The proposed contact–impact model was used to characterise the impact between a legform and the bumper of a vehicle, based on two impact tests conducted at different impact speeds. A single contact–impact definition in MADYMO closely reproduced the dynamics of both tests. The proposed model may be suitable in a wide range of impact conditions where the impact is modelled using multibody techniques and it is practicable to conduct impact tests as part of the modelling process.     

Effect of large elastic strains on cavitation instability predictions for elastic–plastic solids

For an infinite solid containing a void, the cavitation instability limit is defined as the remote stress–and strain state, at which the void grows without bound, driven by the elastic energy stored in the surrounding material. Such cavitation limits have been analysed by a number of authors for metal plasticity as well as for nonlinear elastic solids. The analyses for elastic–plastic solids are here extended to consider the effect of a large initial yield strain, and it is shown how the critical stress value decays for increasing value of the yield strain. Analyses are carried out for remote hydrostatic tension as well as for more general axisymmetric remote stress field, with an initially spherical void. Different levels of strain hardening are considered.     

The development of high-speed particle image velocimetry (20 kHz) for large eddy simulation code refinement in bluff body flows

Flow interaction with a bluff body generates a highly complex flow field and has been the subject of much experimental and theoretical analysis. It has been shown that large eddy simulation (LES) modelling provides a more realistic analysis of the flow for such situations where the large scales of turbulence must be resolved. The inherent small-scale spatial velocity averaging in particle image velocimetry (PIV) is commensurate with the sub-grid scale modelling of LES and, therefore, offers potential as a code refinement technique. To demonstrate this potential, however, PIV must be performed with a temporal resolution of typically kHz and a spatial resolution of sub-mm² to be relevant for the vast majority of flows of engineering interest. This paper reports the development of a high-speed PIV system capable of operating at 20 kHz with a spatial resolution of 0.9 mm². This is the combined highest speed, highest resolution PIV data reported to date. The experiment chosen to demonstrate the system is the study of the steady flow interaction with circular and square cross-section obstacles. A Reynolds number of 3,900 is chosen for the cylinder flow to extend the database used by Breuer M. (1998 Int J Heat Fluid 19:512–521) in his extensive LES modelling of this flow. Data presented include a sequence of two-dimensional velocity and vorticity fields, including flow streamlines. Importantly, the random error, inherent in a PIV measurement, is discussed and a formula presented which allows the error to be estimated and regions of the flow identified where LES comparisons would be uncertain.     

On the Formulation of Continuum Thermodynamic Models for Solids as General Equations for Non-equilibrium Reversible-Irreversible Coupling

The purpose of this work is the comparison of some aspects of the formulation of material models in the context of continuum thermodynamics (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997) with their formulation in the form of a General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC: e.g., Grmela and Öttinger in Phys. Rev. E 56: 6620–6632, 1997; Öttinger and Grmela in Phys. Rev. E 56: 6633–6655, 1997; Öttinger in Beyond equilibrium thermodynamics, Wiley, New York, 2005; Grmela in J. Non-Newton. Fluid Mech. 165: 980–998, 2010). A GENERIC represents a generalization of the Ginzburg-Landau model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to formulate non-equilibrium thermodynamic models for complex fluids, it has recently been applied to anisotropic inelastic solids in a Eulerian setting (Hütter and Tervoort in J. Non-Newton. Fluid Mech. 152: 45–52, 2008; 53–65, 2008; Adv. Appl. Mech. 42: 254–317, 2009) as well as to damage mechanics (Hütter and Tervoort in Acta Mech. 201: 297–312, 2008). In the current work, attention is focused for simplicity on the case of thermoelastic solids with heat conduction and viscosity in a Lagrangian setting (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997, Chaps. 9–12). In the process, the relation of the two formulations to each other is investigated in detail. A particular point in this regard is the concept of dissipation and its model representation in both contexts.     

On the Principal Boundary Value Problem in the Two-Dimensional Anisotropic Theory of Elasticity

By definition, the principal problem of the two-dimensional theory of elasticity consists in solving the equation for the Airy’s stress function in a region with its first order derivatives assigned at a boundary. In this paper, an indirect formulation of this problem based on integral equations with weakly singular kernels is proposed. In a bounded region with a Lyapunov boundary it is reduced to the solution of weakly singular integral equations. Differential properties of its solution are investigated.     

The influences of both offset and mass moment of inertia of a tip mass on the dynamics of a centrifugally stiffened visco-elastic beam

The present study is concerned with the out-of-plane vibrations of a rotating, internally damped (Kelvin-Voigt model) Bernoulli-Euler beam carrying a tip mass. The centroid of the tip mass, possessing also a mass moment of inertia is offset from the free end of the beam and is located along its extended axis. This system can be thought of as an extremely simplified model of a helicopter rotor blade or a blade of an auto-cooling fan. The differential eigenvalue problem is solved by using Frobenius method of solution in power series. The characteristic equation is then solved numerically. The simulation results are tabulated for a variety of the nondimensional rotational speeds, tip mass, tip mass offset, mass moment of inertia and internal damping parameters. These are compared with the results of a conventional finite element modeling as well, and excellent agreement is obtained. Some numerical results are given in graphical form. The numerical results obtained, indicate clearly that the tip mass offset and mass moment of inertia are important parameters on the eigencharacteristics of rotating beams so that they have to be included in the modeling process.     

On the local description of wear-induced volume loss and shape changes for engineering surfaces

The paper examines the possibility to obtain a device-independent representation of the effects induced by wear, compatible with the expectable, future development of a similitude approach to the analysis and comparison of tribological performances. Geometrical coefficients presented in former studies are redefined according to a generalised, local formulation, referred to nominal conditions and not constrained within a particular theoretical model, scale level or measurement technique. Their knowledge (from computation or experiment) allows in principle to evaluate the instantaneous distribution of wear (assessment and prediction of shape changes) and the actual kinematics (deviations from the nominal set-up relating to sliding displacement and clearance increase). First results are presented, from numerical simulations. They are consistent with the available literature; moreover, they suggest that, mainly due to compatibility restrictions, some basic geometrical features of the worn surfaces are fairly invariant with test configuration and wear law, making possible their approximated modelling from a relatively reduced set of data.     

On the determination of the force distribution in overconstrained cable-driven parallel mechanisms

This paper addresses the determination of the force distribution in the cables of a redundantly actuated cable-driven parallel mechanism. First, the static model of cable-driven parallel mechanisms is derived based on the wrench matrix. Then, four performance indices are considered in order to solve the underdetermined problem associated with the distribution of the forces. A simple numerical example is then developed in order to provide insight into the problem, which leads to a geometric interpretation of the results. Based on the presented results, it is proposed to use a p-norm (e.g. a 4-norm) to optimize the distribution of the forces in a cable-driven parallel mechanism in order to minimize the largest deviations from the median forces (or other target values) while maintaining continuity in the solution. A non-iterative polynomial formulation is then proposed for the 4-norm. It is also pointed out that this formulation leads to a unique real solution.     

On Spatial Behavior of the Harmonic Vibrations in Kelvin-Voigt Materials

The present paper deals with the study of the amplitude of the steady-state vibrations in a right finite cylinder made of an isotropic Kelvin-Voigt material. Some exponential decay estimates, similar to those of Saint-Venant type, are obtained for appropriate cross-sectional area measures associated with the amplitude of the steady-state vibrations. It is proved that due to dissipative effects, the estimates in question hold for every value of the frequency of vibrations and for arbitrary values of the elastic coefficients. The results are extended to a semi-infinite cylinder and some alternatives of Phragmèn-Lindelöf type are established.     

An Exact Derivation from Three-Dimensions of the Linear Equations for the Bending of an Elastically Monoclinic Plate and Associated Three-Dimensional Solutions

The exact linear three-dimensional equations for a elastically monoclinic (13 constant) plate of constant thickness are reduced without approximation to a single 4th order differential equation for a thickness-weighted normal displacement plus two auxiliary equations for weighted thickness integrals of a stress function and the normal strain. The 4th order equation is of the same form as in classical (Kirchhoff) theory except the unknown is not the midsurface normal displacement. Assuming a solution of these plate equations, we construct so-called modified Saint-Venant solutions—“modified” because they involve non-zero body and surface loads. That is, solutions of the exact three-dimensional elasticity equations that exhibit no boundary layers and that are subject to a special set of body and surface loads that leave the analogous plate loads arbitrary.     

Parameterization and modelling of large off-road tyres for ride analyses: Part 2 – Parameterization and validation of tyre models

Every mathematical model used in a simulation is an idealization and simplification of reality. Vehicle dynamic simulations that go beyond the fundamental investigations require complex multi-body simulation models. The tyre–road interaction presents one of the biggest challenges in creating an accurate vehicle model. Many tyre models have been proposed and developed but proper validation studies are less accessible. These models were mostly developed and validated for passenger car tyres for application on relatively smooth roads. The improvement of ride comfort, safety and structural integrity of large off-road vehicles, over rough terrain, has become more significant in the development process of heavy vehicles. This paper investigates whether existing tyre models can be used to accurately describe the vertical behaviour of large off road tyres while driving over uneven terrain. [1] Presented an extensive set of experimentally determined parameterization and validation data for a large off-road tyre. Both laboratory and field test are performed for various loads, inflation pressures and terrain inputs. The parameterization process of four tyre models or contact models are discussed in detail. The parameterized models are then validated against test results on various hard but rough off-road terrain and the results are discussed.     

Application of digital image correlation to the study of planar anisotropy of sheet metals at large strains

The aim of this work is to measure and model the planar anisotropy of thin steel sheets. The experimental data have been collected using the digital image correlation technique. This is a powerful tool to measure the strain field on differently shaped specimens subjected to large plastic deformations. In this manner, it is possible to observe the material behaviour under different stress-strain states, from small to large deformation conditions, on the entire specimen surface. The experimental results on smooth and notched samples have been used to characterize the flow stress curve after necking and a nonassociated plastic flow rule is proposed to describe the anisotropic behaviour of the material. To compare the experimental data with the predictions of the adopted constitutive model, a novel method, based on the image correlation results, has been implemented.     

On the averaging principle for one-frequency systems. An application to satellite motions

This paper is related to our previous works (Morosi and Pizzocchero in J. Phys. A, Math. Gen. 39:3673–3702, 2006; Nonlinear Dyn., 2008), on the error estimate of the averaging technique for systems with one fast angular variable. In the cited references, a general method (of mixed analytical and numerical type) has been introduced to obtain precise, fully quantitative estimates on the averaging error. Here, this procedure is applied to the motion of a satellite in a polar orbit around an oblate planet, retaining only the J ₂ term in the multipole expansion of the gravitational potential. To exemplify the method, the averaging errors are estimated for the data corresponding to two Earth satellites; for a very large number of orbits, computation of our estimators is much less expensive than the direct numerical solution of the equations of motion.     

A symbolic algorithm for the automatic computation of multitone-input harmonic balance equations for nonlinear systems

A symbolic algorithm is developed for the automatic generation of harmonic balance equations for multitone input for a class of nonlinear differential systems with polynomial nonlinearities. Generalized expressions are derived for the construction of balance equations for a defined multitone signal form. Procedures are described for determining combinations for a given output frequency from the desired set obtained from box truncated spectra and their permutations to automate symbolic algorithm. An application of method is demonstrated using the well-known Duffing–Van der Pol equation. Then the obtained analytical results are compared with numerical simulations to show the accuracy of the approach. The computation times for both the numerical solutions of equations versus the number of frequency components and the symbolic generation of the equations versus the power of nonlinearity are also investigated.