Strategies for solving index one DAE with non-negative constraints: Application to liquid-liquid extraction

JOURNAL OF COMPUTATIONAL PHYSICS

卷:231

Article

Metivier, Ludovic¹  Montarnal, Philippe¹ 

[1] CEA Saclay, DEN DANS SFME LSET DM2S, F-91191 Gif Sur Yvette, France

关键词: Differential algebraic system;    Non-negative constraints;    Variable stepsize and order BDF methods;    Nonlinear solver;    Conservation of mass;    Trust-regions;    Constrained dogleg methods;   

DOI  :  10.1016/j.jcp.2011.12.039

来源: Elsevier

【 摘 要 】

Liquid-liquid extraction modeling leads to solve an index one DAE system. For the sake of robustness, it is desirable to account for non-negative constraints. Based on the DASSL architecture (a classical index one DAE solver) we propose and compare three different strategies to implement these bound constraints. Each of these strategies corresponds to a different Newton modification: clipping, damping, or interior point method. The comparisons are made on two test cases: the Robertson ODE problem, and an example from liquid-liquid extraction modeling. (C) 2012 Elsevier Inc. All rights reserved.

【 授权许可】

Free