1. A method for generating a discretized physics-based model of a subsurface region having a shape, the method comprising: selecting a pre-solved model related to at least a majority by volume of the subsurface region to be modeled by the discretized physics-based model; applying a solution of potential field lines to the pre-solved model to create a mesh the pre-solved model, the mesh including a first set of domain boundaries and the mesh and the pre-solved model forming an intermediate model; defining the shape of the subsurface region to be modeled by the discretized physics-based model, the shape including a second set of domain boundaries representing subsurface features not reflected in the intermediate model such that the first set of domain boundaries of the intermediate model is not fully aligned with the shape of the subsurface region represented by the second set of domain boundaries; and using a computer processor accessing non-transitory computer processor-readable medium to transform the domain boundaries of the intermediate model to the shape of the subsurface region represented by the second set of domain boundaries of the subsurface region by stretching the mesh of the intermediate model to the shape of the subsurface region defined by the second set of domain boundaries to generate the discretized physics-based model, wherein at least a majority by number of intersections of orthogonal elements within the intermediate model retain orthogonality within 75-105 degrees in the discretized physics-based model.