In this research, a CFD solver is developed for solving the 2D/3D compressible flowproblem: the finite volume method based on multi-block structural grids is used to solve thecompressible Reynolds averaged Navier–Stokes equations (RANS). Included in the methodology aremultiple high-order reconstruction schemes, such as the 3rd-order MUSCL (Monotone UpstreamcenteredSchemes for Conservation Laws), 5th-order WENO (Weight Essentially Non-Oscillatory),and 5th-order MP (Monotonicity-Preserving) schemes. Of the variety of turbulence models that areembedded, this solver is mainly based on the shear stress transport model (SST), which is compatiblewith OpenMP/MPI parallel algorithms. This research uses the CFD solver to conduct steady-stateflow simulation for a two-dimensional supersonic inlet/isolator, incorporating these high-precisionreconstruction schemes to accurately capture the shock wave/expansion wave interaction and shockwave/turbulent boundary layer interaction (SWTBLI), among other effects. By comparing the 2D/3Dcomputation results of the same inlet configuration, it is found that the 3D effects of the side wallcannot be ignored due to the existing strong lateral flow near the corner. To obtain a more refinedturbulence simulation, the commercial software ANSYS Fluent 18.0 is used to carry out the detachededdy simulation (DES) and the large eddy simulation (LES) of the same supersonic inlet, so as toreveal the flow details near the separation area and boundary layers.