We adapt the viscosity method introduced by Rivière in [37] to the free boundary case. Namely, given a compact oriented surface Σ, possibly with boundary, a closed ambient Riemannian manifold (M m , g) and a closed embedded submanifold N n ⊂ M, we study the asymptotic behavior of (almost) critical maps Φ for the functionalon immersions Φ : Σ → M with the constraint Φ(∂Σ) ⊆ N , as σ → 0, assuming an upper bound for the area and a suitable entropy condition.As a consequence, given any collection F of compact subsets of the space of smooth immersions (Σ, ∂Σ) → (M, N ), assuming F to be stable under isotopies of this space we show that the min-max valueis the sum of the areas of finitely many branched minimal immersions Φ (i) : Σ (i) → M with ∂ν Φ (i) ⊥ T N along ∂Σ (i) , whose (connected) domains Σ (i) can be different from Σ but cannot have a more complicated topology.We adopt a point of view which exploits extensively the diffeomorphism invariance of Eσ and, along the way, we simplify several arguments from the original work [37]. Some parts generalize to closed higher-dimensional domains, for which we get a rectifiable stationary varifold in the limit.