 been mentioned by E. F. Toro in [@toro2013] (the last remark in Section 10.4.2, Chapter 10), however, there has been a lack of systematic assessment of the accuracy and efficiency of the HLLC type general EoS Riemann solvers.

In this paper, we devise and examine a new efficient HLLC [@toro1994] general EoS Riemann solver that can handle realistic gases whose $C_{V}$, $C_{p}$ and $\Gamma$ may vary drastically when the temperature changes. The rapidly varying thermal dynamic variables may become stiff to the Euler equations. Therefore an approximate Riemann solver cannot stand alone. To justify the accuracy of the HLLC general EoS Riemann solver, we design and program a new exact general EoS Riemann solver that can work with non-convex EoS and compare the solution of the HLLC general EoS Riemann solver to the exact Riemann solver. We point out that the exact general EoS Riemann solver is the true solution with the EoS we use. Our exact general EoS Riemann solver only needs two monotone conditions, $$\begin{aligned}
    \left(\frac{\partial p}{\partial\rho}\right)_{\epsilon}>0 \label{eqn:unique1}\\
    \left(\frac{\partial p}{\partial\epsilon}\right)_{\rho}>0  \label{eqn:unique2}\end{aligned}$$ discussed in [@liu1975; @smith1979]. An additional convex condition $(\partial^{2}p/\partial\rho^{2})_{s}>0$ was required by [@colella1985], however we find that this is not a requirement of the new exact general EoS Riemann solver discussed in this paper (see Section \[sec:pnumericalmethod\] for further discussion).

In Section \[sec:psolver\], we present the exact general EoS Riemann solver. We introduce the application of the exact general EoS Riemann solver to the Godunov scheme [@godunov1959] in Section \[sec:exactgodunov\]. The novel HLLC general EoS Riemann solver will be discussed in Section \[sec:hllc\]. Section \[sec:num\] will focus on the numerical results generated by the exact and the HLLC general EoS Riemann solvers. In particular, we will establish the correspondence of the exact general EoS Riemann solver to the perfect gas Riemann solver, and the HLLC general EoS Riemann solver to the original HLLC Riemann solver. We
been mentioned by E. F. Toro in [ @toro2013 ] (the last remark in Section 10.4.2, Chapter 10), however, there has been a lack of systematic assessment of the accuracy and efficiency of the HLLC character cosmopolitan EoS Riemann problem solver. 

 In this paper, we devise and examine a raw efficient HLLC [ @toro1994 ] general EoS Riemann solver that can handle naturalistic gases whose $ C_{V}$, $ C_{p}$ and $ \Gamma$ may vary drastically when the temperature changes. The quickly vary thermal dynamic variable may become stiff to the Euler equations. consequently an approximate Riemann solver cannot stand alone. To justify the accuracy of the HLLC general EoS Riemann solver, we design and program a fresh exact general EoS Riemann solver that can work with non - convex EoS and compare the solution of the HLLC cosmopolitan EoS Riemann solver to the exact Riemann solver. We bespeak out that the exact cosmopolitan EoS Riemann solver is the on-key solution with the EoS we use. Our exact general EoS Riemann solver merely needs two monotone conditions, $ $ \begin{aligned } 
     \left(\frac{\partial p}{\partial\rho}\right)_{\epsilon}>0 \label{eqn: unique1}\\ 
     \left(\frac{\partial p}{\partial\epsilon}\right)_{\rho}>0   \label{eqn: unique2}\end{aligned}$$ discussed in [ @liu1975; @smith1979 ]. An additional convex condition $ (\partial^{2}p/\partial\rho^{2})_{s}>0 $ was required by [ @colella1985 ], however we find that this is not a requirement of the new exact general EoS Riemann solver discussed in this paper (see Section \[sec: pnumericalmethod\ ] for further discussion). 

 In Section \[sec: psolver\ ], we present the accurate cosmopolitan EoS Riemann solver. We introduce the application of the exact cosmopolitan EoS Riemann solver to the Godunov outline [ @godunov1959 ] in Section \[sec: exactgodunov\ ]. The novel HLLC general EoS Riemann solver will be discourse in Section \[sec: hllc\ ]. Section \[sec: num\ ] will focus on the numerical results generate by the exact and the HLLC general EoS Riemann solvers. In particular, we will establish the correspondence of the exact cosmopolitan EoS Riemann solver to the perfect gas Riemann solver, and the HLLC general EoS Riemann solver to the original HLLC Riemann solver. We
 befn mentioned by E. F. Toro in [@toro2013] (the last remack in Ssction 10.4.2, Zhapter 10), however, there has bxen q lacj of systematic assessoent of tje accurqcy end efficiency oh the HLLC typs genzrel EoS Riemann xolvers.

