Trond Kvamsdal

Lars H. Odsæter, Mary F. Wheeler, Trond Kvamsdal et al.

A conservative flux postprocessing algorithm is presented for both steady-state and dynamic flow models. The postprocessed flux is shown to have the same convergence order as the original flux. An arbitrary flux approximation is projected into a conservative subspace by adding a piecewise constant correction that is minimized in a weighted $L^2$ norm. The application of a weighted norm appears to yield better results for heterogeneous media than the standard $L^2$ norm which has been considered in earlier works. We also study the effect of different flux calculations on the domain boundary. In particular we consider the continuous Galerkin finite element method for solving Darcy flow and couple it with a discontinuous Galerkin finite element method for an advective transport problem.

Lars H. Odsæter, Trond Kvamsdal, Mats G. Larson

Accurate simulation of fluid flow and transport in fractured porous media is a key challenge in subsurface reservoir engineering. Due to the high ratio between its length and width, fractures can be modeled as lower dimensional interfaces embedded in the porous rock. We apply a recently developed embedded finite element method (EFEM) for the Darcy problem. This method allows for general fracture geometry, and the fractures may cut the finite element mesh arbitrarily. We present here a velocity model for EFEM and couple the Darcy problem to a transport problem for a passive solute. The main novelties of this work is a locally conservative velocity approximation derived from the EFEM solution, and the development of a lowest order upwind finite volume method for the transport problem. This numerical model is compatible with EFEM in the sense that the same computational mesh may be applied, so that we retain the same flexibility with respect to fracture geometry and meshing. Hence, our coupled solution strategy represents a simple approach in terms of formulation, implementation and meshing. We demonstrate our model by some numerical examples on both synthetic and realistic problems, including a benchmark study for single-phase flow. Despite the simplicity of the method, the results are promising.

Florian Stadtman, Adil Rasheed, Trond Kvamsdal et al.

This article presents a comprehensive overview of the digital twin technology and its capability levels, with a specific focus on its applications in the wind energy industry. It consolidates the definitions of digital twin and its capability levels on a scale from 0-5; 0-standalone, 1-descriptive, 2-diagnostic, 3-predictive, 4-prescriptive, 5-autonomous. It then, from an industrial perspective, identifies the current state of the art and research needs in the wind energy sector. The article proposes approaches to the identified challenges from the perspective of research institutes and offers a set of recommendations for diverse stakeholders to facilitate the acceptance of the technology. The contribution of this article lies in its synthesis of the current state of knowledge and its identification of future research needs and challenges from an industry perspective, ultimately providing a roadmap for future research and development in the field of digital twin and its applications in the wind energy industry.

Sindre Stenen Blakseth, Adil Rasheed, Trond Kvamsdal et al.

Upcoming technologies like digital twins, autonomous, and artificial intelligent systems involving safety-critical applications require models which are accurate, interpretable, computationally efficient, and generalizable. Unfortunately, the two most commonly used modeling approaches, physics-based modeling (PBM) and data-driven modeling (DDM) fail to satisfy all these requirements. In the current work, we demonstrate how a hybrid approach combining the best of PBM and DDM can result in models which can outperform them both. We do so by combining partial differential equations based on first principles describing partially known physics with a black box DDM, in this case, a deep neural network model compensating for the unknown physics. First, we present a mathematical argument for why this approach should work and then apply the hybrid approach to model two dimensional heat diffusion problem with an unknown source term. The result demonstrates the method's superior performance in terms of accuracy, and generalizability. Additionally, it is shown how the DDM part can be interpreted within the hybrid framework to make the overall approach reliable.

Eivind Fonn, Harald van Brummelen, Trond Kvamsdal et al.

