Some recent results on positivity-preserving discretisations for the convection-diffusion equation
|
|
- Thomas Johnston
- 6 years ago
- Views:
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
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 informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationSingularly 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 informationALGEBRAIC 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationAlgebraic 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 informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationOn the design of higher-order FEM satisfying the discrete maximum principle
On the design of higher-order FEM satisfying the discrete maximum principle Dmitri Kuzmin Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationOn spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P1 and Q1 finite elements
Volker John, Petr Knobloch On spurious oscillations at layers diminishing (SOLD) methods for convection-diffusion equations: Part II - Analysis for P and Q finite elements Preprint No. MATH-knm-27/4 27.
More informationOn 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 informationMultigrid 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 informationFinite 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 informationHigh-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,
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 informationThe 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 informationCompressed 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 informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationFEM-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 informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationCholesky 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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationWeierstrass 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 informationUne 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 156 A comparison of spurious oscillations at layers diminishing (SOLD) methods for convection
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationAlgorithms 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 informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationKasetsart 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 informationIntroduction 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 informationExplicit 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 informationOn 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 informationOn 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 informationSolution 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 informationA 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 informationFinite 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 informationImplicit 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 informationGeneralised 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 informationNumerical 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 informationConvergence 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 informationReduction of Finite Element Models of Complex Mechanical Components
Reduction of Finite Element Models of Complex Mechanical Components Håkan Jakobsson Research Assistant hakan.jakobsson@math.umu.se Mats G. Larson Professor Applied Mathematics mats.larson@math.umu.se Department
More informationAMS 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 informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationIntrinsic 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 informationUnified 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 informationComparative 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 informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationECE 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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationDiscretization 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 informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationON 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 informationConsistent 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 informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationRelevant 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 informationMarch 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 informationRANK 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 informationLinearization 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 informationMAT 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 informationarxiv: 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 informationResidual 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 informationOn angle conditions in the finite element method. Institute of Mathematics, Academy of Sciences Prague, Czech Republic
On angle conditions in the finite element method Michal Křížek Institute of Mathematics, Academy of Sciences Prague, Czech Republic Joint work with Jan Brandts (University of Amsterdam), Antti Hannukainen
More informationOn 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 informationLINEAR 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 informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationAdaptive 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 informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 SUPG reduced order models for convection-dominated convection-diffusion-reaction
More informationA 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 informationDiscrete Maximum Principle for a 1D Problem with Piecewise-Constant Coefficients Solved by hp-fem
The University of Texas at El Paso Department of Mathematical Sciences Research Reports Series El Paso, Texas Research Report No. 2006-10 Discrete Maximum Principle for a 1D Problem with Piecewise-Constant
More informationMoments, 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 informationEINDHOVEN 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 informationA 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 informationOn 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 informationElectronic Transactions on Numerical Analysis Volume 32, 2008
Electronic Transactions on Numerical Analysis Volume 32, 2008 Contents 1 On the role of boundary conditions for CIP stabilization of higher order finite elements. Friedhelm Schieweck. We investigate the
More informationSUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations
SUPG Reduced Order Models for Convection-Dominated Convection-Diffusion-Reaction Equations Swetlana Giere a,, Traian Iliescu b,1, Volker John a,c, David Wells b,1, a Weierstrass Institute for Applied Analysis
More informationStabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems
INTERNATIONAL CENTER FOR NUMERICAL METHODS IN ENGINEERING Stabilization Techniques for Finite Element Analysis of Convection-Diffusion Problems E. Oñate M. Manzan Publication CIMNE Nº-183, February 000
More informationAn 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 informationMatching 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 informationNUMERICAL 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 informationAn 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 informationFirst-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 informationBounds 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 informationNumerical 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 information30 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 informationAspects 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 informationFailsafe 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 informationarxiv: 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 informationOn 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 informationNumerical 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 informationGraduate 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 informationUndergraduate 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 informationJACOBI 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 informationSuboptimal 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 informationTime-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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationEntropy-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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More information