hp-adaptive Finite Elements in the C++ library Concepts
|
|
- Patrick Branden Austin
- 6 years ago
- Views:
Transcription
1 hp-adaptive Finite Elements in the C++ library Concepts Kersten Schmidt Group POEMS, INRIA, France, Seminar for Applied Mathematics, ETH Zurich, Switzerland 8th July 9, NanoOptics Workshop, Zürich
2 Outline Space discretisation by h-, p- and hp-fem Maxwell s equations Variational formulation Linear system of equations Mesh refinement or increasing the (polynomial) order Resolving material interfaces essential p-fem gives exponential convergence for smooth interfaces hp-fem = suitable combination of h- and p-refinement exponential convergence for interface corners Applications Implementation in Concepts Overview Basis functions with local support Hierarchical basis based on Jacobi polynomials on reference cell Numerical quadrature on reference cell for curved physical cells General Algorithm for Matrix Assembling FEM, DG-FEM and BEM K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
3 Maxwell s equations Electromagnetic wave phenomena are described by Maxwell source problem in time domain permeability µ(x),permittivity ε(x, t), incident E inc (x, t) and initial field E init (x) µ (x) E(x, t) + ε(x, t) t E(x, t) = E inc (x, t) E(x, ) = E init (x) Looking for E(x, t) Maxwell source problem in frequency domain permeability µ(x),permittivity ε(x, ω), incident E inc (x), frequency ω µ (x) E(x) ε(x, ω)ω E(x) = E inc (x) Looking for E(x) Maxwell eigenvalue problem permeability µ(x),permittivity ε(x, ω) Looking for modes (E(x), ω) µ (x) E(x) ε(x, ω)ω E(x) = K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
4 Space discretisation Ansatz with finite functions to approximate exact solution E(x) E N (x) = NX i= c i E i (x) i.e. E E N < ɛ(n) N Finite Differences (FDM) Finite Elements (here p = ) Mesh refinement Expect algebraic convergence ɛ(n) N α Plane Wave Expansion (PWM) Finite Elements (fixed mesh) Increasing Order Expect exponential convergence ɛ(n) exp( βn α ) K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
5 Space discretisation of Maxwell s equation E.g. Maxwell eigenvalue problem µ (x) E(x) ε(x, ω) ω E(x) = ( ) Ansatz with finite functions to approximate exact solution E(x) E N (x) = X N i= c i E i (x) ( ) Finite Differences (FDM) Taylor expansion of E(x)around each grid point to derive stencils for the differential operators. Standard FDTD: Yee algorithm equivalence with variational form, where ansatz functions are Z-splines and test functions are Dirac distributions, J.T. Becerra Sagredo, SAM Report 3-, ETH Zurich. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
6 Space discretisation of Maxwell s equation E.g. Maxwell eigenvalue problem µ (x) E(x) ε(x, ω) ω E(x) = ( ) Ansatz with finite functions to approximate exact solution E(x) E N (x) = X N i= c i E i (x) ( ) Finite Differences (FDM) Taylor expansion of E(x)around each grid point to derive stencils for the differential operators. Standard FDTD: Yee algorithm equivalence with variational form, where ansatz functions are Z-splines and test functions are Dirac distributions, J.T. Becerra Sagredo, SAM Report 3-, ETH Zurich. Galerkin method for Finite Elements and Plane Wave Expansion multiply ( ) with test function F (x) integrate over the domain Ω, use integration by parts formula and insert boundary conditions insert ansatz ( ) and each basis function E i (x) as test function Z µ (x) E N (x) E i (x) ε(x, ω) ω E N (x) E i (x)dx = Ω Same procedure for Maxwell source problem in time or frequency domain. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
7 Space discretisation of Maxwell s equation E.g. Maxwell eigenvalue problem µ (x) E(x) ε(x, ω) ω E(x) = ( ) Ansatz with finite functions to approximate exact solution E(x) E N (x) = X N i= c i E i (x) ( ) Finite Differences (FDM) Taylor expansion of E(x)around each grid point to derive stencils for the differential operators. Standard FDTD: Yee algorithm equivalence with variational form, where ansatz functions are Z-splines and test functions are Dirac distributions, J.T. Becerra Sagredo, SAM Report 3-, ETH Zurich. Galerkin method for Finite Elements and Plane Wave Expansion multiply ( ) with test function F (x) integrate over the domain Ω, use integration by parts formula and insert boundary conditions insert ansatz ( ) and each basis function E i (x) as test function Z µ (x) E N (x) E i (x) ε(x, ω) ω E N (x) E i (x)dx = Ω Same procedure for Maxwell source problem in time or frequency domain. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
8 Space discretisation of Maxwell s equation E.g. Maxwell eigenvalue problem µ (x) E(x) ε(x, ω) ω E(x) = ( ) Ansatz with finite functions to approximate exact solution E(x) E N (x) = X N Galerkin method for Finite Elements and Plane Wave Expansion multiply ( ) with test function F (x) i= c i E i (x) ( ) integrate over the domain Ω, use integration by parts formula and insert boundary conditions insert ansatz ( ) and each basis function E i (x) as test function Z µ (x) E N (x) E i (x) ε(x, ω) ω E N (x)e i (x)dx = Ω Same procedure for Maxwell source problem in time or frequency domain. Strong theory of variational problems, i.e. stability, consistency and convergence coming from functional analysis and mathematical physics for the continuous and numerical analysis for the discrete spaces. Plane wave expansion originately used by mathematical physicists (Schrödinger equation), where Finite Elements come from civil engineering. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
9 Eigenvalue problem in D 6 TE mode at k =, for ε = 9 in [, ] 4 (x) ε(x) h = ω h(x) in [, ] with first non-zero eigenvalue ω = arccos(/4).5.5 ε(x) p = p = p = 3 8 x h = h = / h = / K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
10 Eigenvalue problem in D TE mode at k =, for ε = 9 in [, ] (x) ε(x) h = ω h(x) in [, ] with first non-zero eigenvalue ω = arccos(/4) ωh ω /ω h-fem p = h-fem p = h-fem p = 3 p = p = p = 3 dof h = h = / h = / K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
11 Eigenvalue problem in D TE mode at k =, for ε = 9 in [, ] (x) ε(x) h = ω h(x) in [, ] with first non-zero eigenvalue ω = arccos(/4) p = p = p = 3 ωh ω /ω h-fem p = h-fem p = h-fem p = 3 p-fem dof h = = p-fem with exponential convergence h = / h = / K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
12 Eigenvalue problem in D TE mode at k =, for ε = 9 in [, ] (x) ε(x) h = ω h(x) in [, ] with first non-zero eigenvalue ω = arccos(/4) p = p = p = 3 ωh ω /ω h-fem p = h-fem p = h-fem p = 3 p-fem dof h = = p-fem with exponential convergence h = / h = / ωh ω /ω 3 5 h-fem p = h-fem p = h-fem p = 3 7 p-fem dof K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
13 Eigenvalue problem in D Smooth, periodic dielectricity Plane Wave expansion ε(x).5.5 ωh ω /ω Jumping dielectricity ε(x) x dof 8 4 ε(x) 6 4 ωh ω /ω x dof Global plane wave expansions exponential convergence for smooth material coefficients convergence only algebraic for jumping material coefficients R. Norton, PhD thesis, University of Bath, Sept. 8. Need to resolve material interfaces spectral methods like p-fem conserves exponential convergence K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 6 / 35
14 Eigenvalue problem in D Smooth, periodic dielectricity Plane Wave expansion ε(x).5.5 ωh ω /ω Jumping dielectricity ε(x) x dof Plane Wave expansion versus h- and p-fem 8 4 ε(x) x Global plane wave expansions ωh ω /ω 8 6 PWM h-fem p = h-fem p = h-fem p = 3 p-fem exponential convergence for smooth material coefficients convergence only algebraic for jumping material coefficients R. Norton, PhD thesis, University of Bath, Sept. 8. Need to resolve material interfaces spectral methods like p-fem conserves exponential convergence K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 6 / 35 dof
15 Resolving material interfaces in D FDM / PWM Level Level Level Level 3 h-fem Level Level Level Level 3 Resolving by mesh refinement p-fem p = p = p = 3 p = 4 Resolving by curved elements K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 7 / 35
16 Eigenvalue problems in D Photonic crystal bandstructure dispersive material Photonic crystal composed of cylinders of a polymer with embedded quantum dots in Silicon D. Hermann et. al., Phys. Rev. B 77, 8. Comparison of different FEM and quadratic EWP in k FreeFEM++ with p = without curved cells Code of M. Richter (Univ. Karlsruhe) with p = and curved cells Concepts p-fem with curved cells 3 εp ηα(ω) ε (ω) = ε p + να(ω) α(ω) = ε d (ω) ε p ε d (ω)+ ε p ε d (ω) = + ω p ω ω iγω Relative error (log ) (p ref.) p-fem h-fem h-fem (curved) N (log ) C. Engström, C. Hanfer and K. Schmidt, CTN 6, p , 9 K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 8 / 35
17 Eigenvalue problems in D Photonic crystal wave-guide Supercell approach in a finite cell around a Photonic Cyrstal W-waveguide Joint work with Roman Kappeler (IfE, ETH) S.Soussi, SIAM Journal on Numerical Analysis 43, p. 75, 5. a a a Floquet-transformation in periodicity direction Periodic boundary conditions on horizontal boundaries Eigenvalue problems linear in λ = ω or quadratic in k. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 9 / 35
18 Eigenvalue problems in D Photonic crystal wave-guide Comparison with other codes MPB MIT Photonic-Bands (Group of St.G.Johnson, MIT) Expansion in plane waves, mostly used code for band structure computation. St.G.Johnson and J.D.Joannopoulos, Optics Express 8, p. 73 9,. Comsol FE solver with broad class of applications Linear eigenvalue problem for ω(k) 4 4 error error Concepts mesh for W wave-guide with 5 rows of holes (r =.3a) and ε(x) = degrees of freedom time [s] uniformly sampled on D Brillouin zone (5 points, 3 bands) Sun Fire 46, Dual Core AMD Opteron.53 GHz, Linux K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
19 Eigenvalue problems in D Photonic crystal bandstructure different pattern Contrast ε =, Diameter of cylinder or length of square.9 a, Convergence of p-fem at k = (π/, )/a ωh ω /ω 4 6 ωh ω /ω cylindrical holes dielectric veins curved dielectric rods ( dof) /3 8 cylindrical holes dielectric veins curved dielectric rods ( dof) /3 K.Schmidt, P.Kauf, CMAME 98, p , Mar 9. Weaker convergence behaviour for p-fem in case of material corners Theory : we can expect in general only algebraic convergence K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
20 hp-adaptive refinement in D hp-adaptive FEM Combination of cell refinement and polynomial order enrichment Geometric mesh Exponential convergence also with sharp corners expected λ λ N C exp( βn /3 ). Layers of constant polynomial order around interface corner Polynomial order is raised away from interface corner Cells are refined close to interface corner Occurrence of hanging nodes, possibly multiple hanging nodes per edge for different refinement in two materials K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
21 hp-adaptive refinement in D Maxwell eigenvalue problem BENCHMAX, M. Dauge Benchmark computations for Maxwell equations for the approximation of highly singular solutions TE transmission problem in a square approximated with hp-fem Z Z curl D e(x) curl D e (x)dx = ω ε(x)e(x) e (x)dx Ω Ω Geometry DomE Mesh Convergence of third eigenvalue ε ε ε ε.5.5 ε =, ε =.5.5 relative error - - weighted reg., s = edge elements equivalent problem ndof Ph. Frauenfelder, PhD thesis, ETH Zurich, 4. Approximation of highly singular eigenfunctions with exponential convergence K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
22 hp-adaptive refinement in D Photonic crystal bandstructure Dielectric veins of thickness. a with ε = 8.9 and air holes (ε = ) Coarse Mesh with 9 cells p-fem with algebraic convergence, as expected hp-fem with exponential convergence, as expected relative errors TE mode p-fem N relative errors 5 TM mode N K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
23 hp-adaptive refinement in D Photonic crystal bandstructure Dielectric veins of thickness. a with ε = 8.9 and air holes (ε = ) Coarse Mesh with 9 cells p-fem with algebraic convergence, as expected hp-fem with exponential convergence, as expected relative errors TE mode p-fem hp-fem start with p= N relative errors 5 TM mode N K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
24 Space discretisation by h-, p- and hp-fem Approximation of exact solution E(x) by a finite number of functions P N i= c i E i (x) h-fem, Finite Differences,... mesh refinement p-fem, Global Plane Wave Expansion,... increasing order hp-fem contains both, each cell can be refined or function order can be increased take best choice (a-priori or a-posteriori) in D and 3D material interfaces has to be resolved curved cells hp-fem is able to compute with asymptotically exponential convergence K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
25 Space discretisation by h-, p- and hp-fem Approximation of exact solution E(x) by a finite number of functions P N i= c i E i (x) h-fem, Finite Differences,... mesh refinement p-fem, Global Plane Wave Expansion,... increasing order hp-fem contains both, each cell can be refined or function order can be increased take best choice (a-priori or a-posteriori) in D and 3D material interfaces has to be resolved curved cells hp-fem is able to compute with asymptotically exponential convergence Wave propagation Waves have to be resolved ( 3 dof per wavelength) before we have convergence towards right solution (independent of discretisation scheme). Then, p-fem is more effective (better increase p = 64 7 on a one-cell mesh than refinement of mesh with 64 cells (p = ) three times ( 5 cells). 3 3 K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 6 / 35
26 Space discretisation by h-, p- and hp-fem Approximation of exact solution E(x) by a finite number of functions P N i= c i E i (x) h-fem, Finite Differences,... mesh refinement p-fem, Global Plane Wave Expansion,... increasing order hp-fem contains both, each cell can be refined or function order can be increased take best choice (a-priori or a-posteriori) in D and 3D material interfaces has to be resolved curved cells hp-fem is able to compute with asymptotically exponential convergence Resolving material interfaces Geometries with small details small cells? smallness in or directions anisotropic cells with anisotropic function order domain decomposition and parallel computing Small periodic structures (finite photonic crystals) basis functions with micro-structure (generalised FEM) Yesterday s talk of H. Brandsmeier, joint project with Prof. Ch. Schwab. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 7 / 35
27 Space discretisation by h-, p- and hp-fem Approximation of exact solution E(x) by a finite number of functions P N i= c i E i (x) h-fem, Finite Differences,... mesh refinement p-fem, Global Plane Wave Expansion,... increasing order hp-fem contains both, each cell can be refined or function order can be increased take best choice (a-priori or a-posteriori) in D and 3D material interfaces has to be resolved curved cells hp-fem is able to compute with asymptotically exponential convergence Resolving material interfaces Geometries with small details small cells? smallness in or directions anisotropic cells with anisotropic function order domain decomposition and parallel computing Small periodic structures (finite photonic crystals) basis functions with micro-structure (generalised FEM) Yesterday s talk of H. Brandsmeier, joint project with Prof. Ch. Schwab. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 7 / 35
28 Classification of FEM Finite Element Methods Continuous FEM BEM DG-FEM R Ω µ E E ε ω E E dx = weak solution E H(curl, Ω) FE space is conforming V N H(curl, Ω) (continuous) FEM of boundary integral eq. meshing of surfaces smaller matrices double integral with kernel matrices are dense Fast BEM singular integrals due to singular kernel non-conforming, discontinuous spaces discontinuity over facettes vanishes for N larger matrices, less entries for large p additional terms over facettes or Lagrange multiplicator Edge elements Existing h-, p- and hp-version of Continuous FEM, BEM and DG-FEM Generalised FEM use of problem-adapted basis functions, e.g. local plane wave, singular functions near corners, functions with micro-structure K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 8 / 35
29 Implemenation in Concepts Numerical C++ class library Concepts Started in 995 and further developed at Seminar for Applied Mathematics of ETH Zurich, mainly by Ch. Lage, Ph. Frauenfelder, A.-M. Matache, G. Schmidlin, K. Schmidt, P. Kauf, H. Brandsmeier and severals students. Modular programming Use of interfaces defined in base classes. Re-usability in the development, e.g. various bilinearforms on curved cells. Generality, e.g. meshes with curved cells of mixed type. Concept Oriented Design using mathematical principles. Ph.Frauenfelder and Ch.Lage, Math. Model. Numer. Anal. 36, p ,. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 9 / 35
30 Implemenation in Concepts Numerical C++ class library Concepts hp-adaptive FEM on tensor product elements (quadrilateral, hexahedron) with anisotropic polynomial degree and arbitrary number of hanging nodes per facette for scalar problems in H (Ω) in D, D and 3D, for vectorial problems in (H (Ω)) d in D and 3D, for vectorial problems in H(curl, Ω) in D, for mixed problems in H(curl, Ω) H (Ω) in D. Fast BEM (panel clustering and wavelets) in 3D on triangular surface meshes. Fast writeable sparse matrix format based on hash tables. Indirect linear solver CG, GMRES (intern). Direct linear solver SuperLU, Umfpack, Pardiso (interfaces to libraries). Eigenvalue solver ARPACK, Trilinos (interfaces to libraries). Various time integration schemes, e.g. Euler, Newmark, Runge-Kutta. Solution output or function of the solution, like its gradient, to Matlab, Tecplot, Gnuplot. Post-processing of the solution, e.g. ohmic power loss by integrating Im(ε)E. Class documentation in Doxygen. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
31 Implemenation in Concepts Mesh Import in D and 3D Coord.dat Elm.dat Attrib.dat Radia.dat Concepts Input Generator in D Joint work with R. Kappeler, S. Belfanti, B. Jin (IfE, ETH) Python library for scripting of geometry input Graphical User Interface with Qt Pre-defined pattern, like this K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
32 Implemenation in Concepts Input files HexaWguide-n5-r3-withAirL4.cig string inputpath "../cig/hexa_w_guides" string coordfile "Hexa_W_guide_n5_r3_withAirL4_Coord.dat" string elemfile "Hexa_W_guide_n5_r3_withAirL4_Elms.dat" string attribfile "Hexa_W_guide_n5_r3_withAirL4_Attr.dat" string radiafile "Hexa_W_guide_n5_r3_withAirL4_EdgRadia.dat" string edgecorrfile "Hexa_W_guide_n5_r3_withAirL4_EdgCorr.dat" array double epsilon { } string mode "TE" string refinement "p3, p in 3" string nonrefl_bdattrib "9" double epsarpack e- bool graphics true int graphicpoints 7 Concepts mesh for super-cell of W wave-guide Call for PC bandstructure calculation: PCEVsolve -f HexaWguide-n5-r3-withAirL4.cig K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. / 35
33 Implemenation in Concepts Main C++ file. Reading the input file. Importing the mesh M 3. Defining the FE space on M and perform (adaptive) refinements in h and p 4. Defining the material formula, e.g. for µ(x) and ε(x) 5. Defining the involved bilinearforms with their material or the linearforms for sources. 6. Assembling of the system matrices and vector of the rhs. 7. Solving Time integration scheme with possibly several time solving of linear systems for Maxwell source problem in time domain. Solving a linear system for Maxwell source problem in time domain. Solving an matrix eigenvalue problem for Maxwell eigenvalue problem. 8. Post-processing Graphical output of the solution or function of the solution. Computation of physical quantities like Ohms losses or transmission coefficients. Flexibility by writing problem adapted main files. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
34 Implemenation in Concepts Scalar continuous basis functions (subspace of H (Ω)) vertex function edge function inner or bubble function Representation by local shape function on reference cell [, ] mapped to physical cell.5.5 F K Φ i (x) K = φ k,l (F K x), Continuity is assured by T- or connectivity matrix for each element K Φ i (x) K = X [T K ] kl,i φ kl (F K x), K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
35 Implemenation in Concepts Scalar continuous basis functions (subspace of H (Ω)) vertex function edge function inner or bubble function Representation by local shape function on reference cell [, ] mapped to physical cell.5.5 F K Φ i (x) K = ±φ k,l (F K x), Continuity is assured by T- or connectivity matrix for each element K Φ i (x) K = X [T K ] kl,i φ kl (F K x), K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
36 Implemenation in Concepts Scalar continuous basis functions (subspace of H (Ω)) vertex function edge function inner or bubble function Representation by local shape function on reference cell [, ] mapped to physical cell.5.5 F K Φ i (x) K = φ k,l (F K x), Continuity is assured by T- or connectivity matrix for each element K Φ i (x) K = X [T K ] kl,i φ kl (F K x), K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 4 / 35
37 Implemenation in Concepts Hierarchical basis of H (Ω) We use Ainsworth Coyle s basis with anisotropic polynomial degrees on quadrilateral or hexahedral elements Definition on the reference element η N i,j (ξ, η) = P i (ξ)p j (η) i =,..., p +, j =,..., p + ˆK - ξ P (ξ) = ξ, P (ξ) = +ξ, P i (ξ) = ξ +ξ P, i (ξ) P, i (ξ)... Jacobi polynomials # shape functions = (p + )(p + ) = (p + ) for p = p Vertex function Occuring of P ( ) or P ( ) in both directions K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
38 Implemenation in Concepts Hierarchical basis of H (Ω) We use Ainsworth Coyle s basis with anisotropic polynomial degrees on quadrilateral or hexahedral elements Definition on the reference element η N i,j (ξ, η) = P i (ξ)p j (η) i =,..., p +, j =,..., p + ˆK - ξ P (ξ) = ξ, P (ξ) = +ξ, P i (ξ) = ξ +ξ P, i (ξ) P, i (ξ)... Jacobi polynomials # shape functions = (p + )(p + ) = (p + ) for p = p Edge function Occuring of only one of P ( ) or P ( ) in one directions K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
39 Implemenation in Concepts Hierarchical basis of H (Ω) We use Ainsworth Coyle s basis with anisotropic polynomial degrees on quadrilateral or hexahedral elements Definition on the reference element η N i,j (ξ, η) = P i (ξ)p j (η) i =,..., p +, j =,..., p + ˆK - ξ P (ξ) = ξ, P (ξ) = +ξ, P i (ξ) = ξ +ξ P, i (ξ) P, i (ξ)... Jacobi polynomials # shape functions = (p + )(p + ) = (p + ) for p = p Inner or bubble function Occuring of P i ( ) with i > K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 5 / 35
40 Implemenation in Concepts Hierarchical basis of H(curl, Ω) We use Ainsworth Coyle s basis with anisotropic polynomial degrees on quadrilateral elements Definition on the reference element η ˆK N i,j, (ξ, η) = P i+ (ξ)p (η)` i =,..., p, j j =,..., p + N i,j, (ξ, η) = P i (ξ)p j+ (η)` i =,..., p +, j =,..., p - ξ P (ξ) = ξ, P (ξ) = +ξ, P i (ξ) = ξ +ξ P, i (ξ) P, i (ξ)... Jacobi polynomials # shape functions = (p + )(p + ) + (p + )(p + ) = (p + )(p + ) for p = p Edge function Occuring of P ( ) or P ( ) K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 6 / 35
41 Implemenation in Concepts Hierarchical basis of H(curl, Ω) We use Ainsworth Coyle s basis with anisotropic polynomial degrees on quadrilateral elements Definition on the reference element η ˆK N i,j, (ξ, η) = P i+ (ξ)p (η)` i =,..., p, j j =,..., p + N i,j, (ξ, η) = P i (ξ)p j+ (η)` i =,..., p +, j =,..., p - ξ P (ξ) = ξ, P (ξ) = +ξ, P i (ξ) = ξ +ξ P, i (ξ) P, i (ξ)... Jacobi polynomials # shape functions = (p + )(p + ) + (p + )(p + ) = (p + )(p + ) for p = p Inner or bubble function Occuring of P i ( ) with i > K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 6 / 35
42 Implemenation in Concepts Evaluation on the physical element Nodal element Edge element N i,j (F K x) J N i,j,k (F K x) Element mass matrix for nodal elements Z Z N i,j (F K x)n i,j (FK x) dx = N i,j (ξ, η)n i K ˆK,j (ξ, η) det J dξdη Element stiffness matrix for nodal elements Z Z N i,j (F K x) N i,j (FK x) dx = ( ˆ N i,j (ξ, η)) J J ˆ N i K ˆK,j (ξ, η) det J dξdη Element mass matrix for edge elements Z Z N i,j,k (F K x) J J N i,j,k (FK x) dx = K ˆK N i,j,k (ξ, η) J J N i,j,k (ξ, η) det J dξdη Element stiffness matrix for edge elements Z ( J N i,j,k (F K x)) J N i,j,k (FK x) dx K Z = det J ˆ N i,j,k (ξ, η) det J ˆ N i ˆK,j,k (ξ, η) det J dξdη Numerical quadrature on the reference cell ˆK K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 7 / 35
43 Implemenation in Concepts Evaluation on the physical element Nodal element Edge element N i,j (F K x) J N i,j,k (F K x) Element mass matrix for nodal elements Z Z N i,j (F K x)n i,j (FK x) dx = N i,j (ξ, η)n i K ˆK,j (ξ, η) det J dξdη Element stiffness matrix for nodal elements Z Z N i,j (F K x) N i,j (FK x) dx = ( ˆ N i,j (ξ, η)) J J ˆ N i K ˆK,j (ξ, η) det J dξdη Element mass matrix for edge elements Z Z N i,j,k (F K x) J J N i,j,k (FK x) dx = K ˆK N i,j,k (ξ, η) J J N i,j,k (ξ, η) det J dξdη Element stiffness matrix for edge elements Z ( J N i,j,k (F K x)) J N i,j,k (FK x) dx K Z = ˆ N i,j,k (ξ, η) ˆ N i ˆK,j,k (ξ, η) det J dξdη Numerical quadrature on the reference cell ˆK K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 7 / 35
44 Implemenation in Concepts Sparsity pattern for a rectangle curl D -curl D element matrix identity element matrix K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 8 / 35
45 Implemenation in Concepts Sparsity pattern for a parallelogramme curl D -curl D element matrix identity element matrix K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 9 / 35
46 Implemenation in Concepts Sparsity pattern for a trapezoid curl D -curl D element matrix identity element matrix Sum factorisation reduces computational effort from O(p 6 ) to O(p 5 ) K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
47 Implemenation in Concepts Condition number of A + M in dependence of polynomial degree 8 condition number shape functions p A... element stiffness matrix M... element mass matrix K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
48 Implemenation in Concepts Curved cells η x 3 (ξ) p ˆK ξ F K p 3 x 4 (η) p x (ξ) p x (η) Blending Mapping as linear extension of edge parametrisation into the interior F K ξ = ( η) x (ξ) + ( + ξ) x (η) + ( + η) x 3 (ξ) + ( ξ) x 4 (η) 4 ( ξ)( η) p 4 ( + ξ)( η) p 4 ( + ξ)( + η) p... Jacobian is linear combination of derivates of edge parametrisations and itself «x x J = = «( η) x ξ η (ξ) + x (η) +... x (η) + ( + ξ) x (η) +... Blending technique needs formula for edge parametrisation and derivates Currently implemented: Circular, Ellipsoid, Parabel, Parallel of other one K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 3 / 35
49 Implemenation in Concepts General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) [B] ij = b(φ j, Φ i ) = X X b KK (Φ j, Φ K i K ) K K FEM : K = K DG-FEM : K, K are neighbours BEM : all combination of K and K K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 33 / 35
50 Implemenation in Concepts General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) [B] ij = b(φ j, Φ i ) = X X b KK (Φ j, Φ K i K ) K K = X dim(k) dim(k X X ) X b [T K ] m,j φ K m, [T K ] n,i φ K A n K K m= n= FEM : K = K DG-FEM : K, K are neighbours BEM : all combination of K and K K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 33 / 35
51 Implemenation in Concepts General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) [B] ij = b(φ j, Φ i ) = X X b KK (Φ j, Φ K i K ) K K FEM : K = K = X K X K DG-FEM : K, K are neighbours dim(k) X m= dim(k ) X [T K ] m,j [T K ] n,i b KK n= BEM : all combination of K and K φ K m, φk n K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 33 / 35
52 Implemenation in Concepts General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) [B] ij = b(φ j, Φ i ) = X X b KK (Φ j, Φ K i K ) K K FEM : K = K = X K X K dim(k) X m= DG-FEM : K, K are neighbours dim(k ) X n= BEM : all combination of K and K [T K ] m,j [T K ] n,i {z } connectivity local to global [B KK ] nm, {z } element matrix K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 33 / 35
53 Implemenation in Concepts General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) [B] ij = b(φ j, Φ i ) = FEM : K = K = X K X K dim(k) X m= DG-FEM : K, K are neighbours dim(k ) X n= BEM : all combination of K and K [T K ] m,j [T K ] n,i {z } connectivity local to global [B KK ] nm, {z } element matrix,, over of over of Get by calling with, Multiply with, and copy result into K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 33 / 35
54 Implemenation in Concepts Main C++ file. Reading the input file. Importing the mesh M 3. Defining the FE space on M and perform (adaptive) refinements in h and p 4. Defining the material formula, e.g. for µ(x) and ε(x) 5. Defining the involved bilinearforms with their material or the linearforms for sources. 6. Assembling of the system matrices and vector of the rhs. 7. Solving Time integration scheme with possibly several time solving of linear systems for Maxwell source problem in time domain. Solving a linear system for Maxwell source problem in time domain. Solving an matrix eigenvalue problem for Maxwell eigenvalue problem. 8. Post-processing Graphical output of the solution or function of the solution. Computation of physical quantities like Ohms losses or transmission coefficients. Flexibility by writing problem adapted main files. K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 34 / 35
55 Conclusion Space discretisation by h-, p- and hp-fem Mesh refinement or increasing the (polynomial) order Variational formulation Linear system of equations Resolving material interfaces essential p-fem gives exponential convergence for smooth interfaces hp-fem = suitable combination of h- and p-refinement expon. convergence for interface corners Resolution of number of wave lengthes after coarse resolution an increasing p is more efficient Geometries with small details anisotropic cells, domain decomposition, gfem with micro-structure Applications Bandstructure of disperse photonic crystal (with C. Engström, C. Hafner) Modes in Photonic Crystal Waveguides (with R. Kappeler) Benchmark for highly singular Maxwell solutions (with P. Frauenfelder) PC Bandstructure for polygonial interfaces (with P. Kauf) Implementation in Concepts Modular Programming, Concepts oriented design hp-fem in D and 3D with hanging nodes (edge elements in D) Basis functions with local support Hierarchical basis based on Jacobi polynomials on reference cell Numerical quadrature on reference cell considering shape of curved cell (Blending technique) General Algorithm for Matrix Assembling FEM, DG-FEM and BEM (hp, hanging nodes) K.Schmidt 8th July 9, NanoOptics Workshop, Zürich p. 35 / 35
Modelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationFinite Element Methods for Optical Device Design
DFG Research Center Matheon mathematics for key technologies and Zuse Institute Berlin Finite Element Methods for Optical Device Design Frank Schmidt Sven Burger, Roland Klose, Achim Schädle, Lin Zschiedrich
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
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 informationComputation of the band structure of two-dimensional Photonic Crystals with hp Finite Elements
Computation of the band structure of two-dimensional Photonic Crystals with hp Finite Elements K. Schmidt, P. Kauf Seminar for Applied Mathematics, ETH Zurich, Rämisstrasse, 892 Zürich, Switzerland Abstract
More informationNonlinear Eigenvalue Problems: An Introduction
Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction
More informationParallel Discontinuous Galerkin Method
Parallel Discontinuous Galerkin Method Yin Ki, NG The Chinese University of Hong Kong Aug 5, 2015 Mentors: Dr. Ohannes Karakashian, Dr. Kwai Wong Overview Project Goal Implement parallelization on Discontinuous
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More 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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
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 informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
More informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationTrefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations
Trefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with
More informationGypsilab : a MATLAB toolbox for FEM-BEM coupling
Gypsilab : a MATLAB toolbox for FEM-BEM coupling François Alouges, (joint work with Matthieu Aussal) Workshop on Numerical methods for wave propagation and applications Sept. 1st 2017 Facts New numerical
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More 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 informationEffects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling
Effects of tool eccentricity on wave dispersion properties in borehole acoustic logging while drilling Yibing Zheng and M. Nafi Toksöz Earth Resources Laboratory Dept. of Earth, Atmospheric, and Planetary
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
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 informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
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 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 informationEfficient hp-finite elements
Efficient hp-finite elements Ammon Washburn July 31, 2015 Abstract Ways to make an hp-finite element method efficient are presented. Standard FEMs and hp-fems with their various strengths and weaknesses
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationMethod of Finite Elements I
Method of Finite Elements I PhD Candidate - Charilaos Mylonas HIL H33.1 and Boundary Conditions, 26 March, 2018 Institute of Structural Engineering Method of Finite Elements I 1 Outline 1 2 Penalty method
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationAdding a new finite element to FreeFem++:
Adding a new finite element to FreeFem++: the example of edge elements Marcella Bonazzoli 1 2, Victorita Dolean 1,3, Frédéric Hecht 2, Francesca Rapetti 1, Pierre-Henri Tournier 2 1 LJAD, University of
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationJournal of Computational Physics
Journal of Computational Physics 23 (211) 8573 862 Contents lists available at SciVerse ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp An eigen-based high-order
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 informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
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 informationFEM: Domain Decomposition and Homogenization for Maxwell s Equations of Large Scale Problems
FEM: and Homogenization for Maxwell s Equations of Large Scale Problems Karl Hollaus Vienna University of Technology, Austria Institute for Analysis and Scientific Computing February 13, 2012 Outline 1
More informationL e c t u r e. D r. S a s s a n M o h a s s e b
The Scaled smteam@gmx.ch Boundary Finite www.erdbebenschutz.ch Element Method Lecture A L e c t u r e A1 D r. S a s s a n M o h a s s e b V i s i t i n g P r o f e s s o r M. I. T. C a m b r i d g e December
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationWaveguide Propagation Modes and Quadratic Eigenvalue Problems
Waveguide Propagation Modes and Quadratic Eigenvalue Problems Prof. André Nicolet, Institut Fresnel, UMR CNRS 6133, Université Paul Cézanne, Marseille, France. Prof. Christophe Geuzaine, Mathematics Department,
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationCourse Notes: Week 1
Course Notes: Week 1 Math 270C: Applied Numerical Linear Algebra 1 Lecture 1: Introduction (3/28/11) We will focus on iterative methods for solving linear systems of equations (and some discussion of eigenvalues
More informationPOLYNOMIAL APPROXIMATIONS OF REGULAR AND SINGULAR VECTOR FIELDS WITH APPLICATIONS TO PROBLEMS OF ELECTROMAGNETICS
POLYNOMIAL APPROXIMATIONS OF REGULAR AND SINGULAR VECTOR FIELDS WITH APPLICATIONS TO PROBLEMS OF ELECTROMAGNETICS Alex Besalov School of Mathematics The University of Manchester Collaborators: Norbert
More informationNumerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material
Chapter 30 Numerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material Hiroshi Maeda Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50555
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
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 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 informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationComputational Methods for Engineering Applications II. Homework Problem Sheet 4
S. Mishra K. O. Lye L. Scarabosio Autumn Term 2015 Computational Methods for Engineering Applications II ETH Zürich D-MATH Homework Problem Sheet 4 This assignment sheet contains some tasks marked as Core
More informationJasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015
Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern
More informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationA 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 informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationGeoPDEs. An Octave/Matlab software for research on IGA. R. Vázquez. IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca
GeoPDEs An Octave/Matlab software for research on IGA R. Vázquez IMATI Enrico Magenes, Pavia Consiglio Nazionale della Ricerca Joint work with C. de Falco and A. Reali Supported by the ERC Starting Grant:
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationContinuum mechanics of beam-like structures using one-dimensional finite element based on Serendipity Lagrange cross-sectional discretisation, Mayank Patni, Prof. Paul Weaver, Dr Alberto Pirrera Bristol
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationComputational Electromagnetics Definitions, applications and research
Computational Electromagnetics Definitions, applications and research Luis E. Tobón Pontificia Universidad Javeriana Seminario de investigación Departamento de Electrónica y Ciencias de la Computación
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
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 informationS-matrix approach for calculations of the optical properties of metallic-dielectric photonic crystal slabs
S-matrix approach for calculations of the optical properties of metallic-dielectric photonic crystal slabs N. I. Komarevskiy1,2, T. Weiss3, and S. G. Tikhodeev2 1 Faculty of Physics, Lomonosov Moscow State
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationApproximate tensor-product preconditioners for very high order discontinuous Galerkin methods
Approximate tensor-product preconditioners for very high order discontinuous Galerkin methods Will Pazner 1 and Per-Olof Persson 2 1 Division of Applied Mathematics, Brown University, Providence, RI, 02912
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationarxiv:physics/ v1 [physics.optics] 8 Apr 2005
FEM Modelling of 3D Photonic Crystals and Photonic Crystal Waveguides Sven Burger ab, Roland Klose a, Achim Schädle a, Frank Schmidt ab, and Lin Zschiedrich ab a Zuse Institute Berlin, Takustraße 7, D
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
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 informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More information