Abstract: In some aspects, a number of subsystem models is accessed by a computer system. Each subsystem model represents dynamic attributes of a distinct physical subsystem in a subterranean region. At least one of the number of subsystem models represents dynamic attributes of a mechanical subsystem in the subterranean region. A discrete fracture network (DFN) model representing a fracture network in the subterranean region is accessed at the computer system. The DFN model includes junction models. Each junction model represents interactions between a respective set of subsystem models associated with the junction model. Junction variables of the junction model can be defined based on dynamic attributes of the respective set of subsystem models associated with the junction model. A stimulation treatment for the subterranean region can be simulated by operating the DFN model including the junction variables.

Abstract: In some aspects, a flow model represents a flow path for well system fluid in a subterranean region. A first system of equations from a finite element (FE) discretization of a flow model is accessed by a computer system. A second system of equations that imposes a constraint on errors between an exact solution and an FE solution to the flow model is accessed by the computer system. The FE solution includes nodes representing locations along the flow path for well system fluid in the subterranean region. Locations of the nodes can be computed by operations of the computer system in an iterative manner by computing locations of new nodes based on locations of existing nodes and refining the locations of the new nodes and the locations of the existing nodes based on the first and second systems of equations.

Abstract: In some aspects, a number of subsystem models is accessed by a computer system. Each subsystem model represents dynamic attributes of a distinct physical subsystem in a subterranean region. At least one of the number of subsystem models represents dynamic attributes of a mechanical subsystem in the subterranean region. A discrete fracture network (DFN) model representing a fracture network in the subterranean region is accessed at the computer system. The DFN model includes junction models. Each junction model represents interactions between a respective set of subsystem models associated with the junction model. Junction variables of the junction model can be defined based on dynamic attributes of the respective set of subsystem models associated with the junction model. A stimulation treatment for the subterranean region can be simulated by operating the DFN model including the junction variables.

Abstract: In some aspects, a flow model represents a flow path for well system fluid in a subterranean region. A first system of equations from a finite element (FE) discretization of a flow model is accessed by a computer system. A second system of equations that imposes a constraint on errors between an exact solution and an FE solution to the flow model is accessed by the computer system. The FE solution includes nodes representing locations along the flow path for well system fluid in the subterranean region. Locations of the nodes can be computed by operations of the computer system in an iterative manner by computing locations of new nodes based on locations of existing nodes and refining the locations of the new nodes and the locations of the existing nodes based on the first and second systems of equations.