Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature. Masayuki Yano
|
|
- Marcia Young
- 5 years ago
- Views:
Transcription
1 Model reduction of parametrized aerodynamic flows: stability, error control, and empirical quadrature Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony Patera; NSERC WCCM 2018 New York, United States 26 July 2018
2 Motivation: parametrized aerodynamic flows Goal: rapid and reliable solution of parametrized conservation laws in many-query scenarios. Conservation laws: Euler and Reynolds-averaged Navier-Stokes (RANS) equations Scenarios: 1. parameter sweep (flight-envelope/design study) 2. uncertainty quantification 1
3 Parametrized steady conservation laws µpde: given µ D R P, find u V s.t. ( F (u; µ) + K(u; µ) u) = S(u, u; µ) }{{}}{{}}{{} advection flux diffusion flux source B(u, n K(u; µ) u; µ) = 0 on Ω }{{} boundary operator in Ω µpde (weak): given µ D, find u V s.t. O( ) ops r(u, v; µ) v F (u; µ) + dx = 0 v V. }{{} Ω residual form 2
4 Finite element (FE) method FE space: N -dim space e.g., piecewise polynomial V h V FE integral: K h -pt quadrature Q h {(x k, ρ k )} K h k=1 s.t. e.g., Gauss quadrature [ ]dx Ω (x k,ρ k ) Q h ρ k [ ] xk FE: given µ D, find u h V h s.t. O(N, K h ) ops r h (u h, v; µ) ρ }{{} k [ v F (u h ; µ) + ] xk = 0 v V h. FE residual form (x k,ρ k ) Q h 3
5 Reduced basis (RB) method RB space: N N -dim empirical space V N V h V for parametric manifold {u(µ) µ D}. RB (FE integ.): given µ D, find u N V N s.t. r h (u N, v; µ) }{{} FE residual form O(N, K h ) ops (x k,ρ k ) Q h ρ k [ v N F (u N ; µ) + ] xk = 0 v N V N. 4
6 RB method with empirical quadrature [Ryckelynck; An et al; Farhat et al;... ] EQP: K ν K h -pt quadrature for empirical (parametric) functions (x ν k, ρ ν k) Q ν s.t. [( )(µ)]dx ρ ν k[( )(µ)] x ν k Ω (x ν k,ρν k ) Qν RB-EQP: given µ D, find u ν N V N s.t. O(N, K ν ) ops r ν (u ν N, v; µ) ρ ν }{{} k[ v N F (u ν N; µ)+ ] x ν k = 0 v N V N. EQP residual form (x ν k,ρν k ) Qν 5
7 RB-EQP example: NACA0012 Euler (µ = (α, M )) RB space: V N span{φ i } N i=1 φ 1 φ 2 φ 3 φ 4 EQP: with element-wise grouping N = 4, K ν = 26 N = 10, K ν = 59 6
8 Objectives Goal: stable, accurate, and systematic hyperreduction 7
9 Objectives Goal: stable, accurate, and systematic hyperreduction 1. ensure stability of hyperreduced system energy-stability for linear hyperbolic systems symmetry and positivity for linear diffusion systems r ν (u ν N, v) = 0 v V N can be solved. 7
10 Objectives Goal: stable, accurate, and systematic hyperreduction 1. ensure stability of hyperreduced system energy-stability for linear hyperbolic systems symmetry and positivity for linear diffusion systems r ν (u ν N, v) = 0 v V N can be solved. 2. quantitatively control solution error due to hyperreduction u N u ν N V δ. 7
11 Objectives Goal: stable, accurate, and systematic hyperreduction 1. ensure stability of hyperreduced system energy-stability for linear hyperbolic systems symmetry and positivity for linear diffusion systems r ν (u ν N, v) = 0 v V N can be solved. 2. quantitatively control solution error due to hyperreduction u N u ν N V δ. 3. train hyperreduced system (i.e., identify RB and EQP) in systematic manner r ν (, ) is found systematically. 7
12 Formulation: stable decomposition Overview Hyperbolic systems Diffusion systems (digest) Summary
13 Formulation: stable decomposition Overview Hyperbolic systems Diffusion systems (digest) Summary
14 Discontinuous Galerkin (DG) method DG: given a discontinuous FE space [Reed & Hill; Cockburn & Shu;... ] V h = {v L 2 (Ω) v P p (κ), κ T h }, find u h V h such that, v V h, r h (u h, v) = v F (u h ) dx + v ˆF (u h ) ds = 0. κ T κ h }{{} σ Σ σ h }{{} element integral facet integral Features: stability for conservation laws unstructured meshes hp flexibility u h σ Σ h κ T h 8
15 DG reduced basis (RB) method DG-RB (FE integ.): given a RB space for {u h (µ) µ D}, V N V h, find u N V N such that, v V N, r h (u N, v) = v F (u N ) dx + v ˆF (u N ) ds = 0. κ T κ h }{{} σ Σ σ h }{{} element integral facet integral Hyperreduction: find a hyperreduced form r ν (u N, v N ) r h (u N, v N ) with reduced number of element and facet integrals. 9
16 Hyperreduction of DG residual Step 1: find an element-wise decomposition stability r h (w, v) = η }{{} κ (w, v) }{{} κ T DG residual h element-wise residual such that for any (sparse) weights {ρ ν κ 0} κ Th r ν (w, v) ρ ν }{{} κη κ (w, v) κ T DG-EQP residual h is stable. Step 2. find EQP weights {ρ ν κ} that provides 1. sparsity: nnz{ρ ν κ} T h speed 2. accuracy: u N u ν N V δ accuracy 10
17 Formulation: stable decomposition Overview Hyperbolic systems Diffusion systems (digest) Summary
18 Linear hyperbolic system PDE: for Ω R d and I (0, T ], u t Energy balance: u 2 t 2dx Ω }{{} + Au = 0 in Ω I, (n A) u = (n A) u b on Ω I, u(t = 0) = u 0 in Ω {0}, change in energy = 2 u(n A) u b ds Ω }{{} net energy entering Ω F (u) = Au, A sym. 11
19 DG discretization Steady DG residual: r h (w, v) v F (w)dx + [v] + ˆF (w +, w ; n + )ds Ω Σ }{{}} i {{} volume term interior numerical flux + v + ˆF b (w + ; n + )ds, u Σ h } b {{} boundary numerical flux Energy balance: u 2 t 2dx + Ω }{{} change in energy Σ i [u] + n A [u] + ds }{{} facet jump dissipation: 0 = 2 σ Σ h κ T h u(n A) u b ds Ω }{{} net energy entering Ω 12
20 DG discretization Steady DG residual: r h (w, v) v F (w)dx + [v] + ˆF (w +, w ; n + )ds Ω Σ }{{}} i {{} volume term interior numerical flux + v + ˆF b (w + ; n + )ds, u Σ h } b {{} boundary numerical flux Energy balance: u 2 t 2dx Ω }{{} change in energy σ Σ h κ T h 2 u(n A) u b ds Ω }{{} net energy entering Ω Key: numerical flux provides energy stability. 12
21 Element-wise decomposition Naive decomp.: apply elemental mask to discrete DG residual η κ (w, v) r h (w, v κ ) v F (w)dx + v + ˆF b (w + ; n + )ds κ κ }{{}} b {{} volume term boundary numerical flux + v + ˆF (w +, w ; n + )ds. κ } i {{} one-sided interior numerical flux Element-wise energy: in general, η κ (v, v) 0 no energy stability. 13
22 Energy-stable elemental decomposition Energy-stable decomp.: redistribute interior facet residuals η κ (w, v) v F (w)dx + v + ˆF b (w + ; n + )ds κ κ }{{}} b {{} volume term boundary numerical flux + 1 [vn F (w; µ)] ± 4 ds + 1 [v] + κ 2 ˆF (w +, w ; n + ; µ)ds } i κ {{}} i {{} energy redistribution term split interior numerical flux Element-wise energy: for any v V h, η κ (v, v) 0 (modulo BC) energy stability. 14
23 DG-EQP discretization Steady DG-EQP residual: for weights {ρ ν κ 0} κ Th, r ν (w, v) κ T h ρ ν κη κ (w, v). Energy balance: r ν (u ν N, uν N ) 0 (modulo BC) and hence ( ) u ν t N 2 2dx 2 ρ ν κ u ν + N (n A) u b ds. κ T h κ T h ρ ν κ κ }{{} EQP approx of change in energy σ κ Σ b h } {{ } EQP approx of net energy entering Ω σ 15
24 Formulation: stable decomposition Overview Hyperbolic systems Diffusion systems (digest) Summary
25 Linear diffusion system PDE: for Ω R d Σ i h K u = 0 in Ω, } {{ } interior flux consistency term [c.f., Farhat et al, Carlberg et al.] u = u b on Ω. DG residual r h (w, v) v K wdx + ds Ω Σ }{{} b h }{{} volume term BC term { v K} w + v {K w}ds v {Kl σ ( w )}ds.. SPD: the bilinear form r h (, ) is 1. symmetric: r h (w, v) = r h (v, w) Σ i h } {{ } lifting term 2. positive: r h (v, v) 0. 16
26 Element-wise decomposition of DG residual Naive decomp.: apply elemental mask to discrete residual η κ (w, v) r h (w, v κ ). In general, η κ (, ) is neither symmetric nor non-negative. Structure-preserving decomp.: redistribute interior facet residuals η κ (w, v) = v K wdx + ds κ κ }{{}} b {{} volume term BC term 1 v K w + v K wds 1 v Kl σ ( w )ds 2 κ 2 } i κ {{}} i {{} redistributed interior flux consistency term redistributed lifting term For any w, v V h, η κ (, ) is symmetric and non-negative. 17
27 Formulation: stable decomposition Overview Hyperbolic systems Diffusion systems (digest) Summary
28 Summary: stable decomposition Stable decomp.: for any (sparse) weights {ρ ν κ 0} κ Th the hyperreduced system is r ν (w, v) ρ ν }{{} κη κ (w, v), κ T DG-EQP residual h 1. energy stable for linear hyperbolic systems 2. symmetric and positive for linear diffusion systems Key: naive elemental-mask decomp. stable decomp. 18
29 Formulation: hyperreduction error control
30 RB approximation: discrete form Reduced basis: Z N [φ 1,..., φ N ] RB residual: r N ( w }{{} state RB-EQP residual: ; µ }{{} i = param.) κ T h energy-stable elem. res. {}}{ η κ (Z N w, φ i ; µ) }{{} DG res. form: r h (Z N w,φ i ;µ) r ν N(w; µ) κ T h ρ ν κη N,κ (w; µ), i = 1,..., N = κ T h η N,κ (w; µ) i where EQP weight {ρ ν κ} κ Th achieves 1. sparsity: nnz{ρ ν κ} κ Th T h 2. accuracy: u N (µ) u ν N (µ) V δ, µ D note: control accuracy of solution not quadrature 19
31 RB-EQP: linear program (LP) RB-EQP: set δ δ/ N R >0. Find basic feasible solution ρ = arg min ρ subject to κ T h ρ κ 1. non-negativity: ρ κ 0, κ T h 2. constant accuracy: Ω κ T h ρ κ κ δ 3. manifold accuracy: define g N,κ (µ) J N (u train N (µ); µ) 1 η N,κ (u train N (µ); µ) ; }{{} Jacobian-preconditioned element residual for Ξ train J {µ train j } J j=1 and UJ train {u train N g N (µ) κ T h ρ κ g N,κ (µ) δ (µ)} µ Ξ train, J µ Ξ train J 20
32 RB-EQP: linear program (LP) matrix-vector form RB-EQP: set δ R 0. Find basic feasible solution subject to 1. non-negativity: ρ 0 ρ = arg min ρ R K h 1 T ρ 2. constant accuracy: { κ } T κ T h ρ < > Ω ± δ 3. manifold accuracy: define g N,κ (µ) J N (u train N (µ); µ) 1 η N,κ (u train N (µ); µ) ; }{{} Jacobian-preconditioned element residual g N,1 (µ train 1 )... g N,Kh (µ train 1 ) ρ < > g N,1 (µ train J )... g N,Kh (µ train J ) ρ Kh }{{} (N J) K h g N (µ train 1 ). g N (µ train J ) ± δ 21
33 RB-EQP: properties 1. simple: solved by (dual) simplex method 2. sparse: nnz{ρ ν κ} κ Th T h intuition: l 1 minimization 3. error control: under mild assumptions, if J N (µ) 1 r N (µ) J N (µ) 1 r ν N(µ)) }{{} δ, manifold accuracy constraints I J N (µ) 1 J ν N(µ) max ɛ, then u N (µ) u ν N(µ) V δ. Proof: application of Brezzi-Rappaz-Raviart (BRR). 22
34 RB-EQP: properties 1. simple: solved by (dual) simplex method 2. sparse: nnz{ρ ν κ} κ Th T h intuition: l 1 minimization 3. error control: under mild assumptions, if J N (µ) 1 r N (µ) J N (µ) 1 r ν N(µ)) }{{} δ, manifold accuracy constraints I J N (µ) 1 J ν N(µ) max ɛ, then u N (µ) u ν N(µ) V δ. Proof: application of Brezzi-Rappaz-Raviart (BRR). 22
35 Formulation: automated training
36 Simultaneous FE, RB, and EQP greedy training Input: training set Ξ train J Output: RB V N ; EQP weights {ρ ν κ} κ Th D; EQP tolerance δ; FE tolerance ɛ FE For N = 1,..., N max 1. Find µ N that maximizes the truth dual-norm µ N = arg sup r h (u ν N 1(µ); ; µ) (V h ) µ Ξ train J 2. Solve truth problem: u h (µ N ); enrich FE space as necessary. 3. Update reduced basis: Z N = {Z N 1, u h (µ N )} 4. Update EQP rule repeat N smooth times bootstrap Set U train J {u train N (µ) uν N (µ)} µ Ξ train J Solve LP-EQP for {Ξ train J, UJ train } 23
37 Computational cost per iteration 1. truth residual sampling and J Ξ train J 4. EQP update by bootstrapping J RB-EQP solutions O(JN ) J truth residual evaluations O(JN ) 2. truth snapshot calculation 1 truth solution evaluation O(N ) Note: solution evaluation residual evaluation 24
38 DG-FEM anisotropic mesh adaptation Employ Solve Estimate Aniso Mark Refine. Solve: high-order DG-FEM Estimate: dual-weighted residual (DWR) method Anisotropy Mark: error-to-dof optimization via local solves Refine: anisotropic hanging-node refinement 25
39 Related work
40 Most relevant related work (1/2) Nonlinear MOR: interpolate-then-integrate Gappy POD, MPE, EIM, BPIM, GNAT,... [Everson & Sirovich; Astrid et al; Barrault et al; Nguyen et al; Carlberg et al;... ] Nonlinear MOR: direct integration Hyperreduction [Ryckelynck] Optimal cubature [An et al] Energy-conservative sampling and weighting (ECSW) [Farhat et al] Empirical cubature method [Hernández] EQP for continuous Galerkin Empirical quadrature (for polynomials on arbitrary domains): LP framework [Ryu & Boyd; DeVore et al] 26
41 Most relevant related work (2/2) Structure-preserving/stability-aware MOR: ECSW for solids [Farhat et al] Matrix gappy POD for solids [Carlberg et al] Stable Galerkin ROM for fluids [Barone et al; Kalashinova et al] MOR for parametrized aerodynamics: Euler [LeGresley and Alonso] Euler and RANS [Washabaugh et al] Euler [Zimmermann and Görtz] DG-RB-EQP: hyperreduction with (i) energy stability, (ii) direct error control, (iii) simple LP, and (iv) adaptive snapshots. 27
42 NACA0012 Euler ONERA M6 Euler RAE2822 RANS MDA high-lift RANS Example
43 NACA0012 Euler ONERA M6 Euler RAE2822 RANS MDA high-lift RANS Example
44 NACA0012 Euler Equation: Euler (entropy variables) Parameters.: 1. angle of attack: α [0, 3 ] 2. Mach number: M [0.3, 0.5] (M loc,max > 0.8) 28
45 DG-RB-EQP setup: DG setup: goal-oriented adaptive P 2 DG output (c l ) tolerance: δ DG = 10 4 RB-EQP setup: Ξ train J : 5 5 uniform grid. L 2 (Ω) tolerance: δ EQP = 10 4 Note: DG controls output error but RB-EQP controls L 2 error; full goal-oriented error control is ongoing work. 29
46 NACA0012 Euler: mesh (µ = µ centroid ) Initial: P 2, dof = 13440, c l (5.2% error) Adapted: P 2, dof = 47112, c l counts (0.1% error) 30
47 c l error NACA0012 Euler: spatial error convergence 10-1 best-practice (p=1) ani-h adapt (p=2) % error % error dof For 1% error level, 2nd-order method on best-practice mesh: 400k dof p = 2 anisotropic-h adaptation: 16k dof ( 25 reduction) 31
48 NACA0012 Euler: DG-RB-EQP Reduced basis 1 4 (x-momentum) φ 1 φ 2 φ 3 φ 4 Integration weights: N = 4, K = 26 N = 10, K = 59 32
49 NACA0012 Euler: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) { RB vs RB-EQP }} { N K u N u ν N L 2 (Ω)
50 NACA0012 Euler: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) N K
51 NACA0012 Euler: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) FE vs RB-EQP {}}{ J h J ν h N K
52 NACA0012 Euler: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) FE vs RB-EQP {}}{ J h J ν h N K Online timing (1 thread): N = 10: t rb = 0.23 sec ( 160 speedup wrt adaptive DG) reduced residual, linear solve, & Newton iterations 33
53 NACA0012 Euler ONERA M6 Euler RAE2822 RANS MDA high-lift RANS Example
54 ONERA M6 Euler Equation: Euler (entropy variables) Parameters.: 1. angle of attack: α [0, 3 ] 2. Mach number: M [0.3, 0.5] 34
55 ONERA M6 Euler: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) FE vs RB-EQP {}}{ J h J ν h N K Online timing (24 threads ): N = 10: t rb = 0.81 sec ( 130 speedup) DG-RB-EQP implementation is not correctly load balanced. 35
56 NACA0012 Euler ONERA M6 Euler RAE2822 RANS MDA high-lift RANS Example
57 RAE RANS UQ: problem setup Goal: combined model and discretization error quantification Equation: RANS equations with Spalart-Allmaras turbulence model Parameters: [Spalart & Allmaras; Schaefer et al] 1. turbulent Prandtl number: σ [0.60, 1.00] 2. Kármán constant: κ [0.38, 0.42] 3. 2nd wall destruction constant: c w2 [0.10, 0.35] 4. 3rd wall destruction constant: c w3 [1.75, 2.50] Flow cond.: M = 0.3, Re c = , α =
58 RANS-SA equations RANS mean-flow equations: (ρu) = 0 (ρu u + pi) = (2(µ + µ t )τ) (ρuh) = (c p ( µ + µt ) T + τu) Pr Pr for τ = 1 2 ( u + ut ) 1 3 tr( u)i, H =.... SA model: µ t = ρ νf v1 for f v1 = χ3 χ 3 +c 3 v1 (ρu ν) = ρ(c b1 (1 f t2 ) S ν c w1 f w ( ν d )2 ), χ = ν ν, where 1 (ρ(ν + ν) ν) c b2 σ σ ρ( ν)2 + 1 (ν + ν) ρ ν σ and S = S + ν κ 2 d 2 (1 χ 1+χf v1 ) f w = g( 1+c w3 6 g 6 +c w3 6 ) 1/6, g = r + c w2 (r 6 r), r = ν Sκ 2 d 2 c w1 = c b1 κ c b2 σ 37
59 RAE RANS: mesh (µ = µ centroid ) Initial: P 2, dof = 15180, c d counts (27% error) Adapted: P 2, dof = 41250, c d 90.1 counts (0.3% error) 38
60 c d error RAE RANS: spatial error convergence best-practice (p=1) ani-h adapt (p=2) % error % error dof For 0.3% error level, 2nd-order method on best-practice mesh: 1M dof p = 2 anisotropic-h adaptation: 40k dof ( 25 reduction) 39
61 RAE RANS: DG-RB-EQP Reduced basis 1 4 (x-momentum) Integration weights: N = 4, K = 24 N = 7, K = 55 40
62 RAE RANS: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) FE vs RB-EQP {}}{ J h J ν h / J h N K Online timing (1 thread): N = 7: t rb = 0.77 sec ( 25 speedup wrt adaptive DG) 41
63 RAE RANS UQ: mean c d estimation Drag coefficient distribution (1000 MC solves) c d (counts) E[c d ] 90.4 counts online timing: 13 min (for 1000 solves) error estimate: 0.3% FE }{{} spatial + 0.4% RB-EQP }{{} model reduction % MC }{{} statistical 42
64 NACA0012 Euler ONERA M6 Euler RAE2822 RANS MDA high-lift RANS Example
65 MDA high-lift RANS Equation: RANS equations with Spalart-Allmaras turbulence model Parameters: [Spalart & Allmaras; Schaefer et al] 1. Kármán constant: κ [0.38, 0.42] 2. freestream turbulence intensity χ ν/ν [3, 30] Condition: M = 0.2, Re c = , α = 16 (M loc,max > 0.75) 43
66 MDA high-lift RANS: mesh (µ = µ centroid ) Initial: P 2, dof = , c d counts (47% error) Adapted: P 2, dof = , c d counts (0.1% error) 44
67 c d error MDA high-lift RANS: spatial error convergence 10 0 best-practice (p=1) ani-h adapt (p=2) 1% error 0.1% error dof For 1% error level, 2nd-order method on best-practice mesh: 25M dof p = 2 anisotropic-h adaptation: 250k dof ( 100 reduction) 45
68 MDA high-lift RANS: DG-RB-EQP Reduced basis 1 4 (x-momentum) Integration weights: L 2 (Ω) control output control N = , K =
69 MDA high-lift RANS: DG-RB-EQP convergence Assessment: max error over Ξ test Ξ train of Ξ test = 30 (δ = 10 4 ) RB vs RB-EQP {}}{ u N u ν N L 2 (Ω) FE vs RB-EQP {}}{ u h u ν N L 2 (Ω) FE vs RB-EQP {}}{ J h J ν h / J h N K Online timing (24 threads ): N = 10: t rb = 1.8 sec ( 70 speedup wrt adaptive DG) DG-RB-EQP implementation is not correctly load balanced. UQ: c d = counts 47
70 Summary
71 Summary A model reduction framework for conservation laws based on adaptive DG + RB + empirical quadrature (EQP) DG-RB-EQP provides energy stability (linear hyperbolic systems) and symmetry & positivity (linear diffusive systems) solution error control systematic offline training RB, EQP, and FE errors Ongoing work online & goal-oriented error estimates entropy stability treatment of parameter-dependent singularities 48
Model reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure. Masayuki Yano
Model reduction of parametrized aerodynamic flows: discontinuous Galerkin RB empirical quadrature procedure Masayuki Yano University of Toronto Institute for Aerospace Studies Acknowledgments: Anthony
More informationAn Optimization-based Framework for Controlling Discretization Error through Anisotropic h-adaptation
An Optimization-based Framework for Controlling Discretization Error through Anisotropic h-adaptation Masayuki Yano and David Darmofal Aerospace Computational Design Laboratory Massachusetts Institute
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
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 informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More informationSome recent developments on ROMs in computational fluid dynamics
Some recent developments on ROMs in computational fluid dynamics F. Ballarin 1, S. Ali 1, E. Delgado 2, D. Torlo 1,3, T. Chacón 2, M. Gómez 2, G. Rozza 1 1 mathlab, Mathematics Area, SISSA, Trieste, Italy
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 informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Methods for Nonlinear Systems Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu 1 / 65 Outline 1 Nested Approximations 2 Trajectory PieceWise Linear (TPWL)
More informationNonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems
Nonlinear Model Reduction for Uncertainty Quantification in Large-Scale Inverse Problems Krzysztof Fidkowski, David Galbally*, Karen Willcox* (*MIT) Computational Aerospace Sciences Seminar Aerospace Engineering
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
16th International Conference on Finite Elements in Flow Problems DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology
More informationHigh Order Discontinuous Galerkin Methods for Aerodynamics
High Order Discontinuous Galerkin Methods for Aerodynamics Per-Olof Persson Massachusetts Institute of Technology Collaborators: J. Peraire, J. Bonet, N. C. Nguyen, A. Mohnot Symposium on Recent Developments
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 informationProgress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD
Progress in Mesh-Adaptive Discontinuous Galerkin Methods for CFD German Aerospace Center Seminar Krzysztof Fidkowski Department of Aerospace Engineering The University of Michigan May 4, 2009 K. Fidkowski
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationRICE UNIVERSITY. Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut
RICE UNIVERSITY Nonlinear Model Reduction via Discrete Empirical Interpolation by Saifon Chaturantabut A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy
More informationA reduced basis method and ROM-based optimization for batch chromatography
MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max
More informationEmpirical Interpolation of Nonlinear Parametrized Evolution Operators
MÜNSTER Empirical Interpolation of Nonlinear Parametrized Evolution Operators Martin Drohmann, Bernard Haasdonk, Mario Ohlberger 03/12/2010 MÜNSTER 2/20 > Motivation: Reduced Basis Method RB Scenario:
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
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 informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More 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 informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationGreedy control. Martin Lazar University of Dubrovnik. Opatija, th Najman Conference. Joint work with E: Zuazua, UAM, Madrid
Greedy control Martin Lazar University of Dubrovnik Opatija, 2015 4th Najman Conference Joint work with E: Zuazua, UAM, Madrid Outline Parametric dependent systems Reduced basis methods Greedy control
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More information1st Committee Meeting
1st Committee Meeting Eric Liu Committee: David Darmofal, Qiqi Wang, Masayuki Yano Massachusetts Institute of Technology Department of Aeronautics and Astronautics Aerospace Computational Design Laboratory
More informationApplication of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen
Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations Xiaoyue Shen Supervisor: Prof. David W. Zingg A thesis submitted in conformity with the requirements
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 informationOn Dual-Weighted Residual Error Estimates for p-dependent Discretizations
On Dual-Weighted Residual Error Estimates for p-dependent Discretizations Aerospace Computational Design Laboratory Report TR-11-1 Masayuki Yano and David L. Darmofal September 21, 2011 Abstract This report
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 informationReduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling
Reduced basis method for the reliable model reduction of Navier-Stokes equations in cardiovascular modelling Toni Lassila, Andrea Manzoni, Gianluigi Rozza CMCS - MATHICSE - École Polytechnique Fédérale
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 informationInterior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations
Interior penalty tensor-product preconditioners for high-order discontinuous Galerkin discretizations Will Pazner Brown University, 8 George St., Providence, RI, 9, U.S.A. Per-Olof Persson University of
More informationCFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki
CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
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 informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
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 informationMultigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations
Multigrid Algorithms for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations Khosro Shahbazi,a Dimitri J. Mavriplis b Nicholas K. Burgess b a Division of Applied
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationEntropy-Based Drag Error Estimation and Mesh. Adaptation in Two Dimensions
Entropy-Based Drag Error Estimation and Mesh Adaptation in Two Dimensions Krzysztof J. Fidkowski, Marco A. Ceze, and Philip L. Roe University of Michigan, Ann Arbor, MI 48109 This paper presents a method
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 informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationChapter 1. Introduction and Background. 1.1 Introduction
Chapter 1 Introduction and Background 1.1 Introduction Over the past several years the numerical approximation of partial differential equations (PDEs) has made important progress because of the rapid
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationarxiv: v1 [math.na] 9 Oct 2018
PRIMAL-DUAL REDUCED BASIS METHODS FOR CONVEX MINIMIZATION VARIATIONAL PROBLEMS: ROBUST TRUE SOLUTION A POSTERIORI ERROR CERTIFICATION AND ADAPTIVE GREEDY ALGORITHMS arxiv:1810.04073v1 [math.na] 9 Oct 2018
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
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 informationModel reduction for multiscale problems
Model reduction for multiscale problems Mario Ohlberger Dec. 12-16, 2011 RICAM, Linz , > Outline Motivation: Multi-Scale and Multi-Physics Problems Model Reduction: The Reduced Basis Approach A new Reduced
More informationThe hybridized DG methods for WS1, WS2, and CS2 test cases
The hybridized DG methods for WS1, WS2, and CS2 test cases P. Fernandez, N.C. Nguyen and J. Peraire Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics, MIT 5th High-Order
More informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationParameterized Partial Differential Equations and the Proper Orthogonal D
Parameterized Partial Differential Equations and the Proper Orthogonal Decomposition Stanford University February 04, 2014 Outline Parameterized PDEs The steady case Dimensionality reduction Proper orthogonal
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 informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationEfficient BDF time discretization of the Navier Stokes equations with VMS LES modeling in a High Performance Computing framework
Introduction VS LES/BDF Implementation Results Conclusions Swiss Numerical Analysis Day 2015, Geneva, Switzerland, 17 April 2015. Efficient BDF time discretization of the Navier Stokes equations with VS
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationReduced basis method for efficient model order reduction of optimal control problems
Reduced basis method for efficient model order reduction of optimal control problems Laura Iapichino 1 joint work with Giulia Fabrini 2 and Stefan Volkwein 2 1 Department of Mathematics and Computer Science,
More informationMinimal stabilization techniques for incompressible flows
Minimal stabilization tecniques for incompressible flows G. Lube 1, L. Röe 1 and T. Knopp 2 1 Numerical and Applied Matematics Georg-August-University of Göttingen D-37083 Göttingen, Germany 2 German Aerospace
More informationGerman Aerospace Center (DLR)
German Aerospace Center (DLR) AEROGUST M30 Progress Meeting 23-24 November 2017, Bordeaux Presented by P. Bekemeryer / J. Nitzsche With contributions of C. Kaiser 1, S. Görtz 2, R. Heinrich 2, J. Nitzsche
More informationEfficient and Flexible Solution Strategies for Large- Scale, Strongly Coupled Multi-Physics Analysis and Optimization Problems
University of Colorado, Boulder CU Scholar Aerospace Engineering Sciences Graduate Theses & Dissertations Aerospace Engineering Sciences Spring 1-1-2016 Efficient and Flexible Solution Strategies for Large-
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More information( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
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 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 informationOutput-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results
Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics: Overview and Recent Results Krzysztof J. Fidkowski University of Michigan David L. Darmofal Massachusetts Institute of
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 informationReduced-Order Greedy Controllability of Finite Dimensional Linear Systems. Giulia Fabrini Laura Iapichino Stefan Volkwein
Universität Konstanz Reduced-Order Greedy Controllability of Finite Dimensional Linear Systems Giulia Fabrini Laura Iapichino Stefan Volkwein Konstanzer Schriften in Mathematik Nr. 364, Oktober 2017 ISSN
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationAnisotropic mesh refinement for discontinuous Galerkin methods in two-dimensional aerodynamic flow simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2008; 56:2111 2138 Published online 4 September 2007 in Wiley InterScience (www.interscience.wiley.com)..1608 Anisotropic
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
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 informationC1.2 Ringleb flow. 2nd International Workshop on High-Order CFD Methods. D. C. Del Rey Ferna ndez1, P.D. Boom1, and D. W. Zingg1, and J. E.
C. Ringleb flow nd International Workshop on High-Order CFD Methods D. C. Del Rey Ferna ndez, P.D. Boom, and D. W. Zingg, and J. E. Hicken University of Toronto Institute of Aerospace Studies, Toronto,
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 informationContinuous adjoint based error estimation and r-refinement for the active-flux method
Continuous adjoint based error estimation and r-refinement for the active-flux method Kaihua Ding, Krzysztof J. Fidkowski and Philip L. Roe Department of Aerospace Engineering, University of Michigan,
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 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 informationThe Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation
The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation Y. Maday (LJLL, I.U.F., Brown Univ.) Olga Mula (CEA and LJLL) G.
More informationDiscontinuous Galerkin methods for compressible flows: higher order accuracy, error estimation and adaptivity
Discontinuous Galerkin methods for compressible flows: higher order accuracy, error estimation and adaptivity Ralf Hartmann Institute of Aerodynamics and Flow Technology German Aerospace Center (DLR) Lilienthalplatz
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 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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More information