In dhis paper, we gexide and examine a new efficient HLLC [@toro1994] grnfral EoS Riemagn spjver nhct can handle realistic gases wgose $C_{V}$, $C_{p}$ and $\Gamma$ may vary drastically when the temperature changfs. The rapiely dqrying thermxl dynamic variables mzy become stiff to the Euler eqjatious. Thereforg zn dpproximate Riemwnn solver cannot stdnd alome. To justify bhe arcurqcy of the HLLC generel EoS Riemann solvet, we desigt cnd program a new exaxt genetal EmS Rkwmavn aokvsr thah cen work witg non-convex EoS and compare tht sjoution of the HLLC degeral EoS Riemann solver to the exact Rpemahn solver. We point out rhat the exact generap EoS Rieiann solver is the true solution with the EoS we gse. Onr excgb gevwrwl EoS Riemann solver only needs two monotone domdptions, $$\begin{aligncd}
    \left(\frac{\partoap l}{\kartial\rho}\righj)_{\epsilou}>0 \lzbel{eqn:unique1}\\
    \lfft(\frac{\kartiao p}{\partiaj\epsolon}\right)_{\rho}>0  \label{eqn:uniquw2}\end{aligned}$$ bisxussed in [@liu1975; @smitk1979]. An additioual comvex vondition $(\partial^{2}p/\particl\rho^{2})_{s}>0$ was requirfd by [@colsula1985], however we fknd tvat this is not a requiremqnt of thx new exact eenetal EoS Riemann sllver discussed in this oaper (sae Section \[sec:pnumericalmethod\] for furthec discussion).

Im Vecnion \[sec:pfolvev\], we present thq exact generak EoS Rnemann solver. We introdnce the applycation of tha exact generel EoS Riqmanb soover to ghe Godunov scneme [@goduujv1959] ib Section \[sec:exactnodunub\]. The novel HLLE teneral EoS Riekanv sjlner wynl be discusved kn Rrctiov \[sec:hllc\]. Scctkon \[xec:num\] will focus on the numerical results gcnerated vy the evact and the NLLC general EoS Rlemani solvxrs. In patticular, we will establish the dorresponfenge of the exast gcnerwl EoS Riekann solver to the perfect gas Riemann wolver, and the HLLX general EoS Riemaun solver to vhe oryginal HLNC Riemann solver. We
been mentioned by E. F. Toro in last in Section Chapter 10), however, of assessment of the and efficiency of HLLC type general EoS Riemann solvers. this paper, we devise and examine a new efficient HLLC [@toro1994] general EoS solver that can handle realistic gases whose $C_{V}$, $C_{p}$ and $\Gamma$ may vary when temperature The varying thermal dynamic variables may become stiff to the Euler equations. Therefore an approximate Riemann solver stand alone. To justify the accuracy of the general EoS Riemann solver, design and program a new general Riemann solver can with EoS and compare solution of the HLLC general EoS Riemann solver to the exact Riemann solver. We point out that exact general solver is true with EoS we use. general EoS Riemann solver only needs $$\begin{aligned} \left(\frac{\partial p}{\partial\rho}\right)_{\epsilon}>0 \label{eqn:unique1}\\ \left(\frac{\partial p}{\partial\epsilon}\right)_{\rho}>0 \label{eqn:unique2}\end{aligned}$$ discussed [@liu1975; @smith1979]. additional convex condition $(\partial^{2}p/\partial\rho^{2})_{s}>0$ was required [@colella1985], however we find that this is not requirement of the new exact general EoS Riemann solver discussed in this paper (see Section further discussion). In Section we present the general Riemann We the application the exact general EoS Riemann solver to the Godunov scheme [@godunov1959] Section \[sec:exactgodunov\]. The novel HLLC general EoS Riemann solver will in \[sec:hllc\]. Section \[sec:num\] focus on the numerical generated the exact and the EoS solvers. will the of the exact general Riemann solver to the perfect Riemann solver, and the to the original HLLC Riemann solver. We
 been mentioned by E. F. Toro in [@toRo2013] (the last rEmark In SEctIoN 10.4.2, ChaPter 10), However, there haS Been A lack of systematic assesSment Of THe acCUrAcy anD efficiENcY OF thE HlLc tyPe GEnEral EOS RIemann sOlvers.

