Goaloriented adaptivity for GMsFEM
Abstract
In this paper we develop two goaloriented adaptive strategies for a posteriori error estimation within the generalized multiscale finite element framework. In this methodology, one seeks to determine the number of multiscale basis functions adaptively for each coarse region to efficiently reduce the error in the goal functional. Our first error estimator uses a residual based strategy where local indicators on each coarse neighborhood are the product of local indicators for the primal and dual problems, respectively. In the second approach, viewed as the multiscale extension of the dual weighted residual method (DWR), the error indicators are computed as the pairing of the local residual of the primal problem weighed by a projection into the primal space of the dual solution from an enriched space, over each coarse neighborhood. In both of these strategies, the goaloriented indicators are then used in place of a standard residualbased indicator to mark coarse neighborhoods of the mesh for further enrichment in the form of additional multiscale basis functions. The method is demonstrated on highcontrast problems with heterogeneous multiscale coefficients, and is seen to outperform the standard residual based strategy with respect to efficient reduction of error in the goal function.
1 Introduction
Many practical problems are multiscale in nature, including flow in porous media, seismic wave propagation, and physical processes in perforated media. These problems are described by partial differential equations (PDEs) with potentially high contrast multiscale coefficients. Direct computation of high resolution discrete solutions to these problems can be very expensive. Typically, some type of model reduction techniques are used to solve multiscale problems. Established techniques include numerical homogenization methods [18, 10, 9] and multiscale methods [32, 24, 19, 16, 8, 29, 33, 11, 17, 2, 1, 22]. In numerical homogenization methods, the upscaled media properties are computed over coarse grid blocks, each of which is much larger than a characteristic length scale. In multiscale methods, local multiscale basis functions determined by local finescale problems are constructed in each element of the coarse grid. Generally, one uses a few multiscale basis functions in each coarse element to approximate the global solution by solving a coarse mesh problem over the entire domain. Multiscale basis functions are constructed in an offline step before the coarse mesh problem is solved, after which some type of adaptivity is needed to choose multiscale basis function appropriately.
In this paper we will use multiscale methods, where multiscale basis functions are constructed in each coarse region, as illustrated schematically in Figure 1. To be more specific, we consider a multiscale problem
where and seek the solution in the form
In each coarse region , we construct a set of multiscale basis functions , . These multiscale basis functions, as described below, will be constructed in the offline stage using the generalized multiscale finite element method (GMsFEM) and represent the local heterogeneities of the solution space. In earlier works on the multiscale finite element method (MsFEM) [32], the authors sought one multiscale basis function per coarse element. However as it was later argued, one may need additional basis functions over each coarse element for a sufficiently high fidelity approximation. It is further shown that the optimal number of multiscale basis functions in each region depends on the heterogeneities in the solution space. Typically, some adaptive criteria based on a posteriori error estimation is used to determine how many basis functions to choose in each coarse region . For exampe, in [16], the authors develop an error indicator based on the norm of the residual to determine the number of basis functions to add in each region over each adaptive iteration. Adaptive multiscale methods follow traditional adaptivity concepts as in [41, 5, 7, 37, 30, 35, 40]; however, the multiscale indicators contain information about local heterogeneities. Earlier approaches to goaloriented adaptive methods for multiscale problems include numerical regularization and numerical homogenization with adaptive mesh refinement as in [38, 34]. To the authors’ knowledge, the current presentation is the first to develop a goaloriented enrichment strategy within the general GMsFEM framework.
For many practical problems, one is interested in approximating some function of the solution, known as the quantity of interest, rather than the solution itself. Examples include an average or weighted average of the solution over a particular subdomain, or some localized solution response. In these cases, goaloriented adaptive methods yield a more efficient approximation than standard adaptivity, as the enrichment of degrees of freedom is focused on the local improvement of the quantity of interest rather than across the entire solution [39, 4, 36, 26, 27, 28, 3, 6, 31]. In this paper, we study goaloriented adaptivity for multiscale methods, and in particular the design of error indicators to drive the adaptive enrichment based on the goal function. In multiscale methods, goaloriented adaptivity can play an important role in the efficient approximation of the quantity of interest as heterogeneities in the coefficients may require standard adaptive methods to add degrees of freedom in regions with limited influence on the goal function. In this paper, we develop a goaloriented approach for multiscale methods within GMsFEM framework. In the proposed approach, we increase the accuracy of the approximation by enriching the space rather than refining the mesh by choosing multiscale basis functions computed in the offline stage.
For multiscale basis construction, we use GMsFEM. The construction of multiscale basis functions uses local snapshot spaces and requires solving local spectral problems over each coarse element. The local snapshot functions represents the solution space in each coarse region, and they can include all possible local finegrid functions or harmonic functions. In the snapshot space, we perform a local spectral decomposition and select multiscale basis functions which correspond to the dominant eigenvalues. The multiscale basis functions are constructed by multiplying the dominant eigenmodes by a partition of unity function, e.g., multiscale partition of unity function. In [16], we developed an adaptive approach and developed a posteriori error indicators, which include the information from local spectral problems, e.g., the value of the eigenvalue corresponding to the first eigenvector not included in the coarse space. We derived error estimates and presented numerical results which demonstrate the improved efficiency of the adaptive approach, guided these indicators. For goaloriented problems, we now design goaloriented error indicators, which are different from those developed earlier [16] for multiscale problems, by the additional consideration of a dual problem to direct the adaptivity towards the approximation of the quantity of interest.
In this paper we develop two goaloriented adaptive strategies for a posteriori error estimation. Our first error estimator uses an idea similar to a standard residual based adaptive method, and can be seen as the multiscale extension of the adaptive method presented in [6]. In this case the elementwise indicator is formed by the product of local residual indicators for the primal and dual problems, respectively. In the second approach, viewed as the multiscale extenstion of the dual weighted residual method (DWR), the error indicators are computed as the pairing of the local residual of the primal problem weighed by a projection into the primal space of the dual solution from an enriched space. In both of these strategies, the goaloriented indicators are then used in place of a standard residualbased indicator to mark coarse elements of the mesh for further enrichment in the form of additional multiscale basis functions.
The remainder of the paper is organized as follows. In Section 2 we present and overview of multiscale methods and adaptivity. In Section 3 we give a detailed description of the construction of the multiscale basis functions, review some results on residualbased adaptivity for GMSFEM. In Section 4 we introduce two goaloriented a posteriori error indicators and present an algorithm for goaloriented adaptivity. Finally in Section 5 we present numerical results demonstrating the efficiency of the proposed goaloriented error indicators.
2 Overview of Concepts
In this paper, we consider second order multiscale elliptic problems of the form
(1) 
where is the computational domain, and is a scalar valued heterogeneous coefficient with multiple scales and high contrast. The problem (1) can be solved by many classical numerical techniques, such as the conforming finite element method, but with extremely high computational complexity due to the fact that a very fine mesh is necessary to resolve the multiscale nature of the solution. Thus, some multiscale model reductions are needed to compute an accurate solution efficiently. In the following, we give a brief overview of GMsFEM and its basis enrichment techniques as applied to problem (1). Let be the true solution satisfying
(2) 
where , and . Define the energy norm on by .
To introduce the GMsFEM for the problem (1), we first give the notion of fine and coarse grids. We let be a standard conforming triangulation of the computational domain into finite elements, which can be triangular, rectangular or some other polygons. We refer to this partition as the coarse grid. Subordinate to the coarse grid, we define the fine grid partition, denoted by , by refining each coarse element into a connected union of fine grid blocks. We assume the above refinement is performed such that is a conforming partition of . We let be the number of interior coarse grid nodes, and let be the set of coarse grid nodes or vertices of the coarse mesh . Moreover, we define the coarse neighborhood of the node by
(3) 
which is the union of all coarse elements which have the node as a vertex. See Figure 1 for an illustration of the coarse elements and coarse neighborhoods within the coarse grid. We emphasize the use of to denote a coarse neighborhood, and to denote a coarse element throughout the paper.
Next, we briefly overview the continuous Galerkin (CG) formulation of GMsFEM, a generalization of the classical MsFEM [32]. For each coarse node , we define a set of basis functions supported on the neighborhood . We denote the th basis function supported on the coarse neighborhood by , We remark that in the GMsFEM, we will use multiple basis functions per coarse neighborhood, and the index represents the local numbering of these basis functions. These multiscale basis functions are constructed from a local snapshot space and a local spectral decomposition defined on that snapshot space. The snapshot space contains a collection of many basis functions that can be used to capture most of the fine features of the solution, and the multiscale basis functions are constructed by selecting the dominant modes of a local spectral problem. Using the these multiscale basis functions, the CG solution is represented as . Once the basis functions are identified, the CG global coupling is given through the variational form
(4) 
where is the space spanned by the basis functions , and is the usual bilinear form corresponding to (1). We remark that one can use other formulations, such as the discontinuous Galerkin formulation (see e.g., [12, 14, 20]), the mixed formulation (see e.g., [13, 8]) or the hybridized discontinuous Galerkin formulation (see e.g., [25]) to couple the multiscale basis functions.
In using GMsFEM, it is desirable to determine the number of basis functions per coarse neighhorbood adaptively based on the heterogeneities of the coefficient in order to obtain an efficient representation of the solution. In [16], a residual based a posteriori error indicator is derived and an adaptive basis enrichment algorithm is developed under the CG formulation. In particular, it is shown that
where is the residual of the solution on the coarse neighborhood . Thus, local residuals of the multiscale solution give indicators to the error of the solution in the energy norm, and one can add basis functions to the coarse neighborhoods when the residuals are large. Convergence of this adaptive basis enrichment algorithm is also shown in [16]. On the other hand, for some applications one needs to adaptively construct new basis functions in the online stage in order to capture distant effects. In [15], such online adaptivity is proposed and mathematically analyzed. More precisely, when the local residual is large, one can construct a new basis function in the online stage by solving
where is the restriction of in with zero trace on . Numerical results in [15] show that a couple of these online basis functions can help to reduce the error by a large amount.
The adaptivity procedures discussed above are designed with the aim of reducing the error in the energy norm. In some applications, one may be more interested in reducing error measured by some function of the solution other than a norm. For example, in flow applications, one needs to obtain a good approximation of the pressure in locations where the wells are situated. Therefore, we now consider goaloriented adaptivity within GMsFEM. Specifically, we define a linear functional , referred to as the goal functional. In goaloriented adaptivity, one wants to adaptively enrich the approximation space in order to reduce the goal error defined by . In the construction of goaloriented adaptivity for GMsFEM, we use local indicators based on the solution of a dual problem: finding such that
(5) 
For a primal problem based on bilinear form , the dual form is the formal adjoint of the primal, satisfying , and in the current symmetric case, is identical to the primal. Formally, the primaldual equivalence follows for the solution to the primal problem (2) and the solution to the dual problem (5)
(6) 
Error estimates for the quantity of interest follow from (6) and Galerkin orthogonality with respect to the discrete problems and their respective solutions. Forming error indicators based on both primal and dual problems and these estimates, we add multiscale basis functions to coarse neighborhoods when the values of the corresponding indicators are large. Our numerical examples show that the goaloriented approach performs better than the residual approach for highcontrast problems when the error is measured by the goalfunctional, .
3 The GMsFEM and residualbased adaptivity
In this section, we will give a detailed description of the GMsFEM (see for example [19, 21]) and it’s residual based adaptivity (see for example [16]).
3.1 Local basis functions
We first present the construction of the multiscale basis functions. This construction is performed in the offline stage; that is, basis functions are precomputed before the actual solve of the problem. The construction starts with a snapshot space. This space contains a relatively large set of basis functions which can be used to capture most features of the finescale solution. The next step is to perform a local dimension reduction to obtain a lower dimensional subspace that can still be used to approximate the solution with good accuracy. The local dimension reduction is performed by solving a spectral problem, and the dominant eigenfunctions are used as the multiscale basis functions.
First, we define a snapshot space , where the functions in are supported in . The snapshot space can be the space of all finescale basis functions or the solutions of some local problems with various choices of boundary conditions. For example, we can use the following harmonic extensions to form a snapshot space. Specifically, let , index the set of finegrid vertices that lie on the boundary of each coarse neighborhood, . Define the unit source functions for each . Then construct the snapshot function by solving
(7) 
The snapshot space corresponding to the region , then contains functions
We define the corresponding change of variable matrix
where are considered as the columns of the matrix.
We next determine a set of dominant modes from , and the resulting lower dimensional space is called the offline space . To construct the offline space , we perform a dimension reduction of the space of snapshots using an auxiliary spectral decomposition. The analysis in [23] motivates the following generalized eigenvalue problem for eigenvalues and eigenfunctions in the space of snapshots:
(8) 
where
and
where and denote analogous finescale stiffness and mass matrices as defined by
where is the finescale basis function for . We will give the definition of later on. To generate the offline space we then select the smallest eigenvalues from Equation (8) and form the corresponding eigenvectors in the space of snapshots by setting (for ), where are the coordinates of the vector , and is the number of eigenvectors chosen to span the offline space. We will use the set of local basis functions to form the approximation space in the next section.
3.2 CG formulation
In this section we create an appropriate solution space and the variational formulation for a continuous Galerkin approximation of Equation (1). The idea is to use the basis set , , to form the approximation space, called the offline space, and apply the standard continuous Galerkin formulation. We begin with an initial coarse space , where we recall denotes the number of coarse neighborhoods corresponding to interior coarse nodes. Here, are the standard multiscale partition of unity functions which are supported in and are defined by
(9)  
for all coarse elements , where is a continuous function on which is linear on each edge of . Based on the analysis in [23], the summed, pointwise energy required for the eigenvalue problems (8) is defined as
where denotes the coarse mesh size.
The partition of unity functions are then multiplied by the eigenfunctions to construct the multiscale basis functions
(10) 
where we recall denotes the number of offline eigenvectors that are chosen for each coarse node . We note the construction in Equation (10) yields continuous basis functions due to the multiplication of offline eigenvectors with the initial (continuous) partition of unity. Next, we define the continuous Galerkin spectral multiscale space as
(11) 
Using this offline space, we obtain the GMsFEM as in (4).
3.3 Residualbased adaptivity
In this section, we give a brief review of the adaptive basis enrichment algorithm proposed in [16]. After solving the coarse mesh problem and computing error indicators, the standard Dörfler marking strategy is applied with respect to the neighborhoods indexed by vertices , rather than the elements . The marked neighborhoods are then enriched with additional basis functions. The spectral problem (8) gives a natural ordering of basis functions in the local snapshot space with respect to the eigenvalues, in increasing order. The analysis of [16] then suggests an adaptive procedure to add basis functions based on a local error indicator. The error indicator is defined by a norm based residual, which gives a robust indicator with good performance for cases with high contrast media.
Let be a coarse neighborhood and write . For a given multiscale solution , we define a linear functional on by
(12) 
The norm of is defined as
(13) 
where . The norm gives a measure on how well the solution satisfies the variational problem (2) restricted to . The norm can be obtained by solving an auxiliary problem of finding such that for all and then setting . One can solve this auxiliary problem in the snapshot space to reduce the computational cost. Moreover, in [16] the following a posteriori error bound is proved
(14) 
where is a uniform constant independent of the contrast of , and is the th eigenvalue of the eigenvalue problem (8) for the coarse neighborhood . In particular, is the first, i.e., smallest eigenvalue from the spectral problem (8) for which the corresponding eigenvector is not included in the construction of the offline space.
Using the above error bound (14), a convergent adaptive enrichment algorithm is developed in [16]. We now describe this algorithm. The algorithm is an iterative process, and basis functions are added in each iteration/level based on the magnitudes of local residuals. We use to index the enrichment level and let be the solution space at the th iteration. For each coarse region, let be the number of eigenfunctions used at the enrichment level for the coarse region .
3.4 Adaptive algorithm
The adaptive enrichment algorithm for GMsFEM is now summarized below. Choose a fixed marking parameter . Choose also an initial offline space by specifying a fixed number of basis functions for each coarse neighborhood, and this number is denoted by , for each . Then, generate a sequence of spaces and a sequence of multiscale solutions obtained by solving (4). Specifically, for each , perform the following calculations:

