BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Eric Chung

Eric Chung contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.APComputational Engineering, Finance, and Sciencemath-phmath.MP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

A deep learning based reduced order modeling for stochastic underground flow problems

In this paper, we propose a deep learning based reduced order modeling method for stochastic underground flow problems in highly heterogeneous media. We aim to utilize supervised learning to build a reduced surrogate model from the stochastic parameter space that characterizes the possible highly heterogeneous media to the solution space of a stochastic flow problem to have fast online simulations. Dominant POD modes obtained from a well-designed spectral problem in a global snapshot space are used to represent the solution of the flow problem. Due to the small dimension of the solution, the complexity of the neural network is significantly reduced. We adopt the generalized multiscale finite element method (GMsFEM), in which a set of local multiscale basis functions that can capture the heterogeneity of the media and source information are constructed to efficiently generate globally defined snapshot space. Rigorous theoretical analyses are provided and extensive numerical experiments for linear and nonlinear stochastic flows are provided to verify the superior performance of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A discontinuous Galerkin based multiscale method for heterogeneous elastic wave equations

In this paper, we develop a local multiscale model reduction strategy for the elastic wave equation in strongly heterogeneous media, which is achieved by solving the problem in a coarse mesh with multiscale basis functions. We use the interior penalty discontinuous Galerkin (IPDG) to couple the multiscale basis functions that contain important heterogeneous media information. The construction of efficient multiscale basis functions starts with extracting dominant modes of carefully defined spectral problems to represent important media feature, which is followed by solving a constraint energy minimization problems. Then a Petrov-Galerkin projection and systematization onto the coarse grid is applied. As a result, an explicit and energy conserving scheme is obtained for fast online simulation. The method exhibits both coarse-mesh and spectral convergence as long as one appropriately chose the oversampling size. We rigorously analyze the stability and convergence of the proposed method. Numerical results are provided to show the performance of the multiscale method and confirm the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraint energy minimizing generalized multiscale finite element method for convection diffusion equation

In this paper we present and analyze a constraint energy minimizing generalized multiscale finite element method for convection diffusion equation. To define the multiscale basis functions, we first build an auxiliary multiscale space by solving local spectral problems motivated by analysis. Then constraint energy minimization performed in oversampling domains is exploited to construct the multiscale space. The resulting multiscale basis functions have a good decay property even for high contrast diffusion and convection coefficients. Furthermore, if the number of oversampling layer is chosen properly, we can prove that the convergence rate is proportional to the coarse mesh size. Our analysis also indicates that the size of the oversampling domain weakly depends on the contrast of the heterogeneous coefficients. Several numerical experiments are presented illustrating the performances of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Generalized multiscale finite element method for highly heterogeneous compressible flow

In this paper, we study the generalized multiscale finite element method (GMsFEM) for single phase compressible flow in highly heterogeneous porous media. We follow the major steps of the GMsFEM to construct permeability dependent offline basis for fast coarse-grid simulation. The offline coarse space is efficiently constructed only once based on the initial permeability field with parallel computing. A rigorous convergence analysis is performed for two types of snapshot spaces. The analysis indicates that the convergence rates of the proposed multiscale method depend on the coarse meshsize and the eigenvalue decay of the local spectral problem. To further increase the accuracy of multiscale method, residual driven online multiscale basis is added to the offline space. The construction of online multiscale basis is based on a carefully design error indicator motivated by the analysis. We find that online basis is particularly important for the singular source. Rich numerical tests on typical 3D highly heterogeneous medias are presented to demonstrate the impressive computational advantages of the proposed multiscale method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A local-global generalized multiscale finite element method for highly heterogeneous stochastic groundwater flow problems

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 characterization 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Iterated numerical homogenization for multi-scale elliptic equations with monotone nonlinearity

Nonlinear multi-scale problems are ubiquitous in materials science and biology. Complicated interactions between nonlinearities and (nonseparable) multiple scales pose a major challenge for analysis and simulation. In this paper, we study the numerical homogenization for multi-scale elliptic PDEs with monotone nonlinearity, in particular the Leray-Lions problem (a prototypical example is the p-Laplacian equation), where the nonlinearity cannot be parameterized with low dimensional parameters, and the linearization error is non-negligible. We develop the iterated numerical homogenization scheme by combining numerical homogenization methods for linear equations, and the so-called "quasi-norm" based iterative approach for monotone nonlinear equation. We propose a residual regularized nonlinear iterative method, and in addition, develop the sparse updating method for the efficient update of coarse spaces. A number of numerical results are presented to complement the analysis and valid the numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A pressure robust staggered discontinuous Galerkin method for the Stokes equations

