Transactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
|
|
- Sarah Lindsey Curtis
- 5 years ago
- Views:
Transcription
1 Application of the panel clustering method to the three-dimensional elastostatic problem K. Hayami* and S. A. Sauter** ^Department of Mathematical Engineering and Information Physics, School ofengineering, University of Tokyo, 113 Tokyo, Japan t. u-tokyo. ac.jp **Lehrstuhlfur Praktische Mathematik, Mathematisches Seminar Bereichll, Christian-Albrechts-Universitat zu Kiel, D Kiel, Germany Abstract A formulation of the panel clustering method for the three-dimensional elastostatic problem is presented, in order to reduce the computational work and memory required for the boundary element method. In order to make the necessary polynomial expansions for the clustering simple, the fundamental solution for the displacement and traction components can be expressed as a linear combination of partial derivatives of the distance r between the observation point and thefieldpoint. However, this requires a higher order expansion compared to directly expanding the usual expression for the kernels. 1 Introduction Although the Boundary Element Method (BEM) enjoys the advantage of the boundary only discretization, a serious computational difficulty arises due to its dense matrix formulation, particularly for large scale threedimensional problems arising in engineering. This is because the method requires O(A^) memory and O(A^) computational work using the conventional approach, where TV is the number of unknowns. The situation is even worse for the three-dimensional elastostatic problem, where the number of unknowns is three times that of the potential problem, since the unknowns at each node is a vector instead of a scalar. Rokhlin[8] and Hackbusch and Nowak[5] independently proposed the
2 626 Boundary Elements multipole method and the panel clustering method, respectively, in order to overcome this difficulty. The main idea is to approximate the far field using multipole or polynomial expansions around a centre of a cluster of panels or boundary elements, thus reducing the O(A^) dense matrix vector multiplication to a 0(7V(log AT)^) sparse matrix vector multiplication for each iteration of the iterative linear solver, where d is the dimension of the space. Yam ad a and Hayami[ll] proposed a multipole boundary element method for two-dimensional elastostatics. In this paper, we will present a formulation of the panel clustering method for the three-dimensional elastostatic problem [6]. 2 The boundary integral equation for 3-D elastostatics The boundary integral equation for the three-dimensional (linear, isotropic) elastostatic problem is given by [1] (x, yx(y)df(y) = (/== 1,2,3), (1) where we have used Einstein's convention for the summation over the repeated index k 1,2,3, F is the boundary of the domain under consideration, Uk,pk are the displacement and traction components, respectively, c/jk(x) = ^Sik when F is smooth at x, and the body force term has been neglected. u*j.(x,y) is the fundamental solution corresponding to the A:-th component of the displacement at y due to a unit point load in the /-direction at x, which is given by where JJL is the shear modulus, v is the Poisson's ratio, r = \ y x and r,& = <9r/<%, where y = (2/1,2/2,2/3)^. p*j.(x,y) is the traction component at y corresponding to u*^,(x,y), which is given by k - «, + 3r,r k} + (1-2i/)(n/r,* - OTT (1 V ' JT (3) where n\ is the component of the unit outward normal vector at y G F. (2)
3 3 The boundary element discretization Next, the boundary integral equation (1) is discretized by discretizing the boundary F into boundary elements F^, (a = l,...,n). If constant elements are used, and x* is a point representing F^, we obtain t4(x,y)dr(y), (a = l,...,n, ( = 1,2,3) (4) or where iad 3 a 8 8 / t o o\ /_\ %<=^%, (a=l,...,^, / = 1,2,3), (5) 2 ik t4(x,y)dr(y), T/3 and Einstein's convention is now also applied for the summation over the elements f3 = 1,..., n. Equation (5) can be rearranged as where JI IK k J Ik i k ' \ ) ' * * 5 7 "> ") } i \ ) for u], u^ : unknown boundary displacement component, for uf = uf : given boundary displacement component, ~a8 _ a8 -aft _ ^ 9ik = 9ik > 9ik = v for pi pi : unknown boundary traction component, and -a/3 _ n -a/3 _ a/3 9 ik = u, 9ik = 9 ik for pf = pf : given boundary traction component. Equation (6) is a system of linear equations for the unknown boundary displacement and traction components uf and %. Since its matrix is dense and nonsymmetric, it is usually solved using LU-decomposition, which leads to O(A^) computational work and O(A^) memory, where TV = 3n is the number of unknowns. Alternatively, one could use Krylov subspace iterative solvers for nonsymmetric matrices, such as the GMRES(Generalized Minimal Residual)
4 628 Boundary Elements method[9]. The iterates of this method is guaranteed to converge to the true solution within TV iterations if exact arithmetic is used. More important, the system of equations (6) arising from boundary integral equations such as equation (1) is better conditioned compared to its finite difference and finite element counterparts. That is, for the same mesh (element) size A, the condition number for BEM is 0(1) or 0(ft~*), whereas for the FEM, it is 0(/i~^)[3]. Hence, if a suitable preconditioner is used, the method should converge within M «N iterations for the BEM. Note also that it makes sense to do the full orthogonalization in the GMRES without restarting, since the system is better conditioned and smaller for BEM compared to FEM. In this method, as with all iterative methods, the dominant part of the computation is the dense matrix-vector multiplication (inner product): p" t4(x«,y)dr(y), (a =!,...,«, /= 1,2,3) (8) which costs O(^) each. Here, v% and q% are the components of the iteration vector corresponding to the unknown components itf and pf, respectively. The summations over j3' and /?" are taken for the elements F^/ and F^// where the components Uk and p& are unknown, respectively. Hence, the amount of computational work is reduced from O(TV^) to O(MN^}, but the memory required is still O(TV^), and it is this memory bottle-neck that hinders the solution of large scale problems using the boundary element method. 4 The panel clustering method The reason why the matrix-vector product of equation (8) is dense is because the observation point x" on the element F^ is related to all the elements F/? on the boundary through the kernels p*^.(x",y) and u*j.(x",y). In this paper, we will apply the panel clustering method [5, 10] in order to reduce the required memory and computational work for computing the above mentioned matrix-vector product. The method makes use of polynomial expansions to approximate the integral kernels for cluster of elements which are sufficiently far from the observation point, thus reducing the amount of computation and required memory, and at the same time achieving the required accuracy which is consistent with the error introduced by discretizing the original boundary integral equation (1) into the boundary element equations (4).
5 Boundary Elements 629 Let 0=1 (9) represent the integral operator or inner product corresponding to the second or third term of the right hand side of equation (8). (u(y) is constant over each element F/? in this case.) Following [5], in the farfield (i.e. when x y is sufficiently large) we will replace fc(x,y) by a piecewise (elementwise) polynomial fc(x,y) with respect to y, which may be chosen as an expansion around a suitable centre y GR3: 6(x,y)= ^KXx;y )$Xy), (10) i lm where Im is an index set of size jj/ < Ciw?. The integer m is the order of the expansion. The functions $,(y) must be independent of x and the centre y. Furthermore, integrals of $,(y) over single boundary elements F/3 must be easily computable. Usually, 4>i(y) is a (piecewise) polynomial over each element F/?. The error of the expansion of equation (10) depends on the order m and on the distance y y from the centre. The exact requirements on k are as follows. Let 770 G (0,1) be given. Then, there are constants d, d such that for all x,y G R^, all 0 < rj < TJQ < 1, and all m G N there are expansions k of the form (10) satisfying fc(x,y)-fc(x,y) < Ci(C^r I *(x,y) for all y-y < 7? x- This inequality provides an estimate of the relative error of k. Next, "clusters" r, which are unions of several boundary elements F/?, are introduced. For each observation point x, the boundary F can be represented by a certain number of elements and clusters: F = FiUr2U...UFpUriUr2U...Ur, (F/j : elements, TJ : clusters). (11) Since the clusters TJ are unions of many elements, the sum p-\- c can be considerably smaller than the number of elements n. The integral of equation (9) can now be expressed as The first term corresponding to the "near field" FI U Y^ U... U Fp will be evaluated directly. The clusters TJ correspond to the "far field", where the
6 630 Boundary Elements expansion (10) around a centre y = z^ of TJ can be exploited, i.e. we can approximate the integral over TJ by replacing k by k to obtain /fc(x,y)«(y)dr(y) = / K,(x;y )$,(y)u(y)dr(y) **' ^1 < /m ' $,(y)u(y)dr(y). (12) Since the quantities y;.(u) = / *i(y)u(y)dr(y), (j = l,...,c) (13)./TJ are independent of x, they will be computed in the first phase for all indices i and clusters T. Then the evaluation of equation (12) can be performed for all element nodes x", 1 < a < n, by which has Jj/m = 0(rn?) terms independent of the size of the cluster TJ, and JJr.(w) can be shared among different x*. By taking a hierarchy of clusters, the number of all possible clusters (consisting of more than one element) for all x", 1 < a < n can be kept under the number of elements n [5]. The error of approximating the integral /r.fc(x,y)u(y)df(y) can be controlled by keeping the size of the clusters sufficiently small compared to the distance from x. A partition of F by equation (11) is called "admissible" with respect to x when the clusters are sufficiently small in this sense. The admissible covering with the smallest number of members is used. Apparently, the initial phase computations of the integrals for J^ (u] in equation (13) have to be done for every iteration of the GMRES method, since the iteration vector v% and % in the matrix-vector product of equation (8) is updated for every iteration. However, this can be avoided. In our case, constant elements are used. Hence, /' («) = / $,(y)«(y)dr(y) J Tj = «/,/ $,(y)dr(y) where ^ " 1 *'
7 Boundary Elements 631 This means that the integrals defining Jf, (1 < j3 < n) need be computed only once, and «/ (u) can be computed for different iteration vectors up for each GMRES iteration by just performing the linear combination of equation (14). The same technique for general high order elements is discussed in [5]. In order to approximate the kernels to an accuracy which is consistent with the error introduced by discretizing the original boundary integral equation by the boundary element method, the order m for the approximating polynomials in (10) should depend on the mesh (element) size A, and hence on n. Since it is sufficient to take m ~ O(logn), the computational work and memory required for one matrix-vector multiplication approximating (9) for all observation points x",(a = l,...,n), is O(nlog^n) and O(nlog^n), respectively [5]. The evaluation of the right hand side of the system of linear equations (7) can also be done using the panel clustering technique to save memory and computational work. This time, the matrix-vector multiplication need be performed only once for the given vectors wf and jof coming from the boundary data. It is important to note that, when treating domains with complex geometry, the distance along the boundary does not necessarily reflect the distance in 3-D Euclidean space. Therefore, it makes more sense to partition the whole domain into a hierarchy of sub-domains and to cluster them, instead of clustering the boundary elements along the boundary. This will also prove useful when evaluating the displacement or stress at an internal point using the panel (domain) clustering method. 5 A modified expression for the elastostatic fundamental solution In order to obtain the polynomial expansion of (10) to the ra-th order for the elastostatic kernel %^ of (2) and p^ of (3), one could expand ^ to the m- th order and r to the (m + l)-th order and combine them by multiplication and linear combination. Alternatively, as in [7], if we recall the derivation of the fundamental solution for the elastostatic problem using Galerkin's vector[l] or note that where d, again, is the dimension of space (d 3 in our case), (2) gives
8 632 Boundary Elements Here, we denote rjk = Q Q etc. Since the strain corresponding to u*^ is Hooke's law: 4% = 2^^'^ + %W' where A = ^ ^ is the volumetric component, gives 1 Then, the traction p*^ CT^-HJ corresponding to itj^, where nj is the unit outward normal vector at y P, is given by P* = ^ (l^7 <"* + r,^n, + r _,«A ' J^W,') (^ instead of equation (3). Thus, it suffices to expand r(x,y) = x y around the centre of cluster y, and not ^, since the polynomial expansions for u*k and p*^ can be obtained from the linear combination of the partial derivatives r^ and r^&?, respectively. However, in this formulation, it should be noted [4] that, in order to obtain the ra-th order expansion for u^ and p^, one needs to expand r to the (m + 2)-th and (m + 3)-th order, respectively, instead of the (m + l)-th order when expanding the standard expression in (2) and (3). Hence, a careful comparison must be done between the two strategies for the required computational work and memory to obtain the same consistent accuracy. 6 Evaluation of displacement and stress at an internal point The displacement u/(x) and stress <T^(X) at a point in the domain (internal point) can also be computed using the panel clustering technique using and %'(*) = E^ /?%(x,y)dr(y) - IX (3=1 ^p ff=l where (3=1 07T(1 I/)
9 Boundary Elements 633 and V where we have used -j~ = J^- r,,-. Now the panels (or sub-domains) are clustered according to the distance from the internal point x. 7 Polynomial expansion for the panel clustering method Finally, the (ra-th order) polynomial expansion f(y) of r(y) = y x can be obtained by first taking Taylor expansions of r around the centre of the cluster y = (y, 3/2, 2/3)"^, and then expressing them in terms of polynomials of yi instead of hi = %/, y, i.e., f(y) = f(y + h) m-1 E ^,,,,,3(x;y )2/r^^ + ^m(y,h) (17) =0 where h = (h^h^h^.s = 1/1 + i/i + IA, and «(y,h) = The computation of the coefficients «/i^% j ^3 in (17) can be done recurrently as follows. Let r,, = ^-(y ) etc. and r,- = yf %,-, then, r,ija: := --(r,i r,,;w := --(r,^ etc. Once the polynomial expansions (17) of r are obtained, the expansions for r,/ and then, r,/* -> r^m -> u^ and then, r^j -> r^^nf -> P*k can be easily obtained by differentiating by y^ etc. according to (15) and (16). The polynomial expansion for ^ can be obtained similarly to (17) if the standard formulation of (2) and (3) is preferred.
10 634 Boundary Elements Acknowledgement The first author would like to thank Dr. T. Fukui for enlightening remarks. References [1] Brebbia, C.A. and Dominguez, J., Boundary Elements An Introductory Course, Second Edition, Computational Mechanics Publications and McGraw-Hill Book Company, [2] Greengard, L., The Rapid Evaluation of Potential Fields in Particle Systems, MIT Press, Cambridge, [3] Hackbusch, W. Elliptic Differential Equations, Springer- Verlag, 1992, p.207. [4] Hackbusch, W. (private communication), [5] Hackbusch, W. and Nowak, Z.P. On the fast matrix multiplication in the boundary element method by panel clustering, Numerische Mathematik, Vol. 54, pp , [6] Hay ami, K. and Sauter, S. A formulation of the panel clustering method for the three-dimensional elastostatic problem, Proceedings of the JASon BEM, Tokyo, pp , [7] Pierce, A. P. and Napier, J.A.L. A spectral multipole method for efficient solution of large-scale boundary element models in elastostatics, International Journal for Numerical Methods in Engineering, Vol. 38, pp , [8] Rokhlin, V., Rapid solution of integral equations of classical potential theory, Journal of Computational Physics, Vol. 60, pp , [9] Saad, Y. and Schultz, H., GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SI AM Journal on Scientific and Statistical Computing, Vol. 7, No. 3, pp , [10] Sauter, S., The panel clustering method in 3-d BEM, Wave Propagation in Complex Media, G. Papanicolau ed., IMA- Volumes in Mathematics and its Applications, Vol. 96, pp , Springer- Verlag, [11] Yam ad a, Y. and Hay ami, K., A multipole boundary element method for two dimensional elastostatics, in W. Hackbusch and G. Wittum eds., Boundary Elements: Implementation and Analysis of Advanced Algorithms, Proceedings of the 12th GAMM-Seminar Kiel, January 19 to 21, 1996, Notes on Numerical Fluid Mechanics, Vol. 54, Vie weg- Verlag, Braunschweig, Wiesbaden, Germany, pp , 1996
Transactions on Modelling and Simulation vol 19, 1998 WIT Press, ISSN X
Cost estimation of the panel clustering method applied to 3-D elastostatics Ken Hayami* & Stefan A. Sauter^ * Department of Mathematical Engineering and Information Physics, Graduate School of Engineering,
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
Application of the fast multipole method to the 3-D BEM analysis of electron guns T. Nishida and K. Hayami Department of Mathematical Engineering and Information Physics, School of Engineering, University
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
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 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 informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
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 informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationSOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS. Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA
1 SOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA 2 OUTLINE Sparse matrix storage format Basic factorization
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig H 2 -matrix approximation of integral operators by interpolation by Wolfgang Hackbusch and Steen Borm Preprint no.: 04 200 H 2 -Matrix
More informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationAPPLICATION OF THE SPARSE CARDINAL SINE DECOMPOSITION TO 3D STOKES FLOWS
F. Alouges, et al., Int. J. Comp. Meth. and Ep. Meas., Vol. 5, No. 3 (07) 387 394 APPLICATION OF THE SPARSE CARDINAL SINE DECOMPOSITION TO 3D STOKES FLOWS F. ALOUGES, M. AUSSAL, A. LEFEBVRE-LEPOT, F. PIGEONNEAU
More informationPreconditioned inverse iteration and shift-invert Arnoldi method
Preconditioned inverse iteration and shift-invert Arnoldi method Melina Freitag Department of Mathematical Sciences University of Bath CSC Seminar Max-Planck-Institute for Dynamics of Complex Technical
More informationA new fast multipole boundary element method for solving large-scale two-dimensional elastostatic problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2006; 65:863 88 Published online 7 September 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 0.002/nme.474
More informationTopics. The CG Algorithm Algorithmic Options CG s Two Main Convergence Theorems
Topics The CG Algorithm Algorithmic Options CG s Two Main Convergence Theorems What about non-spd systems? Methods requiring small history Methods requiring large history Summary of solvers 1 / 52 Conjugate
More informationA fast multipole boundary element method for 2D multi-domain elastostatic problems based on a dual BIE formulation
Comput Mech (2008) 42:76 773 DOI 0.007/s00466-008-0274-2 ORIGINAL PAPER A fast multipole boundary element method for 2D multi-domain elastostatic problems based on a dual BIE formulation Y. J. Liu Received:
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 HIGH FREQUENCY ACOUSTIC SIMULATIONS VIA FMM ACCELERATED BEM
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 HIGH FREQUENCY ACOUSTIC SIMULATIONS VIA FMM ACCELERATED BEM PACS: 43.20.Fn Gumerov, Nail A.; Duraiswami, Ramani; Fantalgo LLC, 7496
More informationConjugate gradient method. Descent method. Conjugate search direction. Conjugate Gradient Algorithm (294)
Conjugate gradient method Descent method Hestenes, Stiefel 1952 For A N N SPD In exact arithmetic, solves in N steps In real arithmetic No guaranteed stopping Often converges in many fewer than N steps
More informationAPPROXIMATING GAUSSIAN PROCESSES
1 / 23 APPROXIMATING GAUSSIAN PROCESSES WITH H 2 -MATRICES Steffen Börm 1 Jochen Garcke 2 1 Christian-Albrechts-Universität zu Kiel 2 Universität Bonn and Fraunhofer SCAI 2 / 23 OUTLINE 1 GAUSSIAN PROCESSES
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationAlgebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov
More informationThe flexible incomplete LU preconditioner for large nonsymmetric linear systems. Takatoshi Nakamura Takashi Nodera
Research Report KSTS/RR-15/006 The flexible incomplete LU preconditioner for large nonsymmetric linear systems by Takatoshi Nakamura Takashi Nodera Takatoshi Nakamura School of Fundamental Science and
More informationKrylov subspace projection methods
I.1.(a) Krylov subspace projection methods Orthogonal projection technique : framework Let A be an n n complex matrix and K be an m-dimensional subspace of C n. An orthogonal projection technique seeks
More information4.8 Arnoldi Iteration, Krylov Subspaces and GMRES
48 Arnoldi Iteration, Krylov Subspaces and GMRES We start with the problem of using a similarity transformation to convert an n n matrix A to upper Hessenberg form H, ie, A = QHQ, (30) with an appropriate
More informationITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS
ITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS YOUSEF SAAD University of Minnesota PWS PUBLISHING COMPANY I(T)P An International Thomson Publishing Company BOSTON ALBANY BONN CINCINNATI DETROIT LONDON MADRID
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationFast Methods for Integral Equations. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simulation - Lecture 23 Fast Methods for Integral Equations Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Solving Discretized Integral Equations Using
More informationThe Kinematic Equations
The Kinematic Equations David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 0139 September 19, 000 Introduction The kinematic or strain-displacement
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationParallel Iterative Methods for Sparse Linear Systems. H. Martin Bücker Lehrstuhl für Hochleistungsrechnen
Parallel Iterative Methods for Sparse Linear Systems Lehrstuhl für Hochleistungsrechnen www.sc.rwth-aachen.de RWTH Aachen Large and Sparse Small and Dense Outline Problem with Direct Methods Iterative
More informationFEM/FMBEM coupling for acoustic structure interaction and acoustic design sensitivity analysis with sound-absorbing materials
Boundary Elements and Other Mesh Reduction Methods XXXVIII 113 FEM/FMBEM coupling for acoustic structure interaction and acoustic design sensitivity analysis with sound-absorbing materials Y. M. Xu, H.
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationFinding normalized and modularity cuts by spectral clustering. Ljubjana 2010, October
Finding normalized and modularity cuts by spectral clustering Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu Ljubjana 2010, October Outline Find
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccs Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary c 2008,2010
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationFEM and sparse linear system solving
FEM & sparse linear system solving, Lecture 9, Nov 19, 2017 1/36 Lecture 9, Nov 17, 2017: Krylov space methods http://people.inf.ethz.ch/arbenz/fem17 Peter Arbenz Computer Science Department, ETH Zürich
More informationConjugate Gradient Method
Conjugate Gradient Method direct and indirect methods positive definite linear systems Krylov sequence spectral analysis of Krylov sequence preconditioning Prof. S. Boyd, EE364b, Stanford University Three
More informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationNotes on Numerical Fluid Mechanics
Notes on Numerical Fluid Mechanics Volume 31 Parallel Algorithm,s for Partial Differential Equations Edited by Wolfgang Hackbusch al V1eweg Parallel Algorithms for Partial Differential Equations Proceedings
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this
More informationSummary of Iterative Methods for Non-symmetric Linear Equations That Are Related to the Conjugate Gradient (CG) Method
Summary of Iterative Methods for Non-symmetric Linear Equations That Are Related to the Conjugate Gradient (CG) Method Leslie Foster 11-5-2012 We will discuss the FOM (full orthogonalization method), CG,
More informationAn H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages
PIERS ONLINE, VOL. 6, NO. 7, 2010 679 An H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages Haixin Liu and Dan Jiao School of Electrical
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH
ON ORTHOGONAL REDUCTION TO HESSENBERG FORM WITH SMALL BANDWIDTH V. FABER, J. LIESEN, AND P. TICHÝ Abstract. Numerous algorithms in numerical linear algebra are based on the reduction of a given matrix
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
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 informationFast Multipole Methods for Incompressible Flow Simulation
Fast Multipole Methods for Incompressible Flow Simulation Nail A. Gumerov & Ramani Duraiswami Institute for Advanced Computer Studies University of Maryland, College Park Support of NSF awards 0086075
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationTechnische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.4: Krylov Subspace Methods 2 / 51 Krylov Subspace Methods In this chapter we give
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationMaster Thesis Literature Review: Efficiency Improvement for Panel Codes. TU Delft MARIN
Master Thesis Literature Review: Efficiency Improvement for Panel Codes TU Delft MARIN Elisa Ang Yun Mei Student Number: 4420888 1-27-2015 CONTENTS 1 Introduction... 1 1.1 Problem Statement... 1 1.2 Current
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More informationInstitute of Applied Mechanics. Institut für Baumechanik. Preprint Series. Graz University of Technology
Institute of Applied Mechanics Institut für Baumechanik Preprint Series Graz University of Technology An accelerated symmetric time-domain boundary element formulation for elasticity Matthias Messner,
More informationLecture 11: CMSC 878R/AMSC698R. Iterative Methods An introduction. Outline. Inverse, LU decomposition, Cholesky, SVD, etc.
Lecture 11: CMSC 878R/AMSC698R Iterative Methods An introduction Outline Direct Solution of Linear Systems Inverse, LU decomposition, Cholesky, SVD, etc. Iterative methods for linear systems Why? Matrix
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 informationOn the minimum of several random variables
On the minimum of several random variables Yehoram Gordon Alexander Litvak Carsten Schütt Elisabeth Werner Abstract For a given sequence of real numbers a,..., a n we denote the k-th smallest one by k-
More informationAN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
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 informationSolving integral equations on non-smooth domains
Solving integral equations on non-smooth domains Johan Helsing Lund University Talk at the IMA Hot Topics workshop, August 2-5, 2010 σ yy A B C σ yy Acknowledgement The work presented in this talk has
More informationADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS
J. OPERATOR THEORY 44(2000), 243 254 c Copyright by Theta, 2000 ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS DOUGLAS BRIDGES, FRED RICHMAN and PETER SCHUSTER Communicated by William B. Arveson Abstract.
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 19: Computing the SVD; Sparse Linear Systems Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationAn Adaptive Hierarchical Matrix on Point Iterative Poisson Solver
Malaysian Journal of Mathematical Sciences 10(3): 369 382 (2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal An Adaptive Hierarchical Matrix on Point
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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 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 informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationFast iterative BEM for high-frequency scattering
Fast iterative BEM for high-frequency scattering elastodynamic problems In 3D. Marion Darbas LAMFA UMR CNRS 7352, Amiens, France S. Chaillat (POEMS, CNRS-INRIA-ENSTA) and F. Le Louër (LMAC, UTC) Conference
More informationLecture 18 Classical Iterative Methods
Lecture 18 Classical Iterative Methods MIT 18.335J / 6.337J Introduction to Numerical Methods Per-Olof Persson November 14, 2006 1 Iterative Methods for Linear Systems Direct methods for solving Ax = b,
More informationAn accelerated predictor-corrector scheme for 3D crack growth simulations
An accelerated predictor-corrector scheme for 3D crack growth simulations W. Weber 1 and G. Kuhn 2 1,2 1 Institute of Applied Mechanics, University of Erlangen-Nuremberg Egerlandstraße 5, 91058 Erlangen,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationTwo-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues
Two-level multiplicative domain decomposition algorithm for recovering the Lamé coefficient in biological tissues Si Liu and Xiao-Chuan Cai Department of Applied Mathematics, University of Colorado at
More informationResearch Article Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations
Antennas and Propagation Volume 212, Article ID 756259, 9 pages doi:1.1155/212/756259 Research Article Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations Han Guo, Jun
More informationWeaker assumptions for convergence of extended block Kaczmarz and Jacobi projection algorithms
DOI: 10.1515/auom-2017-0004 An. Şt. Univ. Ovidius Constanţa Vol. 25(1),2017, 49 60 Weaker assumptions for convergence of extended block Kaczmarz and Jacobi projection algorithms Doina Carp, Ioana Pomparău,
More informationA Crash Course on Fast Multipole Method (FMM)
A Crash Course on Fast Multipole Method (FMM) Hanliang Guo Kanso Biodynamical System Lab, 2018 April What is FMM? A fast algorithm to study N body interactions (Gravitational, Coulombic interactions )
More informationRetrospectives on My Studies of Solid Mechanics (II)
Retrospectives on My Studies of Solid Mechanics (II) - a new energy method based on the stationary principle of total energy By Tadahiko Kawai + and Etsu Kazama ++ ABSTRACT A new energy method is proposed
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationParallel Numerics, WT 2016/ Iterative Methods for Sparse Linear Systems of Equations. page 1 of 1
Parallel Numerics, WT 2016/2017 5 Iterative Methods for Sparse Linear Systems of Equations page 1 of 1 Contents 1 Introduction 1.1 Computer Science Aspects 1.2 Numerical Problems 1.3 Graphs 1.4 Loop Manipulations
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More information