Next Article in Journal
Particulate Matter and COVID-19 Disease Diffusion in Emilia-Romagna (Italy). Already a Cold Case?
Previous Article in Journal
Studying a Tumor Growth Partial Differential Equation via the Black–Scholes Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mixed Generalized Multiscale Finite Element Method for a Simplified Magnetohydrodynamics Problem in Perforated Domains

1
Multiscale Model Reduction Laboratory, North-Eastern Federal University, 677007 Yakutsk, Russia
2
Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Key Laboratory of Intelligent Computing & Information Processing of Ministry of Education, School of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
3
Institute for Scientific Computation, Texas A&M University, College Station, TX 77843, USA
4
Department of Mathematics, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong 999077, China
5
Department of Mathematics, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Submission received: 22 May 2020 / Revised: 11 June 2020 / Accepted: 17 June 2020 / Published: 23 June 2020

Abstract

:
In this paper, we consider a coupled system of equations that describes simplified magnetohydrodynamics (MHD) problem in perforated domains. We construct a fine grid that resolves the perforations on the grid level in order to use a traditional approximation. For the solution on the fine grid, we construct approximation using the mixed finite element method. To reduce the size of the fine grid system, we will develop a Mixed Generalized Multiscale Finite Element Method (Mixed GMsFEM). The method differs from existing approaches and requires some modifications to represent the flow and magnetic fields. Numerical results are presented for a two-dimensional model problem in perforated domains. This model problem is a special case for the general 3D problem. We study the influence of the number of multiscale basis functions on the accuracy of the method and show that the proposed method provides a good accuracy with few basis functions.

1. Introduction

The magnetohydrodynamic (MHD) system describes the interactions of electrically conducting incompressible flows in the presence of a magnetic field. MHD flows are usually involved in the contexts: liquid metal cooling of nuclear reactors, electromagnetic casting of metals, MHD generators, accelerators, and MHD ion propulsion, see [1,2,3]. The equations that model the MHD are the coupling of incompressible Navier–Stokes equations and the Maxwell’s equations via the Lorentz force and Ohm’s law. There are some research works devoted to the numerical approximation of MHD equations. For treating the nonlinear terms effectively, three classical iterative methods (Stokes-type, Newtown and Oseen-type ones) within finite element approximation for the steady MHD equations have been developed in [4,5]. Especially, the theoretical analysis and numerical examples in [4,6] show that the Picard iteration (called Oseen-type iteration by the authors in the reference) is suitable for high Reynolds number problems. The local and parallel finite element algorithms based on two-grid methods were given in [7,8,9]. The stabilized finite element schemes that capture physical solutions were investigated in [10]. The mixed finite element methods with divergence-free velocities and magnetics fields are developed [11,12,13]. However, analysis of the above literature mainly focuses on the singly-connected computational domains. In this paper, our goal is to develop multiscale methods for MHD flows on complex domains with perforations.
In this paper, we consider a simplified MHD problem in perforated domains, and propose a multiscale method for constructing the reduced order model on the coarse grid. In particular, we consider a special 2D case of MHD equations and the full MHD problem will be considered in the future. We start with the mixed formulation for both equations of the simplified MHD problem in a perforated domain. For the solution of the nonlinear system, we use a Picard iteration to decouple the equations. We construct a fine grid that resolves perforations on the grid level and write an approximation using a mixed finite element method. Problems in perforated domains have a multiscale nature and approximations on the fine grid lead the large system of equations. To reduce the size of the discrete system, multiscale methods or homogenization techniques are needed [14,15,16,17,18,19]. Homogenization techniques are constructed to solve problems with scale separation, i.e., when perforations are of similar sizes. When the perforated media do not have scale separation, multiscale methods are used [20,21,22,23,24]. In the papers [20,21], the authors presented the multiscale finite element method for problems in perforated domains, but these approaches use a limited number of degrees of freedom per coarse element. In the papers [18,25], the general approach for constructing the multiscale basis functions using local spectral characteristics was presented in [18,25]. Continuous and Discontinuous Galerkin Generalized Multiscale Finite Element Methods (GMsFEM) for coarse grid approximation of the problems in perforated domains was presented in [26,27]. The GMsFEM is a systematic approach to identify multiscale basis functions via local snapshots and local spectral problems. The local snapshots are constructed by solving local problems and contain the information about local geometry structure (perforations). By performing local spectral decomposition, the method identifies multiscale basis functions. For accurately approximating the fluxes, we use a Mixed Generalized Multiscale Finite Element method and construct multiscale basis functions for them [28,29,30,31]. In Mixed GMsFEM, we construct the snapshot space for each coarse edge, perform the local spectral decomposition in the snapshot space and select the dominant modes as multiscale basis functions for the flux. The presented work is based on our previous papers [30,31], where we presented a Mixed GMsFEM for solutions of the problems in perforated domains. In this paper, we apply and study the method for solving the incompressible magnetohydrodynamics (MHD) problem in perforated domains. We present numerical results for some perforated geometries 2D. These numerical results demonstrate one can achieve a good accuracy with fewer basis functions.
The paper is organized as follows. In Section 2, we describe the problem formulation. Next, we decouple the system of equations using Picard iterations. In Section 3 and Section 4, we present a fine grid and multiscale approximations for the magnetic field and flow problems, respectively. Fine grid and coarse grid approximations are presented using mixed formulation. Construction of the multiscale basis functions are presented. In Section 5, we present numerical results for the two-dimensional problem in a perforated domain. The paper ends with a conclusion.