In this paper we propose a pressure robust staggered discontinuous Galerkin method for the Stokes equations on general polygonal meshes by using piecewise constant approximations. We modify the right hand side of the body force in the discrete formulation by exploiting divergence preserving velocity reconstruction operator, which is the crux for pressure independent velocity error estimates. The optimal convergence for velocity gradient, velocity and pressure are proved. In addition, we are able to prove the superconvergence of velocity approximation by the incorporation of divergence preserving velocity reconstruction operator in the dual problem, which is also an important contribution of this paper. Finally, several numerical experiments are carried out to confirm the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A uniformly robust staggered DG method for the unsteady Darcy-Forchheimer-Brinkman problem

In this paper we propose and analyze a uniformly robust staggered DG method for the unsteady Darcy-Forchheimer-Brinkman problem. Our formulation is based on velocity gradient-velocity-pressure and the resulting scheme can be flexibly applied to fairly general polygonal meshes. We relax the tangential continuity for velocity, which is the key ingredient in achieving the uniform robustness. We present well-posedness and error analysis for both the semi-discrete scheme and the fully discrete scheme, and the theories indicate that the error estimates for velocity are independent of pressure. Several numerical experiments are presented to confirm the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive multiscale model reduction for nonlinear parabolic equations using GMsFEM

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 is 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive staggered DG method for Darcy flows in fractured porous media

Modeling flows in fractured porous media is important in applications. One main challenge in numerical simulation is that the flow is strongly influenced by the fractures, so that the solutions typically contain complex features, which require high computational grid resolutions. Instead of using uniformly fine mesh, a more computationally efficient adaptively refined mesh is desirable. In this paper we design and analyze a novel residual-type a posteriori error estimator for staggered DG methods on general polygonal meshes for Darcy flows in fractured porous media. The method can handle fairly general meshes and hanging nodes can be simply incorporated into the construction of the method, which is highly appreciated for adaptive mesh refinement. The reliability and efficiency of the error estmator are proved. The derivation of the reliability hinges on the stability of the continuous setting in the primal formulation. A conforming counterpart that is continuous within each bulk domain for the discrete bulk pressure is defined to facilitate the derivation of the reliability. Finally, several numerical experiments including multiple non-intersecting fractures are carried out to confirm the proposed theories.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of the CEM-GMsFEM for Stokes flows in heterogeneous perforated domains

In this paper, we consider the incompressible Stokes flow problem in a perforated domain and employ the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) to solve this problem. The proposed method provides a flexible and systematical approach to construct crucial divergence-free multiscale basis functions for approximating the displacement field. These basis functions are constructed by solving a class of local energy minimization problems over the eigenspaces that contain local information on the heterogeneities. These multiscale basis functions are shown to have the property of exponential decay outside the corresponding local oversampling regions. By adapting the technique of oversampling, the spectral convergence of the method with error bounds related to the coarse mesh size is proved.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Learning Algorithms for Coarsening Uncertainty Space and Applications to Multiscale Simulations

In this paper, we investigate and design multiscale simulations for stochastic multiscale PDEs. As for the space, we consider a coarse grid and a known multiscale method, the Generalized Multiscale Finite Element Method (GMsFEM). In order to obtain a small dimensional representation of the solution in each coarse block, the uncertainty space needs to be partitioned (coarsened). This coarsening collects realizations that provide similar multiscale features as outlined in GMsFEM (or other method of choice). This step is known to be computationally demanding as it requires many local solves and clustering based on them. In this paper, we take a different approach and learn coarsening the uncertainty space. Our methods use deep learning techniques in identifying clusters(coarsening) in the uncertainty space. We use convolutional neural networks combined with some techniques in adversary neural networks. We define appropriate loss functions in the proposed neural networks, where the loss function is composed of several parts that includes terms related to clusters and reconstruction of basis functions. We present numerical results for channelized permeability fields in the examples of flows in porous media.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mixed GMsFEM for linear poroelasticity problems in heterogeneous porous media

Accurate numerical simulations of interaction between fluid and solid play an important role in applications. The task is challenging in practical scenarios as the media are usually highly heterogeneous with very large contrast. To overcome this computational challenge, various multiscale methods are developed. In this paper, we consider a class of linear poroelasticity problems in high contrast heterogeneous porous media, and develop a mixed generalized multiscale finite element method (GMsFEM) to obtain a fast computational method. Our aim is to develop a multiscale method that is robust with respect to the heterogeneities and contrast of the media, and gives a mass conservative fluid velocity field. We will construct decoupled multiscale basis functions for the elastic displacement as well as fluid velocity. Our multiscale basis functions are local. The construction is based on some suitable choices of local snapshot spaces and local spectral decomposition, with the goal of extracting dominant modes of the solutions. For the pressure, we will use piecewise constant approximation. We will present several numerical examples to illustrate the performance of our method. Our results indicate that the proposed method is able to give accurate numerical solutions with a small degree of freedoms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Online conservative generalized multiscale finite element method for flow models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Staggered DG method with small edges for Darcy flows in fractured porous media

