Description |
As hypersonic aerospace vehicles are designed to increased performance specifications utilizing lighter weight, higher strength materials, fluid-structural interaction (FSI) effects become increasingly important to model, especially considering the increasing use of numerical models in many phases of design. When a fluid flows over a solid, a force is imparted on the solid and the solid deforms. This deformation, in turn, causes a change in the fluid flow field which modifies the force distribution on the structure. This FSI induced deformation is a primary area of study within the field of aeroelasticity. To further complicate the matter, thermodynamic and chemical effects are vitally important to model in the hypersonic flow regime. Traditionally, two separate numerical models are utilized to model the fluid and solid phases and a coupling algorithm accomplishes the task of modeling FSI. Coupling between the two solvers introduces numerical inaccuracies, inefficiencies, and for many mesh-based solvers, large deformations cannot be modeled. For this research, a combined Eulerian grid-based and Lagrangian particle-based solver known as the Material Point Method (MPM) is introduced and defined from prior research by others, and the particular MPM numerical code utilized in this research is outlined. The code combines the two separate solvers into a single numerical algorithm with separate constitutive relations for the fluid and solid phase, thereby allowing FSI modeling within a single computational framework. A limiter is applied to reduce numerical noise and oscillations around shock and expansion waves and exhibits a large reduction in oscillation amplitude and frequency. A Fourier's Law of Conduction heat transfer algorithm is implemented for heat transfer at a fluid-structure interface. The results from this heat transfer algorithm are compared with an independently developed numerical code for the single ramp case and experimental data for the double cone case. Finally, a reacting flow model is exhibited, the results are compared to other numerical solutions for verification and recommendations are made for further research. |