CoChat
Sign Up
PreprintCited by 58Advanced Numerical Methods in Computational MathematicsOptimization and Variational AnalysisAdvanced Control Systems Optimization

A framework for automated PDE-constrained optimisation

Open access
Funke, S. W., Farrell, P. E.·2013·arXiv
DOI 10.48550/arxiv.1302.3894arXiv 1302.3894
Open this paper in CoChat — free accountOpen PDF

Abstract

A generic framework for the solution of PDE-constrained optimisation problems based on the FEniCS system is presented. Its main features are an intuitive mathematical interface, a high degree of automation, and an efficient implementation of the generated adjoint model. The framework is based upon the extension of a domain-specific language for variational problems to cleanly express complex optimisation problems in a compact, high-level syntax. For example, optimisation problems constrained by the time-dependent Navier-Stokes equations can be written in tens of lines of code. Based on this high-level representation, the framework derives the associated adjoint equations in the same domain-specific language, and uses the FEniCS code generation technology to emit parallel optimised low-level C++ code for the solution of the forward and adjoint systems. The functional and gradient information so computed is then passed to the optimisation algorithm to update the parameter values. This approach works both for steady-state as well as transient, and for linear as well as nonlinear governing PDEs and a wide range of functionals and control parameters. We demonstrate the applicability and efficiency of this approach on classical textbook optimisation problems and advanced examples.

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

DataCite · 10.48550/arxiv.1302.3894

Cites

Statistics

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • Previous computational analyses of cell migration have identified diffusivity values from approximately 50 to 3000 μm2/hr.

    “Previous computational analyses of cell migration have identified a wide range of diffusivity values, from approximately 50 to 3000μm2/hr9,40,41.”

    From Inference of weak-form partial differential equations describing migration and proliferation mechanisms in wound healing experiments on cancer cells (2025) · cites this paper as [41]

    The Introduction (p. 2) discusses PDE-constrained optimisation, and the wetting-and-drying example (pp. 21–23) concerns tsunami modelling. Searches for cell migration, diffusivity, and the stated range found no passage reporting cell-migration diffusivity values. (LLM, conf=0.76)

    Experimental — accuracy not yet measured

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

The citation lists haven’t been fetched yet

References 55 · Cited by 58

