A finite volume approach to geometrically non-linear stress analysis

In this past decade finite volume (FV) methods have increasingly been used for the solution of solid mechanics problems. This contribution describes a cell vertex finite volume discretisation approach to the solution of geometrically nonlinear (GNL) problems. These problems, which may well have linear material properties, are subject to large deformation. This requires a distinct formulation, which is described in this paper together with the solution strategy for GNL problem. The competitive performance for this procedure against the conventional finite element (FE) formulation is illustrated for a three dimensional axially loaded column.

A finite volume unstructured mesh approach to dynamic fluid-structure interaction: an assessment of the challenge of flutter analysis

Computational modelling of dynamic fluid-structure interaction (DFSI) is problematical since conventionally computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. Hence, progress in modelling the emerging multi-physics problem of dynamic fluid-structure interaction in a consistent manner is frustrated and significant problems in computation convergence may be encountered in transferring and filtering data from one mesh and solution procedure to another, unless the fluid-structure coupling is either one way, very weak or both. This paper sets out the solution procedure for modelling the multi-physics dynamic fluid-structure interaction problem within a single software framework PHYSICA, using finite volume, unstructured mesh (FV-UM) procedures and will focus upon some of the problems and issues that have to be resolved for time accurate closely coupled dynamic fluid-structure flutter analysis.

Dynamic fluid-structure interactions using finite volume unstructured mesh procedures

A three-dimensional finite volume, unstructured mesh (FV-UM) method for dynamic fluid–structure interaction (DFSI) is described. Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, wind response of buildings, flows in elastic pipes and blood vessels. It involves the coupling of fluid flow and structural mechanics, two fields that are conventionally modelled using two dissimilar methods, thus a single comprehensive computational model of both phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply. More recently, strategies for solving the full coupling between the fluid and solid mechanics behaviour have been developed. A key contribution has been made by Farhat et al. [Int. J. Numer. Meth. Fluids 21 (1995) 807] employing FV-UM methods for solving the Euler flow equations and a conventional finite element method for the elastic solid mechanics and the spring based mesh procedure of Batina [AIAA paper 0115, 1989] for mesh movement. In this paper, we describe an approach which broadly exploits the three field strategy described by Farhat for fluid flow, structural dynamics and mesh movement but, in the context of DFSI, contains a number of novel features: a single mesh covering the entire domain, a Navier–Stokes flow, a single FV-UM discretisation approach for both the flow and solid mechanics procedures, an implicit predictor–corrector version of the Newmark algorithm, a single code embedding the whole strategy.

Dynamic fluid-structure interaction using finite volume unstructured mesh procedures

A three dimensional finite volume, unstructured mesh method for dynamic fluid-structure interation is described. The broad approach is conventional in that the fluid and structure are solved sequentially. The pressure and viscous stresses from the flow algorithm provide load conditions for the solid algorithm, whilst at the fluid structure interface the deformed structure provides boundary condition from the structure to the fluid. The structure algorithm also provides the necessary mesh adaptation for the flow field, the effect of which is accounted for in the flow algorithm. The procedures described in this work have several novel features, namely: * a single mesh covering the entire domain. * a Navier Stokes flow. * a single FV-UM discretisation approach for both the flow and solid mechanics procedures. * an implicit predictor-corrector version of the Newmark algorithm. * a single code embedding the whole strategy. The procedure is illustrated for a three dimensional loaded cantilever in fluid flow.

Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, PHYSICA, for later application to full dynamic fluid structure interaction.

A vertex-based finite volume method applied to non-linear material problems in computational solid mechanics

A vertex-based finite volume (FV) method is presented for the computational solution of quasi-static solid mechanics problems involving material non-linearity and infinitesimal strains. The problems are analysed numerically with fully unstructured meshes that consist of a variety of two- and threedimensional element types. A detailed comparison between the vertex-based FV and the standard Galerkin FE methods is provided with regard to discretization, solution accuracy and computational efficiency. For some problem classes a direct equivalence of the two methods is demonstrated, both theoretically and numerically. However, for other problems some interesting advantages and disadvantages of the FV formulation over the Galerkin FE method are highlighted.

A finite volume unstructured mesh approach to dynamic fluid–structure interaction: an assessment of the challenge of predicting the onset of flutter

Computational modelling of dynamic fluid–structure interaction (DFSI) is a considerable challenge. Our approach to this class of problems involves the use of a single software framework for all the phenomena involved, employing finite volume methods on unstructured meshes in three dimensions. This method enables time and space accurate calculations in a consistent manner. One key application of DFSI simulation is the analysis of the onset of flutter in aircraft wings, where the work of Yates et al. [Measured and Calculated Subsonic and Transonic Flutter Characteristics of a 45° degree Sweptback Wing Planform in Air and Freon-12 in the Langley Transonic Dynamic Tunnel. NASA Technical Note D-1616, 1963] on the AGARD 445.6 wing planform still provides the most comprehensive benchmark data available. This paper presents the results of a significant effort to model the onset of flutter for the AGARD 445.6 wing planform geometry. A series of key issues needs to be addressed for this computational approach. • The advantage of using a single mesh, in order to eliminate numerical problems when applying boundary conditions at the fluid-structure interface, is counteracted by the challenge of generating a suitably high quality mesh in both the fluid and structural domains. • The computational effort for this DFSI procedure, in terms of run time and memory requirements, is very significant. Practical simulations require even finer meshes and shorter time steps, requiring parallel implementation for operation on large, high performance parallel systems. • The consistency and completeness of the AGARD data in the public domain is inadequate for use in the validation of DFSI codes when predicting the onset of flutter.

Combined vertex-based - cell-centred finite volume method for flows in complex geometries

