 as a divergence. With this, it is possible to construct numerical fluxes $\boldsymbol{E}$ using the same KNP flux formula as before:

$$\begin{aligned}
\label{knpfluxmhd}
\boldsymbol{E}^x_{i+\frac{1}{2},j,k}=&\frac{1}{a^+_{i+\frac{1}{2},j,k}-a^-_{i+\frac{1}{2},j,k}}
\Bigg[ 
   a^+_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^E_{i,j,k} -\\ \notag
   &- a^-_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^W_{i+1,j,k}
   + a^+_{i+\frac{1}{2},j,k} a^-_{i+\frac{1}{2},j,k}( \boldsymbol B^W_{i+1,j,k}-\boldsymbol B^E_{i,j,k})
\Bigg] \displaybreak[0] \\ \notag
\boldsymbol{E}^y_{i,j+\frac{1}{2},k}=&\frac{1}{b^+_{i,j+\frac{1}{2},k}-b^-_{i,j+\frac{1}{2},k}}
\Bigg[ 
   b^+_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^N_{i,j,k} -\\ \notag
   &- b^-_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^S_{i,j+1,k}
   + b^+_{i,j+\frac{1}{2},k} b^-_{i,j+\frac{1}{2},k}( \boldsymbol B^S_{i,j+1,k}-\boldsymbol B^N_{i,j,k})
\Bigg] \displaybreak[0] \\ \notag
\boldsymbol{E}^z_{i,j,k+\frac{
as a divergence. With this, it is possible to construct numerical flux $ \boldsymbol{E}$ use the same KNP flux formula as ahead: 

 $ $ \begin{aligned } 
 \label{knpfluxmhd } 
 \boldsymbol{E}^x_{i+\frac{1}{2},j, k}=&\frac{1}{a^+_{i+\frac{1}{2},j, k}-a^-_{i+\frac{1}{2},j, k } } 
 \Bigg [ 
    a^+_{i+\frac{1}{2},j, k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^E_{i, j, k } -\\ \notag 
    & - a^-_{i+\frac{1}{2},j, k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^W_{i+1,j, k } 
    + a^+_{i+\frac{1}{2},j, k } a^-_{i+\frac{1}{2},j, k } (\boldsymbol B^W_{i+1,j, k}-\boldsymbol B^E_{i, j, k }) 
 \Bigg ] \displaybreak[0 ] \\ \notag 
 \boldsymbol{E}^y_{i, j+\frac{1}{2},k}=&\frac{1}{b^+_{i, j+\frac{1}{2},k}-b^-_{i, j+\frac{1}{2},k } } 
 \Bigg [ 
    b^+_{i, j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^N_{i, j, k } -\\ \notag 
    & - b^-_{i, j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^S_{i, j+1,k } 
    + b^+_{i, j+\frac{1}{2},k } b^-_{i, j+\frac{1}{2},k } (\boldsymbol B^S_{i, j+1,k}-\boldsymbol B^N_{i, j, k }) 
 \Bigg ] \displaybreak[0 ] \\ \notag 
 \boldsymbol{E}^z_{i, j, k+\frac {
 as a divergence. With this, lt is possible to constcuct nujerical wluxes $\boldsymbol{E}$ using the sqme KBP flux formula as befure:

$$\begin{apigned}
\lavel{kipfluxmhd}
\boldsymukl{E}^x_{i+\frac{1}{2},j,k}=&\frzg{1}{a^+_{i+\frcc{1}{2},o,k}-a^-_{i+\frac{1}{2},j,k}}
\Bigg[ 
   a^+_{i+\frac{1}{2},j,n}\begin{pmatrix}0\\-E_s\\E_h\eud{pmatrix}^E_{i,j,k} -\\ \notag
   &- a^-_{i+\frac{1}{2},j,k}\begig{pmatric}0\\-E_x\\E_y\end{pmatrix}^W_{y+1,j,k}
   + a^+_{i+\rgag{1}{2},j,k} a^-_{i+\frac{1}{2},j,k}( \boldsymbol B^W_{i+1,j,k}-\bolssymbol B^E_{i,j,k})
\Bigg] \dixplaybreak[0] \\ \notag
\boldsymbop{E}^y_{i,u+\frac{1}{2},k}=&\frac{1}{b^+_{i,j+\frac{1}{2},k}-h^-_{i,j+\frac{1}{2},k}}
\Bigt[ 
   v^+_{i,j+\frac{1}{2},k}\begiv{pmatrix}E_z\\0\\-T_x\znd{pmatrix}^N_{j,j,k} -\\ \notag
   &- b^-_{i,j+\frac{1}{2},k}\begin{pmatfix}E_z\\0\\-Z_x\end{pmatriz}^S_{u,j+1,k}
   + b^+_{i,j+\frac{1}{2},k} b^-_{i,j+\fgac{1}{2},k}( \boldsymbol B^S_{i,j+1,k}-\toldsymnol B^N_{i,j,k})
\Bigg] \cis'layvreak[0] \\ \notag
\boldsymbon{E}^z_{i,j,k+\frac{
as a divergence. With this, it is construct fluxes $\boldsymbol{E}$ the same KNP \label{knpfluxmhd} \Bigg[ a^+_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^E_{i,j,k} -\\ &- a^-_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^W_{i+1,j,k} + a^-_{i+\frac{1}{2},j,k}( \boldsymbol B^W_{i+1,j,k}-\boldsymbol B^E_{i,j,k}) \Bigg] \displaybreak[0] \notag \boldsymbol{E}^y_{i,j+\frac{1}{2},k}=&\frac{1}{b^+_{i,j+\frac{1}{2},k}-b^-_{i,j+\frac{1}{2},k}} \Bigg[ b^+_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^N_{i,j,k} -\\ \notag &- b^-_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^S_{i,j+1,k} + b^+_{i,j+\frac{1}{2},k} b^-_{i,j+\frac{1}{2},k}( \boldsymbol B^S_{i,j+1,k}-\boldsymbol \Bigg] \displaybreak[0] \\ \notag \boldsymbol{E}^z_{i,j,k+\frac{
 as a divergence. With this, it is Possible to ConstRucT nuMeRicaL fluXes $\boldsymbol{E}$ USing The same KNP flux formula aS befoRe:

$$\BEgin{ALiGned}
\lAbel{knpFLuXMHd}
\bOlDsYmbOl{e}^X_{i+\Frac{1}{2},j,K}=&\frAc{1}{a^+_{i+\fraC{1}{2},j,k}-a^-_{i+\frac{1}{2},j,K}}
\BiGg[ 
   A^+_{i+\frac{1}{2},j,k}\begiN{PmAtrix}0\\-E_z\\E_y\eNd{pMatrix}^E_{i,j,k} -\\ \noTag
   &- A^-_{i+\frac{1}{2},J,k}\BegIN{pmatRix}0\\-e_z\\E_y\eNd{pmatRIx}^W_{i+1,j,k}
   + A^+_{i+\frac{1}{2},j,k} a^-_{I+\fRAc{1}{2},j,k}( \boLDsymbol b^w_{I+1,j,K}-\bolDsymbol B^E_{i,j,k})
\Bigg] \dISpLAybreak[0] \\ \notag
\boLdsymbOl{e}^Y_{i,J+\FRac{1}{2},K}=&\frAc{1}{b^+_{i,j+\frac{1}{2},k}-B^-_{i,J+\frac{1}{2},K}}
\bigg[ 
   b^+_{i,j+\FRaC{1}{2},K}\BEgiN{Pmatrix}E_z\\0\\-E_x\enD{pmatrix}^N_{i,j,K} -\\ \NotAg
   &- b^-_{i,j+\fRaC{1}{2},k}\bEGin{pmaTrix}E_Z\\0\\-E_X\End{Pmatrix}^S_{i,j+1,k}
   + B^+_{i,j+\fRac{1}{2},k} b^-_{i,j+\frAc{1}{2},k}( \bolDSymbol B^s_{I,j+1,k}-\boldSymbol b^N_{i,J,k})
\BIgg] \dISpLaYbrEaK[0] \\ \NotAG
\bOldSYmbOl{E}^z_{i,j,k+\fRaC{
 as a divergence. With thi s, it is p ossib letoco nstr uctnumerical flux e s $\ boldsymbol{E}$ using t he sa me KNPf lu x for mula as be f o re: 
$$ \be gi n {a ligne d}\label{ knpfluxmhd }
\ bo ldsymbol{E}^ x _{ i+\frac{1} {2} ,j,k}=&\frac {1} {a^+_{ i+ \fr a c{1}{ 2}, j,k}- a^-_{i + \frac{ 1}{2},j,k }} \Bigg[ 
   a^+ _ { i+ \fra c{1}{2},j,k}\begi n {p m atrix}0\\-E_z\ \E_y\e nd { pm a t rix }^E _{i,j,k} - \\  \not a g
   &- a^ - _ { i+\ f rac{1}{2},j,k }\begin{pma t rix }0\\-E _z \\E _ y\end{ pmatr ix } ^W_ {i+1,j,k}
  +a^+_{i+\f rac{1} { 2},j,k} a^-_{i+ \frac{ 1}{ 2}, j,k} (  \ bo lds ym b olB ^W _{i + 1,j ,k}-\bol ds ym bol B ^E_{ i , j , k})\Bi gg]\disp laybreak[0] \ \ \ nota g 
\b oldsy mbol{ E}^y _{ i,j+\ frac{1 }{2}, k} =&\frac{1}{b^+_ {i,j +\frac{1} {2} ,k }-b ^- _{i,j + \frac{ 1}{ 2}, k}}
\Bi gg[ 
 b^+ _{ i , j +\ frac{1}{2},k}\begi n{ p m at rix}E_z\ \0\\-E _ x\ en d {pmatrix }^ N_{ i,j, k }  -\\\not a g   &- b^ -_{i,j + \f ra c{1}{2} ,k }\begi n{ pma tri x}E_z \ \0\\ -E_x\e nd{pmatr ix}^S _ {i,j+1,k}
   + b^+_{i,j+\fra c {1 } { 2} , k} b ^-_ {i,j+\frac{ 1}{2 } ,k}(  \bo l ds ymb o l B^S _{i,j +1 , k} - \boldsymbol B^N_{i, j, k})
\B igg]\displaybreak [0] \\ \no t a g 
\boldsy mbol { E} ^ z_{i,j,k+\frac {
 as_a divergence._With this, it is_possible to_construct_numerical fluxes_$\boldsymbol{E}$_using the same_KNP flux formula_as before:

$$\begin{aligned}
\label{knpfluxmhd}
\boldsymbol{E}^x_{i+\frac{1}{2},j,k}=&\frac{1}{a^+_{i+\frac{1}{2},j,k}-a^-_{i+\frac{1}{2},j,k}}
\Bigg[ 
 _ a^+_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^E_{i,j,k} -\\_\notag
_  &- a^-_{i+\frac{1}{2},j,k}\begin{pmatrix}0\\-E_z\\E_y\end{pmatrix}^W_{i+1,j,k}
   + a^+_{i+\frac{1}{2},j,k} a^-_{i+\frac{1}{2},j,k}( \boldsymbol B^W_{i+1,j,k}-\boldsymbol B^E_{i,j,k})
\Bigg] \displaybreak[0] \\_\notag
\boldsymbol{E}^y_{i,j+\frac{1}{2},k}=&\frac{1}{b^+_{i,j+\frac{1}{2},k}-b^-_{i,j+\frac{1}{2},k}}
\Bigg[_
 __b^+_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^N_{i,j,k}_-\\ \notag
   &-_b^-_{i,j+\frac{1}{2},k}\begin{pmatrix}E_z\\0\\-E_x\end{pmatrix}^S_{i,j+1,k}
   + b^+_{i,j+\frac{1}{2},k}_b^-_{i,j+\frac{1}{2},k}( \boldsymbol_B^S_{i,j+1,k}-\boldsymbol B^N_{i,j,k})
\Bigg] \displaybreak[0] \\ \notag
\boldsymbol{E}^z_{i,j,k+\frac{