2. Problem Formulation

The full steady incompressible MHD equations in Ω R 3 is a complex system which link the velocity field u = ( u 1 , u 2 , u 3 ) , magnetic field B = ( B 1 , B 2 , B 3 ) and pressure p:
R e 1 Δ u + u · u + p S c curl B × B = f , x Ω , R e m 1 curl ( curl B ) curl ( u × B ) = g , x Ω , · u = 0 , x Ω , · B = 0 , x Ω ,
where R e is the hydrodynamic Reynolds number, R e m the magnetic Reynolds number, S c the coupling number. We refer the readers to [32] for details on the existence, uniqueness, and finite element approximation of solutions of (1).
This article considers a magnetic field B = ( 0 , 0 , B ) perpendicular to the velocity field u = ( u 1 , u 2 , 0 ) of (1) in two-dimensional perforated domain Ω , then the simplified MHD equations without any body force and hydrodynamic convective terms are described as follows:
· ( u B ) · ( R e m 1 B ) = 0 , x Ω , R e 1 Δ u + p S c B B = 0 , x Ω , · u = 0 , x Ω .
Let D = R e m 1 , q = D B be the flux and Ω is the perforated domain (see Figure 1). We write Equation (2) in the mixed form
D 1 q + B = 0 , x Ω , · q + · ( u B ) = 0 , x Ω , R e 1 Δ u + p + S c D 1 B q = 0 , x Ω , · u = 0 , x Ω ,
and consider with the following boundary conditions
u = 0 , x Γ P , ( u p I ) · n = 0 , x Ω / Γ P ,
B = 0 , x Γ P , B = g , x Γ 1 , q · n = 0 , x Γ 2 ,
where Γ P is the boundary of perforations, Γ 1 Γ 2 Γ P = Ω , n is the unit outward normal vector on Ω and I is the n × n identity matrix.
For linearization, we use Picard iterations and obtain the following algorithm
  • Find ( q k + 1 , B n + 1 ) from
    D 1 q k + 1 + B k + 1 = 0 , q k + 1 + ( u k B k + 1 ) = 0 ,
  • Find ( u k + 1 , p n + 1 ) from
    R e 1 Δ u k + 1 + S c p k + 1 + S c D 1 B k + 1 q k + 1 = 0 , u k + 1 = 0 .
Here k is the nonlinear iteration.

3. Magnetic Field Problem

We have the following variational formulations: find ( q k + 1 , B n + 1 ) V q × Q B such that
m ( q k + 1 , v ) + g ( B k + 1 , v ) = f q ( v ) v V q , g ( q k + 1 , r ) + c ( B k + 1 , r ) = 0 , r Q B ,
where k is the nonlinear iteration and
m ( q , v ) = Ω D 1 q · v d x , g ( B , q ) = Ω B q d x ,
c ( B , r ) = Ω ( u k B ) · r d x , f q ( v ) = Γ 1 g v · n d s .
We can rewrite fine grid approximation in the matrix form
M G T G C k q k + 1 B k + 1 = F q 0
where C k = C ( u k ) . The fine grid problem is solved using the lowest-order Raviart-Thomas element.

Multiscale Approximation