Find the multiscale solution in the current space. That is, find such that
(15) 
Compute the local residual. For each coarse region , compute
where
consistent with (12), and the norm is defined in (13) respectively. Next, reenumerate the coarse neighborhoods so the above local residuals are arranged in decreasing order . That is, in the new enumeration, the coarse neighborhood has the largest residual and the coarse neighborhood has the least residual .
Remark 3.1.
An alternate approach to avoid the complexity of the full sort is the standard binning or heapifying strategy [36]. Let , and consider only that satisfy . Let , and perform a partial sort of the remaining indicators collecting or binning the indices for which , for , where is the smallest integer to satisfy

Find the coarse regions where enrichment is needed. Choose the smallest integer such that
(16) The coarse neighborhoods are then enriched with additional basis functions. If the partial sort of Remark 3.1 is used in place of the sort, elements are marked by emptying the first bin, those indicators with , and then continuing on to the second bin, and so forth until (16) is satisfied. As elements within bins are not sorted, this yields a quasioptimal marked set, i.e., the marked set may not be the set of leastcardinality to satisfy (16) as in the full sort, but it is within a factor of two of the least cardinality.

Enrich the space. For each , add basis function for the region according to the following rule. Let be the positive integer such that is large enough compared with (see Remark 3.2). Then include the eigenfunctions in the construction of the basis functions. The resulting space is denoted as . Mathematically, the space is defined as
where and , with , denote the new basis functions obtained by the eigenfunctions . In addition, we set .
4 Goaloriented adaptivity
In this section, we present a goaloriented adaptive enrichment algorithm for GMsFEM. The goaloriented variant of the adaptive method requires the solution of a dual problem in addition to the primal at each iteration. The indicators are computed with both the primal residual and either the dual residual or a projection of an enriched dual solution into the primal space. These indicators predict which neighborhoods to enrich to increase the quality of the approximation of the quantity of interest. After introduction of the discrete dual problem, the finite dimensional analogue of (5), we propose two error indicators for goaloriented enrichment.
The dual problem plays a vital role in goaloriented adaptivity as the vehicle for introducing the goal functional into the adaptive process. Given a goal functional , we define the discrete dual problem on approximation space as: find such that
(17) 
As in (5), the discrete dual form is identical to the primal for symmetric problems. The dual problem, however, features the goal functional as the source. The discrete dual solution, may now be used to define goaloriented error indicators.
based goaloriented indicator
Our first goaloriented indicator is similar in form to the residual based indicator described in the previous section. To motivate this indicator, we introduce the local bilinear form , and the induced localized energy norm . Let be the solution to (2), be the solution to (4), the solution to (5), and the solution to (17). Using Galerkin orthogonality and the relation between the primal and dual problems, the error in the quantity of interest satisfies
(18) 
Decomposing the global integration into neighborhoods by the partition of unity functions given by (9)
(19) 
where we used the fact that is supported in and the fact that . Instead of using the norm of local residual for the primal problem defined in (13) to be our indicator, (19) suggests using the product of norms of local residuals for the primal and dual problems, posed in the same discrete space, . The local dual residual is defined by
(20) 
where is the solution to (17). Analogous to (13), the norm of is defined as
(21) 
The local version of (14), applied to both primal and dual residuals, namely
(22) 
motivates the local error indicator, , defined as
(23) 
where and are defined in (13) and (21), respectively. Applying (23) and (22) to (19) bounds the error in the goal function by
(24) 
In summary, the goalerror over the global domain is bounded by the product of energy errors of the primal and dual problems, which is in turn bounded by the sum of the indicators given by (23), modified by the partition of unity functions. This upper bounds suggests the adequacy of the indicators in reducing the error. The efficiency and a formal convergence analysis are however not addressed here. As in [6] where a similar indicator is used for refinement, this indicator displays similar behavior to the DWRtype indicator, as shown in the numerical experiments; however, it is more amenable to analysis. This indicator has the added advantage of reduced computational cost as compared to the DWRtype indicator described below, as both primal and dual problems are solved over the same discrete spaces, whereas for the DWRtype method, the dual problem has greater computational complexity than the primal.
DWRtype goaloriented indicator
The next error indicator is similar to DWR error indicator. For the primal problem solved in discrete space , the DWR indicator is motived by the following residual equation. For the solution to (5), the solution to (2), and the solution to (4)
(25) 
where in is arbitrary and the global residual . We let be the component of spanned by the basis functions corresponding to the coarse neighborhood . Localizing (25) by the partition of unity functions ,
(26) 
As the exact solution is unavailable, one generally instead replaces by , a discrete solution from a more enriched space than the primal, essentially neglecting the last term of (26). The function is the component of spanned by the basis functions corresponding to the coarse neighborhood . In standard finite element methods the global residual is then solved elementwise and used as an indicator. By Galerkin orthogonality, then function may be taken as any function in but in practice is taken as the projection of the enriched dual solution into .
In this case, the dual problem is solved in the enriched space, called . The space is obtained by adding more basis functions to each coarse neighborhood. Recalling the construction of basis functions in (10), the enriched space is constructed with more than basis functions per coarse neighborhood, specifically
(27) 
where basis functions are added for each . The span of these basis is our . Let be the solution for the dual problem in , that is, satisfies
(28) 
The DWRtype error estimator is defined as
(29) 
In the above definition, is the component of spanned by the basis functions , , corresponding to the coarse neighborhood . Moreover is the component of spanned by the basis functions , , in the offline space. Comparison with (26) yields an heuristic bound, modulo the error term created by replacing by .
Goaloriented adaptive enrichment algorithm
Choose a fixed marking parameter . Choose also an initial offline space by specifying a fixed number of basis functions for each coarse neighborhood, and this number is denoted by . Then, generate a sequence of spaces and a sequence of multiscale solutions obtained by solving (4). Specifically, for each , perform the following calculations:

