Sparse Tensor Methods for PDEs with Stochastic Data II: Convergence Rates of gpcfem
|
|
- Bernadette Blake
- 5 years ago
- Views:
Transcription
1 Sparse Tensor Methods for PDEs with Stochastic Data II: Convergence Rates of gpcfem Christoph Schwab Seminar for Applied Mathematics ETH Zürich Albert Cohen (Paris, France) & Ron DeVore (Texas A& M) R. Andreev, M. Bieri and G.J. Gittelson ( UQ-Group ) (SAM,ETH) Workshop Tutorial Computation with Uncertainty IMA, Minneapolis, MN, USA, October 16 & 17, 2010 ERC Project STAHDPDE and Swiss National Science Foundation
2 Outline 1. Elliptic BVP with stochastic data 2. infinite - Dimensional Deterministic BVP 3. Regularity 4. Best N-term Convergence Rates of sgfem 5. Localization of active gpc-modes. Complexity. 6. Discretization in D. Sparse Tensor sgfem. 7. Parabolic and Hyperbolic Problems 8. Conclusions & References 1
3 Elliptic BVP with stochastic data Given: probability space (Ω,Σ,P) on data space X(D) L (D), V H 1 (D), random diffusion coefficient a(x,ω) L (Ω,dP;X(D)), deterministic source term f H 1 (D) = (H0 1(D)) =: V, (sbvp) Find u(x,ω) L 2 (Ω,dP;V) = L 2 (Ω,dP;H0 1(D)) L2 (Ω,dP) V such that [ˆ ] [ˆ ] E a(x, ) x u(x, ) x v(x, )dx = E f(x)v(x, )dx D for all v L 2 (Ω,dP;V) a L (Ω,dP;X(D)) and a max a(, ) a min > 0!u L 2 (Ω,dP;H0 1(D)). D 2
4 Karhunen-Loève expansion - separation of deterministic and stochastic variables - Example 1 (Karhunen-Loève Expansion) If a L 2 (Ω,dP;L (D)) then in L 2 (Ω,dP;L 2 (D)), a(x,ω) = E[a](x)+ m 1 ψ m (x)y m (ω) = E[a](x)+ m 1 λm ϕ m (x)y m (ω), (λ m,ϕ m ) m 1 eigensequence of (cpt., s.a.) covariance operator ˆ C[a] : L 2 (D) L 2 (D) (C[a]v)(x) := C a (x,x )v(x )dx v L 2 (D), D C a = E[(a E[a]) (a E[a])] Y m (ω) := 1 λmˆ D (a(x,ω) E[a](x))ϕ m (x)dx : Ω Γ m R m = 1,2,... 3
5 Karhunen-Loève expansion - convergence - a(x,ω) = E a (x)+ m 1 λm ϕ m (x)y m (ω) KL expansion converges in L 2 (D Ω), not necessarily in L (D Ω) To ensure L (D Ω) convergence, must - estimate decay rate of KL eigenvalues λ m : Schwab & Todor JCP (2006), - bound ϕ m L (D): Todor Diss ETH (2005), SINUM (2006) - assume: bounds for Y m L (Ω) 4
6 Karhunen-Loève expansion - eigenvalue estimates - Regularity of C a ensures decay of KL-eigenvalue sequence (λ m ) m 1 C a (x,x ) : D D R is piecewise analytic on D D if ex. smoothness partition D = {D j } J j=1 of D into a finite sequence of simplices D j such that D = J j=1 D j (1) and C a (x,x ) analytic in open neighbourhood of D j D j for any (j,j ). piecewise H t,t on D D if V a Hpw t,t (D D) := L 2 (D i,h t (D j )) i,j J 5
7 Karhunen-Loève expansion - eigenvalue estimates - (H,, ) separable Hilbert space, C K(H) compact, s.a., eigenpair sequence (λ m,φ m ) m 1. If C m B(H) is any operator of rank at most m, λ m+1 C C m B(H). (2) (e.g. Pinkus 1985: n-widths in Approx. Theory). 6
8 Proposition 2 (KL-eigenvalue decay) Karhunen-Loève expansion - eigenvalue estimates - (> exponential KL decay: Gaussian C a (x,x )) C a (x,x ) := σ 2 exp( γ x x 2 ) = 0 λ m c(γ,σ)/m! m 1 (exponential KL decay: Piecewise analytic C a (x,x )) C a pw analytic on D D = c > 0 0 λ m cexp( bm 1/d ) m 1 (Algebraic KL-eigenvalue decay for p.w. H t (D)-kernels) C a H t,t pw(d D) (t d/2) = 0 λ m cm t/d m 1 7
9 Karhunen-Loève expansion - eigenvalue estimates - 10 Largest eigenvalues of exp( x x 2 ), x ( 1,1) observed predicted Largest eigenvalues of exp( x x 1+δ ), x ( 1,1) δ=0.00, observed δ=0.00, predicted δ=0.25, observed δ=0.25, predicted δ=0.50, observed δ=0.50, predicted δ=0.75, observed δ=0.75, predicted
10 Karhunen-Loève expansion - eigenfunction estimates - Regularity of C a ensures L bounds for L 2 -scaled eigenfunctions (ϕ m ) m 1 Theorem 3 ( Schwab & Todor JCP 2006) Assume C a H t,t pw(d D) for t > d. Then Hence: δ > 0 ex.c(δ) > 0 s.t. m 1 : ϕ m L (D) C(δ)λ δ m. b m := λ1/2 m ϕ m L (D) inf x D E[a] C(δ)λ 1/2 δ m C(δ)m t/2d δ 9
11 Karhunen-Loève expansion - convergence rate - Conclusion: KL expansion of a(x,ω) L 2 (Ω,dP;L (D)) converges uniformly and exponentially on D Ω if - C a (x,x ) piecewise analytic - (Y m (ω)) m 1 uniformly bounded on Ω (e.g. Y m (ω) U( 1,1)) 10
12 Karhunen-Loève expansion - convergence rate - λ m 10 1 Largest 2000 Eigenvalues of Factorizable Kernel exp( 10 x x 2 ), x ( 1,1) computed eigenvalues λ m estimated eigenvalues 1/ Γ(m 1/3 /2) relative L 2 norm of remainder after truncation at level M L 2 Norm Decay of KL Series Remainder for Kernel e 10 x y 2 in 3D m M 11
13 Karhunen-Loève expansion - truncation from infinite to finite dimension M - > M N KL-truncation order a M (x,ω) := E[a](x)+ M m 1 (SBVP) with stochastic coefficient a(x,ω) div(a(x,ω) x u(x,ω)) = f(x) λm ϕ m (x)y m (ω) in L 2 (Ω,dP;H 1 (D)) (SBVP) M with truncated stochastic coefficient a M (x,ω) div(a M (x,ω) x u M (x,ω)) = f(x) in L 2 (Ω,dP;H 1 (D)) Theorem 4 If C a pw analytic and (Y m ) m 1 uniformly bounded, then δ > 0 ex. b,c(δ),m 0 > 0 such that (SBVP) M well-posed for M M 0 and u u M L2 (Ω;H 1 0(D)) { Cexp( bm 1/d ) M M 0 if C a pw analytic C(δ)M t/2d+1 δ M M 0 if C a H t,t pw(d D) 12
14 High dimensional deterministic bvp a M : D Ω R, a M (x,ω) = E[a](x)+ M m 1 λm ϕ m (x)y m (ω) Assumption (Y m ) m 1 independent, uniformly bounded family of rv s (e.g. Y m uniformly distributed in Γ m = I = ( 1/2,1/2), m = 1,2,3,...) Random variable Y m Parameter y m I m, m = 1,2,... (Y 1,Y 2,...,Y M ) y = (y 1,y 2,...,y M ) I M := I 1... I M dp(ω) ρ(y)dy = m 1ρ m (y m )dy m M ã M : D I M R, ã M (x,y) = E[a](x)+ λm ϕ m (x)y m m 1 13
15 Parametric Deterministic BVP a(x,ω) = E[a](x)+ j 1 ψ j (x)y j (ω) a(x,y) := E[a](x)+ j 1ψ j (x)y j, y = (y 1,y 2,...) U = ( 1,1) N, P(dω) = ρ(dy) = ρ(y)dy = j 1 ρ j (y j )dy j Parametric, deterministic problem, pointwise form: for every y U = ( 1,1) N, find u(x,y) such that x (a(x,y) x u(x,y)) = f(x) in D, u D = 0. 14
16 Parametric, deterministic problem, variational form: for every y U = ( 1,1) N, find U u u(y) V = H0 1 (D) such that ˆ ˆ v V : a(,y) u vdx = f(x)v(x)dx. } D {{} D =:b(y;u,v) Parametric deterministic problem, stochastic Galerkin form: Find u L 2 (U,dµ;V) such that ˆ ˆ ˆ ˆ v L 2 (U,dµ;V) : a(,y) u vdxµ(dy) = y U x D U D }{{} =:B(u,v)= y U b(y;u,v)µ(dy) f(x)v(x, y)dxµ(dy).
17 sgfem F: set of all ν = (ν j ) j 1 N N 0 such that only finitely many ν j are non-zero. Notation: ν F : ν := j 1 ν j = ν l 1 <, ν! = j 1 ν j! <, y ν = j 1 y ν j j y U. Tensorized (L 2 ( 1,1;1/2dx)-normalized) Legendre polynomials: L ν (y) = j 1L νj (y j ). Proposition 5: Assume Y j U( 1,1), i.e. ρ j (dy j ) = µ j (dy j ) := 1/2dy j, µ := j 1 µ j(dy j ). Then (L ν ) ν F is a complete orthonormal system in L 2 (U,dµ). 15
18 each v L 2 (U,dρ;V) has a representation ˆ v = v ν L ν, where v ν = v(,y)l ν (y)µ(dy) V ν F and there holds Parseval s equation v L2 (U,dµ;V) = ( v ν V ) l2 (F). U
19 sgfem For any finite subset Λ F, define X Λ := {v Λ (x,y) = ν Λv ν (x)l ν (y) ; v ν V} L (U,dµ;V) L 2 (U,dµ;V) where {L ν } ν F is the basis of Legendre polynomials. Galerkin approximation: u Λ = ν Λ u νl ν X Λ s.t. Cea s lemma: where C 1 := amax a min. u Λ X Λ such that B(u Λ,v Λ ) = F(v Λ ) v Λ X Λ. u u Λ L2 (U,dµ;V) C 1 inf v Λ X Λ u v Λ L2 (U,dµ;V), 16
20 sgfem Legendre expansion of u(x, y): u(x,y) = ν Fv ν (x)l ν (y) where v ν := ˆ U u(,y)l ν (y)µ(dy) V. u u Λ L2 (U,dµ;V) C 1 ( ν/ Λ v ν 2 V Choice of Λ? Summability of sequence α ν := v ν V, ν F. Lemma 6 (Stechkin): Let 0 < p q and assume α = (α ν ) ν F l p (F). If F N is the set of indices corresponding to N largest α ν, ( α ν q 1 ) q α lp (F)N r, ν/ F N where r := 1 p 1 q 0. )
21 Regularity b = (b j ) j=1, b j := ψ j L (D) a min, j = 1,2,... Theorem 7 (A-priori estimates of gpc coefficients) Assume ψ j L (D) b l1 (N) = κ < 1. a min j 1 Then v ν V ( ) ν! f V ν! bν, b ν = b ν 1 a min 1 bν , ν F. Estimate of v ν V by a real variable bootstrap argument. 18
22 Regularity Theorem 8 (Sequence Summability) For 0 < p 1, ( ν! ν! bν ) ν F l p (F) iff b l1 (N) < 1 b l p (N). 19
23 Convergence rates - semi-discretization in stochastic variable - Focus on approximation in L 2 (mean square, maximum likelihood) and L (worst case). Tool: Classical Legendre basis (P n ) n 0 : P n L ([ 1,1]) = P n (1) = 1. (3) L 2 normalized Legendre Sequence L n (t) = 2n+1P n (t), Note ˆ 1 1 L n (t) 2dt 2 = 1. L 0 = P 0 = 1. For ν F, recall tensorized polynomials P ν (y) := j 1 P νj (y j ) and L ν (y) := j 1L νj (y j ). (4) Note: 20
24 (L ν ) ν F orthonormal basis of L 2 (U,dµ), µ = j dy j every u L (U,dµ;V) L 2 (U,dµ;V) admits gpc expansions u(y) = ν F u ν P ν (y) = ν Fv ν L ν (y). (5) These gpc expansions converge in L 2 (U,dµ;V) with ˆ v ν := u(y)l ν (y)µ(dy) and u ν := j ) j 1(1+2ν U 1/2 v ν. (6) Estimates of u ν V, v ν V through complex analysis. Note v ν V u ν V.
25 Uniform Ellipticity Assumption: UEA(r, R) 0 < r R < such that x D and y U Example: 1-term KL expansion. 0 < r a(x,y) R <. (7) a(x,y) = ā(x)+yψ(x), ess inf x Dā(x) a min, ψ L (D) κa min = r = a min (1 κ).
26 Define A δ = {z C N : δ R(a(x,z)) a(x,z) 2R for every x D.} (8) Polydiscs: let ρ = (ρ j ) j 1 be a sequence of radii ρ j > 0 and define: U U ρ := j C : z j ρ j } C j 1{z N. (9) When a and ψ j are real valued, UEA(r,R) = x D and z U, 0 < r R(a(x,z)) a(x,z) 2R, (10) ρ = (ρ j ) j 1 is δ-admissible iff x D : ρ j ψ j (x) R(ā(x)) δ. (11) j 1 If ρ = (ρ j ) j 1 is δ-admissible, then U ρ A δ.
27 Lemma 9: Assume UEAC(r,R) for some 0 < r R <, ρ = (ρ j ) j 1 δ-admissible for some 0 < δ < r with ρ j > 1 for all j such that ν j 0. Then for any ν F hold a-priori bounds on gpc coeff s v ν V u ν V f V φ(ρ j )(2ν j +1)ρ ν j j, (12) δ where φ(t) := πt 2(t 1) for t > 1. j 1,ν j 0 Choice of δ-admissible weight sequence (ρ j ) j 1?
28 Legendre coefficient estimate in Lemma 6 and ν-dependent δ-admissible ρ : fix 1 < κ 2 such that (κ 1) j 1 ψ j L (D) r 8. (13) Choose J 0 as smallest integer such that j>j 0 ψ j L (D) r(κ 1) 18πκ, (14) Then define ρ j := Then ρ is δ admissible with δ = r 2 and κ if 1 j J 0, rν 2+ j 4 ν F ψ j L (D), j J 0 such that ν j 0, 1 j J 0 such that ν j = 0. v ν V u ν V ρ ν ν F. (15)
29 Theorem 10: Assume a(x,z) satisfies UEAC(r,R) for some 0 < r R < ( ψ j L (D)) j 1 l p (N) for some p < 1, Then for the same value of p ( u ν V ) ν F, ( v ν V ) ν F l p (F). The Legendre expansions (5) converge in L (U,dµ;V): i.e. S ΛN u(y) := ν Λ N u ν (x)p ν (y) = ν Λ N v ν (x)l ν (y) converge lim sup u(y) S ΛN u(y) V = 0 (16) N + y U for (Λ N ) N 1 F any sequence of finite sets which exhausts F. Moreover,...
30 1. L ( worst case ) convergence: If Λ N is the set of ν F corresponding to indices of N largest u ν V, then holds the convergence estimate sup y U u(y) S ΛN u(y) V ( u ν V ) lp (F)N r, r = 1 p 1 0. (17) 2. L 2 ( maximum likelihood ) convergence: If Λ N is the set of ν F corresponding to indices of N largest v ν V, then u S ΛN u L2 (U,dµ;V) ( v ν V ) lp (F)N r, r = 1 p (18)
31 Finding Λ l - I Problem: Theorem does not give a constructive way of finding sets Λ 1 Λ 2... Λ l...f 1. Adaptive sgfem, or 2. A-priori selection of the Λ l. Consider 2: given ρ = (ρ 1,ρ 2,...) > 0 s.t. 1 < ρ 1 ρ 2... ρ m and ε > 0, define Λ ε (ρ) := { ν F ρ ν ε } = {ν F ρ ν 1/ε} If our bounds on gpc coeff s v ν V are sharp, Λ ε (ρ) will contain essentially the #Λ ε (ρ) largest coefficients v ν V. Monotonicity: for any ρ, ρ Λ as above and any ε,ε > 0, ε ε implies Λ ε (ρ) Λ ε (ρ), ρ ρ implies Λ ε (ρ) Λ ε ( ρ). 21
32 Finding Λ l - II Consider comparison sequences ν σ l (N) such that Note: ρ 1 ν σ := {(m+1) σ : m = 1,2,...}, σ > 0. ρ 1 ν σ ν F : ρ ν νσ ν Scaling: given σ > 0, find { Λ ε (ν σ ) = Λ ε 1/σ(ν } 1 ) = ν N N 0 : (m+1) ν m ε 1/σ, ε 0. m N ν Λ ε (ν σ ) iff ν is a multiplicative partition of some integer x < ε 1/σ. 22
33 Complexity (finding Λ l ) Proposition 11 (R. Andreev 08) 1. Λ ε (ν σ ) can be localized in work and memory growing log-linearly in #Λ ε (ν σ ). 2. As ε 0, e 2 logx #Λ ε (ν σ ) x 2 with x = ε 1/σ, π(logx) 3/4 E.R. Canfield, Paul Erdös and C. Pomerance: On a Problem of Oppenheim concerning Factorisatio Numerorum J. Number Theory 17 (1983) 1 ff., G. Szekeres and P. Turán: Über das zweite Hauptproblem der Factorisatio Numerorum Acta Litt. Sci. Szeged 6 (1933)
34 Smoothness in D: Norms: Discretization in D W = {v V : v L 2 (D)} V. v W := v V + v W, Note: D convex = W = H 2 (D) V. v W = v L2 (D). Theorem 12 (Regularity) Assume f L 2 (D) and a(x,z) satisfies UEAC(r,R) for 0 < r R <. If for some 0 < p < 1 ( ψ j L (D)) j 1 l p (N) and ( ψ j L (D)) j 1 l p (N) then ( u ν W ) ν F,( v ν W ) ν F l p (F). 24
35 Discretization in D (Nonadaptive(!)) Finite Element Approximation in D: D Convex = inf w v h V Ch w W CM 1 d w W. v h V h Non-convex polyhedra D: W H 2 (D). Convergence as M = dim(v h ) reduced to inf v h V h w v h V C t M t w W some 0 < t < 1 d. See Brenner & Scott: Intro FEM, Babuška, Kellogg and Pitkäranta (1979) 25
36 Discretization in D (Nonadaptive(!)) Finite Element Approximation in D: Theorem 13: Assume FE spaces (V h ) h>0 in D R d have approximation order 0 < t 1/d, a satisfies UEAC(r,R) and ( ψ j W 1, (D)) j=1 lp (N) some 0 < p 1. Then (i) With Λ N the set of indices corresponding to N largest u ν V, there exist finite element spaces V ν of dimensions M ν, ν Λ N, such that u(y) sup y U ũ ν P ν (y) V CN min{r,t} dof, r := 1 p 1, ν Λ N where N dof = ν Λ N M ν and C = ( C t + ( u ν V ) lp (F)) ( u ν W ) lp (F). 26
37 (ii) With Λ N the set of indices corresponding to N largest v ν V, there exist finite element spaces V ν of dimension M ν, ν Λ N, such that u ṽ ν L ν L2 (U,dµ;V) CN min{r,t} dof, r := 1 p 1 2, ν Λ N where N dof = ν Λ N M ν and C = ( C 2 t + ( v ν V ) 2 l p (F) )1 2 ( v ν W ) lp (F). Here, t ν, ũ ν and ṽ ν are V-projections of t ν, u ν and v ν, respectively, onto V ν.
38 Numerical Example (R. Andreev) D = ( 1,1), E a (x) = 5+x, C a (x,x ) = min{x,x }+1 2 H 1,1 (D D). KL Eigenpairs: λ m = Algebraic KL decay with rate 2. 8 π 2 (2m 1) 2, ψ m(x) = sin((x+1)/ 2λ m ) Bieri, Andreev and CS. SISC (2010) 27
39 Error vs. N total 10 1 error MC, L 2 (H 1 ) FTP, L 2 (H 1 ) STP, L 2 (H 1 ) FTP, L 2 (L 2 ) STP, L 2 (L 2 ) ndof
40 Heat Eqn w. Random Conductivity (w. V.H. Hoang) Q T = I D. In Q T, u t (a(x,ω) u) = g(t,x), u D I = 0, u t=0 = h(x). V = H 1 0(D) H = L 2 (D) H V = H 1 (D) X = L 2 (I;V) H 1 (I;V ) and Y = L 2 (I;V) H. In X and Y norms X and Y, for u X and v = (v 1,v 2 ) Y given by u X = ( u 2 L 2 (I;V) + u 2 H 1 (I;V ) )1/2 and v Y = ( v 1 2 L 2 (I;V) + v 2 2 H )1/2. 28
41 Heat Eqn w. Random Conductivity Parametric Deterministic Problem: u t (t,x,y) x [a(x,y) x u(t,x,y)] = g(t,x) in Q T, u(t,x,y) D I = 0, u t=0 = h(x), where, for y = (y 1,y 2,...) U a(x,y) = ā(x)+ y j ψ j (x). Pointwise parametric (x, t)-variational Formulation [Sc. & R. Stevenson2009]: given f Y, find u(y) : U y X such that j=1 b(y;u,v) = f(v), v = (v 1,v 2 ) Y. where ˆ b(y;u,(v 1,v 2 )) = du ˆ ˆ I dt,v 1(t, ) H dt+ a(x,y) u(t,x) v 1 (t,x)dxdt+ w(0),v 2 H. D I ˆ f(v) = g(t),v 1 (t) H dt+ h,v 2 H, v = (v 1,v 2 ) Y. I 29
42 Heat Eqn w. Random Conductivity Galerkin parametric (x, t) variational formulation: Bochner spaces X = L 2 (U,dρ;X) L 2 (U,dρ) X and Y = L 2 (U,dρ;Y) L 2 (U,dρ) Y. where Find u X such that B(u,v) = F(v) for all v Y B(u,v) = ˆ ˆ b(y, u, v)ρ(dy) and F(v) = U U f(v)ρ(dy). 30
43 Heat Eqn w. Random Conductivity gpc Approximation: Λ F, N = #Λ <. X Λ := {u Λ (t,x,y) = ν Λu ν (t,x)l ν (y) : u ν X} X, and Y Λ := {v Λ (t,x,y) = ν Λv ν (t,x)l ν (y) : v ν X} Y. In the Legendre basis (L ν ) ν F, we write v 1Λ (t,x,y) = ν Λ v 1ν (t,x)l ν (y) and v 2Λ (x,y) = ν Λv 2ν (x)l ν (y), respectively, where v ν = (v 1ν,v 2ν ) Y for all ν F. Galerkin approximation: u Λ X Λ : B(u Λ,v Λ ) = F(v Λ ) v Λ Y Λ. (19) 31
44 Heat Eqn w. Random Conductivity Theorem 14: Λ F, (19) is coupled system of N = #Λ parabolic equations, admits a unique solution u Λ X Λ a-priori error bound: ( u u Λ X c ν/ Λ u ν 2 X) 1/2. u ν X Legendre coefficients of the solution of the parametric deterministic problem, c independent of Λ. ( ψ j L (D)) j 1 l p (N) for some 0 < p 1 = ( u ν ) ν F l p (F) 32
45 and (Λ N ) N N F of index sets with #Λ N such that the solutions u ΛN of the Galerkin semidiscretized problems (19) satisfy u u ΛN X CN r, r = 1 p 1 2. Space-Time (x, t) Discretization: regularity, (x,t) wavelets, optimal error bounds in terms of N dof Analogous results for wave equations in random media
46 Conclusion gpc: Transform SPDE into parametric, deterministic PDE on U = ( 1,1) N ST approximation by gpc Ψ of Regularity: W V smoothness space in (x,t) U y u(y, ) V ( u ν W ) ν F l p (F), 0 < p 1. p = 1: ST-rate equal to MLMC convergence rate (A. Barth, CS, N. Zollinger 2010) Sparse tensorization of D and Ω Galerkin discretizations: Dimension independent convergence rates in terms of N total. 33
47 requires mixed regularity in both, D and U: if { ψ j W 1, (D)} j=1 l p (N) then { u ν W : ν F} l p (F) L 2 (mean square) and L (worst case) convergence rates. Adaptivity essential to avoid curse of dimensionality. Diffusion problems in physical dimension d = 2, with λ m m 2.2 and suppλ {1,...,1500} on PC (R. Andreev, M. Bieri and CS (SISC 2010)) (superior to Smolyak construction). Sparse composite collocation on lattices of points in U = ( 1,1) N, U = R N (M. Bieri Dissertation ETH (2009) and SINUM (2010)) Lognormal Case (C.J. Gittelson (M 3AS 2010)) Heat and Wave Eqn. (C.S. & V. Hoang 2010), Eigenvalue Problems (C.S. & R. Andreev) Adaptive, Convergent Algorithm (w. A. Cohen, R. DeVore 2010)
48 References N. Wiener (Am. J. Math. 1938), Cameron and Martin (Ann. Math. 1947) G. E. Karniadakis et al. ( ), Babuška, Nobile, Tempone etal. (SINUM 2004, 2006, 2008) (CMAME 2005) P. Frauenfelder, R.A. Todor and CS (CMAME 2005) R.A. Todor (SINUM 2005 and Diss. ETH 2005) CS and R.A. Todor (JCP 2006 and IMAJNA 2007) R. Andreev, M. Bieri and CS (SISC 2010) M. Bieri (SINUM 2010) V. H. Hoang and CS (2010) A. Cohen, R. DeVore and CS (2009) (2010) C.J. Gittelson Dissertation ETH (2010). reports 34
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationFoundations of the stochastic Galerkin method
Foundations of the stochastic Galerkin method Claude Jeffrey Gittelson ETH Zurich, Seminar for Applied Mathematics Pro*oc Workshop 2009 in isentis Stochastic diffusion equation R d Lipschitz, for ω Ω,
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationMultilevel accelerated quadrature for elliptic PDEs with random diffusion. Helmut Harbrecht Mathematisches Institut Universität Basel Switzerland
Multilevel accelerated quadrature for elliptic PDEs with random diffusion Mathematisches Institut Universität Basel Switzerland Overview Computation of the Karhunen-Loéve expansion Elliptic PDE with uniformly
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationMulti-Element Probabilistic Collocation Method in High Dimensions
Multi-Element Probabilistic Collocation Method in High Dimensions Jasmine Foo and George Em Karniadakis Division of Applied Mathematics, Brown University, Providence, RI 02912 USA Abstract We combine multi-element
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationSparse Quadrature Algorithms for Bayesian Inverse Problems
Sparse Quadrature Algorithms for Bayesian Inverse Problems Claudia Schillings, Christoph Schwab Pro*Doc Retreat Disentis 2013 Numerical Analysis and Scientific Computing Disentis - 15 August, 2013 research
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationAdaptive Collocation with Kernel Density Estimation
Examples of with Kernel Density Estimation Howard C. Elman Department of Computer Science University of Maryland at College Park Christopher W. Miller Applied Mathematics and Scientific Computing Program
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationarxiv: v2 [math.na] 8 Sep 2017
arxiv:1704.06339v [math.na] 8 Sep 017 A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations Juan Galvis O. Andrés Cuervo September 3, 018 Abstract We use a Monte
More informationSparse polynomial chaos expansions in engineering applications
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Sparse polynomial chaos expansions in engineering applications B. Sudret G. Blatman (EDF R&D,
More informationQuasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient
Quasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient F. Nobile, L. Tamellini, R. Tempone, F. Tesei CSQI - MATHICSE, EPFL, Switzerland Dipartimento di
More informationApproximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity
Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Marcus Sarkis Jointly work with Juan Galvis USC09 Marcus Sarkis (WPI) SPDE without Uniform Ellipticity USC09
More informationJ. Beck 1, F. Nobile 2, L. Tamellini 2 and R. Tempone 1
ESAIM: PROCEEDINGS, October 2011, Vol. 33, p. 10-21 C. Dobrzynski, T. Colin & R. Abgrall, Editors IMPLEMENTATION OF OPTIMAL GALERKIN AND COLLOCATION APPROXIMATIONS OF PDES WITH RANDOM COEFFICIENTS, J.
More informationarxiv: v1 [math.na] 14 Dec 2018
A MIXED l 1 REGULARIZATION APPROACH FOR SPARSE SIMULTANEOUS APPROXIMATION OF PARAMETERIZED PDES NICK DEXTER, HOANG TRAN, AND CLAYTON WEBSTER arxiv:1812.06174v1 [math.na] 14 Dec 2018 Abstract. We present
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationMulti-Level Monte-Carlo Finite Element Methods for stochastic elliptic variational inequalities
Multi-Level Monte-Carlo Finite Element Methods for stochastic elliptic variational inequalities R. Kornhuber and C. Schwab and M. Wolf Research Report No. 2013-12 April 2013 Seminar für Angewandte Mathematik
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More informationSolving the steady state diffusion equation with uncertainty Final Presentation
Solving the steady state diffusion equation with uncertainty Final Presentation Virginia Forstall vhfors@gmail.com Advisor: Howard Elman elman@cs.umd.edu Department of Computer Science May 6, 2012 Problem
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationLow-rank techniques applied to moment equations for the stochastic Darcy problem with lognormal permeability
Low-rank techniques applied to moment equations for the stochastic arcy problem with lognormal permeability Francesca Bonizzoni 1,2 and Fabio Nobile 1 1 CSQI-MATHICSE, EPFL, Switzerland 2 MOX, Politecnico
More informationLecture 1: Center for Uncertainty Quantification. Alexander Litvinenko. Computation of Karhunen-Loeve Expansion:
tifica Lecture 1: Computation of Karhunen-Loeve Expansion: Alexander Litvinenko http://sri-uq.kaust.edu.sa/ Stochastic PDEs We consider div(κ(x, ω) u) = f (x, ω) in G, u = 0 on G, with stochastic coefficients
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 1 / 23 Lecture outline
More informationStochastic methods for solving partial differential equations in high dimension
Stochastic methods for solving partial differential equations in high dimension Marie Billaud-Friess Joint work with : A. Macherey, A. Nouy & C. Prieur marie.billaud-friess@ec-nantes.fr Centrale Nantes,
More informationSTOCHASTIC SAMPLING METHODS
STOCHASTIC SAMPLING METHODS APPROXIMATING QUANTITIES OF INTEREST USING SAMPLING METHODS Recall that quantities of interest often require the evaluation of stochastic integrals of functions of the solutions
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationCollocation methods for uncertainty quantification in PDE models with random data
Collocation methods for uncertainty quantification in PDE models with random data Fabio Nobile CSQI-MATHICSE, EPFL, Switzerland Acknowlegments: L. Tamellini, G. Migliorati (EPFL), R. Tempone, (KAUST),
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationCBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods. Crowder, Adam J and Powell, Catherine E
CBS Constants and Their Role in Error Estimation for Stochastic Galerkin Finite Element Methods Crowder, Adam J and Powell, Catherine E 2017 MIMS EPrint: 2017.18 Manchester Institute for Mathematical Sciences
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationSolving the stochastic steady-state diffusion problem using multigrid
IMA Journal of Numerical Analysis (2007) 27, 675 688 doi:10.1093/imanum/drm006 Advance Access publication on April 9, 2007 Solving the stochastic steady-state diffusion problem using multigrid HOWARD ELMAN
More informationTECHNISCHE UNIVERSITÄT BERLIN
TECHNISCHE UNIVERSITÄT BERLIN Adaptive Stochastic Galerkin FEM with Hierarchical Tensor Representations Martin Eigel Max Pfeffer Reinhold Schneider Preprint 2015/29 Preprint-Reihe des Instituts für Mathematik
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationA sparse grid stochastic collocation method for elliptic partial differential equations with random input data
A sparse grid stochastic collocation method for elliptic partial differential equations with random input data F. Nobile R. Tempone C. G. Webster June 22, 2006 Abstract This work proposes and analyzes
More informationarxiv: v2 [math.na] 8 Apr 2017
A LOW-RANK MULTIGRID METHOD FOR THE STOCHASTIC STEADY-STATE DIFFUSION PROBLEM HOWARD C. ELMAN AND TENGFEI SU arxiv:1612.05496v2 [math.na] 8 Apr 2017 Abstract. We study a multigrid method for solving large
More informationACMAC s PrePrint Repository
ACMAC s PrePrint Repository Analytic Regularity and GPC Approximation for Control Problems Constrained by Linear Parametric Elliptic and Parabolic PDEs Angela Kunoth and Christoph Schwab Original Citation:
More informationhp-dg timestepping for time-singular parabolic PDEs and VIs
hp-dg timestepping for time-singular parabolic PDEs and VIs Ch. Schwab ( SAM, ETH Zürich ) (joint with O. Reichmann) ACMAC / Heraklion, Crete September 26-28, 2011 September 27, 2011 Motivation Well-posedness
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationGiovanni Migliorati. MATHICSE-CSQI, École Polytechnique Fédérale de Lausanne
Analysis of the stability and accuracy of multivariate polynomial approximation by discrete least squares with evaluations in random or low-discrepancy point sets Giovanni Migliorati MATHICSE-CSQI, École
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationA Reduced Basis Approach for Variational Problems with Stochastic Parameters: Application to Heat Conduction with Variable Robin Coefficient
A Reduced Basis Approach for Variational Problems with Stochastic Parameters: Application to Heat Conduction with Variable Robin Coefficient Sébastien Boyaval a,, Claude Le Bris a, Yvon Maday b, Ngoc Cuong
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationAnalysis and finite element approximations of stochastic optimal control problems constrained by stochastic elliptic partial differential equations
Retrospective Theses and issertations Iowa State University Capstones, Theses and issertations 2008 Analysis and finite element approximations of stochastic optimal control problems constrained by stochastic
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationBath Institute For Complex Systems
BICS Bath Institute for Complex Systems Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and its Application to Multilevel Monte Carlo Methods Julia Charrier, Robert Scheichl and
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationCOMPRESSIVE SENSING PETROV-GALERKIN APPROXIMATION OF HIGH-DIMENSIONAL PARAMETRIC OPERATOR EQUATIONS
COMPRESSIVE SENSING PETROV-GALERKIN APPROXIMATION OF HIGH-DIMENSIONAL PARAMETRIC OPERATOR EQUATIONS HOLGER RAUHUT AND CHRISTOPH SCHWAB Abstract. We analyze the convergence of compressive sensing based
More informationSparse Grids. Léopold Cambier. February 17, ICME, Stanford University
Sparse Grids & "A Dynamically Adaptive Sparse Grid Method for Quasi-Optimal Interpolation of Multidimensional Analytic Functions" from MK Stoyanov, CG Webster Léopold Cambier ICME, Stanford University
More informationAnalytic regularity and polynomial approximation of stochastic, parametric elliptic multiscale PDEs
Eidgenössische Technische Hochschule Zürich Ecole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich Analytic regularity and polynomial approximation
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 informationEFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS
EFFICIENT STOCHASTIC GALERKIN METHODS FOR RANDOM DIFFUSION EQUATIONS DONGBIN XIU AND JIE SHEN Abstract. We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We
More informationNumerical Approximation of Stochastic Elliptic Partial Differential Equations
Numerical Approximation of Stochastic Elliptic Partial Differential Equations Hermann G. Matthies, Andreas Keese Institut für Wissenschaftliches Rechnen Technische Universität Braunschweig wire@tu-bs.de
More informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationNovember 7, 2002 GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS
November 7, 00 GALERKIN FINITE ELEMENT APPROXIMATIONS OF STOCHASTIC ELLIPTIC PARTIAL IFFERENTIAL EQUATIONS IVO BABUŠKA, RAÚL TEMPONE AN GEORGIOS E. ZOURARIS Abstract. We describe and analyze two numerical
More informationSparse deterministic approximation of Bayesian inverse problems
Home Search Collections Journals About Contact us My IOPscience Sparse deterministic approximation of Bayesian inverse problems This article has been downloaded from IOPscience. Please scroll down to see
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationarxiv: v3 [math.na] 18 Sep 2016
Infinite-dimensional integration and the multivariate decomposition method F. Y. Kuo, D. Nuyens, L. Plaskota, I. H. Sloan, and G. W. Wasilkowski 9 September 2016 arxiv:1501.05445v3 [math.na] 18 Sep 2016
More informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationLecture I: Asymptotics for large GUE random matrices
Lecture I: Asymptotics for large GUE random matrices Steen Thorbjørnsen, University of Aarhus andom Matrices Definition. Let (Ω, F, P) be a probability space and let n be a positive integer. Then a random
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationOptimal polynomial approximation of PDEs with stochastic coefficients by Galerkin and collocation methods
Optimal polynomial approximation of PDEs with stochastic coefficients by Galerkin and collocation methods Fabio Nobile MOX, Department of Mathematics, Politecnico di Milano Seminaire du Laboratoire Jacques-Louis
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
More informationKarhunen-Loève Approximation of Random Fields by Generalized Fast Multipole Methods
Karhunen-Loève Approximation of Random Fields by Generalized Fast Multipole Methods C. Schwab and R.A. Todor Research Report No. 26- January 26 Seminar für Angewandte Mathematik Eidgenössische Technische
More informationUncertainty Quantification in Discrete Fracture Network Models
Uncertainty Quantification in Discrete Fracture Network Models Claudio Canuto Dipartimento di Scienze Matematiche, Politecnico di Torino claudio.canuto@polito.it Joint work with Stefano Berrone, Sandra
More informationarxiv: v1 [math.na] 3 Apr 2019
arxiv:1904.02017v1 [math.na] 3 Apr 2019 Poly-Sinc Solution of Stochastic Elliptic Differential Equations Maha Youssef and Roland Pulch Institute of Mathematics and Computer Science, University of Greifswald,
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationAnalysis and Computation of Hyperbolic PDEs with Random Data
Analysis and Computation of Hyperbolic PDEs with Random Data Mohammad Motamed 1, Fabio Nobile 2,3 and Raúl Tempone 1 1 King Abdullah University of Science and Technology, Thuwal, Saudi Arabia 2 EPFL Lausanne,
More informationResearch Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based Karhunen-Loève Expansion
Mathematical Problems in Engineering Volume 2012, Article ID 215109, 15 pages doi:10.1155/2012/215109 Research Article Multiresolution Analysis for Stochastic Finite Element Problems with Wavelet-Based
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More informationExperiences with Model Reduction and Interpolation
Experiences with Model Reduction and Interpolation Paul Constantine Stanford University, Sandia National Laboratories Qiqi Wang (MIT) David Gleich (Purdue) Emory University March 7, 2012 Joe Ruthruff (SNL)
More informationA Preconditioned Recycling GMRES Solver for Stochastic Helmholtz Problems
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 2, pp. 342-353 Commun. Comput. Phys. August 2009 A Preconditioned Recycling GMRES Solver for Stochastic Helmholtz Problems Chao Jin 1 and Xiao-Chuan
More informationHeat kernels of some Schrödinger operators
Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Awards Screening Test. February 25, Time Allowed: 90 Minutes Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Awards Screening Test February 25, 2006 Time Allowed: 90 Minutes Maximum Marks: 40 Please read, carefully, the instructions on the following page before you
More informationUncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs.
Uncertainty Quantification and hypocoercivity based sensitivity analysis for multiscale kinetic equations with random inputs Shi Jin University of Wisconsin-Madison, USA Shanghai Jiao Tong University,
More information