Abstract:An immersed boundary method for the fluid--structure--thermal interaction in rarefied gas flow is presented. In this method, the slip model is incorporated with the penalty immersed boundary method to address the velocity and temperature jump conditions at the fluid--structure interface in rarefied gas flow within slip regime. In this method, the compressible flow governed by Navier-Stokes equations are solved by using high-order finite difference method; the elastic solid is solved by using finite element method; the fluid and solid are solved independently and the fluid--structure--thermal interaction are achieved by using a penalty method in a partitioned way. Several validations are conducted including Poiseuille flow in a 2D pipe, flow around a 2D NACA airfoil, moving square cylinder in a 2D pipe, flow around a sphere and moving sphere in quiescent flow. The numerical results from present method show good agreement with the previous published data obtained by other methods, and it confirms the the good ability of the proposed method in handling fluid--structure--thermal interaction for both weakly compressible and highly compressible rarefied gas flow. To overcome the incapability of Navier-Stokes equations at high local Knudsen numbers in supersonic flow, an artificial viscosity is introduced to ease the sharp transition at the shock wave front. Inspired by Martian exploration, the application of proposed method to study the aerodynamics of flapping wing in rarefied gas flow is conducted in both 2D and 3D domains, to obtain some insights for the flapping-wing aerial vehicles operating in Martian environment.
What problem does this paper attempt to address?
This paper attempts to solve the problem of modeling fluid - structure - thermal interaction (FSTI) in rarefied gas flows. Specifically, the author proposes a method based on the Immersed Boundary Method (IBM) to handle the velocity slip and temperature jump conditions in rarefied gas flows within the slip region. This method solves the problems of velocity and temperature jumps at the fluid - structure interface by combining the slip model with the penalty immersed boundary method.
### Main problems
1. **Fluid - structure - thermal interaction in rarefied gas flows**:
- Fluid - structure - thermal interaction (FSTI) in rarefied gas flows is a complex problem, especially in the slip flow regime, that is, the region where the Knudsen number is between 0.001 and 0.1. Traditional Navier - Stokes equations have limitations in dealing with such problems because they assume that the fluid is continuous and do not consider the slip effect.
2. **Velocity slip and temperature jump conditions**:
- In rarefied gas flows, the fluid can slide on the solid surface, and non - uniform temperature can cause the gas to slide from the cold region to the hot region, a phenomenon known as "thermal creep". These effects need to be accurately modeled in numerical simulations to predict fluid - structure - thermal interaction.
3. **Failure of the continuum assumption at high local Knudsen numbers**:
- In supersonic flows, high local Knudsen numbers before and after the shock wave can cause the failure of the continuum assumption, thus affecting the accuracy of the Navier - Stokes equations. To overcome this problem, the author introduces artificial viscosity to smooth the sharp changes in front of the shock wave.
### Solutions
1. **Combination of the slip model and the penalty immersed boundary method**:
- The author combines the slip model with the penalty immersed boundary method (pIB) to handle the velocity slip and temperature jump conditions at the fluid - structure interface. The slip model calculates the slip velocity and temperature jump based on the local velocity and temperature gradients.
2. **High - order finite difference method and finite element method**:
- For the hydrodynamic part, the high - order finite difference method is used to solve the Navier - Stokes equations for compressible fluids; for the elastic solid part, the finite element method is used to solve.
3. **Artificial viscosity**:
- To deal with the problem of the failure of the continuum assumption at high local Knudsen numbers, the author introduces artificial viscosity, which depends linearly on the local Knudsen number. When the local Knudsen number is greater than 0.1, the artificial viscosity will increase, thus smoothing the sharp changes in front of the shock wave.
### Verification and application
- **Verification**:
- The author verifies the accuracy of the proposed method through multiple verification cases, such as Poiseuille flow, flow around NACA airfoil, moving square cylinder, flow around sphere, etc. The numerical results are in good agreement with previous research results, indicating that this method has good ability in dealing with fluid - structure - thermal interaction in rarefied gas flows.
- **Application**:
- Inspired by Mars exploration, the author applies the proposed method to study the dynamics of flapping - wing aircraft in rarefied gas flows to obtain some insights into the operation of flapping - wing aircraft in the Martian environment.
In summary, this paper aims to solve the problem of modeling fluid - structure - thermal interaction in rarefied gas flows by proposing a new numerical method, especially dealing with the velocity slip and temperature jump conditions, and the problem of the failure of the continuum assumption at high local Knudsen numbers.