Abstract
We present a novel partitioned coupling algorithm to solve first-order time-dependent non-linear problems (e.g. transient heat conduction). The spatial domain is partitioned into a set of totally disconnected subdomains. The continuity conditions at the interface are modeled using a dual Schur formulation where the Lagrange multipliers represent the interface fluxes (or the reaction forces) that are required to maintain the continuity conditions. The interface equations along with the subdomain equations lead to a system of differential algebraic equations (DAEs). For the resulting equations a numerical algorithm is developed, which includes choosing appropriate constraint stabilization techniques. The algorithm first solves for the interface Lagrange multipliers, which are subsequently used to advance the solution in the subdomains. The proposed coupling algorithm enables arbitrary numeric schemes to be coupled with different time steps (i.e. it allows subcycling) in each subdomain. This implies that existing software and numerical techniques can be used to solve each subdomain separately. The coupling algorithm can also be applied to multiple subdomains and is suitable for parallel computers. We present examples showing the feasibility of the proposed coupling algorithm.
Original language | English (US) |
---|---|
Pages (from-to) | 1385-1415 |
Number of pages | 31 |
Journal | International Journal for Numerical Methods in Engineering |
Volume | 75 |
Issue number | 12 |
DOIs | |
State | Published - Sep 17 2008 |
Keywords
- Constraint stabilization
- Differential algebraic equation (DAE)
- Domain decomposition methods
- FETI method
- Mixed- and multi-time-step methods
- Non-linear transient analysis
- Subcycling
ASJC Scopus subject areas
- Numerical Analysis
- General Engineering
- Applied Mathematics