Parallel 3-D pseudospectral simulation of seismic wave propagation

Takashi Furumura, B. L N Kennett, Hiroshi Takenaka

Research output: Contribution to journalArticle

68 Citations (Scopus)

Abstract

Three-dimensional pseudospectral modeling for a realistic scale problem is still computationally very intensive, even when using current powerful computers. To overcome this, we have developed a parallel pseudospectral code for calculating the 3-D wavefield by concurrent use of a number of processors. The parallel algorithm is based on a partition of the computational domain, where the field quantities are distributed over a number of processors and the calculation is concurrently done in each subdomain with interprocessor communications. Experimental performance tests using three different styles of parallel computers achieved a fairly good speed up compared with conventional computation on a single processor: maximum speed-up rate of 26 using 32 processors of a Thinking Machine CM-5 parallel computer, 1.6 using a Digital Equipment DEC-Alpha two-CPU workstation, and 4.6 using a cluster of eight Sun Microsystems SPARC-Station 10 (SPARC-10) workstations connected by an Ethernet. The result of this test agrees well with the performance theoretically predicted for each system. To demonstrate the feasibility of our parallel algorithm, we show three examples: 3-D acoustic and elastic modeling of fault-zone trapped waves and the calculation of elastic wave propagation in a 3-D syncline model.

Original languageEnglish
Pages (from-to)279-288
Number of pages10
JournalGeophysics
Volume63
Issue number1
Publication statusPublished - Jan 1998
Externally publishedYes

Fingerprint

Seismic waves
seismic waves
seismic wave
Wave propagation
wave propagation
central processing units
parallel computers
workstations
Parallel algorithms
simulation
trapped wave
interprocessor communication
three-dimensional modeling
elastic wave
syncline
synclines
Computer workstations
Ethernet
fault zone
performance tests

ASJC Scopus subject areas

  • Geochemistry and Petrology
  • Geophysics

Cite this

Parallel 3-D pseudospectral simulation of seismic wave propagation. / Furumura, Takashi; Kennett, B. L N; Takenaka, Hiroshi.

In: Geophysics, Vol. 63, No. 1, 01.1998, p. 279-288.

Research output: Contribution to journalArticle

Furumura, T, Kennett, BLN & Takenaka, H 1998, 'Parallel 3-D pseudospectral simulation of seismic wave propagation', Geophysics, vol. 63, no. 1, pp. 279-288.
Furumura, Takashi ; Kennett, B. L N ; Takenaka, Hiroshi. / Parallel 3-D pseudospectral simulation of seismic wave propagation. In: Geophysics. 1998 ; Vol. 63, No. 1. pp. 279-288.
@article{6d76a64906b14971b1dfe37859479b8b,
title = "Parallel 3-D pseudospectral simulation of seismic wave propagation",
abstract = "Three-dimensional pseudospectral modeling for a realistic scale problem is still computationally very intensive, even when using current powerful computers. To overcome this, we have developed a parallel pseudospectral code for calculating the 3-D wavefield by concurrent use of a number of processors. The parallel algorithm is based on a partition of the computational domain, where the field quantities are distributed over a number of processors and the calculation is concurrently done in each subdomain with interprocessor communications. Experimental performance tests using three different styles of parallel computers achieved a fairly good speed up compared with conventional computation on a single processor: maximum speed-up rate of 26 using 32 processors of a Thinking Machine CM-5 parallel computer, 1.6 using a Digital Equipment DEC-Alpha two-CPU workstation, and 4.6 using a cluster of eight Sun Microsystems SPARC-Station 10 (SPARC-10) workstations connected by an Ethernet. The result of this test agrees well with the performance theoretically predicted for each system. To demonstrate the feasibility of our parallel algorithm, we show three examples: 3-D acoustic and elastic modeling of fault-zone trapped waves and the calculation of elastic wave propagation in a 3-D syncline model.",
author = "Takashi Furumura and Kennett, {B. L N} and Hiroshi Takenaka",
year = "1998",
month = "1",
language = "English",
volume = "63",
pages = "279--288",
journal = "Geophysics",
issn = "0016-8033",
publisher = "Society of Exploration Geophysicists",
number = "1",

}

TY - JOUR

T1 - Parallel 3-D pseudospectral simulation of seismic wave propagation

AU - Furumura, Takashi

AU - Kennett, B. L N

AU - Takenaka, Hiroshi

PY - 1998/1

Y1 - 1998/1

N2 - Three-dimensional pseudospectral modeling for a realistic scale problem is still computationally very intensive, even when using current powerful computers. To overcome this, we have developed a parallel pseudospectral code for calculating the 3-D wavefield by concurrent use of a number of processors. The parallel algorithm is based on a partition of the computational domain, where the field quantities are distributed over a number of processors and the calculation is concurrently done in each subdomain with interprocessor communications. Experimental performance tests using three different styles of parallel computers achieved a fairly good speed up compared with conventional computation on a single processor: maximum speed-up rate of 26 using 32 processors of a Thinking Machine CM-5 parallel computer, 1.6 using a Digital Equipment DEC-Alpha two-CPU workstation, and 4.6 using a cluster of eight Sun Microsystems SPARC-Station 10 (SPARC-10) workstations connected by an Ethernet. The result of this test agrees well with the performance theoretically predicted for each system. To demonstrate the feasibility of our parallel algorithm, we show three examples: 3-D acoustic and elastic modeling of fault-zone trapped waves and the calculation of elastic wave propagation in a 3-D syncline model.

AB - Three-dimensional pseudospectral modeling for a realistic scale problem is still computationally very intensive, even when using current powerful computers. To overcome this, we have developed a parallel pseudospectral code for calculating the 3-D wavefield by concurrent use of a number of processors. The parallel algorithm is based on a partition of the computational domain, where the field quantities are distributed over a number of processors and the calculation is concurrently done in each subdomain with interprocessor communications. Experimental performance tests using three different styles of parallel computers achieved a fairly good speed up compared with conventional computation on a single processor: maximum speed-up rate of 26 using 32 processors of a Thinking Machine CM-5 parallel computer, 1.6 using a Digital Equipment DEC-Alpha two-CPU workstation, and 4.6 using a cluster of eight Sun Microsystems SPARC-Station 10 (SPARC-10) workstations connected by an Ethernet. The result of this test agrees well with the performance theoretically predicted for each system. To demonstrate the feasibility of our parallel algorithm, we show three examples: 3-D acoustic and elastic modeling of fault-zone trapped waves and the calculation of elastic wave propagation in a 3-D syncline model.

UR - http://www.scopus.com/inward/record.url?scp=0031699665&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0031699665&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0031699665

VL - 63

SP - 279

EP - 288

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 1

ER -