Homotopy, polynomial equations, gross boundary data, and small Helmholtz systems
- Univ. of Cambridge (United Kingdom)
Inverse problems of the boundary measurement type appear in several geophysical contexts including DC resistivity, electromagnetic induction, and groundwater flow. The objective is to determine a spatially varying coefficient in a partial differential equation from incomplete knowledge of the dependent variable and its normal gradient at the boundary. Equivalent 2D discrete inverse problems based on the Helmholtz or modified Helmholtz equation reduce to systems of polynomial equations indicating that there are only a finite number of exact solutions, excluding certain pathological cases. A homotopy procedure decides whether real, positive solutions exist and, if so, generates the entire list. The computational complexity of the algorithm scales as M{sup M/2}, where M is the number of model parameters to be found. Measurement errors are accomodated by oversampling the boundary data at additional frequencies. For test Helmholtz and modified Helmholtz inverse problems based on (1) perfect and (2) noisy data the author generates the full list of exact solutions. The homotopy approach applies to large scale, multidimensional geophysical inverse problems but at present is practical only for small systems, up to M = 9. Recent advances in homotopy theory should, however, reduce the complexity, making larger problems tractable in the future. 18 refs., 4 figs., 2 tabs.
- OSTI ID:
- 377024
- Journal Information:
- Journal of Computational Physics, Vol. 124, Issue 2; Other Information: PBD: 15 Mar 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
Numerical Polynomial Homotopy Continuation Method and String Vacua
Parallel Ellipsoidal Perfectly Matched Layers for Acoustic Helmholtz Problems on Exterior Domains