TY - JOUR
T1 - On dual Schur domain decomposition method for linear first-order transient problems
AU - Nakshatrala, K. B.
AU - Prakash, A.
AU - Hjelmstad, Keith
N1 - Funding Information:
The first author (KBN) would like to thank Professor Albert Valocchi and Professor Daniel Tortorelli for their valuable discussions and encouragement. The research reported herein was partly supported (KBN) by the Department of Energy through a SciDAC-2 Project (Grant No. DOE DE-FC02-07ER64323). This support is gratefully acknowledged. The opinions expressed in this paper are those of the authors and do not necessarily reflect that of the sponsor.
PY - 2009/11/20
Y1 - 2009/11/20
N2 - This paper addresses some numerical and theoretical aspects of dual Schur domain decomposition methods for linear first-order transient partial differential equations. The spatially discrete system of equations resulting from a dual Schur domain decomposition method can be expressed as a system of differential-algebraic equations (DAEs). In this work, we consider the trapezoidal family of schemes for integrating the ordinary differential equations (ODEs) for each subdomain and present four different coupling methods, corresponding to different algebraic constraints, for enforcing kinematic continuity on the interface between the subdomains. Unlike the continuous formulation, the discretized formulation of the transient problem is unable to enforce simultaneously the continuity of both the primary variable and its rate along the subdomain interface (except for the backward Euler method). Method 1 (d-continuity) is based on the conventional approach using continuity of the primary variable and we show that this method is unstable for a lot of commonly used time integrators including the mid-point rule. To alleviate this difficulty, we propose a new Method 2 (modified d-continuity) and prove its stability for coupling all time integrators in the trapezoidal family (except the forward Euler). Method 3 (v-continuity) is based on enforcing the continuity of the time derivative of the primary variable. However, this constraint introduces a drift in the primary variable on the interface. We present Method 4 (Baumgarte stabilized) which uses Baumgarte stabilization to limit this drift and we derive bounds for the stabilization parameter to ensure stability. Our stability analysis is based on the "energy" method, and one of the main contributions of this paper is the extension of the energy method (which was previously introduced in the context of numerical methods for ODEs) to assess the stability of numerical formulations for index-2 differential-algebraic equations (DAEs). Finally, we present numerical examples to corroborate our theoretical predictions.
AB - This paper addresses some numerical and theoretical aspects of dual Schur domain decomposition methods for linear first-order transient partial differential equations. The spatially discrete system of equations resulting from a dual Schur domain decomposition method can be expressed as a system of differential-algebraic equations (DAEs). In this work, we consider the trapezoidal family of schemes for integrating the ordinary differential equations (ODEs) for each subdomain and present four different coupling methods, corresponding to different algebraic constraints, for enforcing kinematic continuity on the interface between the subdomains. Unlike the continuous formulation, the discretized formulation of the transient problem is unable to enforce simultaneously the continuity of both the primary variable and its rate along the subdomain interface (except for the backward Euler method). Method 1 (d-continuity) is based on the conventional approach using continuity of the primary variable and we show that this method is unstable for a lot of commonly used time integrators including the mid-point rule. To alleviate this difficulty, we propose a new Method 2 (modified d-continuity) and prove its stability for coupling all time integrators in the trapezoidal family (except the forward Euler). Method 3 (v-continuity) is based on enforcing the continuity of the time derivative of the primary variable. However, this constraint introduces a drift in the primary variable on the interface. We present Method 4 (Baumgarte stabilized) which uses Baumgarte stabilization to limit this drift and we derive bounds for the stabilization parameter to ensure stability. Our stability analysis is based on the "energy" method, and one of the main contributions of this paper is the extension of the energy method (which was previously introduced in the context of numerical methods for ODEs) to assess the stability of numerical formulations for index-2 differential-algebraic equations (DAEs). Finally, we present numerical examples to corroborate our theoretical predictions.
KW - Coupling algorithms
KW - Differential-algebraic equations
KW - Dual Schur domain decomposition method
KW - Generalized trapezoidal family
UR - http://www.scopus.com/inward/record.url?scp=69949092617&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=69949092617&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2009.07.016
DO - 10.1016/j.jcp.2009.07.016
M3 - Article
AN - SCOPUS:69949092617
SN - 0021-9991
VL - 228
SP - 7957
EP - 7985
JO - Journal of Computational Physics
JF - Journal of Computational Physics
IS - 21
ER -