

To examine the Riemann problem for a realistic fluid and determine the accuracy of the HLLC general EoS Riemann solver, we implemented an exact general EoS Riemann solver based on [@colella1985]. We improved the robustness of the root-finding algorithm. The existence and uniqueness of the solution are guaranteed if the EoS tables and the interpolated ones satisfy the two monotone condition in Equation \[eqn:unique1\]-\[eqn:unique2\]. The new exact general EoS Riemann solver can be used even when a convex condition $\left(\frac{\partial^{2}p}{\partial\rho^{2}}\right)_{s}>0$ is not satisfied. We applied the exact general EoS Riemann solver to the Godunov scheme in Section \[sec:exactgodunov\]. To examine the correctness of our new algorithm, we compared our exact general EoS Riemann solver to the exact ideal gas Riemann solver for shock tube tests of perfect gases with $\gamma=\Gamma=1.05$ and $\gamma=\Gamma=1.667$. The solutions are comparable, confirming that our exact general EoS Riemann reduces to the original exact ideal gas Riemann solver when using the corresponding perfect gas EoS. We also compared our HLLC general EoS Riemann solver to the original HLLC Riemann solver by running tests with perfect gas EoS, resulting in comparable solutions. Therefore, our HLLC general EoS Riemann solver reduces to the original HLLC Riemann solver in the case of a perfect gas.

In Section \[sec:hllcvsexact\], we compare the solution of the HLLC general EoS Riemann solver with two different EoS solvers to the exact solution calculated by the algorithm in Section \[sec:psolver\]. One EoS solver is called the tabulated EoS solver; it calculates the thermodynamic variables by interpolating 2D tables with a bi-linear interpolation method (Section \[sec:hllceossolver\]). Another EoS solver is called the analytic EoS solver; it calculates the thermodynamic variables by solving EoS at run-time (see Appendix \[sec:eos\] for details).

We find reasonably consistent results between the approximate and exact Riemann solvers for the shock-tube and two-shock tests. The HLLC solver can resolve the contact discontinuity, which is different from the HLL solver [@kannan2019]. The location and the strength of the shock can be well captured. Although discrepancy is found in the two-rarefaction-wave test
To examine the Riemann problem for a realistic fluid and determine the accuracy of the HLLC cosmopolitan EoS Riemann problem solver, we implement an exact general EoS Riemann problem solver based on [ @colella1985 ]. We improved the robustness of the root - determination algorithm. The universe and uniqueness of the solution are guarantee if the EoS tables and the interpolated one satisfy the two monotone circumstance in Equation \[eqn: unique1\]-\[eqn: unique2\ ]. The new exact general EoS Riemann solver can be used even when a convex condition $ \left(\frac{\partial^{2}p}{\partial\rho^{2}}\right)_{s}>0 $ is not satisfied. We apply the exact general EoS Riemann solver to the Godunov scheme in Section \[sec: exactgodunov\ ]. To test the correctness of our new algorithm, we compared our exact cosmopolitan EoS Riemann solver to the exact ideal natural gas Riemann solver for shock tube test of perfect gases with $ \gamma=\Gamma=1.05 $ and $ \gamma=\Gamma=1.667$. The solutions are comparable, confirming that our exact cosmopolitan EoS Riemann reduce to the original exact ideal gas Riemann solver when using the corresponding perfect gas EoS. We also compare our HLLC general EoS Riemann solver to the original HLLC Riemann solver by running examination with arrant gas EoS, resulting in comparable solutions. Therefore, our HLLC general EoS Riemann solver reduces to the original HLLC Riemann solver in the case of a arrant gas. 

 In Section \[sec: hllcvsexact\ ], we compare the solution of the HLLC general EoS Riemann solver with two different EoS solvers to the exact solution calculated by the algorithm in Section \[sec: psolver\ ]. One EoS solver is called the tabulated EoS solver; it calculates the thermodynamic variables by interpolating 2D table with a bi - analogue interjection method acting (Section \[sec: hllceossolver\ ]). Another EoS solver is called the analytic EoS solver; it calculates the thermodynamic variables by solving EoS at run - time (examine Appendix   \[sec: eos\ ] for details). 

 We find reasonably consistent resultant role between the approximate and exact Riemann problem solver for the shock - tube and two - shock tests. The HLLC solver can decide the contact discontinuity, which is unlike from the HLL solver [ @kannan2019 ]. The location and the forte of the shock can be well captured. Although discrepancy is discover in the two - rarefaction - wave test


