Search

We propose a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution., Comment: 32 pages, 8 figures

Scientific and engineering problems often involve parametric partial differential equations (PDEs), such as uncertainty quantification, optimizations, and inverse problems. However, solving these PDEs repeatedly can be prohibitively expensive, especially for large-scale complex applications. To address this issue, reduced order modeling (ROM) has emerged as an effective method to reduce computational costs. However, ROM often requires significant modifications to the existing code, which can be time-consuming and complex, particularly for large-scale legacy codes. Non-intrusive methods have gained attention as an alternative approach. However, most existing non-intrusive approaches are purely data-driven and may not respect the underlying physics laws during the online stage, resulting in less accurate approximations of the reduced solution. In this study, we propose a new non-intrusive bi-fidelity reduced basis method for time-independent parametric PDEs. Our algorithm utilizes the discrete operator, solutions, and right-hand sides obtained from the high-fidelity legacy solver. By leveraging a low-fidelity model, we efficiently construct the reduced operator and right-hand side for new parameter values during the online stage. Unlike other non-intrusive ROM methods, we enforce the reduced equation during the online stage. In addition, the non-intrusive nature of our algorithm makes it straightforward and applicable to general nonlinear time-independent problems. We demonstrate its performance through several benchmark examples, including nonlinear and multiscale PDEs.

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced-order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size alone and that only a few oversampling layers, as well as basis functions, are needed., Comment: 22 pages, 7 figures, 4 tables, submitted to Journal of Computational Physics, fixed some typos and notation

Multiscale elliptic equations with scale separation are often approximated by the corresponding homogenized equations with slowly varying homogenized coefficients (the G-limit). The traditional homogenization techniques typically rely on the periodicity of the multiscale coefficients, thus finding the G-limits often requires sophisticated techniques in more general settings even when multiscale coefficient is known, if possible. Alternatively, we propose a simple approach to estimate the G-limits from (noisy-free or noisy) multiscale solution data, either from the existing forward multiscale solvers or sensor measurements. By casting this problem into an inverse problem, our approach adopts physics-informed neural networks (PINNs) algorithm to estimate the G-limits from the multiscale solution data by leveraging a priori knowledge of the underlying homogenized equations. Unlike the existing approaches, our approach does not rely on the periodicity assumption or the known multiscale coefficient during the learning stage, allowing us to estimate homogenized coefficients in more general settings beyond the periodic setting. We demonstrate that the proposed approach can deliver reasonable and accurate approximations to the G-limits as well as homogenized solutions through several benchmark problems.

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization's effective coefficients are computed through solving the arising cell problems. To tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method., Comment: 45 pages, 4 figures, 2 tables, major revision. This is the accepted manuscript by Journal of Computational and Applied Mathematics (2021). The published journal article is available at https://doi.org/10.1016/j.cam.2021.113648 (2021)

We study homogenization of a locally periodic two-scale dual-continuum system where each continuum interacts with the other. Equations for each continuum are written separately with interaction terms (exchange terms) added. The homogenization limit depends strongly on the scale of this continuum interaction term with respect to the microscopic scale. In J. S. R. Park and V. H. Hoang, {\it Hierarchical multiscale finite element method for multicontinuum media}, arXiv:1906.04635, we study in details the case where the interaction terms are scaled as $O(1/\ep^2)$ where $\ep$ is the microscale of the problem. We establish rigorously homogenization limit for this case where we show that in the homogenization limit, the dual-continuum structure disappears. In this paper, we consider the case where this term is scaled as $O(1/\ep)$. This case is far more interesting and difficult as the homogenized problem is a dual-continuum system which contains features that are not in the original two scale problem. In particular, the homogenized dual-continuum system contains extra convection terms and negative interaction coefficients while the interaction coefficient between the continua in the original two scale system obtains both positive and negative values. We prove rigorously the homogenization convergence. We also derive rigorously a homogenization convergence rate. Homogenization of dual-continuum system of this type has not been considered before.

We consider in this paper a challenging problem of simulating fluid flows, in complex multiscale media possessing multi-continuum background. As an effort to handle this obstacle, model reduction is employed. In \cite{rh2}, homogenization was nicely applied, to find effective coefficients and homogenized equations (for fluid flow pressures) of a dual-continuum system, with new convection terms and negative interaction coefficients. However, some degree of multiscale still remains. This motivates us to propose the generalized multiscale finite element method (GMsFEM), which is coupled with the dual-continuum homogenized equations, toward speeding up the simulation, improving the accuracy as well as clearly representing the interactions between the dual continua. In our paper, globally, each continuum is viewed as a system and connected to the other throughout the domain. We take into consideration the flow transfers between the dual continua and within each continuum itself. Such multiscale flow dynamics are modeled by the GMsFEM, which systematically generates either uncoupled or coupled multiscale basis (to carry the local characteristics to the global ones), via establishing local snapshots and spectral decomposition in the snapshot space. As a result, we will work with a system of two equations coupled with some interaction terms, and each equation describes one of the dual continua on the fine grid. Convergence analysis of the proposed GMsFEM is accompanied with the numerical results, which support the favorable outcomes., Comment: 35 pages, 6 figures, 4 tables, submitted to Journal of Computational and Applied Mathematics

Simulation in media with multiple continua where each continuum interacts with every other is often challenging due to multiple scales and high contrast. One needs some types of model reduction. One of the approaches is multi-continuum technique, where every process in each continuum is modeled separately and an interaction term is added. Direct numerical simulation in multi scale media is usually not practicable. For this reason, one constructs the corresponding homogenized equations. The paper develops a hierarchical approach for solving cell problems at a dense network of macroscopic points with an essentially optimal computation cost. The method employs the fact that neighboring representative volume elements (RVEs) share similar features; and effective properties of the neighboring RVEs are close to each other. The hierarchical approach reduces computation cost by using different levels of resolution for cell problems at different macroscopic points. The method requires a hierarchy of macroscopic grid points and a corresponding nested approximation spaces with different levels of resolution. Each level of macroscopic points is assigned to an approximation finite element (FE) space which is used to solve the cell problems at the macroscopic points in that level. We prove rigorously that this hierarchical method achieves the same level of accuracy as that of the full solve where cell problems at every macroscopic point are solved using the FE spaces with the highest level of resolution, but at the essentially optimal computation cost. Numerical implementation that computes effective permeabilities of a two scale multicontinuum system via the numerical solutions of the cell problems supports the analytical results. Finally, we prove the homogenization convergence for our multiscale multi-continuum system., Comment: 26 pages, 7 figures