Sensitivity analysis of a nonlinear Newton-Krylov solver for heat transfer with phase change.
- Dana A.
- Douglas B.
- Bryan R.
Development of a complex metal-casting computer model requires information about how varying the problem parameters affects the results (metal flow and solidification). For example, we would like to know how the last point to solidify or the cooling rate at a given location changes when the physical properties of the metal, boundary conditions, or mold geometry are changed. As a preliminary step towards a complete sensitivity analysis of a three-dimensional casting simulation, we examine a one-dimensional version of a metal-alloy phase-change conductive-heat-transfer model by means of Automatic Differentiation (AD). This non-linear 'Jacobian-free' method is a combination of an outer Newton-based iteration and an inner conjugate gradient-like (Krylov) iteration. The implicit solution algorithm has enthalpy as the dependent variable from which temperatures are determined. We examine the sensitivities of the difference between an exact analytical solution for the final temperature and that produced by this algorithm to the problem parameters. In all there are 17 parameters (12 physical constants such as liquid density, heat capacity, and thermal conductivity, 2 initial and boundary condition parameters, the final solution time, and 2 algorithm tolerances). We apply AD in the forward and reverse mode and verify the sensitivities by means of finite differences. In general, the finite-difference method requires at least N+1 computer runs to determine sensitivities for N problem parameters. By forward and reverse, we mean the direction through the solution and in time and space in which the derivative values are obtained. The forward mode is typically more efficient for determining the sensitivity of many responses to one or a few parameters, while the reverse mode is better suited for sensitivities of one or a few responses with respect to many parameters. The sensitivities produced by all the methods agreed to at least three significant figures. The forward and reverse AD code run times were similar and were approximately 34% faster than those of the finite-difference sensitivities. Real problems in three dimensions will certainly have many more parameters describing mold geometry and pouring conditions. If the trend seen here holds true reverse mode AD is favored since the computational time increases only slightly for additional parameters.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Organization:
- USDOE
- OSTI ID:
- 976161
- Report Number(s):
- LA-UR-02-2478; TRN: US201009%%577
- Resource Relation:
- Conference: Submitted to: WCCM V - Fifth World Congress on Computational Mechnanics, Vienna, Austria, July 7-12, 2002
- Country of Publication:
- United States
- Language:
- English
Similar Records
NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations
AUTOMATIC DIFFERENTIATION OF AN EULERIAN HYDROCODE