Find the multiscale solution in the current space. That is, find such that
(30) 
Find the multiscale dual solution in the current space or an enriched space. That is, find such that
(31) where

Compute the local residual. For each coarse region , compute
consistent with (12) and (20); and the norm is defined in (13) and (21) respectively. Next, reenumerate the coarse neighborhoods so that the above local residuals are arranged in decreasing order . That is, in the new enumeration, the coarse neighborhood has the largest residual . As in Remark 3.1 the full sort of the estimators can be replaced by a partial sort for a marked set of quasioptimal cardinality.

Find the coarse regions where enrichment is needed. Choose the smallest integer such that
(32) The coarse neighborhoods are then enriched with additional basis functions. Alternately, a set of neighborhoods based on the binning strategy for a partial sort can be chosen as described in Step 3 of the Adaptive algorithm 3.4.

Enrich the space. For each , add basis functions for the region according to the following rule. Let be the smallest positive integer such that is large enough compared with . Then include the eigenfunctions in the construction of the basis functions. The resulting space is denoted as . Mathematically, the space is defined as
where and , with , denote the new basis functions obtained by the eigenfunctions . In addition, set .
In the next section, we demonstrate the efficiency of the goaloriented adaptive algorithm defined above on a problem with highcontrast multiscale coefficients. The results are compared with the standard residualbased adaptive method defined in the previous section. We note both the increased efficiency in the reduction in goalerror, , and the decreased reduction in the energynorm error with the goaloriented methods. These two observations suggest the method does what it was designed to do: focus the adaptive enrichment towards reduction in goalerror without resolving the solution where it has limited influence on the goalerror. We also note similarity in performance between the two indicators, both demonstrating errors with a similar observed rate of convergence.
5 Numerical Results
In this section, we present two numerical examples for multiscale problems with highcontrast coefficients and compare the performance of the two indicators defined in the previous section. For our simulations, we take the domain , and the inflowoutflow source term , where , and . We consider the problem of finding for the solution to (1), namely
(33) 
The goal functional
is the average value of on the outflow region . In practice, is the location of the wells, and it is important that the average pressure on is accurate. The two examples differ by the highcontrast coefficients , shown in Figure 2. Shown on the left, features a highconductivity channel crossing the domain separating the inflow and outflow; on the right, is a similar coefficient without the channel. These coefficents are visualized with the blue region indicating the value and the red region indicating the contrasts. In each example, we consider two different contrast strengths, and . We note the invariance in the relative performance of the indicators with respect to the contrast strengths.
Higher contrast (). Bottom left: , Bottom right: .
Higher contrast(), Bottom left: , Bottom right: .
For the first example, Figure 3 shows a comparison in error reduction between the three indicators for coefficient . In the two figures on the right, we compare the logarithms of the energy norm errors against the number of unknowns for three types of adaptive enrichment algorithms; namely, the residual based method described in Section 3.4 (denoted in blue in Figure 3), the residualbased goaloriented method (denoted in red in Figure 3), and the DWR type goaloriented method (denoted in black in Figure 3) as described in Section 4. From these results, we see the two types of goaloriented methods behave similarly and outperform the standard adaptive method with an improved rate of goalerror reduction. We note a more stable decrease in error reduction for the goaloriented residualtype method, but slightly improved, if less predictable error reduction for the DWRtype. This last observation is to be expected, as the DWRtype indicator does not account for the error created by using the enriched solution in place of the exact dual solution , as in (26).
On the left of Figure 3, we see the standard residual based method outperforms the goaloriented methods for the energy norm error, . This confirms that the goaloriented methods are driving the adaptivity toward a more efficient evaluation of the quantity of interest without expending additional computational effort resolving features of the solution with limited influence on the goal.
For the second example, we consider finding for that satisfies (33), with given by shown on the right of Figure 2, with no highconductivity channel. As seen in Figure 4, the results are qualitatively similar to the results of the first example. In summary, the plots on the right show the goalerror reduction for the higher and respectively lower contrast cases. The residualtype and DWRtype goaloriented indicators achieve a better rate of error reduction than the standard adaptive method, with the DWRtype showing generally the lowest error, with the leaststeady decrease. The plots on the right of Figure 4 show the reduction of energy error for the three indicators. As in the first example, the standard residual based adaptive method designed to reduce the energy error shows the best performance here, whereas the goaloriented methods yield steady error reduction but are focused on localized error reduction in the region of the goalfunctional, rather than across the entire domain.
These results demonstrate the importance of goaloriented adaptivity in GMsFEM, particularly in cases where the global domain is significantly larger than the region of infuence for the quantity of interest. In partiuclar, problems with many localized features only some of which significantly influence the quantity of interest will benefit from a goaloriented adaptive strategy.
6 Conclusion
In this paper we defined two types of error indicators that can be used in an adaptive algorithm for multiscale problems with highcontrast coefficients. The goaloriented adaptive algorithm fits within the framework of GMsFEM, and focuses the adaptivity on reducing the error in the quantity of interest, rather than in global norm. We first reviewed the general ideas of GMsFEM for highcontrast problems, then gave a detailed overview of the construction of multiscale basis functions, and a residual based adaptive algorithm designed to reduce the energy norm error. We stated the dual problem, then motivated and introduced two goaloriented error indicators, and described their use in a a goaloriented adaptive algorithm. Finally, we demonstrated the efficiency of the goaloriented algorithm and estimators compared with the standard adaptive method introduced earlier. We found for both indicators, the goaloriented method reduced the error in the goalfunction at a better rate than the standard method. We also found the two indicators perform similarly, with some increase in error reduction seen in the DWRtype indicator, but at the cost of solving the dual problem in a more enriched space, increasing the computational complexity. The residualbased indicator on the other hand may be more amenable to convergence analysis, as may be investigated in future work. The current results indicate the goaloriented strategy increases the efficiency of GMsFEM when a function of the solution rather than the solution in its entirety is of interest.
References
 [1] J.E. Aarnes, S. Krogstad, and K.A. Lie. A hierarchical multiscale method for twophase flow based upon mixed finite elements and nonuniform grids. SIAM J. Multiscale Modeling and Simulation, 5(2):337–363, 2006.
 [2] T. Arbogast, G. Pencheva, M.F. Wheeler, and I. Yotov. A multiscale mortar mixed finite element method. Multiscale Model. Simul., 6(1):319–346 (electronic), 2007.
 [3] W. Bangerth and R. Rannacher. Adaptive finite element methods for differential equations. Birkhauser, Boston, 2003.
 [4] R. Becker and R. Rannacher. An optimal control approach to a posteriori error estimation in finite element methods. Acta Numerica, 10:1–102, 2001.
 [5] P. Binev, W. Dahmen, and R. DeVore. Adaptive finite element methods with convergence rates. Numer. Math., 97(2):219–268, 2004.
 [6] M. Buerg and M. Nazarov. Goaloriented adaptive finite element methods for elliptic problems revisited. J. Comput. Appl. Math., 287:125 – 147, 2015.
 [7] J. M. Cascon, C. Kreuzer, R. H. Nochetto, and K. G. Siebert. Quasioptimal convergence rate for an adaptive finite element method. SIAM J. Numer. Anal., 46(5):2524–2550, 2008.
 [8] H. Y. Chan, E. Chung, and Y. Efendiev. Adaptive mixed gmsfem for flows in heterogeneous media. arXiv preprint arXiv:1507.01659, 2015.
 [9] Y. Chen and L. Durlofsky. An ensemble level upscaling approach for efficient estimation of finescale production statistics using coarsescale simulations. In SPE Reservoir Simulation Symposium, Houston, Texas, U.S.A., 2 2007. Society of Petroleum Engineers.
 [10] Y. Chen, L. Durlofsky, M. Gerritsen, and X. Wen. A coupled localglobal upscaling approach for simulating flow in highly heterogeneous formations. Advances in Water Resources, 26:1041–1060, 2003.
 [11] E. Chung and Y. Efendiev. Reducedcontrast approximations for highcontrast multiscale flow problems. SIAM J. Multiscale Modeling and Simulation, 8:1128–1153, 2010.
 [12] E. Chung, Y. Efendiev, and R. Gibson. An energyconserving discontinuous multiscale finite element method for the wave equation in heterogeneous media. Advances in Adaptive Data Analysis, 3:251–268, 2011.
 [13] E. Chung, Y. Efendiev, and C. S. Lee. Mixed generalized multiscale finite element methods and applications. Multiscale Modeling & Simulation, 13(1):338–366, 2015.
 [14] E. Chung, Y. Efendiev, and W. T. Leung. Generalized multiscale finite element methods for wave propagation in heterogeneous media. Multiscale Modeling & Simulation, 12(4):1691–1721, 2014.
 [15] E. Chung, Y. Efendiev, and W. T. Leung. Residualdriven online generalized multiscale finite element methods. arXiv preprint arXiv:1501.04565, 2015.
 [16] E. Chung, Y. Efendiev, and G. Li. An adaptive gmsfem for highcontrast flow problems. Journal of Computational Physics, 273:54–76, 2014.
 [17] E. Chung, Y. Efendiev, G. Li, and M. Vasilyeva. Generalized multiscale finite element methods for problems in perforated heterogeneous domains. Applicable Analysis, to appear, 2015.
 [18] L.J. Durlofsky. Numerical calculation of equivalent grid block permeability tensors for heterogeneous porous media. Water Resour. Res., 27:699–708, 1991.
 [19] Y. Efendiev, J. Galvis, and T. Hou. Generalized multiscale finite element methods. Journal of Computational Physics, 251:116–135, 2013.
 [20] Y. Efendiev, J. Galvis, R Lazarov, M Moon, and M. Sarkis. Generalized multiscale finite element method. symmetric interior penalty coupling. Journal of Computational Physics, 255:1–15, 2013.
 [21] Y. Efendiev, J. Galvis, G. Li, and M. Presho. Generalized multiscale finite element methods: Oversampling strategies. International Journal for Multiscale Computational Engineering, 12(6), 2014.
 [22] Y. Efendiev, J. Galvis, and P.S. Vassilevski. Spectral element agglomerate algebraic multigrid methods for elliptic problems with highcontrast coefficients. In Domain decomposition methods in science and engineering XIX, volume 78 of Lect. Notes Comput. Sci. Eng., pages 407–414. Springer, Heidelberg, 2011.
 [23] Y. Efendiev, J. Galvis, and X.H. Wu. Multiscale finite element methods for highcontrast problems using local spectral basis functions. Journal of Computational Physics, 230:937–955, 2011.
 [24] Y. Efendiev, T. Hou, and X.H. Wu. Convergence of a nonconforming multiscale finite element method. SIAM J. Numer. Anal., 37:888–910, 2000.
 [25] Y. Efendiev, R. Lazarov, and K. Shi. A multiscale hdg method for second order elliptic equations. part i. polynomial and homogenizationbased multiscale spaces. SIAM Journal on Numerical Analysis, 53(1):342–369, 2015.
 [26] D. Estep, M. Holst, and M. Larson. Generalized green’s functions and the effective domain of influence. SIAM J. Sci. Comput, 26:1314–1339, 2002.
 [27] M.B. Giles and E. Süli. Adjoint methods for PDEs: a posteriori error analysis and postprocessing by duality. Acta Numerica, 11:145–236, 2003.
 [28] T. Grätsch and K.J. Bathe. A posteriori error estimation techniques in practical finite element analysis. Comput. Struct., 83(45):235 – 265, 2005.
 [29] H. Hajibeygi, D. Kavounis, and P. Jenny. A hierarchical fracture model for the iterative multiscale finite volume method. Journal of Computational Physics, 230(4):8729–8743, 2011.
 [30] M. Holst. Adaptive numerical treatment of elliptic systems on manifolds. Adv. Comput. Math., 15(1–4):139–191, 2001. Available as http://arxiv.org/abs/1001.1367 arXiv:1001.1367 [math.NA].
 [31] M. Holst, S. Pollock, and Y. Zhu. Convergence of goaloriented adaptive finite element methods for semilinear problems. Comp. Vis. Sci., 17(1):43–63, 2015.
 [32] T. Hou and X.H. Wu. A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys., 134:169–189, 1997.
 [33] P. Jenny, S.H. Lee, and H. Tchelepi. Multiscale finite volume method for elliptic problems in subsurface flow simulation. J. Comput. Phys., 187:47–67, 2003.
 [34] C. Jhurani and L. Demkowicz. Multiscale modeling using goaloriented adaptivity and numerical homogenization. part i: Mathematical formulation and numerical results. Computer Methods in Applied Mechanics and Engineering, 213âÃÃ¬216:399 – 417, 2012.
 [35] K. Mekchay and R. Nochetto. Convergence of adaptive finite element methods for general second order linear elliptic PDE. SIAM J. Numer. Anal., 43(5):1803–1827, 2005.
 [36] M. S. Mommer and R. Stevenson. A goaloriented adaptive finite element method with convergence rates. SIAM J. Numer. Anal., 47(2):861–886, 2009.
 [37] R. H. Nochetto, K. G. Siebert, and A. Veeser. Theory of adaptive finite element methods: an introduction, pages 409 – 542. Springer, 2009.
 [38] J. T. Oden and K. Vemaganti. Estimation of local modeling error and goaloriented adaptive modeling of heterogeneous materials: I. error estimates and adaptive algorithms. J. Comput. Phys., 164:22–47, 2000.
 [39] S. Prudhomme and J. T. Oden. On goaloriented error estimation for elliptic problems: application to the control of pointwise errors. Comput. Method Appl. M., 176(14):313–331, 1999.
 [40] R. Stevenson. Optimality of a standard adaptive finite element method. Found. Comput. Math., 7(2):245–269, April 2007.
 [41] R. Verfürth. A review of a posteriori error estimation and adaptive meshrefinement techniques. Advances in numerical mathematics. WileyTeubner, 1996.