To fxamine the Riemann probuem for a realistic flnid and determive the accuracy of the HLLC jeneeal EiS Riemann solver, we ioplementef an exaxt gtneral EoS Riemani solver based ln [@cmoella1985]. We improyed the robgstness of the ruoc-finding algorithm. The existence and uniquemeds of the solujion swe gharanteed if the EoS tables and ths interkolated ones satisgy the two monotone conditlon ln Equation \[eqn:unieue1\]-\[eqn:uniqug2\]. Thq new exact gdneral EoS Riemann solber can be used even when a conxex cpndition $\lgyr(\frww{\partial^{2}p}{\pactial\rro^{2}}\right)_{s}>0$ is kpt sathsfied. Ee applied the exect teneral EoS Riemann smlver to the Godunjv scheme iu Section \[sec:exactgodyniv\]. To exakine rhe codrxcthess ov onr new algodithm, we conpared our exact gemewql EoS Riemanh solvqr to the exact ideal gas Riemann solver xor shock tube tests of peefect gases with $\gammw=\Gamma=1.05$ anq $\gamma=\Gamma=1.667$. The solutions are comparable, confirmhng tiag obv dzaft general EoS Riemann reduces to the originaj ecabt ideal gas Riemcnn solver when islnb the correspovding 'srrect gas EoS. We aldo compwred iur HLLC denetal EoS Riemann solver to tye original KLLX Riemann solver bv running terts eith lerfect gas EoS, resultiug in domparable dolutions. Gherefore, our HLUC nenaral EoS Riemann solver requces to vhe oxiginal FLLC Riemagn solver ln thc case of a perfect gas.

Iu Secdion \[sec:hlpcvsexact\], we compare the solution of the HLLC ganegal EoS Rnemann solver with ewo different GoS solvexs to ghe exact aolutioi calculated by the algorhjhm in Sectioi \[sec:psolder\]. Ine WoS solxdr is called tne tabulaneb EoS solcer; it calculates bhe tfsrmodynamic varncboes by interpolstivg 2Q nabkef with a bi-litear ingrrpolxtion mtchob (Rectoon \[sec:hllceossolver\]). Anofher EoS solver is cwlled thg analytis EoS solver; ot calculates the uhermovynamir varisblgs by solving EoS at run-time (ses Appendid \[seg:eos\] for detayls).

Wc figd reasonayly consistent results between the approeimate and exact Riemanb solvers for the skogk-tube and txo-shoch tests. Tve HLLC solver can rwsolve the contacu discontinuity, which ia diffarent from the HLL solver [@kannan2019]. The location and the strength of the shock cab be wxlj captured. Zlthpugh giferepancr is hound in the two-rsrefaction-wave test
To examine the Riemann problem for a and the accuracy the HLLC general an general EoS Riemann based on [@colella1985]. improved the robustness of the root-finding The existence and uniqueness of the solution are guaranteed if the EoS tables the interpolated ones satisfy the two monotone condition in Equation \[eqn:unique1\]-\[eqn:unique2\]. The new general Riemann can used even when a convex condition $\left(\frac{\partial^{2}p}{\partial\rho^{2}}\right)_{s}>0$ is not satisfied. We applied the exact general EoS solver to the Godunov scheme in Section \[sec:exactgodunov\]. examine the correctness of new algorithm, we compared our general Riemann solver the ideal Riemann solver for tube tests of perfect gases with $\gamma=\Gamma=1.05$ and $\gamma=\Gamma=1.667$. The solutions are comparable, confirming that our exact EoS Riemann the original ideal Riemann when using the gas EoS. We also compared our Riemann solver to the original HLLC Riemann solver running tests perfect gas EoS, resulting in comparable Therefore, our HLLC general EoS Riemann solver reduces the original HLLC Riemann solver in the case of a perfect gas. In Section \[sec:hllcvsexact\], the solution of the general EoS Riemann with different solvers the exact calculated by the algorithm in Section \[sec:psolver\]. One EoS solver is the tabulated EoS solver; it calculates the thermodynamic variables by tables a bi-linear interpolation (Section \[sec:hllceossolver\]). Another EoS is the analytic EoS solver; the variables at (see \[sec:eos\] for details). We reasonably consistent results between the and exact Riemann solvers tests. The HLLC solver can resolve the contact which is different from the HLL solver The location and the strength of the shock can be well captured. discrepancy is the two-rarefaction-wave test