In this paper, we present and analyze a staggered discontinuous Galerkin method for Darcy flows in fractured porous media on fairly general meshes. A staggered discontinuous Galerkin method and a standard conforming finite element method with appropriate inclusion of interface conditions are exploited for the bulk region and the fracture, respectively. Our current analysis weakens the usual assumption on the polygonal mesh, which can integrate more general meshes such as elements with arbitrarily small edges into our theoretical framework. We prove the optimal convergence estimates in $L^2$ error for all the variables by exploiting the Ritz projection. Importantly, our error estimates are shown to be fully robust with respect to the heterogeneity and anisotropy of the permeability coefficients. Several numerical experiments including meshes with small edges and anisotropic meshes are carried out to confirm the theoretical findings. Finally, our method is applied in the framework of unfitted mesh.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A new staggered DG method for the Brinkman problem robust in the Darcy and Stokes limits

In this paper we propose a novel staggered discontinuous Galerkin method for the Brinkman problem on general quadrilateral and polygonal meshes. The proposed method is robust in the Stokes and Darcy limits, in addition, hanging nodes can be automatically incorporated in the construction of the method, which are desirable features in practical applications. There are three unknowns involved in our formulation, namely velocity gradient, velocity and pressure. Unlike the original staggered DG formulation proposed for the Stokes equations in \cite{KimChung13}, we relax the tangential continuity of velocity and enforce different staggered continuity properties for the three unknowns, which is tailored to yield an optimal $L^2$ error estimates for velocity gradient, velocity and pressure independent of the viscosity coefficient. Moreover, by choosing suitable projection, superconvergence can be proved for $L^2$ error of velocity. Finally, several numerical results illustrating the good performances of the proposed method and confirming the theoretical findings are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Constraint energy minimizing generalized multiscale finite element method for nonlinear poroelasticity and elasticity

In this paper, we apply the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM) to first solving a nonlinear poroelasticity problem. The arising system consists of a nonlinear pressure equation and a nonlinear stress equation in strain-limiting setting, where strains keep bounded while stresses can grow arbitrarily large. After time discretization of the system, to tackle the nonlinearity, we linearize the resulting equations by Picard iteration. To handle the linearized equations, we employ the CEM-GMsFEM and obtain appropriate offline multiscale basis functions for the pressure and the displacement. More specifically, first, auxiliary multiscale basis functions are generated by solving local spectral problems, via the GMsFEM. Then, multiscale spaces are constructed in oversampled regions, by solving a constraint energy minimizing (CEM) problem. After that, this strategy (with the CEM-GMsFEM) is also applied to a static case of the above nonlinear poroelasticity problem, that is, elasticity problem, where the residual based online multiscale basis functions are generated by an adaptive enrichment procedure, to further reduce the error. Convergence of the two cases is demonstrated by several numerical simulations, which give accurate solutions, with converging coarse-mesh sizes as well as few basis functions (degrees of freedom) and oversampling layers.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Generalized multiscale finite element method for a strain-limiting nonlinear elasticity model

In this paper, we consider multiscale methods for nonlinear elasticity. In particular, we investigate the Generalized Multiscale Finite Element Method (GMsFEM) for a strain-limiting elasticity problem. Being a special case of the naturally implicit constitutive theory of nonlinear elasticity, strain-limiting relation has presented an interesting class of material bodies, for which strains remain bounded (even infinitesimal) while stresses can become arbitrarily large. The nonlinearity and material heterogeneities can create multiscale features in the solution, and multiscale methods are therefore necessary. To handle the resulting nonlinear monotone quasilinear elliptic equation, we use linearization based on the Picard iteration. We consider two types of basis functions, offline and online basis functions, following the general framework of GMsFEM. The offline basis functions depend nonlinearly on the solution. Thus, we design an indicator function and we will recompute the offline basis functions when the indicator function predicts that the material property has significant change during the iterations. On the other hand, we will use the residual based online basis functions to reduce the error substantially when updating basis functions is necessary. Our numerical results show that the above combination of offline and online basis functions is able to give accurate solutions with only a few basis functions per each coarse region and updating basis functions in selected iterations.

0 citations0 reviewsNo code linkedOpen access