Singularly Perturbed Partial Differential Equations
|
|
- Roland Reeves
- 5 years ago
- Views:
Transcription
1 WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics, Prague, Czech Republic. bstract. We present a summary of basic numerical methods for solving singularly perturbed partial differential equations. First the effect of the presence of a singular perturbation is introduced, then two different approaches to the elimination of this effect are discussed. Stabilization techniques are demonstrated on several test examples. daptive mesh-refinement methods are taken into account afterwards. Introduction We call a partial differential equation singularly perturbed when its behaviour and properties signifficantly change while one (or more) of its coefficients tends to some critical value. In this paper the singular perturbation will be always caused by the coefficient ε +. The model equation for our purposes will be the scalar convection-diffusion equation: ε u + b u = f in Ω = (,) d, u = u b on Ω. () Here d = or 2, ε > is a constant coefficient of diffusion, b W, (Ω) d is a convective field satisfying the incompressibility condition divb = and f L 2 (Ω) is a given outer force. We also assume that the Dirichlet boundary condition satisfies u b H /2 ( Ω). Since it is usually impossible to find an analytic solution of the problem () for general data, we have to solve it numerically. It is especially difficult when the convection dominates diffusion, i.e. ε b. Then we have a convection dominated problem and the equation becomes singularly perturbed. In this case the solution of () usually possesses interior and boundary layers. These are narrow regions where the solution changes rapidly. If the mesh size is much larger then the width of these regions, the layers cannot be resolved properly and thus spurious (nonphysical) oscillations occur in the numerical solution. To diminish these oscillations one usually uses upwinding or its multi-dimensional analogy SUPG (Streamline upwind Petrov/Galerkin) method. Despite the fact that the SUPG method significantly improves the quality of the discrete solution, it does not remove all spurious oscillations at layers. For this reason another stabilization terms are added to the discrete problem. These are the so-called SOLD (Spurious oscillation at layers diminishing) methods. One can also use some adaptive mesh-refinement algorithm and refine the mesh along layers. However, it does not always bring the desired effect since the mesh width in layer regions should be extremly small. Weak formulation and finite element discretization To employ the finite element method one has to start with a weak formulation of the problem (). Let ũ b H (Ω) be an extension of u b to the whole domain Ω, then the weak formulation of the problem () reads: Find u H (Ω) such that u ũ b H (Ω) and where a(u,v) = (f,v) v H (Ω), (2) a(u,v) = ε( u, v) + (b u,v). (3) Here and throughout this paper (, ) denotes an inner product in L 2 (Ω). Under the assumptions mentioned above the Lax-Milgram theorem gives the existence and uniqueness of the weak solution u H (Ω) [Lamač, 28]. 54
2 LMČ: SINGULRLY PERTURED PRTIL DIFFERENTIL EQUTIONS The starting point of the finite element discretization is the definition of a triangulation T h of the domain Ω. It consists of a finite number of open polygonal elements K (intervals for d = ) and satisfies Ω = K T h K with diam K h, K T h. The intersection of closures of any two different elements K is either empty or a common vertex or a common edge (d = 2). When d = 2 we confine ourselves just to triangular elements with some additional shape-regularity assumption (the minimum angle assumption or bounded ratios of triangle s diameter to its inscribed circle). Now we obtain a finite element discretization of equation (2) replacing the space H (Ω) by a finite-dimensional finite element subspace V h. We say that u h H (Ω) is a discrete solution obtained by the finite element method if u h ũ bh V h and where ũ bh is a finite element interpolate of ũ b. Stabilization methods a(u h,v h ) = (f,v h ) v h V h, (4) Let us start our observations with a D example. Example. Consider equation () with d =, b const. >, f and u b, i.e. εu + bu = in (,), u() = u() =. (5) The analytical solution of the problem (5) is ( u(x) = b x exp(b ε x) exp( b ε ) ). (6) Solving the problem (5) by continuous piecewise linear finite elements leads to the central finite difference scheme and consequently to unwanted spurious oscillations. The usual remedy is using an upwind difference discretization of the convective term. Using the finite element method one can reach the upwind effect by adding an artificial diffusion ε u = 2bh to the diffusion constant ε: ε u i+ 2u i + u i h 2 + b u i u i h ( = ε + bh ) ui+ 2u i + u i 2 h 2 + b u i+ u i. (7) 2h Here and in what follows x i = ih and u i u(x i ), for i =,,...,N and h = /N, N N. Denote α = bh/2ε the local Péclet number and ξ(α) = coth α /α. Then we can obtain a nodally exact solution [Lamač, 28] of the equation (5) by adding the artifical diffusion ε opt = ε u ξ(α) < ε u to the diffusion constant ε in the central difference scheme (see Figure ) Figure. Details of the discrete () and the exact () solution of Example (with b ) at a layer near x = obtained using the finite element method with added artifical diffusion ε u (left) and ε opt (right). nother way how to acquire the stabilization effect is a change of the weighting (test) functions in the discrete problem (4). The standard piecewise linear weighting functions are v i = x x i χ h [xi,x i ] + x i+ x χ h [xi,x i+ ], i =,2,...,N, (8) where χ I denotes a characteristic function of the interval I. 55
3 LMČ: SINGULRLY PERTURED PRTIL DIFFERENTIL EQUTIONS We change these functions by adding some constants on both intervals [x i,x i ] and [x i,x i+ ]. Thus the modified weighting functions have a form (see Figure 2) v i [c L,c R ] = (v i + c L )χ [xi,x i ] + (v i + c R )χ [xi,x i+ ]. (9) If we choose c L = /2 and c R = /2, then the solution obtained using the finite element method is the same as that obtained using the finite difference method with the full upwind discretization of the convective term. If we choose c L = ξ(α)/2 and c R = ξ(α)/2 then we obtain again the nodally exact solution of the equation (5). v i v U i x i+ x i x i x i+ x i x i Figure 2. The standard weighting function v i (left) and the weighting function v Ui = v i [/2, /2] (right) which causes the upwind effect. If we then apply a similar approach of changed weighting functions in higher dimensions (d > ), we obtain the SUPG method [rooks, 982]. Denoting the residual R(u) = ε u + b u f, the SUPG method reads: Find u h H (Ω) such that u h ũ bh V h and a(u h,v h ) + (R(u h ),τb v h ) = (f,v h ) v h V h. () The function τb v h is (in our case) discontinuous and corresponds to the choice of the constants c L,c R in one dimension. The parameter τ is the so-called stabilization parameter and, in D, we obtain the nodally exact solution for τ = h 2bξ(α) [Lamač, 28]. The choice of the stabilization parameter τ L (Ω) is in higher dimensions [John, 27] recommended as: τ K = τ K = µ K h K 2 b ξ(α K) with α K = ν K b h K 2ε. () Here h K is a diameter of the element K in the direction of b and µ K = ν K = /k, where k is the order of approximation of V h K with respect to the norm,k. The SUPG method does not diminish all the spurious oscillations in the discrete solution, particularly at layers. That is why we add another stabilization term to the left hand side of the equation (). These techniques belong to the so-called SOLD (Spurious oscillations at layers diminishing) methods. The most efficient [John, 27] SOLD methods are the do Carmo-Galeão method [Carmo, 99], the urman-ern method [urman, 22] and mainly the Mizukami- Hughes method [Mizukami, 985]. When u h, the do Carmo-Galeão method adds to the left hand side of () the term ( ε u h, v h ), where ε = σ R(u h) 2 and σ = τ max{, b u h u h 2 R(u h ) }. nother SOLD method term was derived by urman and Ern. When b it adds the term ( ε D u h, v h ) to the left hand side of (). Here D = I (b b)/ b 2 is a projection into the plane perpendicular to the vector b and ε = τ b 2 R(u h ) /( b u h + R(u h ) ). While the last two methods add another stabilization term to the equation (), the Mizukami-Hughes method completely changes the weighting functions in (4). It changes the standard piecewise linear weighting functions v i by adding constants Ci K on elements K suppv i : ϕ i = v i + Ci K χ K, v i V h. (2) K supp v i Thus, for each K, we have three constants Ci K corresponding to the vertices of K and it is assumed that their sum vanishes (see Figure 3). The functions ϕ i are then applied to the discrete problem (4). 56
4 LMČ: SINGULRLY PERTURED PRTIL DIFFERENTIL EQUTIONS a) b ) c ) d ) 2/3. v i P i supp v i b.23 2/3 /6 /6 2/3 ϕ i Figure 3. The construction of the weighting function in the Mizukami-Hughes method: for each element K from the support (b) of the standard piecewise linear weighting function v i (a) one computes three constants Cj K coresponding to the vertices of K (c). The constants Ci K corresponding to the vertex P i are then on each element K added to the function v i (d). proper choice of constants Ci K depends not only on the direction of the vector b, but also on the solution u h itself. Thus, as well as the two SOLD methods mentioned above, the Mizukami-Hughes method is nonlinear. etter properties, in particular at layers, has the improved Mizukami-Hughes method (IMH) introduced by Knobloch [Knobloch, 26]. The following example shows a result of an application of both the SUPG and the improved Mizukami- Hughes method in 2D. Example 2. Let us consider the problem () with d = 2, ε = 3, b = (/2, 3/2) T, f and a boundary condition u b = for x = or y.7 and u b = otherwise. If we apply the SUPG method on an isotropic mesh, we obtain a solution with spurious oscillations at layers. s we can see in Figure 4 the improved Mizukami-Hughes method provides a remedy Figure 4. Solutions of Example 2 solved on an isotropic mesh (left, 8 triangles) using the SUPG (center) and the improved Mizukami-Hughes (right) method. different approach that leads to the stabilizing effect is a stabilization via local projection [Roos, 28]. It introduces a projection π h : L 2 (Ω) D h into a second finite element space D h and then defines a stabilization term based on the fluctuation κ h (b u h ) of b u h, where κ h = id L 2 (Ω) π h. The space D h is generally defined on the different triangulation M h in such a way that several assumptions on its dimension and order of approximation are satisfied. In our case of continuous piecewise linear finite elements we define the space D h on the same triangulation as the space V h, but it consists only of piecewise constant functions. The space V h is also changed. It is enriched element by element with bubble functions β K defined as the product of the barycentric coordinates. Then we say that u h H (Ω) is a solution obtained by the local projection stabilization method, if u h ũ bh V h and where the stabilizing term is given by S h (u h,v h ) = a(u h,v h ) + S h (u h,v h ) = (f,v h ) v h V h, (3) K T h τ K (κ h (b u h ),κ h (b v h )) K. (4) 57
5 LMČ: SINGULRLY PERTURED PRTIL DIFFERENTIL EQUTIONS If we choose τ K = τ opt = h 6 (/ξ(α) 3/α) < h 6 on each element K T h, then the obtained solution is in D for constant data (and an equidistant mesh with diam K = h for all K T h ) nodally exact. The result of an application of this metod is depicted in Figure C C C Figure 5. n application of the local projection method to Example. The coefficient τ K is chosen as, τ opt and h/6, respectively. The solution obtained by the local projection method () is compared with the exact solution (C). In order to remove the oscillations bubble functions are omitted afterwards (). daptively refined meshes The strategy of using adaptively refined meshes consists in a localization of those parts of the computational domain Ω where the discrete solution differs the most from the exact solution. In these regions the mesh is refined in order to allow better approximation of the solution. To localize the problematic regions one uses a posteriori error estimators. They usually depend not only on the coefficients of the equation but also on the applied triangulation or computed solution. While local error estimators (indicators) determine regions where the highest error appears, a global error estimator determines the error in the whole computational domain Ω. It can be derived as a sum of local error estimators. Let us assume that we can assign a local error estimator η K to each element K T h and let the global error estimator η be a sum of local error estimators. Then a general mesh-refining algorithm has following steps:. Construction of an initial triangulation T h. 2. Solution of the discrete problem on T hk. 3. Computation of local error estimators η K. 4. Evaluation of a global error estimate (if sufficient then stop). 5. Selection of elements to refine taking η K into account. 6. Construction of a new triangulation, go to 2. very convenient algorithm for the construction of new triangulations (in step 6) is the redgreen algorithm. It divides selected elements in such a way that the shape-regularity (minimum angle) condition is always satisfied, i.e. the angels in each triangle of each triangulation do not degenerate. nother possibility how to costruct new triangulation is moving nodes in the direction of the increasing error. Combining both approaches one can reach satisfactory results Figure 6. Meshes obtained by the red-green algorithm (left) and NGENER (right). They consist of 7833 and 9859 elements, respectively. 58
6 LMČ: SINGULRLY PERTURED PRTIL DIFFERENTIL EQUTIONS The previous pictures show meshes obtained by the red-green algorithm and a program NGENER (see references) applied to Example 2. During the red-green algorithm the a posteriori error estimator introduced by Kunert [Kunert, 22] was used. The program NGENER not only refines the mesh but also moves nodes. Nevertheless the program uses an error estimator which works just with the computed discrete solution and does not depend directly on data of the equation (ε,b,f). Finally we can evaluate the global error estimate (by Kunert) for all types of methods we used while solving the Example 2. The following table confirms our expectations that the global error estimate (computed using the error estimator introduced by Kunert) decreases much more rapidly when mesh-refinement algorithms are applied. Meshes from Figures 4 and 6 were used. Method SUPG IMH RED-GREEN NGENER NGENER+IMH Error estimate Discussion We have presented several efficient stabilization techniques for solving singularly perturbed partial differential equations and compared them on test examples. If we use only stabilization methods we cannot achieve a sufficiently small global error (unless we use large amount of nodes). On the other hand in a strongly convection dominated case we cannot diminish all spurious oscillations at layers using only adaptive mesh-refinement since the layers are too narrow. In higher dimensions it is probably important to construct such triangulations whose elements are thinner in the direction where the oscillations are observed. Due to an added artificial diffusion, stabilization methods can also cause smearing of the solution near layers. Conclusion ll experiments confirmed the fact that singularly perturbed partial differential equations are numerically very difficult to solve. If we want to solve them successfully, we have to appropriately combine stabilization and adaptive methods in order not only to diminish all spurious oscillations at layers but also acquire sufficiently accurate solution. References NGENER 3., mesh generator and mesh adaptor, version 3., rooks,. N., T. J. R. Hughes, Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations, Computer Methods In pplied Mechanics nd Engineering, 32, , 982. urman, E., E. Ern, Nonlinear diffusion and discrete maximum principle for stabilized Galerkin approxiations of the convection-diffusion-reaction equation, Computer Methods In pplied Mechanics nd Engineering, 9, , 22. Carmo, E. G. D. do,. C. Galeão, Feedback Petrov-Galerkin methods for convection-dominated problems, Computer Methods In pplied Mechanics nd Engineering, 88, -6, 99. John, V., P. Knobloch, comparison of spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equation: Part I - review, Computer Methods In pplied Mechanics nd Engineering, 96, , 27. Knobloch, P., Improvements of the Mizukami-Hughes method for convection-diffusion equations, Computer Methods In pplied Mechanics nd Engineering, 96, , 26. Kunert, G., posteriori error estimation for convection-dominated problems on anisotropic meshes, Preprint- Reihe des Chemnitzer Sonderforschungbereich 393/2-4, Technische Universität Chemnitz, March, 22. Lamač, J., Numerical solution of convection-diffusion equations using stabilization and adaptive methods, Master s Thesis, Charles University in Prague, 28. Mizukami,., T. J. R. Hughes, Petrov-Galerkin finite element method for convection-dominated flows: n accurate upwinding technique for satisfying the maximum principle, Computer Methods In pplied Mechanics nd Engineering, 5, 8-93, 985. Roos, H.-G., M. Stynes, L. Tobiska, Robust Numerical Methods for Singularly Perturbed Differential Equations, Springer-Verlag erlin Heidelberg,
NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationOn spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P1 and Q1 finite elements
Volker John, Petr Knobloch On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P and Q finite elements Preprint No. MATH-knm-27/4 27.
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 156 A comparison of spurious oscillations at layers diminishing (SOLD) methods for convection
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationSTABILIZED GALERKIN FINITE ELEMENT METHODS FOR CONVECTION DOMINATED AND INCOMPRESSIBLE FLOW PROBLEMS
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 29 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 1994 STABILIZED GALERIN FINITE ELEMENT METHODS FOR CONVECTION
More informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationSuperconvergence analysis of the SDFEM for elliptic problems with characteristic layers
Superconvergence analysis of the SDFEM for elliptic problems with characteristic layers S. Franz, T. Linß, H.-G. Roos Institut für Numerische Mathematik, Technische Universität Dresden, D-01062, Germany
More informationNumerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems
DOI 10.1007/s10915-009-9337-6 Numerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems Pengtao Sun Long Chen Jinchao Xu Received: 17 April 2009 / Revised:
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 informationNUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE
IJRRAS 5 () November NUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE. K. Sharath Babu,* & N. Srinivasacharyulu Research Scholar, Department of Mathematics, National Institute
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationReview Article Robust Numerical Methods for Singularly Perturbed Differential Equations: A Survey Covering
International Scholarly Research Network ISRN Applied Mathematics Volume 2012, Article ID 379547, 30 pages doi:10.5402/2012/379547 Review Article Robust Numerical Methods for Singularly Perturbed Differential
More informationA stabilized finite element method for the convection dominated diffusion optimal control problem
Applicable Analysis An International Journal ISSN: 0003-6811 (Print) 1563-504X (Online) Journal homepage: http://www.tandfonline.com/loi/gapa20 A stabilized finite element method for the convection dominated
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationA NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES
A NOTE ON CONSTANT-FREE A POSTERIORI ERROR ESTIMATES R. VERFÜRTH Abstract. In this note we look at constant-free a posteriori error estimates from a different perspective. We show that they can be interpreted
More informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationJURJEN DUINTJER TEBBENS
Proceedings of ALGORITMY 005 pp. 40 49 GMRES ACCELERATION ANALYSIS FOR A CONVECTION DIFFUSION MODEL PROBLEM JURJEN DUINTJER TEBBENS Abstract. When we apply the GMRES method [4] to linear systems arising
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationAN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY SINGULARLY PERTURBED PROBLEM
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 11, Number 1, Pages 24 33 c 2014 Institute for Scientific Computing and Information AN OPTIMAL UNIFORM A PRIORI ERROR ESTIMATE FOR AN UNSTEADY
More informationA unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations
Noname manuscript No. (will be inserted by the editor) A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations Gabriel R. Barrenechea Volker John Petr Knobloch Richard
More informationA VOLUME MESH FINITE ELEMENT METHOD FOR PDES ON SURFACES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2012) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 10-14, 2012 A VOLUME MESH FINIE ELEMEN MEHOD FOR
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 informationA Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods
ScientifiComputing A Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods Thomas-Peter Fries, Hermann G. Matthies Institute of Scientific Computing Technical University
More informationFinite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems
Finite Element Superconvergence Approximation for One-Dimensional Singularly Perturbed Problems Zhimin Zhang Department of Mathematics Wayne State University Detroit, MI 48202 Received 19 July 2000; accepted
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationSome remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes
Some remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes Stefano Micheletti Simona Perotto Marco Picasso Abstract In this paper we re-address the anisotropic recipe
More informationNitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem
Nitsche XFEM with Streamline Diffusion Stabilization for a Two Phase Mass Transport Problem Christoph Lehrenfeld and Arnold Reusken Bericht Nr. 333 November 2011 Key words: transport problem, Nitsche method,
More informationGunar Matthies 1, Piotr Skrzypacz 2 and Lutz Tobiska 2. Introduction
A UNIFIED CONVERGENCE ANALYSIS FOR LOCAL PROJECTION STABILISATIONS APPLIED TO THE OSEEN PROBLE Gunar atthies 1, Piotr Skrzypacz 2 and Lutz Tobiska 2 Abstract. The discretisation of the Oseen problem by
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationLINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 445 461 c World Scientific Publishing Company LINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationSUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations
SUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations Swetlana Giere a,, Traian Iliescu b,1, Volker John a,c, David Wells b,1, a Weierstrass Institute for Applied Analysis
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 SUPG reduced order models for convection-dominated convection-diffusion-reaction
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 information^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil
A space-timefiniteelement formulation for the shallow water equations F.L.B. Ribeiro,' A.C. Galeao," L. Landau* "Programa de Engenharia Civil, COPPE / Universidade Federal do Rio de Janeiro, Caixa Postal
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
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 informationInterpolation in h-version finite element spaces
Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationSuperconvergence Using Pointwise Interpolation in Convection-Diffusion Problems
Als Manuskript gedruckt Technische Universität Dresden Herausgeber: Der Rektor Superconvergence Using Pointwise Interpolation in Convection-Diffusion Problems Sebastian Franz MATH-M-04-2012 July 9, 2012
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 informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationConvergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements
Convergence on Layer-adapted Meshes and Anisotropic Interpolation Error Estimates of Non-Standard Higher Order Finite Elements Sebastian Franz a,1, Gunar Matthies b, a Department of Mathematics and Statistics,
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationBenchmarking high order finite element approximations for one-dimensional boundary layer problems
Benchmarking high order finite element approximations for one-dimensional boundary layer problems Marcello Malagù,, Elena Benvenuti, Angelo Simone Department of Engineering, University of Ferrara, Italy
More informationA Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation
University of Colorado at Denver and Health Sciences Center A Petrov-Galerkin Enriched Method: A Mass Conservative Finite Element Method For The Darcy Equation Gabriel R. Barrenechea, Leopoldo P. Franca,
More informationANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION
Proceedings of EQUADIFF 2017 pp. 127 136 ANALYSIS OF THE FEM AND DGM FOR AN ELLIPTIC PROBLEM WITH A NONLINEAR NEWTON BOUNDARY CONDITION MILOSLAV FEISTAUER, ONDŘEJ BARTOŠ, FILIP ROSKOVEC, AND ANNA-MARGARETE
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationSuitability of LPS for Laminar and Turbulent Flow
Suitability of LPS for Laminar and Turbulent Flow Daniel Arndt Helene Dallmann Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics VMS 2015 10th International Workshop on
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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
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 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 informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationNeumann-Type A Posteriori Error Estimation For Steady Convection-Diffusion Equation
INTRNATIONAL JOURNAL OF TCHNOLOGY NHANCMNTS AND MRGING NGINRING RSARCH, VOL 3, ISSU 8 8 Neumann-Type A Posteriori rror stimation For Steady Convection-Diffusion quation G. Temesgen Meuria, J. Anand Rao
More informationMODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS. YIN-TZER SHIH AND HOWARD C. ELMAN y
MODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS YIN-TZER SHIH AND HOWARD C. ELMAN y Abstract. We consider the design of robust and accurate nite element approximation methods for
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationProjection Methods for Rotating Flow
Projection Methods for Rotating Flow Daniel Arndt Gert Lube Georg-August-Universität Göttingen Institute for Numerical and Applied Mathematics IACM - ECCOMAS 2014 Computational Modeling of Turbulent and
More informationStabilized and Coupled FEM/EFG Approximations for Fluid Problems
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM4, Sept. 5-, 24, Beijing, China c 24 Tsinghua University Press & Springer-Verlag Stabilized and Coupled FEM/EFG Approximations for Fluid Problems
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationStabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems
INTERNATIONAL CENTER FOR NUMERICAL METHODS IN ENGINEERING Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems E. Oñate M. Manzan Publication CIMNE Nº-183, February 000
More informationAdaptivity and Variational Stabilization for Convection-Diffusion Equations
Adaptivity and Variational Stabilization for Convection-Diffusion Equations Albert Cohen, Wolfgang Dahmen, Gerrit Welper January 31, 2011 Abstract In this paper we propose and analyze stable variational
More informationA STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION
A STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION Eugenio Oñate International Centre for Numerical Methods in Engineering Universidad Politécnica
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationOn the choice of a Stabilizing Subgrid for Convection-Diffusion Problems
On the choice of a tabilizing ubgrid for Convection-Diffusion Problems F. Brezzi,2, L. D. Marini,2, and A. Russo,2 Abstract UPG and Residual-Free Bubbles are closely related methods that have been used
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More 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 information