A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

Time-dependent fracture toughness of cornea.

The fracture and time-dependent properties of cornea are very important for the development of corneal scaffolds and prostheses. However, there has been no systematic study of cornea fracture; time-dependent behavior of cornea has never been investigated in a fracture context. In this work, fracture toughness of cornea was characterized by trouser tear tests, and time-dependent properties of cornea were examined by stress-relaxation and uniaxial tensile tests. Control experiments were performed on a photoelastic rubber sheet. Corneal fracture resistance was found to be strain-rate dependent, with values ranging from 3.39±0.57 to 5.40±0.48kJm(-2) over strain rates from 3 to 300mmmin(-1). Results from stress-relaxation tests confirmed that cornea is a nonlinear viscoelastic material. The cornea behaved closer to a viscous fluid at small strain but became relatively more elastic at larger strain. Although cornea properties are greatly dependent on time, the stress-strain responses of cornea were found to be insensitive to the strain rate when subjected to tensile loading.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

Vibration of pressurised, elastic, spherical shells

Balloons are one example of pressurised, elastic, spherical shells. Whilst analytical solutions exist for the vibration of pressurised spheres, these models only incorporate constant tension in the membrane. For elastic shells, changes in curvature will result in restoring forces that are proportional to the elasticity in the membrane; hence the assumption of constant tension is not valid. This paper describes an analytical solution for the natural frequencies of an elastic spherical shell subject to internal pressure. When the membrane tension is set to zero, the results are shown to converge to the analytical solution for a spherical shell, and when the skin elasticity is neglected, the results converge to the constant-tension solution. This analytical solution is used to predict the natural frequencies of a small balloon, based on a value for the elastic modulus that is determined using biaxial tensile testing. These predictions are compared to experimental measurements of balloon vibrations using impact hammer testing, and good agreement is seen.

Wetness loss prediction for a low pressure steam turbine using computational fluid dynamics

Two-phase computational fluid dynamics modelling is used to investigate the magnitude of different contributions to the wet steam losses in a three-stage model low pressure steam turbine. The thermodynamic losses (due to irreversible heat transfer across a finite temperature difference) and the kinematic relaxation losses (due to the frictional drag of the drops) are evaluated directly from the computational fluid dynamics simulation using a concept based on entropy production rates. The braking losses (due to the impact of large drops on the rotor) are investigated by a separate numerical prediction. The simulations show that in the present case, the dominant effect is the thermodynamic loss that accounts for over 90% of the wetness losses and that both the thermodynamic and the kinematic relaxation losses depend on the droplet diameter. The numerical results are brought into context with the well-known Baumann correlation, and a comparison with available measurement data in the literature is given. The ability of the numerical approach to predict the main wetness losses is confirmed, which permits the use of computational fluid dynamics for further studies on wetness loss correlations. © IMechE 2013 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav.

Hopping robot based on free vibration of an elastic curved beam

This study presents a novel approach to the design of low-cost and energy-efficient hopping robots, which makes use of free vibration of an elastic curved beam. We found that a hopping robot could benefit from an elastic curved beam in many ways such as low manufacturing cost, light body weight and small energy dissipation in mechanical interactions. A challenging problem of this design strategy, however, lies in harnessing the mechanical dynamics of free vibration in the elastic curved beam: because the free vibration is the outcome of coupled mechanical dynamics between actuation and mechanical structures, it is not trivial to systematically design mechanical structures and control architectures for stable locomotion. From this perspective, this paper investigates a case study of simple hopping robot to identify the design principles of mechanics and control. We developed a hopping robot consisting of an elastic curved beam and a small rotating mass, which was then modeled and analyzed in simulation. The experimental results show that the robot is capable of exhibiting stable hopping gait patterns by using a small actuation with no sensory feedback owing to the intrinsic stability of coupled mechanical dynamics. Furthermore, an additional analysis shows that, by exploiting free vibration of the elastic curved beam, cost of transport of the proposed hopping locomotion can be in the same rage of animals' locomotion including human running. © 2011 IEEE.

A lattice dynamical treatment for the total potential energy of single-walled carbon nanotubes and its applications: relaxed equilibrium structure, elastic properties, and vibrational modes of ultra-narrow tubes

