Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Fluid-structure interaction simulation of an inflatable aerodynamic tension-cone decelerator

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. In this paper, we focus on an inflatable tension cone design that has potential advantages over other geometries. A computational fluid-structure interaction model of a tension cone is employed to investigate the behavior of the inflatable aeroshell at supersonic speeds for conditions matching recent experimental results. A parametric study is carried out to investigate the deflections of the tension cone as a function of inflation pressure of the torus at a Mach of 2.5. Comparison of the behavior of the structure, amplitude of deformations, and determined loads are reported. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

International Workshop on Fluid-Structure Interaction. Theory, Numerics and Applications

... Fluid-Structure Interaction Simulation of Viscoplastic and Fracturing Thin- Shells Subjected to Underwater Shock Loading Ralf Deiterding, Fehmi Cirak, ...

Fluid-structure interaction with mean flow: Over-scattering and unstable resonance growth

It is known theoretically [1-3] that infinitely long fluid loaded plates in mean flow exhibit a range of unusual phenomena in the 'long time' limit. These include convective instability, absolute instability and negative energy waves which are destabilized by dissipation. However, structures are necessarily of finite length and may have discontinuities. Moreover, linear instability waves can only grow over a limited number of cycles before non-linear effects become dominant. We have undertaken an analytical and computational study to investigate the response of finite, discontinuous plates to ascertain if these unusual effects might be realized in practice. Analytically, we take a "wave scattering" [2,4] - as opposed to a "modal superposition" [5] - view of the fluttering plate problem. First, we solve for the scattering coefficients of localized plate discontinuities and identify a range of parameter space, well outside the convective instability regime, where over-scattering or amplified reflection/transmission occurs. These are scattering processes that draw energy from the mean flow into the plate. Next, we use the Wiener-Hopf technique to solve for the scattering coefficients from the leading and trailing edges of a baffled plate. Finally, we construct the response of a finite, baffled plate by a superposition of infinite plate propagating waves continuously scattering off the plate ends and solve for the unstable resonance frequencies and temporal growth rates for long plates. We present a comparison between our computational results and the infinite plate theory. In particular, the resonance response of a moderately sized plate is shown to be in excellent agreement with our long plate analytical predictions. Copyright © 2010 by ASME.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

Fluid-structure interaction of three-dimensional magnetic artificial cilia

A numerical model is developed to analyse the interaction of artificial cilia with the surrounding fluid in a three-dimensional setting in the limit of vanishing fluid inertia forces. The cilia are modelled using finite shell elements and the fluid is modelled using a boundary element approach. The coupling between both models is performed by imposing no-slip boundary conditions on the surface of the cilia. The performance of the model is verified using various reference problems available in the literature. The model is used to simulate the fluid flow due to magnetically actuated artificial cilia. The results show that narrow and closely spaced cilia create the largest flow, that metachronal waves along the width of the cilia create a significant flow in the direction of the cilia width and that the recovery stroke in the case of the out-of-plane actuation of the cilia strongly depends on the cilia width. © 2012 Cambridge University Press.

Multi-objective optimization of a fluid structure interaction benchmarking

The integration and application of a new multi-objective tabu search optimization algorithm for Fluid Structure Interaction (FSI) problems are presented. The aim is to enhance the computational design process for real world applications and to achieve higher performance of the whole system for the four considered objectives. The described system combines the optimizer with a well established FSI solver which is based on the fully implicit, monolithic formuFlation of the problem in the Arbitrary Lagrangian-Eulerian FEM approach. The proposed solver resolves the proposed uid-structure interaction benchmark which describes the self-induced elastic deformation of a beam attached to a cylinder in laminar channel ow. The optimized ow characteristics of the aforementioned geometrical arrangement illustrate the performance of the system in two dimensions. Special emphasis is given to the analysis of the simulation package, which is of high accuracy and is the core of application. The design process identifies the best combination of ow features for optimal system behavior and the most important objectives. In addition, the presented methodology has the potential to run in parallel, which will significantly speed-up the elapsed time. Finite Element Method (FEM), Fluid-Structure Interaction (FSI), Multi-Ojective Tabu search (MOTS2). Copyright © 2013 Tech Science Press.

Fluid-structure interaction simulations of a tension-cone inflatable aerodynamic decelerator

