Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Application of domain decomposition and transmission line modelling techniques to 2D, time-domain, finite element problems

In this paper a recently published finite element method, which combines domain decomposition with a novel technique for solving nonlinear magnetostatic finite element problems is described. It is then shown how the method can be extended to, and optimised for, the solution of time-domain problems. © 1999 IEEE.

Active control of high-frequency vibration: Optimisation using the hybrid modelling method

This work presents active control of high-frequency vibration using skyhook dampers. The choice of the damper gain and its optimal location is crucial for the effective implementation of active vibration control. In vibration control, certain sensor/actuator locations are preferable for reducing structural vibration while using minimum control effort. In order to perform optimisation on a general built-up structure to control vibration, it is necessary to have a good modelling technique to predict the performance of the controller. The present work exploits the hybrid modelling approach, which combines the finite element method (FEM) and statistical energy analysis (SEA) to provide efficient response predictions at medium to high frequencies. The hybrid method is implemented here for a general network of plates, coupled via springs, to allow study of a variety of generic control design problems. By combining the hybrid method with numerical optimisation using a genetic algorithm, optimal skyhook damper gains and locations are obtained. The optimal controller gain and location found from the hybrid method are compared with results from a deterministic modelling method. Good agreement between the results is observed, whereas results from the hybrid method are found in a significantly reduced amount of time. © 2012 Elsevier Ltd. All rights reserved.

Theoretical simulation studies of pulsed field magnetisation of (RE)BCO bulk superconductors

The magnetisation of bulk high temperature superconductors (HTS), such as RE-Ba-Cu-O [(RE)BCO, where RE is a rare earth element or Y], by a practical technique is essential for their application in high field, permanent magnet-like devices. Research to-date into the pulsed field magnetisation (PFM) of these materials, however, has been limited generally to experimental techniques, with relatively little progress in the development of theoretical models. This is because not only is a multi-physics approach needed to take account of the heating of the samples but also the high electric fields generated are well above the regime in which there are reliable experimental results. This paper describes a framework of theoretical simulation using the finite element method (FEM) that is applicable to both single- and multi-pulse magnetisation processes of (RE)BCO bulk superconductors. The model incorporates the heat equation and provides a convenient way of determining the distribution of trapped field, current density and temperature change within a bulk superconductor at each stage of the magnetisation process. An example of the single-pulse magnetisation of a (RE)BCO bulk is described. Potentially, the model may serve as a cost-effective tool for the optimisation of the bulk geometry and the magnetisation profile in multi-pulse magnetisation processes. © 2010 IOP Publishing Ltd.

Elastodynamic erosion of thermal barrier coatings

The finite element method is used to analyze the elastodynamic response of a columnar thermal barrier coating due to normal impact and oblique impact by an erosive particle. An assessment is made of the erosion by crack growth from preexisting flaws at the edge of each column: it is demonstrated that particle impacts can be sufficiently severe to give rise to columnar cracking. First, the transient stress state induced by the normal impact of a circular cylinder or a sphere is calculated in order to assess whether a 2D calculation adequately captures the more realistic 3D behavior. It is found that the transient stress states for the plane strain and axisymmetric models are similar. The sensitivity of response to particle diameter and to impact velocity is determined for both the cylinder and the sphere. Second, the transient stress state is explored for 2D oblique impact by a circular cylindrical particle and by an angular cylindrical particle. The sensitivity of transient tensile stress within the columns to particle shape (circular and angular), impact angle, impact location, orientation of the angular particle, and to the level of friction is explored in turn. The paper concludes with an evaluation of the effect of inclining the thermal barrier coating columns upon their erosion resistance. © 2011 The American Ceramic Society.

Hydrodynamic interaction of two unsteady model microorganisms.

The study of pair-wise interactions between swimming microorganisms is fundamental to the understanding of the rheological and transport properties of semi-dilute suspensions. In this paper, the hydrodynamic interaction of two ciliated microorganisms is investigated numerically using a boundary-element method, and the microorganisms are modeled as spherical squirmers that swim by time-dependent surface deformations. The results show that the inclusion of the unsteady terms in the ciliary propulsion model has a large impact on the trajectories of the interacting cells, and causes a significant change in scattering angles with potential important consequences on the diffusion properties of semi-dilute suspensions. Furthermore, the analysis of the shear stress acting on the surface of the microorganisms revealed that the duration and the intensity of the near-field interaction are significantly modified by the presence of unsteadiness. This observation may account for the hydrodynamic nature of randomness in some biological reactions, and supersedes the distinction between intrinsic randomness and hydrodynamic interactions, adding a further element to the understanding and modeling of interacting microorganisms.

