Research Statement. Tonatiuh Sánchez-Vizuet
|
|
- Darrell Farmer
- 5 years ago
- Views:
Transcription
1 I am especially interested in devising, rigorously analyzing and implementing methods for the numerical simulation of transient wave phenomena. During my doctoral studies I have focused on Boundary Integral Equation techniques for the simulation of time domain linear wave scattering and transmission problems. Prior to that I was involved in the simulation of nonlinear wave propagation by shock-capturing Godunov-style methods. Currently my main interest is the study and computational simulation of elastic wave propagation and its interaction with other kinds of waves and physical processes (acoustic, electric, poroelastic, fluid flow, etc.). I work in the design and analysis of numerical methods for their simulation using for instance boundary element, continuous and discontinuous Galerkin, or finite volume methods and the potential application of these techniques to direct and inverse problems of applied sciences and engineering. Boundary Integral Equations in the Time Domain and Convolution Quadrature Consider a linear time-independent elliptic partial differential operator L with an associated fundamental solution G(x,y). The solution u to the partial differential equation Lu = 0 in Ω R d, satisfying appropriate boundary conditions on Ω (and radiation conditions as x if Ω is unbounded) admits an integral representation in terms of the fundamental solution of the form u(x) = G(x,y)λ(y)dy ν(y) G(x,y)ϕ(y)dy, Ω where λ and ϕ are unknown density functions defined on Ω. The boundary conditions of the problem are then used to set up integral equations on the boundary from which the density functions can be retrieved. Boundary integral formulations exploit the possibility to transfer the volume problem onto the boundary of the region, effectively reducing the dimensionality of the problem by one. Moreover, once the equations have been discretized, all computations are done solely on the boundary, making these methods particularly well suited for the treatment of exterior problems in unbounded domains with homogeneous physical properties. In the time domain an analogous procedure can be followed resulting in boundary integral equations that have a convolutional structure. For instance, in R 3, the single layer potential for the wave equation has the following expression (S λ)(x,t) := Γ Ω λ(y,t c 1 x y ) dγ(y). 4π x y In this case, Convolution Quadrature [11, 12] (CQ) constitutes a very powerful tool for discretization. It relies on the observation that the convolution of a causal function g with a causal Laplace-transformable kernel f can be written in terms of its Laplace transform, F, and the solution to a first order evolution equation whose inhomogeneous term is exactly g. Formally, for some σ > 0 (f g)(t) = 1 2πi σ+i t σ i where y(t;s) is the unique solution of 0 F(s)e s(t τ) g(τ)dτ ds = 1 2πi ẏ(t) sy(t) = g(t), y(0) = 0. σ+i σ i F(s)y(t;s)ds, This differential equation can be approximately solved by any ODE-solving procedure using time-domain information from g. The numerical solution can then be used to approximate the complex contour integral using the generally more tractable Laplace-domain kernel F. Currently it is standard to use low order backward differentiation, trapezoidal rule or Runge-Kutta based CQ implementations for which 1
2 the behaviour of the error and convergence properties are well understood. In the near future I plan to explore the posibility of developing higher order schemes involving deferred correction to attain enhanced performance in applications of CQ to wave propagation. Traditionally, the stability and convergence analysis of the solutions obtained following this procedure depended on carefully bounding the operators of the Calderón projector associated to the PDE in terms of the Laplace parameter s and then applying results related to the inversion of the Laplace transform. In [8] we use a novel approach that sidesteps the Laplace-domain analysis using the machinery of abstract evolution equations and semi-group theory to obtain time-domain estimates and regularity results directly for the equations of acoustics. For wave-structure interaction, the presence of time derivatives in the coupling conditions adds difficulty to the problem. Nevertheless, we are currently working on the extension of these analysis techniques to such problems, in particular to the acoustic/elastic and acoustic/piezoelectric interaction [2]. Discretization of Boundary Integral Equations with reduced quadrature (deltabem) For smooth parametrizable scatterers in 2D we devised and implemented an inexpensive third-order method based on an improvement of ǫ-quadrature, a technique dating back to the work of Saranen and Vainikko in the 90 s [15]. The method, which we refer to as deltabem, is built on the idea of collocating the equation on equispaced points in parametric space and then approximating the integral operators using a one-point quadrature rule in carefully shifted symmetric grids as shown in Figure 1. The technique can be understood as a non-conforming Petrov-Galerkin discretization using piecewise constants to approximate H 1/2 (Γ) and Dirac deltas for H 1/2 (Γ). The shape of the test functions is controlled by the shifting parameter ǫ and is optimized for third order convergence. In collaboration with Francisco Javier Sayas, Matthew Hassell, Sijiang Lu and Tianyu Qiu, we implemented this idea for the operators associated with Laplace, Helmholtz and Navier-Lamé equations. The code is written in Matlab and is freely available in Generalizing the ideas behind deltabem to 3D scattering is the subject of ongoing research. Trial Test i 2 i 1 i i+1 i+2 i 2 i 1 i i+1 i+2 i 2 i 1 i i+1 i+2 i 2 i 1 i i+1 i+2 Figure 1: Left: The kite-shaped figure is a cartoon of the staggered grids on the boundary of a scatterer. Open dots represent collocation points and colored crosses represent the symmetrically shifted quadrature points. Center: the shapes for the optimal trial and test functions used as approximations of the trace space (bottom row) and its dual (top row). Right: A possible generalization for 3D scatterers; the equation is collocated at the circles while the red crosses are the location of the optimized quadrature points. Elastic wave propagation and deltabem The analysis of deltabem was recently developed for the operators related to the Helmholtz equation resulting in a one-parameter family of test functions[3]. Nevertheless, for the integral operators associated to the Navier-Lamé resolvent equation (µ( u+ u )+λ ui ) = s 2 u, the extension we developed in [4] was not straightforward due to the two difficulties listed below. 2
3 First of all, in the elastic case, the double layer operator and its adjoint contain a strongly singular operator (a perturbation of the periodic Hilbert transform) that makes the operators of the second kind much more difficult to handle. In collaboration with Francisco Javier Sayas and Víctor Domínguez [4], we have provided mathematical and numerical evidence that a particular symmetric average of the shiftings ǫ = ±1/6,±5/6 yields third order convergence and is in fact the only choice with this property. The second difficulty was the lack of an appropriate regularization of the hypersingular operator compatible with our method for discretization. We implemented a regularization due to Frangi and Novati [7] which results in an integro-differential operator that is consistent with the deltabem discretization. This provided a complete discretization of the operators of the Calderón Projector that can be used to solve well-posed frequency domain problems or can be combined with Convolution Quadrature for timedomain problems. Figure 2 shows an example of a simulation obtained using this method. Future lines of research in this direction include the analysis of a deltabem discretization of the Calderón projector associated to Stokes s resolvent equations to discretize integral formulations like the one outlined in [1]. Figure 2: A plane pressure wave interacts with a rigid body, the total wave satisfies homogeneous Dirichlet boundary conditions. After the wave is reflected it decomposes into pressure and shear components according to the angle of incidence with respect to the normal. The simulation uses deltabem for space discretization and Trapezoidal Rule Convolution Quadrature for time evolution. Boundary Integral Equation and coupled schemes for wave-structure interaction The study of the scattering of acoustic waves by obstacles with various physical properties and the corresponding reaction of the scatterer is of great importance in the applied sciences. Applications include structural analysis, noise reduction, seismic imaging, and many others. The analysis and simulation of this processes in the time domain can greatly enhance the understanding of the underlying physics by providing additional information related to the transient stages of the problem. In [10, 14, 9] we have addressed the time-domain scattering of acoustic waves by bounded linearly elastic solids, piezoelectric solids and thermoelastic solids. These problems require the coupling of PDE s of hyperbolic, eliptic and parabolic type with non-standard transmission conditions. Table 1 shows the concrete system of PDE s modeling each interaction. In the acoustic/elastic case we have developed two well posed formulations; one based solely on Boundary Integral Equations for homogeneous scatterers, and one coupled Boundary Integral Equation/Variational formulation for general non-homogeneous anisotropic scatterers. For the piezoelectric and thermoelastic cases we have developed symmetric coupling schemes that allow for boundary element treatment of the acoustic wave and a finite element solution for the elastic, electric and thermal unknowns in the solid. The analysis is done in the Laplace domain assuming Lipschitz scatterers and sufficiently smooth incident waves. Well-posedness, at the continuous and semi-discrete-in-space levels, is established by equivalent transmission and variational problems from which stability bounds in terms of the Laplace parameter are derived. Once the systems are discretized in time with Convolution Quadrature, the stability bounds can be translated into time-domain estimates with explicit time dependence. In the short term my goal in this area is to develop integro-differential and boundary-integral formulations of this type for the fluid-structure interaction of elastic bodies of different sorts immersed in 3
4 a Stokes s flow. Such well posed formulations would allow for the simulation of a wide variety of fluidsolid interactions and, in combination with the computational and analytical tools like deltabem and CQ, would provide a computational toolbox for inexpensive time-domain simulation of a wide array of multiphysics problems. Acoustic/Elastic Acoustic/Piezoelectric Acoustic/Thermoelastic v =c 2 v tt in Ω + v =c 2 v tt in Ω + v =c 2 σ =ρ Σ u tt in Ω σ ζ θ =ρ Σ u tt in Ω v tt in Ω + D =0 in Ω κ θ =(θ +κη u) t in Ω σ =ρ Σ u tt in Ω u t ν = ν(v +v inc u t ν = ν(v +v inc ) on Γ (σ ζθi)ν = ρ f (v +v inc ) tν on Γ ) on Γ σν = ρ f (v +v inc σν = ρ f (v +v inc ) tν on Γ u t ν = ν(v +v inc ) on Γ ) tν on Γ D ν =η d on Γ N νθ =η d on Γ N ψ =µ d on Γ D θ =µ d on Γ D Table 1: Wave-structure Interaction: The equations are posed for t 0 and all the unknowns satisfy homogeneous initial conditions. Here Ω + denotes the unbounded exterior of the scatterer, Ω the interior of the scatterer and Γ the interface which is split into Dirichlet and Neumann parts Γ = Γ D Γ N in the piezoelectric and thermoelastic problems. The density of the elastic body is represented by ρ Σ and the speed of sound in the acoustic domain is c, also v is the unknown scattered acoustic wave, u is the elastic displacement wave and ψ is the electric potential. The stress tensor σ is a function of the displacement u. In piezoelectric solids both σ and the electric displacement vector D depend on a combination of u and ψ, thus coupling the system. The temperature variation with respect to equilibrium is denoted by θ, the thermal diffusivity κ and the thermal expansion coefficient ζ are considered constant. Finite volume methods for non-linear conservation laws In the past I have worked on the numerical treatment of first order non-linear conservation laws. The non-linearity of the problem implies that singularities and discontinuities can develop in the solutions for finite times even for smooth initial data. The classic Godunov-type methods solve this problem by discretizing the initial conditions (for instance averaging on a control volume or cell) and approximating them with piecewise constants or polynomials. The equations are then solved exactly for this approximate initial conditions which requires solving or approximating so-called Riemann problems. Methods of this family are particularly well suited to capture developing shocks in the solution, avoid Gibbs oscillations, and preserve the conservative nature of the continuous problem. For my masters thesis [13] I implemented a FORTRAN code for the solution of Euler s equations of gas dynamics in vacuum t ρ ρu ρv ρw E + x ρu ρu 2 +p ρvu ρwu u(e +p) + y ρv ρuv ρv 2 +p ρwv v(e +p) + z ρw ρuw ρvw ρw 2 +p w(e +p) = 0. Here u,v,w are the components of the velocity, ρ is the gas density, E is the energy and p is the pressure. The 3D implementation was based on dimensional splitting and required the handling of the interfaces with vacuum. It also explored different choices of flux limiters to increase the order of the approximation in the regions far from the shock profiles. Figure 3 shows snapshots of a simulation of a gas sphere colliding with a wall of gas. Recent work has been done [5, 6] on the coupling of Boundary Elements with Finite Volumes for time-independent problems. In the near future I plan to undertake the analysis and implementation of the coupling of these two methods in the time domain. Such couplings would potentially be applied to the time-domain simulation of phenomena such as crack propagation, where the behavior of the fluid in the cracks is traditionally computed on the time-independent regime using finite volumes and the elastic effects are handled with Boundary Elements. 4
5 Figure 3: Finite volume simulation of a sphere of gas colliding with a wall of gas. The plot shows density iso-surfaces at four different times. The method used dimensional splitting and first order explicit time-stepping. It involved solving a nonlinear Riemann problem at the inteface of each cell for every time step and verifying a positivity condition to handle the generation of vacuum. References [1] C. Bacuta, M. E. Hassell, G. C. Hsiao, and F.-J. Sayas. Boundary integral solvers for an evolutionary exterior stokes problem. SIAM Journal on Numerical Analysis, 53(3): , [2] T. Brown, T. Sánchez-Vizuet, and F.-J. Sayas. Evolution of a semidiscrete system modeling the scattering of acoustic waves by a piezoelectric solid. Submitted, [3] V. Domínguez, S. Lu, and F.-J. Sayas. A Nyström flavored Calderón calculus of order three for two dimensional waves, time-harmonic and transient. Comput. Math. Appl., 67(1): , [4] V. Domínguez, T. Sánchez-Vizuet, and F.-J. Sayas. A fully discrete Calderón calculus for the two-dimensional elastic wave equation. Comput. Math. Appl., 69(7): , [5] C. Erath. Comparison of two couplings of the finite volume method and the boundary element method. In J. Fuhrmann, M. Ohlberger, and C. Rohde, editors, Finite Volumes for Complex Applications VII-Methods and Theoretical Aspects, volume 77 of Springer Proceedings in Mathematics & Statistics, pages Springer International Publishing, [6] C. Erath, G. Of, and F.-J. Sayas. A non-symmetric coupling of the finite volume method and the boundary element method. submitted, [7] A. Frangi. and G. Novati. Regularized symmetric Galerkin BIE formulations in the Laplace transform domain for 2D problems. Computational Mechanics, 22(1):50 60, [8] M. Hassell, T. Qiu, T. Sánchez-Vizuet, and F.-J. Sayas. A new and improved analysis of the time domain boundary integral operators for acoustics. To appear in the Journal of Integral Equations and Applications, [9] G. Hsiao, T. Sánchez-Vizuet, F.-J. Sayas, and R. Weinacht. A time-dependent wave-thermoelastic solid interaction. Submitted, [10] G. C. Hsiao, T. Sánchez-Vizuet, and F.-J. Sayas. Boundary and coupled boundary-finite element methods for transient wave-structure interaction. IMA J. Numer. Anal., 37(1): , [11] C. Lubich. Convolution quadrature and discretized operational calculus. I. Numer. Math., 52(2): , [12] C. Lubich. Convolution quadrature and discretized operational calculus. II. Numer. Math., 52(4): , [13] T. Sánchez-Vizuet. Numerical solution of Euler s equations of gas dynamics in 3D. Master s thesis, Universidad Nacional Autónoma de México. Instituto de Matemáticas, June (In Spanish). [14] T. Sánchez-Vizuet and F.-J. Sayas. Symmetric boundary-finite element discretization of time dependent acoustic scattering by elastic obstacles with piezoelectric behavior. Journal of Scientific Computing, 70(3): , [15] J. Saranen and G. Vainikko. Periodic Integral and Pseudodifferential Equations with Numerical Approximation. Springer Monographs in Mathematics. Springer-Verlag, Berlin,
Discontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationMath 126 Winter 2012: Rough lecture notes
Math 126 Winter 2012: Rough lecture notes Alex Barnett, Dartmouth College March 1, 2012 1 Introduction Numerical mathematics is at the intersection of analysis (devising and proving theorems), computation
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationMETHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS
METHODS FOR SOLVING MATHEMATICAL PHYSICS PROBLEMS V.I. Agoshkov, P.B. Dubovski, V.P. Shutyaev CAMBRIDGE INTERNATIONAL SCIENCE PUBLISHING Contents PREFACE 1. MAIN PROBLEMS OF MATHEMATICAL PHYSICS 1 Main
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
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 informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
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 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 informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
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 informationApplied Numerical Mathematics. High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces
Applied Numerical Mathematics 111 (2017) 64 91 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum High-order numerical schemes based on difference potentials
More informationA fully discrete Calderón Calculus for the two-dimensional elastic wave equation
A fully discrete Calderón Calculus for the two-dimensional elastic wave equation Víctor Domínguez, Tonatiuh Sánchez-Vizuet & Francisco Javier Sayas July 16, 18 arxiv:149.7653v1 [math.na] 6 Sep 14 Abstract
More informationSOME APPLICATIONS OF INTEGRAL EQUATIONS TO THE SOLUTION OF TRANSIENT PARTIAL DIFFERENTIAL EQUATIONS. Matthew E. Hassell
SOME APPLICATIONS OF INTEGRAL EQUATIONS TO THE SOLUTION OF TRANSIENT PARTIAL DIFFERENTIAL EQUATIONS by Matthew E. Hassell A dissertation submitted to the Faculty of the University of Delaware in partial
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationA fully discrete BEM-FEM scheme for transient acoustic waves
A fully discrete BEM-FEM scheme for transient acoustic waves Matthew E. Hassell & Francisco Javier Sayas Department of Mathematical Sciences, University of Delaware, USA {mhassell,fjsayas}@udel.edu arxiv:1601.08248v1
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationOn the Solution of the Elliptic Interface Problems by Difference Potentials Method
On the Solution of the Elliptic Interface Problems by Difference Potentials Method Yekaterina Epshteyn and Michael Medvinsky Abstract Designing numerical methods with high-order accuracy for problems in
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
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 informationNEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEABED
Georgian Mathematical Journal Volume 4 2007, Number, 09 22 NEAR FIELD REPRESENTATIONS OF THE ACOUSTIC GREEN S FUNCTION IN A SHALLOW OCEAN WITH FLUID-LIKE SEAED ROERT GILERT AND MIAO-JUNG OU Abstract. In
More informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
More informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
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 informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationR. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant
R. Courant and D. Hilbert METHODS OF MATHEMATICAL PHYSICS Volume II Partial Differential Equations by R. Courant CONTENTS I. Introductory Remarks S1. General Information about the Variety of Solutions.
More informationNumerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry
Numerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry V. S. Ajaev and J. Tausch 2 Southern Methodist University ajaev@smu.edu 2 Southern Methodist University tausch@smu.edu
More informationDetermination of thin elastic inclusions from boundary measurements.
Determination of thin elastic inclusions from boundary measurements. Elena Beretta in collaboration with E. Francini, S. Vessella, E. Kim and J. Lee September 7, 2010 E. Beretta (Università di Roma La
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationFactorization method in inverse
Title: Name: Affil./Addr.: Factorization method in inverse scattering Armin Lechleiter University of Bremen Zentrum für Technomathematik Bibliothekstr. 1 28359 Bremen Germany Phone: +49 (421) 218-63891
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationPartial Differential Equations
Partial Differential Equations Analytical Solution Techniques J. Kevorkian University of Washington Wadsworth & Brooks/Cole Advanced Books & Software Pacific Grove, California C H A P T E R 1 The Diffusion
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationFrequency functions, monotonicity formulas, and the thin obstacle problem
Frequency functions, monotonicity formulas, and the thin obstacle problem IMA - University of Minnesota March 4, 2013 Thank you for the invitation! In this talk we will present an overview of the parabolic
More informationFast iterative BEM for high-frequency scattering
Fast iterative BEM for high-frequency scattering elastodynamic problems In 3D. Marion Darbas LAMFA UMR CNRS 7352, Amiens, France S. Chaillat (POEMS, CNRS-INRIA-ENSTA) and F. Le Louër (LMAC, UTC) Conference
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationOn uniqueness in the inverse conductivity problem with local data
On uniqueness in the inverse conductivity problem with local data Victor Isakov June 21, 2006 1 Introduction The inverse condictivity problem with many boundary measurements consists of recovery of conductivity
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 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 informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
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 informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
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 informationGeneralized Functions Theory and Technique Second Edition
Ram P. Kanwal Generalized Functions Theory and Technique Second Edition Birkhauser Boston Basel Berlin Contents Preface to the Second Edition x Chapter 1. The Dirac Delta Function and Delta Sequences 1
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
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 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 informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
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 informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationCausal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases
Causal Dissipation for the Relativistic Fluid Dynamics of Ideal Gases Heinrich Freistühler and Blake Temple Proceedings of the Royal Society-A May 2017 Culmination of a 15 year project: In this we propose:
More informationApplications of Integral Equation Methods to a Class of Fundamental Problems in Mechanics and Mathematical Physics
Outline Applications of Integral Equation Methods to a Class of Fundamental Problems in Mechanics and Mathematical Physics George C. Hsiao Department of Mathematical Sciences University of Delaware Newark,
More informationThe Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations
The Embedded Boundary Integral Method (EBI) for the Incompressible Navier-Stokes equations George Biros Lexing Ying and Denis Zorin Courant Institute of Mathematical Sciences, New York University, NY 10012,
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationASYMPTOTIC THEORY FOR WEAKLY NON-LINEAR WAVE EQUATIONS IN SEMI-INFINITE DOMAINS
Electronic Journal of Differential Equations, Vol. 004(004), No. 07, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) ASYMPTOTIC
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
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 recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationDRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS
DRIFT-DIFFUSION SYSTEMS: VARIATIONAL PRINCIPLES AND FIXED POINT MAPS FOR STEADY STATE SEMICONDUCTOR MODELS Joseph W. Jerome Department of Mathematics Northwestern University Evanston, IL 60208 Abstract
More informationTHE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS
THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. BY DR. LOTHAR COLLATZ
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More information