Simulation of nonisothermal, multiphase flow through fractured geothermal reservoirs involves the solution of a system of strongly nonlinear algebraic equations. The Newton-Raphson method used to solve such a nonlinear system of equations requires the evaluation of a Jacobian matrix. In this paper we discuss automatic differentiation (AD) as a method for analytically computing the Jacobian matrix of derivatives. Robustness and efficiency of the AD-generated derivative codes are compared with a conventional derivative computation approach based on first-order finite differences.

Press the Back button in your browser, or search again.

Copyright 2002, Stanford Geothermal Program: Readers who download papers from this site should honor the copyright of the original authors and may not copy or distribute the work further without the permission of the original publisher.