Some recent results on positivity-preserving discretisations for the convection-diffusion equation

Size: px
Start display at page:

Download "Some recent results on positivity-preserving discretisations for the convection-diffusion equation"

Transcription

1 Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics, University of Strathclyde, Scotland 2 WIAS Institute, Berlin, Germany 3 Charles University, Prague, Czech Republic Recent advances in Discontinuous Galerkin Methods, University of Reading, Reading, September 11-12, 2014

2 Introduction: The discrete maximum principle The continuous maximum principle : Theorem Let u be the solution of the problem ε u + b u = f in Ω, and u = 0 on Ω. Then, if f 0 in Ω, then u 0 in Ω, and attains its minimum at the boundary. G.R. Barrenechea (Strathclyde) Reading, September / 24

3 The convection-diffusion equation The DMP : Theorem Let u h P 1 (Ω) be the solution of the problem ε ( u h, v h ) Ω + (b u h, v h ) Ω = (f, v h ) Ω v h P 1 (Ω). Then, if f 0 in Ω, the mesh is acute, and b h 2ε < 1, then u h 0 in Ω, and attains its minimum at the boundary. Remark : Under these hypothesis, the matrix [ε( λ j, λ i ) Ω + (b λ j, λ i ) Ω ] is an M-matrix. This is, it is invertible, all the diagonal elements are positive, and the off-diagonal ones are non-positive. G.R. Barrenechea (Strathclyde) Reading, September / 24

4 The convection-diffusion equation The DMP : Theorem Let u h P 1 (Ω) be the solution of the problem ε ( u h, v h ) Ω + (b u h, v h ) Ω = (f, v h ) Ω v h P 1 (Ω). Then, if f 0 in Ω, the mesh is acute, and b h 2ε < 1, then u h 0 in Ω, and attains its minimum at the boundary. Remark : Under these hypothesis, the matrix [ε( λ j, λ i ) Ω + (b λ j, λ i ) Ω ] is an M-matrix. This is, it is invertible, all the diagonal elements are positive, and the off-diagonal ones are non-positive. G.R. Barrenechea (Strathclyde) Reading, September / 24

5 Some early solutions Artificial diffusion : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + s(u h, v h ) = (f, v h ) Ω v h P 1 (Ω). Bad news : The linear schemes, such as the artificial diffusion, have two main drawbacks: - their consistency error leads to a convergence of O( h); - they produce results which are extremely diffusive. G.R. Barrenechea (Strathclyde) Reading, September / 24

6 Some early solutions Artificial diffusion : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + α h ( u h, v h ) Ω = (f, v h ) Ω v h P 1 (Ω). Bad news : The linear schemes, such as the artificial diffusion, have two main drawbacks: - their consistency error leads to a convergence of O( h); - they produce results which are extremely diffusive. G.R. Barrenechea (Strathclyde) Reading, September / 24

7 Some early solutions Artificial diffusion : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + α h ( u h, v h ) Ω = (f, v h ) Ω v h P 1 (Ω). Bad news : The linear schemes, such as the artificial diffusion, have two main drawbacks: - their consistency error leads to a convergence of O( h); - they produce results which are extremely diffusive. G.R. Barrenechea (Strathclyde) Reading, September / 24

8 Some early solutions Artificial diffusion : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + α h ( u h, v h ) Ω = (f, v h ) Ω v h P 1 (Ω). Bad news : The linear schemes, such as the artificial diffusion, have two main drawbacks: - their consistency error leads to a convergence of O( h); - they produce results which are extremely diffusive. G.R. Barrenechea (Strathclyde) Reading, September / 24

9 Solution: nonlinear schemes Idea : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + N(u h ; u h, v h ) = (f, v h ) Ω v h P 1 (Ω). Main features : - N is a continuous form, may depend on the residual, or not. - In some cases (not that many!), the maximum principle can be proved (cf. Burman & Ern). - Optimal convergence can be proved in most cases. A more recent alternative (D. Kuzmin) : Algebraic Flux Correction schemes. These work at the matrix level, and have provided very convincing numerical results. G.R. Barrenechea (Strathclyde) Reading, September / 24

