Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
|
|
- Grace Hortense Malone
- 5 years ago
- Views:
Transcription
1 Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM) WADG 6/12/18 1 / 27
2 Collaborators and acknowledgements Prof. T. Warburton: Virginia Tech, Department of Mathematics Dr. Russell Hewett: TOTAL Research and Technology USA Prof. Maarten de Hoop: Rice University, Department of Computational Mathematics Khemraj Shukla: Oklahoma State University, Dept. of Geophysics. Chan (CAAM) WADG 6/12/18 2 / 27
3 High order DG methods for wave propagation Unstructured (tetrahedral) meshes for geometric flexibility. Low numerical dissipation and dispersion. High order approximations: more accurate per unknown. Explicit time stepping: high performance on many-core. Figure courtesy of Axel Modave. Chan (CAAM) WADG 6/12/18 3 / 27
4 High order DG methods for wave propagation Unstructured (tetrahedral) meshes for geometric flexibility. Low numerical dissipation and dispersion. High order approximations: more accurate per unknown. Explicit time stepping: high -0.8 performance on many-core Fine linear approximation. Chan (CAAM) WADG 6/12/18 3 / 27
5 High order DG methods for wave propagation Unstructured (tetrahedral) meshes for geometric flexibility. Low numerical dissipation and dispersion. High order approximations: more accurate per unknown. Explicit time stepping: high -0.8 performance on many-core Coarse quadratic approximation. Chan (CAAM) WADG 6/12/18 3 / 27
6 High order DG methods for wave propagation Unstructured (tetrahedral) meshes for geometric flexibility. Low numerical dissipation and dispersion. High order approximations: more accurate per unknown. Explicit time stepping: high performance on many-core. Max errors vs. dofs. Chan (CAAM) WADG 6/12/18 3 / 27
7 High order DG methods for wave propagation Unstructured (tetrahedral) meshes for geometric flexibility. Low numerical dissipation and dispersion. High order approximations: more accurate per unknown. Explicit time stepping: high performance on many-core. Graphics processing units (GPU). Chan (CAAM) WADG 6/12/18 3 / 27
8 Outline Weight-adjusted DG methods: acoustics 1 Weight-adjusted DG methods: acoustics 2 Extension to elastic wave propagation 3 Acoustic-elastic coupling, poroelasticity Chan (CAAM) WADG 6/12/18 3 / 27
9 Outline Weight-adjusted DG methods: acoustics 1 Weight-adjusted DG methods: acoustics 2 Extension to elastic wave propagation 3 Acoustic-elastic coupling, poroelasticity Chan (CAAM) WADG 6/12/18 3 / 27
10 Weight-adjusted DG methods: acoustics High order methods and arbitrary heterogeneous media (a) Mesh and exact c 2 (b) Piecewise const. c 2 (c) High order c 2 Efficient implementations on triangular or tetrahedral meshes assume piecewise constant wavespeed c 2. High order vs constant approx. of media: spurious reflections. This talk: generalized mass lumping for DG with provable stability, high order accuracy, simple numerical fluxes (central flux + penalty). Chan (CAAM) WADG 6/12/18 4 / 27
11 Weight-adjusted DG methods: acoustics Energy stable discontinuous Galerkin formulations Model problem: acoustic wave equation (pressure-velocity system) 1 p c 2 (x) t = u, u t = p Local formulation on element D k 1 p D k c 2 (x) t q = uq + 1 ([[u]] n + τ p [[p]]) q D k 2 D k u t v = p v + 1 ([[p]] + τ u [[u]] n) v D k 2 D k D k Energy stability (penalty terms weakly enforce continuity conditions) ( ) Dk p 2 t c 2 (x) + u 2 = τ p [[p]] 2 + τ u [[u n]] 2 0. k k D k Chan (CAAM) WADG 6/12/18 5 / 27
12 Weight-adjusted DG methods: acoustics Semi-discrete formulation and weighted mass matrices Weighted mass matrix M 1/c 2: SPD, induces weighted norm on p ( M1/c 2) ij = 1 D k c 2 (x) φ j(x)φ i (x) dx. Semi-discrete form: face mass and weak derivative matrices M f, S i. dp M 1/c 2 dt = Must build and factorize M 1/c 2 d i=1 S i u i + faces M du i dt = S ip + M f F u. faces c 2 constant and ( M 1/c 2) 1 = c 2 M 1 ). M f F p, separately over each element (unless Chan (CAAM) WADG 6/12/18 6 / 27
13 Weight-adjusted DG methods: acoustics Weight-adjusted DG: generalized mass lumping Weight-adjusted DG (WADG): M ( M 1/w ) 1 M is SPD and an energy stable approximation of weighted mass matrix. du M w dt M ( ) 1 du M 1/w M dt Bypasses inverse of weighted matrix (M w ) 1 = right hand side. ( M ( M 1/w ) 1 M ) 1 = M 1 M 1/w M 1. Low storage matrix-free application of M 1 M 1/w using quadrature. M = V T q WV q, V q = eval. at quad. pts, W = quad. weights (M) 1 M 1/w RHS = M 1 V T q W }{{} P q diag (1/w) V q (RHS). Chan (CAAM) WADG 6/12/18 7 / 27
14 Weight-adjusted DG methods: acoustics WADG and high order accuracy Generates norm with same equivalence constants w min u T Mu u T M w u w max u T Mu High order accuracy of weighted projection P w, WADG projection P w u P w u C w h N+1 w L 2 W N+1, u W N+1,2 P w u P w u C w,nh N+2 w L 2 W N+1, u W N+1,2. WADG retains high order accuracy for moments: if v P M v T M w u v T M 1 M 1/w M 1 u C w w W N+1, h 2N+2 M u W N+1,2 v L 2. Chan (CAAM) WADG 6/12/18 8 / 27
15 Weight-adjusted DG methods: acoustics WADG and high order accuracy Generates norm with same equivalence constants w min u T Mu u T M 1 M 1/w M 1 u w max u T Mu High order accuracy of weighted projection P w, WADG projection P w u P w u C w h N+1 w L 2 W N+1, u W N+1,2 P w u P w u C w,nh N+2 w L 2 W N+1, u W N+1,2. WADG retains high order accuracy for moments: if v P M v T M w u v T M 1 M 1/w M 1 u C w w W N+1, h 2N+2 M u W N+1,2 v L 2. Chan (CAAM) WADG 6/12/18 8 / 27
16 Weight-adjusted DG methods: acoustics WADG and high order accuracy Generates norm with same equivalence constants w min u T Mu u T M 1 M 1/w M 1 u w max u T Mu High order accuracy of weighted projection P w, WADG projection P w u P w u C w h N+1 w L 2 W N+1, u W N+1,2 P w u P w u C w,nh N+2 w L 2 W N+1, u W N+1,2. WADG retains high order accuracy for moments: if v P M v T M w u v T M 1 M 1/w M 1 u C w w W N+1, h 2N+2 M u W N+1,2 v L 2. Chan (CAAM) WADG 6/12/18 8 / 27
17 Weight-adjusted DG methods: acoustics WADG and high order accuracy Generates norm with same equivalence constants w min u T Mu u T M 1 M 1/w M 1 u w max u T Mu High order accuracy of weighted projection P w, WADG projection P w u P w u C w h N+1 w L 2 W N+1, u W N+1,2 P w u P w u C w,nh N+2 w L 2 W N+1, u W N+1,2. WADG retains high order accuracy for moments: if v P M v T M w u v T M 1 M 1/w M 1 u C w w W N+1, h 2N+2 M u W N+1,2 v L 2. Chan (CAAM) WADG 6/12/18 8 / 27
18 Weight-adjusted DG methods: acoustics Acoustic wave equation: variable coefficients (a) Weight w = c 2 (x, y) (b) Standard DG Figure: Standard vs. weight-adjusted DG with spatially varying c 2 containing both smooth variations and a discontinuity. Chan (CAAM) WADG 6/12/18 9 / 27
19 Weight-adjusted DG methods: acoustics Acoustic wave equation: variable coefficients (a) Weight w = c 2 (x, y) (b) Weighted-adjusted DG Figure: Standard vs. weight-adjusted DG with spatially varying c 2 containing both smooth variations and a discontinuity. Chan (CAAM) WADG 6/12/18 9 / 27
20 Weight-adjusted DG methods: acoustics Low storage WADG: efficiency on GPUs N = 1 N = 2 N = 3 N = 4 N = 5 N = 6 N = 7 Mw WADG Speedup Time (ns) per element: storing/applying M 1 1/w vs WADG (deg. 2N quadrature). Efficiency on GPUs: reduce memory accesses and data movement. (Tuned) low storage WADG faster than storing and applying M 1 1/w! Chan (CAAM) WADG 6/12/18 10 / 27
21 Outline Extension to elastic wave propagation 1 Weight-adjusted DG methods: acoustics 2 Extension to elastic wave propagation 3 Acoustic-elastic coupling, poroelasticity Chan (CAAM) WADG 6/12/18 10 / 27
22 Extension to elastic wave propagation Matrix-valued weights and elastic wave propagation Symmetric hyperbolic system: velocity v, stress σ = vec(s) with σ = (σ xx, σ yy, σ zz, σ yz, σ xz, σ xy ) T. ρ v t = d i=1 A T i σ x i, 1 σ C t = d i=1 A i v x i. Factoring out the constitutive stiffness tensor C results in simple and spatially constant matrices A i µ + λ λ λ C = λ 2µ + λ λ λ λ 2µ + λ, A 1 = µi }{{} for isotropic media Hughes and Marsden Classical elastodynamics as a linear symmetric hyperbolic system. Chan (CAAM) WADG 6/12/18 11 / 27
23 Extension to elastic wave propagation DG formulation for elasticity: energy stability D k D k Analogous to acoustics: numerical fluxes independent of media! ( ) 1 σ T d v C q = A i q + 1 ( ) A n [[v]] + τ σ A n A T n [[σ]] q, t x i 2 ( ρ v ) T w = t D k D k i=1 d i=1 A T i Energy method: take q = σ, w = v D k ( C 1 σ D k D k σ w + 1 ( ) A T n [[σ]] + τ v A T n A n [[v]] w. x i 2 D k ) T σ = t t D ( k ρ v ) T v = t t ( ) σ(x) T C 1 (x)σ(x) D k ρ(x) v 2. Chan (CAAM) WADG 6/12/18 12 / 27
24 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) d v A i q + 1 ( ) A n [[v]] + τ σ A n A T n [[σ]] q, x i 2 D k i=1 d A T i D k i=1 σ x i w D k D k ( ) A T n [[σ]] + τ v A T n A n [[v]] w Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
25 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) d v A T q ( i + A n {{v}} + τ ) σ x i 2 A na T n [[σ]] q, D k i=1 d A T i D k i=1 σ x i w D k D k ( ) A T n [[σ]] + τ v A T n A n [[v]] w Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
26 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) d v A T σ ( i + A n {{v}} + τ ) σ x i 2 A na T n [[σ]] σ, D k i=1 d A T i D k i=1 σ x i v D k D k ( ) A T n [[σ]] + τ v A T n A n [[v]] v Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
27 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) D k Sum volume terms, which cancel to zero + ) σ ( A n {{v}} + τ σ 2 A na T n [[σ]] ( A T n [[σ]] + τ v A T n A n [[v]] ) v Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
28 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) D k Sum volume terms, which cancel to zero + ) σ ( A n {{v}} + τ σ 2 A na T n [[σ]] ( A T n [[σ]] + τ v A T n A n [[v]] ) v Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
29 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) Sum over all elements: central flux terms cancel + τ σ 2 A na T n [[σ]] σ + τ v 2 AT n A n [[v]] v k D k Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
30 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) Sum over all elements: combine penalty flux terms+ τ σ 2 AT n [[σ]] A T n σ + τ v 2 A n [[v]] A n v k D k Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
31 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) Sum over all elements: combine penalty flux terms+ τ σ 2 AT n [[σ]] A T n [[σ]] + τ v 2 A n [[v]] A n [[v]] k D k Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
32 Extension to elastic wave propagation DG formulation for elasticity: energy stability, cont. A i constant, integration by parts if quadrature exact for degree (2N 1) Sum over all elements: combine penalty flux terms+ τ σ 2 AT n [[σ]] A T n [[σ]] + τ v 2 A n [[v]] A n [[v]] k D k... and done! Energy stability (penalty weakly enforces continuity conditions): 1 ρ v 2 + σ T C 1 σ = 2 t D k D k τ v faces f 2 A n [[v]] 2 + τ σ A T n [[σ]] τ v, τ σ > 0 penalizes normal jumps [[v]] n 0, [[S]] n 0. Chan (CAAM) WADG 6/12/18 13 / 27
33 Extension to elastic wave propagation Semi-discrete system: matrix-valued weighted Matrix-weighted mass matrix: let W R d d be SPD with entries w ij M w11... M w1d M W =..... M wd1... M wdd Semi-discrete DG formulation involves C 1 -weighted mass matrix Σ M C 1 t = M ρi V t = d (A i S i ) V + ( ) I M f F σ, faces d ) (A T i S i Σ + ( ) I M f F v. faces i=1 i=1 Chan (CAAM) WADG 6/12/18 14 / 27
34 Extension to elastic wave propagation Weight-adjusted DG: matrix-valued weights Matrix-weighted mass matrix large, hard to invert M C 1... M 11 C 1 1d M C 1 =..... M C 1 d1... M C 1 dd Weight-adjusted approximation for C 1 decouples each component M 1 C 1 ( I M 1) M C ( I M 1 ). Evaluate RHS components at quadrature points, apply C(x i ) to component vectors at quadrature points, project back to polynomials. Recovers Kronecker product M 1 C 1 = C M 1 for constant C 1. Chan (CAAM) WADG 6/12/18 15 / 27
35 Extension to elastic wave propagation Energy stability and DG spectra (a) Central flux (τ v, τ σ = 0) (b) Penalty flux (τ v, τ σ = 1) Figure: DG spectra with random heterogeneities at each quadrature point. Guaranteed energy stability for both energy conservative and energy dissipative numerical interface fluxes. CFL can be improved by setting τ σ 1/ C, τ v ρ Chan (CAAM) WADG 6/12/18 16 / 27
36 Extension to elastic wave propagation Elastic wave propagation: convergence Convergence for harmonic oscillation, Rayleigh, Lamb, and Stoneley waves: between O(h N+1 ) and O(h N+1/2 ). σ error grows as C 1 (e.g. incompressible limit λ/µ ). Relative L 2 error Mesh size h (a) L 2 errors (Stoneley wave) v error σ error λ/µ (b) C 1, N = 3, h = 1/8. Chan (CAAM) WADG 6/12/18 17 / 27
37 Extension to elastic wave propagation Elastic wave propagation: stiff inclusion Figure: tr(σ) and σ xy for stiff inclusion with N = 5, h 1/50. Chan (CAAM) WADG 6/12/18 18 / 27
38 Extension to elastic wave propagation Elastic wave propagation: anisotropy Simple implementation for anisotropy - fluxes independent of C. (a) Vertical velocity, t = 30µs (b) Vertical velocity, t = 60µs Anisotropic heterogeneous media: transverse isotropy (x < 0), isotropy (x > 0). Komatitsch, Barnes, Tromp Simulation of anisotropic wave propagation based upon a spectral element method. Chan (CAAM) WADG 6/12/18 19 / 27
39 Extension to elastic wave propagation Curved meshes and heterogeneous media Map curved elements D k to reference element ˆD, introduce determinant of Jacobian of mapping J (varies of ˆD). RHS terms: discretize skew-symmetric form d v A T i ˆD i=1 d A T i ˆD i=1 q ( J + A n {{v}} + τ ) σ x i D k 2 A na T n [[σ]] q, σ x i wj D k ( ) A T n [[σ]] + τ v A T n A n [[v]] w. Time-derivative terms: J incorporated into matrix weighting. D k ( C ) 1 σ T ( q = JC t ˆD ) 1 σ T q. t Chan (CAAM) WADG 6/12/18 20 / 27
40 Extension to elastic wave propagation Curved meshes and random heterogeneous media (b) Central flux (τ v, τ σ = 0) (c) Central flux (τ v, τ σ = 1) A skew-symmetric discretization guarantees discrete energy stability. Chan (CAAM) WADG 6/12/18 21 / 27
41 Extension to elastic wave propagation Elastic wave propagation: 3D isotropic media (a) Computational mesh (b) Homogeneous isotropic media Figure: tr(σ) with µ(x) = 1 + H(y) cos(3πx) cos(3πy) cos(3πz), N = 5. Chan (CAAM) WADG 6/12/18 22 / 27
42 Extension to elastic wave propagation Elastic wave propagation: 3D isotropic media (a) Computational mesh (b) Piecewise constant C(x) Figure: tr(σ) with µ(x) = 1 + H(y) cos(3πx) cos(3πy) cos(3πz), N = 5. Chan (CAAM) WADG 6/12/18 22 / 27
43 Extension to elastic wave propagation Elastic wave propagation: 3D isotropic media (a) Computational mesh (b) High order C(x) Figure: tr(σ) with µ(x) = 1 + H(y) cos(3πx) cos(3πy) cos(3πz), N = 5. Chan (CAAM) WADG 6/12/18 22 / 27
44 Outline Acoustic-elastic coupling, poroelasticity 1 Weight-adjusted DG methods: acoustics 2 Extension to elastic wave propagation 3 Acoustic-elastic coupling, poroelasticity Chan (CAAM) WADG 6/12/18 22 / 27
45 Acoustic-elastic coupling, poroelasticity Acoustic-elastic coupling Coupling conditions for fluid (acoustic) and solid (elastic) media: S n = pn, v n = u n. Replace [[p]] and [[σ]] with residuals of coupling conditions [[p]] = n (S + n ) p, A T n [[σ]] = [[S]] n = p + n S n, (acoustic side) (elastic side). Energy stable for arbitrary acoustic and elastic media. Chan (CAAM) WADG 6/12/18 23 / 27
46 Acoustic-elastic coupling, poroelasticity Acoustic-elastic coupling: arbitrary heterogeneous media (a) Low order c 2 (x), µ(x) (b) tr(σ) Figure: Acoustic-elastic waves from a Ricker pulse (N = 10, h = 1/16). Chan (CAAM) WADG 6/12/18 24 / 27
47 Acoustic-elastic coupling, poroelasticity Acoustic-elastic coupling: arbitrary heterogeneous media (a) High order c 2 (x), µ(x) (b) tr(σ) Figure: Acoustic-elastic waves from a Ricker pulse (N = 10, h = 1/16). Chan (CAAM) WADG 6/12/18 24 / 27
48 Acoustic-elastic coupling, poroelasticity Poroelasticity: low frequency Biot s system Biot s system adds pore pressure p and relative fluid velocity q. Symmetric first order system with augmented stress, velocity. σ = (σ xx, σ yy, σ zz, σ yz, σ xz, σ xy, p) T, ṽ d E v t = A T σ i Dṽ, x i i=1 ṽ = (v x, v y, v z, q x, q y, q z ) T E σ σ t = d i=1 A i ṽ x i. E s, E v SPD, A i spatially const. + sparse (a single ±1 per row). ρ ρ f ( E v = C ρ f m 1, E 1 C s = 1 ) α α T C 1 1 M + αt C 1, α Lemoine, Ou, LeVeque High-resolution finite volume modeling of wave propagation in orthotropic poroelastic media. Chan (CAAM) WADG 6/12/18 25 / 27
49 Acoustic-elastic coupling, poroelasticity Numerical results for Biot (with de Hoop, Shukla) (a) Orthotropic media, 1.56 ms (b) Epoxy-glass, 1.8 ms (c) Isotropic interface, 230 ms Inviscid Biot solution (mass particle velocity) showing fast P, slow S, and slow P waves. Simulation utilizes N = 3 and a uniform mesh of 128 elements per side. Systematic derivation of WADG formulations: can symmetrize a system by defining an appropriate convex entropy (energy function). Carcione Wave propagation in anisotropic, saturated porous media: plane wave theory and numerical simulation. Lemoine, Ou, LeVeque High-resolution finite volume modeling of wave propagation in orthotropic poroelastic media. Chan (CAAM) WADG 6/12/18 26 / 27
50 Acoustic-elastic coupling, poroelasticity Summary and acknowledgements Weight-adjusted DG (WADG) for acoustic and elastic wave propagation in heterogeneous media and on curved meshes. Generalized mass lumping: energy stability and high order accuracy. Significantly simpler formulations using symmetric forms of PDEs. This work is supported by DMS and DMS Thank you! Questions? Chan, Hewett, Warburton Weight-adjusted DG methods: wave propagation in heterogeneous media (SISC). Chan Weight-adjusted DG methods: matrix-valued weights and elastic wave prop. in heterogeneous media (IJNME). Chan (CAAM) WADG 6/12/18 27 / 27
51 Acoustic-elastic coupling, poroelasticity Additional slides Chan (CAAM) WADG 6/12/18 27 / 27
52 Acoustic-elastic coupling, poroelasticity Computational costs at high orders of approximation Note: WADG at high orders becomes expensive! Runtime (s) WADG runtimes Volume Surface Update Large dense matrices: O(N 6 ) work per tet. High orders usually use tensor-product elements: O(N 4 ) vs O(N 6 ) cost, but less geometric flexibility Degree N Idea: choose basis such that matrices are sparse. WADG runtimes for 50 timesteps, elements. Chan (CAAM) WADG 6/12/18 27 / 27
53 Acoustic-elastic coupling, poroelasticity BBDG: efficient volume, surface kernels Nodal DG Bernstein-Bezier DG 6 Volume Surface 6 Volume Surface Runtime (s) Degree N Degree N Bernstein-Bezier speedup requires constant coefficient RHS evaluation! Chan (CAAM) WADG 6/12/18 27 / 27
54 Acoustic-elastic coupling, poroelasticity BBDG: efficient volume, surface kernels BBDG speedup over nodal DG 8 Volume Surface 6 Speedup Degree N Bernstein-Bezier speedup requires constant coefficient RHS evaluation! Chan (CAAM) WADG 6/12/18 27 / 27
55 Acoustic-elastic coupling, poroelasticity BBWADG: polynomial multiplication and projection (a) Exact c 2 (b) M = 0 approximation (c) M = 1 approximation WADG reuses fast Bernstein volume and surface kernels. Fast O(N 3 ) Bernstein algorithms for polynomial multiplication: represent c 2 P M, p(x) P N, and construct c 2 (x)p(x) P M+N. Fast O(N 4 ) polynomial L 2 projection P N : P M+N P N. Chan (CAAM) WADG 6/12/18 27 / 27
56 Acoustic-elastic coupling, poroelasticity BBWADG: computational runtime for M = 1 Runtime (s) BBWADG BBWADG 10 6 N 4 N 4 WADG WADG 10 7 N 6 N (a) Degree N Degree N (b) Chan (CAAM) WADG 6/12/18 27 / 27
57 Acoustic-elastic coupling, poroelasticity BBWADG: update kernel speedup over WADG (acoustics) N = 3 N = 4 N = 5 N = 6 N = 7 N = 8 WADG 1.60e e e e e e-6 BBWADG 2.20e e e e e e-7 Speedup For N 8, quadrature (and WADG) becomes much more expensive. (a) N = 7 quadrature (b) N = 8 quadrature Chan (CAAM) WADG 6/12/18 27 / 27
58 Acoustic-elastic coupling, poroelasticity BBWADG: approximating c 2 and accuracy L 2 error piecewise constant c 2 2 piecewise linear c 2 4 quadratic c 2 5 WADG (N = 4) Mesh size h M = 0 M = 1 M = 2 Approximating smooth c 2 (x) using L 2 projection: O(h 2 ) for M = 0, O(h 4 ) for M = 1, O(h M+3 ) for 0 < M N 2. Chan (CAAM) WADG 6/12/18 27 / 27
arxiv: v2 [math.na] 8 Dec 2017
Weight-adjusted discontinuous Galerkin methods: matrix-valued weights and elastic wave propagation in heterogeneous media Jesse Chan arxiv:70.0025v2 [math.na] 8 Dec 207 Abstract Weight-adjusted inner products
More informationPoroelasticity simulation in CLAWPACK for geophysical and biological applications
Poroelasticity simulation in CLAWPACK for geophysical and biological applications October 27, 2012 Joint work with Randall J. LeVeque (University of Washington) and Yvonne Ou (University of Delaware) Supported
More informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationarxiv: v2 [math.na] 20 Oct 2017
On the penalty stabilization mechanism for upwind discontinuous Galerkin formulations of first order hyperbolic systems Jesse Chan a,, T. Warburton b a Department of Computational and Applied Mathematics,
More informationarxiv: v2 [math.na] 21 Aug 2016
GPU-ACCELERATED BERNSTEIN-BEZIER DISCONTINUOUS GALERKIN METHODS FOR WAVE PROBLEMS JESSE CHAN AND T WARBURTON arxiv:15165v [mathna] 1 Aug 16 Abstract We evaluate the computational performance of the Bernstein-Bezier
More informationOperator Upscaling for the Wave Equation
Operator Upscaling for the Wave Equation Tetyana Vdovina Susan E. Minkoff UMBC), Oksana Korostyshevskaya Department of Computational and Applied Mathematics Rice University, Houston TX vdovina@caam.rice.edu
More informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More informationA nodal discontinuous Galerkin finite element method for the poroelastic wave equation
Comput. Geosci., manuscript No. (will be inserted by the editor) A nodal discontinuous Galerkin finite element method for the poroelastic wave equation Nodal DG for poroelastic system Khemraj hukla ( )
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationApproximation SEM-DG pour les problèmes d ondes elasto-acoustiques
Approximation SEM-DG pour les problèmes d ondes elasto-acoustiques Helene Barucq, Henri Calandra, Aurélien Citrain, Julien Diaz, Christian Gout To cite this version: Helene Barucq, Henri Calandra, Aurélien
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationLASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE
LASER GENERATED THERMOELASTIC WAVES IN AN ANISOTROPIC INFINITE PLATE H. M. Al-Qahtani and S. K. Datta University of Colorado Boulder CO 839-7 ABSTRACT. An analysis of the propagation of thermoelastic waves
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationDefmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states
Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationWe briefly discuss two examples for solving wave propagation type problems with finite differences, the acoustic and the seismic problem.
Excerpt from GEOL557 Numerical Modeling of Earth Systems by Becker and Kaus 2016 1 Wave propagation Figure 1: Finite difference discretization of the 2D acoustic problem. We briefly discuss two examples
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationMéthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel
Méthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel Julien Diaz, Équipe Inria Magique 3D Séminaire de Mathématiques Appliquées, Laboratoire de Mathématiques Nicolas Oresme
More informationWave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media
SECTION I. SEISMIC EXPLORATION Volume 38 WAVE FIELDS IN REAL MEDIA: Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media (SECOND EDITION, REVISED AND EXTENDED) by Jose M. CARCIONE
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationApplication of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling
1 Application of Nodal Discontinuous Glaerkin Methods in Acoustic Wave Modeling Xin Wang ABSTRACT This work will explore the discontinuous Galerkin finite element method (DG-FEM) for solving acoustic wave
More informationNotes: Outline. Diffusive flux. Notes: Notes: Advection-diffusion
Outline This lecture Diffusion and advection-diffusion Riemann problem for advection Diagonalization of hyperbolic system, reduction to advection equations Characteristics and Riemann problem for acoustics
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationTHE SPECTRAL ELEMENT METHOD (SEM): FORMULATION AND IMPLEMENTATION OF THE METHOD FOR ENGINEERING SEISMOLOGY PROBLEMS
THE SPECTRAL ELEMENT METHOD (SEM): FORMULATION AND IMPLEMENTATION OF THE METHOD FOR ENGINEERING SEISMOLOGY PROBLEMS Kristel C. Meza-Fajardo 1 and Apostolos S. Papageorgiou 2 1 Professor, Depto. de Ingenieria
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationA lumped nodal DGTD-PIC method to ensure charge conservation for the 3D Vlasov-Maxwell system on nonconforming grids
A lumped nodal DGTD-PIC method to ensure charge conservation for the 3D Vlasov-Maxwell system on nonconforming grids B. Fornet Nuclétudes Wednesday, October 13 1 Some Physics 2 DGTD-PIC : why and how 3
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationOn propagation of Love waves in an infinite transversely isotropic poroelastic layer
Journal of Physics: Conference Series PAPER OPEN ACCESS On propagation of Love waves in an infinite transversely isotropic poroelastic layer To cite this article: C Nageswara Nath et al 2015 J. Phys.:
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationMethod in 1D. Introduction of the Basic Concepts and Comparison to Optimal FD Operators
The Spectral Element Method in 1D Introduction of the Basic Concepts and Comparison to Optimal FD Operators Seismology and Seismics Seminar 17.06.2003 Ludwigs-Maximilians-Universität München by Bernhard
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More information12. Thermal and scalar field analysis
2. THERMAL AND SCALAR FIELD ANALYSIS... 2. INTRODUCTION... 2.2 GENERAL FIELD PROBLEM... 2 2.3 COMMON FLUX COMPONENTS... 5 2. GALERKIN INTEGRAL FORM... 6 2.5 GALERKIN INTEGRAL FORM*... 7 2.6 ORTHOTROPIC
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationHydraulic properties of porous media
PART 5 Hydraulic properties of porous media Porosity Definition: Void space: n V void /V total total porosity e V void /V solid Primary porosity - between grains Secondary porosity - fracture or solution
More informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationA model order reduction technique for speeding up computational homogenisation
A model order reduction technique for speeding up computational homogenisation Olivier Goury, Pierre Kerfriden, Wing Kam Liu, Stéphane Bordas Cardiff University Outline Introduction Heterogeneous materials
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationTrefftz-DG solution to the Helmholtz equation involving integral equations
Trefftz-DG solution to the Helmholtz equation involving integral equations H. Barucq, A. Bendali, M. Fares, V. Mattesi, S. Tordeux Magique 3D Inria Bordeaux Sud Ouest LMA UMR CNRS 5142 INSA of Toulouse
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationCME342 Parallel Methods in Numerical Analysis. Matrix Computation: Iterative Methods II. Sparse Matrix-vector Multiplication.
CME342 Parallel Methods in Numerical Analysis Matrix Computation: Iterative Methods II Outline: CG & its parallelization. Sparse Matrix-vector Multiplication. 1 Basic iterative methods: Ax = b r = b Ax
More informationq is large, to ensure accuracy of the spatial discretization 1 In practice, q = 6 is a good choice.
Time-stepping beyond CFL: a locally one-dimensional scheme for acoustic wave propagation Leonardo Zepeda-Núñez 1*), Russell J. Hewett 1, Minghua Michel Rao 2, and Laurent Demanet 1 1 Dept. of Mathematics
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
More informationMATH 304 Linear Algebra Lecture 19: Least squares problems (continued). Norms and inner products.
MATH 304 Linear Algebra Lecture 19: Least squares problems (continued). Norms and inner products. Orthogonal projection Theorem 1 Let V be a subspace of R n. Then any vector x R n is uniquely represented
More informationNumerical modeling of dynamic rupture propagation
IX. Slovenská geofyzikálna konferencia, 22. - 23. jún 2011 Fakulta matematiky, fyziky a informatiky UK, Bratislava Numerical modeling of dynamic rupture propagation Martin Galis Peter Moczo Jozef Kristek
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationThe Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
More informationMODELING OF SEISMIC WAVE SCATTERING FROM ROUGH FRACTURE SURFACES
Proceedings of the Project Review, Geo-Mathematical Imaging Group (Purdue University, West Lafayette IN), Vol. 1 (2011) pp. 243-252. MODELING OF SEISMIC WAVE SCATTERING FROM ROUGH FRACTURE SURFACES CHRISTOPHER
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationMaxwell s Equations:
Course Instructor Dr. Raymond C. Rumpf Office: A-337 Phone: (915) 747-6958 E-Mail: rcrumpf@utep.edu Maxwell s Equations: Physical Interpretation EE-3321 Electromagnetic Field Theory Outline Maxwell s Equations
More informationPhase Plane Behavior of Solitary Waves in Nonlinear Layered Media
Phase Plane Behavior of Solitary Waves in Nonlinear Layered Media Randall J. LeVeque 1 and Darryl H. Yong 2 1 Department of Applied Mathematics, University of Washington, Box 35242, Seattle, WA 98195-242
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationA new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation
A new second order absorbing boundary layer formulation for anisotropic-elastic wavefeld simulation Junxiao Li, Kris Innanen and Bing Wang University of Calgary China University of Petroleum-Beijing Summary
More informationImplementation of the dg method
Implementation of the dg method Matteo Cicuttin CERMICS - École Nationale des Ponts et Chaussées Marne-la-Vallée March 6, 2017 Outline Model problem, fix notation Representing polynomials Computing integrals
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationA 3 D finite element model for sound transmission through a double plate system with isotropic elastic porous materials
Acoustics and Vibrations Group Université de Sherbrooke, QC CANADA Département génie mécanique Université de Sherbrooke Sherbrooke, QC CANADA Tel.: (819) 821-7157 Fax: (819) 821-7163 A 3 D finite element
More informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
More information