Reduced-basis methods (RB methods or RBMs) form one of the most promising techniques to deliver numerical solutions of parametrized PDEs in real-time performance with reasonable accuracy. For incompressible flow problems, RBMs based on LBB stable velocity-pressure spaces do not generally inherit the stability of the underlying high-fidelity model and, instead, additional stabilization techniques must be introduced. One way of bypassing the loss of LBB stability in the RBM is to inflate the velocity space with supremizer modes. This however deteriorates the performance of the RBM in the performance-critical online stage, as additional DOFs must be introduced to retain stability, while these DOFs do not effectively contribute to accuracy of the RB approximation. In this work we consider a velocity-only RB approximation, exploiting a solenoidal velocity basis. The solenoidal reduced basis emerges directly from the high-fidelity velocity solutions in the offline stage. By means of Piola transforms, the solenoidality of the velocity space is retained under geometric transformations, making the proposed RB method suitable also for the investigation of geometric parameters. To ensure exact solenoidality of the high-fidelity velocity solutions that constitute the RB, we consider approximations based on divergence-conforming compatible B-splines. We show that the velocity-only RB method leads to a significant improvement in computational efficiency in the online stage, and that the pressure solution can be recovered a posteriori at negligible extra cost. We illustrate the solenoidal RB approach by modeling steady two-dimensional Navier-Stokes flow around a NACA0015 airfoil at various angles of attack.

Sindre Stenen Blakseth, Adil Rasheed, Trond Kvamsdal et al.

In this work, we introduce, justify and demonstrate the Corrective Source Term Approach (CoSTA) -- a novel approach to Hybrid Analysis and Modeling (HAM). The objective of HAM is to combine physics-based modeling (PBM) and data-driven modeling (DDM) to create generalizable, trustworthy, accurate, computationally efficient and self-evolving models. CoSTA achieves this objective by augmenting the governing equation of a PBM model with a corrective source term generated using a deep neural network. In a series of numerical experiments on one-dimensional heat diffusion, CoSTA is found to outperform comparable DDM and PBM models in terms of accuracy -- often reducing predictive errors by several orders of magnitude -- while also generalizing better than pure DDM. Due to its flexible but solid theoretical foundation, CoSTA provides a modular framework for leveraging novel developments within both PBM and DDM. Its theoretical foundation also ensures that CoSTA can be used to model any system governed by (deterministic) partial differential equations. Moreover, CoSTA facilitates interpretation of the DNN-generated source term within the context of PBM, which results in improved explainability of the DNN. These factors make CoSTA a potential door-opener for data-driven techniques to enter high-stakes applications previously reserved for pure PBM.

Omer San, Adil Rasheed, Trond Kvamsdal

Most modeling approaches lie in either of the two categories: physics-based or data-driven. Recently, a third approach which is a combination of these deterministic and statistical models is emerging for scientific applications. To leverage these developments, our aim in this perspective paper is centered around exploring numerous principle concepts to address the challenges of (i) trustworthiness and generalizability in developing data-driven models to shed light on understanding the fundamental trade-offs in their accuracy and efficiency, and (ii) seamless integration of interface learning and multifidelity coupling approaches that transfer and represent information between different entities, particularly when different scales are governed by different physics, each operating on a different level of abstraction. Addressing these challenges could enable the revolution of digital twin technologies for scientific and engineering applications.

Suraj Pawar, Omer San, Burak Aksoylu et al.

Recent applications of machine learning, in particular deep learning, motivate the need to address the generalizability of the statistical inference approaches in physical sciences. In this letter, we introduce a modular physics guided machine learning framework to improve the accuracy of such data-driven predictive engines. The chief idea in our approach is to augment the knowledge of the simplified theories with the underlying learning process. To emphasise on their physical importance, our architecture consists of adding certain features at intermediate layers rather than in the input layer. To demonstrate our approach, we select a canonical airfoil aerodynamic problem with the enhancement of the potential flow theory. We include features obtained by a panel method that can be computed efficiently for an unseen configuration in our training procedure. By addressing the generalizability concerns, our results suggest that the proposed feature enhancement approach can be effectively used in many scientific machine learning applications, especially for the systems where we can use a theoretical, empirical, or simplified model to guide the learning module.

Yared W. Bekele, Eivind Fonn, Trond Kvamsdal et al.

Pressure oscillations at small time steps have been known to be an issue in poroelasticity simulations. A review of proposed approaches to overcome this problem is presented. Critical time steps are specified to alleviate this in finite element analyses. We present a mixed isogeometric formulation here with a view to assessing the results at very small time steps. Numerical studies are performed on Terzaghi's problem and consolidation of a layered porous medium with a very low permeability layer for varying polynomial degrees, continuities across knot spans and spatial discretizations. Comparisons are made with equal order simulations.