To construct the multiscale space for the coarse scale approximation, we use a Mixed GMsFEM [23,28,29,30,31]. Let T H be a coarse grid of domain Ω and E H be the set of all faces of the coarse grid, where H is the coarse mesh size and N e be the total number of faces. We define the neighborhood of the face E i E H by
ω i = j K j T H | E i K j .
and ω i is a union of two coarse grid blocks when E i lies in the interior of the domain Ω (see Figure 2).
We construct a multiscale space for the flux
V q , H = span { ψ i } i = 1 N q , Q B , H = { r L 2 ( Ω ) : r | K P 0 ( K ) , K T H }
and for the B, we use the space of piecewise constant functions over the coarse triangulation T H .
Here ψ i are the multiscale basis functions that supported in a coarse neighborhood ω i , N q = dim ( V q , H ) is the number of basis functions and N B = dim ( Q B , H ) is equal to the number of coarse grid cells.
For the coarse grid approximation, we use a Mixed GMsFEM and have the following variational formulation: find ( q H k + 1 , B H k + 1 ) V q , H × Q B , H such that
m ( q H k + 1 , v ) + g ( B H k + 1 , v ) = f ( v ) v V q , H , g ( q H k + 1 , r ) + c ( B H k + 1 , r ) = 0 , r Q B , H ,
We can write our system in matrix form for each nonlinear iteration
M H G H G H T C H q H B H = F H 0 ,
where
M H = R q M R q T , G H = R q G R B T , C H = R B C R B T , F H = R q F q .
Here R q and R B are the projection matrices
R q = ψ 1 , , ψ N q , R B = η 1 , , η N B .
We remark that { ψ i } i = 1 N q are the multiscale basis functions for flux and { η i } i = 1 N B are the basis functions for B.
For construction of the multiscale space for the flux, we start with the construction of the snapshot space that contains an extensive set of basis functions formed by the solution of local problems with all possible boundary conditions up to the fine grid resolution. After that, we solve a spectral problem to select dominant modes in the snapshot space.
Snapshot space. We solve the following problem on the coarse neighborhood ω i that corresponds to the coarse-grid edge E i E H : Find ( ϕ l i , η ) such that
m ( ϕ l i v ) + g ( η , v ) = 0 , v V h i , g ( ϕ l i , r ) = ω i c r d x , r Q h i ,
with boundary condition ϕ l i · n i = 0 on ω i and additional boundary condition ϕ l i · n i = δ l i on coarse edge E i , where n i is the outward unit-normal vector on ω i . The local problem (10) is solved separately on each coarse-grid element K j i ω i using the fine-grid defined in ω i by the lowest-order Raviart–Thomas element. Here E i = l = 1 J i e l , where e l are the fine-grid edges on E i / Γ P , and J i is the number of the fine grid edges on E i / Γ P . The delta function δ l i is a piecewise constant function defined on E i / Γ P such that it has value 1 on e l and value 0 on the other fine-grid edges. The constant c in (10) is chosen by compatibility condition, c = 1 | K j i | E i ϕ l i · n i d s , j = 1 , 2 .
The collection of the solutions of the above local problems generates the snapshot space in ω i
V i , snap = { ϕ l i : 1 l J i } , R i , snap = ϕ 1 i , , ϕ J i i .
Multiscale space. For reduction on the snapshot space, we solve the following local spectral problem in ω i
A i , snap ψ k i , snap = λ k S i , snap ψ k i , snap ,
where
A i , snap = [ a m n snap ] = a i ( ϕ m i , ϕ n i ) = R i , snap A i R i , snap T , S i , snap = [ m m n snap ] = s i ( ϕ m i , ϕ n i ) = R i , snap S i R i , snap T ,
and
a i ( q , v ) = E i ( q · n i ) ( v · n i ) d s , s i ( q , v ) = ω i q v d x + ω i ÷ q ÷ v d x .
To form a multiscale space, we arrange the eigenvalues in increasing order, and choose the first M i eigenvalues and take the corresponding eigenvectors ψ k i = R i , snap ψ k i , snap as the basis functions ( k = 1 , , M i )
V q , H = span { ψ k i : 1 i N e , 1 k M i } , R q = ψ 1 1 , , ψ M 1 1 , , ψ 1 N e , , ψ M N e N e .
This space will be used as the approximation space for the flux.

4. Flow Problem

For the approximation of the flow problem on the fine grid, we use a Discontinuous Galerkin method [33,34,35,36,37,38]. Let T h be a fine-grid partition of the domain Ω , E h be the set of facets in T h and E h = E i n t h E o u t h ( E i n t h is the set of interior facets and E o u t h is the set of boundary facets). We use the notations K and E to denote cell and facet in the fine grid T h . We define the jump [ u ] and the average { u } of a function u on interior facet
[ u ] E = u | K + u | K , { u } E = u | K + + u | K 2 ,
where n is the unit normal vector on E, K + and K are the two cells sharing the facet E. For E E o u t h , we define
[ u ] E = u | E , { u } E = u | E .
We have the following variational formulations for the DG approach: find ( u k + 1 , p n + 1 ) V u × Q p such that
a DG ( u k + 1 , v ) + b DG ( p k + 1 , v ) = f u ( v ) , v V u , b DG ( u k + 1 , r ) = 0 , r Q p ,
where
a DG ( u , v ) = Ω R e 1 u · v d x E E h E { R e 1 u · n } · [ v ] + { R e 1 v · n } · [ u ] E γ h R e 1 [ u ] · [ v ] d s ,
b DG ( u , p ) = K T H K p u d x + E E h E p [ u ] · n d s , f u ( v ) = Ω ( S c D 1 B k + 1 q k + 1 ) · v d x .
The fine scale velocity space V u = { v L 2 ( Ω ) : v | K ( P 1 ( T ) ) 2 , K T h } contains functions that are piecewise linear in each fine-grid element K and are continuous along the fine-grid edges, but are discontinuous across coarse grid edges. For the pressure, we use the space of piecewise constant functions.
In the matrix form, we have
A G T G 0 u k + 1 p k + 1 = F u k + 1 0
where F u k + 1 = F u ( B k + 1 , q k + 1 ) .
Next, we present a multiscale method for solutions of the Stokes problem. Similarly, to the Mixed GMsFEM to obtain multiscale basis functions, we solve problems in local domains with various boundary conditions to form a snapshot space and use a spectral problem to perform a dimension reduction.

Multiscale Approximation

