DG-GMsFEM for Problems in Perforated Domains with Non-Homogeneous Boundary Conditions

: Problems in perforated media are complex and require high resolution grid construction to capture complex irregular perforation boundaries leading to the large discrete system of equations. In this paper, we develop a multiscale model reduction technique based on the Discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM) for problems in perforated domains with non-homogeneous boundary conditions on perforations. This method implies division of the perforated domain into several non-overlapping subdomains constructing local multiscale basis functions for each. We use two types of multiscale basis functions, which are constructed by imposing suitable non-homogeneous boundary conditions on subdomain boundary and perforation boundary. The construction of these basis functions contains two steps: (1) snapshot space construction and (2) solution of local spectral problems for dimension reduction in the snapshot space. The presented method is used to solve different model problems: elliptic, parabolic, elastic, and thermoelastic equations with non-homogeneous boundary conditions on perforations. The concepts for coarse grid construction and deﬁnition of the local domains are presented and investigated numerically. Numerical results for two test cases with homogeneous and non-homogeneous boundary conditions are included, as well. For the case with homogeneous boundary conditions on perforations, results are shown using only local basis functions with non-homogeneous boundary condition on subdomain boundary and homogeneous boundary condition on perforation boundary. Both types of basis functions are needed in order to obtain accurate solutions, and they are shown for problems with non-homogeneous boundary conditions on perforations. The numerical results show that the proposed method provides good results with a signiﬁcant reduction of the system size. solution as a reference solution. To reduce the size of the ﬁne grid system, we present a Discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM) with the construction of two types of multiscale basis functions for subdomain boundary and perforation boundary, separately. Presented perforation boundary basis functions are used to approximate non-homogeneous boundary conditions on perforations. Numerical investigation of the presented method was performed for four model problems: (1) elliptic, (2) parabolic, (3) elasticity, and (4) thermoelasticity problems in the perforated domain. The results for two cases with homogeneous and non-homogeneous boundary conditions are given. For the case with homogeneous boundary conditions on perforations, the given results are obtained using only subdomain boundary basis functions. However, for a non-homogeneous boundary condition, both subdomain and perforation boundary basis functions should be used. Numerical results are presented for different concepts of the coarse grid construction (structured, quasi-structured, and unstructured coarse grids). Numerical results show that the proposed method can provide good results and give a signiﬁcant reduction of the system size with appropriate choosing construction of the multiscale basis functions for problems with homogeneous and non-homogeneous perforation boundary conditions. In future works, we plan to consider oversampling techniques for multiscale basis function construction, construct a multiscale solver for three-dimensional problems, and consider construction of the accurate multiscale techniques for nonlinear problems in heterogeneous perforated media.