CFD modelling of 'real-life' processes often requires solutions in complex three dimensional geometries, which can often result in meshes where aspects of it are badly distorted. Cell-centred finite volume methods, typical of most commercial CFD tools, are computationally efficient, but can lead to convergence problems on meshes which feature cells with high non-orthogonal shapes. The vertex-based finite volume method handles distorted meshes with relative ease, but is computationally expensive. A combined vertex-based - cell-centred (VB-CC) technique, detailed in this paper, allows solutions on distorted meshes that defeat purely cell-centred physical models to be employed in the solution of other transported quantities. The VB-CC method is validated with benchmark solutions for thermally driven flow and turbulent flow. An early application of this hybrid technique is to three-dimensional flow over an aircraft wing, although it is planned to use it in a wide variety of processing applications in the future.

A cell-centred finite volume method for modelling viscoelastic flow

An unstructured cell-centred finite volume method for modelling viscoelastic flow is presented. The method is applied to the flow through a planar channel and the 4:1 planar contraction for creeping flow of an Oldroyd-B fluid. The results are presented for a range of Weissenberg numbers. In the case of the planar channel results are compared with analytical solutions. For the 4:1 planar contraction benchmark problem the convection terms in the constitutive equations are approximated using both first and second order differencing schemes to compare the techniques and the effect of mesh refinement on the solution is investigated. This is the first time that a fully unstructured, cell-centredfinitevolume technique has been used to model the Oldroyd-B fluid for the test cases presented in this paper.

Muscle activation patterns in shoulder impingement patients

Introduction: Shoulder impingement is one of the most common presentations of shoulder joint problems 1. It appears to be caused by a reduction in the sub-acromial space as the humerus abducts between 60o -120o – the 'painful arc'. Structures between the humeral head and the acromion are thus pinched causing pain and further pathology 2. Shoulder muscle activity can influence this joint space but it is unclear whether this is a cause or effect in impingement patients. This study aimed to observe muscle activation patterns in normal and impingement shoulder patients and determine if there were any significant differences. Method: 19 adult subjects were asked to perform shoulder abduction in their symptomatic arm and non-symptomatic. 10 of these subjects (age 47.9 ± 11.2) were screened for shoulder impingement, and 9 subjects (age 38.9 ± 14.3) had no history of shoulder pathology. Surface EMG was used to collect data for 6 shoulder muscles (Upper, middle and lower trapezius, serratus anterior, infraspinatus, middle deltoids) which was then filtered and fully rectified. Subjects performed 3 smooth unilateral abduction movements at a cadence of 16 beats of a metronome set at 60bpm, and the mean of their results was recorded. T-tests were used to indicate any statistical significance in the data sets. Significance was set at P<0.05. Results: There was a significant difference in muscle activation with serratus anterior in particular showing a very low level of activation throughout the range when compared to normal shoulder activation patterns (<30%). Middle deltoid recruitment was significantly reduced between 60-90o in the impingement group (30:58%).Trends were noted in other muscles with upper trapezius and infraspinatus activating more rapidly and erratically (63:25%; 60:27% respectively), and lower trapezius with less recruitment (13:30%) in the patient group, although these did not quite reach significance. Conclusion: There appears to be some interesting alterations in muscle recruitment patterns in impingement shoulder patients when compared against their own unaffected shoulders and the control group. In particular changes in scapula control (serratus anterior and trapezius) and lateral rotation (infraspinatus), which have direct influence on the sub-acromial space, should be noted. It is still not clear whether these alterations are causative or reactionary, but this finding gives a clear indication to the importance of addressing muscle reeducation as part of a rehabilitation programme in shoulder impingement patients.

The influence of crank configuration on muscle activity and torque production during arm crank ergometry

This study investigated the effect of crank configuration on muscle activity and torque production during submaximal arm crank ergometry. Thirteen non-specifically trained male participants volunteered. During the research trials they completed a warm-up at 15 W before two 3-min exercise stages were completed at 50 and 100 W; subjects used either a synchronous or asynchronous pattern of cranking. During the final 30-s of each submaximal exercise stage electromyographic and torque production data were collected. After the data had been processed each parameter was analysed using separate 2-way ANOVA tests with repeated measures. The activity of all muscles increased in line with external workload, although a shift in the temporal pattern of muscle activity was noted between crank configurations. Patterns of torque production during asynchronous and synchronous cranking were distinct. Furthermore, peak, minimum and delta (peak-minimum) torque values were different (P < 0.05) between crank configurations at both workloads. For example, at 100 W, peak torque using synchronous [19.6 (4.3) Nm] cranking was higher (P < 0.05) compared to asynchronous [16.8 (1.6) Nm] cranking. In contrast minimum torque was lower (P < 0.05) at 100 W using synchronous [4.8 (1.7) Nm] compared to asynchronous [7.3 (1.2) Nm] cranking. There was a distinct bilateral asymmetry in torque production during asynchronous cranking with the dominant transmitting significantly more force to the crank arm. Taken together, these preliminary data demonstrate the complex nature of muscle activity during arm crank ergometry performed with an asynchronous or synchronous crank set-up. Further work is required to determine how muscle activity (EMG activity) and associated patterns of torque production influence physiological responses and functional capacity during arm crank ergometry.

Explicit finite volume modeling of aircraft anti-icing and de-icing

In this paper, a runback water and ice prediction model is extended to anti-icing and thermal de-icing situations. The resulting coupled equations that govern thin-film flow, ice accretion, and heat conduction in the multilayered system substrate-ice-water are solved using an explicit finite volume approach. The procedure is implemented in the three-dimensional icing code ICECREMO2, and both structured and unstructured grids can be considered. Numerical results are presented to compare the present code simulations to some data provided by other ice prediction codes and to show the capabilities of the present numerical tool.