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.

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.

The interaction of zinc oxide-based dental cements with aqueous solutions of potassium fluoride

The ability of zinc oxide-based dental cements (zinc phosphate and zinc polycarboxylate) to take up fluoride from aqueous solution has been studied. Only zinc phosphate cement was found to take up any measurable fluoride after 5 h exposure to the solutions. The zinc oxide filler of the zinc phosphate also failed to take up fluoride from solution. The key interaction for this uptake was thus shown to involve the phosphate groups of the set cement. However, whether this took the form of phosphate/fluoride exchange, or the formation of oxyfluoro-phosphate groups was not clear. Fluoride uptake followed radicaltime kinetics for about 2 h in some cases, but was generally better modelled by the Elovich equation, dq(t)/dt = alpha exp(-beta q(t)). Values for alpha varied from 3.80 to 2.48 x 10(4), and for beta from 7.19 x 10(-3) to 0.1946, though only beta showed any sort of trend, becoming smaller with increasing fluoride concentration. Fluoride was released from the zinc phosphate cements in processes that were diffusion based up to M(t)/M(infinity) of about 0.4. No further release occurred when specimens were placed in fresh volumes of deionised water. Only a fraction of the fluoride taken up was re-released, demonstrating that most of the fluoride taken up becomes irreversibly bound within the cement.