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

Sheen, Dong-Hoon
Tuncay, Kağan
Baag, Chang-Eob
Ortoleva, Peter J.
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.


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 ...
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...
Wave Propagation in Poroelastic Media Saturated by Two Fluids
Tuncay, Kağan (ASME International, 1997-6-1)
<jats:p>A theory of wave propagation in isotropic poroelastic media saturated by two immiscible Newtonian fluids is presented. The macroscopic constitutive relations, and mass and momentum balance equations are obtained by volume averaging the microscale balance and constitutive equations and assuming small deformations. Momentum transfer terms are expressed in terms of intrinsic and relative permeabilities assuming the validity of Darcy’s law. The coefficients of macroscopic constitutive relations are expr...
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: