An asymptotic model for computing the scattering of waves by small heterogeneities in the time domain
|
|
- Polly Cornelia Hodges
- 5 years ago
- Views:
Transcription
1 An asymptotic model for computing the scattering of waves by small heterogeneities in the time domain Simon Marmorat EPI POems INRIA Saclay / CEA List Joint work with Patrick Joly (EPI POems INRIA/ENSTA/CNRS) Xavier Claeys (LJLL, Paris) Simon MARMORAT - Junior Seminar November, 19th 2013
2 Motivation and context 2 «Non-destructive testing (NDT) is a wide group of analysis techniques used in science and industry to evaluate the properties of a material, component or system without causing damage.» Wikipedia Applications: Medicine (ultrasonography, EEG, )! Geophysics (ground exploration)! Industrial maintenance (planes, nuclear plants, )
3 Non destructive testing 3 Study ultrasonic wave propagation with applications to Non-Destructive Testing. Emitter / receiver From the knowledge of the incident and the scattered waves, recover information (shape, size, physical properties, ) on the object of interest.
4 Non destructive testing 3 Study ultrasonic wave propagation with applications to Non-Destructive Testing. Emitter / receiver From the knowledge of the incident and the scattered waves, recover information (shape, size, physical properties, ) on the object of interest.
5 Non destructive testing 3 Study ultrasonic wave propagation with applications to Non-Destructive Testing. Emitter / receiver From the knowledge of the incident and the scattered waves, recover information (shape, size, physical properties, ) on the object of interest.
6 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t)
7 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Equation is said to be well-posed
8 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic
9 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic How to represent the phenomenon anyway?! Thanks to simulation!
10 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic How to represent the phenomenon anyway?! Thanks to simulation! Represent by a finite number of points, compute an approximation of the unknown on each points.
11 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic Mesh How to represent the phenomenon anyway?! Thanks to simulation! Represent by a finite number of points, compute an approximation of the unknown on each points.
12 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic Mesh Since equation time dependent, we also discretize time How to represent the phenomenon anyway?! Thanks to simulation! Represent by a finite number of points, compute an approximation of the unknown on each points. t
13 What is Numerical Analysis? 4 Equation modeling some physical phenomenon occurring on, with unknown u(x, t) Prove the equation has a unique solution Usually no analytical (no nice formula) solution Equation is said to be well-posed When the equation or the geometry! are not simplistic Mesh Since equation time dependent, we also discretize time How to represent the phenomenon anyway?! Thanks to simulation! Represent by a finite number of points, compute an approximation of the unknown on each points. Leads to a matrix system: solvable by a computer! t
14 What is Numerical Analysis? 5 Key parameter: the number of points on the mesh number of degrees of freedom n Mesh
15 What is Numerical Analysis? 5 Key parameter: the number of points on the mesh number of degrees of freedom n More degrees of freedom Better approximation Higher computational cost inverting a matrix of size has complexity O(n 3 ) n Mesh
16 What is Numerical Analysis? 5 Key parameter: the number of points on the mesh number of degrees of freedom n More degrees of freedom Better approximation Higher computational cost inverting a matrix of size has complexity O(n 3 ) n Mesh We note x the maximum diameter of mesh cells! and the time step. t
17 What is Numerical Analysis? 5 Key parameter: the number of points on the mesh number of degrees of freedom n More degrees of freedom Better approximation Higher computational cost inverting a matrix of size has complexity O(n 3 ) n Mesh We note x the maximum diameter of mesh cells! and the time step. t Remark: n ' 1 x
18 PhD motivations 6 Non destructive testing for ultrasonic waves in media containing small heterogeneities (for instance concrete with small aggregate in it). Propagation of waves inside smooth media disturbed by small defects with respect to the wavelength. Classical simulation tools struggles to deal with multi scale problems.
19 PhD motivations 7 Simulation of an elastic wave propagating inside a media disturbed by small inclusions Finite element method Simulation 35 time (µs) CPU time 3527 (s) Degrees of freedom t (µs)
20 PhD motivations 7 Simulation of an elastic wave propagating inside a media disturbed by small inclusions Finite element method Simulation 35 time (µs) CPU time 3527 (s) Degrees of freedom t (µs) CPU time for simulation in a medium of same size without defects is about 5 minutes!!!
21 PhD motivations 7 Simulation of an elastic wave propagating inside a media disturbed by small inclusions Finite element method Simulation 35 time (µs) CPU time 3527 (s) Degrees of freedom t (µs) CPU time for simulation in a medium of same size without defects is about 5 minutes!!!
22 PhD motivations 8 Homogeneous medium Medium with defects
23 PhD motivations 8 Homogeneous medium Medium with defects about degrees of freedom about degrees of freedom
24 PhD motivations 8 Homogeneous medium Medium with defects about degrees of freedom about degrees of freedom
25 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " <<
26 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy?
27 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy? Homogeneous medium x '
28 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy? Homogeneous medium Medium with defects x ' x ' "
29 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy? Homogeneous medium Medium with defects x ' x ' " Very bad since we have to choose (stability condition) t apple C x
30 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy? Homogeneous medium Medium with defects x ' x ' "
31 PhD motivations 9 We deal with time domain waves propagating in a smooth media with defects of characteristic size central wavelength! Assumption: " << How to choose x to reach satisfying accuracy? Homogeneous medium Medium with defects x ' x ' " take " What we want: x ' -scale into account via special numerical feature based on asymptotics
32 A model problem 10 Wave acoustic equation: look for u " (x, t) such that x 2 R d, t > 0.! "
33 A model problem 10 Wave acoustic equation: look for u " (x, t) such that x 2 R d, t > 0.! " rf = d 1 C A r ~g = dx i
34 A model problem 10 Wave acoustic equation: look for u " (x, t) such that x 2 R d, t > 0. The physical parameters ", µ " admit a jump on the defect boundary.! " Construction of an approximate model: defect taken into account via auxiliary sources.
35 MATCHED ASYMPTOTIC METHOD Simon MARMORAT - Junior Seminar November, 19th 2013
36 Matched asymptotics 12 Look for the solution u " of the following model problem: (P " ) 2 u 2 u " (x, 0) r (µ " ru " )=f, (x, 0) = 0, x 2 R d, t > 0, x 2 R d.! " Asymptotic behavior of u " as "! 0?
37 Matched asymptotics 12 Look for the solution u " of the following model problem: (P " ) 2 u 2 u " (x, 0) r (µ " ru " )=f, (x, 0) = 0, x 2 R d, t > 0, x 2 R d. Asymptotic behavior of as? We would like something like u " "! 0 u " (x, t) =! " 1X n=0 and compute the u n (independent of " ). " n u n (x, t)
38 Matched asymptotics 12 Look for the solution u " of the following model problem: (P " ) 2 u 2 u " (x, 0) r (µ " ru " )=f, (x, 0) = 0, x 2 R d, t > 0, x 2 R d.! " Asymptotic behavior of as? u " "! 0 No uniform development of u " in function of "! Need to distinguish a far-field zone and a near-field zone inside which development holds.
39 Matched asymptotics 13 Look for the solution u " of the following model problem: (P " ) 2 u 2 u " (x, 0) r (µ " ru " )=f, (x, 0) = 0, x 2 R d, t > 0, x 2 R d.! "
40 Matched asymptotics 13 Look for the solution u " of the following model problem: (P " ) 2 u 2 u " (x, 0) r (µ " ru " )=f, (x, 0) = 0, x 2 R d, t > 0, x 2 R d. We denote! " " (x) = (x) if x 2!" + (x) if x 2 R d \! " µ " (x) = µ (x) if x 2!" µ + (x) if x 2 R d \! "
41 Matched asymptotics 14 Model problem 2 u 2 r (µ " ru " )=f,! "
42 Matched asymptotics 14 Model problem 2 u 2 r (µ " ru " )=f, p "! 2 " p " Near field zone Matching zone Far field zone
43 Ansatz and total field 15 Select a cutoff function,! rp scale it "(r) =, "! and represent the total field as the sum (r) u " (x, t) = " ( x ) U N " (x, t)+(1 " ( x )) u F " (x, t)
44 Ansatz and total field 15 Select a cutoff function,! rp scale it "(r) =, "! and represent the total field as the sum (r) u " (x, t) = " ( x ) U N " (x, t)+(1 " ( x )) u F " (x, t) = when x apple p " 2
45 Ansatz and total field 15 Select a cutoff function,! rp scale it "(r) =, "! and represent the total field as the sum (r) u " (x, t) = " ( x ) U N " (x, t)+(1 " ( x )) u F " (x, t) = when x p "
46 Matching principle 16 Far field equations u 2 r µ + ru 0 = f, Limit field u 2 r µ + ru n =0, n 1. Near field equations r (µru n )=0, n =0, 1 r (µru n U n 2, n 2
47 Matching principle 16 Far field equations u 2 r µ + ru 0 = f, Limit field Near field equations r (µru n )=0, n =0, 1 u 2 r µ + ru n =0, n 1. NX Matching principle: NX The partial sums " n u n and " n U n must coïncide in the matching zone, n=0 n=0 r (µru n U n 2, n 2 8N 0.
48 Matching principle 16 Far field equations u 2 r µ + ru 0 = f, Limit field Near field equations r (µru n )=0, n =0, 1 u 2 r µ + ru n =0, n 1. NX Matching principle: NX The partial sums " n u n and " n U n must coïncide in the matching zone, n=0 n=0 r (µru n U n 2, n 2 8N 0. Equations + matching principle the U n s can be computed.
49 TOWARDS AN ASYMPTOTIC METHOD Simon MARMORAT - Junior Seminar November, 19th 2013
50 An asymptotic model 18 From the initial model problem (P " ) 2 u 2 r (µ " ru " )=f, x 2 R d, t > 0,
51 An asymptotic model 18 From the initial model problem (P " ) 2 u 2 r (µ " ru " )=f, x 2 R d, t > 0, seeing ", µ " as perturbations of +, µ + we derive the following approximate model:
52 An asymptotic model 18 From the initial model problem (P " ) 2 u 2 r (µ " ru " )=f, x 2 R d, t > 0, seeing ", µ " as perturbations of +, µ + we derive the following approximate model: auxiliary source ( e P " ) u 2 r µ + ru " = f + S" N v ", v " = S N "? u", x 2 R d, t > 0, coupling the auxiliary source with the unknown
53 An asymptotic model 18 From the initial model problem " - wave operator, not suited for discretization (P " ) 2 u 2 r (µ " ru " )=f, x 2 R d, t > 0, seeing ", µ " as perturbations of +, µ + we derive the following approximate model: ( e P " ) u 2 r µ + ru " = f + S" N v ", v " = S N "? u", " - independent wave operator, well suited for discretization S" N is an operator involving the first N near-field functions N is the order of approximation of the method: ksol (P )" sol ( e P )" kapplec" N x 2 R d, t > 0,
54 An asymptotic model 18 From the initial model problem " - wave operator, not suited for discretization (P " ) 2 u 2 r (µ " ru " )=f, x 2 R d, t > 0, seeing ", µ " as perturbations of +, µ + we derive the following approximate model: ( e P " ) u 2 r µ + ru " = f + S" N v ", v " = S N "? u", " - independent wave operator, well suited for discretization S" N is an operator involving the first N near-field functions N is the order of approximation of the method: ksol (P )" sol ( e P )" kapplec" N x 2 R d, t > 0,
55 1D Numerical example 19 Physical parameters: + =1 = 10 µ = 10 µ + =1 " =0.04
56 1D Numerical example 19 Physical parameters: + =1 = 10 µ = 10 µ + =1 " =0.04 Computation of a reference solution via classical numerical method x =0.01
57 1D Numerical example 19 Physical parameters: + =1 = 10 µ = 10 µ + =1 " =0.04 Computation of a reference solution via classical numerical method x =0.01 Computation of an approximate solution via asymptotic numerical method x =0.04
58 1D Numerical example 20 Total field Diffracted field CPU time (s) Reference solution h ref = Asymptotic solution h as = Difference Relative L 2 error < 1%
59 1D Numerical example 20 Total field Diffracted field CPU time (s) Reference solution h ref = Asymptotic solution h as = Difference Relative L 2 error < 1%
60 Numerical example 21 Case of a unique defect Number of degrees of freedom CPU time 5 minutes!
61 Numerical example 21 Case of a unique defect Number of degrees of freedom CPU time 5 minutes!
62 Numerical example 21 Case of a unique defect Number of degrees of freedom CPU time 5 minutes!
63 Numerical example 21 Case of a unique defect Number of degrees of freedom CPU time 5 minutes!
64 Numerical example 22 Case of a network of defects Number of degrees of freedom CPU time 5 minutes!
65 Numerical example 22 Case of a network of defects Number of degrees of freedom CPU time 5 minutes!
66 Conclusion and prospects 23 Efficient and accurate numerical method for computing the scattering of waves by small defects. Rigorous error estimates.! To be done: Adapt the method to the 3D case. Extension to defects of arbitrary shape. Extension to the elastodynamic wave equation.
67 Thank you for your attention! Simon MARMORAT - Junior Seminar November, 19th 2013
MODELING OF ACOUSTIC PROCESSES IN SOLIDS BASED ON PARTICLE INTERACTION
MATEC Web of Conferences 55, 44 (8) IME&T 7 https://doi.org/.5/matecconf/85544 MODELING OF ACOUSTIC PROCESSES IN SOLIDS BASED ON PARTICLE INTERACTION Angela Kuzovova,, Timur Muksunov Tomsk State University,
More informationULTRASONIC WAVE PROPAGATION IN DISSIMILAR METAL WELDS APPLICATION OF A RAY-BASED MODEL AND COMPARISON WITH EXPERIMENTAL RESULTS
ULTRASONIC WAVE PROPAGATION IN DISSIMILAR METAL WELDS APPLICATION OF A RAY-BASED MODEL AND COMPARISON WITH EXPERIMENTAL RESULTS Audrey GARDAHAUT 1, Hugues LOURME 1, Frédéric JENSON 1, Shan LIN 2, Masaki
More informationWe G Elastic One-Return Boundary Element Method and Hybrid Elastic Thin-Slab Propagator
We G05 0 Elastic One-Return Boundary Element Method and Hybrid Elastic Thin-lab Propagator R.. Wu* (University of California) & Z. Ge (Beijing University, China) UMMARY We developed the theory and algorithm
More informationWeierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings
Weierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings Johannes Elschner & Guanghui Hu Mohrenstrasse 39 10117 Berlin Germany
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationEFIT SIMULATIONS FOR ULTRASONIC NDE
EFIT SIMULATIONS FOR ULTRASONIC NDE René Marklein, Karl-Jörg Langenberg, Klaus Mayer (University of Kassel, Department of Electrical and Computer Engineering, Electromagnetic Field Theory, Wilhelmshöher
More informationEVALUATION OF DAMAGES DUE TO ALKALI-SILICA REACTION WITH ACOUSTICS TECHNIQUES. DEVELOPMENT OF A NEW NONLINEAR METHOD.
EVALUATION OF DAMAGES DUE TO ALKALI-SILICA REACTION WITH ACOUSTICS TECHNIQUES. DEVELOPMENT OF A NEW NONLINEAR METHOD. Apedovi S. Kodjo (1, 2), Patrice Rivard (1), Frederic Cohen-Tenoudji (3) and Jean-Louis
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationOn the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar
NDT&E International 33 (2000) 401 407 www.elsevier.com/locate/ndteint On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar T.-T. Wu*, J.-H. Sun, J.-H.
More informationSimulation of an Ultrasonic Immersion Test for the Characterization of Anisotropic Materials
DIPARTIMENTO DI SCIENZE DELL INGEGNERIA CIVILE E DELL ARCHITETTURA POLITECNICO DI BARI Simulation of an Ultrasonic Immersion Test for the Characterization of Anisotropic Materials A. Castellano*, P. Foti,
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationSIMULATION OF ULTRASONIC NDT IN COMPOSITE RADIUS
SIMULATION OF ULTRASONIC NDT IN COMPOSITE RADIUS N. Dominguez 1, O. Grellou 2, S. Van-der-Veen 2 1 European Aeronautic Defense and Space Company (EADS), Innovation Works Dept., 1 rue Marius Terce, 325
More informationP137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients
P137 Our Experiences of 3D Synthetic Seismic Modeling with Tip-wave Superposition Method and Effective Coefficients M. Ayzenberg (StatoilHydro), A. Aizenberg (Institute of Petroleum Geology and Geophysics),
More informationEnhanced resolution in structured media
Enhanced resolution in structured media Eric Bonnetier w/ H. Ammari (Ecole Polytechnique), and Yves Capdeboscq (Oxford) Outline : 1. An experiment of super resolution 2. Small volume asymptotics 3. Periodicity
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationLECTURE NO. 4-5 INTRODUCTION ULTRASONIC * PULSE VELOCITY METHODS
LECTURE NO. 4-5 ULTRASONIC * PULSE VELOCITY METHODS Objectives: To introduce the UPV methods To briefly explain the theory of pulse propagation through concrete To explain equipments, procedures, calibrations,
More information9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR'S ACRONYM(S)
REPORT DOCUMENTATION PAGE 1 AFRL-SR-AR-TR-05_ The public reporting burden for this collection of information is estimated to average 1 hour per response, including tlde time gathering and maintaining the
More informationDoppler echocardiography & Magnetic Resonance Imaging. Doppler echocardiography. History: - Langevin developed sonar.
1 Doppler echocardiography & Magnetic Resonance Imaging History: - Langevin developed sonar. - 1940s development of pulse-echo. - 1950s development of mode A and B. - 1957 development of continuous wave
More informationApplication of Ultra-Sonic Pulse Velocity Techniques for Concrete Quality Measurements
Application of Ultra-Sonic Pulse Velocity Techniques for Concrete Quality Measurements Uday Bhise, Engineer. NDT; Durocrete Engineering Services Pvt. Ltd. Introduction This method of testing was originally
More informationAPPLICATION-DIRECTED MODELING OF RADIATION AND PROPAGATION OF ELASTIC WAVES IN ANISOTROPIC MEDIA: GPSS AND OPOSSM
APPLICATION-DIRECTED MODELING OF RADIATION AND PROPAGATION OF ELASTIC WAVES IN ANISOTROPIC MEDIA: GPSS AND OPOSSM M. Spies, F. Walte Fraunhofer-Institute for Nondestructive Testing (IzfP) 66123 Saarbriicken,
More informationTraveling Waves. Why there is reflection? The one-dimensional (1D) case. A traveling wave is the propagation of motion (disturbance) in a medium.
The one-dimensional (1D) case Traveling Waves A traveling wave is the propagation of motion (disturbance) in a medium. Reflection Why there is reflection? The perturbation propagates on. Traveling Wave
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-17 HARMONIC PLANE WAVES Introduction In this lecture, we discuss propagation of 1-D planar waves as they travel
More informationNon-Destructive Testing of Concrete Based on Analysis of Velocity Dispersion of Laser Ultrasonics
ECNDT 26 - Poster 223 Non-Destructive Testing of Concrete Based on Analysis of Velocity Dispersion of Laser Ultrasonics Kenichiro TSUYUKI, Ryuta KATAMURA, Satoru MIURA, Kajima Technical Research Institute,
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
More informationS. Marmorat (POEMS - INRIA Rocquencourt) F. Collino (CERFACS) J.-D. Benamou (POEMS - INRIA Rocquencourt)
Numerical MicroLocal Analysis (NMLA) S. Marmorat (POEMS - INRIA Rocquencourt) F. Collino (CERFACS) J.-D. Benamou (POEMS - INRIA Rocquencourt) Given frequency domain wave data, the proposed algorithm gives
More informationCRACK-TIP DIFFRACTION IN A TRANSVERSELY ISOTROPIC SOLID. A.N. Norris and J.D. Achenbach
CRACK-TIP DIFFRACTION IN A TRANSVERSELY ISOTROPIC SOLID A.N. Norris and J.D. Achenbach The Technological Institute Northwestern University Evanston, IL 60201 ABSTRACT Crack diffraction in a transversely
More informationNEW RESULTS ON TRANSMISSION EIGENVALUES. Fioralba Cakoni. Drossos Gintides
Inverse Problems and Imaging Volume 0, No. 0, 0, 0 Web site: http://www.aimsciences.org NEW RESULTS ON TRANSMISSION EIGENVALUES Fioralba Cakoni epartment of Mathematical Sciences University of elaware
More informationDOING PHYSICS WITH MATLAB WAVE MOTION
DOING PHYSICS WITH MATLAB WAVE MOTION THE [1D] SCALAR WAVE EQUATION THE FINITE DIFFERENCE TIME DOMAIN METHOD Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY
More informationThe failure of the sounding assumption in electroseismic investigations
The failure of the sounding assumption in electroseismic investigations F.D. Fourie, J.F. Botha Institute for Groundwater Studies, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
More informationStudy and design of a composite acoustic sensor to characterize an heterogeneous media presenting a complex matrix
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 Study and design of a composite acoustic sensor to characterize an heterogeneous media presenting a complex matrix PACS: 43.58.-e Georges,
More informationA PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging
A PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging Ali Fathalian and Kris Innanen Department of Geoscience, University of Calgary Summary The constant-q
More informationFinite element computation of leaky modes in straight and helical elastic waveguides
Finite element computation of leaky modes in straight and helical elastic waveguides Khac-Long Nguyen, Fabien Treyssede, Christophe Hazard, Anne-Sophie Bonnet-Ben Dhia To cite this version: Khac-Long Nguyen,
More informationSnell s law in transversely isotropic media using linearized group velocities and related quantities
Snell's law using group angles and velocities Snell s law in transversely isotropic media using linearized group velocities and related quantities P.F. Daley ABSTRACT Using a linearized approximation for
More informationWhen is the error in the h BEM for solving the Helmholtz equation bounded independently of k?
BIT manuscript No. (will be inserted by the editor) When is the error in the h BEM for solving the Helmholtz equation bounded independently of k? I. G. Graham M. Löhndorf J. M. Melenk E. A. Spence Received:
More informationElectromagnetic Acoustic Transducers for In and Out of plane Ultrasonic Wave Detection
7th World Conference on Nondestructive Testing, 5-8 Oct 8, Shanghai, China Electromagnetic Acoustic Transducers for In and Out of plane Ultrasonic Wave Detection Xiaoming JIAN, Steve DIXON, Karl QUIK Phoenix
More informationDetermining the effective electromagnetic properties of negative-refractive-index metamaterials from internal fields
Determining the effective electromagnetic properties of negative-refractive-index metamaterials from internal fields Bogdan-Ioan Popa* and Steven A. Cummer Department of Electrical and Computer Engineering,
More informationP125 Method for Calibrating Seismic Imaging Velocities
P125 Method for Calibrating Seismic Imaging Velocities P.C. Docherty* (Fairfield Industries Inc.) SUMMARY Anisotropy can cause large depth errors in images obtained with isotropic velocity analysis and
More informationInverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing
Inverse Scattering Theory: Transmission Eigenvalues and Non-destructive Testing Isaac Harris Texas A & M University College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with: F. Cakoni, H.
More informationGraduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto, Japan.
On relationship between contact surface rigidity and harmonic generation behavior in composite materials with mechanical nonlinearity at fiber-matrix interface (Singapore November 2017) N. Matsuda, K.
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 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 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 informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
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 informationDiscrete (Difference) Equations
Discrete (Difference) Equations A discrete, or difference, equation expresses a relationship between the elements of a sequence, {y n }, where n N 0 {0, 1, 2,...}. For example, the trivial difference equation
More informationQuantum Field Theory and Condensed Matter Physics: making the vacuum concrete. Fabian Essler (Oxford)
Quantum Field Theory and Condensed Matter Physics: making the vacuum concrete Fabian Essler (Oxford) Oxford, June 2013 Lev Landau This work contains many things which are new and interesting. Unfortunately,
More informationSeismic Waves and Earthquakes A Mathematical Overview
Seismic Waves and Earthquakes A Mathematical Overview The destruction caused by earthquakes is caused by the seismic waves propagating inside and in particular, near the surface of earth. We describe the
More informationApproximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL
147th Annual Meeting: Acoustical Society of America New York City, May 25, 2004 Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets Speaker: James G. Berryman University
More informationAdd-on unidirectional elastic metamaterial plate cloak
Add-on unidirectional elastic metamaterial plate cloak Min Kyung Lee *a and Yoon Young Kim **a,b a Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro, Gwanak-gu, Seoul,
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 informationOn double integrals over spheres
Inverse Problems 4 (1988) 1-10, Printed in the UK On double integrals over spheres R Burridge and G Beylkin Schlumberger-Doll Research, Old Quarry Road, Ridgefield, CT 06877-4108, USA Received 7 April
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationLamb Wave Behavior in Bridge Girder Geometries
Lamb Wave Behavior in Bridge Girder Geometries I. J. Oppenheim a*, D. W. Greve b, N. L. Tyson a a Dept. of Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA 15213 b Dept.
More informationUNIVERSITY OF SOUTHAMPTON
UNIVERSITY OF SOUTHAMPTON PHYS2023W1 SEMESTER 1 EXAMINATION 2016-2017 WAVE PHYSICS Duration: 120 MINS (2 hours) This paper contains 9 questions. Answers to Section A and Section B must be in separate answer
More informationHomogenization of the Transmission Eigenvalue Problem for a Periodic Media
Homogenization of the Transmission Eigenvalue Problem for a Periodic Media Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work
More informationThe Helically Reduced Wave Equation as a Symmetric Positive System
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2003 The Helically Reduced Wave Equation as a Symmetric Positive System Charles G. Torre Utah State University Follow this
More informationImprovement of Low Strain Pile Integrity Test
Improvement of Low Strain Pile Integrity Test Wenzhang Luo 1, Fan Chen 2, and Junling Hu 1 1 Deparment of Mechanical Engineering, University of Bridgeport, Bridgeport, CT 06604 2 National Center for Quality
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 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 informationThe effect of anticlines on seismic fracture characterization and inversion based on a 3D numerical study
The effect of anticlines on seismic fracture characterization and inversion based on a 3D numerical study Yungui Xu 1,2, Gabril Chao 3 Xiang-Yang Li 24 1 Geoscience School, University of Edinburgh, UK
More informationElastodynamic response of a layered anisotropic plate : comparative approaches. Towards the development of hybrid methods
Journées Ondes du Sud-Ouest 10 mars 2016 Elastodynamic response of a layered anisotropic plate : comparative approaches. Towards the development of hybrid methods PhD thesis of Pierric MORA, successfully
More informationULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS
ULTRASONIC REFLECTION BY A PLANAR DISTRIBUTION OF SURFACE BREAKING CRACKS A. S. Cheng Center for QEFP, Northwestern University Evanston, IL 60208-3020 INTRODUCTION A number of researchers have demonstrated
More informationInterface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel
Interface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel Michael Schoenberg Schlumberger-Doll Research, Old Quarry
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 informationPinning and depinning of interfaces in random media
Pinning and depinning of interfaces in random media Patrick Dondl joint work with Nicolas Dirr and Michael Scheutzow March 17, 2011 at Université d Orléans signal under crossed polarizers is measured and
More informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde, Steve Arden Funded by: Leverhulme Trust University of Maryland, September 2005
More informationTheoretical and Experimental Investigation of Rayleigh Waves on Spherical and Cylindrical Surfaces
1st International Symposium on Laser Ultrasonics: Science, Technology and Applications July 16-18 2008, Montreal, Canada Theoretical and Experimental Investigation of Rayleigh Waves on Spherical and Cylindrical
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationA model to predict modal radiation by finite-sized sources in semi-infinite isotropic plates
Journal of Physics: Conference Series PAPER OPEN ACCESS A model to predict modal radiation by finite-sized sources in semi-infinite isotropic plates To cite this article: M Stévenin et al 207 J. Phys.:
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More information(1) (2) (3) Main Menu. Summary. reciprocity of the correlational type (e.g., Wapenaar and Fokkema, 2006; Shuster, 2009):
The far-field approximation in seismic interferometry Yingcai Zheng Yaofeng He, Modeling Imaging Laboratory, University of California, Santa Cruz, California, 95064. Summary Green s function retrieval
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 informationSpatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Spatio-Temporal Characterization of Bio-acoustic Scatterers in Complex Media Karim G. Sabra, School of Mechanical Engineering,
More informationA Novel Sensor Design for Generation and Detection of Shear-Horizontal Waves Based on Piezoelectric Fibres
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic A Novel Sensor Design for Generation and Detection of Shear-Horizontal Waves Based on Piezoelectric
More informationIntegral equations for crack systems in a slightly heterogeneous elastic medium
Boundary Elements and Other Mesh Reduction Methods XXXII 65 Integral equations for crack systems in a slightly heterogeneous elastic medium A. N. Galybin & S. M. Aizikovich Wessex Institute of Technology,
More informationSummary. Introduction
: a philosophical view on seismic interferometry Kees Wapenaar*, Delft University of Technology, and Roel Snieder, Colorado School of Mines Summary We discuss the phenomenon of turning noise into signal
More informationSurvey of Inverse Problems For Hyperbolic PDEs
Survey of Inverse Problems For Hyperbolic PDEs Rakesh Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA Email: rakesh@math.udel.edu January 14, 2011 1 Problem Formulation
More informationDetailed Program of the Workshop
Detailed Program of the Workshop Inverse Problems: Theoretical and Numerical Aspects Laboratoire de Mathématiques de Reims, December 17-19 2018 December 17 14h30 Opening of the Workhop 15h00-16h00 Mourad
More information= (G T G) 1 G T d. m L2
The importance of the Vp/Vs ratio in determining the error propagation and the resolution in linear AVA inversion M. Aleardi, A. Mazzotti Earth Sciences Department, University of Pisa, Italy Introduction.
More informationLamb Wave Propagation Modeling Using Cellular Automata
6th European Workshop on Structural Health Monitoring - We.3.E.2 More info about this article: http://www.ndt.net/?id=14171 Lamb Wave Propagation Modeling Using Cellular Automata P. KLUSKA 1, W. J. STASZEWSKI
More informationPart 1. The simple harmonic oscillator and the wave equation
Part 1 The simple harmonic oscillator and the wave equation In the first part of the course we revisit the simple harmonic oscillator, previously discussed in di erential equations class. We use the discussion
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
More informationINVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN. Leonid Lyubchyk
CINVESTAV Department of Automatic Control November 3, 20 INVERSE MODEL APPROACH TO DISTURBANCE REJECTION AND DECOUPLING CONTROLLER DESIGN Leonid Lyubchyk National Technical University of Ukraine Kharkov
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationSignal Loss. A1 A L[Neper] = ln or L[dB] = 20log 1. Proportional loss of signal amplitude with increasing propagation distance: = α d
Part 6 ATTENUATION Signal Loss Loss of signal amplitude: A1 A L[Neper] = ln or L[dB] = 0log 1 A A A 1 is the amplitude without loss A is the amplitude with loss Proportional loss of signal amplitude with
More informationAcoustic scattering : high frequency boundary element methods and unified transform methods
Acoustic scattering : high frequency boundary element methods and unified transform methods Book or Report Section Accepted Version Chandler Wilde, S.N. and Langdon, S. (2015) Acoustic scattering : high
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More information2D Modeling of Elastic Wave Propagation in Solids Containing Closed Cracks with Friction
2D Modeling of Elastic Wave Propagation in Solids Containing Closed Cracks with Friction S. Delrue 1, V. Aleshin 2, O. Bou Matar 2, K. Van Den Abeele 1 1 Wave Propagation & Signal Processing Research Group,
More informationNanoscale Energy Conversion and Information Processing Devices - NanoNice - Photoacoustic response in mesoscopic systems
Nanoscale Energy Conversion and Information Processing Devices - NanoNice - Photoacoustic response in mesoscopic systems Photonics group W. Claeys, S. Dilhair, S. Grauby, JM. Rampnoux, L. Patino Lopez,
More informationa~ a~ a~ a~ DEVELOPMENTS IN THE FINITE ELEMENT MODELING OF ULTRASONIC NDT PHENOMENA R. Ludwig and w. Lord Department of Electrical Engineering
DEVELOPMENTS N THE FNTE ELEMENT MODELNG OF ULTRASONC NDT PHENOMENA R. Ludwig and w. Lord Department of Electrical Engineering Colorado State University, Fort Collins. CO 8523 NTRODUCTON Numerical modeling
More informationLIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS. Suba Thomas Harry G. Kwatny Bor-Chin Chang
LIE TRANSFORM CONSTRUCTION OF EXPONENTIAL NORMAL FORM OBSERVERS Suba Thomas Harry G Kwatny Bor-Chin Chang Department of Mechanical Engineering and Mechanics, Drexel University, Philadelphia, PA 1914, USA
More information3 component laboratory experiments by laser interferometry: anisotropy estimations using polarization of quasi P-waves and S waves
3 component laboratory experiments by laser interferometry: anisotropy estimations using polarization of quasi P-waves and S waves Maxim Lebedev, Andrej Bóna, Roman Pevzner, Boris Gurevich 19.10.2011 Department
More informationOn convergence of Approximate Message Passing
On convergence of Approximate Message Passing Francesco Caltagirone (1), Florent Krzakala (2) and Lenka Zdeborova (1) (1) Institut de Physique Théorique, CEA Saclay (2) LPS, Ecole Normale Supérieure, Paris
More informationPoisson Regression. Ryan Godwin. ECON University of Manitoba
Poisson Regression Ryan Godwin ECON 7010 - University of Manitoba Abstract. These lecture notes introduce Maximum Likelihood Estimation (MLE) of a Poisson regression model. 1 Motivating the Poisson Regression
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 information