Introduction
Problems in perforated domains are of great interest additionally proposing many real-world applications. Take, for example, diffusion in perforated domains, mechanical processes in granular media, pore-scale flows in porous media, and so on [1][2][3][4]. Problems in perforated domains with non-homogeneous boundary conditions on perforations have great importance for a lot of applications in physics, biology, geology, and chemistry [5][6][7][8].
The main characteristic of perforated domains is the multiscale nature of the underlying processes. The solution techniques for these problems require high resolution in grid construction to capture complex irregular boundaries of perforations. Using direct numerical methods to solve these problems is computationally challenging and expensive. In order to reduce the size of the system with accurate approximation, we use homogenization techniques and multiscale methods.
Reduction techniques are widely used today for problems in heterogeneous perforated media to reduce computational cost. In recent years, many methods have been developed focused on obtaining solutions on a coarse mesh. The homogenization method can be used for problems in perforated domains with scale separation [8][9][10]. For example, in Reference [11], the author considers the periodic homogenization problem in perforated domains that are formed by removing a periodic array of small holes from a fixed open bounded and connected domain with regular boundary. In the numerical homogenization method, the approximation of the solution on a coarse grid is constructed by calculating effective properties for coarse grid cells [12]. In Reference [13], the authors construct the machine learning methods for a fast calculation of effective characteristic for domains with random inclusions. In addition, for the perforated media without scale separation, the multiscale techniques are widely used [14,15]. Multiscale methods construct multiscale basis functions in each coarse-grid cell and couple these basis functions in a global formulation. Several studies have shown employment of multiscale methods as solution to problems in perforated domains [16][17][18][19]. There are also studies presenting the multiscale approach based on Crouzeix-Raviart coupling of multiscale finite element basis [16,17]. The generalization of the heterogeneous multiscale finite element method for elliptic problems in perforated domains is presented in work of Henning and Ohlberger (2009) [18]. In order to avoid a limited number of degrees of freedom per coarse element, this paper considers the Generalized Multiscale Finite Element Method (GMsFEM) for solution problems in perforated domains [19][20][21][22]. The GMsFEM is a general multiscale procedure, wherein the model reduction is based on some local multiscale basis functions that are constructed using local spectral decomposition. In our earlier study [19], we presented an accurate multiscale approximation using GMsFEM with several basis functions in each local domain. The problems with non-homogeneous boundary conditions on perforations are studied in References [23,24], where authors evaluate construction of the additional basis functions for perforation boundary. The multiscale model reduction technique, based on Discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM), is presented in Reference [22], where we consider the solution of the problems with homogeneous boundary conditions. Additional points of multiscale methods can distinguish methods, such as upscaling method using non-local multicontinuum method (NLMC) [24,25], which can be effective for heterogeneous. The main part of the method is the construction of suitable local basis functions with the capability of capturing multiscale features and non-local effects. Another method that deserves attention is the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) [26,27]. The main purpose of this method is that the convergence of the method is independent of contrast and decreases linearly with respect to the cell size if the oversampling size is chosen appropriately. In addition, for punctured tasks, the oversampling method can be especially useful and improve accuracy and convergence [28][29][30]. Besides, the oversampling strategy is used to reduce the mismatching effects of boundary conditions imposed artificially in the construction of snapshot basis functions.
This paper presents an extension of the DG-GMsFEM for problems with non-homogeneous boundary conditions on perforations. The presented approach is based on the construction of a separate multiscale space for coarse cell interfaces and for local perforation boundary. We present a unified approach for different types of problems: elliptic, parabolic, elastic, and thermoelastic problems. Using the DG-GMsFEM for coarse grid coupling, we construct multiscale basis functions to generate an accurate lower-dimensional model on a coarse grid. The paper starts with the construction of the snapshot space that contains a set of local solutions in each coarse cell with various boundary conditions. In order to perform a dimension reduction of the snapshot space and define multiscale basis functions, we solve a local spectral problem in the snapshot space. The construction of two types of multiscale basis functions for subdomain boundary and perforation boundary are presented separately. The presented perforation boundary basis functions are used to approximate non-homogeneous boundary conditions on perforations. The numerical results for two-dimensional problems with homogeneous and non-homogeneous boundary conditions on perforations are presented. We investigate a different number of subdomain and perforation boundary multiscale basis functions and show that an accurate solution can be obtained with a large reduction of the system size with a sufficient number of multiscale basis functions. The coarse grid construction and definition of the local domains are presented and investigated numerically for each of the considered equations.
This paper is organized as follows. Section 2 contains description of the problem formulation. The fine-scale approximation using the interior penalty discontinuous Galerkin (IPDG) method is given in Section 3. The multiscale method based on DG-GMsFEM is presented in Section 4 with the description of the construction of the multiscale basis functions. Section 5 presents numerical results for the two-dimensional problem in a perforated domain with homogeneous and non-homogeneous boundary conditions on perforations. Results are shown for the different number of multiscale basis functions, as well as results for different types of problems (elliptic, parabolic, elastic, and thermoelastic). The results for quasi-structured and unstructured coarse grid are presented. Finally, the conclusion is given in Section 6.

