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 generic time and space accurate numerical approach to closely coupled fluid-structure interaction problems

Fluid structure interaction, as applied to flexible structures, has wide application in diverse areas such as flutter in aircraft, flow in elastic pipes and blood vessels and extrusion of metals through dies. However a comprehensive computational model of these multi-physics phenomena is a considerable challenge. Until recently work in this area focused on one phenomenon and represented the behaviour of the other more simply even to the extent in metal forming, for example, that the deformation of the die is totally ignored. More recently, strategies for solving the full coupling between the fluid and soild mechanics behaviour have developed. Conventionally, the computational modelling of fluid structure interaction is problematical since computational fluid dynamics (CFD) is solved using finite volume (FV) methods and computational structural mechanics (CSM) is based entirely on finite element (FE) methods. In the past the concurrent, but rather disparate, development paths for the finite element and finite volume methods have resulted in numerical software tools for CFD and CSM that are different in almost every respect. Hence, progress is frustrated in modelling the emerging multi-physics problem of fluid structure interaction in a consistent manner. Unless the fluid-structure coupling is either one way, very weak or both, transferring and filtering data from one mesh and solution procedure to another may lead to significant problems in computational convergence. Using a novel three phase technique the full interaction between the fluid and the dynamic structural response are represented. The procedure is demonstrated on some challenging applications in complex three dimensional geometries involving aircraft flutter, metal forming and blood flow in arteries.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Computational modeling of mold filling and related free-surface flows in shape casting: an overview of the challenges involved

Accurate representation of the coupled effects between turbulent fluid flow with a free surface, heat transfer, solidification, and mold deformation has been shown to be necessary for the realistic prediction of several defects in castings and also for determining the final crystalline structure. A core component of the computational modeling of casting processes involves mold filling, which is the most computationally intensive aspect of casting simulation at the continuum level. Considering the complex geometries involved in shape casting, the evolution of the free surface, gas entrapment, and the entrainment of oxide layers into the casting make this a very challenging task in every respect. Despite well over 30 years of effort in developing algorithms, this is by no means a closed subject. In this article, we will review the full range of computational methods used, from unstructured finite-element (FE) and finite-volume (FV) methods through fully structured and block-structured approaches utilizing the cut-cell family of techniques to capture the geometric complexity inherent in shape casting. This discussion will include the challenges of generating rapid solutions on high-performance parallel cluster technology and how mold filling links in with the full spectrum of physics involved in shape casting. Finally, some indications as to novel techniques emerging now that can address genuinely arbitrarily complex geometries are briefly outlined and their advantages and disadvantages are discussed.

Modelling the tilt-casting process for the tranquil filling of titanium alloy turbine blades

The tilt-casting method is used to achieve tranquil filling of gamma-TiAl turbine blades. The reactive alloy is melted in a cold crucible using an induction coil and then the complete crucible-mould- running system assembly is rotated through 180degrees to transfer the metal into the mould. The induction current is ramped down gradually as the rotation starts and the mould is preheated to maintain superheat. The liquid metal then enters the mould and the gas within it (argon) escapes through the inlet aperture and through auxiliary vents. Solidification starts as soon the metal enters the mould and it is important to account for this effect to predict and prevent misruns. The rotation rate has to be controlled carefully to allow sufficient time for gas evacuation, but at the same time preserve superheat. This 3-phase system is modelled using the FV method, with a fast implicit numerical scheme used to capture the transient liquid free surface. The enthalpy method is used to model solidification and predict defects such as trapped bubbles, macro-porosity or surface connected porosity. Modeling is used to support an experimental program for the development of a production method for gamma-TiAl blades, with a target length of 40cm. The experiments provide validation for the model and the model in turn optimizes the tilt-casting process. The work is part of the EU project IMPRESS.

Reduction in photosystem II efficiency during a virus-controlled Emiliania huxleyi bloom

During viral infection of Emiliania huxleyi, laboratory studies have shown that photo-system (PS) II efficiency declines during the days post-infection and is thought to be associated with viral-induced interruption of electron transport rates between photosystems. However,measuring the impact of viral infection on PSII function in E. huxleyi populations from natural,taxonomically diverse phytoplankton communities is difficult, and whether this phenomenon occurs in nature is presently unknown. Here, chlorophyll fluorescence analysis was used to assess changes in PSII efficiency throughout an E. huxleyi bloom during a mesocosm experiment off the coast of Norway. Specifically, we aimed to determine whether a measurable suppression of the efficiency of PSII photochemistry could be observed due to viral infection of the natural E. huxleyi populations. During the major infection period prior to bloom collapse, there was a significant reduction in PSII efficiency with an average decrease in maximum PSII photochemical efficiency (Fv/Fm) of 17% and a corresponding 75% increase in maximum PSII effective absorption cross section(σPSII); this was concurrent with a significant decrease in E. huxleyi growth rates and an increase in E. huxleyi virus (EhV) production. As E. huxleyi populations dominated the phytoplankton community and potentially contributed up to 100% of the chlorophyll a pool, we believe that the variable chlorophyll fluorescence signal measured during this period was derived predominantly from E. huxleyi and, thus, reflects changes occurring within E. huxleyi cells. This is the first demonstration of suppression of PSII photochemistry occurring during viral infection of natural coccolithophore populations.

