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.

Linear stability of miscible two-fluid flow down an incline

We show that miscible two-layer free-surface flows of varying viscosity down an inclined substrate are different in their stability characteristics from both immiscible two-layer flows, and flows with viscosity gradients spanning the entire flow. New instability modes arise when the critical layer of the viscosity transport equation overlaps the viscosity gradient. A lubricating configuration with a less viscous wall layer is identified to be the most stabilizing at moderate miscibility (moderate Peclet numbers). This also is in contrast with the immiscible case, where the lubrication configuration is always destabilizing. The co-existence that we find under certain circumstances, of several growing overlap modes, the usual surface mode, and a Tollmien-Schlichting mode, presents interesting new possibilities for nonlinear breakdown. © 2013 AIP Publishing LLC.

Design and modelling of a fluid inerter

Mechanical spring-damper network performance can often be improved by the inclusion of a third passive component called the inerter. This ideally has the characteristic that the force at the terminals is directly proportional to the relative acceleration between them. The fluid inerter presented here has advantages over mechanical ball screw devices in terms of simplicity of design. Furthermore, it can be readily adapted to implement various passive network layouts. Variable orifices and valves can be included to provide series or parallel damping. Test data from prototypes with helical tubes have been compared with models to investigate parasitic damping effects of the fluid. © 2013 Copyright Taylor and Francis Group, LLC.

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.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

Wetness loss prediction for a low pressure steam turbine using computational fluid dynamics

Two-phase computational fluid dynamics modelling is used to investigate the magnitude of different contributions to the wet steam losses in a three-stage model low pressure steam turbine. The thermodynamic losses (due to irreversible heat transfer across a finite temperature difference) and the kinematic relaxation losses (due to the frictional drag of the drops) are evaluated directly from the computational fluid dynamics simulation using a concept based on entropy production rates. The braking losses (due to the impact of large drops on the rotor) are investigated by a separate numerical prediction. The simulations show that in the present case, the dominant effect is the thermodynamic loss that accounts for over 90% of the wetness losses and that both the thermodynamic and the kinematic relaxation losses depend on the droplet diameter. The numerical results are brought into context with the well-known Baumann correlation, and a comparison with available measurement data in the literature is given. The ability of the numerical approach to predict the main wetness losses is confirmed, which permits the use of computational fluid dynamics for further studies on wetness loss correlations. © IMechE 2013 Reprints and permissions: sagepub.co.uk/journalsPermissions.nav.

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea

Based on the comprehensive interpretation and study of the Neogene fracture system and diapiric structure, it can be concluded that the diapiric structures, high-angle fractures and vertical fissure system are the main gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, northern South China Sea. The Neogene fractures widely developed in the study area may be classed into two groups: NW (NNW)-trending and NE (NNE)-trending. The first group was active in the Late Miocene, while the second one was active since the Pliocene. The NE (NNE)-trending fractures were characterized by lower activity strength and larger scale, and cut through the sediment layers deposited since the Pliocene. Within the top sediment layers, the high-angle fracture and vertical fissure system was developed. The diapiric structures display various types such as a turtle-back-like arch, weak piercing, gas chimney, and fracture (or crack, fissure). On the seismic profile, some diapiric structures show the vertical chimney pathway whose top is narrow and the bottom is wide, where some ones extend horizontally into pocket or flower-shaped structures and formed the seismic reflection chaotic zones. Within the overlying sediment layers of the diapiric structure, the tree branch, flower-shaped high-angle fractures and vertical fissures were developed and became the pathway and migration system of the gas-bearing fluid influx. In the study area, the diapiric structures indicate a high temperature/over pressure system ever developed. Closely associated and abundant bright-spots show the methane-bearing fluid influx migrated vertically or horizontally through the diapiric structures, high-angle fractures and vertical fissures. In the place where the temperature and pressure conditions were favor for the formation of gas hydrate, the hydrate reservoir deposition sub-system was developed.

Investigation of Interfacial Tensions for Carbon Dioxide Aqueous Solutions by Perturbed-Chain Statistical Associating Fluid Theory Combined with Density-Gradient Theory

The perturbed-chain statistical associating fluid theory and density-gradient theory are used to construct an equation of state (EOS) applicable for the phase behaviors of carbon dioxide aqueous solutions. With the molecular parameters and influence parameters respectively regressed from bulk properties and surface tensions of pure fluids as input, both the bulk and interfacial properties of carbon dioxide aqueous solutions are satisfactorily correlated by adjusting the binary interaction parameter (k(ij)). Our results show that the constructed EOS is able to describe the interfacial properties of carbon dioxide aqueous solutions in a wide temperature range, and illustrate the influences of temperature, pressure, and densities in each phase on the interfacial properties.

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.