Adjoint method for an inverse problem of CCPF model

Abstract

The problem for determining the exchange rate function of 2D CCPF model by measurements on the partial boundary is considered and solved as one PDE-constraint optimization problem. The optimal variant is the minimum of a cost functional that quantifies the difference between the measurements and the exact solutions. Gradient-based algorithm is used to solve this optimization problem. At each step, the derivative of the cost functional with respect to the exchange rate function is calculated and only one forward solution and one adjoint solution are needed. One method based on the adjoint equation is developed and implemented. Numerical examples show the efficiency of the adjoint method.

Keywords

Adjoint method / Inverse problem / CCPF model / PDE-constraint optimization

References

[1] Adams R A. Sobolev Spaces, 1975, New York: Academic Press

[2] Water Resour. Res., 2003, 39 10

[3] Bear J, Verrruijt A. Modeling Groundwater Flow and Pollution, D. Reidel, 1987, Mass: Pub. Co., Norwell

[4] Birk S, Liedl R, Sauter M, Teutsch G. Hydraulic boundary conditions as a controllong factor in karst gensis: A numerical modeling study on artesian conduit development in gypsum. Water Resour. Res., 2003, 39(1): SBH2.1-SBH2.14

[5] Bobok E. Fluid Mechanics for Petroleum Engineers, 1993, New York: Elservier Sci

[6] Bonnans J F, Gilbert J C, Lemarechal C, Sagastizabal C A. Numerical Optimization: Theoretical and Practical Aspects, 2006 2nd Edition Berlin: Springer-Verlag

[7] Brandenburg C, Lindemann F, Ulbrich M, Ulbrich S. A continuous adjoint approach to shape optimization for Navier-Stokes flow, Optimal Control of Coupled Systems of Partial Differential Equations. International Series of Numerical Mathematics, 2009, 158: 35-56

[8] Cao Y, Gunzburger M, Hua F, Wang X. Analysis and finite element approximation of a coupled, continuum pipe-flow/Darcy model for flow in porous media with embedded conduits. Numer. Math. PDE, 2010, 27(5): 1243-1252

[9] Chen N, Gunzburger M, Hu B Calibrating the exchange coefficient in the modified coupled continuum pipe-flow model for flows in karst aquifers. J. Hydrol, 2012, 414–415: 294-301

[10] Chen W, Cheng J, Lin J, Wang L. A level set method to reconstruct the interface of discontinuity in the conductivity. Science in China Series A: Mathematics, 2009, 52(1): 29-44

[11] Feijóo G R, Malhotra M, Oberai A A Shape sensitivity calculations for exterior acoustics problems. Eng. Comput., 2001, 18(3/4): 376-391

[12] Feijóo G R, Oberai A A, Pinsky P M. An application of shape optimization in the solution of inverse acoustic scattering problems. Inverse Problems, 2004, 20(1): 199-228

[13] Hinze M, Pinnau R, Ulbrich M, Ulbrich S. Optimization with PDE Constraints, 2009, Berlin: Springer-Verlag

[14] Hua F. Modeling, analysis and simulation of Stokes-Darcy system with Beavers-Joseph interface condition, 2009, Tallahassee: Florida State University

[15] Liedl R, Sauter M, Hckinghaus D Simulation of the development of Karst aquifers using a coupled continuum pipe flow model. Water Resour. Res., 2003, 39(3): SBH6.1-SBH6.11

[16] Lu S, Chen N, Hu B, Cheng J. On the inverse problems for the coupled continuum pipe flow model for flows in karst aquifers. Inverse Problems, 2012, 28(6): 065003

[17] Narasimhan T N. Multidimensional numerical simulation of fluid flow in fractured porous media. Water Resour. Res., 1982, 18(4): 1235-1247

[18] Vogel C R. Computational methods for inverse problems, 2002, Philadelphia: SIAM

[19] Wang X. On the coupled continuum pipe flow model (CCPF) for flows in Karst aquifer. DCDS-B, 2010, 13(2): 489-501