In this work, we present an upscaled model for mixed dimensional coupled flow problem in fractured porous media. We consider both embedded and discrete fracture models (EFM and DFM) as fine scale models which contain coupled system of equations. For fine grid discretization, we use a conservative finite-volume approximation. We construct an upscaled model using the non-local multicontinuum (NLMC) method for the coupled system. The proposed upscaled model is based on a set of simplified multiscale basis functions for the auxiliary space and a constraint energy minimization principle for the construction of multiscale basis functions. Using the constructed NLMC-multiscale basis functions, we obtain an accurate coarse grid upscaled model. We present numerical results for both fine-grid models and upscaled coarse-grid models using our NLMC method. We consider model problems with (1) discrete fracture fine grid model with low and high permeable fractures;(2) embedded fine grid model for two types of geometries with differnet fracture networks and (3) embedded fracture fine grid model with heterogeneous permeability. The simulations using the upscaled model provide very accurate solutions with significant reduction in the dimension of the problem.

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high

Mortar methods are widely used techniques for discretizations of partial differential equations and preconditioners for the algebraic systems resulting from the discretizations. For problems with high contrast and multiple scales, the standard mortar spaces are not robust, and some enrichments are necessary in order to obtain an efficient and robust mortar space. In this paper, we consider a class of flow problems in high contrast heterogeneous media, and develop a systematic approach to obtain an enriched multiscale mortar space. Our approach is based on the constructions of local multiscale basis functions. The multiscale basis functions are constructed from local problems by following the framework of the Generalized Multiscale Finite Element Method (GMsFEM). In particular, we first create a local snapshot space. Then we select the dominated modes within the snapshot space using an appropriate Proper Orthogonal Decomposition (POD) technique. These multiscale basis functions show better accuracy than polynomial basis for multiscale problems. Using the proposed multiscale mortar space, we will construct a multiscale finite element method to solve the flow problem on a coarse grid and a preconditioning technique for the fine scale discretization of the flow problem. In particular, we develop a multiscale mortar mixed finite element method using the mortar space. In addition, we will design a two-level additive preconditioner and a two-level hybrid preconditioner based on the proposed mortar space for the iterative method applied to the fine scale discretization of the flow problem. We present several numerical examples to demonstrate

The objective of this paper is to design novel multi-layer neural network architectures for multiscale simulations of flows taking into account the observed data and physical modeling concepts. Our approaches use deep learning concepts combined with local multiscale model reduction methodologies to predict flow dynamics. Using reduced-order model concepts is important for constructing robust deep learning architectures since the reduced-order models provide fewer degrees of freedom. Flow dynamics can be thought of as multi-layer networks. More precisely, the solution (eg, pressures and saturations) at the time instant n+ 1 depends on the solution at the time instant n+ 1 and input parameters, such as permeability fields, forcing terms, and initial conditions. One can regard the solution as a multi-layer network, where each layer, in general, is a nonlinear forward map and the number of layers relates to the internal time steps. We will rely on rigorous model reduction concepts to define unknowns and connections for each layer. In each layer, our reduced-order models will provide a forward map, which will be modified (" trained") using available data. It is critical to use reduced-order models for this purpose, which will identify the regions of influence and the appropriate number of variables. Because of the lack of available data, the training will be supplemented with computational data as needed and the interpolation between data-rich and data-deficient models. We will also use deep learning algorithms to train the elements of the reduced model discrete system. We will present main ingredients of our approach and numerical results. Numerical

We develop a PetrovGalerkin stabilization method for multiscale convectiondiffusion transport systems. Existing stabilization techniques add a limited number of degrees of freedom in the form of bubble functions or a modified diffusion, which may not be sufficient to stabilize multiscale systems. We seek a local reduced-order model for this kind of multiscale transport problems and thus, develop a systematic approach for finding reduced-order approximations of the solution. We start from a PetrovGalerkin framework using optimal weighting functions. We introduce an auxiliary variable to a mixed formulation of the problem. The auxiliary variable stands for the optimal weighting function. The problem reduces to finding a test space (a dimensionally reduced space for this auxiliary variable), which guarantees that the error in the primal variable (representing the solution) is close to the projection error of the full