Next, we present the construction of the multiscale space for the coarse scale approximation of velocity. Note that, in the DG multiscale method we construct a multiscale basis in each coarse block [26,27,39,40,41]. Let T H be a coarse-grid partition of the domain Ω with mesh size H and E H be the set of facets in T H , E H = E i n t H E o u t H . For the pressure approximation, we use the piecewise constant function space Q p , H over the coarse cells. We define V u , H as the multiscale velocity space, which contains a set of basis functions supported in each coarse block K (see Figure 3).
We construct a multiscale space for the velocity T H
V u , H = span { ψ i } i = 1 N u , Q p , H = { r L 2 ( Ω ) : r | K P 0 ( K ) , K T H }
and for the pressure, we use the space of piecewise constant functions over the coarse triangulation, N u = dim ( V u , H ) is the number of basis functions and N p = dim ( Q p , H ) is equal to the number of coarse grid cells.
For the coarse grid approximation, we use a DG approach and have the following variational formulation: find ( u H k + 1 , p H k + 1 ) V u , H × Q p , H such that
a DG ( u H k + 1 , v ) + g DG ( p H k + 1 , v ) = f ( v ) v V u , H , g ( u H k + 1 , r ) = 0 , r Q p , H ,
We can write the system in matrix form for each nonlinear iteration
A H G H G H T 0 u H p H = F H 0 ,
where
A H = R u A R u T , G H = R u G R p T , F H = R u F u .
Here R u and R p are the projection matrices
R u = ψ 1 , , ψ N u , R p = η 1 , , η N p ,
where { ψ i } i = 1 N u are the multiscale basis functions for velocity.
For construction of the multiscale space for the velocity, we start with the construction of the snapshot space that contains an extensive set of basis functions formed by the solution of local problems with all possible boundary conditions up to the fine grid resolution. After that, we solve a spectral problem to select dominant modes of the snapshot space.
Snapshot space. We construct local snapshot basis in each coarse block K i , ( i = 1 , , N ) , where N is the number of coarse blocks in Ω . The local snapshot space consists of functions which are solutions ϕ l i V u ( K i ) of
R e 1 Δ ϕ l i + η = 0 , x K i , · ϕ l i = c , x K i ,
with ϕ l i = δ i l on K i ( l = 1 , , J i ), where J i is the number of fine grid nodes on the boundary of K i , and δ i l is the discrete delta function defined on K i . This problem is solved on the fine mesh using some appropriate approximation spaces (DG). Here constant c is chosen by the compatibility condition, c = 1 | K i | K i δ i l · n d s .
Using local solutions, we form a local snapshot space in K i
V i , snap = { ϕ l i : 1 l J i } , R i , snap = ϕ 1 i , , ϕ J i i .
Multiscale space. To reduce the size of the snapshot space, we solve the following local spectral problem in the snapshot space for K i
A i , snap ψ k i , snap = λ k S i , snap ψ k i , snap ,
where
A i , snap = R i , snap A i R i , snap T , S i , snap = R i , snap S i R i , snap T ,
and A i is the matrix representation of the bilinear form a i ( u , v ) and S i is the matrix representation of the bilinear form s i ( u , v )
a i ( u , v ) = K i R e 1 u · v d x , s i ( u , v ) = K i u · v d x .
We remark that the integral in s i ( u , v ) is defined on the boundary of the coarse block. In this case, the number of the spectral problem equals the number of coarse blocks.
We arrange the eigenvalues in increasing order and choose the first eigenvectors corresponding to the first M i the smallest eigenvalues ψ k i = R i , snap ψ k i , snap as the basis functions ( k = 1 , , M i )
V i u , H = span { ψ k i : 1 i N e , 1 k M i } , R u = ψ 1 1 , , ψ M 1 1 , , ψ 1 N c , , ψ M N c N c .
This space will be used as the approximation space for the velocity. The global multiscale space V u , H is the combination of the local ones, i.e.,
V u , H = span { ϕ k i : 1 k L i , 1 i N } .

5. Numerical Results

