View source for PARTITIONED COUPLING FOR SOLVING STOKES-DARCY SYSTEMS

= PARTITIONED COUPLING FOR SOLVING STOKES-DARCY SYSTEMS =

PID (if applicable): [https://arxiv.org/pdf/2108.13229.pdf arxiv:2108.13229]

== Problem Statement ==

Instationary, coupled Stokes-Darcy two-domain problem: coupled systems of free flow adjacent to permeable (porous) media

=== Object of Research and Objective ===

The motivation for this is to avoid problems when using simple direct solvers (sparse direct solvers for the (linearized) subproblems): bad parallel scaling, untrustworthy solution with bad conditioning.

=== Procedure ===

Splitting the problem into two subproblems (into free flow and porous media flow) and coupling (with DuMux and preCICE): use of specialized solvers for the subdomains, coupling with preCICE.

[[File:https://i.imgur.com/BFmhwG0.jpg]]

=== Involved Disciplines ===

Environmental Systems, Mathematics

=== Data Streams ===

== Model ==

Stokes flow in the free-flow domain Darcy’s law for the porous domain

=== Discretization ===

* Time: first-order backward Euler scheme
* Space: finite volumes
** Porous Domain (Darcy): mit two-point flux approximation for pressure
** Free Flow domain (Stokes): staggered grid for pressure and velocity, upwind scheme for approximation of fluxes

=== Variables ===

{|
! Name
! Unit
! Symbol
|-
| Pressure (Dirichlet pressure)
| -
|  <math>p</math>
|-
| Velocity (Neumann velocity)
| -
|  <math>v</math>
|}

== Process Informationen ==

=== Process Steps ===

{|
! Name
! Description
! Input
! Output
! Method
! Parameter
! Environment
! Mathematical Area
|-
| Solving
| Numerical Solution of the Problem
| LGS; <math>S^{ff}</math>: <math>v^{pm}*n</math>; <math>S^{pm}</math>: <math>p^{ff}</math>
|
| Solver (UMFPACK, PD-GMRS, Bi-CGSTAB), Preconditioner (Uzawa, AMG), Coupling (preCICE, Picard iteration; inverse least-squares interface quasi-Newton)
|
| DuMux, preCICE, DUNE
| Numerical Mathematics
|}

=== Verwendete Methoden ===

=== Applied Methods ===

{|
! ID
! Name
! Process Step
! Parameter
! implemented by
|-
| wikidata:Q7001954
| Dirichlet-Neumann coupling
| Coupling
|
|
|-
| wikidata:Q1683631
| Picard iteration mit fixed-point iteration
| Coupling
|
|
|-
| wikidata:Q25098909
| inverse least-squares interface quasi-Newton
| Coupling
|
|
|-
| wikidata:Q1069090
| block-Gauss-Seidel method
| Preconditioner
|
|
|-
| wikidata:Q2467290
| Umfpack
| Solver
|
|
|-
| wikidata:Q56564057
| PD-GMRES
| Solver
| k (subiteration parameter, determined automatically), tolerance: relative residual...
|
|-
| wikidata:Q56560244
| Bi-CGSTAB
| Solver
| tolerance: relative residual...
|
|-
| wikidata:Q1471828
| AMG method
| Preconditioner
|
|
|-
| wikidata: Q17144437
| Uzawa-iterations
| Preconditioner
|
|
|-
| wikidata:Q1654069
| ILU(0) factorization
| Preconditioner
|
|
|}

=== Verwendete Software ===

=== Software used ===

{|
! ID
! Name
! Description
! Version
! Programming Language
! Dependencies
! versioned
! published
! documented
|-
| sw:8713
| preCICE
| Library for coupling simulations
| v2104.0
| Core library in C++; Bindings for Fortran, Python, C, Matlab; Adaptors dependant on simulation code (Fortran, Python, C, Matlab)
| Linux, Boost, MPI, ...
| [https://github.com/precice/precice https://github.com/precice/precice]
| [doi:10.18419/darus-2125]
| [https://precice.org/docs.html https://precice.org/docs.html]
|-
| sw:14293
| DuMux
| DUNE for Multi-{Phase, Component, Scale, Physics, …} flow and transport in porous media
|
| C++, python-bindings, utility-skripts in python
| Linux, DuNE (C++ Framework), cmake (module chains), package-config, compiler, build-essentials, dpg
| [https://git.iws.uni-stuttgart.de/dumux-repositories/dumux https://git.iws.uni-stuttgart.de/dumux-repositories/dumux]
| [https://zenodo.org/record/5152939#.YQva944zY2w https://zenodo.org/record/5152939#.YQva944zY2w]
| [https://dumux.org/docs/ https://dumux.org/docs/]
|-
| sw:18749
| ISTL
| Iterative Solver Template Library” (ISTL) which is part of the “Distributed and Unified Numerics Environment” (DUNE).
|
| C++
| Linux, DuNE (C++ Framework)
| [https://gitlab.dune-project.org/core/dune-istl https://gitlab.dune-project.org/core/dune-istl]
| [doi:10.1007/978-3-540-75755-9_82]
| [https://www.dune-project.org/modules/dune-istl/ https://www.dune-project.org/modules/dune-istl/]
|}

=== Hardware ===

{|
! ID
! Name
! Processor
! Compiler
! #Nodes
! #Cores
|-
|
|
| AMD EPYC 7551P CPU
|
| 1
| 1
|}

=== Input Data ===

{|
! ID
! Name
! Size
! Data Structure
! Format Representation
! Format Exchange
! binary/text
! proprietary
! to publish
! to archive
|-
|
| LGS
| O($2*10^6$) Matrix size
| Data structure in DUNE/DuMux
|
|
| numbers
| open
| ?
| ?
|}

=== Output Data ===

{|
! ID
! Name
! Size
! Data Structure
! Format Representation
! Format Exchange
! binary/text
! proprietary
! to publish
! to archive
|-
|
| Solution Vector
|
| Data Structure in DUNE/DuMux
|
|
| numbers
| open
| ?
| ?
|-
|
| Runtime Behavior
|
| Data Structure in DUNE/DuMux
|
|
| numbers
| open
| ?
| ?
|-
| wikidata:Q18812775
| vtk-image files, from the daten structure
| ?
|
| ASCII or binary
| open
| ?
| ?
|
|
|}

== Reproducibility ==

=== Mathematical Reproducibility ===

yes, by all parameters

=== Runtime Reproducibility ===

yes, for same input samples

=== Reproducibility of Results ===

Due to floating point arithmetic, no bitwise reproducibility

=== Reproducibility on original Hardware ===

=== Reproducibility on other Hardware ===

a) Serial Computation

b) Parallel Computation

=== Transferability to ===

a) similar model parameters (other initial and boundary values)

b) other models

= Legend =

The following abbreviations are used in the document to indicate/resolve IDs:

doi: DOI / [https://dx.doi.org/ https://dx.doi.org/]

sw: swMATH / [https://swmath.org/software/ https://swmath.org/software/]

wikidata: [https://www.wikidata.org/wiki/ https://www.wikidata.org/wiki/]
[[Category:Workflow]]