Main navigation | Main content

HOME » PROGRAMS/ACTIVITIES » Annual Thematic Program

PROGRAMS/ACTIVITIES

Annual Thematic Program »Postdoctoral Fellowships »Hot Topics and Special »Public Lectures »New Directions »PI Programs »Math Modeling »Seminars »Be an Organizer »Annual »Hot Topics »PI Summer »PI Conference »Applying to Participate »

Talk Abstract

Seminar on Industrial Problems

Gaussian Spectral Rules for Second Order Finite-difference Schemes

Seminar on Industrial Problems

Gaussian Spectral Rules for Second Order Finite-difference Schemes

December 3, 1999

**Presented
by:
**

Schlumberger-Doll Research

druskin@ridgefield.sdr.slb.com

10:10 am

The subject of this talk is targeted grid optimization for elliptic
and time-domain problems arising in remote sensing (geophysics,
computed topography, etc.), where the solution is needed only
at few receiver points.

The optimization can be viewed as an extension of the conception of the Gaussian quadratures rules to the second order finite-difference schemes. A standard Gaussian k-point quadrature for numerical integration is chosen to be exact for 2k polynomials, and an optimal grid with k nodes is chosen to match the impedance at the receiver points for some 2k frequencies. To solve this problem we employ methods of rational approximation, linear algebra and inverse problem theory. The optimization yields exponential convergence of the impedance, i.e., the standard second order scheme with the three-point stencil exhibits spectral superconvergence.

The optimized scheme is applied to two- and three- dimensional problems in electromagnetic and acoustic well logging. Our numerical experiments exhibit exponential superconvergence at prescribed points (receivers), where the cost per grid node is close to that of the standard second order finite-difference scheme. We observe more than one order speedup for practically important problems.

Collaborators: Sergey Asvadurov (SLB), David Ingerman (Princeton-MIT), Shari Moskow (UFL) and Leonid Knizhnerman (CGE).

**Materials
Used for the Talk**

- Gaussian
spectral rules for the three-point second differences: I.
A two-point positive problem in a semiinfinite domain
pdf (387KB)
- Application
of the difference Gaussian rules to solution of hyperbolic
problems pdf (637KB)
- Three-point
finite difference schemes. Padé and the spectral Galerkin
method. I. One-sided impedance approximation pdf (557KB)
- Gaussian
spectral rules for second order finite-differ pdf (70KB)
- pcolorwave.pdf
(182KB)
- pcolorwaven.pdf (196KB)

Back to 1999-2000 Seminars on Industrial Problems