A DEIM driven reduced basis method for the diffuse Stokes/Darcy model coupled at parametric phase-field interfaces

Computational Geosciences - Tập 26 - Trang 1465-1502 - 2022
Stein K.F. Stoter1, Etienne Jessen1, Viktor Niedens1, Dominik Schillinger1
1Institute of Mechanics and Computational Mechanics, Leibniz Universität Hannover, Hannover, Germany

Tóm tắt

In this article, we develop a reduced basis method for efficiently solving the coupled Stokes/Darcy equations with parametric internal geometry. To accommodate possible changes in topology, we define the Stokes and Darcy domains implicitly via a phase-field indicator function. In our reduced order model, we approximate the parameter-dependent phase-field function with a discrete empirical interpolation method (DEIM) that enables affine decomposition of the associated linear and bilinear forms. In addition, we introduce a modification of DEIM that leads to non-negativity preserving approximations, thus guaranteeing positive-semidefiniteness of the system matrix. We also present a strategy for determining the required number of DEIM modes for a given number of reduced basis functions. We couple reduced basis functions on neighboring patches to enable the efficient simulation of large-scale problems that consist of repetitive subdomains. We apply our reduced basis framework to efficiently solve the inverse problem of characterizing the subsurface damage state of a complete in-situ leach mining site.

Tài liệu tham khảo

Quarteroni, A., Manzoni, A., Negri, F.: Reduced basis methods for partial differential equations: an introduction (2015)

Hesthaven, S.J., Rozza, G., Stamm, B.: Certified reduced basis methods for parametrized partial differential equations (2015)

Maury, B.: A fat boundary method for the Poisson problem in a domain with holes. J. Sci. Comput. 16(3), 319–339 (2001)

Li, X., Lowengrub, J., Rȧtz, A., Voigt, A.: Solving PDEs in complex geometries: a diffuse domain approach. Commun. Math. Sci. 7(1), 81 (2009)

Oden, J.T., Hawkins, A., Prudhomme, S.: General diffuse-interface theories and an approach to predictive tumor growth modeling. Math. Models Methods Appl. Sci. 20(3), 477–517 (2010)

Nguyen, L.H., Stoter, S.K.F., Baum, T.J., Kirschke, S., Ruess, M., Yosibash, Z., Schillinger, D.: Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures. Int. J. Numer. Methods Biomed. Eng. https://doi.org/10.1002/cnm.2880 (2017)

Nguyen, L.H., Stoter, S.K.F., Ruess, M., Sanchez Uribe, M.A., Schillinger, D.: The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries. Int. J. Numer. Methods Eng. 113(4), 601–633 (2018). https://doi.org/10.1002/nme.5628

Boyer, F.: Mathematical study of multi-phase flow under shear through order parameter formulation. Asymptot. Anal. 20(2), 175–212 (1999)

Alnæs, M.S., Blechta, J., Hake, J., Johansson, A., Kehlet, B., Logg, A., Richardson, C., Ring, J., Rognes, M.E., Wells, G.N.: The fenics project version 1.5. Arch. Numer. Softw. 3(100), 9–23 (2015)

Balay, S, Abhyankar, S, Adams, M.F, Brown, J, Brune, P, Buschelman, K, Dalcin, L, Eijkhout, V, Gropp, W.D, Kaushik, D, Knepley, M.G, McInnes, L.C, Rupp, K, Smith, B.F, Zampini, S, Zhang, H.: PETSc Web page. Available: https://petsc.org/ (2021)

In Situ Leach Mining (ISL) of Uranium - World Nuclear Association: https://world-nuclear.org/information-library/nuclear-fuel-cycle/mining-of-uranium/in-situ-leach-mining-of-uranium.aspx, Accessed 01 June 2021

Honeymoon project: https://www.epa.sa.gov.au/files/4771306_honeymoon_application.pdf, Accessed 01 June 2021