We present numerical results for the model problem in the perforated domain. We use the Generalized Multiscale Finite Element method (GMsFEM) to construct a coarse grid approximation. The presented method based on the calculations of local multiscale basis functions in order to form a projection matrix that is used to construct a lower-dimensional approximation. After obtaining the solution of the reduced-order system, we reconstruct a fine grid resolution using a transposed projection matrix. For error calculations, we use a solution of the problem with traditional finite element approximation on the fine grid as the reference solution.
A computational perforated domain is presented in Figure 4 and has dimensions Ω = [ 0 , 1 ] × [ 0 , 1 ] . For the approximation of the perforated domain, we use an unstructured fine grid with triangular cells that resolves perforations on the grid level. This can lead to the large discrete system of equations for the traditional finite element approximations. To reduce the size of the system, we present a multiscale solver that used to construct a reduced-order approximation on the coarse grid. The coarse grid can be either structured or unstructured. In this work, for simplicity, we consider a 10 × 10 structured grid to illustrate the multiscale method. In many applications, it is easier to use structured coarse grids and resolve the subgrid effects as we do in our simulations. We construct an unstructured fine mesh that contains 61,912 cells, 93,668 facets and 31,727 vertices. The coarse grid is uniform, which contains 100 cells, 220 facets and 121 vertices. We set R e 1 = 1 , D = 10 , S c = 1 and g = 1 .
In Figure 5, Figure 6 and Figure 7, we present the results of the proposed multiscale method. The distribution of the magnitude fields (B and q) are presented in Figure 5 and Figure 6. In Figure 7, we show the velocity field. At the top of the figures, we presented the fine-scale solutions. Multiscale solutions are present on the bottom figure with D O F c = 980 for q and D O F c = 4100 for velocity field u. Note that, we use the piecewise constant functions for multiscale solution of B. We use a Picard iteration and perform three iterations for all cases with ϵ B = 10 9 .
We calculate following errors:
e u = Ω | u u m s | 2 d x Ω u 2 d x · 100 % , e q = Ω | q q m s | 2 d x Ω q 2 d x · 100 % ,
e L 2 p ¯ = Ω | p ¯ p H | 2 d x Ω p ¯ 2 d x , · 100 % , e L 2 B ¯ = Ω | B ¯ B H | 2 d x Ω B ¯ 2 d x · 100 % ,
where ( q , B ) and ( u , p ) are the reference solutions using traditional fine grid solvers, B ¯ and p ¯ are the coarse cell average for reference solutions, ( q H , B H ) and ( u H , p H ) are the coarse grid solutions using presented multiscale solvers, q m s = R q T q H and u m s = R u T u H are the reconstructed multiscale solutions.
In Table 1, we show the relative error for q and B. We presented errors for different velocities: (1) multiscale velocities with D O F c = 1100 (10 basis), (2) multiscale velocities with D O F c = 2100 (20 basis), (3) multiscale velocities with D O F c = 3100 (30 basis), (4) multiscale velocities with D O F c = 4100 (40 basis), (5) multiscale velocities with D O F c = 6100 (60 basis) and (6) fine-scale solution with D O F f = 433,384. We see the relative error for q is reduced from 4.1 % to 0.8 % , and the coarse cell average error for B is around one percent.
In Table 2 we present the relative errors for velocity and for pressure field. We presented errors for u and p with different B: (1) multiscale magnetic field with D O F c = 320 (1 basis), (2) multiscale magnetic field with D O F c = 540 (2 basis), (3) multiscale magnetic field with D O F c = 760 (3 basis), (4) multiscale magnetic field with D O F c = 980 (4 basis) and (5) fine-scale solution with D O F f = 155,580. We see the relative velocity error reduces from 60 % to 8 % with a larger number of multiscale basis functions. We use a discontinuous Galerkin Generalized Multiscale Finite Element Method (DG-GMsFEM) for the flow problem and Mixed Generalized Multiscale Finite Element Method (Mixed GMsFEM) for the magnetic field problem. The construction of basis functions in both methods follows the general framework of GMsFEM by using local snapshots and local spectral problems. The errors can vary depending on a choice of discretization. In order to further error reduction for the multiscale approximation of the flow problem, we will apply oversampling techniques in future work [39].

6. Conclusions

We considered the magnetohydrodynamics (MHD) problem in a perforated domain. A mathematical model was described using a coupled system of equations for the magnetic field and for the flow problem. We constructed a fine grid that resolved perforations on the grid level in order to use a traditional approximation. For the solution on the fine grid, we constructed approximation using the mixed finite element method. To reduce the size of the fine grid system with an accurate solution on the coarse grid, a Mixed Generalized Multiscale Finite Element Method (Mixed GMsFEM) was presented. Numerical results were presented for a two-dimensional model problem in a perforated domain and we studied the influence of the number of multiscale basis functions on the method accuracy. In the multiscale solver, the main computational gain was due to fewer multiscale basis functions, which are used to reduce the size of the system. For the traditional fine grid solutions, we have D O F f = 155,580 for the magnetic field problem and D O F f = 433,384 for the flow problem. For multiscale solvers, we have 0.9 % error for q with three multiscale basis functions with D O F c = 760 ( 0.48 % of D O F f ) and 9% error for u with 40 multiscale basis functions with D O F c = 4100 ( 0.9 % of D O F f ). In the future, we plan to consider oversampling techniques and construct a multiscale solver for three-dimensional problems.

Author Contributions

Writing—original draft, M.V.; Writing—review & editing, V.A., Q.T., E.T.C. and Y.E. All authors have read and agreed to the published version of the manuscript.

Funding

