PROGRESS ON PRACTICAL METHODS OF ERROR ESTIMATION FOR ENGINEERING CALCULATIONS
|
|
- Josephine Newton
- 5 years ago
- Views:
Transcription
1 ECCM-2001 European Conference on Computational Mechanics June 26-29, 2001 PROGRESS ON PRACTICAL METHODS OF ERROR ESTIMATION FOR ENGINEERING CALCULATIONS J. Tinsley Oden and Serge Prudhomme, The Texas Institute for Computational and Applied Mathamatics The University of Texas at Austin, Austin, TX 78712, U.S.A. Tim Westermann and Jon Bass Altair Engineering Inc. Austin, Texas, U.S.A. Key words: Goal-oriented error estimation, Displacement error, Error in stress, Upper and lower bounds, Dual problem, Residual method. Abstract. In the present paper, the goal-oriented error estimation method is applied to several complex engineering analysis of three-dimensional elastic solids. In the examples considered, the errors are estimated with respect to displacement or stress components at a given point in the solid. The numerical results indicate that the method can be used effectively for complex engineering applications.
2 J.T. Oden, S. Prudhomme, T. Westermann, J. Bass 1 Introduction The use of a posteriori error estimation as a means for verification of accuracies of computer simulations has been a topic of research for nearly two decades. Until recently, such error estimates have been global in nature, giving bounds on global approximation errors in energy-type norms. Extensions of global estimation methods to procedures for estimating errors in local quantities of interest has been described in recent papers [I, 2]. However, applications of these new techniques have been confined to very simple model problems involving scalar-valued functions on one- and two-dimensional domains. In the present investigation, we review the theory of goal-oriented error estimation and discuss applications to several complex engineering analysis of three-dimensional elastic solids. These examples indicate that the goal-oriented methods can be used effectively in quite complex engineering applications. 2 Model problem and approximation error Let 0 be a bounded open domain in lr d, d = 2,3, with Lipschitz continuous boundary ao. The region D is occupied by a linearly elastic material body in static equilibrium under the action of body forces f E (L 2 (D))d and surface tractions g E (L 2 (r N ))d, where r N C ao. The portion of the boundary on which displacements are prescribed is denoted r D, such that rdurn = ao, and rdnrn = 0. The stress 0', related to the displacement u by the constitutive equation (T = ~E : (Vu + Vu T ), satisfies inside the linearly elastic the momentum equations: -V u = f, in D. (1) Here E is the elasticity tensor, which is assumed to satisfy the standard symmetry and uniform ellipticity conditions. Moreover, the displacement and stress are subjected to the boundary conditions: u = 0, { (2) u n=g, Let s be a positive integer and 8 an open region in jrd (8 may stand for the whole domain 0, for a subdomain of 0 or for a portion of the boundary). Denote by HS(8) the Sobolev space of functions in L2(8) with distributional derivatives of order ~ s in L2(8) as well. The space of admissible displacements V is therefore defined as: (3) the boundary values being understood in the sense of traces of Hl-functions. The classical variational form associated with the elasticity problem consists in finding the displacement u E V such that: B(u, v) = F(v), VVEV, (4) 2
3 ECCM-2001, where B(u, v) F(v) = InVv : EVu dx, 1 f.v d+1g' v ds. n fn The bilinear form B(,.) is symmetric positive-definite on Y x Y and induces a norm, the so-called energy norm, Ilvlle =.jb(v, v). (6) Let Ph denote a partition of n, Le. Ph is a finite collection of Ne open elements K i, i = 1,2,..., Ne. The partition is assumed to satisfy 0 = UKiE'Ph K i, and Ki n Kj = 0 whenever i is different from j. The size h K of an element K is measured in terms of its diameter. We associate with the partition the parameter h = maxke'ph hk. In order to approximate shape functions such as the solution u, one constructs a finite element space of hierarchical v h = {v E H1(O); V = 0 on 80 and VIK = v 0 Fj(l, v E PPK(K), VK E Ph} (7) where FK is the affine mapping from the master element K to the element K in the partition, and PpK(K) is the space of polynomial functions defined on K of degree at most PK. We employ here hierarchical shape functions based on the integrated Legendre polynomials [3, 4]. These shape functions can be classified into nodal, edge,face and interior bubble functions in three dimensions. We associate with the partition Ph the parameter P = minke'phpk. Then, the finite element space yh,p, or, for the sake of simplicity in the notation, yh, is defined as: Using the classical Galerkin method, a finite element approximation Uh E yh of U is given as the solution of B(Uh, v) = F(v), (9) The numerical error in the approximation shall be defined as e = U - Uh, e E Y. (5) (8) 3 Goal-oriented error estimates Rather than estimating the numerical error in the energy norm, we aim here at evaluating the error with respect to a quantity of interest. The quantity of interest represents the goal of the computational analysis, and is supposed to be defined mathematically as a bounded linear functional on the space of trial displacements Y. Many quantities of interest actually fall into that category. Let L denote such a functional. Therefore, our goal is to predict the quantity L(u). The numerical error in which we are now interested in is given by the quantity e L E lr, which, by linearity of L, reduces to: (10) 3
4 J.T. Oden, S. Prudhomme, T. Westermann, J. Bass We briefly explain how to estimate L(e). The reader is referred to [2,5] for more details. The starting point is to find a function w which relates the quantity L(e) to the residual 'R obtained from the momentum equation. This function, the so-called influence function, is shown to be the solution of the dual problem: B(v, w) = L(v), VVEV, (11) and is approximated in Vh by the dinite element solution wh E V h of the discrete dual problem: B(v, Wh) = L(v), (12) As for the error in the solution, we introduce the error in the approximate influence function as g = W - Who Then follows the remarkable result which establishes the relation between the error quantity L(e) and the energy norms of linear combinations of e and g. In Equation (13), s defines a scaling factor in R Several methods to obtain bounds on errors in the energy norm are now available in the literature. In the following numerical experiments, we shall use a simple method described in [6]. In order to derive an estimate of and bounds on L( e), we suppose that the following global error estimates are available: 1J~w ~ lise + s-lglle ~ 1J~p' 1J1~w ~ lise - s-lglle ~ 1J;;µp' Then, it is straightforward to show that the quantities 1Jtw and 1J!;µp, defined as: (14) L _ 1( +)2 1( - )2 1J/ow - 4 1J/ow - 4 1Jupp, L _ 1( +)2 1( - )2 1JIIPP - 4 1JIIPP - 4 1J/ow. (15) produce respectively a lower and an upper bounds on L( e), i.e. 1Jtw ~ L( e) ~ 1J~P' (16) An estimate of L(e) can be derived by taking the average of 1Jtw and 1J/~P' that is: L 1( L L 1Jest = "2 1J/ow + 1Jupp)' (17) Finally, the value of s is chosen so that the quantities Jlselleand Jls-lglle have same amplitudes, i.e. IIseile = Iis-lglle, which implies that: 8= (18) Such a choice of s is justified because it minimizes the quantities Jlse+ g/sll~ and lise - g/sll~. 4
5 ECCM-2001, Figure 1: Mesh topology and boundary conditions for cantelever beam. 4 Numerical experiments The main purpose of the present paper is to provide an account on the performance of goaloriented error estimation for three-dimensional elasticity problems. Our estimates and error bounds are verified by comparing the results with a refined and/or enriched mesh which is believed to provide a nearly "exact solution". We consider here as quantities of interest pointwise displacements and stresses. By poinwise quantities, we actually mean to take local averages using mollifying kernels defined in the form: ifr < c otherwise (19) The constant C corresponds to a normalizing factor and is chosen such that In ke;(x) dx = 1. It follows that the quantity of interest corresponding to the displacement at a point Xo E 0 will be represented by the linear functional: L(u) = Inke(lx - xol)u(x)dx ~ u(xo). (20) Likewise, for the stress a xx at a given point Xo, we shall consider the quantity of interest: L(u) = Inke;(lx - xol)axxdx ~ axx(xo). (21) The parameter c can be prescribed by the user and varies according to the problem and quantity of interest. In most cases, c is chosen of the order of the mesh size h. 4.1 The cantilever beam In the first example, we consider a cantilever beam with dimension L = (10,1,2) subjected to a uniform load. The beam is fixed at x = O. The mesh topology and boundary conditions are shown in Fig. 1. The material properties are set to E = 200 OPa and v = 0.3 where E is the Young's modulus and v the Poisson's ratio. 5
6 J.T. Oden, S. Prudhomme, T. Westermann, J. Bass Mesh Estimated error "Exact error" Elements p Estimate Lower Upper L(u - Uh) 4x1x x1x x 2 x x 4 x Table 1: Error estimates and bounds for the quantity of interest (displacement Uz at Xo (10.0,0.5, 1.0) in the cantilever beam). InIIueru:e fwlction for W iu x direction ~~~ InOueuco function for W in \' direction Ih1)(rSolve l.ot4s......,!"......,..,... Figure 2: Influence function corresponding to the displacement Uz at Xo = (10.0,0.5,1.0) (left) in x-direction (right) in y-direction Displacement at a point We estimate here the error in the displacement Uz (in direction z) at the free end-point Xo = (10.0,0.5,1.0). We select a sequence of meshes which are uniformly refined by dividing each element into two new elements in one or several directions. Estimates and bounds on the error in the displacement are shown in Table 1. The "exact errors" are computed here by substracting the discrete displacement from the displacement obtained on the refined mesh 32 x 4 x 8 with spectral order 2. The displacement at (10.0,0.5,1.0) on the mesh of reference reads We observe in Table 1 that the bounds and estimates, although very accurate, all slightly underestimate the exact errors. We attribute this behavior to the fact that the global estimates of the solution and influence function in the energy norm use edge bubble functions only and hence are rather coarse. Finally, we show in Fig. 2 the x- and y-component of the corresponding influence function Stress at a point In the next example, we estimate the error in the stress component a xx at Xo = (5.0,0.5,2.0). This point is located exactly at the center of the upper face of the beam, as it can be observed from the shape of the influence function shown in Fig. 3. As before we collect in Table 2 the estimates of the errors in the quantity of interest. This time, the estimates are rather accurate, but 6
7 ECCM-2001, Mesh Estimated error "Exact error" h p Estimate Lower Upper L(u - Uh) 4x1x x1x2 1 -] x 2 x x 4 x Table 2: Error estimates and bounds for the quantity of interest (Stress O'xx at Xo = (5.0,0.5,2.0) in the cantilever beam). Hu)erSoh"C' HvperSolv~ 02 Figure 3: Influence function corresponding to the stress O'xx at Xo = (5.0,0.5,2.0) (left) in x-direction (right) in z-direction. the bounds have the tendency to diverge from the exact error. This behavior was also observed in [2] for two-dimensional problems. We finally show in Fig. 4 the distribution of the elementwise contributions to the error in the displacement and stress component to be used as refinement indicators. We clearly observe the error in the displacement at the free-end is sensitive to the inaccuracy of the solution at the fixed-end, as this can be expected. On the other hand, we notice that the error in the stress has a more local origin. 4.2 The "hoop" problem The "hoop" problem is concerned with the simulation of a portion of a tire under a specific loading, as shown in Fig. 5. The blue arrows show the direction in which the hoop is held fixed while the orange arrows indicate a uniformly distributed load. The elements of the mesh are of order 2 along the surface of the tire and linear in the transversal direction and the provided mesh consists of 200 elements. Here, we analyze the stress component O'xx at Xo = (16,61,37) and show the errors for a sequence of meshes. In this case, as shown in Table 3, the bounds for the quantity of interest 7
8 J.T. Oden, S. Prudhomme, T. Westermann, J. Bass Contributions to the error ingol for W ~'i'ig,i~~e h.o.,... 1.&. 0" eo.dbuljom 10lhe enoc!nqoi for SIOl ,0.+00 rf' t.o.o.t I~ & Figure 4: Distribution of the contributions to the error in the quantity of interest (left) for the displacement at Xo = (10.0,0.5, 1.0) (right) for the stress at Xo = (5.0,0.5,2.0). Mesh and Loads lit_solve Figure 5: Mesh topology and boundary conditions for the "hoop" problem. are more accurate as compared to the beam problem. We also show in Fig. 6 the corresponding influence function and the distribution of the contributions to the error in the stress component. 4.3 The control arm In the final example, we consider a control arm whose geometry and loading are shown in Fig. 7. Once again, we are interested here in the stress component a xx whose distribution we show in Fig. 7 as well. The goal is to estimate the error in the computed stress at location Xo = (150,25,0) on the given mesh. We found the following estimate and bounds: rl:s~'" = CT",,,, 2 T/lolV = T/I~;; = to be compared with the "exact" error in the stress e(a xx ) = obtained by using a refined mesh with higher spectral order. We also show a plot of the influence function and the 8
9 ECCM-2001, Mesh Estimated error "Exact error" Elements p Estimate Lower Upper Initial mesh (1,2) Refine in (zeta) (1,2) Refine in (xi, zeta) (1,2) Table 3: Error estimates and bounds for the quantity of interest (Stress O'xx at Xo = (16,61,37) for the hoop problem). IIVJ'l(rSolve ~lothcem:cinq..~r[(lhellfui ~~I,. ",. I.' I ID )p I " "D Figure 6: Distribution of the contributions to the error in the quantity of interest (left) for the displacement at Xo = (57,23,37) (right) for the stress at Xo = (16,61,37). distribution of the error in the stress component in Fig. 8. On the latter, we observe that the contributions are not just local and that remote locations do contribute as well. An adapted mesh based on these elementwise contributions would then control the error in the quantity of interest by taking care of all contributing sources, near and remote. 5 Conclusions In this paper, we have applied the goal-oriented error estimation technique to the case of threedimensional linear elasticity problems. We have concentrated on estimating pointwise errors in displacement and stress components. The results obtained on the various examples shown here are very promising. In all cases, the error was predicted with very good accuracy. We would like to point out that relatively simple estimates of the error in the energy norm were used in those numerical experiments, which, to our own surprise, provided us with very reasonable estimates and bounds on the error. Acknowledgements This research was supported by a contract with the General Motors Research and Development Center. This support and the encouragement of Dr. Mark Botkin of GM are greatfully acknowledged. 9
10 J.T. Oden, S. Prudhomme, T. Westermann, J. Bass Control ami JbP<'rSolve Stress Sxx M~W,'~e ,;;;- "'",0> 100,... '''''' '00,(. - '00..., Figure 7: The control arm problem: (left) Geometry and loading (right) Distribution of the stress component axx. Influence function " 11 -~. u (.". Figure 8: The control arm problem: (left) Influence function in x-direction corresponding to the pointwise stress axx atxo = (150,25,0) (right) Distribution of the error in the stress component axx at Xo = (150,25,0). References [1] S. Prudhomme and J.T. aden. On goal-oriented error estimation for elliptic problems: Application to the control of pointwise errors. Compo Meth. in Appl. Mech. and Eng., 176, , (1999). [2] J.T. aden and S. Prudhomme. Goal-oriented error estimation and adaptivity for the finite element method. Computers Math. Applic., 41(5-6), , (2001). [3] B. Szab6 and J. Babuska. Finite Element Analysis. John Wiley & Sons, (1991). [4] J.T. aden. Optimal h-p finite element methods. Compo Meth. in Appl. Mech. and Eng., 112, , (1994). 10
11 ECCM-2001, [5] M. Ainsworth and J.T. Oden. A Posteriori Error Estimation in Finite Element Analysis. John Wiley & Sons, (2000). [6] S. Prudhomme and J.T. Oden. Simple techniques to improve the reliability of a posteriori error estimates for finite element approximations. ECCN 2001, (2001). 11
SIMPLE TECHNIQUES TO IMPROVE THE RELIABILITY OF APPROXIMATIONS
ECCM-2001 European Conference on Computational Mechanics 511 June 26-29, 2001 Cracow, Poland SIMPLE TECHNIQUES TO IMPROVE THE RELIABILITY OF A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS
More informationFirst Edition: Theory of Global Estimation of Modeling Error J. Tinsley Oden, 1 Serge Prudhomme, 1 Daniel C. Hammerand, 2 and Mieczyslaw S. Kuczma 3 1
First Edition: Theory of Global Estimation of Modeling Error J. Tinsley Oden, 1 Serge Prudhomme, 1 Daniel C. Hammerand, 2 and Mieczyslaw S. Kuczma 3 1 Texas Institute for Computational and Applied Mathematics
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More 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 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 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 informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationJ. Tinsley Oden, Kumar S. Vemaganti approaches. Preliminary results on extensions of these adaptive approaches to a class of models depicting material
Estimation of Local Modeling Error and Goal-Oriented Adaptive Modeling of Heterogeneous Materials; Part I : Error Estimates and Adaptive Algorithms J. Tinsley Oden and Kumar S. Vemaganti Texas Institute
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationA Posteriori Existence in Adaptive Computations
Report no. 06/11 A Posteriori Existence in Adaptive Computations Christoph Ortner This short note demonstrates that it is not necessary to assume the existence of exact solutions in an a posteriori error
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More 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 informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
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 informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
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 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 informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationRELIABILITY, UNCERTAINTY, ESTIMATES VALIDATION AND VERIFICATION
RELIABILITY, UNCERTAINTY, ESTIMATES VALIDATION AND VERIFICATION Transcription of I. Babuska s presentation at the workshop on the Elements of Predictability, J. Hopkins Univ, Baltimore, MD, Nov. 13-14,
More informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationModeling error and adaptivity in nonlinear continuum mechanics
Comput. Methods Appl. Mech. Engrg. 190 2001) 6663±6684 www.elsevier.com/locate/cma Modeling error and adaptivity in nonlinear continuum mechanics J. Tinsley Oden a, Serge Prudhomme a, *, Daniel C. Hammerand
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationSolving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells
Copyright 2010 Tech Science Press CMES, vol.57, no.2, pp.109-135, 2010 Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells E. J. Sellountos 1,
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationA domain decomposition algorithm for contact problems with Coulomb s friction
A domain decomposition algorithm for contact problems with Coulomb s friction J. Haslinger 1,R.Kučera 2, and T. Sassi 1 1 Introduction Contact problems of elasticity are used in many fields of science
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationError Control for Molecular Statics Problems
Error Control for Molecular Statics Problems Serge Prudhomme 1, Paul T. Bauman 2, and J. Tinsley Oden 3 Institute for Computational Engineering and Sciences The University of Texas at Austin Austin, Texas
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
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 informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 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 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 informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
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 informationAn example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction
An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationError Control for Molecular Statics Problems
International Journal for Multiscale Computational Engineering, 4(5&6)647 662(26) Error Control for Molecular Statics Problems Serge Prudhomme, Paul T. Bauman & J. Tinsley Oden Institute for Computational
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
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 informationModule 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method
Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without
More informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationA Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations
A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations M. Chiumenti y, Q. Valverde Λ, C. Agelet de Saracibar y and M. Cervera y ydept. de Resistència
More informationCOMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES
COMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES P. M. Mohite and C. S. Upadhyay** Department of Aerospace Engineering, IIT Kanpur 0806, INDIA, e-mail: mohite@iitk.ac.in Assistant Professor,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationDYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,
More informationA spacetime DPG method for acoustic waves
A spacetime DPG method for acoustic waves Rainer Schneckenleitner JKU Linz January 30, 2018 ( JKU Linz ) Spacetime DPG method January 30, 2018 1 / 41 Overview 1 The transient wave problem 2 The broken
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationFRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
Proc. 2 nd ESIS TC4 Conference on Polymers and Composites, in press, 1999 Author prepared reprint FRACTURE MECHANICS OF COMPOSITES WITH RESIDUAL STRESSES, TRACTION-LOADED CRACKS, AND IMPERFECT INTERFACES
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
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 informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationLimit Analysis with the. Department of Mathematics and Computer Science. Odense University. Campusvej 55, DK{5230 Odense M, Denmark.
Limit Analysis with the Dual Ane Scaling Algorithm Knud D. Andersen Edmund Christiansen Department of Mathematics and Computer Science Odense University Campusvej 55, DK{5230 Odense M, Denmark e-mail:
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 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 informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
More informationOn Mixed Methods for Signorini Problems
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 30, 2003, Pages 45 52 ISS: 1223-6934 On Mixed Methods for Signorini Problems Faker Ben Belgacem, Yves Renard, and Leila Slimane Abstract. We
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More information