To examine the Riemann probleM for a realiStic fLuiD anD dEterMine The accuracy of tHE HLLc general EoS Riemann solvEr, we iMpLEmenTEd An exaCt generAL EOs rieMaNn SolVeR BaSed on [@ColElla1985]. We iMproved the RobUsTness of the roOT-fInding algoRitHm. The existenCe aNd uniqUeNesS Of the SolUtion Are guaRAnteed If the EoS tAbLEs and tHE interpOLAtEd onEs satisfy the two moNOtONe condition in EQuatioN \[eQN:uNIQue1\]-\[Eqn:Unique2\]. The nEw Exact GEneral EOs RIEMAnn SOlver can be useD even when a cONveX condiTiOn $\lEFt(\frac{\PartiAl^{2}P}{\ParTial\rho^{2}}\righT)_{s}>0$ is Not satisfIed. We aPPlied thE Exact geNeral EOS RIemAnn sOLvEr To tHe gOduNOv SchEMe iN Section \[SeC:eXactgOdunOV\]. tO ExamIne The cOrrecTness of our new AlgOritHM, we CompaRed ouR exaCt GenerAl EoS RIemanN sOlver to the exact IdeaL gas RiemaNn sOlVer FoR shocK Tube teSts Of pErfect gAses witH $\GamMa=\gAMMa=1.05$ And $\gamma=\Gamma=1.667$. The soLuTIOnS are compArable, COnFiRMing that OuR exAct gENEral EOS RiEMaNn reduceS to the ORiGiNal exacT iDeal gaS RIemAnn SolveR When Using tHe corresPondiNG perfect gas EoS. wE also compared OUr hllC GEnerAl EOS Riemann soLver TO the OrigINaL HLlc RiemAnn soLvER bY Running tests with perFeCt gas EOS, resUlting in compaRable solutIONS. TherefoRe, ouR hLlc general EoS RieMann sOlver reducES to the orIginaL HLLC RieMann solveR IN the case Of a PerFecT gaS.

iN SEction \[sec:hllcVSExacT\], wE comparE thE solutiOn oF thE HLlC gEnEral EoS RiEmann solVeR wItH tWo dIfferENt EoS solVeRs tO tHe eXact sOLution CalcuLateD bY tHE alGorithm IN SECTion \[SeC:pSolvEr\]. ONe eoS soLver IS caLled the Tabulated eoS SOlveR; iT cAlculatEs the thermodyNaMic variablEs By iNterpoLATing 2D tabLes with a bi-linear interpoLAtion meThoD (SectIon \[sEc:hllceosSolVer\]). AnoTheR eoS solVer is cAlled ThE anALYtic EOs SoLveR; iT calculateS THe tHermoDyNamiC variabLes by solving EoS at rUN-tiMe (see Appendix \[Sec:Eos\] fOR DeTaiLS).

WE FinD rEAsoNABly consistent reSults betweEn THe ApproximatE And ExAct RiemAnn solvErs foR The shocK-tube and tWo-shock teStS. The hllC sOlver can reSolve the Contact diSContiNUiTy, whiCh iS diffeReNt fRom thE HLL soLVer [@KannaN2019]. The loCaTion anD the sTrEngth of tHe shock can be well captureD. AlthoUgh diScrEpancy is fOunD In tHe two-rareFactIon-wave tesT


To examine the Riemann p roblem for  a re ali sti cflui d an d determine th e  acc uracy of the HLLC gene ral E oS Riem a nn  solv er, wei mp l e men te danex a ct  gene ral  EoS Ri emann solv erba sed on [@col e ll a1985]. We  im proved the r obu stness  o f t h e roo t-f indin g algo r ithm.The exist en c e andu niquene s s  o f th e solution are gu a ra n teed if the Eo S tabl es an d the  in terpolated  o nes s a tisfy t h et w o  mo n otone conditi on in Equat i on\[eqn: un iqu e 1\]-\[ eqn:u ni q ue2 \]. The new  exa ct genera l EoSR iemanns olver c an beuse d e venw he na c on v exc on dit i on$\left(\ fr ac {\par tial ^ { 2 } p}{\ par tial \rho^ {2}}\right)_{ s}> 0$ i s  no t sat isfie d. W eappli ed the  exac tgeneral EoS Rie mann  solver t o t he  Go du nov s c heme i n S ect ion \[s ec:exac t god un o v \ ].  To examine the co rr e c tn ess of o ur new al go r ithm, we  c omp ared o ur ex actg en eral EoS  Riema n nso lver to  t he exa ct  id eal  gasR iema nn sol ver forshock tube tests ofp erfect gasesw it h $\ g amma =\G amma=1.05$and$ \gam ma=\ G am ma= 1 .667$ . The  s o lu t ions are comparable ,confir mingthat our exac t generalE o S  Riemann  red u ce s  to the origin al ex act idealg as Riema nn so lver whe n using t h e  corresp ond ing  pe rfe c t  g as EoS. We al s o  com pa red our  HL LC gene ral  Eo S R iem an n solverto the o ri gi na lHLL C Rie m ann solv er  by  r unn ing t e sts wi th pe rfec tga s  Eo S, resu l ti n g  inco mp arab leso lutio ns.T her efore,our HLLCgen e ralEo SRiemann  solver reduc es  to the or ig ina l HLLC R iemann s olver in the case of ap erfectgas .