Contribution of the UV component of natural sunlight to photoinhibition of photosynthesis in six species of subtidal brown and red seaweeds.

Field-collected specimens of three species of Laminaria and three species of subtidal red algae (Delesseria sanguinea, Plocamium cartilagineum and Phyllophora pseudoceranoides) were exposed to natural summer sunlight on Helgoland (southern North Sea) for up to 4 h at 15 °C. Dark-adapted variable fluorescence (Fv : Fm) was measured immediately after these treatments, and following 6, 24 and 48 h of recovery in moderate irradiances of white light. The response of plants to the full spectrum of natural sunlight was compared with that to PAR alone, UV-A + visible, UV-A + UV-B, or UV-A alone. The Fv : Fm values of all species were reduced to minimal values after 4 h in all of these treatments, but those of the more resistant species (Laminaria spp. and P. pseudoceranoides) were higher after shorter exposures to UV radiation alone than to PAR with or without UV. The recovery of Fv : Fm in all species was also more rapid in the two treatments that contained UV radiation alone than in those that included PAR. These results suggest that it is the high irradiances of PAR in natural sunlight which are responsible for the photoinhibition of photosynthesis of subtidal seaweeds and that the current ambient irradiances of UV radiation (either UV-B or UV-A) in northern temperate latitudes would not contribute significantly to this photoinhibition.

The effect of increased dietary fruit and vegetable consumption on endothelial activation, inflammation and oxidative stress in hypertensive volunteers

Background and aims

Public health campaigns recommend increased fruit and vegetable (FV) consumption as an effective means of cardiovascular risk reduction. During an 8 week randomised control trial among hypertensive volunteers, we noted significant improvements in endothelium-dependent vasodilatation with increasing FV consumption. Circulating indices of inflammation, endothelial activation and insulin resistance are often employed as alternative surrogates for systemic arterial health. The responses of several such biomarkers to our previously described FV intervention are reported here.
Methods and results

Hypertensive volunteers were recruited from medical outpatient clinics. After a common 4 week run-in period during which FV consumption was limited to 1 portion per day, participants were randomised to 1, 3 or 6 portions daily for 8 weeks. Venous blood samples for biomarker analyses were collected during the pre and post-intervention vascular assessments. A total of 117 volunteers completed the 12 week study. Intervention-related changes in circulating levels of high sensitivity C-reactive protein (hsCRP), soluble intracellular adhesion molecule-1 (sICAM-1), soluble vascular cell adhesion molecule-1 (sVCAM-1), von Willebrand factor (vWF) and plasminogen activator inhibitor-1 (PAI-1) did not differ significantly between FV groups. Similarly, there were no significant between group differences of change in homeostasis model assessment (HOMA) scores.
Conclusions

Despite mediating a significant improvement in acetylcholine induced vasodilatation, increased FV consumption did not affect a calculated measure of insulin resistance or concentrations of the circulating biomarkers measured during this study. Functional indices of arterial health such as endothelium-dependent vasomotion are likely to provide more informative cardiovascular end-points during short-term dietary intervention trials.

Effect of fruit and vegetable intake on oxidative stress and inflammation in COPD: a randomised controlled trial

Epidemiological evidence supports a positive relationship between fruit and vegetable (FV) intake, lung function and chronic obstructive pulmonary disease (COPD). Increasing FV intake may attenuate the oxidative stress and inflammation associated with COPD.

An exploratory randomised controlled trial to examine the effect of increased consumption of FV on oxidative stress and inflammation in moderate-to-severe COPD was conducted. 81 symptomatically stable patients with a habitually low FV intake (two or fewer portions of FV per day) were randomised to the intervention group (five or more portions of FV per day) or the control group (two or fewer portions of FV per day). Each participant received self-selected weekly home deliveries of FV for 12 weeks.

75 participants completed the intervention. There was a significant between-group change in self-reported FV intake and biomarkers of FV intake (zeaxanthin (p=0.034) and ß-cryptoxanthin (p=0.015)), indicating good compliance; post-intervention intakes in intervention and control groups were 6.1 and 1.9 portions of FV per day, respectively. There were no significant changes in biomarkers of airway inflammation (interleukin-8 and myeloperoxidase) and systemic inflammation (C-reactive protein) or airway and systemic oxidative stress (8-isoprostane).

This exploratory study demonstrated that patients with moderate-to-severe COPD were able to comply with an intervention to increase FV intake; however, this had no significant effect on airway or systemic oxidative stress and inflammation.