The direct field boundary impedance of two-dimensional periodic structures with application to high frequency vibration prediction.

Large sections of many types of engineering construction can be considered to constitute a two-dimensional periodic structure, with examples ranging from an orthogonally stiffened shell to a honeycomb sandwich panel. In this paper, a method is presented for computing the boundary (or edge) impedance of a semi-infinite two-dimensional periodic structure, a quantity which is referred to as the direct field boundary impedance matrix. This terminology arises from the fact that none of the waves generated at the boundary (the direct field) are reflected back to the boundary in a semi-infinite system. The direct field impedance matrix can be used to calculate elastic wave transmission coefficients, and also to calculate the coupling loss factors (CLFs), which are required by the statistical energy analysis (SEA) approach to predicting high frequency vibration levels in built-up systems. The calculation of the relevant CLFs enables a two-dimensional periodic region of a structure to be modeled very efficiently as a single subsystem within SEA, and also within related methods, such as a recently developed hybrid approach, which couples the finite element method with SEA. The analysis is illustrated by various numerical examples involving stiffened plate structures.

Hydrodynamic entrapment of bacteria swimming near a solid surface.

The near-surface motility of bacteria is important in the initial formation of biofilms and in many biomedical applications. The swimming motion of Escherichia coli near a solid surface is investigated both numerically and experimentally. A boundary element method is used to predict the hydrodynamic entrapment of E. coli bacteria, their trajectories, and the minimum separation of the cell from the surface. The numerical results show the existence of a stable swimming distance from the boundary that depends only on the shape of the cell body and the flagellum. The experimental validation of the numerical approach allows one to use the numerical method as a predictive tool to estimate with reasonable accuracy the near-wall motility of swimming bacteria of known geometry. The analysis of the numerical database demonstrated the existence of a correlation between the radius of curvature of the near-wall circular trajectory and the separation gap. Such correlation allows an indirect estimation of either of the two quantities by a direct measure of the other without prior knowledge of the cell geometry. This result may prove extremely important in those biomedical and technical applications in which the near-wall behavior of bacteria is of fundamental importance.

Modelling of equipment loaded panels for robust control investigations

This paper develops a modelling technique for equipment load panels which directly produces (adequate) models of the underlying dynamics on which to base robust controller design/evaluations. This technique is based on the use of the Lagrange's equations of motion and the resulting models are verified against those produced by a finite Element Method Model.

A cohesive approach to thin-shell fracture and fragmentation

We develop a finite-element method for the simulation of dynamic fracture and fragmentation of thin-shells. The shell is spatially discretized with subdivision shell elements and the fracture along the element edges is modeled with a cohesive law. In order to follow the propagation and branching of cracks, subdivision shell elements are pre-fractured ab initio and the crack opening is constrained prior to crack nucleation. This approach allows for shell fracture in an in-plane tearing mode, a shearing mode, or a bending of hinge mode. The good performance of the method is demonstrated through the simulation of petalling failure experiments in aluminum plates. © 2005 Elsevier B.V. All rights reserved.

The elastic-plastic indentation response of a columnar thermal barrier coating

Thermal barrier coatings with a columnar microstructure are prone to erosion damage by a mechanism of surface cracking upon impact by small foreign particles. In order to explore this erosion mechanism, the elastic indentation and the elastic-plastic indentation responses of a columnar thermal barrier coating to a spherical indenter were determined by the finite element method and by analytical models. It was shown that the indentation response is intermediate between that of a homogeneous half-space and that given by an elastic-plastic mattress model (with the columns behaving as independent non-linear springs). The sensitivity of the indentation behaviour to geometry and to the material parameters was explored: the diameter of the columns, the gap width between columns, the coefficient of Coulomb friction between columns and the layer height of the thermal barrier coating. The calculations revealed that the level of induced tensile stress is sufficient to lead to cracking of the columns at a depth of about the column radius. It was also demonstrated that the underlying soft bond coat can undergo plastic indentation when the coating comprises parallel columns, but this is less likely for the more realistic case of a random arrangement of tapered columns. © 2009 Elsevier B.V.