In  Sec tion \[se c:h llcvse xac t \], we  compa re th esol u t ion o f th e H LL C generalE o S R ieman nsolv er with  two different EoS sol vers to the e xac t so l u ti onc al c ula te d  by t he algorithm in  Section \ [s e c: psolver\]. One  E oS solv er is c alled the tab ulated Eo S solver;  i t ca l c ula tes the th ermodyna mic varia b les b y  i nterp ola ting 2 Dtab les w ith ab i-l inear  inter po lation  meth od  (Sectio n \[sec:hllceossolver\] ). Ano therEoS  solver i s c a lle d the ana lyti c EoS solv er;  it  calc ula t es th e th e rm ody n amicvari a bles by s o lv ing E oS  at run-tim e ( see  Appe ndi x  \[sec :eos \] for details).We find reason ably c ons ist e nt r es ults between t heap p r oximatean d exact Rie mann sol ve r s for  the s hock-t ube and t wo - shocktest s.The HLLCsol ve r  can re so lv e  the c onta ct  disco ntinui t y, w h i ch is differentfromt h e HLL sol ver [ @k annan20 1 9].The locati on and thestreng th o f the  shockca n be w ell  c aptured. A l though di screp ancy is  f ound  in  the t wo-r a r efact ion- wa vetest


To examine_the Riemann_problem for a realistic_fluid and_determine_the accuracy_of_the HLLC general_EoS Riemann solver,_we implemented an exact_general EoS Riemann_solver_based on [@colella1985]. We improved the robustness of the root-finding algorithm. The existence and_uniqueness_of the_solution_are_guaranteed if the EoS tables_and the interpolated ones satisfy_the two_monotone condition in Equation \[eqn:unique1\]-\[eqn:unique2\]. The new exact_general_EoS Riemann solver_can be used even when a convex condition $\left(\frac{\partial^{2}p}{\partial\rho^{2}}\right)_{s}>0$_is not satisfied. We applied the_exact general EoS_Riemann_solver_to the Godunov scheme_in Section \[sec:exactgodunov\]. To examine the_correctness of our new algorithm, we_compared our exact general EoS Riemann solver_to the exact ideal gas Riemann_solver for shock tube tests_of perfect_gases with $\gamma=\Gamma=1.05$ and $\gamma=\Gamma=1.667$._The solutions are_comparable, confirming_that our exact_general EoS Riemann reduces to the_original exact ideal_gas Riemann solver when using the_corresponding_perfect gas EoS._We_also_compared our_HLLC general EoS_Riemann_solver to_the_original HLLC Riemann solver by running_tests_with perfect gas EoS, resulting in comparable_solutions. Therefore, our HLLC_general_EoS Riemann solver reduces_to the original HLLC Riemann_solver in the case of a_perfect gas.

In_Section \[sec:hllcvsexact\],_we compare the solution of the HLLC general EoS Riemann solver_with two different EoS solvers to_the exact solution calculated_by the_algorithm_in Section \[sec:psolver\]._One_EoS solver_is called the tabulated EoS solver; it_calculates the_thermodynamic variables by interpolating 2D tables_with a bi-linear interpolation_method_(Section \[sec:hllceossolver\]). Another EoS solver is_called the analytic EoS solver; it_calculates the thermodynamic variables by_solving_EoS_at run-time (see Appendix \[sec:eos\] for_details).

We find reasonably consistent results between_the approximate and_exact Riemann solvers for the shock-tube and_two-shock_tests. The HLLC solver can resolve_the_contact discontinuity, which is different from_the_HLL_solver [@kannan2019]. The location and_the strength of the shock can_be well captured. Although discrepancy is found in the_two-rarefaction-wave test