Articles


Document Type
Journal article (JA)
Title
A new upwind flux for a jump boundary condition applied to 3D viscous fracture modeling
Author
Zhan, Qiwei(1,2); Sun, Qingtao(1); Zhuang, Mingwei(1,4); Mao, Yiqian(1); Ren, Qiang(3); Fang, Yuan(1); Huang, Wei-Feng(1); Liu, Qing Huo(1)
Address
(1) Department of Electrical and Computer Engineering, Duke University, Durham; NC; 27708, United States; (2) Department of Civil and Environmental Engineering, Duke University, Durham; NC; 27708, United States; (3) School of Electronics and Information Engineering, Beihang University, Beijing; 100191, China; (4) Institute of Electromagnetics and Acoustics, Xiamen University, Xiamen, Fujian; 100191, China
RPAddress
Email
ResearchID
ORCID
Journal
Computer Methods in Applied Mechanics and Engineering
Publisher
Elsevier B.V.
ISSN
0045-7825
Published
2018-04-01, 331:456-473.
JCR
2
ImpactFactor
3.949
ISBN
Fund_Code
HYMC
HYDD
HYKSRQ
HYJSRQ
HYLWLB
HYJB
Keywords
Boundary conditions - Concrete beams and girders - Domain decomposition methods - Elastic waves - Fracture - Numerical methods - Wave propagation
Abstract
We present a discontinuous Galerkin (DG) algorithm with nonconformal meshes to simulate 3D elastic wave propagation in heterogeneous media with arbitrary discrete fractures. In our method, the fractures are not limited to be planar, single, and lossless, but can be curved, intersecting, and viscous. In contrast to the exact volumetric modeling for the extremely thin layer, explicitly treating an individual fracture as a geometry surface (i.e., an imperfect contact interface) requires the jump condition for displacement/velocity, but the continuity of traction vector on the fracture interface. A new upwind flux is proposed to weakly impose this jump boundary condition in the DG framework. This flux guarantees the stability and accuracy of the DG schemes to model arbitrary fractures. Unlike conventional Riemann solvers applied to continuous media, this solution involves an evolutionary update on the Godunov states. Besides this, no extra computational cost is added. In addition, we can extend the fracture interface into a perfectly matched layer to mimic an infinitely large fracture. Quantitative comparisons of the waveforms between our algorithm and an independent finite element code demonstrate the accuracy and efficiency of our algorithm. ? 2017 Elsevier B.V.
WOS Categories
Engineering, Multidisciplinary; Mathematics, Interdisciplinary Applications; Mechanics
Accession Number
WOS:000425737200020
EI收录号
20175204576409
DOI
10.1016/j.cma.2017.11.002
ESI_Type
COMPUTER
Collection
SCIE, EI

Back to List