Problem Formulation
Let Ω be a perforated domain and Γ p be the perforation boundary. In this work, we consider the elliptic equation in perforated domain and non-homogeneous boundary conditions on perforation boundary where f is a given source vector, L denotes a linear differential operator, B is a normal derivative operator, and ∂Ω = Γ g ∪ Γ p .
• For the Laplace operator, we have with where q(u) is the flux, k is the diffusion coefficient, and n is the outward unit normal on ∂Ω. • For the elasticity operator, we have with where u is the displacement field, ε(u) is the strain tensor, σ(u) is the stress tensor, u Tr is the transpose of u, and λ and µ are the Lamé coefficients.

Fine-Grid Approximation
This section shows an approximation using the interior penalty discontinuous Galerkin method (IPDG) on the fine grid that resolves all perforations on the grid level.
Let T h be a fine-grid partition of the domain Ω given by where N h cell is the number of fine grid cells. We use E h to denote the set of facets in T h with are the set of interior and boundary facets (see Figure 1).
where u + = u| K + , u − = u| K − , and K + and K − are the two cells sharing the facet E. For boundary facets, we have [u] = u| E and {u} = u| E , E ∈ E h b . For IPDG approximation, the following variational formulation of the problem is used: where: • for the Laplace operator, we have • for the elasticity operator, we have Here, γ f is the penalty parameter. The above systems can be written in the matrix form as follows:

Multiscale Method
This section describes the construction of multiscale approximation on the coarse grid using the Discontinuous Galerkin Generalized Multiscale Finite Element method (DG-GMsFEM) [22,31,32].
Let T H be a coarse-grid partition of the domain Ω with coarse mesh size H (see Figure 2).
where N H cell is the number of coarse grid cells, and K i is coarse cell (local domain). Use E H to denote the set of facets in T H with E H = E H o ∪ E H b . Note that, in DG-GMsFEM, the multiscale basis functions are supported in each coarse cell K i .
Defining V H as the multiscale space parabolic, elastic, and thermoelastic N u = dim(V H ) is the number of basis functions. For the coarse grid approximation, we use a DG approach and have the following variational formulation: find u H ∈ V H such that Note that the coarse-scale system can be formed by projecting the fine-scale system onto the coarse grid. The projection matrix can be assembled using the multiscale basis functions.
To construct the multiscale space, we start with the construction of the snapshot space that contains a set of basis functions formed by the solution of local problems with all possible boundary conditions up to the fine grid resolution in each coarse cell K i (local domain) for i = 1, · · · , N, parabolic, elastic, and thermoelastic N is the number of coarse blocks in Ω. After that, we solve a spectral problem to select dominant modes of the snapshot space.
Next, we present the details of constructing the multiscale basis functions considering two types of basis functions related to the two boundaries: (1) subdomain boundaries Γ i g and (2) perforation boundary Γ i p . Subdomain boundary muliscale basis functions. In the local snapshot space consisting of functions u i l , which are solutions to the following local problem with the following boundary condition on subdomain boundaries, and, on perforation boundary, we set homogeneous boundary condition related to (3) where l = 1, · · · , L g i . Local problems are solved using IPDG approximation on fine mesh T h (K i ). For Laplace operator, we have L g i = J g i local problems, where J g i is the number of fine grid facets on Γ i g , and g l i = δ l i is the Kronecker delta function that has value 1 if i = l and value 0 else. For elasticity operator, we have L g i = d · J g i local problems, where d is the dimension, and g l i = (δ l i , 0) and (0, δ l i ) for d = 2. The collection of the solutions of the above local problems generates the snapshot space in the local domain K i To reduce the size of the snapshot space, we solve the following local spectral problem in the snapshot space V i,snap gÃ K i gψ i whereÃ Here, A K i h and S K i h are the matrix representation of the bilinear forms a K i (u, v) and For Laplace operator: • For elasticity operator: It should be noted that the integral in s K i (u, v) is defined on the boundary of the coarse block K i and based on the definition of the snapshot space used to extract dominant modes related to the outer boundary of the local domain.
Next, we arrange the eigenvalues in increasing order and choose the first eigenvectors corresponding to the first smallest eigenvalues ψ i g,k = (R i,snap g ) Trψi g,k as the basis functions (k = 1, . . . , M i g ) V g H = span{ψ i g,k : 1 ≤ i ≤ N H cell , 1 ≤ k ≤ M i g }. The first five eigenvectors for some local domains K i are depicted in Figure 3 for Laplace operator and in Figure 4 for elasticity operator.

