Parallel implementation of a velocity-stress staggered-grid finite-difference method for 2-D poroelastic wave propagation

Date

2006-10-01

Author

Sheen, Dong-Hoon
Tuncay, Kağan
Baag, Chang-Eob
Ortoleva, Peter J.

Metadata

Show full item record

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Item Usage Stats

96
views

0
downloads

Numerical simulation of wave propagation in poroelastic media demands significantly more computational capability compared to elastic media simulation. Use of serial codes in a single scientific workstation limits the size of problem. To overcome this difficulty, a parallel velocity-stress staggered-grid finite-difference method is developed for efficient simulation of wave propagation in 2-D poroelastic media. The finite difference formulation of Biot's theory has the properties of fourth order accuracy in space and second order accuracy in time. The model is decomposed into small subdomains for each processor. After each processor updates wavefields within its domain, the processors exchange the wavefields via message passing interface (MPI). The parallel implementation reduces the computational time and also allows one to study larger problems. From our numerical experiment, consistent with other I-D experiments, it is found that the presence of heterogeneity of porous medium can produce significant P-wave attenuation in the seismic frequency range.

Subject Keywords

Biot's theory, Poroelastic media, Wave propagation, Finite difference, Message passing interface (MPI)

URI

https://hdl.handle.net/11511/39763

Journal

COMPUTERS & GEOSCIENCES

DOI

https://doi.org/10.1016/j.cageo.2005.11.001

Collections

Department of Civil Engineering, Article

Suggestions

OpenMETU
Core

Fast Algorithms for Digital Computation of Linear Canonical Transforms Koc, Aykut; Öktem, Sevinç Figen; Ozaktas, Haldun M.; Kutay, M. Alper (2016-01-01) Fast and accurate algorithms for digital computation of linear canonical transforms (LCTs) are discussed. Direct numerical integration takes O.N-2/time, where N is the number of samples. Designing fast and accurate algorithms that take O. N logN/time is of importance for practical utilization of LCTs. There are several approaches to designing fast algorithms. One approach is to decompose an arbitrary LCT into blocks, all of which have fast implementations, thus obtaining an overall fast algorithm. Another a...
Assessment of alternative simulation techniques in nonlinear time history analyses of multi-story frame buildings: A case study Karim Zadeh Naghshineh, Shaghayegh; Askan Gündoğan, Ayşegül; Yakut, Ahmet (2017-07-01) In regions with sparse ground motion data, simulations provide alternative acceleration time series for evaluation of the dynamic response of a structure. Different ground motion simulation methods provide varying levels of goodness of fit between observed and synthetic data. Before using the seismologically acceptable synthetic records for engineering purposes, it is critical to investigate the efficiency of synthetics in predicting observed seismic responses of structures. For this purpose, in this study ...
Column Parallel Incremental Zoom ADC for Uncooled Imaging Applications Beyenir, Aycan; Akın, Tayfun; Department of Electrical and Electronics Engineering (2022-8-11) This thesis presents a new column-parallel analog to digital converter (ADC) integrated circuit with high precision (20-bit), high signal-to-noise ratio (119.5 dB), low power incremental Zoom ADC architecture and demonstrates its integration to a 384x288 pixel format analog output microbolometer readout integrated circuit (ROIC) with 25µm pixel pitch. The Zoom ADC is a hybrid of a low-resolution successive approximation register (SAR) ADC and a high-resolution Delta-Sigma ADC. The main aim of this design is...
Investigation of Alternative Array Configurations of Nanowires for Maximum Power Transmission at Optical Frequencies Karaosmanoglu, Bariscan; Satana, Hasan Aykut; Dikmen, Fatih; Ergül, Özgür Salih (2017-07-14) We present a computational analysis and comparison of nanowire arrays for efficient transmission of electromagnetic power at optical frequencies. A reliable simulation environment involving the multilevel fast multipole algorithm for fast iterative solutions of surface integral equations is used to accurately compute field intensity and power density values in the vicinity of nanowires in diverse configurations. The results demonstrate favorable properties of honeycomb structures that can provide efficient ...
Simulation of Directional Microphones in Digital Waveguide Mesh-Based Models of Room Acoustics Hacıhabiboğlu, Hüseyin; Günel Kılıç, Banu (2010-02-01) Digital waveguide mesh (DWM) models are time-domain numerical methods providing computationally simple solutions for wave propagation problems. They have been used in various acoustical modeling and audio synthesis applications including synthesis of musical instrument sounds and speech, and modeling of room acoustics. A successful model of room acoustics should be able to account for source and receiver directivity. Methods for the simulation of directional sources in DWM models were previously proposed. T...

Citation Formats

D.-H. Sheen, K. Tuncay, C.-E. Baag, and P. J. Ortoleva, “Parallel implementation of a velocity-stress staggered-grid finite-difference method for 2-D poroelastic wave propagation,” COMPUTERS & GEOSCIENCES, pp. 1182–1191, 2006, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/39763.