INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
|
|
- Kenneth Fitzgerald
- 6 years ago
- Views:
Transcription
1 INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station, TX Supervisor: Dr. Andrzej J. Safjan June 20, 2003 Support of this work by NSF grant # and by ONR grant #N is greatfully acknowledged. Typeset by AMS-TEX
2 Motivating Problem: Scattering of Sonar Waves from a Submarine Hull. TRUNCATION BOUNDARY ASPECT RATIO = 7: Fig. 3: Surrounding a Submarine with an Infinite Element Prolate Spheroid. 2
3 Objectives. To obtain highly-accurate computational solutions to Helmholtz equation problems posed on unbounded domains using infinite elements. 2. To formulate infinite elements which permit a large number of degrees of freedom to be utilized in the radial (infinite) direction. 3. To formulate infinite elements which have numerically well-conditioned stiffness matrices. 4. To formulate infinite elements with h-p-r adaptivity in the radial direction. 5. To formulate infinite elements with the ability to surround boundaries of various aspect ratios. TRUNCATION BOUNDARY 2b 2a Fig. 2: Elliptical Boundary with Aspect Ratio=b/a. 3
4 Three-Dimensional Helmholtz Equation We are interested in solving the following three-dimensional Helmholtz problem posed on the unbounded domain Ω with boundary Ω: u k 2 u =0, u n = g, lim r r Ω x Ω ) ( u r + iku =0 where r is the radial distance from the origin. TRUNCATION BOUNDARY Ω Ω x r ELASTIC SOLID j 0 i n ACOUSTIC FLUID Fig. : Unbounded Domain and Boundary. 4
5 3-D Spherical Multipole Expansion Recall that in spherical coordinates x = r sin θ(cos φi + sin φj)+r cos θk where r 0, 0 θ π, 0 φ<2π, the outward traveling wave solutions of the 3-D Helmholtz equation admit the following eigenfunction expansion u( r, θ, φ )= n n=0 m=0 c nm h (2) n (kr)pn m (cos θ)e imφ where h (2) n (x) are the spherical Hankel functions of the second kind and Pn m (x) are the associated Legendre functions of the first kind. The spherical Hankel functions of the second kind are given by h (2) n (x) = i ( ) x p n e ix, for n =0,, 2,... x where p n (x) are polynomials which satisfy the recurrence relation p n+ (x) =(2n +)xp n (x) p n (x), for n =, 2,... and p 0 (x) =, p (x) =x + i. 5
6 Thus follows the classical multipole expansion theroem by Atkinson and Wilcox [Atkinson, 23], [Wilcox, 24,25]. Theorem. Let u(x) be a radiation function for the region exterior to a sphere x x 0 = c, and let ( r, θ, φ ) be the spherical coordinates for x relative to an origin at x 0. Then u( r, θ, φ )=e ikr n= G n ( θ, φ ) r n where the series converges for r>cand converges absolutely and uniformly in r, θ and φ in any region r c + ɛ>c. The series may be differentiated term by term with respect to r, θ and φ any number of times and the resulting series all converge absolutely and uniformly. 3-D Spheroidal Multipole Expansions [Holford, ] In [] Holford extended the multipole expansion theorem u( ρ, ˆθ, ˆφ )=e ikρ n= G n ( ˆθ, ˆφ ) ρ n of Atkinson and Wilcox to Prolate Spheroidal coordinates x = ρ 2 f 2 sin ˆθ(cos ˆφi + sin ˆφj)+ρcos ˆθk and Oblate Spheroidal coordinates x = ρ sin ˆθ(cos ˆφi + sin ˆφj)+ ρ 2 f 2 cos ˆθk where f 0, ρ f, 0 ˆθ π, 0 ˆφ <2π. 6
7 Weak Formulation of the 3-D Helmholtz Problem Weighted and Conjugated Weak Formulation. where a: U V C π 2π a( u, v )= 0 0 ρ ρ 2 Find u U such that a( u, v )=l(v), v V is the sesquilinear form and l : V C l(v) = ( u v k 2 uv 2 ) u gρρ ρ ρ v J dρd ˆφdˆθ Ω ρ 2 gvds is the antilinear functional. Here U = V = H w(ω) where H w(ω) is the Sobolev space endowed with the norm u H w (Ω) = π 0 2π 0 ρ ρ 2 ( u 2 + u 2 ) J dρd ˆφdˆθ. On the boundary Ω weconsider the function space L 2 ( Ω) endowed with the norm u L 2 ( Ω) = Ω u 2 ds. 7
8 For Prolate Spheroidal coordinates J = g ρρ = (ρ 2 f 2 ρ 2 f ρ 2 f 2 ( + cos 2ˆθ), ) ( + cos 2ˆθ) sin ˆθ >0. Infinite Element Approximation. where Find u h U h such that a( u h,v h )=l(v h ), v h V h U h U, V h V dim U h <, dim V h < 8
9 3-D Multipole Based Infinite Elements [Burnett, 7 0] In [7 0] Burnett developed 3-D infinite elements based on the Prolate and Oblate Spheroidal multipole expansions. The ϕ j ( ρ, ξ, η ), infinite element basis functions are the tensor product of the multipole functions ψ n (ρ) = e ikρ ρ n, for n =, 2,..., ρ ( ρ, ) in the radial direction with χ i ( ξ,η ), Lagrange polynomial shape functions in the angular directions where the transformation from local coordinates to global angular coordinates is given by ˆθ( ξ,η )= i ˆθ i χ i ( ξ,η ), ˆφ( ξ,η )= The corresponding approximate solution is i ˆφ i χ i ( ξ,η ). u h ( ρ, ξ, η )= j u j ϕ j ( ρ, ξ, η ). 9
10 Model Problem. Using the 3-D Global Multipole infinite element we compute the solution to the 3-D Helmholtz equation whose analytic solution is u( r, θ, φ )=h (2) 0 (kr)p 0(cos θ). The domain is taken to be the region outside the sphere Ω={ ( r, θ, φ ) r> }, Ω={ ( r, θ, φ ) r =}. This model problem is solved for k =. INFINITE ELEMENT MESH Fig. 4: Spherical Multipole Infinite Element Mesh. 0
11 CONDITION NUMBER.E+2.E+9.E+7 MULTIPOLE INFINITE ELEMENT.E+5.E+3.E+.E+9.E+7 NE= Number of Infinite Elements in mesh=00 NORMALIZED ERROR e Ω = u u h H w (Ω) u H w (Ω).E+5.E-.E-2.E-3 MULTIPOLE INFINITE ELEMENT NUMBER OF RADIAL DEGREES OF FREEDOM.E-4.E-5 Fig. 5: ILL-Conditioning of the Spherical Multipole Infinite Element.
12 Remedies for the Numerical Ill-Conditioning. A preconditioning method based on a Gram-Schmidt-like process applied to the global multipole functions. 2. Replace the global multipole functions with piecewise multipole functions which have compact support. Preconditioning Method Given the multipole functions ψ n (ρ) = e ikρ ρ n, for n =, 2,..., ρ ( ρ, ) and the sesquilinear form a: V V C associated with the weak formulation of the Helmholtz problem, consider the Gram matrix Γ C n n defined by a( ψ,ψ ) a( ψ,ψ 2 ) a( ψ,ψ n ) a( ψ 2,ψ ) a( ψ 2,ψ 2 ) a( ψ 2,ψ n ) Γ = a( ψ n,ψ ) a( ψ n,ψ 2 ) a( ψ n,ψ n ) Let ψ C n be the vector ψ(ρ) = 2 ψ (ρ) ψ 2 (ρ). ψ n (ρ)
13 and consider a linear transformation of ψ induced by an n n invertible matrix X ˆψ(ρ) =Xψ(ρ). The Gram matrix Γ is transformed into the Gram matrix ˆΓ = XΓX H. Consider now the following decomposition of Γ Γ = LDU = LŨ where D C n n is a diagonal matrix, L C n n is a lower triangular matrix with unit diagonal elements, and U C n n is an upper triangular matrix with unit diagonal elements. Moreover, L = LD 2, Ũ = D 2 U, and D 2 is the positive square root of D computed using the following formula for the positive square root of a complex number: { 2 a + ib = 2 ( a2 + b 2 + a + i a2 + b 2 a), b ( a2 + b 2 + a i. a2 + b 2 a), b < 0 The Gram-Schmidt transformation consists of taking the transformation matrix X to be the lower triangular matrix X = L = D 2 L. Consequently, the transformed Gram matrix ˆΓ = XΓX H H = Ũ L is upper triangular with diagonal elements having unit complex modulus. In other words, the Gram-Schmidt transformation yields transformed fields ˆψ which satisfy the following conditions { a( ˆψ i, ˆψ 0 i<j j ) = i = j. 3
14 The preconditioned infinite element basis functions ˆϕ j ( ρ, ξ, η ) are the tensor product of the Gram-Schmidt transformed multipole functions ˆψ n (ρ) inthe radial direction with the Lagrange polynomial shape functions χ i ( ξ,η )inthe angular directions. The approximate solution is u h ( ρ, ξ, η )= j û j ˆϕ j ( ρ, ξ, η ). 4
15 infinity Ω Ω j 0 i e α3 e α2 e Ω e α Fig. 6: Infinite Element Mesh with Typical Element e. 5
16 CONDITION NUMBER NE=70.E+7 NE=60 NE=50 NE=40.E+6 NE=30 NE=20 NE=0 NORMALIZED ERROR.E+5 e Ω = u u h u.e NE= Number of Infinite Elements in mesh NUMBER OF MULTIPOLE TERMS NE=0.E-5 NE=20 NE=30 NE=40 NE=50 NE=60 NE=70 Fig. 7: Assembled Stiffness Matrix Condition Numbers and Convergence Curves of the Preconditioned Multipole Infinite Element. 6
17 3-D Piecewise Multipole Infinite Element In formulating the 3-D Piecewise Multipole infinite element we replace the global multipole functions ψ n (ρ) = e ikρ ρ n, for n =, 2,..., ρ ( ρ, ) with Lagrange interpolation functions σ p (ρ) inthe variable ρ. Again, the infinite element basis functions ϕ j ( ρ, ξ, η ) are the tensor product of the piecewise multipole functions σ p (ρ) in the radial direction with Lagrange polynomial shape functions χ i ( ξ,η )inthe angular directions. -D Meshing of the Unbounded Radial Domain. We consider a -D meshing of the unbounded radial domain ( ρ, )into element subdomains. Let n e =number of elements in the mesh n =number of nodes per element n n = total number of nodes in the mesh Each element for e =,...,n e isbounded and has n nodes while the element e = n e is unbounded and has n nodes. The total number of radial nodes for this mesh is n n =(n )n e. 7
18 ρ ρ ρ 2 ρ n infinity e e Τ ζ - Fig. 8: -D Meshing of the Unbounded Radial Domain. 8
19 Radial Shape Functions For e =,...,n e the radial shape functions are n ( ) σ p (ρ) = q= q p n ( ρ ρ q )e ikρ, for p =, 2,...,n q= q p ρ p ρ q where ρ,ρ 2,...,ρ n are the nodal coordinates. To perform numerical integrations we introduce the transformation e T :(, ) ( ρ,ρ n ) e T (ζ) = 2 ( ρ n ρ ) ζ + 2 ( ) ρ n + ρ with Jacobian e J(ζ) = ( 2 ( 2 ( ) ρ n ρ ρ n ρ ) ζ + 2 ( )) 2. ρ n + ρ 9
20 For e = n e we take the limit as ρ n. Therefore the shape functions become σ p (ρ) = ρ n q= q p n ( ( ρ ρ q ) ) e ikρ, for p =, 2,...,n. ρ p q= q p ρ p ρ q The transformation e T :(, ) ( ρ,ρ n )becomes and the Jacobian becomes e T (ζ) = 2ρ ζ e J(ζ) = 2ρ ( ζ) 2. 20
21 Distribution of the Radial Nodes in 3-D The distribution of the radial nodes for the 3-D Piecewise Multipole infinite element is based on the transformation T :(, 0) ( ρ, ) T (z) = ρ z.. Uniform Distribution of the Radial Nodes. A uniform partition of the interval (, 0)into n n subintervals yields the radial nodal coordinates ρ q = T ( + q ), for q =, 2,...,n n. n n 2. Non-Uniform Distribution of the Radial Nodes. Using a non-uniform arithmetic partition of the interval (, 0) into n n subintervals yields the radial nodal coordinates ρ q = T ( + ) q( q), for q =, 2,...,n n. n n (n n +) 2
22 TYPICAL INFINITE ELEMENT Fig. 29: Elliptical Mesh of Piecewise Multipole Infinite Elements. 22
23 Model Problem 3. 3-D Numerical Examples Using the Multipole and Piecewise Multipole infinite elements we compute the solution to the 3-D Helmholtz equation whose analytic solution is u( ρ, ˆθ, ˆφ )=h (2) (kr( ρ, ˆθ, ˆφ ))P (cos θ( ρ, ˆθ, ˆφ )) where the transformation from Prolate Spheroidal to Spherical coordinates is given by r = ρ 2 f 2 2 ( cos 2ˆθ), cos θ = ρ cos ˆθ. ρ 2 f 2 2 ( cos 2ˆθ) The domain is taken to be the region outside the spheroid defined by Ω={ ( ρ, ˆθ, ˆφ ) ρ> }, Ω={ ( ρ, ˆθ, ˆφ ) ρ =}. This model problem is solved for k =and aspect ratio of 7:. 23
24 ASPECT RATIO = 7: INFINITE ELEMENT MESH Fig. 2: Prolate Piecewise Multipole Infinite Element Mesh. 24
25 NORMALIZED BOUNDARY ERROR e Ω = u u h L2 ( Ω) u L2 ( Ω).E- GLOBAL MULTIPOLE INFINITE ELEMENT NUMBER OF RADIAL DEGREES OF FREEDOM.E-2 P=.E-3 P=2.E-4.E-5 P= Degree of Radial Shape Functions. NE= Number of Infinite Elements in mesh=00 P=3.E-6 ASPECT RATIO = 7: P=4 P=5 Fig. 5: Comparison of the Boundary Convergence Curves of the Global Multipole and Piecewise Multipole Infinite Elements. Weighted and Conjugated Formulation. Aspect Ratio=7:. Non-uniform Radial Mesh. 25
26 CONDITION NUMBER.E+2.E+9.E+7 GLOBAL MULTIPOLE INFINITE ELEMENT.E+5.E+3.E+ P= Degree of Radial Shape Functions. NE= Number of Infinite Elements in mesh=00 ASPECT RATIO = 7:.E+9.E+7 P= P=2 P=3 P=4 P=5 NORMALIZED ERROR e Ω = u u h H w (Ω) u H w (Ω).E+5.E- GLOBAL MULTIPOLE INFINITE ELEMENT NUMBER OF RADIAL DEGREES OF FREEDOM P=.E-2 P=2.E-3 P=3.E-4 P=4 P=5 Fig. 6: Comparison of the Assembled Stiffness Matrix Condition Numbers and Domain Convergence Curves of the Global Multipole and Piecewise Multipole Infinite Elements. Weighted and Conjugated Formulation. Aspect Ratio=7:. Non-uniform Radial Mesh. 26
27 Unweighted and Unconjugated Weak Formulation. An Unweighted and Unconjugated weak formulation of the 3-D Helmholtz problem is a( u, v )=l(v), v V, u U where a: U V C π 2π ρ ( a( u, v )= lim u v k 2 uv ) J dρd ˆφdˆθ ρ 0 0 ρ lim ρ (ρ2 f 2 ) π 2π is the bilinear form and l : V C l(v) = gv ds 0 Ω 0 u v sin ˆθdˆφdˆθ ρ is the linear functional. Here U = V = H w(ω) where H w(ω) is the Sobolev space endowed with the norm u H w (Ω) = π 0 2π 0 ρ ρ 2 ( u 2 + u 2 ) J dρd ˆφdˆθ. On the boundary Ω weconsider the function space L 2 ( Ω) endowed with the norm u L 2 ( Ω) = Ω u 2 ds. 27
28 NORMALIZED BOUNDARY ERROR e Ω = u u h L2 ( Ω) u L2 ( Ω).E-.E-2 GLOBAL MULTIPOLE INFINITE ELEMENT P=2.E-3 P=3.E-4 P=4.E-5 P= Degree of Radial Shape Functions. P=5.E-6.E-7 NE= Number of Infinite Elements in mesh=00 ASPECT RATIO = 7:.E NUMBER OF RADIAL DEGREES OF FREEDOM Fig. 23: Comparison of the Boundary Convergence Curves of the Global Multipole and Piecewise Multipole Infinite Elements. Unweighted and Unconjugated Formulation. Aspect Ratio=7:. Non-uniform Radial Mesh. 28
29 NORMALIZED ERROR e Ω = u u h H w (Ω) u H w (Ω) GLOBAL MULTIPOLE INFINITE ELEMENT.E+4.E+3.E+2 P= Degree of Radial Shape Functions. NE= Number of Infinite Elements in mesh=00 ASPECT RATIO = 7:.E+.E-.E P=4 P=3 P=2 P=5 NUMBER OF RADIAL DEGREES OF FREEDOM.E-3.E-4.E-5 Fig. 24: Comparison of the Domain Convergence Curves of the Global Multipole and Piecewise Multipole Infinite Elements. Unweighted and Unconjugated Formulation. Aspect Ratio=7:. Non-uniform Radial Mesh. 29
30 CONDITION NUMBER.E+2.E+9 GLOBAL MULTIPOLE INFINITE ELEMENT P= Degree of Radial Shape Functions..E+7 NE= Number of Infinite Elements in mesh=00.e+5 ASPECT RATIO = 7:.E+3 P=4.E+.E+9 P=3 P=2 P=5.E+7.E NUMBER OF RADIAL DEGREES OF FREEDOM Fig. 25: Comparison of the Assembled Stiffness Matrix Condition Numbers of the Global Multipole and Piecewise Multipole Infinite Elements. Unweighted and Unconjugated Formulation. Aspect Ratio=7:. Non-uniform Radial Mesh. 30
31 Effect of Frequency Variations Using the 2-D Piecewise Multipole infinite element we compute the solution to model problem 4 for k =,k =0and k = 00. CONDITION NUMBER.E+5 K= P= Degree of Infinite Element Shape Functions=2.E+4 NE= Number of Infinite Elements in mesh=0 K=0.E+3 K=00 NORMALIZED ERROR e Ω = u u h H w (Ω) u H w (Ω).E NUMBER OF RADIAL DEGREES OF FREEDOM K=.E-5 K=0 K=00 Fig. 32: Assembled Stiffness Matrix Condition Numbers and Convergence Curves for Various Wave Numbers. 3
32 Conclusions. A new piecewise multipole infinite element formulation has been developed. By replacing the global multipole functions with piecewise multipole functions the numerical conditioning of the 2-D and 3-D infinite element global stiffness matrices is greatly improved, enabling an effective use of more than 00 degrees of freedom in the radial direction. The 2-D and 3-D piecewise multipole infinite elements performed well in solving the model problems considered here. 2. The piecewise multipole formulation seems to be insensitive to the value of ka for the same mesh. For ka =, ka =0and ka = 00 the condition number and relative error essentially stay the same. 3. The piecewise multipole formulation gives convergence on the domain Ω and on the boundary Ω when implemented with the UWUC weak formulation of the 3-D Helmholtz equation problem. This contrasts with the global multipole formulation which diverges on Ω when implemented with the UWUC weak formulation. 32
33 Publications Complete details of this research can be found in: () A. Safjan and M. Newman, The Ill-Conditioning of Infinite Element Stiffness Matrices, Computers and Mathematics with Applications, Volume 4, Numbers 0-, May-June 200, pp (2) A. Safjan and M. Newman, On Two-Dimensional Infinite Elements Utilizing Basis Functions with Compact Support, Computers Methods in Applied Mechanics and Engineering, Volume 90, Number 48, September 200, pp (3) A. Safjan and M. Newman, On Three-Dimensional Infinite Elements Utilizing Basis Functions with Compact Support, Computers and Mathematics with Applications, Volume 43, Numbers 8-9, April-May 2002, pp Acknowledgements The support of the National Science Foundation (NSF) grant # and the Office of Naval Research (ONR) grant #N is greatly appreciated. 33
University of Kentucky, Lexington, Kentucky Bachelor of Science, August 1983 to May 1990 Major: Mechanical Engineering
Michael Newman Department of Aerospace Engineering 734 H.R. Bright Building College Station, TX 77843-3141 Home: (979)268-8335 Work: (979)845-0750 OBJECTIVE EDUCATION An entry level faculty position to
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 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 informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationStudia Geophysica et Geodaetica. Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids
Meshless BEM and overlapping Schwarz preconditioners for exterior problems on spheroids Journal: Studia Geophysica et Geodaetica Manuscript ID: SGEG-0-00 Manuscript Type: Original Article Date Submitted
More informationComputation of the scattering amplitude in the spheroidal coordinates
Computation of the scattering amplitude in the spheroidal coordinates Takuya MINE Kyoto Institute of Technology 12 October 2015 Lab Seminar at Kochi University of Technology Takuya MINE (KIT) Spheroidal
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
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 informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationSolving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method
Solving the Dirichlet-to-Neumann map on an oblate spheroid by a mesh-free method A. Costea a, Q.T. Le Gia b,, D. Pham c, E. P. Stephan a a Institut für Angewandte Mathematik and QUEST (Centre for Quantum
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationThe Atkinson Wilcox theorem in ellipsoidal geometry
J. Math. Anal. Appl. 74 00 88 845 www.academicpress.com The Atkinson Wilcox theorem in ellipsoidal geometry George Dassios Division of Applied Mathematics, Department of Chemical Engineering, University
More informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Discretization of Boundary Conditions Discretization of Boundary Conditions On
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationSome negative results on the use of Helmholtz integral equations for rough-surface scattering
In: Mathematical Methods in Scattering Theory and Biomedical Technology (ed. G. Dassios, D. I. Fotiadis, K. Kiriaki and C. V. Massalas), Pitman Research Notes in Mathematics 390, Addison Wesley Longman,
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationSound radiation from the open end of pipes and ducts in the presence of mean flow
Sound radiation from the open end of pipes and ducts in the presence of mean flow Ray Kirby (1), Wenbo Duan (2) (1) Centre for Audio, Acoustics and Vibration, University of Technology Sydney, Sydney, Australia
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
More informationNonlinear, Transient Conduction Heat Transfer Using A Discontinuous Galerkin Hierarchical Finite Element Method
Nonlinear, Transient Conduction Heat Transfer Using A Discontinuous Galerkin Hierarchical Finite Element Method by Jerome Charles Sanders B.S. in Physics, May 2002 The College of New Jersey A Thesis submitted
More informationScattering. March 20, 2016
Scattering March 0, 06 The scattering of waves of any kind, by a compact object, has applications on all scales, from the scattering of light from the early universe by intervening galaxies, to the scattering
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
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 informationAn adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
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 informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationA family of closed form expressions for the scalar field of strongly focused
Scalar field of non-paraxial Gaussian beams Z. Ulanowski and I. K. Ludlow Department of Physical Sciences University of Hertfordshire Hatfield Herts AL1 9AB UK. A family of closed form expressions for
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationChapter 2 Acoustical Background
Chapter 2 Acoustical Background Abstract The mathematical background for functions defined on the unit sphere was presented in Chap. 1. Spherical harmonics played an important role in presenting and manipulating
More informationHIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST SQUARES METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 4-5, Pages 604 626 c 2017 Institute for Scientific Computing and Information HIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationNull Cones to Infinity, Curvature Flux, and Bondi Mass
Null Cones to Infinity, Curvature Flux, and Bondi Mass Arick Shao (joint work with Spyros Alexakis) University of Toronto May 22, 2013 Arick Shao (University of Toronto) Null Cones to Infinity May 22,
More informationInverse Obstacle Scattering
, Göttingen AIP 2011, Pre-Conference Workshop Texas A&M University, May 2011 Scattering theory Scattering theory is concerned with the effects that obstacles and inhomogenities have on the propagation
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE HELMHOLTZ EQUATION
MATHEMATICS OF COMPUTATION Volume 68, Number 227, Pages 945 953 S 0025-5718(99)01064-9 Article electronically published on February 15, 1999 ON THE COUPLING OF BEM AND FEM FOR EXTERIOR PROBLEMS FOR THE
More informationH 2 -matrices with adaptive bases
1 H 2 -matrices with adaptive bases Steffen Börm MPI für Mathematik in den Naturwissenschaften Inselstraße 22 26, 04103 Leipzig http://www.mis.mpg.de/ Problem 2 Goal: Treat certain large dense matrices
More informationPART II : Least-Squares Approximation
PART II : Least-Squares Approximation Basic theory Let U be an inner product space. Let V be a subspace of U. For any g U, we look for a least-squares approximation of g in the subspace V min f V f g 2,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationEigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator.
PHYS208 spring 2008 Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. 07.02.2008 Adapted from the text Light - Atom Interaction PHYS261 autumn 2007 Go to list of topics
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
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 informationThe Coupling Method with the Natural Boundary Reduction on an Ellipse for Exterior Anisotropic Problems 1
Copyright 211 Tech Science Press CMES, vol.72, no.2, pp.13-113, 211 The Coupling Method with the Natural Boundary Reduction on an Ellipse for Exterior Anisotropic Problems 1 Quan Zheng 2, Jing Wang 2 and
More informationMath 575-Lecture Fundamental solutions for point force: Stokeslet, stresslet
Math 575-Lecture 2 Fundamental solutions for point force: tokeslet, stresslet Consider putting a point singular force with strength F in a tokes flow of infinite extent. We then need to solve P µ U = δ(x)f,
More informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
More informationNumerical Methods for Two Point Boundary Value Problems
Numerical Methods for Two Point Boundary Value Problems Graeme Fairweather and Ian Gladwell 1 Finite Difference Methods 1.1 Introduction Consider the second order linear two point boundary value problem
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationChapter 1. Condition numbers and local errors in the boundary element method
Chapter 1 Condition numbers and local errors in the boundary element method W. Dijkstra, G. Kakuba and R. M. M. Mattheij Eindhoven University of Technology, Department of Mathematics and Computing Science,
More informationPoint estimates for Green s matrix to boundary value problems for second order elliptic systems in a polyhedral cone
Maz ya, V. G., Roßmann, J.: Estimates for Green s matrix 1 ZAMM Z. angew. Math. Mech. 00 2004 0, 1 30 Maz ya, V. G.; Roßmann, J. Point estimates for Green s matrix to boundary value problems for second
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
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 informationThe Discontinuous Galerkin Method for Hyperbolic Problems
Chapter 2 The Discontinuous Galerkin Method for Hyperbolic Problems In this chapter we shall specify the types of problems we consider, introduce most of our notation, and recall some theory on the DG
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationFINITE ELEMENT SOLUTION OF SCATTERING IN COUPLED FLUID-SOLID SYSTEMS by Mirela O. Popa B.S., Harvey Mudd College, Claremont CA, 1995 M.S.
FINITE ELEMENT SOLUTION OF SCATTERING IN COUPLED FLUID-SOLID SYSTEMS by Mirela O. Popa B.S., Harvey Mudd College, Claremont CA, 995 M.S., University of Colorado at Denver, 997 A thesis submitted to the
More informationFast domain decomposition algorithm for discretizations of 3-d elliptic equations by spectral elements
www.oeaw.ac.at Fast domain decomposition algorithm for discretizations of 3-d elliptic equations by spectral elements V. Korneev, A. Rytov RICAM-Report 2006-21 www.ricam.oeaw.ac.at Fast domain decomposition
More informationOn spherical-wave scattering by a spherical scatterer and related near-field inverse problems
IMA Journal of Applied Mathematics (2001) 66, 539 549 On spherical-wave scattering by a spherical scatterer and related near-field inverse problems C. ATHANASIADIS Department of Mathematics, University
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationUniversity of Illinois at Chicago Department of Physics
University of Illinois at Chicago Department of Physics Electromagnetism Qualifying Examination January 4, 2017 9.00 am - 12.00 pm Full credit can be achieved from completely correct answers to 4 questions.
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
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 informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More information