Perforation boundary muliscale basis functions.
To handle non-homogeneous boundary conditions on the perforation boundaries, we construct an additional multiscale basis functions. The snapshot space is constructed by solution of the following problem in local domain K i that contains perforations with the homogeneous boundary condition on subdomain boundaries, and, on the perforation boundary, we set the following boundary condition related to (3) We perform a dimension reduction in the snapshot space using the local spectralproblem The eigenvalues are arranged in increasing order, and, by choosing the first eigenvectors corresponding to the first smallest eigenvalues, we define perforation boundary multiscale basis functions ) Trψi g,k for k = 1, . . . , M i p , and N cell,p is the number of local domains with perforations. It should be noted that the first eigenvector was not taken because constant valued vector already exists in outer boundary multiscale space. The perforation boundary of multiscale basis functions are presented in Figure 5 for Laplace problem and in Figure 6 for elasticity problem.    Interior multiscale basis functions. We add one interior basis functions [22,32] constructed by solution of the following local problem: with f = 1 and homogeneous boundary condition on subdomain and perforation boundaries, Finally, there is the following multiscale space: where N p = N cell,p , and N g = N cell Coarse scale system. To construct the coarse grid system, we generate a projection matrix using multiscale basis functions Using projection matrix, we have the following coarse grid system in matrix form: where After the solution of the coarse-scale system, the fine-scale solution for displacement can be recovered U ms = R Tr U H .

Numerical Results
We present numerical results for the model problems in the perforated domain Ω = [0, L x ] × [0, L y ] with L x = L y = 1. The computational domain with coarse and fine grids is presented in Figure 7. The fine grid contains 14,648 vertices and 28,410 cells. The coarse grid contains 121 vertices and 100 cells (local domains). To construct the computational domain with computational mesh, we use the mesh generator Gmsh [33]. To investigate the presented multiscale method for problems solving in perforated domains, the following tests are considered: (1) elliptic equation and (2) elasticity equation. The extension of the method for solution of the parabolic and thermoelasticity problems are presented in Appendixes A and B . The numerical implementation of the model problems is based on the FEniCS computing platform [34]. There are two test cases for each problem: • Case 1. homogeneous boundary conditions on perforations; • Case 2. non-homogeneous boundary conditions on perforations.
We calculate the following relative errors in L 2 and energy norm (semi-norm H 1 ) between multiscale and fine-grid solution: where u and u ms are the fine-scale and multiscale solutions, respectively. To start, there are results for elliptic and elastic equations. Whereat there is evaluation of the unstructured coarse grids and heterogeneous perforated domains. The results for parabolic and thermoelasticity problems are given in Appendixes A and B .

Results for Elliptic Equation
We consider the elliptic equation in the perforated domain Ω with k = 1 and the following boundary conditions: • Case 1. Homogeneous boundary conditions on perforations: • Case 2. Non-homogeneous boundary conditions on perforations: We consider three different α = 1, 25, and 100. Tables 1 and 2 cell is the total number of coarse grid cells, and N p = N H cell,p is the number of coarse grid cells with perforations. For the problem with homogeneous boundary conditions (Case 1), we use multiscale space without perforation basis functions (M p = 0). Results for multiscale spaces with and without perforation boundary basis functions (M p = M) are presented for Case 2 with non-homogeneous boundary conditions. For Case 1, we observe that the relative error in L 2 norm is reduced to 1% for α = 1, 25, and 100, when we take the sufficient number of multiscale basis functions (M g = 12 and M p = 0). For Case 2 of the boundary conditions, the multiscale method gives a solution with a large error when perforation boundary multiscale basis functions is not taken (M p = 0). However, we obtain good results after adding perforation boundary bases (M p = M). For example, the relative L 2 error for M = 12 multiscale basis functions is 12% for M p = 0 and reduced to 1.2% for M p = M (α = 25). For the test problem with larger α = 100, L 2 error from 21% to 0.67% is reduced. The graphical representation of the Tables 1 and 2 is given in Figure 8.       For Case 1 with homogeneous boundary conditions on perforations and Case 2 with non-homogeneous boundary conditions on perforations, based on results, it could be said that the presented multiscale method shows good results with small errors and demonstrates high accuracy for the both cases. Using non-homogeneous boundary conditions on perforations, the multiscale method using perforation boundary basis functions greatly improves results. We consider the elasticity equation with λ = µ = 1 with the following boundary conditions:

