×

The independent set perturbation method for efficient computation of sensitivities with applications to data assimilation and a finite element shallow water model. (English) Zbl 1391.76322

Summary: An adjoint model for a 2D Galerkin/Petrov-Galerkin finite element (FE) shallow water (S-W) model is developed using the Independent Set Perturbation (ISP) sensitivity analysis. Its performance in a full 4-D Var setup with a limited area shallow water equations model is assessed by comparing with the adjoint model derived by the automatic differentiation approach (TAMC), where it is used for optimising the initial conditions. It is shown that the ISP sensitivity analysis provides a very simple approach of forming the adjoint code/gradients/differentiation of discrete forward models (even complex governing equations, discretization methods and non-linear parameterizations) and is realised using a graph colouring approach combined with a perturbation method. Importantly, the adjoint is automatically updated as the forward code continues to be developed. In the test cases, it is shown that the adjoint model using the ISP sensitivity analysis can achieve the accuracy of traditional adjoint models derived by the automatic differentiation method (TAMC). Further comparison shows that the CPU time required for running the adjoint model using the ISP sensitivity analysis is much less than that required for the automatic differentiation derived adjoint model since the ISP derived adjoint CPU time scales linearly with the problem size. { } The ISP sensitivity analysis is further applied to a highly non-linear Petrov-Galerkin FE model. The perturbation size used in deriving the tangent linear model with the ISP sensitivity analysis method is then optimised and the resulting approach used to assimilate both sparse (more realistic) and dense observational data for optimising the initial conditions. A simple first order formula is developed to calculate the perturbation size for each variable, at each node and time level. By applying the ISP sensitivity method to an intermediate complexity model (a shallow water model) this paper outlines steps towards applying the approach to data assimilation (DA) problems involving realistic complex models.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
76B15 Water waves, gravity waves; dispersion and scattering, nonlinear interaction
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Cacuci, D. G.; Weber, C. F.; Oblow, E. M.; Marable, J. H., Sensitivity theory for general systems of nonlinear equations, Nucl Sci Eng, 75, 88-110, (1980)
[2] Sandu, A.; Daescu, D.; Carmichael, G. R., Direct and adjoint sensitivity analysis of chemical kinetic systems with KPP: I - theory and software tools, Atmos Environ, 37, 5083-5096, (2003)
[3] Bennett, A. F.; Chua, B. S.; Leslie, L. M., Generalized inversion of a global numerical weather prediction model, Meteor Atmos Phys, 60, 165-178, (1996)
[4] Courtier, P.; Talagrand, O., Variational assimilation of meteorological observations with the direct and adjoint shallow-water equations, Tellus, 42A, 531-549, (1990)
[5] Langland, R. H.; Elsberry, R. L.; Errico, R. M., Evaluation of physical process in an idealized extratropical cyclone using adjoint sensitivity, Quart J Roy Meteor Soc, 121, 1349-1386, (1995)
[6] Ghil M, Ide K, Bennett A, Courtier P, Kimoto M, Nagata M, et al. (Eds.), Data assimilation in meteorology and oceanography. Meteorological Society of Japan, 300 pp.
[7] Farrell, B. F., Optimal excitation of neutral Rossby waves, J Atmos Sci, 45, 163-172, (1988)
[8] Anderson, W. K.; Venkatakrishnan, V., Aerodynamic design optimization on unstructured grids with a continuous adjoint formulation, Comput Fluids, 28, 4, 443-480, (1999) · Zbl 0968.76074
[9] Barhen, J.; Cacuci, D. G.; Wagschal, J. J.; Bjerke, M. A.; Mullins, C. B., Uncertainty analysis of time-dependent nonlinear systems: theory and application to transient thermal-hydraulics, Nucl Sci Eng, 81, 23-44, (1982)
[10] Ionescu-Bujor M, Cacuci DG. Adjoint Sensitivity Analysis of the Two-fluid Model in RELAP5/MOD3.2, Invited Paper, M and C’99, International conference on mathematics and computation, reactor physics and environmental analysis in nuclear applications, vol. 1. Madrid, Spain, September 27-30; 1999. p. 745-54 [Senda Editorial, Madrid, Spain].
[11] Marta, A. C.; Mader, C. A.; Martins, J. R.R. A.; Van der Weide, E.; Alonso, J. J., A methodology for the development of discrete adjoint solvers using automatic differentiation tools, Int J Comput Fluid Dyn, 21, 9-10, 307-327, (2007) · Zbl 1184.76832
[12] Griewank, A.; Walther, A., Implementation of checkpointing for the reverse or adjoint mode of computational differentiation, ACM Trans Math Software, 26, 1, 19-45, (2000) · Zbl 1137.65330
[13] Griewank, A.; Walther, A., Evaluating derivatives: principles and techniques of algorithmic differentiation, (2008), SIAM Philadelphia · Zbl 1159.65026
[14] Griewank, A.; Juedes, D.; Utke, J., ADOL-C, a package for the automatic differentiation of algorithms written in C/C++, ACM Trans Math Software, 22, 131-167, (1996) · Zbl 0884.65015
[15] Nielsen, E. J.; Kleb, W. L., Efficient construction of discrete adjoint operators on unstructured grids using complex variables, AIAA J, 44, 4, 827-836, (2006)
[16] Giering, R.; Kaminski, T.; Slawig, T., Generating efficient derivative code with TAF: adjoint and tangent linear Euler flow around an airfoil, Future Gener Comput Syst, 21, 8, 1345-1355, (2005)
[17] Hascoet L, TAPENADE: a tool for Automatic Differentiation of programs. In: Proceedings of the ECCOMAS conference. Jyvaskyla, Finland, July 2004; 2004. p. 14.
[18] Rostaing, N.; Dalmas, S.; Galligo, A., Automatic differentiation in odyssée, Tellus A, 45, 558-568, (1993)
[19] Giering, R.; Kaminski, T., Recipes for adjoint code construction, Trans Math Software, 24, 437-474, (1998) · Zbl 0934.65027
[20] Gockenbach MS. Understanding code generated by TAMC, IAAA Paper TR00-29 2000. Department of Computational and Applied Mathematics, Rice University: Texas, USA; 2000.
[21] Vermeulen, P. T.M.; Heemink, A. W., Model-reduced variational data assimilation, Mon Weather Rev, 134, 2888-2899, (2006)
[22] Altaf, M. U.; Heemink, A. W.; Verlaan, M., Inverse shallow-water flow modeling using model reduction, Int J Multiscale Comput Eng, 7, 577-596, (2009)
[23] Cao, Y.; Zhu, J.; Navon, I. M.; Luo, Z., A reduced order approach to four-dimensional variational data assimilation using proper orthogonal decomposition, Int J Numer Meth Fluids, 53, 10, 1571-1583, (2007) · Zbl 1370.86002
[24] Daescu, D. N.; Navon, I. M., A dual-weighted approach to order reduction in 4D-VAR data assimilation, Mon Weather Review, 136, 3, 1026-1041, (2008)
[25] Fang, F.; Pain, C. C.; Navon, I. M.; Piggott, M. D.; Gorman, G. J.; Farrell, P. E., POD reduced order 4D-VAR adaptive mesh Ocean modelling, Int J Numer Meth Fluids, 60, 7, 709-732, (2009) · Zbl 1163.86002
[26] Bischof, C.; Carles, A.; Corliss, G.; Griewank, A.; Hovland, P., ADIFOR - generating derivative codes from Fortran programs, Sci Program, 1, 11-29, (1992)
[27] Carle A, Fagan M. ADIFOR 3.0 overview, Technical report CAAM-TR-00-02. Rice University; 2000.
[28] Charpentier, I., Checkpointing schemes for adjoint codes: application to the meteorological model meso-NH, SIAM J Sci Comput, 22, 6, 2135-2151, (2001) · Zbl 0989.86003
[29] (Bischof, C. H.; Bucker, H. M.; Hovland, P.; Naumann, U.; Utke, J., Lecture notes in computational science and engineering, Advances in automatic differentiation, vol. 64, (2008), Springer-Verlag), 370, doi:10.1007/978-3-540-68942-3
[30] Christakopoulos F, Jones D, Muller JD. Time-stepping for adjoint CFD codes from automatic differentiation. In: Pereira JC, Sequeira A. editors, ECCOMAS CFD 2010. Lisbon, Portugal; 2010.
[31] Zhu, K.; Navon, I. M.; Zou, X., Variational data assimilation with a variable resolution finite- element shallow-water equations model, Mon Weather Rev, 122, 5, 946-965, (1994)
[32] Navon, I. M., Finite-element simulation of the shallow-water equations model on a limited area domain, Appl Math Modell, 3, 1, 337-348, (1979) · Zbl 0438.76017
[33] Navon, I. M., FEUDX: a two-stage, high-accuracy, finite-element FORTRAN program for solving shallow-water equations, Comput Geosci, 13, 3, 255-285, (1987)
[34] Tan, W. Y., Shallow-water hydrodynamics: mathematical theory and numerical solution for a two-dimensional system of shallow-water equations, (1992), Elsevier Oceanography Series Beijing, p. 434
[35] Vreugdenhil, C. B., Numerical methods for shallow-water flow, (1994), Kluwer Academic Publishers Boston, p. 276
[36] Galewsky, J.; Scott, R. K.; Polvani, L. M., An initial-value problem for testing numerical models of the global shallow-water equations, Tellus, 56, 5, 429-440, (2004)
[37] Wang, H. H.; Halpern, P.; Douglas, J.; Dupont, T., Numerical solutions of the one-dimensional primitive equations using Galerkin approximation with localized basic functions, Mon Weather Rev, 100, 10, 738-746, (1972)
[38] Douglas, J.; Dupont, T., Galerkin method for parabolic problems, SIAM J Numer Anal, 7, 575-626, (1970) · Zbl 0224.35048
[39] Grammeltvedt, A., A survey of finite-difference schemes for the primitive equations for a barotropic fluid, Mon Weather Rev, 97, 384-404, (1969)
[40] Fang, F.; Pain, C. C.; Navon, I. M.; Gorman, G. J.; Piggott, M. D.; Allison, P. A., The independent set perturbation adjoint method: a new method of differentiating mesh-based fluids models, J Numer Meth Fluids, 66, 976-999, (2011) · Zbl 1285.76015
[41] Zienkiewicz, O. C.; Taylor, R. L.; Nithiarasu, P., The finite element method for fluid dynamics, (2005), Butterworth-Heinemann, p. 400 · Zbl 1278.76006
[42] Horwedel, J. E., GRESS, a preprocessor for sensitivity analysis of Fortran programs, (Griewank, A., Automatic differentiation of algorithms: theory, implementation, and application, (1991), SIAM Philadelphia), 243-250 · Zbl 0782.68121
[43] Liu, D. C.; Nocedal, J., On the limited memory BFGS method for large scale optimization mathematical, Programming, 45, 503-528, (1989) · Zbl 0696.90048
[44] Navon, I. M.; Zou, X.; Derber, J.; Sela, J., Variational data assimilation with an adiabatic version of the NMC spectral model, Mon Weather Rev, 120, 55-79, (1992)
[45] Le Dimet, F. X.; Talagrand, O., Variational algorithms for analysis and assimilation of meteorological observations: theoretical aspects, Tellus A, 38A, 2, 97-110, (1986)
[46] Daley, R., Atmospheric data analysis, (1991), Cambridge University Press
[47] Miller, R. N.; Zaron, E. O.; Bennett, A. F., Data assimilation in models with convective adjustment, Mon Weather Rev, 122, 2607-2613, (1994)
[48] Ide, K.; Courtier, P.; Ghil, M.; Lorenc, A. C., Unified notation for data assimilation: operational, sequential and variational, J Meteorologcial Soc Jpn, 75, 1B, 181-189, (1997)
[49] Brezillon, P.; Staub, J. F.; Perault-Staub, A. M.; Milhaud, G., Numerical estimation of the first order derivative: approximate evaluation of an optimal step, Comput Maths Appls, 7, 333-347, (1981) · Zbl 0472.65015
[50] Donea, J.; Huerta, A., Finite element methods for flow problems, (2003), John Wiley & Sons, Ltd., ISBN: 0-471-49666-9
[51] Errico, R. M., What is an adjoint model?, Bull Am Meteorol Soc, 78, 11, 2577-2591, (1997)
[52] Pires, C.; Vautard, R.; Talagrand, O., On extending the limits of variational assimilation in nonlinear chaotic systems, Tellus A, 48, 96-121, (1996)
[53] Maffray, Frédéric., On the coloration of perfect graphs, (Reed, Bruce A.; Sales, Cláudia L., Recent advances in algorithms and combinatorics, CMS books in mathematics, vol. 11, (2003), Springer-Verlag), 65-84, doi:10.1007/0-387-22444-0_3 · Zbl 1023.05052
[54] Gebremedhin, A. H.; Manne, F.; Pothen, A., What color is your Jacobian? graph coloring for computing derivatives, SIAM Rev, 47, 4, 629-705, (2005) · Zbl 1076.05034
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.