Research works printed but not linked: 40

  1. ALNÆS, M. S. 2011. UFL: A finite element form language. In Automated Solution of Differential Equations by the Finite Element Method , A. Logg, K. A. Mardal, and G. N. Wells, Eds. Springer-Verlag, Berlin, Heidelberg, New York, Chapter 17, 299-334.
  2. ALNÆS, M. S., LOGG, A., ØLGAARD, K. B., ROGNES, M. E., AND WELLS, G. N. 2012. Unified Form Language: A domain-specific language for weak formulations of partial differential equations. arXiv:1211.4047 [cs.MS]. · arXiv
  3. ANDERSSON, J., ˚ AKESSON, J., AND DIEHL, M. 2012. CasADi - a symbolic package for automatic differentiation and optimal control. In Recent Advances in Algorithmic Differentiation , S. Forth, P. Hovland, E. Phipps, J. Utke, and A. Walther, Eds. Lecture Notes in Computational Science and Engineering Series, vol. 87. Springer-Verlag, Berlin, Heidelberg, New York.
  4. BEN ELGHALI, S. E., BALME, R., LE SAUX, K., EL HACHEMI BENBOUZID, M., CHARPENTIER, J. F., AND HAUVILLE, F. 2007. A simulation model for the evaluation of the electrical power potential harnessed by a marine current turbine. IEEE Journal of Oceanic Engineering 32, 4, 786-797.
  5. BIEGLER, L. T., NOCEDAL, J., AND SCHMID, C. 1995. A reduced Hessian method for large-scale constrained optimization. SIAM Journal on Optimization 5, 2, 314-347.
  6. BIROS, G. AND GHATTAS, O. 2000. Parallel Lagrange-Newton-Krylov-Schur methods for PDE-constrained optimization. Part I: The Krylov-Schur solver. SIAM Journal on Scientific Computing 27 , 2005.
  7. CORLISS, G. F. AND GRIEWANK, A. 1993. Operator overloading as an enabling technology for automatic differentiation. In Proceedings of the Annual Object-Oriented Numerics Conference . Sunriver, Oregon.
  8. DIVETT, T., VENNELL, R., AND STEVENS, C. 2011. Optimisation of multiple turbine arrays in a channel with tidally reversing flow by numerical modelling with adaptive mesh. In Proceedings of the European Wave and Tidal Energy Conference 2011 . Southampton, UK.
  9. FARRELL, P. E., HAM, D. A., FUNKE, S. W., AND ROGNES, M. E. 2012. Automated derivation of the adjoint of high-level transient finite element programs. SIAM Journal on Scientific Computing . Accepted.
  10. GRIEWANK, A. AND WALTHER, A. 2000. Algorithm 799: revolve: An implementation of checkpointing for the reverse or adjoint mode of computational differentiation. ACM Transactions on Mathematical Software 26, 1, 19-45.
  11. HINZE, M., PINNAU, R., ULBRICH, M., AND ULBRICH, S. 2009. Optimization with PDE constraints . Mathematical Modelling: Theory and Applications Series, vol. 23. Springer-Verlag, Berlin, Heidelberg, New York.
  12. JAMESON, A. 1988. Aerodynamic design via control theory. Journal of Scientific Computing 3, 3, 233-260. JONES, E., OLIPHANT, T., PETERSON, P., ET AL. 2001. SciPy: Open source scientific tools for Python.
  13. KAIPIO, J. AND SOMERSALO, E. 2004. Statistical and computational inverse problems . Applied Mathematical Sciences Series, vol. 160. Springer-Verlag, Berlin, Heidelberg, New York.
  14. KARUSH, W. 1939. Minima of functions of several variables with inequalities as side constraints. M.S. thesis, Department of Mathematics, University of Chicago.
  15. KIRBY, R. C. AND LOGG, A. 2006. A compiler for variational forms. ACM Transactions on Mathematical Software 32, 3, 417-444.
  16. KIRBY, R. C. AND LOGG, A. 2007. Efficient compilation of a class of variational forms. ACM Transactions on Mathematical Software 33, 3, 93-138.
  17. KUHN, H. W. AND TUCKER, A. W. 1951. Nonlinear programming. In Second Berkeley symposium on mathematical statistics and probability . Berkeley, California, 481-492.
  18. LIONS, J. L. 1971. Optimal control of systems governed by partial differential equations . Grundlehren der mathematischen Wissenschaften Series, vol. 170. Springer-Verlag, Berlin, Heidelberg, New York.
  19. LOGG, A., MARDAL, K.-A., AND WELLS, G. N., Eds. 2012. Automated Solution of Differential Equations by the Finite Element Method . Lecture Notes in Computational Science and Engineering Series, vol. 84. Springer-Verlag, Berlin, Heidelberg, New-York.
  20. LOGG, A. AND WELLS, G. N. 2010. DOLFIN: Automated finite element computing. ACM Transactions on Mathematical Software 37, 2.
  21. LOGG, A., WELLS, G. N., AND HAKE, J. E. 2011. DOLFIN: A C++/Python finite element library. In Automated Solution of Differential Equations by the Finite Element Method , A. Logg, K.-A. Mardal, and G. N. Wells, Eds. Springer-Verlag, Berlin, Heidelberg, New-York, Chapter 10.
  22. MARKALL, G. R., SLEMMER, A., HAM, D. A., KELLY, P. H. J., CANTWELL, C. D., AND SHERWIN, S. J. 2012. Finite element assembly strategies on multi-core and many-core architectures. International Journal for Numerical Methods in Fluids 71, 1, 80-97.
  23. MATSUYAMA, M. AND TANAKA, H. 2001. An experimental study of the highest runup height in the 1993 Hokkaido Nansei-oki earthquake tsunami. In National Tsunami Hazard Mitigation Program Review and International Tsunami Symposium (ITS 2001) . 879-889.
  24. MEDEIROS, S. C. AND HAGEN, S. C. 2013. Review of wetting and drying algorithms for numerical tidal flow models. International Journal for Numerical Methods in Fluids 71, 4, 473-487.
  25. NASH, S. 1984. Newton-type minimization via the Lanczos method. SIAM Journal on Numerical Analysis 21, 4, 770-788.
  26. NAUMANN, U. 2012. The Art of Differentiating Computer Programs: An Introduction to Algorithmic Differentiation . Software, Environments, and Tools Series, vol. 24. SIAM, Philadelphia.
  27. NELDER, J. A. AND MEAD, R. 1965. A simplex method for function minimization. The Computer Journal 7, 4, 308-313.
  28. NOCEDAL, J. AND WRIGHT, S. J. 2006. Numerical optimization 2 Ed. Springer Series in Operations Research and Financial Engineering. Springer-Verlag, Berlin, Heidelberg, New-York.
  29. ØLGAARD, K. B. AND WELLS, G. N. 2010. Optimizations for quadrature representations of finite element tensors through automated code generation. ACM Transactions on Mathematical Software 37, 1, 8-1.
  30. OROZCO, C. E. AND GHATTAS, O. N. 1997. A reduced SAND method for optimal design of non-linear structures. International Journal for Numerical Methods in Engineering 40, 15, 2759-2774.
  31. RABIER, F., J ¨ ARVINEN, H., KLINKER, E., MAHFOUF, J. F., AND SIMMONS, A. 2000. The ECMWF operational implementation of four-dimensional variational assimilation. I: Experimental results with simplified physics. Quarterly Journal of the Royal Meteorological Society 126, 564, 1143-1170.
  32. SCH¨ OBERL, J. AND ZULEHNER, W. 2007. Symmetric indefinite preconditioners for saddle point problems with applications to PDE-constrained optimization problems. SIAM Journal on Matrix Analysis and Applications 29, 3, 752-773.
  33. SONG, L., ZHOU, J., LI, Q., YANG, X., AND ZHANG, Y. 2011. An unstructured finite volume model for dambreak floods with wet/dry fronts over complex topography. International Journal for Numerical Methods in Fluids 67, 8, 960-980.
  34. STUMM, P. AND WALTHER, A. 2009. Multistage approaches for optimal offline checkpointing. SIAM Journal on Scientific Computing 31, 3, 1946-1967.
  35. TR´ EMOLI ` ERES, R., LIONS, J.-L., AND GLOWINSKI, R. 1981. Numerical analysis of variational inequalities . Studies in Mathematics and its Applications Series, vol. 8. North-Holland, Amsterdam, New-York, Oxford.
  36. WELLS, G. N., KUHL, E., AND GARIKIPATI, K. 2006. A discontinuous Galerkin method for the CahnHilliard equation. Journal of Computational Physics 218, 2, 860-877.
  37. WESTERINK, J. J., LUETTICH, R. A., FEYEN, J. C., ATKINSON, J. H., DAWSON, C., ROBERTS, H. J., POWELL, M. D., DUNION, J. P., KUBATKO, E. J., AND POURTAHERI, H. 2008. A basin- to channel-scale unstructured grid hurricane storm surge model applied to southern Louisiana. Monthly Weather Review 136, 3, 833-864.
  38. YALCINER, A. C., IMAMURA, F., AND SYNOLAKIS, C. E. 2011. Amplitude evolution and runup of long waves: Comparison of experimental and numerical data on a 3D complex topography. In Advanced Numerical Models For Simulating Tsunami Waves And Runup , P. L. Liu, H. Yeh, and C. Synolakis, Eds. World Scientific Publishing Company, Chapter 9, 243-247.
  39. ZHANG, J., LIN, J., SUN, X., AND CHEN, C. 2009. Application of an improved wetting and drying scheme in POM. In Proceedings of the 2009 International Conference on Energy and Environment Technology . Vol. 2. IEEE Computer Society, Washington, DC, 427-432.
  40. ZHU, C., BYRD, R. H., LU, P., AND NOCEDAL, J. 1997. Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization. ACM Transactions on Mathematical Software 23, 4, 550560.

Authors
Funke, S. W., & Farrell, P. E.
ORCID
Venue
arXiv
Publisher
arXiv
Published
2013
ISSN
2331-8422 (ISSN-L)
DOI
10.48550/arxiv.1302.3894
arXiv
1302.3894
OpenAlex
W1530180322
Cited by
58
References
55
Citations by year
2025: 1 · 2024: 5 · 2023: 7 · 2022: 6 · 2021: 5 · 2020: 9 · 2019: 2 · 2018: 6 · 2017: 7 · 2016: 4 · 2015: 1 · 2014: 2 · 2013: 3
Type
posted content
Language
en
Open access
Licence
arxiv.org/licenses/nonexclusive-distrib/1.0/
Locations
  • arXiv (Cornell University) (submittedVersion)
  • arXiv (Cornell University)
Funders
Engineering and Physical Sciences Research Council (EP/I00405X/1)
Topics
  • Physical Sciences › Engineering › Computational Mechanics › Advanced Numerical Methods in Computational Mathematics
  • Physical Sciences › Computer Science › Computational Theory and Mathematics › Optimization and Variational Analysis
  • Physical Sciences › Engineering › Control and Systems Engineering › Advanced Control Systems Optimization
Keywords
adjoint model, FEniCS, domain-specific languages, variational problems, code generation, Navier-Stokes equations, automatic differentiation, gradient computation, parameter optimization

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review