An Error Minimized Pseudospectral Penalty Direct Poisson Solver
|
|
- Tabitha Lamb
- 5 years ago
- Views:
Transcription
1 An Error Minimized Pseudospectral Penalt Direct Poisson Solver b Horng Tz-Leng Dept. of Applied Mathematics Feng Chia Universit Taichung Taiwan tlhorng123@gmail.com and Teng Chun-Hao Dept. of Applied Mathematics ational Chiao Tung Universit Hsin-Chu Taiwan chunhao.teng@gmail.com
2 Abstract This paper presents a direct Poisson solver based on an error minimized Chebshev pseudospectral penalt formulation for problems defined on rectangular domains. The current method is based on eigenvalueeigenvector matri diagonalization methods developed long time before and elaborated b Chen Su and Shizgal (2000) [1]. However this kind of fast Poisson solvers utilizing pseudospectral discretization is usuall restricted to periodic or Dirichlet boundar conditions onl for efficient implementation. Enforcement of eumann or Robin boundar conditions requests mess row operations that reduces the easiness of operation. The current method can easil accommodate all three tpes of boundar conditions: Dirichlet eumann and Robin boundar conditions. The reason for that is that it emplos penalt method. Besides the penalt parameters are determined analticall such that the discrete L2 error is minimized. 2D and 3D numerical eperiments are conducted and the results show that the penalt scheme computes numerical solutions with better accurac compared to the traditional approaches with boundar conditions enforced strongl. The current method can be etended to generall linear 2nd order elliptic-tpe partial differential boundar value problems with arbitrar coefficients with onl one restriction that those arbitrar coefficients must be able to be separable to coordinate-variable functions. Under this generalization this method can be surel applied to solve famous Helmholtz equation too. As to computing efficienc the asmptotic operation count for current method in 3D case is 2z ( + + z) and its counterpart for 2D case is 2( + ). Obviousl it is far superior to method epressing Laplace operator in tensor product. For an FFT-based method the asmptotic operation count is basicall 2z (log() + log() + log(z)) for 3D case and 2(log() + log()) for 2D case. Indeed the current method which relies on etensive matri-matri multiplications is inferior to FFT-based methods in theor. However this inferiorit also depends on hardware. For moderate and z it is not necessaril inferior in computers nowadas. Of course for large grid resolution FFTbased methods remain the the best choice if handling boundar conditions is not a problem. [1] H. CHE AD Y. SU AD B.D. SHIZGAL. A Direct Spectral Collocation Poisson Solver in Polar and Clinder Coordinates. J. Comput. Phs. 160 (2000)
3 Approimation b Chebshev polnomials 1 Suppose that u is approimated b a truncated series of Chebshev polnomials as u uˆ ktk k 0 where T denotes Chebshev polnomial and uˆ epansion coefficient. k 2 Choosing the Chebshev Gauss-Lobatto collocations points: j j cos j 01 3 We can form the following discrete transform and inverse transform 2 1 u u T j k k j k j k j c k j0 cj u u T k 0 2 j 0 2 k 0 where cj and ck. 1 otherwise 1 otherwise k
4 4 For pseudospectral method interpolating u at the collocation points j 01. j 5 Hence u can be also epressed as u Lk u k 0 where L is the Lagrange interpolating polnomial defined as k L k1 2 l ( 1) (1 ) T( ) k 2. l0 k l ck ( k) lk 2 du d u 6 Then the derivatives and can be approimated as 2 d d 2 2 du dlk d u d Lk j uk j and 2 juk 2 j d d d d j j dl d derivative matri. k0 k0 k where j is usuall referred to the Chebshev collocation k
5 Chebshev Collocation Derivative Matri 1 The entries of first-order Chebshev collocation derivative matri with respect to based on Gauss-Lobatto collocation points are given as below: D D 6 6 j 2 for i 0 or D j 12 1 where jj 2 ci 21 1 otherwise j i j ci 1 D i j and i j 01 2 ij c j i j j here j cos j 012 are Chebshev Gauss-Lobatto collocation points. Also there is a version for avoiding round-off error when is large. (W. Don and S. Solomonoff in SIAM J. Sci. Comp. Vol. 6 pp (1994)) 2 D D 2.
6 Chebshev Gauss-Lobatto collocation mesh
7 Eample: 1D Poisson equation (PE) with Dirichlet BC s 1 Consider an 1D Poisson equation with Dirichlet boundar conditions as follows u f 11 u 1 g and u 1 g. 2 Discretization: ( The first/last row/column of D can due to Dirichlet boundar conditions.) be stripped u 1 1 f D g D 10 g 1 u f 2 2 D g D 20 g 2 D ij u 2 f 2 D g 20 D g 2 u 1 f 1 D g 10 D g 1 where D D i j ij ij (3) Eponential order of convergence.
8 1 Consider a 2D Poisson problem with Dirichlet boundar conditions in a rectangular domain shown as follows u u f u g u g u g u g 2 Discretizing collocation points in the and directions for u Eample: 2D Poisson equation with Dirichlet BC s respectivel and keeping u in the shape of a rectangular matri: t t t ( +1) ( +1) D u D u D u ud f with u f R. 3 The first and last rows and columns of D and D can be stripped due to Dirichlet BCs then we can obtain D D i j D D i j ij ij ij ij g j D g j D g i D g i ( -1) ( -1) t D u ud f where u f R where f f D ij ij i0 i 1 2 1; j ( ) ( ) ( ) ( ) with i j0 j
9 References for Diagonalization Method (Dirichlet BC s) 1 D. B. Haidvogel Thomas Zang The accurate solution of Poisson's equation b epansion in Chebshev polnomials J. Comput. Phs. 30 (1979) (2D diagonalization under spectral method) (2) P. Haldenwang G. Labrosse and S. Abboudi Chebshev 3-D spectral and 2-D pseudospectral solvers for the Helmholtz equation J. Comput. Phs. 55 (1984) Hwar C. Ku Richard S. Hirsh and Thomas D. Talor A pseudospectral method for solution of the three-dimensional incompressible avier-stokes equations J. Comput. Phs. 10 (1987) (earliest 2D) (4) D. Funaro and D. Gottlieb A new method of imposing boundar conditions in pseudospectral approimations of hperbolic equations Mathematics of Computation 51:184 (1988) (penalt method analsis) 5 U. Ehrenstein R. Peret A Chebshev collocation method for the avier-stokes equations with application to double-diffusive convection International Journal for umerical Methods in Fluids 9 (1989) (2D Helmholtz) 6 Henr H. Yang Bernie Shizgal Chebshev pseudospectral multi-domain technique for viscous flow calculation Comput. Methods Appl. Mech. Engrg. 118 (1994) (2D Helmholtz eumann BC)
10 7 S. Zhao and M. J. Yedlin A new iterative Chebshev spectral method for solving the elliptic equation u f J. Comput. Phs. 113 (1994) (2D) (8) H. B. Chen K. andakumar W. H. Finla and H. C. Ku Three-dimensional viscous flow through a rotating channel: a pseudospectral matri method approach 23 (1996): (3D Helmholtz Dirichlet BC) 9 Y. Su Collocation spectral methods in the solution of Poisson equation MS thesis Universit of British Columbia Canada (deatiled review and 2D derivation) 10 Heli Chen Yuhong Su and B. D. Shizgal A direct spectral collocation Poisson solver in polar and clindrical coordinates J. Comput. Phs. 160 (2000) (2D/3D clindrical coordinate and smbolic derivation) 11 Roger Peret Spectral methods for incompressible viscous flow Springer-Verlag ew York (2D p 88) 12 S. guen C. Delcarte A spectral collocation method to solve Helmholtz problems with boundar conditions involving mied tangential and normal derivatives J. Comput. Phs. 200 (2004) (complicated Helmholtz solver single/multi-block)
11 Diagonalization via Eigenvector De where i ie D i Y Y 1 here and e represent the eigenvalue and eigenvector of D and Y e e. 1 e Similarl i t t D e e D X X j i j j 1 2 where. 1
12 D u ud f 1 DYY ux ux fx 1 YY ux ux fx Y ux Y ux Y fx where and 1 1 V V B V Y ux B Y fx t D ux ud X fx t D ux ux fx i.e. v v b i and j v ij i ij ij j ij i b ij j where v ij and b are the elements of V and B respectivel. ij 1 Once V is solved u can be recovered from u YVX. The diagonalization method so far onl works easil for Dirichlet BC's. For non-dirichlet BC's such as eumann and Robin BC's it requests mess row operations.
13 Eample: 1D Poisson equation with BC s involving derivative 1 Consider an 1D Poisson equation with Dirichlet/eumann boundar conditions as follows u f 11 u 1 g and u 1 g. 2 Stripping off first/last row/column of (3) Discretization: D no longer works. u g 0 u 1 f D 1 0 j D ˆ where Dˆ D ij ij u 1 f 1 I j u g 1 1 here I j stands for the last row of the 1 1 identit matri.
14 Eample: 2D Poisson equation with BC s involving derivative 1 Consider a 2D Poisson equation with Dirichlet/eumann boundar conditions as follows u u f u g u g u g u g 2 Following the idea of 1D case discretizing collocation points in and directions respectivel for u and keeping u in the form of a rectangular matri: D u ud fˆ u fˆ R ˆ ˆ t ( +1) ( +1) where Dˆ fˆ D 0 j D 0 j D Dˆ D ij ij I I j j ( +1) ( +1) g or g g g g or g g g fij g g g or g g g g or g ( +1) ( + 1) ( +1) ( +1) i here. j 1 2 1
15 ˆ ˆ t 3 However using 55 matri as an eample D u ud fˆ will become u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u uu u u u u u u uu g or g g g g g or g g f f f g g f f f g g f f f g g or g g g g g or g which is wrong at BC's. How to fi this? Resorting to Kronecker product I Dˆ Dˆ I is both time and memor consuming.
16 Diagonalization with Penalt Method 1 Consider the 2D Poisson equation with all kind of boundar conditions in a rectangular domain i.e. where Bu 1 g ( ) Bu 1 g ( ) B and and are specified constants. At least one of and B must be non-zero for a unique solution to eist. 2 At this time we can not strip first/last row/column of D and D as we did for the case of Dirichlet BCs. We follow the previous idea but modif the substitution b penalt method shown as belows.
17 3 Replace D and D b the following D and D respectivel D D D I0 j D j D D I j D 0 0 j ij D j j 0 j I 0 j D I j D j j D 0 j where and are the penalt weights which are usuall large numbers. ij D
18 t 4 Using 55 matri as an eample D uud f will become u Bu Bu u Bu u Bu u Bu u Bu Bu ubu u u u ubu ubu u u u ubu u Bu u u u u Bu u Bu Bu u Bu u Bu u Bu u Bu Bu f g g f g f g f g f g g f g f f f f g f g f f f f g f g f f f f g f g g f g f g f g f g g which well approimates u f with the specified BC's when 's are large enough.
19 Diagonalization via Eigenvector again D e where i ie D i Y Y here and e represent the eigenvalue and eigenvector of D and Y e e. 0 e 1 2 Similarl j i t t j j i D e e D X X 0 1 where.
20 D uud f 1 DYY ux ux fx 1 YY ux ux fx Y ux Y ux Y fx where and 1 1 V V B V Y ux B Y fx t D ux ud X fx t D ux ux fx i.e. v v b i 01 and j 01 v ij i ij j ij ij i b ij respectivel. j where v and b are the elements of V and B ij 1 Once is solved can be recovered from. V u u YVX ij
21 1 D D 3D Case D 0 j I j D D D I j D 0 0 j ij D j j D 0 j I j D D ij D I j D 0 0 j j j z z z Dzz 0 j I j Dz D z z z Dzz I z j Dz 0 0 j zz zz ij j j z z z z
22 (2) D Y Y Y D Y 1 il lm m im mi il lm m l0 i0 l0 D X X X D X t 1 t lj jn n ln nl lj jn n j0 l0 j0 D Z Z Z D Z. z z z t 1 t zz lk ko o lo ol zz lk ko o k0 l0 k0
23 (3) z t t il ljk ilk lj ijl zz lk ijk l0 l0 l0 D u u D u D F z z D u X Z u X Z il ljk jn ko ilk n ln ko k0 j0 l0 k0 l0 z z u Z X F X Z ijl o lo jn ijk jn ko j0 l0 k0 j0 z D z 1 Y il lmyml uljk X jnzko uilk XlnZko n k0 j0 l0 m0 k0 l0 z u X Z ijl jn lo o ijk jn ko j0 l0 k0 j0 z F X Z
24 z z 1 YimmYml uljk X jnzko uilk XlnZkon k0 j0 l0 k0 l0 z u X Z ijl jn lo o ijk jn ko j0 l0 k0 j0 z F X Z z z 1 1 myml uljk X jnzko Ymi uilk XlnZkon k0 j0 l0 i0 k0 l0 z z 1 1 Ymi uijl X jnzloo Ymi Fijk X jnzko i0 l0 j0 i0 k0 j0
25 z z 1 1 m Yml uljk X jnzko n Ymi uilk XlnZko i0 k0 j0 i0 k0 l0 z z 1 1 o Ymi uijl X jnzlo Ymi Fijk X jnzko i0 l0 j0 i0 k0 j0 Hence W B W m n o mno mno mno z z 1 1 where Wmno Ymi uilk XlnZko Bmno Ymi Fijk X jnzko. i0 k0 l0 i0 k0 j0 m B mno n o mno z 1 1 ijk im mno nj ok m0 o0 n0 u Y W X Z.
26 Error Minimization for Penalt Parameter 2 1 Taking 1D PE as eample u( ) f( ). 2 f( ) 16 sin 4 so eact solution u( ) sin(4 ). 3 BC's: (a) u( 1) 0 u(1) 4 (b) u( 1) 4 u(1) u(1) 4. Figure 1: BC (a) =16. Figure (2): BC (b) =16.
27 Error-Minimized Tau vs. Strongl Enforcement Strongl-enforced BC is equivalen to. Table 1: BC (a) Table 2: BC (b)
28 Error-Minimized Tau Can Be Derived Details see: Tz-Leng Horng and Chun-Hao Teng* 2012 "An error minimized pseudospectral penalt direct Poisson solver" Journal of Computational Phsics 231(6):
29 Error-Minimized Tau Can Be Etended to Higher-Dimension PE Taking 2D PE as an eample u u u f( ) with BC's Bu 1 g ( ) Bu 1 g ( ). There eists a unique set u( ) u( ) f ( ) f ( ) such that u ( ) u( ) u( ) and f( ) f ( ) f ( ) with u f subject to B u 1 g ( ) u f subject to B u 1 g ( ) So we can calculate error-minimized and in each direction alone according to previous formula.
30 umerical eperiments for 2D PE 1 u f( ). 2 2 f( ) 32 sin(4 )sin(4 ) so eact solution is u ( ) sin(4 )sin(4 ). u( 1 ) (3) BC's: (a) u( 1 ) 4 sin(4 ) u ( 1) 4sin(4 ).
31 BC's: (b) u( 1 ) 0 u( 1) 0. CT (Chebshev tau): D.B. Haidvogel T. Zang The accurate solution of Poisson s equation b epansion in Chebshev polnomials J. Comput. Phs. 30 (1979) CC (Chebshev collocation): U. Ehrenstein R. Peret A Chebshev collocation method for the avier Stokes equations with application to double-diffusive convection Int. J. umer. Methods fluids 9 (1989) CG (Chebshev Galerkin): J. Shen Efficient spectral-galerkin method II. Direct solvers of second and fourth order equations b using Chebshev polnomials SIAM J. Sci. Comput. 16 (1995)
32 u(1 ) BC's: (c) u1 0 4sin(4 ) u ( 1) u ( 1) 0 u ( 1) 4sin(4 ).
33 umerical eperiment for 3D PE 1 u f( z). 2 2 f( z) 48 sin(4 )sin(4 )sin(4 z) so eact solution is u ( ) sin(4 )sin(4 )sin(4 z). u(1 z) (3) BC's: u( 1 z) 0 4sin(4 )sin(4 z) u ( 1 z) u ( 1 z) 4sin(4 )sin(4 z) u ( 1 z) 4sin(4 )sin(4 z) u ( 1) u ( 1) u ( 1) 4 sin(4 )sin(4 ) 4sin(4 )sin(4 ). z z
34 1 Asmptotic operation account: Conclusion 2D: 2 3D: 2 compared with FFT-based method z z 2D: 2 log log 3D: 2 log log log. z z (2) This method can be easil etended to the following tpe of elliptic problem subjected to all kinds of BC's in a rectangular domain: where L a L L Lzu f z b c L a b c Lz a3 z b 2 3 z c3 z z z Sure it includes Helmholtz equation. (3) So far analtic formula for error-minimized tau can be onl derived for PE.. 4 For all BC's being eumann in PE which has infinitel man solutions and the difference of an two of them will be a constant the current method will give ou one of the solutions. 5 The basis function is not limited to Chebshev polnomials. It can be replaced b Legendre polnomials or an other basis function as far as there is an associated collocation derivative matri. It can be applied to finite difference too. (6) Sample matlab code is available at
35 Thank ou for our attention. An question?
A FAST POISSON SOLVER BY CHEBYSHEV PSEUDOSPECTRAL METHOD USING REFLEXIVE DECOMPOSITION. Teng-Yao Kuo, Hsin-Chu Chen and Tzyy-Leng Horng*
TAWANESE JOURNAL OF MATHEMATCS Vol. 17, No. 4, pp. 1167-1181, August 213 DO: 1.1165/tjm.17.213.2574 This paper is available online at http://journal.taiwanmathsoc.org.tw A FAST POSSON SOLVER BY CHEBYSHEV
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationKrylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationFast Direct Solver for Poisson Equation in a 2D Elliptical Domain
Fast Direct Solver for Poisson Equation in a 2D Elliptical Domain Ming-Chih Lai Department of Applied Mathematics National Chiao Tung University 1001, Ta Hsueh Road, Hsinchu 30050 Taiwan Received 14 October
More informationAn improved ADI-DQM based on Bernstein polynomial for solving two-dimensional convection-diffusion equations
Mathematical Theor and Modeling ISSN -50 (Paper) ISSN 5-05(Online) Vol 3, No.1, 013 An improved ADI-DQM based on Bernstein polnomial for solving two-dimensional convection-diffusion equations A.S.J. Al-
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationCOMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER PARABOLIC EQUATIONS. Lili Ju. Xinfeng Liu.
DISCRETE AND CONTINUOUS doi:13934/dcdsb214191667 DYNAMICAL SYSTEMS SERIES B Volume 19, Number 6, August 214 pp 1667 1687 COMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
More informationApply mass and momentum conservation to a differential control volume. Simple classical solutions of NS equations
Module 5: Navier-Stokes Equations: Appl mass and momentum conservation to a differential control volume Derive the Navier-Stokes equations Simple classical solutions of NS equations Use dimensional analsis
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationCourse Requirements. Course Mechanics. Projects & Exams. Homework. Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications
Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov What are multipole methods and what is this course about. Problems from phsics, mathematics,
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationMAE 82 Engineering Mathematics
Class otes : First Order Differential Equation on Linear AE 8 Engineering athematics Universe on Linear umerical Linearization on Linear Special Cases Analtical o General Solution Linear Analtical General
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 978 994 c 999 Societ for Industrial and Applied Mathematics A STUDY OF MONITOR FUNCTIONS FOR TWO-DIMENSIONAL ADAPTIVE MESH GENERATION WEIMING CAO, WEIZHANG HUANG,
More informationEfficient Chebyshev Spectral Method for Solving Linear Elliptic PDEs Using Quasi-Inverse Technique
Numer. Math. Theor. Meth. Appl. Vol. 4, No. 2, pp. 197-215 doi: 1.428/nmtma.211.42s.5 Ma 211 Efficient Chebshev Spectral Method for Solving Linear Elliptic PDEs Using Quasi-Inverse Technique Fei Liu, Xingde
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 informationnorges teknisk-naturvitenskapelige universitet Semi-Lagrangian exponential integrators for the incompressible Navier-Stokes equations
norges teknisk-naturvitenskapelige universitet Semi-Lagrangian eponential integrators for the incompressible Navier-Stokes equations b Elena Celledoni, Bawfeh Kingsle Kometa and Olivier Verdier preprint
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationREVISION SHEET FP2 (MEI) CALCULUS. x x 0.5. x x 1.5. π π. Standard Calculus of Inverse Trig and Hyperbolic Trig Functions = + = + arcsin x = +
the Further Mathematics network www.fmnetwork.org.uk V 07 REVISION SHEET FP (MEI) CALCULUS The main ideas are: Calculus using inverse trig functions & hperbolic trig functions and their inverses. Maclaurin
More informationConsider a slender rod, fixed at one end and stretched, as illustrated in Fig ; the original position of the rod is shown dotted.
4.1 Strain If an object is placed on a table and then the table is moved, each material particle moves in space. The particles undergo a displacement. The particles have moved in space as a rigid bod.
More information4 Strain true strain engineering strain plane strain strain transformation formulae
4 Strain The concept of strain is introduced in this Chapter. The approimation to the true strain of the engineering strain is discussed. The practical case of two dimensional plane strain is discussed,
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationAccelerator Physics Statistical and Beam-Beam Effects. G. A. Krafft Old Dominion University Jefferson Lab Lecture 14
Accelerator Phsics Statistical and Beam-Beam Effects G. A. Krafft Old Dominion Universit Jefferson Lab Lecture 4 Graduate Accelerator Phsics Fall 7 Waterbag Distribution Lemons and Thode were first to
More informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More informationMethods for Advanced Mathematics (C3) Coursework Numerical Methods
Woodhouse College 0 Page Introduction... 3 Terminolog... 3 Activit... 4 Wh use numerical methods?... Change of sign... Activit... 6 Interval Bisection... 7 Decimal Search... 8 Coursework Requirements on
More informationA Hybrid Numerical Method for High Contrast. Conductivity Problems. Liliana Borcea, and George C. Papanicolaou y. June 10, 1997.
A Hbrid Numerical Method for High Contrast Conductivit Problems Liliana Borcea, and George C. Papanicolaou June, 997 Abstract We introduce and test etensivel a numerical method for computing ecientl current
More information67. (a) Use a computer algebra system to find the partial fraction CAS. 68. (a) Find the partial fraction decomposition of the function CAS
SECTION 7.5 STRATEGY FOR INTEGRATION 483 6. 2 sin 2 2 cos CAS 67. (a) Use a computer algebra sstem to find the partial fraction decomposition of the function 62 63 Find the area of the region under the
More informationMath 214 Spring problem set (a) Consider these two first order equations. (I) dy dx = x + 1 dy
Math 4 Spring 08 problem set. (a) Consider these two first order equations. (I) d d = + d (II) d = Below are four direction fields. Match the differential equations above to their direction fields. Provide
More informationA high order Padé ADI method for unsteady convection-diffusion equations
Center for Turbulence Research Annual Research Briefs 2005 85 A high order Padé ADI method for unstead convection-diffusion equations B D. You 1. Motivation and objectives The unstead convection-diffusion
More informationSTABILIZED FEM SOLUTIONS OF MHD EQUATIONS AROUND A SOLID AND INSIDE A CONDUCTING MEDIUM
Available online: March 09, 2018 Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. Volume 68, Number 1, Pages 197 208 (2019) DOI: 10.1501/Commua1_0000000901 ISSN 1303 5991 http://communications.science.ankara.edu.tr/index.php?series=a1
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More informationTriple Integrals in Cartesian Coordinates. Triple Integrals in Cylindrical Coordinates. Triple Integrals in Spherical Coordinates
Chapter 3 Multiple Integral 3. Double Integrals 3. Iterated Integrals 3.3 Double Integrals in Polar Coordinates 3.4 Triple Integrals Triple Integrals in Cartesian Coordinates Triple Integrals in Clindrical
More information2: Distributions of Several Variables, Error Propagation
: Distributions of Several Variables, Error Propagation Distribution of several variables. variables The joint probabilit distribution function of two variables and can be genericall written f(, with the
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationPart D. Complex Analysis
Part D. Comple Analsis Chapter 3. Comple Numbers and Functions. Man engineering problems ma be treated and solved b using comple numbers and comple functions. First, comple numbers and the comple plane
More informationA DLM/FD Method for Simulating Particle Motion in Viscoelastic Fluid
A Method for Simulating Particle Motion in Viscoelastic Fluid Tsorng-Wha Pan and Roland Glowinski Department of Mathematics Universit of Houston Houston, Teas Ma 3, 6 International Workshop on Fluid-Structure
More informationTime-Frequency Analysis: Fourier Transforms and Wavelets
Chapter 4 Time-Frequenc Analsis: Fourier Transforms and Wavelets 4. Basics of Fourier Series 4.. Introduction Joseph Fourier (768-83) who gave his name to Fourier series, was not the first to use Fourier
More informationA Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 A Meshless Simulations for D Nonlinear Reaction-diffusion Brusselator Sstem Ahmad Shirzadi, Vladimir Sladek, Jan Sladek 3 Abstract: This paper
More informationNotes 7 Analytic Continuation
ECE 6382 Fall 27 David R. Jackson Notes 7 Analtic Continuation Notes are from D. R. Wilton, Dept. of ECE Analtic Continuation of Functions We define analtic continuation as the process of continuing a
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More information8.3. Integration of Rational Functions by Partial Fractions. 570 Chapter 8: Techniques of Integration
570 Chapter 8: Techniques of Integration 8.3 Integration of Rational Functions b Partial Fractions This section shows how to epress a rational function (a quotient of polnomials) as a sum of simpler fractions,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationDATA STRUCTURES FOR THE NONLINEAR SUBGRID STABILIZATION METHOD
DATA STRUCTURES FOR TE NONLINEAR SUBGRID STABILIZATION METOD Isaac P. Santos Lucia Catabriga Departamento de Informática Universidade Federal do Espírito Santo - UFES Av. Fernando Ferrari, 54 - Vitória,
More informationExact Differential Equations. The general solution of the equation is f x, y C. If f has continuous second partials, then M y 2 f
APPENDIX C Additional Topics in Differential Equations APPENDIX C. Eact First-Order Equations Eact Differential Equations Integrating Factors Eact Differential Equations In Chapter 6, ou studied applications
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 informationClosed form expressions for the gravitational inner multipole moments of homogeneous elementary solids
Closed form epressions for the gravitational inner multipole moments of homogeneous elementar solids Julian Stirling 1,2, and Stephan Schlamminger 1 1 National Institute of Standards and Technolog, 1 Bureau
More informationVibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions
Tamkang Journal of Science and Engineering, Vol. 6, No. 4, pp. 7-6 (003) 7 Vibration Analsis of Isotropic and Orthotropic Plates with Mied Boundar Conditions Ming-Hung Hsu epartment of Electronic Engineering
More information15. Eigenvalues, Eigenvectors
5 Eigenvalues, Eigenvectors Matri of a Linear Transformation Consider a linear ( transformation ) L : a b R 2 R 2 Suppose we know that L and L Then c d because of linearit, we can determine what L does
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationCoordinate geometry. + bx + c. Vertical asymptote. Sketch graphs of hyperbolas (including asymptotic behaviour) from the general
A Sketch graphs of = a m b n c where m = or and n = or B Reciprocal graphs C Graphs of circles and ellipses D Graphs of hperbolas E Partial fractions F Sketch graphs using partial fractions Coordinate
More informationPhysically Based Rendering ( ) Geometry and Transformations
Phsicall Based Rendering (6.657) Geometr and Transformations 3D Point Specifies a location Origin 3D Point Specifies a location Represented b three coordinates Infinitel small class Point3D { public: Coordinate
More informationON THE CONVERGENCE OF THE FOURIER APPROXIMATION FOR EIGENVALUES AND EIGENFUNCTIONS OF DISCONTINUOUS PROBLEMS
SIAM J. UMER. AAL. Vol. 40, o. 6, pp. 54 69 c 003 Society for Industrial and Applied Mathematics O THE COVERGECE OF THE FOURIER APPROXIMATIO FOR EIGEVALUES AD EIGEFUCTIOS OF DISCOTIUOUS PROBLEMS M. S.
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationParticular Solutions
Particular Solutions Our eamples so far in this section have involved some constant of integration, K. We now move on to see particular solutions, where we know some boundar conditions and we substitute
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationPolynomial approximation and Splines
Polnomial approimation and Splines 1. Weierstrass approimation theorem The basic question we ll look at toda is how to approimate a complicated function f() with a simpler function P () f() P () for eample,
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More informationAnswers Investigation 3
Answers Investigation Applications. a., b. s = n c. The numbers seem to be increasing b a greater amount each time. The square number increases b consecutive odd integers:,, 7,, c X X=. a.,,, b., X 7 X=
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationMa 227 Final Exam Solutions 5/8/03
Ma 7 Final Eam Solutions 5/8/3 Name: Lecture Section: I pledge m honor that I have abided b the Stevens Honor Sstem. ID: Directions: Answer all questions. The point value of each problem is indicated.
More informationElliptic Equations. Chapter Definitions. Contents. 4.2 Properties of Laplace s and Poisson s Equations
5 4. Properties of Laplace s and Poisson s Equations Chapter 4 Elliptic Equations Contents. Neumann conditions the normal derivative, / = n u is prescribed on the boundar second BP. In this case we have
More informationAffine transformations. Brian Curless CSE 557 Fall 2014
Affine transformations Brian Curless CSE 557 Fall 2014 1 Reading Required: Shirle, Sec. 2.4, 2.7 Shirle, Ch. 5.1-5.3 Shirle, Ch. 6 Further reading: Fole, et al, Chapter 5.1-5.5. David F. Rogers and J.
More informationRoberto s Notes on Integral Calculus Chapter 3: Basics of differential equations Section 3. Separable ODE s
Roberto s Notes on Integral Calculus Chapter 3: Basics of differential equations Section 3 Separable ODE s What ou need to know alread: What an ODE is and how to solve an eponential ODE. What ou can learn
More informationAffine transformations
Reading Optional reading: Affine transformations Brian Curless CSE 557 Autumn 207 Angel and Shreiner: 3., 3.7-3. Marschner and Shirle: 2.3, 2.4.-2.4.4, 6..-6..4, 6.2., 6.3 Further reading: Angel, the rest
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More informationAdditional Topics in Differential Equations
6 Additional Topics in Differential Equations 6. Eact First-Order Equations 6. Second-Order Homogeneous Linear Equations 6.3 Second-Order Nonhomogeneous Linear Equations 6.4 Series Solutions of Differential
More informationSpectral method for the unsteady incompressible Navier Stokes equations in gauge formulation
Report no. 04/09 Spectral method for the unsteady incompressible Navier Stokes equations in gauge formulation T. W. Tee I. J. Sobey A spectral method which uses a gauge method, as opposed to a projection
More informationDISCONTINUOUS GALERKIN DISCRETIZATION WITH EMBEDDED BOUNDARY CONDITIONS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol.3(23), No., pp.35 58 c 23 Editorial Board of the Journal Computational Methods in Applied Mathematics Joint Co. Ltd. DISCONTINUOUS GALERKIN DISCRETIZATION
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More informationExponential time differencing methods for nonlinear PDEs
XXI Congreso de Ecuaciones Diferenciales Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, 1- septiembre 9 (pp. 1 8) Exponential time differencing methods for nonlinear PDEs F. de la Hoz 1,
More informationThe Control-Volume Finite-Difference Approximation to the Diffusion Equation
The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering gerr@mepdedu ME 448/548: D
More informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More informationOn the Missing Modes When Using the Exact Frequency Relationship between Kirchhoff and Mindlin Plates
On the Missing Modes When Using the Eact Frequenc Relationship between Kirchhoff and Mindlin Plates C.W. Lim 1,*, Z.R. Li 1, Y. Xiang, G.W. Wei 3 and C.M. Wang 4 1 Department of Building and Construction,
More informationFall Exam II. Wed. Nov. 9, 2005
Fall 005 10.34 Eam II. Wed. Nov. 9 005 (SOLUTION) Read through the entire eam before beginning work and budget your time. You are researching drug eluting stents and wish to understand better how the drug
More informationA posteriori error estimates have been an integral part of ever adaptive method, thus for the DGM to be used in an adaptive setting one needs to devel
A Posteriori Discontinuous Finite Element Error Estimation for Two-Dimensional Hperbolic Problems Slimane Adjerid and Thomas C. Masse Department of Mathematics Virginia Poltechnic Institute and State Universit
More informationExact Equations. M(x,y) + N(x,y) y = 0, M(x,y) dx + N(x,y) dy = 0. M(x,y) + N(x,y) y = 0
Eact Equations An eact equation is a first order differential equation that can be written in the form M(, + N(,, provided that there eists a function ψ(, such that = M (, and N(, = Note : Often the equation
More informationMathematics. Section B. CBSE XII Mathematics 2012 Solution (SET 2) General Instructions:
CBSE XII Mathematics 0 Solution (SET ) General Instructions: Mathematics (i) (ii) (iii) (iv) (v) All questions are compulsor. The question paper consists of 9 questions divided into three Sections A, B
More informationThe Basic (Physical) Problem: Experiments. Internal Gravity Waves and Hyperbolic Boundary-Value Problems. Internal Gravity Waves. Governing Equations
The Basic (Phsical) Problem: Eperiments Internal Gravit Waves and Hperbolic Boundar-Value Problems P. A. Martin Colorado chool of Mines Collaborator: tefan Llewelln mith, UC an Diego D. E. Mowbra & B..
More informationSimplified Fast Multipole Methods for Micromagnetic Modeling
implified Fast Multipole Methods for Micromagnetic Modeling D. M. Apalkov and P. B.Visscher Department of Phsics and Astronom and MIT Center The Universit of Alabama This project was supported b F grants
More informationNumerical study of some nonlinear wave equations via Chebyshev collocation method
IJST (203) 37A4: 477-482 Iranian Journal of Science & Technology http://ijsts.shirazu.ac.ir Numerical study of some nonlinear wave equations via Chebyshev collocation method N. Y. A. Elazem and A. Ebaid*
More informationLecture 1: Course Introduction.
Lecture : Course Introduction. What is the Finite Element Method (FEM)? a numerical method for solving problems of engineering and mathematical physics. (Logan Pg. #). In MECH 40 we are concerned with
More informationD u f f x h f y k. Applying this theorem a second time, we have. f xx h f yx k h f xy h f yy k k. f xx h 2 2 f xy hk f yy k 2
93 CHAPTER 4 PARTIAL DERIVATIVES We close this section b giving a proof of the first part of the Second Derivatives Test. Part (b) has a similar proof. PROOF OF THEOREM 3, PART (A) We compute the second-order
More information