Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
|
|
- Justin Martin
- 6 years ago
- Views:
Transcription
1 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 1 INRIA Bordeaux-Sud-Ouest, team-project Magique 3D 2 INRIA Sophia-Antipolis-Méditerranée, team-project Nachos 3 TOTAL Exploration-Production
2 Motivation Examples of seismic applications M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
3 Motivation Imaging method : the full wave inversion Iterative procedure Inverse problem requiring to solve a lot of forward problems M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
4 Motivation Imaging method : the full wave inversion Iterative procedure Inverse problem requiring to solve a lot of forward problems Seismic imaging : time-domain or harmonic-domain? Time-domain : imaging condition complicated but low computational cost Harmonic-domain : imaging condition simple but huge computational cost M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
5 Motivation Imaging method : the full wave inversion Iterative procedure Inverse problem requiring to solve a lot of forward problems Seismic imaging : time-domain or harmonic-domain? Time-domain : imaging condition complicated but low computational cost Harmonic-domain : imaging condition simple but huge computational cost Forward problem of the inversion process Elastic wave propagation in harmonic domain : Helmholtz equation Reduction of the size of the linear system M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
6 Motivation Seismic imaging in heterogeneous complex media Complex topography High heterogeneities M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
7 Motivation Seismic imaging in heterogeneous complex media Complex topography High heterogeneities Use of unstructured meshes with FE methods M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
8 Motivation Seismic imaging in heterogeneous complex media Complex topography High heterogeneities Use of unstructured meshes with FE methods DG method Flexible choice of interpolation orders (p adaptativity) Highly parallelizable method Increased computational cost as compared to classical FEM M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
9 Motivation Seismic imaging in heterogeneous complex media Complex topography High heterogeneities Use of unstructured meshes with FE methods DG method Flexible choice of interpolation orders (p adaptativity) Highly parallelizable method Increased computational cost as compared to classical FEM DOF of classical FEM DOF of DGM M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
10 Motivation Objective of this work Development of an hybridizable DG (HDG) method Comparison with a reference method : a standard nodal DG method Figure : Degrees of freedom of DGM Figure : Degrees of freedom of HDGM M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
11 HDG methods HDG methods B. Cockburn, J. Gopalakrishnan, R. Lazarov Unified hybridization of discontinuous Galerkin, mixed and continuous Galerkin methods for second order elliptic problems, SIAM Journal on Numerical Analysis, Vol. 47 (2009) S. Lanteri, L. Li, R. Perrussel, Numerical investigation of a high order hybridizable discontinuous Galerkin method for 2d time-harmonic Maxwell s equations, COMPEL, Vol. 32 (2013) (time-harmonic domain) N.C. Nguyen, J. Peraire, B. Cockburn, High-order implicit hybridizable discontinuous Galerkin methods for acoustics and elastodynamics, J. of Comput. Physics, Vol. 230 (2011) (time domain for seismic applications) M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
12 2D Helmholtz equations Contents 2D Helmholtz elastic equations Notations and definitions Hybridizable Discontinuous Galerkin method Numerical results Conclusions-Perspectives M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
13 2D Helmholtz equations 2D Helmholtz elastic equations First order formulation of Helmholtz wave equations x = (x, y) Ω R 2, { iωρ(x)v(x) = σ(x) + fs (x) iωσ(x) = C(x) ε(v(x)) Free surface condition : σn = 0 on Γ l Absorbing boundary condition : σn = v p (v n)n + v s (v t)t on Γ a v : velocity vector σ : stress tensor ε : strain tensor M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
14 2D Helmholtz equations 2D Helmholtz elastic equations First order formulation of Helmholtz wave equations x = (x, y) Ω R 2, { iωρ(x)v(x) = σ(x) + fs (x) iωσ(x) = C(x) ε(v(x)) Free surface condition : σn = 0 on Γ l Absorbing boundary condition : σn = v p (v n)n + v s (v t)t on Γ a ρ : mass density C : tensor of elasticity coefficients v p : P-wave velocity v s : S-wave velocity f s : source term, f s L 2 (Ω) M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
15 Definitions Contents 2D Helmholtz elastic equations Notations and definitions Hybridizable Discontinuous Galerkin method Numerical results Conclusions-Perspectives M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
16 Definitions Notations and definitions Notations T h mesh of Ω composed of triangles M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
17 Definitions Notations and definitions Notations T h mesh of Ω composed of triangles F h set of all faces F of T h F M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
18 Definitions Notations and definitions Notations T h mesh of Ω composed of triangles F h set of all faces F of T h n the normal outward vector of an element n M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
19 Definitions Notations and definitions Approximations spaces P p () set of polynomials of degree at most p on M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
20 Definitions Notations and definitions Approximations spaces P p () set of polynomials of degree at most p on V p h = {v ( L 2 (Ω) ) 2 : v V p () = (P p ()) 2, T h } M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
21 Definitions Notations and definitions Approximations spaces P p () set of polynomials of degree at most p on V p h = {v ( L 2 (Ω) ) 2 : v V p () = (P p ()) 2, T h } Σ p h = {σ ( L 2 (Ω) ) 3 : σ Σ p () = (P p ()) 3, T h } M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
22 Definitions Notations and definitions Approximations spaces P p () set of polynomials of degree at most p on V p h = {v ( L 2 (Ω) ) 2 : v V p () = (P p ()) 2, T h } Σ p h = {σ ( L 2 (Ω) ) 3 : σ Σ p () = (P p ()) 3, T h } M h = {η ( L 2 (F h ) ) 2 : η F (P p (F )) 2, F F h } M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
23 Definitions Notations and definitions Definitions Jump [[ ]] of a vector v through F : [[v]] = v + n + + v n = v + n + v n + Jump of a tensor σ through F : [[σ]] = σ + n + + σ n = σ + n + σ n + + n n + M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
24 HDG method Contents 2D Helmholtz elastic equations Notations and definitions Hybridizable Discontinuous Galerkin method Formulation Algorithm Numerical results Conclusions-Perspectives M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
25 HDG method Formulation HDG formulation of the equations Local HDG formulation { iωρv σ = 0 iωσ Cε (v) = 0 M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
26 HDG method Formulation HDG formulation of the equations Local HDG formulation iωρ v w + iωσ : ξ + v σ : w ( C ξ σ n w = 0 ) v C ξ n = 0 σ and v are numerical traces of σ and v respectively on M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
27 HDG method Formulation HDG formulation of the equations Local HDG formulation We define : iωρ v w + iωσ : ξ + v σ : w ( C ξ σ n w = 0 ) v C ξ n = 0 v F = λ F, F F h, σ n = σ n τi ( v λ ), on where τ is the stabilization parameter (τ > 0) M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
28 HDG method Formulation HDG formulation of the equations Local HDG formulation We replace v and ( σ n ) by their definitions into the local equations iωρ v w + σ : w σ n w ( ) + τi v λ w = 0 ( ) iωσ : ξ + v C ξ λ C ξ n = 0 M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
29 HDG method Formulation HDG formulation of the equations Local HDG formulation iωρ v w iωσ : ξ + v ( σ ) w + ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
30 HDG method Formulation HDG formulation of the equations Local HDG formulation iωρ v w iωσ : ξ + v ( σ ) w + ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 We define : W = ( v x, v z, σ xx, σ zz, σ ) T xz Λ = ( ) T Λ F1, Λ F2,..., Λ Fn f, where nf = card(f h ) Discretization of the local HDG formulation A W + C Λ = 0 M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
31 HDG method Formulation HDG formulation of the equations Transmission condition In order to determine λ, the continuity of the normal component of σ is weakly enforced, rendering this numerical trace conservative : [[ σ n]] η = 0 F M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
32 HDG method Formulation HDG formulation of the equations Transmission condition In order to determine λ, the continuity of the normal component of σ is weakly enforced, rendering this numerical trace conservative : [[ σ n]] η = 0 F Replacing ( σ n ) and summing over all faces, the transmission condition becomes : ( ) σ n η ( ) τi v λ η = 0 T h T h M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
33 HDG method Formulation HDG formulation of the equations Transmission condition In order to determine λ, the continuity of the normal component of σ is weakly enforced, rendering this numerical trace conservative : [[ σ n]] η = 0 F Replacing ( σ n ) and summing over all faces, the transmission condition becomes : ( ) σ n η ( ) τi v λ η = 0 T h T h Discretization of the transmission condition [ B W + L Λ ] = 0 T h M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
34 HDG method Formulation HDG formulation of the equations Global HDG formulation iωρ v w iωσ : ξ + T h ( σ ) w + v ( σ n ) η T h ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 τi ( v λ ) η = 0 M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
35 HDG method Formulation HDG formulation of the equations Global HDG formulation iωρ v w iωσ : ξ + T h ( σ ) w + v ( σ n ) η T h ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 τi ( v λ ) η = 0 Global HDG discretization A W + C Λ = 0 [ B W + L Λ ] = 0 T h M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
36 HDG method Formulation HDG formulation of the equations Global HDG formulation iωρ v w iωσ : ξ + T h ( σ ) w + v ( σ n ) η T h ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 τi ( v λ ) η = 0 Global HDG discretization W = (A ) 1 C Λ [ B W + L Λ ] = 0 T h M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
37 HDG method Formulation HDG formulation of the equations Global HDG formulation iωρ v w iωσ : ξ + T h ( σ ) w + v ( σ n ) η T h ( ) C ξ τi ( v λ ) w = 0 λ C ξ n = 0 τi ( v λ ) η = 0 Global HDG discretization [ B (A ) 1 C + L ] Λ = 0 T h M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
38 HDG method Algorithm Main idea of the algorithm using the HDG formulation 1. Construction of the linear system MΛ = 0 with M = T h [ B (A ) 1 C + L ] for = 1 to Nb tri do Compute matrices B, (A ) 1, C and L Construction of the corresponding section of end for M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
39 HDG method Algorithm Main idea of the algorithm using the HDG formulation 1. Construction of the linear system MΛ = 0 2. Construction of the right hand side S M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
40 HDG method Algorithm Main idea of the algorithm using the HDG formulation 1. Construction of the linear system MΛ = 0 2. Construction of the right hand side S 3. Resolution MΛ = S M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
41 HDG method Algorithm Main idea of the algorithm using the HDG formulation 1. Construction of the linear system MΛ = 0 2. Construction of the right hand side S 3. Resolution MΛ = S 4. Computation of the solutions of the initial problem for = 1 to Nb tri do Compute W = (A ) 1 C Λ end for M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
42 Numerical results Contents 2D Helmholtz elastic equations Notations and definitions Hybridizable Discontinuous Galerkin method Numerical results Plane wave in an homogeneous medium Marmousi test-case Conclusions-Perspectives M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
43 Numerical results Plane wave in an homogeneous medium Plane wave m Physical parameters : ρ = 2000kg.m 3 λ = 16GPa µ = 8GPa Plane wave : i(k cos θx+k sin θy) u = e m Computational domain Ω setting where k = ω v p θ = 0 Three meshes : 3000 elements elements elements M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
44 Numerical results Plane wave in an homogeneous medium Plane wave W a W e P 1 P 2 P 3 P h max Convergence order of the HDG scheme M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
45 Numerical results Marmousi test-case Marmousi test-case Computational domain Ω composed of triangles M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
46 Numerical results Marmousi test-case Parallel results for the Marmousi test-case with the HDG-P3 scheme, f = 2Hz M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
47 Numerical results Marmousi test-case Caracteristics of the computing processors used Plafrim platform Hardware specification : 16 nodes, 12 cores by nodes Caracteristics of computing nodes : 2 Hexa-core Westmere Intel R Xeon R X5670 Frequency : 2,93 GHz Cache L3 : 12 Mo RAM : 96 Go Infiniband DDR : 20Gb/s Ethernet : 1Gb/s M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
48 Numerical results Marmousi test-case Efficiency of the parallelism of the global matrix construction 1 Efficiency HDGm IPDGm Nb proc M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
49 Numerical results Marmousi test-case Efficiency of the parallelism of the whole algorithm 1 Efficiency HDGm IPDGm Nb proc M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
50 Numerical results Marmousi test-case Speed up for the global matrix construction 5 TimeIPDG/TimeHDG Nb proc M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
51 Numerical results Marmousi test-case Speed up (Total simulation time) 10 TimeIPDG/TimeHDG Nb proc M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
52 Numerical results Marmousi test-case Memory required (GB) for the simulation 120 Memory (GB) HDGm IPDGm Nb proc M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
53 Conclusion Contents 2D Helmholtz elastic equations Notations and definitions Hybridizable Discontinuous Galerkin method Numerical results Conclusions-Perspectives M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
54 Conclusion Conclusions-Perspectives Conclusions On a same mesh, with the HDG method : Memory gain Computational time gain M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
55 Conclusion Conclusions-Perspectives Conclusions On a same mesh, with the HDG method : Memory gain Computational time gain Perspectives Develop 3D HDG formulation for Helmholtz equations Solution strategy for the HDG linear system M. Bonnasse-Gahot - DG and HDG methods for Helmholtz wave equations March 5, /28
56 Conclusion Thank you!
Mé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 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 informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
More informationHybrid iterative-direct domain decomposition based solvers for the time-harmonic Maxwell equations
Hybrid iterative-direct domain decomposition based solvers for the time-harmonic Maxwell equations 1st workshop of the joint INRIA-UIUC laboratory for Petascale Computing Victorita Dolean 1, Mohamed El
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationJoseph Charles. INRIA Sophia Antipolis-Méditerranée Nachos project-team
Amélioration des performances de méthodes Galerkin discontinues d ordre élevé pour la résolution numérique des équations de Maxwell instationnaires sur des maillages simplexes Joseph Charles INRIA Sophia
More informationComparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping
Comparison of a One and Two Parameter Family of Transmission Conditions for Maxwell s Equations with Damping M. El Bouajaji 1, V. Dolean 2, M. J. Gander 3 and S. Lanteri 1 1 Introduction Transmission conditions
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 informationSTABILIZATION IN RELATION TO WAVENUMBER IN HDG METHODS
STABILIZATION IN RELATION TO WAVENUMBER IN HDG METHODS J. GOPALAKRISHNAN, S. LANTERI, N. OLIVARES, AND R. PERRUSSEL Abstract. Simulation of wave propagation through complex media relies on proper understanding
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 informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
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)
More informationA Multiscale Hybrid-Mixed method for the elastodynamic model in time domain
1/21 A Multiscale Hybrid-Mixed method for the elastodynamic model in time domain Weslley da Silva Pereira 1 Nachos project-team, Inria Sophia Antipolis - Mediterranée (FR) Inria Sophia Antipolis - Mediterranée
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
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 informationA hybridizable discontinuous Galerkin method for time-harmonic Maxwell s equations
A hybridizable discontinuous Galerkin method for time-harmonic Maxwell s equations Stéphane Lanteri, Liang Li, Ronan Perrussel To cite this version: Stéphane Lanteri, Liang Li, Ronan Perrussel. A hybridizable
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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationUtilisation de la compression low-rank pour réduire la complexité du solveur PaStiX
Utilisation de la compression low-rank pour réduire la complexité du solveur PaStiX 26 Septembre 2018 - JCAD 2018 - Lyon Grégoire Pichon, Mathieu Faverge, Pierre Ramet, Jean Roman Outline 1. Context 2.
More informationUnified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems.
Unified Hybridization Of Discontinuous Galerkin, Mixed, And Continuous Galerkin Methods For Second Order Elliptic Problems. Stefan Girke WWU Münster Institut für Numerische und Angewandte Mathematik 10th
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationFROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS
FROM GODUNOV TO A UNIFIED HYBRIDIZED DISCONTINUOUS GALERKIN FRAMEWORK FOR PARTIAL DIFFERENTIAL EQUATIONS TAN BUI-THANH Abstract. By revisiting the basic Godunov approach for linear system of hyperbolic
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
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 informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
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 informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationAn Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems
An Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems M. Baumann, and M.B. van Gijzen Email: M.M.Baumann@tudelft.nl Delft Institute of Applied Mathematics Delft University
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 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 informationCohesive band model: a triaxiality-dependent cohesive model for damage to crack transition in a non-local implicit discontinuous Galerkin framework
University of Liège Aerospace & Mechanical Engineering Cohesive band model: a triaxiality-dependent cohesive model for damage to crack transition in a non-local implicit discontinuous Galerkin framework
More informationComputational Seismology: An Introduction
Computational Seismology: An Introduction Aim of lecture: Understand why we need numerical methods to understand our world Learn about various numerical methods (finite differences, pseudospectal methods,
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationTo CG or to HDG: A Comparative Study
To CG or to HDG: A Comparative Study Robert M. Kirby a,1, Spencer J. Sherwin b,3 Bernardo Cockburn c,2, a School of Computing, Univ. of Utah, Salt Lake City, UT, USA b Department of Aeronautics, Imperial
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
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 informationNumerical Methods in Geophysics. Introduction
: Why numerical methods? simple geometries analytical solutions complex geometries numerical solutions Applications in geophysics seismology geodynamics electromagnetism... in all domains History of computers
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 informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
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 informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationextended Hybridizable Discontinuous Galerkin (X-HDG) for Void and Bimaterial Problems
extended Hybridizable Discontinuous Galerkin (X-HDG) for Void and Bimaterial Problems Ceren Gürkan, Esther Sala-Lardies, Martin Kronbichler, and Sonia Fernández-Méndez Abstract A strategy for the Hybridizable
More informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationHélène Barucq 1, Lionel Boillot 1, Henri Calandra 2 and Julien Diaz 1
ESAIM: PROCEEDINGS AND SURVEYS, September 04, Vol. 45, p. 400-409 J.-S. Dhersin, Editor ABSORBING BOUNDARY CONDITIONS FOR D TILTED TRANSVERSE ISOTROPIC ELASTIC MEDIA Hélène Barucq, Lionel Boillot, Henri
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationA free-vibration thermo-elastic analysis of laminated structures by variable ESL/LW plate finite element
A free-vibration thermo-elastic analysis of laminated structures by variable ESL/LW plate finite element Authors: Prof. Erasmo Carrera Dr. Stefano Valvano Bologna, 4-7 July 2017 Research group at Politecnico
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
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 & Department of Mathematics, University of
More informationThe Finite Cell Method: High order simulation of complex structures without meshing
The Finite Cell Method: High order simulation of complex structures without meshing E. Rank, A. Düster, D. Schillinger, Z. Yang Fakultät für Bauingenieur und Vermessungswesen Technische Universität München,
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationCohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework
University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition
More informationAn Error Estimator for the Finite Element Approximation of Plane and Cylindrical Acoustic Waves.
Copyright 205 Tech Science Press CMES, vol.06, no.2, pp.27-45, 205 An Error Estimator for the Finite Element Approximation of Plane and Cylindrical Acoustic Waves. J. E. Sebold, L. A. Lacerda 2 and J.
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
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 informationOn the diminishing of spurious oscillations in explicit finite element analysis of linear and non-linear wave propagation and contact problems
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16315 On the diminishing of spurious oscillations
More informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
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 informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationFast transient structural FE analysis imposing prescribed displacement condition. by using a model order reduction method via Krylov subspace
Journal of Applied Mechanics Vol.13 Vol.13, (August pp.159-167 21) (August 21) JSCE JSCE Fast transient structural FE analysis imposing prescribed displacement condition by using a model order reduction
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 informationAcoustics Analysis of Speaker ANSYS, Inc. November 28, 2014
Acoustics Analysis of Speaker 1 Introduction ANSYS 14.0 offers many enhancements in the area of acoustics. In this presentation, an example speaker analysis will be shown to highlight some of the acoustics
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationHybrid Discontinuous Galerkin Methods for Fluid Dynamics and Solid Mechanics
Hybrid Discontinuous Galerkin Methods for Fluid Dynamics and Solid Mechanics Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationMaps & Tents. Jay Gopalakrishnan. Portland State University. in collaboration with P. Monk, J. Schöberl, P. Sepúlveda, C.
Maps & Tents Jay Gopalakrishnan Portland State University in collaboration with P. Monk, J. Schöberl, P. Sepúlveda, C. Wintersteiger IMA, Minneapolis, June 2017 Jay Gopalakrishnan 1/46 Who are we celebrating?
More informationA SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS. Jennifer Proft CERMICS, ENPC. J. Proft CERMICS, ENPC
A SIMPLE TUTORIAL ON DISCONTINUOUS GALERKIN METHODS Jennifer Proft Outline Definitions A simple example Issues Historical development elliptic equations hyperbolic equations Discontinuous vs. continuous
More informationParallel Performance Studies for a Numerical Simulator of Atomic Layer Deposition Michael J. Reid
Section 1: Introduction Parallel Performance Studies for a Numerical Simulator of Atomic Layer Deposition Michael J. Reid During the manufacture of integrated circuits, a process called atomic layer deposition
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationA phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber
A phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber Chi Yeung Lam * and Chi-Wang Shu ** * Department of Mathematics, The Chinese University
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
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 informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationStress and Energy Transmission by Inhomogeneous Plane Waves into Dissipative Media
Purdue University Purdue e-pubs Publications of the Ray W. Herrick Laboratories School of Mechanical Engineering 11-6-2015 Stress and Energy Transmission by Inhomogeneous Plane Waves into Dissipative Media
More informationThe Evaluation of the Thermal Behaviour of an Underground Repository of the Spent Nuclear Fuel
The Evaluation of the Thermal Behaviour of an Underground Repository of the Spent Nuclear Fuel Roman Kohut, Jiří Starý, and Alexej Kolcun Institute of Geonics, Academy of Sciences of the Czech Republic,
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 Dept. of
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 informationModel Reduction via Discretization : Elasticity and Boltzmann
Model Reduction via Discretization : Elasticity and Boltzmann Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology IMA
More informationPrediction of High-frequency Vibro-acoustic Coupling in Anechoic Chamber Using Energy Finite Element Method and Energy Boundary Element Method
Copyright 2012 Tech Science ress CMES, vol.85, no.1, pp.65-78, 2012 rediction of High-frequency Vibro-acoustic Coupling in Anechoic Chamber Using Energy Finite Element Method and Energy Boundary Element
More informationFEMxDEM double scale approach with second gradient regularization applied to granular materials modeling
FEMxDEM double scale approach with second gradient regularization applied to granular materials modeling Albert Argilaga Claramunt Stefano Dal Pont Gaël Combe Denis Caillerie Jacques Desrues 16 december
More informationSpring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014)
Spring 2014: Computational and Variational Methods for Inverse Problems CSE 397/GEO 391/ME 397/ORI 397 Assignment 4 (due 14 April 2014) The first problem in this assignment is a paper-and-pencil exercise
More information