A Feed-Back Approach to Error Control in Finite Element. Methods: Basic Analysis and Examples. Roland Becker and Rolf Rannacher 1
|
|
- Morgan Lindsey
- 5 years ago
- Views:
Transcription
1 A Feed-Back Approac to Error Control in Finite Element Metods: Basic Analysis and Examples Roland Becker and Rolf Rannacer 1 Summary. A rened approac to residual-based error control in nite element (FE) discretizations is presented. Te conventional strategies for adaptive mes renement in FE metods are mostly based on a posteriori error estimates in te global energy or L 2 -norm involving local residuals of te computed solution. Te mes renement process ten aims at equilibrating tese local error indicators. Suc estimates reect te approximation properties of te nite element space by local interpolation constants wile te stability properties of te continuous model enter troug a global stability constant, wic may be known explicitly in simple cases. Meses generated on te basis of suc global error estimates may not be appropriate in cases of strongly varying coecients and for te computation of local quantities as, for example, point values or contour integrals. More detailed information about te mecanism of error propagation can be obtained by employing duality arguments specially adapted to te quantity of interest. Tis results in a posteriori error estimates in wic te local information derived from te dual solution is used in te form of weigts multiplied by local residuals. On te basis of suc estimates, a feed-back process in wic te weigts are numerically computed wit increasing accuracy leads to almost optimal meses for various kinds of error functionals. Tis approac is developed ere for a simple model problem, namely te Poisson equation in two dimensions, in order to empasize its basic features. However, te underlying concept is rater universal and as, on a euristic basis, already been successfully applied to muc more complex problems in structural and uid mecanics as well as in astropysics. 1 Introduction Te strategies for mes renement conventionally used in FE metods are mostly based on a posteriori error estimates in global norms, e.g., te energy or te L 2 -norm, involving local residuals of te computed solution. Te mes renement process ten aims at equilibrating tese local error indicators. However, meses generated on te basis of suc global error estimates may not be appropriate in cases of strongly varying coecients and for controlling te accuracy in approximating local quantities as, e.g., point values or contour integrals. For tis one needs more detailed information on te mecanism of error propagation wic can be obtained by employing suitable duality arguments known from te a priori error analysis as te so-called Aubin-Nitsce trick (see, e.g., [11]). Te corresponding dual solutions ten yield te appropriate weigt factors to be used in te a posteriori error estimates. Consider, for example, te Poisson equation? u = f in ; u = 0 ; (1.2) on a bounded domain R 2 wic, for te moment, is assumed to be convex polygonal. Te restriction to two dimensions is only for simplicity as te extension to arbitrary dimensions is 1 Tis work as been supported by te German Researc Foundation (DFG), SFB 359 Reactive Flow, Diusion and Transport, Institute of Applied Matematics, University of Heidelberg, INF 293, D Heidelberg, Germany 1
2 straigtforward. By (; ) we denote te L 2 inner product and by k k te corresponding norm on. Let (1.2) be discretized by a standard nite element Galerkin metod using piecewise polynomial, say, linear or bilinear sape functions on meses T = fg satisfying te usual regularity conditions (see, e.g., [11]). To facilitate local mes renement, anging nodes may be allowed. For eac T, let max = max 2T, were = diam(). Te corresponding nite element subspaces are V V := H 1 0 (). Teir approximation properties are caracterized in terms of local approximation estimates (see ([11]), o max nkv? I vk ; 1=2 kv? I C i; 1+r kr1+r vk ~ ; r 2 f0; 1g ; (1.3) for v 2 V \ H 1+r (), were I v 2 V is some locally dened approximation to v, and k k B denotes te L 2 -norm over a set B. For r = 0, ~ is te union of all neigboring elements of, wile for r = 1, one simply as ~ =. Starting from te variational formulation of (1.2), te nite element metod seeks to determine approximations u 2 V, suc tat (ru ; r' ) = (f; ' ) 8 ' 2 V : (1.4) Subtracting (1.4) from te variational formulation of (1.2) results in te following ortogonality relation for te error e = u? u, (re; r' ) = 0 8 ' 2 V : (1.5) Tis discretization allows for optimal-order a priori estimates in te energy and L 2 -norm, kek + max krek C 2 maxkr 2 uk; (1.6) provided tat te solution u is suciently regular. Corresponding estimates also old wit respect to te maximum norm (for references see, e.g., [11]). For error control in te energy or L 2 -norm, one may proceed as follows. Using Galerkin ortogonality (1.5) and integration by parts on eac element yields (re; rz) = 2T n (f + u ; z? z )? 1 2 (n [ru ]; z? z o ; (1.7) for any z 2 V, were [ru ] denotes te jump of ru across te element boundary, and z 2 V is a suitable approximation of z. On edges along te boundary, we set [ru ] = ru. Ten, using Holder's inequality on eac element, we obtain j(re; rz)j 2T! ; (1.8) wit local residuals and weigts! dened by := kf + u k =2 kn [ru ; (1.9) n o! := max?1 kz? z k ;?1=2 kz? z : (1.10) By virtue of (1.3), te approximation z may be cosen suc tat for r 2 f0; 1g, tere olds! C i; r kr1+r zk ~ : (1.11) 2
3 Consequently, it follows tat j(re; rz)j C i 2T 2r 2 1=2kr 1+r zk ; (1.12) for z 2 H 1+r (). Te interpolation constant C i is usually of size 0:1 C i 1, depending on te sape of te elements. Taking now te supremum over z 2 V \ H 1+r (), one obtains te following a posteriori error bound in te energy norm, for r = 0, or in te L 2 -norm, for r = 1, kr 1?r ek C s C i 2T 2r 2 1=2 =: (u ) : (1.13) Te stability constant C s measures te stability properties of te dual problem in terms of te global a priori estimate (r'; rz) = (r 1?r e; r 1?r ') 8 ' 2 V; (1.14) kr 1+r zk C s kr 1?r ek ; (1.15) wic, of course, is trivial for r = 0. Clearly, by analogous arguments a posteriori error estimates may also be derived wit respect to L q -norms including te limit case q = 1. For earlier work on a posteriori error estimates for FE metods, we refer to te pioneering papers of Babuska and Reinboldt [1], Bank and Weiser [3], and Babuska and Gui [2]. Te underlying approaces ave been surveyed in Verfurt [28] and [30]. Te concept of a posteriori error estimation by duality arguments as been developed by Eriksson, Jonson and teir coworkers (see [13], [21], [15], and te literature cited terein). Heuristically based error indicators for elliptic problems ave been devised, e.g., by ienkiewicz and u [34]. Based on an a posteriori error estimate of te type (1.13) te mes renement process may be organized according to te so-called error per cell strategy. For some prescribed error tolerance T OL, te goal is to reac a most economical mes T on wic (u ) T OL. Accordingly, te mes renement process aims at equilibrating te local error indicators := r by rening (or coarsening) te elements 2 T according to te criterion T OL p NCi C s ; N = #f 2 T g : (1.16) Te renement strategies based on te conventional a posteriori error estimate (1.13) rely on te assumption tat te local error indicators := r properly describe te dependence of te global error on te local mes size. However, tis may not be true in certain situations since te a posteriori error estimate (1.13) contains information about te mecanism of error propagation only troug te global stability constant C s. To overcome tis deciency, it as been proposed in [7] to use te quantities! in te estimate (1.8) as weigt-factors multiplied by te local residuals and to compute tem numerically. Tese weigts contain all information about te local approximation properties of te spaces V, as well as te local stability properties of te underlying continuous problem. Tis correspondence can be used for te mes renement algoritm troug a feed-back process. In te course of te renement te dual solution is calculated on te current mes yielding approximate weigts. On te basis of te resulting 3
4 a posteriori error estimate te mes is rened according to one of te criteria described below. Tis process is repeated yielding more and more accurate weigts, i.e., a posteriori error bounds, until te prescribed stopping criterion is fullled. Tis approac allows one to construct almost optimal meses for various kinds of error quantities, were \optimal" can mean \most economical for acieving a prescribed accuracy T OL" or \most accurate for a given maximal number N max of mes points". 2 Residual based weigted a posteriori error estimates First, we illustrate our concept at a simple model situation. Let = (?1; 1) 2 and te pointvalue u(0) be te quantity to be computed wit best accuracy. It is not likely tat te meses obtained troug energy or global L 1 -error control are optimal for tis purpose. To do better, we ave rst to derive an a posteriori error estimate for te point-error. Tis can be done by a duality argument. Let z be te Green function corresponding to te origin, i.e., te solution of te \dual" problem (r'; rz) = '(0) 8 ' 2 V \ C( ) ; (2.2) were, in tis special situation, z(x) log(jxj). Ten, we ave e(0) = (re; r(z? z )); (2.3) wit te nodal interpolant z 2 V, appropriately regularized at x = 0. Consequently, in view of te error identity (1.7), tere olds je(0)j 2T! ; (2.4) wit residuals and weigts! as dened above. Using te particular form of te dual solution, and te usual local interpolation properties of nite elements, tis error estimator takes te form je(0)j C i 2 2T r 2 were r := maxfdist(; 0); g and C i 1. =: point (u ) ; (2.5) Suppose tat for some prescribed tolerance T OL te renement process as led to a nal mes T opt on wic te local error indicators =! are equilibrated, Tis implies te relation T OL N opt ; N opt = #f T opt g : (2.6) T OL r 2 N opt 1=4 ; (2.7) if we assume tat 2. Te latter assumption means tat on te current mes, tere olds 1=2 kn [ru 2 max jn?1 [ru ]j c(u) 2 ; (2.8)
5 a ypotesis wic is supported by computational experiments (see te te examples, below). As consequence, te number N opt of elements of te nal mes T opt is given by N opt 2?2 = N 1=2 2 opt 2T opt 2T opt T OL 1=2 r : (2.9) Since r?1 is integrable over, it follows tat N opt TOL?1 and te \optimal" mes-size distribution is (TOLr ) 1=2. In contrast, te global energy-error estimator (1.13), for r = 0, would generate meses of complexity N opt T OL?2, in order to guarantee an asymptotic beavior krek TOL. Te corresponding global L 2 -error estimator 1=2 =: 2 (u ) (2.10) kek C s C i 2T 2 2 seems to be more appropriate as, at least, it as te rigt asymptotic beavior as! 0. Alternatively, one may use te L 1 -error estimator wic is obtained by an L 1 =L 1 -duality argument (see te proof of Lemma 4.2 in Section 4). It takes te form max jej C i C s j log 0 j 1=2 max 2T 1; =: 1 (u ) ; (2.11) were 0 is te diameter of te element containing te point x = 0, and 1; := 2 max jf + u j + max jn [ru ]j Tis estimator also refers to a global error norm and does not particularly empasize te point x = 0 at wic te solution u is to be computed. Terefore, it will not be able to produce a most economical mes for tis particular purpose. However, asymptotically, it yields meses of te same complexity in terms of TOL, N opt TOL?1, as te point-value estimator point (). Tis situation canges wen point-values of te gradient ru(0) (\stress values") are to be computed. Ten, te global error bounds require N opt TOL?2 cells to acieve a corresponding error of size TOL, wile te weigted error estimator acieves tis again wit only N opt TOL?1 cells (see Example 1, below). Of course, in more general situations, explicit bounds for te dual solution z are not known and ave to be computed numerically. Tis essential question will be discussed in more detail in Section 4. Te simplest approac is based on te local interpolation estimate (1.11) for r = 1. Estimates for te weigts! may be obtained from te approximate dual solution ~z 2 V, computed on te current mes, by taking appropriate 2nd-order dierence quotients, were x denotes te mid-point of element.! ~! := 2 jr2 ~z (x )j ; (2.12) Heuristically, te meaning of our a posteriori error estimate is as follows. We expect tat te weigts! converge to certain limits, as! 0, in te sense?1=2!! jd 2 z(x )j, were x denotes a point contained in an innite sequence of nested elements. Furter, te quantity?1 n [ru ] can be viewed as a certain second-order dierence quotient of u. Hence, it may be expected tat, for a sequence of properly rened meses, te residual terms in (2.4) also converge to certain limits as! 0, jj?1=2 kf + u k! j(f + u)(x )j = 0; j@j?1=2 k?1 n [ru jd 2 u(x )j: 5
6 In tis sense, te a posteriori error estimate (2.4) asymptotically, for! 0, takes te form je(0)j 2T 4 jd 2 z(x )jjd 2 u(x )j: (2.13) Te local residuals consist of two parts, te domain residual kf + u k and te normal-jump terms 1=2 kn [ru For linear elements, it can be argued tat, in te case of smoot f, te contribution of te normal jumps asymptotically dominates tat of te domain residual, and te latter may terefore be neglected. Tis aspect will be discussed in more detail in Section 4. Hence, wit te above notation, we are led to te practical error estimator were je(0)j C i ~! ~ =: approx 2T point (u ) ; (2.14) ~ := 1=2 kn [ru ; ~! := 3=2 kn [r~z ; (2.15) and te interpolation constant may be set to C i 1. Numerical test. Te weigted error estimator (2.5) and its approximate version (2.14) ave been tested, wit C i = 0:1, for te solution u = 10 sin(2x 1 +x 2 +2) using bilinear elements. For successively reduced tolerances TOL, te resulting balanced meses T opt and tose obtained by te global L 2 -error estimator (2.10) and te L 1 -estimator (2.11) are sown in Figure 1. Te corresponding errors are listed in Table 1. Here, L is te number of renement levels of te nal mes T opt and I ef f := je(0)j=(u ). For te problem considered, te predicted dependence N opt TOL?1 for te weigted a posteriori estimate (2.4) is conrmed. Te sligt tendency of our error estimator to underestimate te true error could be compensated by using a larger interpolation constant C i. 6
7 2 1 point Figure 1: Rened meses (left to rigt) for computing u(0) for u = 10 sin(2x 1 + x 2 + 2), using te L 2 - and L 1 -error estimators 2 and 1, respectively, compared to te weigted error estimator point Te above argument directly generalizes to te case of an arbitrary linear error functional J() dened on te space V, or on a suitable subspace containing te nite element space V and te exact solution u. Oter relevant cases are, for instance, te moments of u over (torsion moment), te value of its derivative ru(x 0 ) at some point x 0 2 (stress values), or te weigted integral of te normal n u over te (total surface tension), J(') = ' dx; J(') = r'(x 0 ); n ' ds : (2.16) 7
8 2 TOL N L je(0)j I ef f e e ? e e ? e e ? e e ? e e ? e e ? e e TOL N L je(0)j I ef f e e ? e e ? e e ? e e ? e e ? e e ? e e point TOL N L je(0)j I e e e ? e e ? e e ? e e ? e e ? e e ? e e approx point TOL N L je(0)j I e e e ? e e ? e e ? e e ? e e ? e e ? e e Table 1: Results obtained for computing u(0) for u = 10 sin(2x 1 + x 2 + 2), using te L 2 - and 1 L - error estimators 2 and 1, respectively, compared to te weigted estimator point and its approximate analogue approx point, all wit C i = 0:1 Ten, following te above approac, we ave te a posteriori error estimate jj(e)j 2T! ; (2.17) wit weigts! and residuals as dened above. We like to empasize tat, in tis particular problem, it is crucial to utilize te approximation property of te nite element space to its full extent also in te evaluation of te weigts!, in order to obtain economical meses. Tis applies even toug tere may be no a priori bound for te resulting norms of te continuous dual solution z (see [24], [10] and [19]). We will come back to tis point in Example 4 at te end of tis paper. In many interesting cases, te functional J() is \singular", i.e., not properly dened on H 1 0 () but only on some subspace containing te solution u and on te nite element space V. Typical examples are te evaluation of point-values or of contour integrals of derivatives. In suc a case, it is advisable to regularize te functional according to te prescribed tolerance, in order to avoid over-renement near te singularities. Te type of regularization depends on te particular situation and sould be suc tat jj(u)? J T OL (u)j TOL : (2.18) For example, te point-value estimator (2.14) may be used in te regularized form je(0)j 2T 2 (r + TOL) 2 : (2.19) 8
9 Furter examples of tis type will be discussed in Section 5. Since some of te examples presented below deal wit domains wit a curved boundary and also iger order nite elements, we ave to extend our a posteriori error analysis to tis situation. Now, te domain is assumed eiter to ave a smoot boundary or to be convex polygonal, in order to ensure tat te usual L 2 =H 2 -sift teorem olds true. Accordingly, along te te elements may be curved (isoparametric elements). For te mes size, we use te same notation, etc., as above. In addition, we now introduce te mes domain := [f 2 T g and assume tat all nodal points lie Tis means tat, in general, 6 for a non-convex domain. In suc a case, one as to cope wit a boundary strip f n g [ f n g of local widt ; m+1. Te nite element spaces are dened by V := fv 2 H 1 0 ( ); v j 2 P ()g ; were P () is a suitable vector space of polynomials over containing P m (). Tey approximate te solution space V = H0(), 1 suc tat local approximation properties (1.3) old for v 2 V \ H r+1 ( [ ) wit 1 r m. We introduce te subset of elements adjacent to te boundary, := f 2 T ; meas(@ ) 6= 0g and accordingly set T 0 := T n To eac element 2 tere naturally belongs a section of te strip f n g [ f n g denoted by S. By a variant of te Poincare inequality, tere olds In tis case te nite element sceme reads kvk S C i; m+1 krvk : (2.20) (ru ; r' ) = ( f; ' ) 8 ' 2 V ; (2.21) were f (if necessary) is a suitable extension of f to [. For te following, we extend also te nite element solution u as well as te dual solution z by zero to [ (in accordance wit omogeneous boundary conditions). Using tis convention, tere olds (re; rz) = (f; z)? (ru ; rz) \? (ru ; rz) n = (f; z) n + ( f; z)? (ru ; rz) = (f; z) n + ( f; z? z )? (ru ; r(z? z )) ; from wic by integration by parts, we conclude te following analogue of te identity (1.7), o (re; rz) = n(f + u ; z? z )? 12 (n [ru ]; z? z (2.22) 2T 0 + o n( f + u ; z? z )? 12 (n [ru ]; z? z + (f; z) S : Ten, te weigted a posteriori estimate (2.17) olds wit te local residuals and weigts as dened above, wit te only modication tat for 2 := kf + u k =2 kn [ru + kfk S ; (2.23) n o! := max?1 kz? z k ;?1=2 kz? z ; kzk S : (2.24) In view of (2.20), we ave! C i; m+1 kzk H m+1 () for 2 For a detailed analysis of te eect of boundary approximation using similar arguments, we refer to [12]. 9
10 3 A feed-back algoritm for adaptive mes renement We briey discuss ow a mes renement process may be organized on te basis of an a posteriori error estimate of te type (2.17) or (2.14). Suppose tat some error tolerance T OL and maximum number N max of mes points are given. Te goal is to nd a most economical mes T on wic jj(e)j (u ) = 2T T OL ; (3.2) wit te local error indicators :=!. Usually, one starts from an initial coarse mes wic is ten successively rened. Tere are essentially tree alternative strategies: 1. Error per cell strategy. Te mes generation aims at equilibrating te local error indicators, by rening (or coarsening) te elements 2 T according to te criterion T OL N ; N = #f 2 T g: Since N depends on te result of te renement decision, tis strategy is implicit and usually needs iteration. It is common practice to work wit a varying value N on eac renement level wic is permanently updated according to te renement process. Te result is a mes on wic (u ) T OL, provided tat N max is not exceeded. 2. Fixed fraction strategy. In eac renement cycle, te elements are ordered according to te size of and eiter a xed portion (say 30%) of te elements wit largest or te portion of elements wic make up for a certain part of te estimator, (u ), is rened. Te appropriate coice of te parameter is crucial and depends very muc on te particular situation. For \regular" functionals, one may coose = 0:6?0:8, wile for \singular" functions a smaller coice = 0:1? 0:2 is advisable, in order to enance local renement. Tis process is repeated until te stopping criterion (u ) T OL is fullled, or N max is exceeded. 3. Tolerance reduction strategy. One works wit a varying tolerance T OL var. If on a mes T a discrete solution u old as been obtained wit corresponding error estimator (u old ), te tolerance is set to T OL var := (u old ) ; wit some reduction factor 2 (0; 1) (usually = 0:5). In te next step, one (or more) cycles of te error per cell strategy are applied wit tolerance T OL var yielding a rened mes T new and te new solution u new wit corresponding error bound (u new ). Ten, te tolerance is reduced again and a new renement cycle begins. Tis process is repeated until T OL var T OL, or N max is exceeded. Te error per cell strategy is very delicate as its performance strongly depends on te parameters in te renement decision. In certain cases it may appen tat an inappropriate coice leads to very slow mes renement and in turn to an inecient overall solution. On te oter and, te xed fraction strategy guarantees tat in eac renement cycle a suciently large number of elements is rened. However, in its pure form, it does not allow for mes coarsening and in certain cases may tend to over-rene te mes. Terefore, for practical computations, we recommend te tolerance reduction strategy wic seems to be te most robust and accurate strategy among te tree. Its use is particularly advisable if te weigt-function!(x) as 10
11 nonintegrable singularities wic would oterwise dominate te renement process (see te examples at te end). However, for simplicity, in all te test calculations presented below, te \optimal" meses ave been generated by te xed fraction strategy based on te simplied estimator (2.14) T opt ~(u ) := C i kn [ru kn [rz (1) 2T ; were z (1) is te bilinear approximation to z obtained on te current mes. Te interpolation constant is usually set C i = 0:1, in order to compensate for te over-estimation in tis bound. Accordingly, te local error indicators are dened as := C i kn [ru kn [rz (1) : (3.3) 4 Model analysis of te weigted a posteriori error estimation We will give some teoretical support for te proposed approac to adaptive error control by looking in more detail at te model problem introduced in te previous section. Te mes renement process is organized as described above using (isoparametric) bilinear sape functions on quadrilateral elements. Eac element is allowed to ave just one anging node on eac of its edges. Te crucial question is weter te a posteriori error bounds obtained by te proposed approac are not only asymptotically optimal (order-optimal) but also asymptotically sarp in te sense jj(e)j I e := lim T OL!0 (u ) = 1 ; possibly depending on te work te user is willing to invest in te evaluation of te weigts!, i.e., in te numerical solution of te dual problem. First, we run some computational experiments in order to see wat can be expected. Unfortunately, te result is rater negative, i.e., wit acceptable eort asymptotic sarpness does not seem to be acievable. Wit all te metods considered te eciency index I e never really tends to one, but in most relevant cases stays well below 2, wat may actually be considered as good enoug. Tere are two separate aspects to be considered: Te sarpness of te global error bound (u ), and te optimality of te local error indicators =! wic are used in te mes renement process. Starting from te error representation (1.7), J(e) = 2T n (f + u ; z? z )? 1 2 (n [ru ]; z? z o ; (4.2) tere are various steps in evaluating te rigt and side. First, one estimates jj(e)j 2T ; (4.3) wit te (positive) local error indicators := (f + u ; z? z )? 1 2 (n [ru ]; z? z ; (4.4) 11
12 wic are used in te mes renement process. By tis step, in general, te asymptotic sarpness of te global error estimate is lost. Tis is seen, for instance, in te case were te exact as well as te approximate solution are anti-symmetric wit respect to te x-axis meaning tat e(0) = 0, but (u ) 6= 0. Terefore, we will analyze various ways of directly evaluating te error identity (4.2). Here, z 2 V may be taken as te bilinear interpolation of te approximation for z on te current mes, or simply be set zero, as its global contributions (if exactly evaluated) cancel out to zero under summation. Te approximate values obtained for J(e) are denoted by (i) (u ). 1. Approximation by iger order metods: Te dual problem is solved by using biquadratic nite elements on te current mes yielding an approximation z (2) to z. Te error estimator obtained by using tis approximation in (4.3) and (4.4) is denoted by (1) (). Computing approximations to z by using a iger order nite element sceme does not appear very economical in estimating te error in te low-order sceme. Hence, one may try to work wit te given metod, i.e., piecewise bilinear approximation, and increase its accuracy troug defect correction wit biquadratic interpolation on a coarser mes. Notice tat simple Ricardson extrapolation using two consecutive meses yields improved approximations only at te nodes of te coarser mes. However, since tis procedure requires to solve at least twice on te current mes, it appears too costly for practical purposes and is not furter considered. 2. Approximation by iger order interpolation: A furter simplication is acieved by simple local biquadratic interpolation of te bilinear approximation z (1) on te current mes yielding an approximation I (2) z(1) to z. Te resulting global error estimator is denoted by (2) (). Tis requires some special care on elements wit anging nodes, in order to preserve te iger order accuracy of te interpolation process. 3. Approximation by iger order dierence quotients: Witout loss of generality, consider a rectangular element oriented in te cartesian directions wit side lengts 1 = ;1 and 2 = ;2. Its vertices are numbered in te counter-clockwise sense as a i = (x i 1; x i 2) (i = 0; :::; 3). On te error in te bilinear interpolation can be expanded like (z? I (1) z)(x) = E (x) 2 + R ( ; x) ; (4.5) wit an expansion coecient E (x) and a remainder term R ( ; x) = O( 3 ). Tere olds E (x) = n'(' o?2? 1)( z 11 (a 0 ) + z 11 (a 3 )) + (? 1)('z 11 (a 1 ) + ' z 11 (a 2 )) ?2 n (? 1)( z 22 (a 0 ) + 'z 22 (a 1 )) + (? 1)( z 22 (a 1 ) + z 22 (a 2 )) were ' = (x 1? x 0 1)=(x 1 1? x 0 1), = (x 1 1? x 1 )=(x 1 1? x 0 1), = (x 2? x 0 2)=(x 1 2? x 0 2), = (x 1 2? x 2 )=(x 1 2? x0 2 ), and jr ( ; x)j 3 max jr 3 zj. Te elementary but tedious computation is omitted. According to (4.5), te local error indicators are estimated in te form 2 (f + u ; E )? 1(n [ru 2 ]; E : (4.6) Te second derivatives of z in te expansion coecient E are ten approximated by suitable central dierence quotients of its bilinear approximation z (1) obtained on te current mes. Again at anging nodes special care is required. Te resulting global error estimator is denoted by (3) (). o ; 12
13 4. Approximation by local problems: Following te idea already used in te energy-error estimator of Bank and Weiser [3], on eac element te local problems (rv ; r' ) = (f + u ; ' )? 1 2 (n [ru ]; ' 8 ' 2 V ; (4.7) are solved, were V = Q 2 () te space of biquadratic polynomials on. Using te resulting solutions, te local residual terms are replaced by ~ := j(rv ; r(z? z )) j : (4.8) Te resulting global error estimator is denoted by (4) (). Te related error indicator based on solving patc-wise Diriclet problems is not considered ere, since it turns out to be too expensive in practice. From te various results obtained in our tests, we quote tose for te point-value computation discussed above, as well as tose for te error measured in te global L 2 -norm. Te observed eciency indices are listed in Tables 2 and 3. Tese results indicate tat even for te simplest model situations, an eciency index I e = 1 is acievable only at te expense of unacceptably ig cost, e.g., by approximating te dual solution using iger-order elements. L N J(e) = e(0) J(e)= (1) J(e)= (2) J(e)= (3) J(e)= (4) e e e e e Table 2: Eciency of various weigted error indicators obtained for te point-error J(e) = je(0)j L N J(e) = kek J(e)= (1) J(e)= (2) J(e)= (3) J(e)= (4) e e e e e Table 3: Eciency of various weigted error indicators obtained for te L 2 -error J(e) = kek To prove tat te approximation strategies described above actually work in constructing \optimal" meses seems very dicult, as for tis one would need precise information about te local (on te element level) approximation properties of te nite element Ritz projection on general unstructured meses. Unfortunately, results of tis type are not available in literature and may not even be true at all. Below, we will prove a result for te point-error estimator wic is not quite wat one would like to ave, but provides some partial justication. It states tat, at least, on an optimally rened mes on te basis of te \exact" estimator (u ), te perturbed estimator ~(u ) using a iger-order approximation to te dual solution z sows te 13
14 rigt beavior in terms of te tolerance TOL. To tis end, we ave to make some reasonable tecnical assumptions. First, suppose tat on te \optimal" meses T, tere olds c 2 n o max jf + u j + max jn r[u ]?1 j C(u) 2 2 T ; wic indicates te (global) stability of te nite element Ritz projection wit respect to te discrete W1 2 -norm on general locally rened meses. On quasi-uniform meses (satisfying te uniform sape and size condition) tis immediately follows by te known error estimates in te W1-norm 1 and te inverse properties of nite elements. For more general meses, tis as to be left as a conjecture. However, (4.9) may be cecked in te course of te mes renement process. Lemma 4.1. For te model problem (1.2) on te square = (?1; 1) 2, let te dual solution z, i.e., te Green function to te point x = 0, be approximated on te current mes by its Ritz projection z (2) into te subspace of piecewise quadratic functions. Furter, suppose tat for a sequence of tolerances, TOL! 0, optimally balanced meses T opt ave been generated on te basis of te exact estimator point (u ), suc tat point (u ) TOL and TOL N opt ; 2 T opt : (4.10) Ten, te approximate error estimator ~ point (u ) := n 2T is asymptotically exact in te sense (f + u ; z (2) )? 1 2 (n [ru ]; z (2) (4.11) j point (u )? ~ point (u )j = o(t OL) : (4.12) Proof. We only give a sketc of te proof wic is based on results from local nite element error analysis (see [11]). Using te abbreviation w := z? z (2), we ave by Galerkin ortogonality point (u )? ~ point (u ) = 2T n (f + u ; w? I w )? 1 2 (n [ru ]; w? I w o ; and, consequently, by te local approximation properties of te interpolation operator I, Tis is furter estimated by point (u )? ~ point (u ) point (u )? ~ point (u ) C i Te second factor on te rigt is bounded by 2T kr 2 w k : 1=2 C i C(u) 2 2 r?4 rkr 4 2 w k 2 2T 2T 2T r 4 kr2 w k 2 1=2 14 r 4 jr 2 w j 2 dx 1=2 ; 1=2 : (4.13)
15 were r(x) := maxfjxj; 0 g, 0 = min being te diameter of te element containing te point x = 0. By standard arguments from local nite element analysis, one can prove tat r 4 jr 2 w j 2 dx 1=2 cmax c max r 4 jr 3 zj 2 dx 1=2 r?2 dx 1=2 cmax jlog( 0 )j 1=2 : We omit te rater tecnical details and refer to [18] and [11] for similar arguments. Under te assumptions (4.10) and (4.9), for te optimal mes T, tere olds N T OL?1 and (T OLr ) 1=2. Terefore, te rst factor on te rigt of (4.13) can be bounded by 1=2 1=2 C(u)T OL C(u) jlog(0 )j 1=2 T OL : 2 2T 2 r?4 r?2 2T Combining te foregoing estimates in (4.13) yields point (u )? ~ point (u ) C i C(u) jlog( 0 )j max T OL : In view of max T OL 1=2, tis implies te assertion. In view of te computational results sown in Table 3, te statement of Lemma 4.1 does not seem to be true for te approximation of te global L 2 -error. Since asymptotically sarp a posteriori control does not seem to be feasible in general, we ave to live wit error estimates involving at least certain interpolation constants C i. Tese constants may usually be assumed to be of moderate size C i = 0:1?1, but certainly represent some factor of uncertainty. In particular situations, one may try to calibrate tese constants in te course of te renement process wic, owever, is a very delicate matter. Te full teoretical justication of te error estimator ~() would require an analysis of te eect of replacing te exact second derivatives of te dual solution z by suitable dierence quotients of some approximation z obtained numerically on te current mes. Yet, tis is strongly supported by numerical evidence, even in te present simple case te proof must be left as an open problem. For te residual-based energy error estimator (1.13) one as te two-sided bound krek energy (u ) fkrek + kf? f kg ; (4.14) were f is a suitable approximation of f. Tis sows asymptotic optimality. An analogous bound cannot old in general for local error quantities like te point-error. Wat we can ope to acieve is only an optimality estimate of te type presented in te following lemma. Lemma 4.2. For te nite element approximation of te model problem (1.2) on te square = (?1; 1) 2, te error estimator point () dened in (2.4) is asymptotically optimal in te sense je(0)j point (u ) c jlog( 0 )j max 2T max were 0 is te diameter of te element containing te point x = 0. jej + 2 max jr2 uj ; (4.15) 15
16 Proof. Starting from (2.4), we use te Holder inequality to obtain point (u ) 2 2 2T r 2 1=2 2T?2 r2!2 1=2; (4.16) were r := maxfdist(; 0); g. Te second factor can be estimated by (see [18] or [11]) 2T?2 r2! 2 c f 0 <jxj<1g jxj 2 jr 2 zj 2 dx + cjlog( 0 )j c jlog( 0 )j : (4.17) Furter, using te usual inverse inequality for nite elements togeter wit te local interpolation estimate (1.3), we conclude in a standard way tat kek =2 kn c?1 kek ~ + c kr 2 uk ~ ; were ~ denotes te union of all cells in T intersecting. Tis implies 2 2T r 2 Combining (4.17) and (4.18), te assertion follows. 2 c jlog( 0 )j max max 2T jej + 2 ~ 2 max ~ uj jr2 : (4.18) For a special situation (uniform tensor-product mes) it as been sown in [32] and [33] tat in te energy-error estimate (1.13) eiter te domain residual terms kf + u k or te jump terms kn [ru may be asymptotically neglected depending on weter te degree of te elements used is odd or even, respectively. Tis result as recently been extended in [31] to general adaptively rened meses for lowest-order elements (odd degree). We will derive a corresponding result for te present situation wic to some extent justies te use of te simplied weigted error estimator ~() dened in (2.14). Lemma 4.3. Suppose tat f 2 H 1 () in te model problem (1.2). Ten, for piecewise bilinear elements, tere olds jj(e)j C i 2T were ~ is te union of all cells in T intersecting. 3 krfk + ~ 3=2 kn [ru kr 2 zk ; (4.19) ~ Proof. For simplicity, we give te proof only for rectangular meses. Te generalization to more general meses requires some more tecnical work. In tis case u j 0 on eac element. Since te mes T is generated from a coarse mes by local cell-wise bisection (i.e., subdivision into four sub-cells), it can be decomposed into macro-cells consisting of four sub-cells eac. Tese macro-cells form a coarser mes named ~ T. Te case of cells wit anging nodes requires obvious modications and will not be discussed in detail. On ~T we dene an interpolation operator ~I : C( )! V, by requiring ~ I v(a) = v(a), for vertices a of ~T, and R ~ ~ I v dx = R ~ v dx, for eac macro-cell ~ 2 ~T. By standard arguments (Bramble-Hilbert lemma) we conclude te estimate kv? ~ I vk + 1=2 kv? ~ I C i; 2 kr2 vk ~ ; (4.20) 16
17 wit certain interpolation constants C i; independent of. Now, we coose z = ~ I z in te error representation (4.2) and obtain J(e) = 2T n (f? ~ f ; z? ~ I z )? 1 2 (n [ru ]; z? ~ I z o ; (4.21) were ~ f denotes te piecewise constant interpolation f on te macro-mes ~T. Hence, in view of (4.20), we conclude tat jj(u )j C i wic proves te assertion. 2T 3 krfk + ~ 3=2 kn [ru kr 2 zk ; (4.22) ~ 5 Applications Te following examples are intended to support some of te claims made above on te basis of teoretical analysis. In te tables sown below, we list te results obtained for a sequence of successively reduced values of TOL, te number of points, N opt, of te nal mes T opt, on wic te error indicators are balanced, te maximum number of renement levels, L, in tis mes, and te eciency index I e := jj(e)j=(u ). Example 1. As rst example, we consider te computation of isolated stress values, i.e., point-values of te gradient ru, in te model Poisson problem (1.2) on te square = (?1; 1) 2. Let, for example, te error functional be cosen as J(u) 1 u(0) : Since te corresponding dual solution does not exist in te sense of H 1 0(), we ave to regularize te functional according to J (u) := jb (0)j?1 B 1 u dx ; were B (0) := fx 2 ; jxj < g, and = TOL. Ten, jj(u)? J (u)j TOL, and te corresponding dual solution z is well dened in H0() 1. Observing tat z (x) x 1 =(jxj + ) 2, we are led to te weigted a posteriori error estimate j@ 1 e(0)j C i 2T 2 (r + TOL) 3 + ctol ; (5.2) were r = maxfdist(; 0); g. We again assume tat te local residuals satisfy 2. Ten, on a mes T opt on wic te local error indicators =! are equilibrated, tere olds T OL N opt ; N opt = #f T opt g ; (5.3) T OL (r + TOL) 3 1=4 ; (5.4) N opt 17
18 Consequently, te number N opt of elements of T opt is given by N opt 2T opt 2?2 = 2T opt 2 N 1=2 opt : (5.5) TOL 1=2 3=2 (r + TOL) Since r?3=2 is integrable over, it follows tat N opt TOL?1, as for te evaluation of te point value u(0). Furter, te mes size is distributed like TOL 1=2 (r + TOL) 3=4 ; and, consequently, min TOL 5=4 and max TOL 1=2. We note tat if te functional J() is not regularized, te same reasoning yields te mes-size distribution TOL 1=2 r 3=4 and, in tis case, min TOL 2, meaning a strong over-renement at x = 0. Tis eect is conrmed in our test computations. In contrast to te weigted estimator (5.2), te global energy-error estimator (1.13), for r = 0, or te analogue of te L 1 -error estimator (2.11), wit 1=2 krek C s C i 2 =: energy (u ) ; (5.6) 2T max j@ 1 ej C i C s j log 0 j 1=2 max 2T 1; =: 1 (u ) ; (5.7) 1; := max jf + u j + max jn [ru ]j would generate meses of complexity N opt T OL?2, in order to guarantee (u ) TOL. Figure 2 sows te balanced mes for TOL = 4?4 and te approximation to te dual solution z, = TOL, computed on tis mes. Te corresponding errors for a sequence of tolerances are listed in Table 4. For te problem considered, te weigted a posteriori estimate (2.4) for te point error is obviously asymptotically optimal and te predicted dependence N opt TOL?1 is conrmed. Example 2. integral As a very particular example, we consider te computation of te contour n un 1 ds ; (5.8) for te solution of te model problem (1.2) on te unit disc, = B 1 = fx 2 R 2 ; jxj < 1g, were n 1 is te x 1 -component of te outer unit vector Error functionals of tis type occur, for instance, in te computation of drag and lift coecients of blunt bodies in viscous ows modeled by te Navier-Stokes equations (see [7]). In view of te pointwise analogue of te a priori error estimate (1.6), we ave J(e) = O( max ), so tat convergence is no problem. Te question is ow te mes renement process sould be organized in order to obtain J(u) wit best accuracy. We will see tat te answer is somewat surprising and seems, at te rst glance, to conict wit wat one is used to from a nite element metod. 18
19 Figure 2: Rened mes and approximate dual solution for 1 u(0) in Example 1, using te weigted error estimator point, wit TOL = 4?4 energy TOL N L j@ 1 e(0)j I e e e ? e e ? e e ;1 TOL N L j@ 1 e(0)j I e e e ? e e ? e e ;point TOL N L j@ 1 e(0)j I e e e ? e e ? e e ? e e ? e e ? e e ? e e regular 1;point TOL N L j@ 1 e(0)j I e e e-2-4? e e-2-4? e e-2-4? e e ? e e ? e e ? e e Table 4: Results for 1 u(0) in Example 1, using te energy- and W 1 -error estimators energy and 1;1, respectively, compared to te weigted estimator 1;point and its regularized variant regular 1;point First, in oder to avoid aving to deal wit measures, we regularize te functional J(). Wit te cut-o function ( 1 ; 0 jxj < 1? ; (x) =?1 (1? jxj) ; 1? jxj < 1 ; we make te following ansatz for te dual solution z (x) =?x 1 (x) ; 19
20 were = TOL. Let ' 2 V [ C 2 ( ) be arbitrary. We split te integration over B 1 according to B 1 = B 1? [ S wit S = B 1 n B 1? and obtain, using polar coordinates (r; ), n o rz r' dx n z ' ds r r ' + r B 1 ' dx 1? By construction, tere n z =?@ r x 1 =?cos 1 n B 1?, and consequently, observing x 1 = r cos(), n z ' ds =? cos()'(1? ; ) ds = O() : 0 Furter, we ave and r r ' dx =? Finally, observing tat we nd 1? x r ' dx? r ' dx =? S S?1 x r ' dx + O() ; S S ' dx =?1?1 S x r ' dx = S r?2 (1? r)sin()@ ' dx = O() 1 n n ' ds + O() ; J(') = (rz ; r') + O() : (5.9) Tis implies tat te approximate weigts ~! kr 2 z k in te a posteriori error estimate (2.17) corresponding to te functional J (u ) := (rz ; r') are all zero on elements outside a boundary strip of widt, i.e., tere is no contribution to te error J(e) from tose interior elements. Tis remains true if we let! 0. Terefore, starting from some coarse mes, te \optimal" renement strategy would be in eac step only to rene tose elements wic are adjacent to te boundary and to leave te oters uncanged. Te resulting meses are maximally rened along te boundary down to a minimal mes size min, wile in te interior of te elements are kept of size max 1. If only local mes bisection is allowed in eac renement step, te mes obtained after L steps as minimal size min 2?L and consists of N L 2 L elements. To reac an error of size jj(e)j T OL, about k log 2 (1=T OL) renement cycles are needed leading to an \optimal" mes wit N 1=T OL elements. However, to realize te full accuracy on te constructed mes, te rigt and side consisting of integrals over te force f multiplied by te nodal basis functions as to be calculated almost exactly also on te coarse elements. Figure 3 sows te te nest mes obtained after 7 renement steps and te approximation to te dual solution z ; = min, computed on tis mes. In Table 5 te results obtained by te usual energy error estimator energy (u ) are compared to tose of te \optimal" renement strategy. Example 3. By te tird example, we demonstrate tat carrying weigts in te a posteriori error estimates can also be advantageous in estimating te error in global norms. Tis applies 20
21 Figure 3: Rened mes and approximate dual solution for Example 2 obtained by te approximate weigted error estimator approx ; after 7 renement steps weigt energy L N J(e) : : : :49? :77? :73? 1 6 memory exausted 7 memory exausted \optimal" strategy L N J(e) : : : :48? :62? :67? :93? :03? 2 Table 5: Results obtained wit te energy norm estimator energy strategy for Example 2 and by te \optimal" renement particularly wen global a priori bounds for te dual solution do not appropriately reect te local stability properties of te underlying problem, e.g., for problems wit strongly varying coecients. Consider, for example te boundary value problem? divfarug = f in ; u = 0 ; (5.10) on te square = (?1; 1) 2 wit a coecient function, a(x) = 0:1 + e 3(x 1+x 2) : (5.11) We discretize tis problem by bilinear nite elements and recall te weigted L 2 -a posteriori error estimates following from (1.8), kek 2T! ; (5.12) were te residuals and weigts! are dened by (1.9) and (1.10), were z is te solution 21
22 of te dual problem,? divfarzg = e=kek in ; z = 0 ; (5.13) First, te weigts ~! are computed by taking second dierence quotients of te discrete dual solution z (1) 2 V on te current mes T, wile on te rigt and side of (5.13) te exact error e is used. Te reference solution is as before u = 10 sin(2x 1 + x 2 + 2). Tis gives te approximate stability constant 1=2 ~C s := : 2T ~! 2 Tis way, we obtain a rst approximate weigted L 2 -error estimator ~ weigt (u ) := 2T ~! as well as te corresponding global L 2 -error estimator ~ 2 (u ) := C i ~ Cs 2T 2 2 te interpolation constant being cosen as C i = 0:2. Next, te unknown error e in (5.13) is approximated by ~e := I (2) u? u, were I (2) u is te patc-wise biquadratic interpolation of u. Te resulting error estimator is denoted by ~ approx weigt. Te results obtained for te above test problem by using tese tree L 2 -error estimators are listed in Table 6. Figure 4 sows te errors on te corresponding \optimal" meses for te global estimator ~ 2 (u ) and te weigted estimator ~ weigt (u ). Tese tests clearly demonstrate te superiority of ~ weigt over ~ 2. ; Figure 4: Errors on balanced meses wit N for Example 3 by te L 2 -error estimator ~ 2 (u ) (left, scaled by 1:10) and te weigted error estimator ~ approx weigt (u ) (rigt, scaled by 1:30) Example 4. Te nal example is intended to support our claim tat, in using iger order nite elements, one sould indeed use te approximation property of te nite element space to its full extent in generating te weigts! in te a posteriori estimate (2.17). First, we 22
23 ~ 2 TOL N L kek I e Cs ~ e e ? e e ? e e ? e e ? e e ? e e ~ weigt TOL N L kek I e e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e ~ approx weigt TOL N L kek I e e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e ? e e Table 6: Results obtained for Example 3 by te global L 2 -error estimator ~ 2, compared to te weigted estimators ~ weigt and ~ approx weigt illustrate tis point by an analytical argument. Consider again te model Poisson problem (1.2) on te unit circle = fx 2 R 2 ; jxj < 1g. Suppose tat te mean value J(u) := u(x) dx as to be computed. In tis case te corresponding dual solution is just te quadratic polynomial z(x) = 1? 1 2 jxj2, wic satises?z = 1 and z j@ = 0. We tink tis problem to be solved by biquadratic nite elements (case p = 2). Ten, te \optimal" a posteriori error bound jj(e)j C i 2T 2 kr 3 zk would indicate tat, similar as in Example 2, te error is concentrated along te boundary, wile tere is no contribution from interior mes cells. In contrast, te sub-optimal a posteriori bound jj(e)j C i 2T kr 2 zk would yield contributions to te error from all mes cells. For furter illustration, we consider again te computation of te point value u(0) in te above model situation. Accordingly, we compare te following two a posteriori error estimates, je(0)j C i 3 2T r 3 23 =: (3) point (u ) ; (5.14)
24 and je(0)j C i 2 2T r 2 Te results sown in Table 7 clearly support our claim. =: (2) point (u ) : (5.15) (3) point TOL N L 4? e-1 4? e-2 4? e-3 4? e-3 4? e-4 4? e-4 4? e-5 4? e-5 4? e-6 4? e-7 (2) point TOL N L 4? e-1 4? e-2 4? e-2 4? e-3 4? e-4 4? e-4 4? e-5 4? e-5 4? e-6 4? e-7 Table 7: Results for computing e(0) by quadratic elements in Example 4, obtained wit te local estimators (3) point (u ) and (2) point (u ) 6 Generalization to nonlinear problems Te approac to residual-based error estimation described above can be extended to general nonlinear problems; for ideas in te same direction see [14], [22], [23], [16], as well as [29] and [26]. Related energy-error estimators ave been derived in [29]. We outline our general concept in an abstract setting. Let V be a Hilbert space wit inner product (; ) and corresponding norm k k, and a(; ) a semi-linear form continuously dened on V V. We seek a solution to te abstract variational problem u 2 V : a(u; ') = 0 8 ' 2 V: (6.2) Tis problem is approximated by a Galerkin metod using a sequence of nite dimensional subspaces V V parameterized by a discretization parameter. Te discrete problems read u 2 V : a(u ; ' ) = 0 8 ' 2 V : (6.3) Wit te derivative a 0 (; ; ) of a(; ), we ave te following ortogonality relation for te error e = u? u : 1 0 a 0 (tu + (1? t)u ; e; ' ) dt = a(u; ' )? a(u ; ' ) = 0 8' 2 V : (6.4) Tis suggests te use of te bilinear form L(u; u ; '; z) = 1 0 a 0 (tu + (1? t)u ; '; z) dt ; (6.5) 24
Order of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
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 informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationThe Complexity of Computing the MCD-Estimator
Te Complexity of Computing te MCD-Estimator Torsten Bernolt Lerstul Informatik 2 Universität Dortmund, Germany torstenbernolt@uni-dortmundde Paul Fiscer IMM, Danisc Tecnical University Kongens Lyngby,
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationarxiv: v1 [math.na] 3 Nov 2011
arxiv:.983v [mat.na] 3 Nov 2 A Finite Difference Gost-cell Multigrid approac for Poisson Equation wit mixed Boundary Conditions in Arbitrary Domain Armando Coco, Giovanni Russo November 7, 2 Abstract In
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More information2 Multi-Dimensional Variational Principles y s n u = α θ pu = β n Ω Figure 3..: Two-dimensional region wit and normal vector n x were
Capter 3 Multi-Dimensional Variational Principles 3. Galerkin's Metod and Extremal Principles Te construction of Galerkin formulations presented in Capters and 2 for one-dimensional problems readily extends
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationMath 2921, spring, 2004 Notes, Part 3. April 2 version, changes from March 31 version starting on page 27.. Maps and di erential equations
Mat 9, spring, 4 Notes, Part 3. April version, canges from Marc 3 version starting on page 7.. Maps and di erential equations Horsesoe maps and di erential equations Tere are two main tecniques for detecting
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More information3. Using your answers to the two previous questions, evaluate the Mratio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More information1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More information7 Semiparametric Methods and Partially Linear Regression
7 Semiparametric Metods and Partially Linear Regression 7. Overview A model is called semiparametric if it is described by and were is nite-dimensional (e.g. parametric) and is in nite-dimensional (nonparametric).
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationConvergence and Descent Properties for a Class of Multilevel Optimization Algorithms
Convergence and Descent Properties for a Class of Multilevel Optimization Algoritms Stepen G. Nas April 28, 2010 Abstract I present a multilevel optimization approac (termed MG/Opt) for te solution of
More informationChapter 10. Function approximation Function approximation. The Lebesgue space L 2 (I)
Capter 1 Function approximation We ave studied metods for computing solutions to algebraic equations in te form of real numbers or finite dimensional vectors of real numbers. In contrast, solutions to
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationOptimization of stress modes by energy compatibility for 4-node hybrid quadrilaterals
INTERNATIONAL JOURNAL FOR NUMERIAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 2004; 59:293 33 (DOI: 0.002/nme.877) Optimization of stress modes by energy compatibility for 4-node ybrid quadrilaterals
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationTest 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =
Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationGRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES
GRID CONVERGENCE ERROR ANALYSIS FOR MIXED-ORDER NUMERICAL SCHEMES Cristoper J. Roy Sandia National Laboratories* P. O. Box 5800, MS 085 Albuquerque, NM 8785-085 AIAA Paper 00-606 Abstract New developments
More informationDerivatives. By: OpenStaxCollege
By: OpenStaxCollege Te average teen in te United States opens a refrigerator door an estimated 25 times per day. Supposedly, tis average is up from 10 years ago wen te average teenager opened a refrigerator
More information