V.A.’s work is supported by the mega-grant of the Russian Federation Government N14.Y26.31.0013 and RFBR N19-31-90076. M.V.’s work is supported by the mega-grant of the Russian Federation Government N14.Y26.31.0013. The research of Qili Tang is partially supported by the National Natural Science Foundation of China (No:r 11701151), the Natural Science Foundation of Hunan Province (No: 2019JJ40279), Excellent Youth Program of Scientific Research Project of Hunan Provincial Department of Education (No: 18B064) and International Scientific and Technological Innovation Cooperation Base of Hunan Province for Computational Science (No: 2018WK4006). The research of Eric Chung is partially supported by the Hong Kong RGC General Research Fund (Project numbers 14304217 and 14302018) and CUHK Faculty of Science Direct Grant 2018-19. Y.E. would like to thank the partial support from NSF 1620318 and 1934904. Y.E. would also like to acknowledge the support of Mega-grant of the Russian Federation Government (N 14.Y26.31.0013)

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moreau, R.J. Magnetohydrodynamics; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2013; Volume 3. [Google Scholar]
  2. Davidson, P.A. An Introduction to Magnetohydrodynamics. Am. J. Phys. 2002, 70, 781. [Google Scholar] [CrossRef]
  3. Gerbeau, J.F.; Le Bris, C.; Lelièvre, T. Mathematical Methods for the Magnetohydrodynamics of Liquid Metals; Clarendon Press: London, UK, 2006. [Google Scholar]
  4. Dong, X.; He, Y.; Zhang, Y. Convergence analysis of three finite element iterative methods for the 2D/3D stationary incompressible magnetohydrodynamics. Comput. Methods Appl. Mech. Eng. 2014, 276, 287–311. [Google Scholar] [CrossRef]
  5. Dong, X.; He, Y. Convergence of some finite element iterative methods related to different Reynolds numbers for the 2D/3D stationary incompressible magnetohydrodynamics. Sci. China Math. 2016, 59, 589–608. [Google Scholar] [CrossRef]
  6. Dong, X.; He, Y. The Oseen Type Finite Element Iterative Method for the Stationary Incompressible Magnetohydrodynamics. Adv. Appl. Math. Mech. 2017, 9, 775–794. [Google Scholar] [CrossRef]
  7. Dong, X.; He, Y.; Wei, H.; Zhang, Y. Local and parallel finite element algorithm based on the partition of unity method for the incompressible MHD flow. Adv. Comput. Math. 2018, 44, 1295–1319. [Google Scholar] [CrossRef]
  8. Tang, Q.; Huang, Y. Local and parallel finite element algorithm based on Oseen-type iteration for the stationary incompressible MHD flow. J. Sci. Comput. 2017, 70, 149–174. [Google Scholar] [CrossRef]
  9. Tang, Q.; Huang, Y. Analysis of Local and Parallel Algorithm for Incompressible Magnetohydrodynamics Flows by Finite Element Iterative Method. Commun. Comput. Phys. 2019, 25, 729–751. [Google Scholar] [CrossRef]
  10. Badia, S.; Codina, R.; Planas, R. On an unconditionally convergent stabilized finite element approximation of resistive magnetohydrodynamics. J. Comput. Phys. 2013, 234, 399–416. [Google Scholar] [CrossRef] [Green Version]
  11. Greif, C.; Li, D.; Schötzau, D.; Wei, X. A mixed finite element method with exactly divergence-free velocities for incompressible magnetohydrodynamics. Comput. Methods Appl. Mech. Eng. 2010, 199, 2840–2855. [Google Scholar] [CrossRef]
  12. Hiptmair, R.; Li, L.; Mao, S.; Zheng, W. A fully divergence-free finite element method for magnetohydrodynamic equations. Math. Models Methods Appl. Sci. 2018, 28, 659–695. [Google Scholar] [CrossRef] [Green Version]
  13. Hu, K.; Ma, Y.; Xu, J. Stable finite element methods preserving ∇ · {B} = 0 exactly for MHD models. Numer. Math. 2017, 135, 371–396. [Google Scholar] [CrossRef]
  14. Jikov, V.V.; Kozlov, S.M.; Oleinik, O.A. Homogenization of Differential Operators and Integral Functionals; Springer: Berlin, Germany, 1991. [Google Scholar]
  15. Bakhvalov, N.; Panasenko, G. Homogenization in Periodic Media, Mathematical Problems of the Mechanics of Composite Materials; Nauka: Moscow, Russia, 1984. [Google Scholar]
  16. Allaire, G.; Brizzi, R. A multiscale finite element method for numerical homogenization. SIAM J. Multiscale Modeling Simul. 2005, 4, 790–812. [Google Scholar] [CrossRef] [Green Version]
  17. Efendiev, Y.; Hou, T. Multiscale Finite Element Methods: Theory and Applications. In Surveys and Tutorials in the Applied Mathematical Sciences; Springer: New York, NY, USA, 2009; Volume 4. [Google Scholar]
  18. Efendiev, Y.; Galvis, J.; Hou, T.Y. Generalized multiscale finite element methods (GMsFEM). J. Comput. Phys. 2013, 251, 116–135. [Google Scholar] [CrossRef] [Green Version]
  19. Ohlberger, M.; Schindler, F. Error control for the localized reduced basis multi-scale method with adaptive on-line enrichment. arXiv 2015, arXiv:1501.05202. [Google Scholar]
  20. Muljadi, B.; Narski, J.; Lozinski, A.; Degond, P. Non-conforming multiscale finite element method for Stokes flows in heterogeneous media. Part I: Methodologies and numerical experiments. arXiv 2014, arXiv:1404.2837. [Google Scholar]
  21. Le Bris, L.; Legoll, F.; Lozinski, A. An MsFEM type approach for perforated domains. Multiscale Model. Simul. 2014, 12, 1046–1077. [Google Scholar] [CrossRef] [Green Version]
  22. Henning, P.; Ohlberger, M. The heterogeneous multiscale finite element method for elliptic homogenization problems in perforated domains. Numer. Math. 2009, 113, 601–629. [Google Scholar] [CrossRef]
  23. Chen, Z.; Hou, T. A mixed multiscale finite element method for elliptic problems with oscillating coefficients. Math. Comput. 2002, 72, 541–576. [Google Scholar] [CrossRef] [Green Version]
  24. Brown, D.L.; Peterseim, D. A multiscale method for porous microstructures. arXiv 2014, arXiv:1411.1944. [Google Scholar] [CrossRef] [Green Version]
  25. Efendiev, Y.; Galvis, J.; Wu, X.H. Multiscale finite element methods for high-contrast problems using local spectral basis functions. J. Comput. Phys. 2011, 230, 937–955. [Google Scholar] [CrossRef]
  26. Chung, E.T.; Efendiev, Y.; Li, G.; Vasilyeva, M. Generalized multiscale finite element methods for problems in perforated heterogeneous domains. Appl. Analys. 2016, 95, 2254–2279. [Google Scholar] [CrossRef] [Green Version]
  27. Chung, E.T.; Efendiev, Y.; Vasilyeva, M.; Wang, Y. A multiscale discontinuous Galerkin method in perforated domains. In Proceedings of the Institute of Mathematics and Mechanics, Maringá, Brazil, 22–23 September 2016; Volume 42, pp. 212–229. [Google Scholar]
  28. Chung, E.T.; Efendiev, Y.; Lee, C.S. Mixed generalized multiscale finite element methods and applications. Multiscale Model. Simul. 2015, 13, 338–366. [Google Scholar] [CrossRef] [Green Version]
  29. Chan, H.Y.; Chung, E.; Efendiev, Y. Adaptive mixed GMsFEM for flows in heterogeneous media. Numer. Math. Theory Methods Appl. 2016, 9, 497–527. [Google Scholar] [CrossRef] [Green Version]
  30. Chung, E.T.; Leung, W.T.; Vasilyeva, M. Mixed GMsFEM for second order elliptic problem in perforated domains. J. Comput. Appl. Math. 2016, 304, 84–99. [Google Scholar] [CrossRef]
  31. Chung, E.T.; Leung, W.T.; Vasilyeva, M.; Wang, Y. Multiscale model reduction for transport and flow problems in perforated domains. J. Comput. Appl. Math. 2018, 330, 519–535. [Google Scholar] [CrossRef]
  32. Gunzburger, M.D.; Meir, A.J.; Peterson, J.S. On the existence, uniqueness, and finite element approximation of solutions of the equations of stationary, incompressible magnetohydrodynamics. Math. Comput. 1991, 56, 523–563. [Google Scholar] [CrossRef]
  33. Arnold, D.; Brezzi, F.; Cockburn, B.; Marini, L. Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal. 2002, 39, 1749–1779. [Google Scholar] [CrossRef]
  34. Ewing, R.; Wang, J.; Yang, Y. A stabilized discontinuous finite element method for elliptic problems. Numer. Linear Algebra Appl. 2003, 10, 83–104. [Google Scholar] [CrossRef]
  35. Lazarov, R.; Pasciak, J.; Schöberl, J.; Vassilevski, P. Almost optimal interior penalty discontinuous approximations of symmetric elliptic problems on non-matching grids. Numer. Math. 2003, 96, 295–315. [Google Scholar] [CrossRef]
  36. Lazarov, R.; Tomov, S.; Vassilevski, P. Interior penalty discontinuous approximations of elliptic problems. Comput. Methods Appl. Math. 2001, 1, 367–382. [Google Scholar] [CrossRef] [Green Version]
  37. Riviere, B. Discontinuous Galerkin Methods for Solving Elliptic And Parabolic Equations: Theory and Implementation; SIAM: Philadelphia, PA, USA, 2008. [Google Scholar]
  38. Girault, V.; Riviére, B.; Wheeler, M. A discontinuous Galerkin method with nonoverlapping domain decomposition for the Stokes and Navier-Stokes problems. Math. Comput. 2005, 74, 53–84. [Google Scholar] [CrossRef] [Green Version]
  39. Chung, E.T.; Vasilyeva, M.; Wang, Y. A conservative local multiscale model reduction technique for Stokes flows in heterogeneous perforated domains. J. Comput. Appl. Math. 2017, 321, 389–405. [Google Scholar] [CrossRef] [Green Version]
  40. Chung, E.T.; Efendiev, Y.; Leung, W.T. An adaptive generalized multiscale discontinuous Galerkin method for high-contrast flow problems. Multiscale Model. Simul. 2018, 16, 1227–1257. [Google Scholar] [CrossRef] [Green Version]
  41. Chung, E.; Leung, W.T. A sub-grid structure enhanced discontinuous Galerkin method for multiscale diffusion and convection-diffusion problems. Commun. Comput. Phys. 2013, 14, 370–392. [Google Scholar] [CrossRef]
