Numerical Methods For Partial Differential Algebraic Systems of Equations
Numerical Methods For Partial Differential Algebraic Systems of Equations
¾Linear solvers
¾Nonlinear solvers
¾Algorithms
¾Examples from mathematical geosciences
Approximation of PDAE’s
relative
speedup
year
D. E. Keyes, Columbia
University
Where to go past O(N) ?
¾ Hence, for instance, algebraic multigrid (AMG),
obtaining O(N) in indefinite, anisotropic, or
inhomogeneous problems
¾ Since O(N) is already optimal, there is nowhere further
“upward” to go in efficiency, but one must extend
optimality “outward”, to more general problems
AMG Framework
n
R
error easily algebraically Choose coarse grids,
damped by
smooth transfer operators, and
pointwise
error smoothers to eliminate
relaxation
these “bad” components
within a smaller
dimensional space, and
D. E. Keyes, Columbia recur
University
Computational Performance
¾ Mixed-form of RE
¾ Arithmetic mean relative permeabilities
¾ Analytical evaluation of closure relations
¾ Low-order finite differences or finite
element methods in space
¾ Backward Euler approximation in time
¾ Modified Picard iteration for nonlinear
systems
¾ Thomas algorithm for linear equation
solution
Algorithm Advancements