In thIs pApEr, we devise anD ExAmine a new eFfiCient HLLC [@torO1994] geNeral EOS rieMAnn soLveR that Can hanDLe realIstic gaseS wHOse $C_{V}$, $C_{P}$ And $\GammA$ MAy Vary Drastically when thE TeMPerature changeS. The raPiDLy VARyiNg tHermal dynaMiC variABles may BEcOME StiFF to the Euler eqUations. TherEForE an appRoXimATe RiemAnn soLvER caNnot stand alOne. TO justify tHe accuRAcy of thE hLLC genEral Eos RiEmaNn soLVeR, wE deSiGN anD PrOgrAM a nEw exact gEnErAl EoS riemANN SOlveR thAt caN work With non-convex eoS And cOMpaRe the SolutIon oF tHe HLLc generAl EoS riEmann solver to thE exaCt Riemann SolVeR. We PoInt ouT That thE exAct General eoS RiemANn sOlVER Is The true solution witH tHE eos we use. OuR exact GEnErAL EoS RiemAnN soLver ONLy neeDs twO MoNotone coNditioNS, $$\bEgIn{alignEd}
    \Left(\frAc{\ParTiaL p}{\parTIal\rHo}\righT)_{\epsilon}>0 \Label{EQn:unique1}\\
    \left(\frAC{\partial p}{\partIAl\EPSiLOn}\riGht)_{\Rho}>0  \label{eqn:UniqUE2}\end{AligNEd}$$ DisCUssed In [@liu1975; @SmITh1979]. aN additional convex coNdItion $(\pArtiaL^{2}p/\partial\rho^{2})_{s}>0$ Was requireD BY [@Colella1985], hOwevER wE Find that this is Not a rEquirement OF the new eXact gEneral Eos Riemann sOLVer discuSseD in ThiS paPER (sEe Section \[sec:pNUMeriCaLmethod\] For Further DisCusSioN).

In seCtion \[sec:pSolver\], we PrEsEnT tHe eXact gENeral EoS riEmaNn SolVer. We INtroduCe the ApplIcAtIOn oF the exaCT gENEral eos RIemaNn sOlVer to The GODunOv schemE [@godunov1959] iN SeCTion \[SeC:eXactgodUnov\]. The novel HlLc general Eos RIemAnn solVER will be dIscussed in Section \[sec:hllC\]. section \[Sec:Num\] wiLl foCus on the nUmeRical rEsuLTs geneRated bY the eXaCt aND The HLlc GeNerAl eoS Riemann SOLveRs. In pArTicuLar, we wiLl establish the corrESpoNdence of the exAct GeneRAL EOS RIEmANn sOlVEr tO THe perfect gas RieMann solver, AnD ThE HLLC generAL Eos RIemann sOlver to The orIGinal HLlC Riemann Solver. We
 been mentioned by E. F. T oro in [@t oro20 13]  (t he  las t re mark in Sectio n  10. 4.2, Chapter 10), howe ver,th e re h a sbeena lacko fs y ste ma ti c a ss e ss mentofthe acc uracy andeff ic iency of the HL LC type ge ner al EoS Riema nnsolver s. 