Figure 1. Illustration of a perforated domain.
Figure 1. Illustration of a perforated domain.
Computation 08 00058 g001
Figure 2. Illustration of a coarse grid and local domain ( ω i ) for the Magnetic field problem.
Figure 2. Illustration of a coarse grid and local domain ( ω i ) for the Magnetic field problem.
Computation 08 00058 g002
Figure 3. Illustration of a coarse grid and local domain ( K i ) for the Flow problem.
Figure 3. Illustration of a coarse grid and local domain ( K i ) for the Flow problem.
Computation 08 00058 g003
Figure 4. Perforated domain with coarse (left) and fine grids (right).
Figure 4. Perforated domain with coarse (left) and fine grids (right).
Computation 08 00058 g004
Figure 5. The distribution of the magnetic field. (Left): Reference (fine grid) solution, D O F f = 155,580. (Center): Cell average for reference (fine grid) solution. (Right): Multiscale solution, D O F c = 980 (0.48% from D O F f ).
Figure 5. The distribution of the magnetic field. (Left): Reference (fine grid) solution, D O F f = 155,580. (Center): Cell average for reference (fine grid) solution. (Right): Multiscale solution, D O F c = 980 (0.48% from D O F f ).
Computation 08 00058 g005
Figure 6. The distribution of the flux. Magnitude of flux, x and y components of flux, q, q x , q y (from left to right). (Top): Reference (fine grid) solution, D O F f = 155,580. (Bottom): Multiscale solution, D O F c = 980 (0.48% from D O F f ).
Figure 6. The distribution of the flux. Magnitude of flux, x and y components of flux, q, q x , q y (from left to right). (Top): Reference (fine grid) solution, D O F f = 155,580. (Bottom): Multiscale solution, D O F c = 980 (0.48% from D O F f ).
Computation 08 00058 g006
Figure 7. The distribution of the velocity. Magnitude of velocity, x and y components of velocity, u, u x u y (from left to right). (Top): Reference (fine grid) solution, D O F f = 433,384. (Bottom): Multiscale solution, D O F c = 4100 (0.9% from D O F f ).
Figure 7. The distribution of the velocity. Magnitude of velocity, x and y components of velocity, u, u x u y (from left to right). (Top): Reference (fine grid) solution, D O F f = 433,384. (Bottom): Multiscale solution, D O F c = 4100 (0.9% from D O F f ).
Computation 08 00058 g007
Table 1. Relative errors for flux and magnetic field with a different number of multiscale basis functions for 10 × 10 meshes. D O F f = 155,580.
Table 1. Relative errors for flux and magnetic field with a different number of multiscale basis functions for 10 × 10 meshes. D O F f = 155,580.
DOF c Flux, e q (%)Magnetic Field, e B ¯ (%)
u-multiscale solution with D O F c = 1100
3204.1701.690
5401.1811.578
7600.9581.578
9800.8351.578
u-multiscale solution with D O F c = 2100
3204.1651.688
5401.1761.577
7600.9531.577
9800.8311.577
u-multiscale solution with D O F c = 3100
3204.1641.688
5401.1761.577
7600.9531.577
9800.8301.577
u-multiscale solution with D O F c = 4100
3204.1641.688
5401.1761.577
7600.9521.577
9800.8301.577
u-multiscale solution with D O F c = 6100
3204.1631.688
5401.1761.577
7600.9521.577
9800.8301.577
u-fine-scale solution with D O F f = 433,384
3204.1611.688
5401.1741.577
7600.9541.577
9800.8321.577
Table 2. Relative errors for velocity and pressure with different number of multiscale basis functions for 10 × 10 meshes. D O F f = 433,384.
Table 2. Relative errors for velocity and pressure with different number of multiscale basis functions for 10 × 10 meshes. D O F f = 433,384.
DOF c Velocity, e u (%)Pressure, e p ¯ (%)
B-multiscale solution with D O F c = 320
110060.50125.800
210016.6892.005
310011.2811.226
41009.7921.138
61008.8431.078
B-multiscale solution with D O F c = 540
110060.47725.798
210016.6941.996
310011.2821.225
41009.7921.138
61008.8431.078
B-multiscale solution with D O F c = 760
110060.47325.797
210016.6951.997
310011.2821.225
41009.7921.138
61008.8431.078
B-multiscale solution with D O F c = 980
110060.47225.797
210016.6941.997
310011.2821.225
41009.7921.138
61008.8431.078
B-fine-scale solution with D O F f = 155,580
110059.42025.300
210016.6151.878
310011.2381.199
41009.7541.112
61008.8091.055

Share and Cite

MDPI and ACS Style

Alekseev, V.; Tang, Q.; Vasilyeva, M.; Chung, E.T.; Efendiev, Y. Mixed Generalized Multiscale Finite Element Method for a Simplified Magnetohydrodynamics Problem in Perforated Domains. Computation 2020, 8, 58. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020058

AMA Style

Alekseev V, Tang Q, Vasilyeva M, Chung ET, Efendiev Y. Mixed Generalized Multiscale Finite Element Method for a Simplified Magnetohydrodynamics Problem in Perforated Domains. Computation. 2020; 8(2):58. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020058

Chicago/Turabian Style

Alekseev, Valentin, Qili Tang, Maria Vasilyeva, Eric T. Chung, and Yalchin Efendiev. 2020. "Mixed Generalized Multiscale Finite Element Method for a Simplified Magnetohydrodynamics Problem in Perforated Domains" Computation 8, no. 2: 58. https://0-doi-org.brum.beds.ac.uk/10.3390/computation8020058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop