In this paper, we propose a local-global multiscale method for highly heterogeneous stochastic groundwater flow problems under the framework of reduced basis method and the generalized multiscale finite element method (GMsFEM). Due to incomplete char
acterization of the medium properties of the groundwater flow problems, random variables are used to parameterize the uncertainty. As a result, solving the problem repeatedly is required to obtain statistical quantities. Besides, the medium properties are usually highly heterogeneous, which will result in a large linear system that needs to be solved. Therefore, it is intrinsically inevitable to seek a computational-efficient model reduction method to overcome the difficulty. We will explore the combination of the reduced basis method and the GMsFEM. In particular, we will use residual-driven basis functions, which are key ingredients in GMsFEM. This local-global multiscale method is more efficient than applying the GMsFEM or reduced basis method individually. We first construct parameter-independent multiscale basis functions that include both local and global information of the permeability fields, and then use these basis functions to construct several global snapshots and global basis functions for fast online computation with different parameter inputs. We provide rigorous analysis of the proposed method and extensive numerical examples to demonstrate the accuracy and efficiency of the local-global multiscale method.
In this paper, we systemically review and compare two mixed multiscale finite element methods (MMsFEM) for multiphase transport in highly heterogeneous media. In particular, we will consider the mixed multiscale finite element method using limited gl
obal information, simply denoted by MMsFEM, and the mixed generalized multiscale finite element method (MGMsFEM) with residual driven online multiscale basis functions. Both methods are under the framework of mixed multiscale finite element methods, where the pressure equation is solved in the coarse grid with carefully constructed multiscale basis functions for the velocity. The multiscale basis functions in both methods include local and global media information. In terms of MsFEM using limited global information, only one multiscale basis function is utilized in each local neighborhood while multiple basis are used in MGMsFEM. We will test and compare these two methods using the benchmark three-dimensional SPE10 model. A range of coarse grid sizes and different combinations of basis functions (offline and online) will be considered with CPU time reported for each case. In our numerical experiments, we observe good accuracy by the two above methods. Finally, we will discuss and compare the advantages and disadvantages of the two methods in terms of accuracy and computational costs.
In this paper, we propose a coupled Discrete Empirical Interpolation Method (DEIM) and Generalized Multiscale Finite element method (GMsFEM) to solve nonlinear parabolic equations with application to the Allen-Cahn equation. The Allen-Cahn equation i
s a model for nonlinear reaction-diffusion process. It is often used to model interface motion in time, e.g. phase separation in alloys. The GMsFEM allows solving multiscale problems at a reduced computational cost by constructing a reduced-order representation of the solution on a coarse grid. In arXiv:1301.2866, it was shown that the GMsFEM provides a flexible tool to solve multiscale problems by constructing appropriate snapshot, offline and online spaces. In this paper, we solve a time dependent problem, where online enrichment is used. The main contribution is comparing different online enrichment methods. More specifically, we compare uniform online enrichment and adaptive methods. We also compare two kinds of adaptive methods. Furthermore, we use DEIM, a dimension reduction method to reduce the complexity when we evaluate the nonlinear terms. Our results show that DEIM can approximate the nonlinear term without significantly increasing the error. Finally, we apply our proposed method to the Allen Cahn equation.
In this paper, we consider an online enrichment procedure using the Generalized Multiscale Finite Element Method (GMsFEM) in the context of a two-phase flow model in heterogeneous porous media. The coefficient of the elliptic equation is referred to
as the permeability and is the main source of heterogeneity within the model. The elliptic pressure equation is solved using online GMsFEM, and is coupled with a hyperbolic transport equation where local conservation of mass is necessary. To satisfy the conservation property, we aim at constructing conservative fluxes within the space of multiscale basis functions through the use of a postprocessing technique. In order to improve the accuracy of the pressure and velocity solutions in the online GMsFEM we apply a systematic online enrichment procedure. The increase in pressure accuracy due to the online construction is inherited by the conservative flux fields and the desired saturation solutions from the coupled transport equation. Despite the fact that the coefficient of the pressure equation is dependent on the saturation which may vary in time, we may construct an approximation space using the initial coefficient where no further basis updates follow. Numerical results corresponding to four different types of heterogeneous permeability coefficients are exhibited to test the proposed methodology.
We consider nonlinear elastic wave equations generalizing Goldbergs five constants model. We analyze the nonlinear interaction of two distorted plane waves and characterize the possible nonlinear responses. Using the boundary measurements of the nonl
inear responses, we solve the inverse problem of determining elastic parameters from the displacement-to-traction map.
For scalar semilinear wave equations, we analyze the interaction of two (distorted) plane waves at an interface between media of different nonlinear properties. We show that new waves are generated from the nonlinear interactions, which might be resp
onsible for the observed nonlinear effects in applications. Also, we show that the incident waves and the nonlinear responses determine the location of the interface and some information of the nonlinear properties of the media. In particular, for the case of a jump discontinuity at the interface, we can determine the magnitude of the jump.