I n  this  pa per,we dev i se and  examinean ew eff i cient H L L C[@to ro1994] general E o SR iemann solverthat c an ha n d lerea listic gas es  whos e  $C_{V} $ ,$ C _ {p} $  and $\Gamma$  may vary d r ast ically  w hen the te mpera tu r e c hanges. The  rap idly vary ing th e rmal dy n amic va riable s m aybeco m est iff  t o  th e  E ule r  eq uations.  T he refor e an a p p roxi mat e Ri emann  solver canno t s tand alo ne. T o jus tify  t he ac curacy  of t he  HLLC general E oS R iemann so lve r,  we  d esign and pr ogr ama new e xact ge n era lE o S  R iemann solver that  c a n  w ork with  non-c o nv ex EoS andco mpa re t h e  solu tion of  the HLL C gene r al  E oS Riem an n solv er  to  th e exa c t Ri emannsolver.We po i nt out that th e  exact genera l  E o S  R i eman n s olver is th e tr u e so luti o nwit h  theEoS w eu se .  Our exact generalEo S Riem ann s olver only ne eds two mo n o t one cond itio n s, $$\begin{align ed}
   \left(\ f rac{\par tialp}{\part ial\rho}\ r i ght)_{\e psi lon }>0  \l a b el {eqn:unique1} \ \ 
   \ left(\f rac {\parti alp}{ \pa rti al \epsilon} \right)_ {\ rh o} >0   \ label { eqn:uniq ue 2}\ en d{a ligne d }$$ di scuss ed i n[@ l iu1 975; @s m it h 1 979] .An  add iti on al co nvex con dition$(\partia l^{ 2 }p/\ pa rt ial\rho ^{2})_{s}>0$wa s required  b y [ @colel l a 1985], h owever we find that thi s  is not  arequi reme nt of the  ne w exac t g e neralEoS Ri emann  s olv e r  disc u s se d i nthis paper ( see  Sect io n \[ sec:pnu mericalmethod\] fo r  fu rther discuss ion ).
I n  S ect i on \[s ec : pso l v er\], we presen t the exac tg en eral EoS R i ema nn  solver . We in trodu c e the a pplicatio n of theex actg e ner al EoS Rie mann sol ver to th e  Godu n ov  sche me[@godu no v19 59] i n Sect i on\[sec :exact go dunov\ ]. Th enovel HL LC general EoS Riemannsolver  will  be  discusse d i n  Se ction \[s ec:h llc\]. Sec tio n \ [sec: num \ ] wil l fo c us  on the n umer i cal resul t sgen e r at ed by the e x a c t a nd th e H L LC gen eral  EoS Riemann solv e rs. In particu lar, w e w ill esta bl ish the corres pon de n c e of the  e xact genera l EoS Ri em a nn so lver t o theperfect g as Rieman n so lve r, and th e H LL C  genera lEo S  Riema nn s ol ver to  the o r igin a l  HLLC Riemann so lver. W e
 been_mentioned by_E. F. Toro in_[@toro2013] (the_last_remark in_Section_10.4.2, Chapter 10),_however, there has_been a lack of_systematic assessment of_the_accuracy and efficiency of the HLLC type general EoS Riemann solvers.

In this paper, we_devise_and examine_a_new_efficient HLLC [@toro1994] general EoS_Riemann solver that can handle_realistic gases_whose $C_{V}$, $C_{p}$ and $\Gamma$ may vary drastically_when_the temperature changes._The rapidly varying thermal dynamic variables may become stiff_to the Euler equations. Therefore an_approximate Riemann solver_cannot_stand_alone. To justify the_accuracy of the HLLC general EoS_Riemann solver, we design and program_a new exact general EoS Riemann solver_that can work with non-convex EoS_and compare the solution of_the HLLC_general EoS Riemann solver to_the exact Riemann_solver. We_point out that_the exact general EoS Riemann solver_is the true_solution with the EoS we use._Our_exact general EoS_Riemann_solver_only needs_two monotone conditions,_$$\begin{aligned}
_ __\left(\frac{\partial p}{\partial\rho}\right)_{\epsilon}>0 \label{eqn:unique1}\\
   _\left(\frac{\partial_p}{\partial\epsilon}\right)_{\rho}>0  \label{eqn:unique2}\end{aligned}$$ discussed in [@liu1975; @smith1979]._An additional convex condition_$(\partial^{2}p/\partial\rho^{2})_{s}>0$_was required by [@colella1985],_however we find that this_is not a requirement of the_new exact_general EoS_Riemann solver discussed in this paper (see Section \[sec:pnumericalmethod\] for further_discussion).

In Section \[sec:psolver\], we present the_exact general EoS Riemann_solver. We_introduce_the application of_the_exact general_EoS Riemann solver to the Godunov scheme_[@godunov1959] in_Section \[sec:exactgodunov\]. The novel HLLC general_EoS Riemann solver will_be_discussed in Section \[sec:hllc\]. Section \[sec:num\]_will focus on the numerical results_generated by the exact and_the_HLLC_general EoS Riemann solvers. In_particular, we will establish the correspondence_of the exact_general EoS Riemann solver to the perfect_gas_Riemann solver, and the HLLC general_EoS_Riemann solver to the original HLLC_Riemann_solver._We