arXiv:2610.36349v1 Announce Type: cross Abstract: Nearly incompressible materials arise in soft tissues, elastomers, gels, and other soft solids. Yet, numerical modeling remains challenging because displacement-based formulations can exhibit volumetric locking near incompressibility. Mixed displacement-pressure formulations address this difficulty but require compatible approximation spaces satisfying the inf-sup condition. In mixed formulations, discontinuous pressure spaces are attractive since pressure degrees of freedom can be condensed locally. In high-order isogeometric analysis, however, their Galerkin stability requires pressure-grid coarsening that increases with polynomial degree. Here, we develop a variational multiscale (VMS) stabilized mixed isogeometric formulation that retains this locality while substantially relaxing the restriction. The key observation is that VMS stabilization controls nonconstant elementwise pressure modes and is unconditionally stable with respect to its parameter, provided the elementwise constant pressure component satisfies the inf-sup condition. Consequently, for maximally continuous displacement splines, a pressure grid coarsened by a factor of two satisfies the constant-mode inf-sup criterion independently of polynomial degree. Local condensation then yields a displacement-only system, facilitating use in existing displacement solvers. The formulation is further extended to immersed isogeometric analysis with extended B-splines, pressure-cell aggregation, subdivision-based cut-cell quadrature, and weak enforcement of essential boundary conditions. Discrete inf-sup studies and convergence analyses in the Kirsch and Cook's membrane benchmarks demonstrate stability and show that the immersed treatment preserves mixed stability properties. The proposed framework therefore provides a local, high-order, and robust approach for nearly incompressible problems in mesh-conforming and immersed settings.