Results for Elasticity Equation
and • Case 1. Homogeneous boundary conditions on perforations: • Case 2. Non-homogeneous boundary conditions on perforations:   Tables 3 and 4 and Figure 13 present relative errors for elasticity problem for Case 1 and Case 2, respectively. Evidently the L 2 error reduce from 117% for M p = 0 to 2.3% for M p = M in Case 2 with M = 24. For Case 1, we have 2% of L 2 error for M g = M = 24 with M p = 0. For both cases, we obtain good results with a sufficient number of multiscale basis functions. In Case 1 with homogeneous boundary conditions for elasticity on perforations, good results with outer basis functions are obtained. In Case 2 with non-homogeneous boundary conditions for elasticity on perforations, it is seen that, for displacement, we should use perforation boundary basis functions in order to approximate the non-homogeneous boundary condition.

Unstructured Coarse Grids
We consider a solution of the elliptic and elasticity problems using the presented multiscale method on quasi-structured and unstructured coarse grids. Figure 14 presents quasi-structured and unstructured coarse grids with an illustration of the local domains. Coarse grids contain 100 local domains. Problems with Case 2 boundary conditions are considered for numerical investigation. Table 5 presents results for elliptic equation with α = 100 for non-homogeneous boundary conditions (Case 2). The results for elasticity equation are presented in Table 6. Numerical results are shown for multiscale space with outer and boundary basis functions (M g = M p = M) on quasi-structured and unstructured coarse grids with 100 local domains. Using 32 multiscale basis functions (M p = M g = 32), similar results are obtained with less than 1% of L 2 error for structured, quasi-structured, and unstructured coarse grids with 100 local domains with DOF H = 3844 (4.5% of DOF f = 85,230). By comparing results with less number of multiscale basis functions (M p = M g = 8), it is observed that the error for structured coarse grid is smaller (1% and 19% for L 2 and H 1 errors) than for unstructured coarse grid (4.9% and 33% for L 2 and H 1 errors). Furthermore, the results show that the quasi-structured grid give a smaller errors than unstructured coarse grid (2% and 23% for L 2 and H 1 errors). For the elasticity problem with M = 32, we have 2% and 13% for L 2 and H 1 errors for structured grid, 3% and 15% for L 2 and H 1 errors for quasi-structured grid, and 5% and 22% for L 2 and H 1 errors for unstructured grid.  For elasticity problem with a smaller number of multiscale basis functions, there less errors for a structured grid than an unstructured grid. By comparison of the quasistructured and unstructured grids, we see less errors for the quasi-structured coarse grid. Good results with less errors for any coarse grids can be obtained with a large number of multiscale basis functions. It should be noted that the main advantage of using unstructured coarse grids is the similar number of cells in each local domain (load balancing). Conforming triangulation of the domain with coarse edges is used in construction of the structured grid, which can be difficult for complex geometries with a large number of perforations. Though the conforming construction is not needed in the quasi-structured grid. Numerical results show that the presented multiscale method works very well with any coarse grid construction concepts.

Heterogeneous Coefficients
Finally, we consider the efficiency of the presented method for the solution of the elliptic and elasticity problems with heterogeneous coefficients (k = k(x) and E = E(x)). The heterogeneous elasticity parameter and heterogeneous diffusion coefficient are presented in Figure 15. The parameters λ and µ are given as follows: is the elasticity parameter, and ν = 0.3 is the constant Poisson's ratio. Figure 15. Heterogeneous coefficient. Table 7 presents the relative L 2 and H 1 errors for elliptic (left) and elasticity (right) problems with heterogeneous coefficients. Case 2 with non-homogeneous boundary conditions for perforation is considered (M g = M p = M). For elliptic problem, the results are presented with α = 100. Numerical results for elliptic and elasticity problems show that the presented multiscale method works well with heterogeneous coefficients.

Conclusions
This paper presents the multiscale method for solution problems in the perforated domain with non-homogeneous boundary conditions on perforations. For approximation on the fine grid that resolved perforations on the grid level, we apply the Discontinuous Galerkin finite element method and use the solution as a reference solution. To reduce the size of the fine grid system, we present a Discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM) with the construction of two types of multiscale basis functions for subdomain boundary and perforation boundary, separately. Presented perforation boundary basis functions are used to approximate non-homogeneous boundary conditions on perforations. Numerical investigation of the presented method was performed for four model problems: (1) elliptic, (2) parabolic, (3) elasticity, and (4) thermoelasticity problems in the perforated domain. The results for two cases with homogeneous and non-homogeneous boundary conditions are given. For the case with homogeneous boundary conditions on perforations, the given results are obtained using only subdomain boundary basis functions. However, for a non-homogeneous boundary condition, both subdomain and perforation boundary basis functions should be used. Numerical results are presented for different concepts of the coarse grid construction (structured, quasistructured, and unstructured coarse grids). Numerical results show that the proposed method can provide good results and give a significant reduction of the system size with appropriate choosing construction of the multiscale basis functions for problems with homogeneous and non-homogeneous perforation boundary conditions. In future works, we plan to consider oversampling techniques for multiscale basis function construction, construct a multiscale solver for three-dimensional problems, and consider construction of the accurate multiscale techniques for nonlinear problems in heterogeneous perforated media.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Parabolic Problem
We consider the parabolic problem in perforated domain Ω: with the following boundary conditions and initial condition T = T 0 in Ω for t = 0. For approximation by time, we use an implicit approximation with a time step τ. We have the following fine grid approximation on T h using the IPDG method: find u h ∈ V h such that 1 andǔ h is the solution from the previous time step. The matrix form is following In the multiscale method, we construct a multiscale space similar to the Laplace problem and generate a projection matrix R = ϕ 1 g,1 , . . . , ϕ Finally, we obtain the following coarse grid system for a parabolic problem and calculation of the coarse-scale solution, we reconstruct the fine-scale solution U ms = R Tr U H . Numerical results. We consider the parabolic problem with k = 1. We set the following boundary conditions: • Case 1. Homogeneous boundary conditions on perforations: and perform simulation for t max = 0.02. • Case 2. Non-homogeneous boundary conditions on perforations: and perform simulation for t max = 0.01. We consider three different α = 1, 25, and 100. We set zero initial condition and perform a simulation with 40 time steps.    The reference and multiscale solutions at the different times (t m with m = 5, 15, and 40) are presented in Figure A1 for Case 1 and in Figure A2 for Case 2. The multiscale solution is obtained using 12 multiscale basis functions (M = 12). In the first row, there is a fine-scale solution with DOF h = 85,230. In the second row, there is a multiscale solution with DOF H = 1300 for M p = 0. In the third row, we give a multiscale solution with DOF H = 1776 for M p = M in Case 2. We observe a bad multiscale approximation near the perforation boundary for Case 2 with e L 2 = 22% at final time, if a perforation boundary basis functions is not used to handle non-homogeneous boundary conditions. For M p = M in Case 2, we observe a good multiscale solution with e L62 = 0.7% at final time. The relative errors for parabolic problem (Case 1) between reference solution and multiscale solution for α = 1, 25, and 100 are shown in Table A1 and in Figure A3 (left picture). We have a good multiscale solution with nearly 1-2% of L 2 error when we take 12 outer boundary multiscale basis functions. Table A2 and Figure A3 (right picture) present relative errors for Case 2. Similar to the previous results for the elliptic equation, we see large errors for M p = 0 and can obtain a good multiscale solution when we take the sufficient number of multiscale basis functions with outer and perforation boundary multiscale basis functions. For example, if 12 multiscale basis functions (M = 12) are taken, we have 7-22% of L 2 errors for M p = 0 and reduce errors to 1% when we add perforation boundary basis functions for α = 1, 25, and 100. Table A3 shows results for non-homogeneous boundary conditions (Case 2) on quasistructured and unstructured coarse grids with 100 local domains. We present results for α = 100. Using 32 multiscale basis functions (M p = M g = 32), we obtain similar results with less than 1% of L 2 error for structured, quasi-structured, and unstructured coarse grids with 100 local domains with DOF H = 3844 (4.5% of DOF f = 85,230).

Appendix B. Thermoelasticity Problem
We consider the thermoelasticity problem in Ω that is described by a system of equations for temperature and displacement: We consider Equation (A2) with the following boundary conditions: and initial condition T = T 0 , u = 0 in Ω for t = 0. For the IPDG method on the fine grid T h , we have the following variational formulation: The matrix form is the following: . For the construction of the multiscale method for the thermoelasticity problem, we construct multiscale spaces for displacement and temperature, separately. Then, we generate projection matrices using multiscale basis functions Finally, we construct the following coarse grid system for the thermoelasticity problem   Figure A4 show relative errors for the thermoelasticity problem between a fine-grid and multiscale solution with different multiscale basis functions at the final time. In presented results, we tested different approaches for choosing M u g , M u p and M T g , M T p for the coupled process. Here, M v are the number of multiscale basis functions for displacement (v = u) and temperature (v = T), DOF v h and DOF v H are the numbers of degrees of freedom for reference (fine-grid) and multiscale solution, and e(v) L 2 and e(v) H 1 are the L 2 and H 1 relative errors for v = u, T. Here, DOF h = DOF u h + DOF T h with DOF u h = 6 · N h cell and DOF T h = 3 · N h cell for fine grid solution. For multiscale solution, we have DOF H = DOF u H + DOF T H with DOF u H = (2 · M u g + 1) · N H cell + 2 · M u p · N H cell,p and DOF T H = (M T g + 1) · N H cell + M T p · N H cell,p . The results for quasi-structured and unstructured coarse grids with 100 local domains are presented in Table A5. Numerical results are shown for multiscale space with outer and boundary basis functions (M g = M p = M) for Case 2. We obtain similar results with 4% and 0.3% for displacement and temperature L 2 errors for all grids for M = 32. Figures A5 and A6 show the results of the fine-grid and multiscale solvers for the thermoelasticity problem at the final time on structured coarse grid. The temperature distribution is presented on the left figure, the displacement X and Y components are shown in the center and right figures. In presented results, we used 24 multiscale basis functions for displacement and temperature (M u = M T = 24, M u g = M u p = M u , and M T g = M T p = M T ). Reference solution with DOF u h = 170,460 for displacement, DOF T h = 85,230 for temperature is shown in the first row. Solution on a coarse grid using outer and perforation boundary multiscale basis functions with DOF u H = 5988 for displacement and DOF T H = 3044 for temperature is shown in the second row.   Case 1 with non-homogeneous boundary conditions for elasticity and homogeneous boundary conditions for temperature on perforations shows that, for displacement, we should use perforation boundary basis functions in order to approximate the nonhomogeneous boundary condition for the elasticity part of the coupled system. For temperature, only outer boundary basis functions can be used because we have homogeneous boundary conditions for temperature. For Case 1 with M u g = M u p = 24, M T g = 24, and M T p = 0, we obtain a good results with 3% and 8% of L 2 errors for displacement and temperature, respectively. When we use only outer boundary basis functions for displacement (M u p = 0), we obtain a large error for both displacement and temperature. In Case 2 with homogeneous boundary conditions for elasticity and non-homogeneous boundary conditions for temperature on perforations, we can obtain good results with outer and perforation boundary basis functions for both displacement and temperature (M u g = M u p = M u and M T g = M T p = M T ). It should be noted that perforation boundary basis functions should be added for displacement due to the coupled process for accurate multiscale approximation even for non-homogeneous boundary conditions for elasticity on perforations.