A series of fluid-structure interaction simulations of an aerodynamic tension-cone supersonic decelerator prototype intended for large mass payload deployment in planetary explorations are discussed. The fluid-structure interaction computations combine large deformation analysis of thin shells with large-eddy simulation of compressible turbulent flows using a loosely coupled approach to enable quantification of the dynamics of the vehicle. The simulation results are compared with experiments carried out at the NASA Glenn Research Center. Reasonably good agreement between the simulations and the experiment is observed throughout a deflation cycle. The simulations help to illuminate the details of the dynamic progressive buckling of the tension-cone decelerator that ultimately results in the collapse of the structure as the inflation pressure is decreased. Furthermore, the tension-cone decelerator exhibits a transient oscillatory behavior under impulsive loading that ultimately dies out. The frequency of these oscillations was determined to be related to the acoustic time scale in the compressed subsonic region between the bow shock and the structure. As shown, when the natural frequency of the structure and the frequency of the compressed subsonic region approximately match, the decelerator exhibits relatively large nonaxisymetric oscillations. The observed response appears to be a fluid-structure interaction resonance resulting from an acoustic chamber (pistonlike) mode exciting the structure. Copyright © 2013 by Christopher Porter, R. Mark Rennie, Eric J. Jumper.

A computational fluid-structure interaction study of inflation cycles of a tension cone decelerator

Inflatable aerodynamic decelerators present potential advantages for planetary entry in missions of robotic and human exploration. The design of these structures face many engineering challenges, including complex deformable geometries, anisotropic material response, and coupled shockturbulence interactions. In this paper, we describe a comprehensive computational fluid-structure interaction study of an inflation cycle of a tension cone decelerator in supersonic flow and compare the simulations with earlier published experimental results. The aeroshell design and flow conditions closely match recent experiments conducted at Mach 2.5. The structural model is a 16-sided polygonal tension cone with seams between each segment. The computational model utilizes adaptive mesh refinement, large-eddy simulation, and shell mechanics with self-contact modeling to represent the flow and structure interaction. This study focuses on the dynamics of the structure as the inflation pressure varies gradually, and the behavior of forces experienced by the flexible and rigid (the payload capsule) structures. © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Instabilities in the oscillatory flow of a complex fluid.

The dynamics of a fluid in a vertical tube, subjected to an oscillatory pressure gradient, is studied experimentally for both a Newtonian and a viscoelastic shear-thinning fluid. Particle image velocimetry is used to determine the two-dimensional velocity fields in the vertical plane of the tube axis, in a range of driving amplitudes from 0.8 to 2.5 mm and of driving frequencies from 2.0 to 11.5 Hz. The Newtonian fluid exhibits a laminar flow regime, independent of the axial position, in the whole range of drivings. For the complex fluid, instead, the parallel shear flow regime exhibited at low amplitudes [Torralba, Phys. Rev. E 72, 016308 (2005)] becomes unstable at higher drivings against the formation of symmetric vortices, equally spaced along the tube. At even higher drivings the vortex structure itself becomes unstable, and complex nonsymmetric structures develop. Given that inertial effects remain negligible even at the hardest drivings (Re < 10(-1)), it is the complex rheology of the fluid that is responsible for the instabilities observed. The system studied represents an interesting example of the development of shear-induced instabilities in nonlinear complex fluids in purely parallel shear flow.

Characteristics of blood flow in microchannels and relevant impact on modelling blood behaviour in biochip separators

This paper reports a perspective investigation of computational modelling of blood fluid in microchannel devices as a preparation for future research on fluid-structure interaction (FSI) in biofluid mechanics. The investigation is carried out through two aspects, respectively on physical behaviours of blood flow in microchannels and appropriate methodology for modelling. The physics of blood flow is targeted to the challenges for describing blood flow in microchannels, including rheology of blood fluid, suspension features of red blood cells (RBCs), laminar hydrodynamic influence and effect of surface roughness. The analysis shows that due to the hyperelastic property of RBC and its comparable dimension with microchannels, blood fluid shows complex behaviours of two phase flow. The trajectory and migration of RBCs require accurate description of RBC deformation and interaction with plasma. Following on a discussion of modelling approaches, i.e. Eulerian method and Lagrangian method, the main stream modelling methods for multiphase flow are reviewed and their suitability to blood flow is analysed. It is concluded that the key issue for blood flow modelling is how to describe the suspended blood cells, modelled by Lagrangian method, and couple them with the based flow, modelled by Eulerian method. The multiphase flow methods are thereby classified based on the number of points required for describing a particle, as follows: (i) single-point particle methods, (ii) mutli-point particle methods, (iii) functional particle methods, and (iv) fluid particle methods. While single-point particle methods concentrate on particle dynamic movement, multipoint and functional particle methods can take into account particle mechanics and thus offer more detailed information for individual particles. Fluid particle methods provide good compromise between two phases, but require additional information for particle mechanics. For furthermore detailed description, we suggest to investigate the possibility using two domain coupling method, in which particles and base flow are modelled by two separated solvers. It is expected that this paper could clarify relevant issues in numerical modelling of blood flow in microchannels and induce some considerations for modelling blood flow using multiphase flow methods. © 2012 IEEE.