In this paper, we propose a lattice dynamic treatment for the total potential energy of single-walled carbon nanotubes (SWCNTs) which is, apart from a parameter for the nonlinear effects, extracted from the vibrational energy of the planar graphene sheet. The energetics, elasticity and lattice dynamics are treated in terms of the same set of force constants, independently of the tube structures. Based upon this proposal, we have investigated systematically the relaxed lattice configuration for narrow SWCNTs, the strain energy, the Young's modulus and Poisson ratio, and the lattice vibrational properties with respect to the relaxed equilibrium tubule structure. Our calculated results for various physical quantities are nicely in consistency with existing experimental measurements. In particular, we verified that the relaxation effect makes the bond length longer and the frequencies of various optical vibrational modes softer. Our calculation provides evidence that the Young's modulus of an armchair tube exceeds that of the planar graphene sheet, and that the large diameter limits of the Young's modulus and Poisson ratio are in agreement with the experimental values of graphite; the calculated radial breathing modes for ultra-narrow tubes with diameters ranging between 2 and 5 angstrom coincide with the experimental results and the existing ab initio calculations with satisfaction. For narrow tubes with a diameter of 20 angstrom, the calculated frequencies of optical modes in the tubule's tangential plane, as well as those of radial breathing modes, are also in good agreement with the experimental measurements. In addition, our calculation shows that various physical quantities of relaxed SWCNTs can actually be expanded in terms of the chiral angle defined for the corresponding ideal SWCNTs.

Structural, elastic, and electronic properties of cubic perovskite BaHfO3 obtained from first principles

We investigated the structural, elastic, and electronic properties of the cubic perovskite-type BaHfO3 using a first-principles method based on the plane-wave basis set. Analysis of the band structure shows that perovskite-type BaHfO3 is a wide gap indirect semiconductor. The band-gap is predicted to be 3.94 eV within the screened exchange local density approximation (sX-LDA). The calculated equilibrium lattice constant of this compound is in good agreement with the available experimental and theoretical data reported in the literatures. The independent elastic constants (C-11, C-12, and C-44), bulk modules B and its pressure derivatives B', compressibility beta, shear modulus G, Young's modulus Y, Poisson's ratio nu, and Lame constants (mu, lambda) are obtained and analyzed in comparison with the available theoretical and experimental data for both the singlecrystalline and polycrystalline BaHfO3. The bonding-charge density calculation make it clear that the covalent bonds exist between the Hf and 0 atoms and the ionic bonds exist between the Ba atoms and HfO3 ionic groups in BaHfO3. (C) 2009 Elsevier B.V. All rights reserved.

Depth dependent elastic strain in ZnO epilayer: combined Rutherford backscattering/channeling and X-ray diffraction

A ZnO layer was grown by metalorganic chemical vapor deposition (MOCVD) on a sapphire (0 0 0 1) substrate. The perpendicular and parallel elastic strain of the ZnO epilayer, e(perpendicular to) = 0.19%, e(parallel to) = -0.29%, respectively, were derived by using the combination of Rutherford backscattering (RBS)/channeling and X-ray diffraction (XRD). The ratio vertical bar e(parallel to)/ e(perpendicular to)vertical bar = 1.5 indicates that ZnO layer is much stiffer in the a-axis direction than in the c-axis direction. By using RBS/C, the depth dependent elastic strain was deduced. The strain is higher at the depth close to the interface and decreases towards the surface. The negative tetragonal distortion was explained by considering the lattice mismatch and thermal mismatch in ZnO thin film. (c) 2004 Elsevier B.V. All rights reserved.

Dependence of elastic strain field on the self-organized ordering of quantum dot superlattices

A systematic investigation of the strain distribution of self-organized, lens-shaped quantum dot in the case of growth direction on (001) substrate was presented. The three-dimensional finite element analysis for an array of dots was used for the strain calculation. The dependence of the strain energy density distribution on the thickness of the capping layer was investigated in detail when the elastic characteristics of the matrix material were anisotropic. It is shown that the elastic anisotropic greatly influences the stress, strain, and strain energy density in the quantum dot structures. The anisotropic ratio of the matrix material and the combination with different thicknesses of the capping layer, may lead to different strain energy density minimum locations on the capping layer surface, which can result in various vertical ordering phenomena for the next layer of quantum dots, i.e. partial alignment, random alignment, and complete alignment.

Chemical composition and elastic strain in AlInGaN quaternary films

High-quality AlInGaN quaternary layers were grown on c-Al2O3 using a thick GaN template. A full width at half maximum of 0.075 degrees from AlInGaN(0004) rocking curve and a minimum yield of 5.6% from Rutherford backscattering/channelling spectrometry (RBS) prove the AlInGaN layer of a comparable crystalline quality with GaN layers. The chemical compositions (both of Al and In contents) of AlInGaN layers are directly obtained from RBS and elastic recoil detection analysis. The lattice parameters both in perpendicular and parallel directions are deduced from X-ray diffraction. The AlInGaN layer is found to process a compressive strain in parallel direction and a tensile strain in perpendicular direction. (c) 2006 Elsevier B.V. All rights reserved.

Elastic deformation of wafer surfaces in bonding

A model has been proposed for describing elastic deformation of wafer surfaces in bonding. The change of the surface shape is studied on the basis of the distribution of the periodic strain field. With the condition of diminishing periodic strain away from the interface, Airy stress function has been found that satisfies the elastic mechanical equilibrium. The result reveals that the wavy interface elastically deforms a spatial wavelength from the interface. (C) 2000 American Institute of Physics. [S0021-8979(00)04219-5].

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.