Eulerian models of multi-fluid flows consist in dividing the probability density function (PDF) of the flow in several groups (called fluids) which undergo separated but coupled transport and evolution. These models involve many features which challenge numerical simulation schemes. Noticeably, non-standard multi-modal PDF, compressibility with contrasted equations of state (EOS), and multi-way couplings at large volume fractions strain the schemes in terms of PDF resolution, entropy consistency, and coupling stiffness. Stiff couplings can notably appear in mass and momentum exchanges and in pressure works. For instance, velocity relaxation (momentum exchange) can be stiff due to large drag forces. Stiffness in the pressure work can appear in energy equations when fluids are of highly contrasted compressibility or volume fractions: global volume changes of the mixture then impact fluids very unevenly. An extreme example of high practical impact is the case of air bubbles in water at small volume fractions. To preserve robustness, numerical schemes generally involve implicit or explicit discretizations for respectively negative or positive stiff terms: for instance relaxation terms such as stiff drag require implicit schemes. Now, the sign of pressure work in the internal energy equations may change depending on expansion and compression phases. In numerous published studies, this constraint is mollified by introducing a finite-time relaxation equation of pressure imbalance. This comes at the price of some serious detrimental and well-known side-effects. For a stiff pressure equilibration (or instantaneous relaxation), schemes based on implicit pressure (possibly exponential) appear much more convenient in order to preserve robustness with thermodynamic consistency in all conditions (for instance forbidding negative densities): high contrasts of volume fractions, densities, and EOS become accessible. We illustrate this by successfully simulating not only the most strenuous mixed water-air and air-water shock tubes published in the literature, but also the two theoretical extremes of expansion to vacuum and shock from vanishing pressure. The simulations are carried out with inter-fluid pressure equilibration embedded in a previously developed Euler-Euler explicit compressible scheme [Vazquez-Gonzalez et al. IJMF 2020] without complex multi-dimensional Newton-Raphson iterations and under normal CFL conditions [Heulhard de Montigny & Llor IJMF 2022].