dougshidong / PHiLiP

Parallel High-Order Library for PDEs through hp-adaptive Discontinuous Galerkin methods
Other
45 stars 37 forks source link

boundary_inflow() does not give correct boundary values for supersonic inlet #237

Open ssaishruthi98 opened 1 year ago

ssaishruthi98 commented 1 year ago

As per the reference cited in the boundary_inflow() function, the boundary should be set to the prescribed density, velocity and pressure for supersonic inflow. However, when running the Shu-Osher test case with the inflow BC (and wall BC on the other end), the value of density initially drops then increases over time. The results obtained from using the boundary_inflow() function are as follows (initial condition top, final time bottom): shuOsher_StrongDG-LaxF_SupersonicInflowBC_time0s shuOsher_StrongDG-LaxF_SupersonicInflowBC_time1 8s

The expected results for the Shu-Osher test case are as follows (initial condition top, final time bottom): shuOsher_StrongDG-LaxF_RiemannBC_time0s shuOsher_StrongDG-LaxF_RiemannBC_time1 8s

This is obtained using the boundary_riemann() function which sets the boundary to the prescribed values, but ideally these results should be obtained when using the boundary_inflow() function as well.

The Shu-Osher test case has not been merged yet, but when it is merged, I will include the command to run the test!