In van Opstal et al. (Comput Mech 50:779---788, 2012) airbag inflation simulations were performed where the flow was approximated by Stokes flow. Inside the intricately folded initial geometry the Stokes assumption is argued to hold. This linearity assumption leads to a boundary-integral representation, the key to bypassing mesh generation and remeshing. It therefore enables very large displacements with near-contact. However, such a coarse assumption cannot hold throughout the domain, where it breaks down one needs to revert to the original model. The present work formalizes this idea. A model adaptive approach is proposed, in which the coarse model (a Stokes boundary-integral equation) is locally replaced by the original high-fidelity model (Navier---Stokes) based on a-posteriori estimates of the error in a quantity of interest. This adaptive modeling framework aims at taking away the burden and heuristics of manually partitioning the domain while providing new insight into the physics. We elucidate how challenges pertaining to model disparity can be addressed. Essentially, the solution in the interior of the coarse model domain is reconstructed as a post-processing step. We furthermore present a two-dimensional numerical experiments to show that the error estimator is reliable.