AUSM⁺-up Flux Scheme

The wellbore simulator uses the AUSM⁺-up (Advection Upstream Splitting Method) flux scheme for computing numerical fluxes at cell interfaces. Based on Liou (2006) with area-weighted extensions from Sacconi & Mahgerefteh (2020).

Interface Mach Number

The numerical flux at each cell interface is determined by a split Mach number. Left and right states each contribute through the $\mathcal{M}_4^{\pm}$ polynomials. A pressure-difference dissipation term scaled by $K_\mathrm{p}$ improves stability at low Mach numbers.

$$M_{i+\frac{1}{2}} = \mathcal{M}_4^{+} (M_\mathrm{L}) \;+\; \mathcal{M}_4^{-} (M_\mathrm{R}) \;-\; K_\mathrm{p} \, \max\bigl(1 - \bar{M}^2,\; 0\bigr) \, \frac{p_\mathrm{R} - p_\mathrm{L}}{\bar{\rho}\,\bar{a}^2} \tag{1}$$

Mach Splitting Polynomials

The fourth-order Mach polynomials smoothly split the interface Mach number into left- and right-running contributions. For supersonic flow ($|M| \geq 1$) they reduce to simple upwinding. The subsonic branch uses $\beta$ to control blending:

$$\mathcal{M}_4^{\pm}(M) = \begin{cases} \dfrac{1}{2}(M \pm |M|) & |M| \geq 1 \\[8pt] \pm\,\dfrac{1}{4}(M \pm 1)^2 \Bigl[1 + 16\,\beta\;\dfrac{1}{4}(M \mp 1)^2\Bigr] & |M| < 1 \end{cases} \tag{2}$$

Interface Pressure

The interface pressure is reconstructed using $\mathcal{P}_5^{\pm}$ polynomials with a velocity-difference dissipation term scaled by $K_\mathrm{u}$:

$$p_{i+\frac{1}{2}} = \mathcal{P}_5^{+}\,p_\mathrm{L} \;+\; \mathcal{P}_5^{-}\,p_\mathrm{R} \;-\; K_\mathrm{u}\;\mathcal{P}_5^{+}\,\mathcal{P}_5^{-}\,(\rho_\mathrm{L} + \rho_\mathrm{R})\,\bar{a}\;\frac{A_\mathrm{R}\,u_\mathrm{R} - A_\mathrm{L}\,u_\mathrm{L}}{\max(A_\mathrm{L},\,A_\mathrm{R})} \tag{3}$$

Pressure Splitting Polynomials

$$\mathcal{P}_5^{\pm}(M) = \begin{cases} \dfrac{\mathcal{M}_1^{\pm}}{M} & |M| \geq 1 \\[8pt] \pm\,\mathcal{M}_2^{\pm} \Bigl(2 \mp M - 16\,\alpha\;M\,\mathcal{M}_2^{\mp}\Bigr) & |M| < 1 \end{cases} \tag{4}$$

Mass Flux

The final convective mass flux uses upwind selection based on the interface Mach number:

$$\dot{m}_{i+\frac{1}{2}} = \bar{a}\left[\frac{(\rho A)_\mathrm{L}}{2}\bigl(M_{i+\frac{1}{2}} + |M_{i+\frac{1}{2}}|\bigr) \;+\; \frac{(\rho A)_\mathrm{R}}{2}\bigl(M_{i+\frac{1}{2}} - |M_{i+\frac{1}{2}}|\bigr)\right] \tag{5}$$

References

• Liou, M.-S. (2006). A sequel to AUSM, Part II: AUSM⁺-up for all speeds. Journal of Computational Physics, 214(1):137–170. doi:10.1016/j.jcp.2005.09.020
• Sacconi, A. & Mahgerefteh, H. (2020). Modelling start-up injection of CO₂ into highly-depleted gas fields. Energy, 191, 116530. doi:10.1016/j.energy.2019.116530