10 Solution: nonlinear schemes Idea : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + N(u h ; u h, v h ) = (f, v h ) Ω v h P 1 (Ω). Main features : - N is a continuous form, may depend on the residual, or not. - In some cases (not that many!), the maximum principle can be proved (cf. Burman & Ern). - Optimal convergence can be proved in most cases. A more recent alternative (D. Kuzmin) : Algebraic Flux Correction schemes. These work at the matrix level, and have provided very convincing numerical results. G.R. Barrenechea (Strathclyde) Reading, September / 24

11 Solution: nonlinear schemes Idea : Find u h P 1 (Ω) such that ε ( u h, v h ) Ω + (b u h, v h ) Ω + N(u h ; u h, v h ) = (f, v h ) Ω v h P 1 (Ω). Main features : - N is a continuous form, may depend on the residual, or not. - In some cases (not that many!), the maximum principle can be proved (cf. Burman & Ern). - Optimal convergence can be proved in most cases. A more recent alternative (D. Kuzmin) : Algebraic Flux Correction schemes. These work at the matrix level, and have provided very convincing numerical results. G.R. Barrenechea (Strathclyde) Reading, September / 24

12 Goals and Outline 1 Goals: Understand the method, and its main features. Give the first steps towards a numerical analysis of it. Study its numerical behaviour. 2 The method for the 1D problem. 3 The discrete maximum principle. 4 Solvability of the linear problems, and the nonlinear one. 5 Concluding remarks. G.R. Barrenechea (Strathclyde) Reading, September / 24

13 Goals and Outline 1 Goals: Understand the method, and its main features. Give the first steps towards a numerical analysis of it. Study its numerical behaviour. 2 The method for the 1D problem. 3 The discrete maximum principle. 4 Solvability of the linear problems, and the nonlinear one. 5 Concluding remarks. G.R. Barrenechea (Strathclyde) Reading, September / 24

14 Algebraic flux correction schemes Starting point : A finite element discretisation of our problem of the form: AU = G. Define: D := (d ij ) where d ij := max{a ij, 0, a ji } for i j, d ii = j i d ij. Remark: The matrix à is an M-matrix. Then, it preserves positivity. G.R. Barrenechea (Strathclyde) Reading, September / 24

15 Algebraic flux correction schemes Starting point : A finite element discretisation of our problem of the form: AU = G. Define: D := (d ij ) where d ij := max{a ij, 0, a ji } for i j, d ii = j i d ij. Remark: The matrix à is an M-matrix. Then, it preserves positivity. G.R. Barrenechea (Strathclyde) Reading, September / 24

16 Algebraic flux correction schemes Starting point : A finite element discretisation of our problem of the form: ( A + D)U = G + DU. Define: D := (d ij ) where d ij := max{a ij, 0, a ji } for i j, d ii = j i d ij. Remark: The matrix à is an M-matrix. Then, it preserves positivity. G.R. Barrenechea (Strathclyde) Reading, September / 24

