Simulation of Parachute Inflation Dynamics Using an Eulerian Computational Framework for Fluid-Structure Interfaces Evolving in High-Speed Turbulent Flows

Huang, Zhengyu; Avery, Philip; Farhat, Charbel; Rabinovitch, Jason; Derkevorkian, Armen; Peterson, Lee D.

doi:10.2514/6.2018-1540

Cited by 34 publications

(20 citation statements)

References 18 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The numerical simulation of multiphysics problems involving multiple physical models or multiple simultaneous physical phenomena is significant in many engineering and scientific applications, e.g., fluid-structure interactions (FSI) in aeroelasticity [3,4,5] or biomechanics [6,7,8], chemical reaction in combustion or subsurface flows [9,10], electricity and magnetism with hydrodynamics in plasma physics [11,12,13], among others. These problems are generally highly nonlinear, feature multiple scales and strong coupling effects, and require heterogeneous discretizations for the various physics subsystems.…”

Section: Introductionmentioning

confidence: 99%

High-order partitioned spectral deferred correction solvers for multiphysics problems

Huang

Pazner

Persson

et al. 2020

Journal of Computational Physics

View full text Add to dashboard Cite

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1,2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1,2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

show abstract

Section: Introductionmentioning

confidence: 99%

High-order partitioned spectral deferred correction solvers for multiphysics problems

Huang

Pazner

Persson

et al. 2020

Journal of Computational Physics

View full text Add to dashboard Cite

show abstract

“…The underlying fluid‐structure computational model must account for material and geometric porosities in the parachute's soft goods, must be able to handle the geometric and material nonlinearities of the structural subsystem, must resolve the shocks and turbulent wakes in the fluid subsystem, and must capture the effects of all interactions between these structural and flow features on various physical instabilities such as flutter and pulsation. Furthermore, the large displacements, rotations, and deformations and, more importantly, the large topological changes and massive self‐contact incurred during the deployment of a parachute call for an explicit‐explicit Eulerian computational framework for FSI equipped with an embedded (or immersed) boundary method for CFD (for example, see the work of Huang et al). To this end, the Eulerian framework of the AERO Suite of codes developed at Stanford University for the simulation of highly nonlinear FSI phenomena with topological changes is chosen to perform the FSI simulations of parachute inflation dynamics (PID) reported herein.…”

Section: Application To the Simulation Of Parachute Inflation Dynamicsmentioning

confidence: 99%

Fast computation of the wall distance in unsteady Eulerian fluid‐structure computations

Grimberg

Farhat

2018

Numerical Methods in Fluids

Self Cite

View full text Add to dashboard Cite

Summary Highly nonlinear, turbulent, dynamic, fluid‐structure interaction problems characterized by large structural displacements and deformations, as well as self‐contact and topological changes, are encountered in many applications. For such problems, the Eulerian computational framework, which is often equipped with an embedded (or immersed) boundary method for computational fluid dynamics, is often the most appropriate framework. In many circumstances, it requires the computation of the time‐dependent distance from each active mesh vertex of the embedding mesh to the nearest embedded discrete surface. Such circumstances include, for example, modeling turbulence using the Spalart‐Allmaras or detached eddy simulation turbulence models and performing adaptive mesh refinement in order to track the boundary layer. Evaluating at each time step the distance to the wall is computationally prohibitive, particularly in the context of explicit‐explicit fluid‐structure time‐integration schemes. Hence, this paper presents two complementary approaches for reducing this computational cost. The first one recognizes that many quantities depending on the wall distance are relatively insensitive to its inaccurate evaluation in the far field. Therefore, it simplifies a state‐of‐the‐art algorithm for computing the wall distance accordingly. The second approach relies on an effective wall distance error estimator to update the evaluation of the wall distance function only when otherwise, a quantity of interest that depends on it would become tainted by an unacceptable level of error. The potential of combining both approaches for dramatically accelerating the computation of the wall distance is demonstrated with the Eulerian simulation of the inflation of a disk‐gap‐band parachute system in a supersonic airstream.

show abstract

“…For body-fitted CFD meshes, a "dressing" approach based on phantom surface elements and massless rigid elements was proposed in [16] to address this issue. While it can also be used in non body-fitted CFD frameworks [3], this approach has computational disadvantages that are identified and discussed in this work. For this reason, an alternative approach is proposed here for computing cable-driven fluid structure interactions.…”

Section: Introductionmentioning

confidence: 99%

An embedded boundary approach for resolving the contribution of cable subsystems to fully coupled fluid‐structure interaction

Huang

Avery

Farhat

2020

Numerical Meth Engineering

Self Cite

View full text Add to dashboard Cite

Cable subsystems characterized by one or more extremely long, slender and flexible structural elements are featured in numerous engineering systems including parachutes, suspension bridges, marine drilling risers, and aerial refueling equipment. In each one of these systems, interaction between the cable and its surrounding fluid is inevitable. However, the nature and consequences of such Fluid-Structure Interactions (FSI) have received relatively little attention in the computational mechanics open literature possibly due to an inherent complexity associated with resolution of the multiple scales present in problems of this type. This work proposes an embedded boundary approach for simulating the FSI of cable subsystems, in which the dynamics of the solid cable is captured by a discretization of the cable centerline using conventional beam elements which are typically available in commercial and open source finite element structural codes, while the geometry of the cable is represented within the fluid domain using a discrete -for instance, triangulated -embedded surface. The proposed approach is built on: master/slave kinematics between beam elements (master) and the embedded surface (slave); a highly accurate algorithm for computing the embedded surface displacement based on the beam displacement; and an energy-conserving method for transferring distributed forces and moments acting on the nodes of the discrete surface to the beam elements. Hence, both the flow-induced forces on the cable and the effect of the structural dynamic response of the cable on the nearby flow are taken into account. Moreover, the proposed model can be readily incorporated into the Eulerian computational framework, which enables handling of the large deformations of the cable subsystem. The effectiveness of the proposed approach is demonstrated using a model aerial refueling system and a challenging supersonic parachute inflation problem.

show abstract

Simulation of Parachute Inflation Dynamics Using an Eulerian Computational Framework for Fluid-Structure Interfaces Evolving in High-Speed Turbulent Flows

Cited by 34 publications

References 18 publications

High-order partitioned spectral deferred correction solvers for multiphysics problems

High-order partitioned spectral deferred correction solvers for multiphysics problems

Fast computation of the wall distance in unsteady Eulerian fluid‐structure computations

An embedded boundary approach for resolving the contribution of cable subsystems to fully coupled fluid‐structure interaction

Contact Info

Product

Resources

About