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.

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.

Ca2+- and volume-sensitive chloride currents are differentially regulated by agonists and store-operated Ca2+ entry.

Using patch-clamp and calcium imaging techniques, we characterized the effects of ATP and histamine on human keratinocytes. In the HaCaT cell line, both receptor agonists induced a transient elevation of [Ca2+]i in a Ca2+-free medium followed by a secondary [Ca2+]i rise upon Ca2+ readmission due to store-operated calcium entry (SOCE). In voltage-clamped cells, agonists activated two kinetically distinct currents, which showed differing voltage dependences and were identified as Ca2+-activated (ICl(Ca)) and volume-regulated (ICl, swell) chloride currents. NPPB and DIDS more efficiently inhibited ICl(Ca) and ICl, swell, respectively. Cell swelling caused by hypotonic solution invariably activated ICl, swell while regulatory volume decrease occurred in intact cells, as was found in flow cytometry experiments. The PLC inhibitor U-73122 blocked both agonist- and cell swelling–induced ICl, swell, while its inactive analogue U-73343 had no effect. ICl(Ca) could be activated by cytoplasmic calcium increase due to thapsigargin (TG)-induced SOCE as well as by buffering [Ca2+]i in the pipette solution at 500 nM. In contrast, ICl, swell could be directly activated by 1-oleoyl-2-acetyl-sn-glycerol (OAG), a cell-permeable DAG analogue, but neither by InsP3 infusion nor by the cytoplasmic calcium increase. PKC also had no role in its regulation. Agonists, OAG, and cell swelling induced ICl, swell in a nonadditive manner, suggesting their convergence on a common pathway. ICl, swell and ICl(Ca) showed only a limited overlap (i.e., simultaneous activation), although various maneuvers were able to induce these currents sequentially in the same cell. TG-induced SOCE strongly potentiated ICl(Ca), but abolished ICl, swell, thereby providing a clue for this paradox. Thus, we have established for the first time using a keratinocyte model that ICl, swell can be physiologically activated under isotonic conditions by receptors coupled to the phosphoinositide pathway. These results also suggest a novel function for SOCE, which can operate as a "selection" switch between closely localized channels.