17 Algebraic flux correction schemes Starting point : A finite element discretisation of our problem of the form: ( A + D )U = G + DU. }{{} =:à Define: D := (d ij ) where d ij := max{a ij, 0, a ji } for i j, d ii = j i d ij. Remark: The matrix à is an M-matrix. Then, it preserves positivity. G.R. Barrenechea (Strathclyde) Reading, September / 24

18 Algebraic flux correction schemes Starting point : A finite element discretisation of our problem of the form: ( A + D )U = G + DU. }{{} =:à Define: D := (d ij ) where d ij := max{a ij, 0, a ji } for i j, d ii = j i d ij. Remark: The matrix à is an M-matrix. Then, it preserves positivity. G.R. Barrenechea (Strathclyde) Reading, September / 24

19 Algebraic flux correction schemes Equivalent system : Ã U = G + DU. From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

20 Algebraic flux correction schemes Equivalent system : Ã U = G + DU. From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

21 Algebraic flux correction schemes Equivalent system : (Ã U)i = g i + From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i j i f ij Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

22 Algebraic flux correction schemes Equivalent system : (Ã U)i = g i + From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i j i f ij Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

23 Algebraic flux correction schemes Equivalent system : (Ã U)i = g i + α ij (U)f ij From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i j i Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

24 Algebraic flux correction schemes Equivalent system : (Ã U)i = g i + α ij (U)f ij From the properties of D it follows that ( ) DU i = f ij where f ij = d ij (u j u i ) are the fluxes. j i j i Goal : To limit the fluxes f ij which are responsible for spurious oscillations. The limiters α ij should satisfy the following: - α ij [0, 1]; - α ij should be as close to 1 as possible; - α ij 1 where the Galerkin solution is smooth. G.R. Barrenechea (Strathclyde) Reading, September / 24

25 Definition of the limiters 1 Compute P + i, P i, Q+ i, Q i in such a way that, for each pair of neighbouring nodes x i, x j with indices such that a ji a ij one performs the updates P + i := P + i + max{0, f ij }, P i := P i max{0, f ji }, Q + i := Q + i + max{0, f ji }, Q i := Q i max{0, f ij }, Q + j := Q+ j + max{0, f ij}, Q j := Q j max{0, f ji}, 2 Set R + i { := min 1, Q+ i P + i } {, R i := min 1, Q i P i }. 3 Finally, α ij = { R + i if f ij > 0, R i if f ij < 0, i, j = 1,..., N. G.R. Barrenechea (Strathclyde) Reading, September / 24

26 The 1D convection-diffusion equation Model problem : ε u + bu = g in (0, 1) u(0) = u(1) = 0, with positive constants ε and b. Galerkin FEM : Equidistant nodes x i = ih, with h = 1/N. Find u h P 1 (0, 1) such that u h (0) = u h (1) = 0 and ε(u h, v h) + (bu h, v h ) = (g, v h ) v h P 1 (0, 1). Difference equation form : Setting u i = u h (x i ), this problem is rewritten as ε u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i i = 1,..., N 1. Assume: P e := bh 2ε > 1. G.R. Barrenechea (Strathclyde) Reading, September / 24

27 The 1D convection-diffusion equation Model problem : ε u + bu = g in (0, 1) u(0) = u(1) = 0, with positive constants ε and b. Galerkin FEM : Equidistant nodes x i = ih, with h = 1/N. Find u h P 1 (0, 1) such that u h (0) = u h (1) = 0 and ε(u h, v h) + (bu h, v h ) = (g, v h ) v h P 1 (0, 1). Difference equation form : Setting u i = u h (x i ), this problem is rewritten as ε u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i i = 1,..., N 1. Assume: P e := bh 2ε > 1. G.R. Barrenechea (Strathclyde) Reading, September / 24

28 The 1D convection-diffusion equation Model problem : ε u + bu = g in (0, 1) u(0) = u(1) = 0, with positive constants ε and b. Galerkin FEM : Equidistant nodes x i = ih, with h = 1/N. Find u h P 1 (0, 1) such that u h (0) = u h (1) = 0 and ε(u h, v h) + (bu h, v h ) = (g, v h ) v h P 1 (0, 1). Difference equation form : Setting u i = u h (x i ), this problem is rewritten as ε u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i i = 1,..., N 1. Assume: P e := bh 2ε > 1. G.R. Barrenechea (Strathclyde) Reading, September / 24

29 The 1D convection-diffusion equation Model problem : ε u + bu = g in (0, 1) u(0) = u(1) = 0, with positive constants ε and b. Galerkin FEM : Equidistant nodes x i = ih, with h = 1/N. Find u h P 1 (0, 1) such that u h (0) = u h (1) = 0 and ε(u h, v h) + (bu h, v h ) = (g, v h ) v h P 1 (0, 1). Difference equation form : Setting u i = u h (x i ), this problem is rewritten as ε u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i i = 1,..., N 1. Assume: P e := bh 2ε > 1. G.R. Barrenechea (Strathclyde) Reading, September / 24

30 The 1D convection-diffusion equation Algebraic problem with limited fluxes: ( AU ) i + j i (1 α ij )f ij = g i with f ij = d ij (u j u i ). For the 1D problem: the system reduces to u 0 = u N = 0, and (ε + β i ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i, i = 1,..., N 1, where 1 if u i+1 u i and β i = 0 otherwise, u i u i 1 u i+1 u i < 1, and ε = b h 2 ε = ε (P e 1). G.R. Barrenechea (Strathclyde) Reading, September / 24

31 The Discrete Maximum Principle Theorem Consider any ε b h/2 ε. Then any solution of the nonlinear problem satisfies the discrete maximum principle, i.e., for any i {1,..., N}, one has g i 0 u i min{u i 1, u i+1 }. Moreover, for any k, l {0, 1,..., N + 1} with k + 1 < l, one has g i 0, i = k + 1,..., l 1 u i min{u k, u l }, i = k,..., l. G.R. Barrenechea (Strathclyde) Reading, September / 24

32 Some numerics and the choice of ε Other possible choices: The artificial diffusion matrix D can be defined using different combinations of the diffusion and convection matrices. For example: (F) ε = b h 2 (C) ε = b h 2. ε = ε(p e 1). ( coth P e 1 P e). (P) ε = b h 2 Data: b = f = 1, N = 16, ε = 0.03, i.e., we solve and u(0) = u(1) = u + u = 1 in (0, 1), G.R. Barrenechea (Strathclyde) Reading, September / 24

33 Some numerics and the choice of ε Figure 1 : from (P). Comparison of the exact solution (green) and discrete solution with ε G.R. Barrenechea (Strathclyde) Reading, September / 24

34 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

35 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

36 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

37 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

38 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

39 Bad news from the numerics - Computations very sensitive to rounding errors. Idea : replace the condition u i < min{u i 1, u i+1 } by u i < min{u i 1, u i+1 } τ. - Not a remedy! Conclusion: The nonlinear problem is not solvable in general! Example: N = 4, ε = 0.03, b = 1, f 1 = 6, f 2 = 6, f 3 = 3, f 4 = 2, and ε from (F). Reminder of the problem: (ε + β i (u) ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i. G.R. Barrenechea (Strathclyde) Reading, September / 24

40 Bad news from the numerics G.R. Barrenechea (Strathclyde) Reading, September / 24

41 Solvability of the linear subproblems Theorem For every choice of ε [ b h 2 ε, ] b h 2, and every possible βi [0, 1], the problem has a unique solution. (ε + β i ε) u i 1 2 u i + u i+1 h 2 + b u i+1 u i 1 2 h = g i, G.R. Barrenechea (Strathclyde) Reading, September / 24

42 Solvability of the nonlinear problem Main remark : The lack of solvability is due to the discontinuity of the coefficients β i Theorem Let us suppose that the functions β i : R N+1 [0, 1], i = 1,..., N 1, are continuous, and let ε be any of the previous choices. Then, the nonlinear FCT scheme has a solution. Proof: Write the method as the fixed point equation M(β(u)) u = g, apply the fact that the determinant is a continuous function of the entries of a matrix, and Brouwer s fixed point Theorem. G.R. Barrenechea (Strathclyde) Reading, September / 24

43 Solvability of the nonlinear problem Main remark : The lack of solvability is due to the discontinuity of the coefficients β i Theorem Let us suppose that the functions β i : R N+1 [0, 1], i = 1,..., N 1, are continuous, and let ε be any of the previous choices. Then, the nonlinear FCT scheme has a solution. Proof: Write the method as the fixed point equation M(β(u)) u = g, apply the fact that the determinant is a continuous function of the entries of a matrix, and Brouwer s fixed point Theorem. G.R. Barrenechea (Strathclyde) Reading, September / 24

44 Graphical representation of the regularisation G.R. Barrenechea (Strathclyde) Reading, September / 24

45 The price to pay: A weak version of the DMP Theorem Let u 0,..., u N+1 be a solution of the modified FCT scheme with any functions β 1,..., β N [0, 1] as described before. Then g i 0 u i min{u i 1, u i+1 } or u i max{u i 1, u i+1 } δ h, for i = 1,..., N. G.R. Barrenechea (Strathclyde) Reading, September / 24

46 Numerical evidence on the violation of the DMP The problem : εu + u = 0 subject to u(0) = 1 and u(1) = 0. We measured MAX := u max h 1; RMAX := max{(u max h 1)/h}; P e RMAX the value of P e for which the maximum RMAX is attained. Table 1 : Violation of the discrete maximum principle for the continuous β i. P e [1, 20) P e [20, ) ε MAX RMAX P e RMAX MAX RMAX P e RMAX no P e no P e G.R. Barrenechea (Strathclyde) Reading, September / 24

47 Some preliminary numerics in 2D: The Hemker problem Data: ε = 10 4, 12, 000 Q 1 elements, discontinuous α ij as before, continuous as follows { R + i = R min{q + i i = min 1,, Q i } } max{p + i, P i, τ}. Figure 2 : Discontinuous α ij, non-symmetric G.R. Barrenechea (Strathclyde) Reading, September / 24

48 Some preliminary numerics in 2D: The Hemker problem Figure 3 : Continuous α ij, non-symmetric G.R. Barrenechea (Strathclyde) Reading, September / 24

49 Conclusions and perspectives 1 Some further insight on FCT schemes. 2 Analysis of a wider class of schemes. 3 Counter-examples of existence of solutions for the original method. 4 A modification that is proved to possess solutions, but satisfies only a weak version of the DMP. Future extensions: Deeper study of the symmetric version in higher dimensions. Maximum principle on general meshes. (Order of) convergence. Time-dependent problems. Coupled nonlinear problems in chemical reactions. G.R. Barrenechea (Strathclyde) Reading, September / 24

50 Conclusions and perspectives 1 Some further insight on FCT schemes. 2 Analysis of a wider class of schemes. 3 Counter-examples of existence of solutions for the original method. 4 A modification that is proved to possess solutions, but satisfies only a weak version of the DMP. Future extensions: Deeper study of the symmetric version in higher dimensions. Maximum principle on general meshes. (Order of) convergence. Time-dependent problems. Coupled nonlinear problems in chemical reactions. G.R. Barrenechea (Strathclyde) Reading, September / 24

A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations

A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations Noname manuscript No. (will be inserted by the editor) A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations Gabriel R. Barrenechea Volker John Petr Knobloch Richard

More information

NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE

NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING

More information

On discontinuity capturing methods for convection diffusion equations

On discontinuity capturing methods for convection diffusion equations On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,

More information

Singularly Perturbed Partial Differential Equations

Singularly Perturbed Partial Differential Equations WDS'9 Proceedings of Contributed Papers, Part I, 54 59, 29. ISN 978-8-7378--9 MTFYZPRESS Singularly Perturbed Partial Differential Equations J. Lamač Charles University, Faculty of Mathematics and Physics,

More information

ALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS

ALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION

More information

Universität des Saarlandes. Fachrichtung 6.1 Mathematik

Universität des Saarlandes. Fachrichtung 6.1 Mathematik Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations

More information

Algebraic flux correction and its application to convection-dominated flow. Matthias Möller

Algebraic flux correction and its application to convection-dominated flow. Matthias Möller Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany

More information

On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion

On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11

More information

On the design of higher-order FEM satisfying the discrete maximum principle

On the design of higher-order FEM satisfying the discrete maximum principle On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany

More information

On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P1 and Q1 finite elements

On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P1 and Q1 finite elements Volker John, Petr Knobloch On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P and Q finite elements Preprint No. MATH-knm-27/4 27.

More information

On solving linear systems arising from Shishkin mesh discretizations

On solving linear systems arising from Shishkin mesh discretizations On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October

More information

Multigrid Methods and their application in CFD

Multigrid Methods and their application in CFD Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential

More information

Finite difference method for elliptic problems: I

Finite difference method for elliptic problems: I Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen

More information

High-resolution finite element schemes for an idealized Z-pinch implosion model

High-resolution finite element schemes for an idealized Z-pinch implosion model High-resolution finite element schemes for an idealized Z-pinch implosion model Matthias Möller a, Dmitri Kuzmin b, John N. Shadid c a Institute of Applied Mathematics (LS III), Dortmund University of

More information

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,

[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of, 269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,

More information

The CG1-DG2 method for conservation laws

The CG1-DG2 method for conservation laws for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity

More information

Compressed Sensing and Neural Networks

Compressed Sensing and Neural Networks and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications

More information

RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS

RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.

More information

FEM-Level Set Techniques for Multiphase Flow --- Some recent results

FEM-Level Set Techniques for Multiphase Flow --- Some recent results FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3

More information

Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?

Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal

More information

Basic Principles of Weak Galerkin Finite Element Methods for PDEs

Basic 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 information

Cholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems

Cholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems Cholesky factorisations of linear systems coming from a finite difference method applied to singularly perturbed problems Thái Anh Nhan and Niall Madden The Boundary and Interior Layers - Computational

More information

DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS

DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several

More information

Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, Berlin, Germany,

Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, Berlin, Germany, Volker John On the numerical simulation of population balance systems Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany, Free University of Berlin, Department

More information

Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé

Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques

More information

Universität des Saarlandes. Fachrichtung 6.1 Mathematik

Universität des Saarlandes. Fachrichtung 6.1 Mathematik Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 156 A comparison of spurious oscillations at layers diminishing (SOLD) methods for convection

More information

On 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 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 information

WELL POSEDNESS OF PROBLEMS I

WELL POSEDNESS OF PROBLEMS I Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according

More information

Algorithms for Scientific Computing

Algorithms for Scientific Computing Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling

More information

In Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint

In Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION

More information

Discontinuous Galerkin Methods

Discontinuous Galerkin Methods Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.

More information

Kasetsart University Workshop. Multigrid methods: An introduction

Kasetsart University Workshop. Multigrid methods: An introduction Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available

More information

Introduction to Partial Differential Equations

Introduction to Partial Differential Equations Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These

More information

Explicit and implicit FEM-FCT algorithms with flux linearization

Explicit and implicit FEM-FCT algorithms with flux linearization Explicit and implicit FEM-FCT algorithms with flux linearization Dmitri Kuzmin Institute of Applied Mathematics (LS III), Dortmund University of Technology Vogelpothsweg 87, D-44227, Dortmund, Germany

More information

On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection

On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection D. Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund,

More information

On the impact of the scheme for solving the higher dimensional equation in coupled population balance systems

On the impact of the scheme for solving the higher dimensional equation in coupled population balance systems INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2010; 82:1450 1474 Published online 8 February 2010 in Wiley InterScience (www.interscience.wiley.com)..2830 On the

More information

Solution Methods. Steady convection-diffusion equation. Lecture 05

Solution Methods. Steady convection-diffusion equation. Lecture 05 Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate

More information

A Hybrid Method for the Wave Equation. beilina

A Hybrid Method for the Wave Equation.   beilina A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)

More information

Finite Elements for Nonlinear Problems

Finite Elements for Nonlinear Problems Finite Elements for Nonlinear Problems Computer Lab 2 In this computer lab we apply finite element method to nonlinear model problems and study two of the most common techniques for solving the resulting

More information

Implicit FEM-FCT algorithms and discrete Newton methods for transient convection problems

Implicit FEM-FCT algorithms and discrete Newton methods for transient convection problems Implicit FEM-FCT algorithms and discrete Newton methods for transient convection problems M. Möller 1,,, D. Kuzmin 1, and D. Kourounis 2 1 Institute of Applied Mathematics (LS III), University of Dortmund

More information

Generalised Summation-by-Parts Operators and Variable Coefficients

Generalised Summation-by-Parts Operators and Variable Coefficients Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation

More information

Numerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements

Numerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna

More information

Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner

Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu

More information

Reduction of Finite Element Models of Complex Mechanical Components

Reduction 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 information

AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends

AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary

More information

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.

Lecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C. Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference

More information

Intrinsic products and factorizations of matrices

Intrinsic products and factorizations of matrices Available online at www.sciencedirect.com Linear Algebra and its Applications 428 (2008) 5 3 www.elsevier.com/locate/laa Intrinsic products and factorizations of matrices Miroslav Fiedler Academy of Sciences

More information

Unified A Posteriori Error Control for all Nonstandard Finite Elements 1

Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Unified A Posteriori Error Control for all Nonstandard Finite Elements 1 Martin Eigel C. Carstensen, C. Löbhard, R.H.W. Hoppe Humboldt-Universität zu Berlin 19.05.2010 1 we know of Guidelines for Applicants

More information

Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems

Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.

More information

Parameter robust methods for second-order complex-valued reaction-diffusion equations

Parameter robust methods for second-order complex-valued reaction-diffusion equations Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,

More information

ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017

ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017 ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM Prof. Peter Bermel February 6, 2017 Outline Static Equilibrium Dynamic Equilibrium Thermal transport mechanisms Thermal transport modeling

More information

SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS

SECOND 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 information

Discretization of Convection Diffusion type equation

Discretization of Convection Diffusion type equation Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty

More information

Lecture 4.2 Finite Difference Approximation

Lecture 4.2 Finite Difference Approximation Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by

More information

ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS

ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS ON THE CONSERVATION OF MASS AND ENERGY IN HYGROTHERMAL NUMERICAL SIMULATION WITH COMSOL MULTIPHYSICS Michele Bianchi Janetti, Fabian Ochs, and Wolfgang Feist,2 Unit for Energy Efficient Buildings, University

More information

Consistent Positive Directional Splitting of Anisotropic Diffusion

Consistent Positive Directional Splitting of Anisotropic Diffusion In: Boštjan Likar (ed.): Proc. of Computer Vision Winter Workshop, Bled, Slovenia, Feb. 7-9, 2001, pp. 37-48. Consistent Positive Directional Splitting of Anisotropic Diffusion Pavel Mrázek and Mirko Navara

More information

PREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO

PREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS

More information

Relevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]

Relevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE] Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential

More information

March 2002, December Introduction. We investigate the facial structure of the convex hull of the mixed integer knapsack set

March 2002, December Introduction. We investigate the facial structure of the convex hull of the mixed integer knapsack set ON THE FACETS OF THE MIXED INTEGER KNAPSACK POLYHEDRON ALPER ATAMTÜRK Abstract. We study the mixed integer knapsack polyhedron, that is, the convex hull of the mixed integer set defined by an arbitrary

More information

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA

RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Discussiones Mathematicae General Algebra and Applications 23 (2003 ) 125 137 RANK AND PERIMETER PRESERVER OF RANK-1 MATRICES OVER MAX ALGEBRA Seok-Zun Song and Kyung-Tae Kang Department of Mathematics,

More information

Linearization at equilibrium points

Linearization at equilibrium points TMA4165 2007 Linearization at equilibrium points Harald Hanche-Olsen hanche@math.ntnu.no This note is about the behaviour of a nonlinear autonomous system ẋ = f (x) (where x(t) R n ) near an equilibrium

More information

MAT 1332: CALCULUS FOR LIFE SCIENCES. Contents. 1. Review: Linear Algebra II Vectors and matrices Definition. 1.2.

MAT 1332: CALCULUS FOR LIFE SCIENCES. Contents. 1. Review: Linear Algebra II Vectors and matrices Definition. 1.2. MAT 1332: CALCULUS FOR LIFE SCIENCES JING LI Contents 1 Review: Linear Algebra II Vectors and matrices 1 11 Definition 1 12 Operations 1 2 Linear Algebra III Inverses and Determinants 1 21 Inverse Matrices

More information

arxiv: v1 [math.na] 22 Nov 2018

arxiv: v1 [math.na] 22 Nov 2018 Asymptotic preserving Deferred Correction Residual Distribution schemes Rémi Abgrall and Davide Torlo arxiv:1811.09284v1 [math.na] 22 Nov 2018 Abstract This work aims to extend the residual distribution

More information

Residual Distribution. basics, recents developments, relations with other techniques

Residual Distribution. basics, recents developments, relations with other techniques Residual Distribution basics, recents developments, relations with other techniques MARIO RICCHIUTO July 19, 2012 INTRODUCTION with historical perspective t u+ F(u) = 0 In the 80 s numerical techniques

More information

On angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic

On angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen

More information

On the diminishing of spurious oscillations in explicit finite element analysis of linear and non-linear wave propagation and contact problems

On the diminishing of spurious oscillations in explicit finite element analysis of linear and non-linear wave propagation and contact problems 11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16315 On the diminishing of spurious oscillations

More information

LINEAR INTERVAL INEQUALITIES

LINEAR INTERVAL INEQUALITIES LINEAR INTERVAL INEQUALITIES Jiří Rohn, Jana Kreslová Faculty of Mathematics and Physics, Charles University Malostranské nám. 25, 11800 Prague, Czech Republic 1.6.1993 Abstract We prove that a system

More information

Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS

Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general

More information

Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems

Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear

More information

Weierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN

Weierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 SUPG reduced order models for convection-dominated convection-diffusion-reaction

More information

A parameter-free smoothness indicator for high-resolution finite element schemes

A parameter-free smoothness indicator for high-resolution finite element schemes Cent. Eur. J. Math. 11(8) 2013 1478-1488 DOI: 10.2478/s11533-013-0254-4 Central European Journal of Mathematics A parameter-free smoothness indicator for high-resolution finite element schemes Research

More information

Discrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem

Discrete 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 information

Moments, Model Reduction and Nonlinearity in Solving Linear Algebraic Problems

Moments, Model Reduction and Nonlinearity in Solving Linear Algebraic Problems Moments, Model Reduction and Nonlinearity in Solving Linear Algebraic Problems Zdeněk Strakoš Charles University, Prague http://www.karlin.mff.cuni.cz/ strakos 16th ILAS Meeting, Pisa, June 2010. Thanks

More information

EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008

EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008 EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,

More information

A Multiscale DG Method for Convection Diffusion Problems

A Multiscale DG Method for Convection Diffusion Problems A Multiscale DG Method for Convection Diffusion Problems Annalisa Buffa Istituto di Matematica Applicata e ecnologie Informatiche - Pavia National Research Council Italy Joint project with h. J.R. Hughes,

More information

On several aspects of finite element approximations of laminar/turbulent transition models

On several aspects of finite element approximations of laminar/turbulent transition models EPJ Web of Conferences 67 23 (24) DOI:.5/ epjconf/ 246723 C Owned by the authors published by EDP Sciences 24 On several aspects of finite element approximations of laminar/turbulent transition models

More information

Electronic Transactions on Numerical Analysis Volume 32, 2008

Electronic Transactions on Numerical Analysis Volume 32, 2008 Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the

More information

SUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations

SUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations SUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations Swetlana Giere a,, Traian Iliescu b,1, Volker John a,c, David Wells b,1, a Weierstrass Institute for Applied Analysis

More information

Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems

Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems INTERNATIONAL CENTER FOR NUMERICAL METHODS IN ENGINEERING Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems E. Oñate M. Manzan Publication CIMNE Nº-183, February 000

More information

An Introduction to Variational Inequalities

An Introduction to Variational Inequalities An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan

More information

Matching moments and matrix computations

Matching moments and matrix computations Matching moments and matrix computations Jörg Liesen Technical University of Berlin and Petr Tichý Czech Academy of Sciences and Zdeněk Strakoš Charles University in Prague and Czech Academy of Sciences

More information

NUMERICAL ASPECTS OF A BACTERIA GROWTH MODEL

NUMERICAL ASPECTS OF A BACTERIA GROWTH MODEL Proceedings of ALGORITMY 2005 pp. 175 184 NUMERICAL ASPECTS OF A BACTERIA GROWTH MODEL ROBERT STRAKA AND ZDENĚK ČULÍK Abstract. We numerically investigate the bacteria growth model proposed by Mimura et

More information

An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84

An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A

More information

First-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012

First-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012 First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via

More information

Bounds on the radius of nonsingularity

Bounds on the radius of nonsingularity Towards bounds on the radius of nonsingularity, Milan Hladík Department of Applied Mathematics, Charles University, Prague and Institute of Computer Science of the Academy of Czech Republic SCAN 2014,

More information

Numerical Solution I

Numerical Solution I Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions

More information

30 crete maximum principle, which all imply the bound-preserving property. But most

30 crete maximum principle, which all imply the bound-preserving property. But most 3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical

More information

Aspects of Multigrid

Aspects of Multigrid Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)

More information

Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model

Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model Failsafe FCT algorithms for strongly time-dependent flows with application to an idealized Z-pinch implosion model Matthias Möller 1 joint work with: Dmitri Kuzmin 2 1 Institute of Applied Mathematics

More information

arxiv: v1 [physics.comp-ph] 10 Aug 2015

arxiv: v1 [physics.comp-ph] 10 Aug 2015 Numerical experiments on the efficiency of local grid refinement based on truncation error estimates Alexandros Syrakos a,, Georgios Efthimiou a, John G. Bartzis a, Apostolos Goulas b arxiv:1508.02345v1

More information

On the interplay between discretization and preconditioning of Krylov subspace methods

On the interplay between discretization and preconditioning of Krylov subspace methods On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy

More information

Numerical Programming I (for CSE)

Numerical Programming I (for CSE) Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let

More information

Graduate Mathematical Economics Lecture 1

Graduate Mathematical Economics Lecture 1 Graduate Mathematical Economics Lecture 1 Yu Ren WISE, Xiamen University September 23, 2012 Outline 1 2 Course Outline ematical techniques used in graduate level economics courses Mathematics for Economists

More information

Undergraduate Mathematical Economics Lecture 1

Undergraduate Mathematical Economics Lecture 1 Undergraduate Mathematical Economics Lecture 1 Yu Ren WISE, Xiamen University September 15, 2014 Outline 1 Courses Description and Requirement 2 Course Outline ematical techniques used in economics courses

More information

JACOBI S ITERATION METHOD

JACOBI S ITERATION METHOD ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes

More information

Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids

Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität

More information

Time-dependent variational forms

Time-dependent variational forms Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION

More information

On Multigrid for Phase Field

On Multigrid for Phase Field On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis

More information

Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing

Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Entropy-stable discontinuous Galerkin nite element method with streamline diusion and shock-capturing Siddhartha Mishra Seminar for Applied Mathematics ETH Zurich Goal Find a numerical scheme for conservation

More information

PREPRINT 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 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 information