Algebraic flux correction and its application to convection-dominated flow. Matthias Möller
|
|
- Darcy Lucas
- 5 years ago
- Views:
Transcription
1 Algebraic flux correction and its application to convection-dominated flow Matthias Möller Institute of Applied Mathematics (LS III) University of Dortmund, Germany 8th Hirschegg Workshop on Conservation Laws 9-15 September 2007, Hirschegg 1
2 Overview 1 Algebraic flux correction Motivation Design principles Flux limiting 2 Nonlinear solution algorithm Choice of preconditioners Construction of Jacobians Numerical examples 3 Recent work Application to compressible flows Grid adaptivity 4 Conclusions 2
3 Motivation Scalar transport equation u t + (vu) = 0 High-order methods: wiggles Standard Galerkin FEM M c du dt + Ku = 0 Low-order methods: smearing Remedy: nonlinear combination of high- and low-order discretizations. 3
4 Algebraic flux correction, Kuz05a High-order scheme M c du dt + Ku = 0 K = {k ij}, k ij > 0, j i Apply artificial diffusion D = {d ij }, d ij = max{k ij, 0, k ji } = d ji Low-order scheme M l du dt + Lu = 0 L = K + D, l ij 0, j i Decompose the residual difference into internodal fluxes r H i r L i = f ij, j i f ij = [ m ij d dt + d ij] (uj u i ) = f ji Nonlinear high-resolution scheme containing limited antidiffusion M l du dt + Lu = f, f i = j i α ij f ij, 0 α ij 1 4
5 Algebraic flux correction, Kuz05a High-order scheme M c du dt + Ku = 0 K = {k ij}, k ij > 0, j i Apply artificial diffusion D = {d ij }, d ij = max{k ij, 0, k ji } = d ji Low-order scheme M l du dt + Lu = 0 L = K + D, l ij 0, j i Decompose the residual difference into internodal fluxes r H i r L i = f ij, j i f ij = [ m ij d dt + d ij] (uj u i ) = f ji Nonlinear high-resolution scheme containing limited antidiffusion M l du dt + Lu = f, f i = j i α ij f ij, 0 α ij 1 4
6 Algebraic flux correction, Kuz05a High-order scheme M c du dt + Ku = 0 K = {k ij}, k ij > 0, j i Apply artificial diffusion D = {d ij }, d ij = max{k ij, 0, k ji } = d ji Low-order scheme M l du dt + Lu = 0 L = K + D, l ij 0, j i Decompose the residual difference into internodal fluxes r H i r L i = f ij, j i f ij = [ m ij d dt + d ij] (uj u i ) = f ji Nonlinear high-resolution scheme containing limited antidiffusion M l du dt + Lu = f, f i = j i α ij f ij, 0 α ij 1 4
7 Algebraic flux correction, Kuz05a High-order scheme M c du dt + Ku = 0 K = {k ij}, k ij > 0, j i Apply artificial diffusion D = {d ij }, d ij = max{k ij, 0, k ji } = d ji Low-order scheme M l du dt + Lu = 0 L = K + D, l ij 0, j i Decompose the residual difference into internodal fluxes r H i r L i = f ij, j i f ij = [ m ij d dt + d ij] (uj u i ) = f ji Nonlinear high-resolution scheme containing limited antidiffusion M l du dt + Lu = f, f i = j i α ij f ij, 0 α ij 1 4
8 Flux limiting Nonlinear high-resolution scheme (lumped-mass version) M l du dt + Lu = f, f i = j i α ij f ij, f ij = d ij (u j u i ) = f ji α ij 0 f = 0 M l du dt + Lu = 0 α ij 1 f = Du M l du dt + Ku = 0 low-order scheme high-order scheme Goal: Determine α ij as close to 1 as possible such that (Lu f ) i = j i q ij (u j u i ), q ij 0 j i 5
9 Multidimensional flux correction 1. Compute the sums of positive/negative antidiffusive fluxes P + i = j i max{0, f ij }, P i = j i min{0, f ij } 2. Define the corresponding upper/lower bounds Q + i = j i q ij max{0, u j u i }, Q i = j i q ij min{0, u j u i } 3. Evaluate the nodal correction factors for positive/negative fluxes R + i = min{1, Q + i /P + i }, R i = min{1, Q i /P i } 4. Apply correction factors to the raw antidiffusive fluxes f i = j i α ij f ij a ij = { min{r + i, R j }, if f ij 0 min{r i, R + j }, otherwise 6
10 Multidimensional flux correction 1. Compute the sums of positive/negative antidiffusive fluxes P + i = j i max{0, f ij }, P i = j i min{0, f ij } 2. Define the corresponding upper/lower bounds Q + i = j i q ij max{0, u j u i }, Q i = j i q ij min{0, u j u i } 3. Evaluate the nodal correction factors for positive/negative fluxes R + i = min{1, Q + i /P + i }, R i = min{1, Q i /P i } 4. Apply correction factors to the raw antidiffusive fluxes f i = j i α ij f ij a ij = { min{r + i, R j }, if f ij 0 min{r i, R + j }, otherwise 6
11 Multidimensional flux correction 1. Compute the sums of positive/negative antidiffusive fluxes P + i = j i max{0, f ij }, P i = j i min{0, f ij } 2. Define the corresponding upper/lower bounds Q + i = j i q ij max{0, u j u i }, Q i = j i q ij min{0, u j u i } 3. Evaluate the nodal correction factors for positive/negative fluxes R + i = min{1, Q + i /P + i }, R i = min{1, Q i /P i } 4. Apply correction factors to the raw antidiffusive fluxes f i = j i α ij f ij a ij = { min{r + i, R j }, if f ij 0 min{r i, R + j }, otherwise 6
12 Multidimensional flux correction 1. Compute the sums of positive/negative antidiffusive fluxes P + i = j i max{0, f ij }, P i = j i min{0, f ij } 2. Define the corresponding upper/lower bounds Q + i = j i q ij max{0, u j u i }, Q i = j i q ij min{0, u j u i } 3. Evaluate the nodal correction factors for positive/negative fluxes R + i = min{1, Q + i /P + i }, R i = min{1, Q i /P i } 4. Apply correction factors to the raw antidiffusive fluxes f i = j i α ij f ij a ij = { min{r + i, R j }, if f ij 0 min{r i, R + j }, otherwise 6
13 Nonlinear solution algorithm Nonlinear algebraic system for 0 < θ 1 N(u) := Lu f M l u n+1 u n t + θn(u n+1 ) + (1 θ)n(u n ) = 0 Fixed-point iteration u (0) := u n u (m+1) = u (m) C 1 F (u (m) ), m = 0, 1,... until F (u (m+1) ) tol or F (u (m+1) ) tol F (u n ). Practical implementation C u (m+1) = F (u (m) ) m = 0, 1,... u (m+1) = u (m) u (m+1) u (0) = u n 7
14 Nonlinear solution algorithm Nonlinear algebraic system for 0 < θ 1 N(u) := Lu f F (u n+1 u n+1 u n ) := M l + θn(u n+1 ) + (1 θ)n(u n ) = 0 t Fixed-point iteration u (0) := u n u (m+1) = u (m) C 1 F (u (m) ), m = 0, 1,... until F (u (m+1) ) tol or F (u (m+1) ) tol F (u n ). Practical implementation C u (m+1) = F (u (m) ) m = 0, 1,... u (m+1) = u (m) u (m+1) u (0) = u n 7
15 Choice of preconditioner Linearized system C u (m+1) = F (u (m) ) Diagonal mass matrix C = M l / t cheap, but only usable for very small time steps Low-order operator C = M l / t + θ L fixed-point defect correction scheme may converge slowly no update needed if the velocity field remains unchanged Jacobian operator C = M l / t + θ N(u) u nonlinear operator N(u) is constructed discretely flux limiter may not be globally differentiable Idea: approximate by divided differences and see what happens 8
16 Jacobian matrix, Moe07a Nodal approximation of Jacobian J = {j ik } D k [N i ] := N i(u + he k ) N i (u he k ) 2h j ik = D k [N i ] + O(h 2 ) D k [N i ] = D k [(Lu f ) i ] = D k [(Ku + Du f ) i ] = D k [ j k ij u j + j i (1 α ij )f ij ] = k ik + j D k [k ij ]u j + j i D k [(1 α ij )f ij ] Averaged coefficient kik := k ik(u + he k ) + k ik (u he k ) 2 9
17 Jacobian matrix, Moe07a Nodal approximation of Jacobian J = {j ik } D k [N i ] := N i(u + he k ) N i (u he k ) 2h j ik = D k [N i ] + O(h 2 ) D k [N i ] = D k [(Lu f ) i ] = D k [(Ku + Du f ) i ] = D k [ j k ij u j + j i (1 α ij )f ij ] = k ik + j D k [k ij ]u j + j i D k [(1 α ij )f ij ] Averaged coefficient kik := k ik(u + he k ) + k ik (u he k ) 2 9
18 Jacobian matrix, Moe07a Nodal approximation of Jacobian J = {j ik } D k [N i ] := N i(u + he k ) N i (u he k ) 2h j ik = D k [N i ] + O(h 2 ) D k [N i ] = D k [(Lu f ) i ] = D k [(Ku + Du f ) i ] = D k [ j k ij u j + j i (1 α ij )f ij ] = k ik + j D k [k ij ]u j + j i D k [(1 α ij )f ij ] Averaged coefficient kik := k ik(u + he k ) + k ik (u he k ) 2 9
19 Sparsity pattern, Moe07a Jacobian coefficients j ik = k ik + j D k[k ij ]u j + j i D k[(1 α ij )f ij ] There exists an edge ij iff the FE basis functions have overlapping supports G := K Structure of the Jacobian Remarks z kk 0, g ii = 1, g ij = 1 ij Z := J z ik 0 ik z jk 0 ik ij Same structure is used for edge-oriented stabilization techniques Sparsity pattern can be generated by multiplication: Z = G 2 10
20 Sparsity pattern, Moe07a Jacobian coefficients j ik = k ik + j D k[k ij ]u j + j i D k[(1 α ij )f ij ] There exists an edge ij iff the FE basis functions have overlapping supports G := K g ii = 1, g ij = 1 ij u he k k Structure of the Jacobian Z := J z kk 0, z ik 0 ik z jk 0 ik ij Remarks Same structure is used for edge-oriented stabilization techniques Sparsity pattern can be generated by multiplication: Z = G 2 10
21 Sparsity pattern, Moe07a Jacobian coefficients j ik = k ik + j D k[k ij ]u j + j i D k[(1 α ij )f ij ] There exists an edge ij iff the FE basis functions have overlapping supports G := K Structure of the Jacobian Remarks z kk 0, g ii = 1, g ij = 1 ij Z := J z ik 0 ik z jk 0 ik ij Same structure is used for edge-oriented stabilization techniques Sparsity pattern can be generated by multiplication: Z = G 2 10
22 Sparsity pattern, Moe07a Jacobian coefficients j ik = k ik + j D k[k ij ]u j + j i D k[(1 α ij )f ij ] There exists an edge ij iff the FE basis functions have overlapping supports G := K Structure of the Jacobian Remarks z kk 0, g ii = 1, g ij = 1 ij Z := J z ik 0 ik z jk 0 ik ij Same structure is used for edge-oriented stabilization techniques Sparsity pattern can be generated by multiplication: Z = G 2 10
23 Sparsity pattern, Moe07a Jacobian coefficients j ik = k ik + j D k[k ij ]u j + j i D k[(1 α ij )f ij ] There exists an edge ij iff the FE basis functions have overlapping supports G := K Structure of the Jacobian Remarks z kk 0, g ii = 1, g ij = 1 ij Z := J z ik 0 ik z jk 0 ik ij Same structure is used for edge-oriented stabilization techniques Sparsity pattern can be generated by multiplication: Z = G 2 10
24 Example: Stationary convection-diffusion Convection-diffusion equation v u d u = 0 v = (cos 10, sin 10 ) FEM-TVD: d = 10 3 Initial guess in Ω = (0, 1) 2 u u (x, y) = Boundary conditions { 1 x if y otherwise Q 1 -elements u(x, 0) = 0, u(1, y) = 0, { u 1 if y 0.5, (x, 1) = 0, u(0, y) = y 0 otherwise. 11
25 Example: τ u + v u d u = 0 d = 10 3, τ = 1.0, BiCGSTAB+ILU, rel lin 10 3 norm of nonlinear residual 10 0 Newton level 6 Newton level Newton level 8 Defcor level 6 Defcor level Defcor level number of nonlinear iterations NVT CPU NN NL Newton s method defect correction Newton s method vs. defect correction total CPU time reduces by factor 2-3 reduction of nonlinear iterations by factor 4-16
26 Example: τ u + v u d u = 0 d = 10 3, τ = 10.0, BiCGSTAB+ILU, rel lin 10 3 norm of nonlinear residual 10 0 Newton level 6 Newton level Newton level 8 Defcor level 6 Defcor level Defcor level number of nonlinear iterations NVT CPU NN NL Newton s method defect correction Newton s method vs. defect correction total CPU time reduces by factor 2-3 reduction of nonlinear iterations by factor 4-16
27 Example: τ u + v u d u = 0 d = 10 4, τ = 1.0, BiCGSTAB+ILU, rel lin 10 3 norm of nonlinear residual 10 0 Newton level 6 Newton level Newton level 8 Defcor level 6 Defcor level Defcor level number of nonlinear iterations NVT CPU NN NL Newton s method defect correction Newton s method vs. defect correction total CPU time reduces by factor 2-3 reduction of nonlinear iterations by factor 4-16
28 Example: τ u + v u d u = 0 d = 10 4, τ = 10.0, BiCGSTAB+ILU, rel lin 10 3 norm of nonlinear residual 10 0 Newton level 6 Newton level Newton level 8 Defcor level 6 Defcor level Defcor level number of nonlinear iterations NVT CPU NN NL Newton s method defect correction Newton s method vs. defect correction total CPU time reduces by factor 2-3 reduction of nonlinear iterations by factor
29 Example: τ u + v u d u = 0 d = 10 4, τ = 10.0, BiCGSTAB+ILU, rel lin 10 3 norm of nonlinear residual 10 0 Newton level 6 Newton level Newton level 8 Defcor level 6 Defcor level Defcor level number of nonlinear iterations NVT CPU NN NL Newton s method defect correction Newton s method vs. defect correction total CPU time reduces by factor 2-3 reduction of nonlinear iterations by factor
30 Example: Burgers equation in space-time Inviscid Burgers equation ( ) u 2 x + t u = 0 2 Discrete upwind: 32,768 P 1 -elements Space-time domain Ω = (0, 1) (0, 0.5) Boundary conditions 1 if 0 x < 0.4 t = if 0.4 x 8 t = 0 u(x, t) = 0 if 0.8 < x 1 t = 0 or x = 0 0 t 0.5 u u h 1 = e 2 u u h 2 = e 2 13
31 Example: Burgers equation in space-time Inviscid Burgers equation ( ) u 2 x + t u = 0 2 FEM-TVD: 32,768 P 1 -elements Space-time domain Ω = (0, 1) (0, 0.5) Boundary conditions 1 if 0 x < 0.4 t = if 0.4 x 8 t = 0 u(x, t) = 0 if 0.8 < x 1 t = 0 or x = 0 0 t 0.5 u u h 1 = e 3 u u h 2 = e 2 13
32 Example: τ u + x ( 1 2 u2 ) + t u = 0 Discrete upwind: 1 nit./ τ, BiCGSTAB+ILU, rel lin 10 3 τ = 0.1 Newton defect correction error NVT CPU NN NL CPU NN NL u u h 1 u u h e e e e e e-2 τ = 1.0 Newton defect correction error NVT CPU NN NL CPU NN NL u u h 1 u u h e e e e e e-2 Discrete Jacobian operator yields robust Newton iteration 14
33 Example: τ u + x ( 1 2 u2 ) + t u = 0 FEM-TVD: up to 10 nit./ τ n, τ n [0.01, 1.0], rel nonl 10 3 BiCGSTAB+ILU, rel lin 10 3 Newton defect correction CPU NSTEP NN NL CPU NSTEP NN NL error u u h e e e-3 u u h e e e-2 Newton s method vs. defect correction total CPU time reduces by factor 3 reduction of outer iterations by factor
34 Example: Swirling flow, Leveque Swirling flow in Ω = (0, 1) 2 u t + (vu) = 0 Time-dependent velocity t (0, T ) v x = sin 2 (πx) sin(2πy)g(t) v y = sin 2 (πy) sin(2πx)g(t) g(t) = cos(πt/t ), T = 1.5 Time step control: PID t [10 3, 10 1 ], un+1 u n u n+1 0.5% 16
35 Example: Swirling flow, contd. Newton defect correction NVT NSTEP CPU NN NL CPU NN NL (70) (77) (66) nodes nodes nodes adaptive time step error reduction 17 CPU improved by factor 4.5; reduction of iterations by factor 7
36 Example: Swirling flow, contd. Newton defect correction NVT NSTEP CPU NN NL CPU NN NL (70) (77) (66) nodes nodes nodes norm of nonlinear residual adaptive time step number of nonlinear iterations 17 CPU improved by factor 4.5; reduction of iterations by factor 7
37 Outlook: Compressible flows, Kuz05b Divergence Form U t + 3 d=1 F d x d = 0 A d = F d U Quasi-linear form U t + 3 U d=1 Ad x d = 0 FEM discretization [ M du ] C dt = i j i a ij(u j u i ) a ij = r ij Λ ij r 1 ij Transformation to local characteristic variables makes it possible to perform upwinding and algebraic flux correction as in the scalar case. M = 3, α = 0, density Mach number 18
38 Grid adaption, Moe triangles triangles triangles triangles Density distribution on the final mesh cutline at y = 0.6 grid1 grid2 grid3 grid
39 Time-dependent grid adaptivity 20
40 Time-dependent grid adaptivity 20
41 Time-dependent grid adaptivity 20
42 Time-dependent grid adaptivity 20
43 Time-dependent grid adaptivity 20
44 Time-dependent grid adaptivity 20
45 Time-dependent grid adaptivity 20
46 Conclusions and Outlook nonlinear high-resolution schemes can be constructed by adding artificial diffusion + limited antidiffusion implicit flux correction schemes can be combined with algebraic Newton methods to speed up convergence Jacobian matrix is approximated by divided differences and assembled efficiently edge-by-edge in a black-box fashion the sparsity pattern of the Jacobian needs to be extended which can be accomplished by standard matrix multiplication Further research: compressible Euler/Navier-Stokes equations time-dependent grid adaptation + goal oriented error estimation 21
47 References Kuz05a D. Kuzmin, M.M. Algebraic flux correction I. Scalar conservation laws. In: Flux-Corrected Transport, Principles, Algorithms, and Applications, D. Kuzmin, R. Löhner, S. Turek (eds.). Springer: Germany, 2005; Kuz05b D. Kuzmin, M.M. Algebraic flux correction II. Compressible Euler equations. In: Flux-Corrected Transport, Principles, Algorithms, and Applications, D. Kuzmin, R. Löhner, S. Turek (eds.). Springer: Germany, 2005; Moe06 M.M., D. Kuzmin. Adaptive mesh refinement for high-resolution finite element schemes. IJNMF 52, 2006; Moe07a M.M. Efficient solution techniques for implicit finite element schemes with flux limiters. IJNMF (in press). Moe07b M.M, D. Kuzmin, D. Kourounis. Implicit FEM-FCT algorithms and discrete Newton methods for transient problems. Tech. Rep. 340, 2007, University of Dortmund. 22
48 PID time step control Normalize relative changes e n = un u n 1 tol u n Reject time step if e n > 1 t n := β t n, 0 < β < 1 Compute new time step t n+1 = ( en 1 e n e n Parameters (Valli, Carey, Coutinho) ) kp ( ) ki ( 1 e 2 ) kd n 1 t n e n e n 2 k P = 0.075, k I = 0.175, k D = 0.01 Back 23
49 Semi-implicit FCT limiter, Part I 1. Initialization P ± i Q ± i R ± i 0 2.Compute positivity-preserving intermediate low-order solution ũ = u n + (1 θ) tm 1 l Lu n 3. Evaluate raw antidiffusive flux f n ij P ± i := P ± i + max min {0, f n = td n ij (un i ij }, P ± j := P ± j + max min 4. Update admissible increments for both nodes Q ± i uj n ) and update {0, f n ij } := max min {Q± i, ũ j ũ i }, Q ± j := max min {Q± j, ũ i ũ j } 5. Limite the raw antidiffusive fluxes f n ij R ± i := m i Q ± i /P ± i f ij := { min{r + i, R j }fij n, min{r i, R + j }fij n, if F ij n > 0, otherwise. 24
50 Semi-implicit FCT limiter, Part II 1. Evaluate the target flux f ij = [m ij + θ td (m) ij ](u (m) i u (m) j ) [m ij (1 θ) tdij n ](ui n uj n ) 2.Constrain each flux by means of f ij { f ij = min{f ij, max{0, f ij }}, if f ij > 0, max{f ij, min{0, f ij }}, otherwise. 3. Insert the limited antidiffusive flux into the right-hand side b(m) i := b (m) i + f ij, b(m) j := b (m) j f ij Back 25
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 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 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 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 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 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 finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More informationA flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds
A flux-corrected RBF-FD method for convection dominated problems in domains and on manifolds Andriy Sokolov Dmitri Kuzmin, Oleg Davydov and Stefan Turek Institut für Angewandte Mathematik (LS3) TU Dortmund
More informationAlgebraic Flux Correction II
Algebraic Flux Correction II Compressible Flow Problems Dmitri Kuzmin, Matthias Möller, and Marcel Gurris Abstract Flux limiting for hyperbolic systems requires a careful generalization of the design principles
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
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 informationEfficient FEM-multigrid solver for granular material
Efficient FEM-multigrid solver for granular material S. Mandal, A. Ouazzi, S. Turek Chair for Applied Mathematics and Numerics (LSIII), TU Dortmund STW user committee meeting Enschede, 25th September,
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
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 informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNUMERICAL SIMULATION OF TURBULENT BUBBLY FLOWS
3rd International Symposium on Two-Phase Flow Modelling and Experimentation Pisa, 22-24 September 2004 NUMERICAL SIMULATION OF TURBULENT BUBBLY FLOWS Dmitri Kuzmin and Stefan Turek Institute of Applied
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 informationNUMERICAL SIMULATION OF TURBULENT BUBBLY FLOWS
3rd International Symposium on Two-Phase Flow Modelling and Experimentation Pisa, 22-24 September 2004 NUMERICAL SIMULATION OF TURBULENT BUBBLY FLOWS Dmitri Kuzmin and Stefan Turek Institute of Applied
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 informationAlgebraic Flux Correction III
Algebraic Flux Correction III Incompressible Flow Problems Stefan Turek and Dmitri Kuzmin Abstract This chapter illustrates the use of algebraic flux correction in the context of finite element methods
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 informationAcceleration of a Domain Decomposition Method for Advection-Diffusion Problems
Acceleration of a Domain Decomposition Method for Advection-Diffusion Problems Gert Lube 1, Tobias Knopp 2, and Gerd Rapin 2 1 University of Göttingen, Institute of Numerical and Applied Mathematics (http://www.num.math.uni-goettingen.de/lube/)
More informationAlgebraic Flux Correction III. Incompressible Flow Problems
Algebraic Flux Correction III. Incompressible Flow Problems Stefan Turek 1 and Dmitri Kuzmin 2 1 Institute of Applied Mathematics (LS III), University of Dortmund Vogelpothsweg 87, D-44227, Dortmund, Germany
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
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 informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationNumerical Simulation of Powder Flow
Numerical Simulation of Powder Flow Stefan Turek, Abderrahim Ouazzi Institut für Angewandte Mathematik, Univ. Dortmund http://www.mathematik.uni-dortmund.de/ls3 http://www.featflow.de Models for granular
More informationNEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS
NEWTON-GMRES PRECONDITIONING FOR DISCONTINUOUS GALERKIN DISCRETIZATIONS OF THE NAVIER-STOKES EQUATIONS P.-O. PERSSON AND J. PERAIRE Abstract. We study preconditioners for the iterative solution of the
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationPALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations
2013 SIAM Conference On Computational Science and Engineering Boston, 27 th February 2013 PALADINS: Scalable Time-Adaptive Algebraic Splitting and Preconditioners for the Navier-Stokes Equations U. Villa,
More informationFine-grained Parallel Incomplete LU Factorization
Fine-grained Parallel Incomplete LU Factorization Edmond Chow School of Computational Science and Engineering Georgia Institute of Technology Sparse Days Meeting at CERFACS June 5-6, 2014 Contribution
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationConvergence and accuracy of sea ice dynamics solvers. Martin Losch (AWI) Annika Fuchs (AWI), Jean-François Lemieux (EC)
Convergence and accuracy of sea ice dynamics solvers Martin Losch (AWI) Annika Fuchs (AWI), Jean-François Lemieux (EC) Sea ice thickness with MITgcm (JPL) Review of 2D momentum equations m Du Dt = mf k
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids 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 informationA Variational Multiscale Stabilized Finite Element. Method for the Solution of the Euler Equations of
A Variational Multiscale Stabilized Finite Element Method for the Solution of the Euler Equations of Nonhydrostatic Stratified Flows M. Moragues S. Marras M. Vázquez G. Houzeaux Barcelona Supercomputing
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationA monolithic FEM solver for fluid structure
A monolithic FEM solver for fluid structure interaction p. 1/1 A monolithic FEM solver for fluid structure interaction Stefan Turek, Jaroslav Hron jaroslav.hron@mathematik.uni-dortmund.de Department of
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
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 informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
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 informationA Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization. Laslo Tibor Diosady
A Linear Multigrid Preconditioner for the solution of the Navier-Stokes Equations using a Discontinuous Galerkin Discretization by Laslo Tibor Diosady B.A.Sc., University of Toronto (2005) Submitted to
More informationA unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations
Noname manuscript No. (will be inserted by the editor) A unified analysis of Algebraic Flux Correction schemes for convection-diffusion equations Gabriel R. Barrenechea Volker John Petr Knobloch Richard
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
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 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 informationC. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2
s: application to chemical vapor deposition C. Vuik 1 S. van Veldhuizen 1 C.R. Kleijn 2 1 Delft University of Technology J.M. Burgerscentrum 2 Delft University of Technology Department of Multi Scale Physics
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSolving linear systems (6 lectures)
Chapter 2 Solving linear systems (6 lectures) 2.1 Solving linear systems: LU factorization (1 lectures) Reference: [Trefethen, Bau III] Lecture 20, 21 How do you solve Ax = b? (2.1.1) In numerical linear
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
More informationRECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES
RECENT DEVELOPMENTS IN VERY HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR INVISCID AND VISCOUS PROBLEMS. R. Abgrall Team Bacchus INRIA Bordeaux Sud Ouest and UniversitédeBordeauxTalence,France Roscoff,
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
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 informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
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 informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu Univ. of Texas at Austin January 12, 2007 1 Introduction 2 Weighted Residuals 3 Poisson Equation 4 Other Examples 5 Essential
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationNonlinear examples. June 21, Solution by fixed-point iteration 3. 5 Exercises 9
Nonlinear examples June 21, 2012 Contents 1 Bratu s equation 1 1.1 Manufacturing a problem to match a solution.................. 2 1.1.1 A manufactured Bratu s problem in 1D.................. 2 1.2 Weak
More informationA DISCRETE PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW WITH CORIOLIS FORCE
A DISCRETE PROJECTION METHOD FOR INCOMPRESSIBLE VISCOUS FLOW WITH CORIOLIS FORCE ANDRIY SOKOLOV, MAXIM A. OLSHANSKII, AND STEFAN TUREK Abstract. The paper presents a new discrete projection method (DPM)
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More information