Implicit Riquier Bases for PDAE and their semi-discretizations

  • Authors:
  • Wenyuan Wu;Greg Reid;Silvana Ilie

  • Affiliations:
  • Applied Mathematics Department, University of Western Ontario, London, N6A 5B7, Canada;Applied Mathematics Department, University of Western Ontario, London, N6A 5B7, Canada;Department of Computer Science, University of Toronto, Toronto, Ontario, M5S 2E4, Canada

  • Venue:
  • Journal of Symbolic Computation
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

Complicated nonlinear systems of pde with constraints (called pdae) arise frequently in applications. Missing constraints arising by prolongation (differentiation) of the pdae need to be determined to consistently initialize and stabilize their numerical solution. In this article we review a fast prolongation method, a development of (explicit) symbolic Riquier Bases, suitable for such numerical applications. Our symbolic-numeric method to determine Riquier Bases in implicit form, without the unstable eliminations of the exact approaches, applies to square systems which are dominated by pure derivatives in one of the independent variables. The method is successful provided the prolongations with respect to a single dominant independent variable have a block structure which is uncovered by Linear Programming and certain Jacobians are nonsingular when evaluated at points on the zero sets defined by the functions of the pdae. For polynomially nonlinear pdae, homotopy continuation methods from Numerical Algebraic Geometry can be used to compute approximations of the points. Our method generalizes Pryce's method for dae to pdae. Given a dominant independent time variable, for an initial value problem for a system of pdae we show that its semi-discretization is also naturally amenable to our symbolic-numeric approach. In particular, if our method can be successfully applied to such a system of pdae, yielding an implicit Riquier Basis, then under modest conditions, the semi-discretized system of dae is also an implicit Riquier Basis.