Finite element approximation of the stochastic heat equation with additive noise
|
|
- Clinton Boyd
- 5 years ago
- Views:
Transcription
1 p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson
2 p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u() = u Abstract framework Mild solution Finite element approximation Strong convergence Weak convergence
3 p. 3/32 Co-workers Yubin Yan Mihály Kovács Fredrik Lindgren Matthias Geissert (Darmstadt)
4 p. 4/32 Stochastic heat equation The stochastic heat equation is u(x,t) t u(x, t) = f(u(x, t)) + g(u(x, t))ẇ(x, t), x Ω Rd, t > u(x, t) =, x Ω, t > u(x, ) = u (x), x Ω I know four approaches to SPDE: Martingale measure approach: Walsh 1986 Variational approach: Pardoux 1972, Krylov and Rozowski 1979 Semigroup approach: Da Prato and Zabczyk 1992 Wick product approach: Øksendal 1996 We use the semigroup approach. It is one way of giving a rigorous meaning to the equation.
5 p. 5/32 Semigroup approach First we change to stochastic notation: Ω D R d, x ξ D R d, u(x, t) X t (ξ), so that the stochastic heat equation can be written in the form Define dx t X t dt = f(x t ) dt + g(x t ) dw t H = L 2 (D), A =, D(A) = H 2 (D) H 1 (D) so that we can write { dx t + AX t dt = f(x t ) dt + g(x t ) dw t X = x
6 p. 6/32 Semigroup approach Let {e ta } t be the analytic semigroup generated by A, that is, u(t) = e ta u = e tλ j (u, φ j )φ j j=1 is the solution of the homogeneous heat equation { u + Au =, t > u() = u The solution of the non-homogeneous equation { u + Au = f, t > u() = u is then given by u(t) = e ta u + e (t s)a f(s) dt
7 p. 7/32 Analytic semigroup e ta v = e tλ j (v, φ j )φ j j=1 can be extended as a holomorphic function of t smoothing property: A β e ta v Ct β v, β A 1/2 e sa v 2 ds C v 2 Theory of semigroups of bounded linear operators
8 p. 8/32 Semigroup approach The stochastic heat equation: { dx t + AX t dt = f(x t ) dt + g(x t ) dw t X = x It is given a rigorous meaning as the integral equation: X t = e ta x + e (t s)a f(x s ) ds + e (t s)a g(x s ) dw s Such a solution is called a mild solution. Other approaches: different formulations and solution concepts. We need to define the stochastic integral.
9 p. 9/32 Semigroup approach Linear equation with additive noise : { dx + AX dt = dw, t > X() = X H = L 2 (D),, (, ), D R d, bounded domain A =, D(A) = H 2 (D) H(D) 1 X(t), H-valued stochastic process on probability space (Ω, F, P) W(t), H-valued Wiener process with filtration F t E(t) = e ta, analytic semigroup generated by A Mild solution (stochastic convolution): X(t) = E(t)X + E(t s) dw(s), t
10 p. 1/32 Q-Wiener process covariance operator Q : H H, self-adjoint, positive definite, bounded, linear operator Qe j = γ j e j, γ j >, {e j } j=1 ON basis β j (t), independent identically distributed, real-valued, Brownian motions W(t) = γ 1/2 j β j (t)e j j=1 Two important cases: Tr(Q) <. W(t) converges in L 2 (Ω, H): E γ 1/2 2 j β j (t)e j = γ j E ( β j (t) 2) = t j=1 j=1 γ j = t Tr(Q) < Q = I, white noise. W(t) is not H-valued, since Tr(I) =, but converges in a weaker sense j=1
11 p. 11/32 Q-Wiener process If Tr(Q) < : W() = continuous paths independent increments Gaussian law P (W(t) W(s)) 1 = N(, (t s)q), s t W(t) generates a filtration F t so that it becomes a square integrable H-valued martingale. We can integrate with respect to W.
12 p. 12/32 Stochastic integral X(t) = E(t)X + E(t s) dw(s), t We can define the stochastic integral E B(s)Q 1/2 2 HS ds < Isometry property E B(s) dw(s) 2 = E B(s)Q 1/2 2 HS ds B(s) dw(s) if Hilbert-Schmidt operator B: B 2 HS = Bϕ l 2 <, {ϕ l } arbitrary orthonormal basis in H l=1 Da Prato and Zabczyk, Stochastic Equations in Infinite Dimensions
13 p. 13/32 Regularity v β = A β/2 v, Ḣ β = D(A β/2 ), β R for example: Ḣ 2 = H 2 (D) H(D), 1 Ḣ 1 = H(D) 1 v 2 L 2 (Ω,Ḣβ ) = E( v 2 β) = A β/2 v 2 dξ dp(ω), β R Ω Theorem. If A (β 1)/2 Q 1/2 HS < for some β, then ( ) X(t) L2 (Ω,Ḣ β ) C X L2 (Ω,Ḣ β ) + A(β 1)/2 Q 1/2 HS D Two cases: If Q 1/2 2 HS = Tr(Q) <, then β = 1 If Q = I, d = 1, A = 2 ξ 2, then A (β 1)/2 HS < for β < 1/2 A (β 1)/2 2 HS = j λ (1 β) j j j (1 β)2/d < iff d = 1, β < 1/2
14 p. 14/32 Proof with X = X(t) 2 L 2 (Ω,Ḣ β ) = E ( = = = C A β/2 E(t s)dw(s) 2 ) A β/2 E(s)Q 1/2 2 HS ds A 1/2 E(s)A (β 1)/2 Q 1/2 2 HS ds k=1 A 1/2 E(s)A (β 1)/2 Q 1/2 φ k 2 ds A (β 1)/2 Q 1/2 φ k 2 k=1 = C A (β 1)/2 Q 1/2 2 HS A 1/2 E(s)v 2 ds C v 2
15 p. 15/32 Regularity of W Theorem. If A (β 1)/2 Q 1/2 HS < for some β, then W(t) L 2 (Ω, Ḣ β 1) ) L 2 (Ω, Ḣ 1 ) Proof. W(t) = dw(s) ( W(t) 2 L 2 (Ω,Ḣ β 1 ) = E A (β 1)/2 dw(s) 2 ) = A (β 1)/2 Q 1/2 2 HS ds = t A (β 1)/2 Q 1/2 2 HS < If Q 1/2 2 HS = Tr(Q) <, then β = 1 and W(t) L 2(Ω, H) If Q = I, d = 1, A = 2 x 2, then A (β 1)/2 HS < for β < 1/2
16 p. 16/32 The finite element method Deterministic heat equation, with u t = u t : { u t u = f, x D, t > u =, x D, t > (u, v) = D uv dξ (u t, v) ( u, v) }{{} =( u, v) = (f, v), v H 1 (D) = Ḣ 1 weak form: { u(t) H 1 (D), u() = u (u t, v) + ( u, v) = (f, v), v H 1 (D)
17 p. 17/32 The finite element method weak form: { u(t) H 1 (D), u() = u (u t, v) + ( u, v) = (f, v), v H 1 (D) triangulations {T h } <h<1, mesh size h function spaces {S h } <h<1, S h H 1 (D) = Ḣ1 { u h (t) S h, (u h (), v) = (u, v), v S h (u h,t, v) + ( u h, v) = (f, v), v S h (u h,t, v) + ( u h, v) }{{} =(A h u h,v) { u h,t + A h u h = P h f, t > u h () = P h u = (f, v), v S h }{{} =(P h f,v) continuous piecewise linear functions the same abstract framework
18 p. 18/32 The finite element method triangulations {T h } <h<1, mesh size h finite element spaces {S h } <h<1 S h H(D) 1 = Ḣ1 S h continuous piecewise linear functions A h : S h S h, discrete Laplacian, (A h ψ, χ) = ( ψ, χ), χ S h P h : L 2 S h, orthogonal projection, (P h f, χ) = (f, χ), χ S h { Xh (t) S h, X h () = P h X dx h + A h X h dt = P h dw P h W(t) is Q h -Wiener process with Q h = P h QP h. Mild solution, with E h (t) = e ta h, X h (t) = E h (t)p h X + E h (t s)p h dw(s)
19 p. 19/32 Error estimates for the deterministic problem { ut + Au =, t > u() = v u(t) = E(t)v { uh,t + A h u h =, t > u h () = P h v u h (t) = E h (t)p h v Denote F h (t)v = E h (t)p h v E(t)v, v β = A β/2 v We have, for β 2, F h (t)v Ch β v β, t ( 1/2 F h (s)v ds) 2 Ch β v β 1, t V. Thomée, Galerkin Finite Element Methods for Parabolic Problems
20 p. 2/32 Strong convergence in L 2 norm Theorem. If A (β 1)/2 Q 1/2 HS < for some β [, 2], then X h (t) X(t) L2 (Ω,H) Ch β( X L2 (Ω,Ḣ β ) + A(β 1)/2 Q 1/2 HS ) Recall: X h (t) X(t) L2 (Ω,H) = ( ) 1/2 E( X h (t) X(t) 2 ) Two cases: If Q 1/2 2 HS = Tr(Q) <, then the convergence rate is O(h). If Q = I, d = 1, A = 2 ξ 2, then the rate is almost O(h 1/2 ). No result for Q = I, d 2.
21 p. 21/32 Strong convergence: proof X(t) = E(t)X + X h (t) = E h (t)p h X + F h (t) = E h (t)p h E(t) E(t s) dw(s) E h (t s) P h dw(s) X h (t) X(t) = F h (t)x + F h (t s) dw(s) = e 1 (t) + e 2 (t) F h (t)x Ch β X β (deterministic error estimate) = e 1 (t) L2 (Ω,H) Ch β X L2 (Ω,Ḣ β )
22 p. 22/32 Strong convergence: proof E B(s) dw(s) 2 t = E B(s)Q 1/2 2 HS ds (isometry) ( 1/2 F h (s)v ds) 2 Ch β v β 1, with v = Q 1/2 ϕ l (deterministic) = e 2 (t) 2 L 2 (Ω,H) = E F h (t s) dw(s) 2 = = l=1 F h (t s)q 1/2 ϕ l 2 ds C F h (t s)q 1/2 2 HS ds h 2β Q 1/2 ϕ l 2 β 1 l=1 = Ch 2β A (β 1)/2 Q 1/2 ϕ l 2 = Ch 2β A (β 1)/2 Q 1/2 2 HS l=1 If Tr(Q) <, we may choose β = 1, otherwise β < 1.
23 p. 23/32 Time discretization { dx + AX dt = dw, t > X() = X The implicit Euler method: k = t, t n = nk, W n = W(t n ) W(t n 1 ) { X n h S h, X h = P h X X n h X n 1 h + ka h X n h = P h W n, X n h = E kh X n 1 h + E kh P h W n, E kh = (I + ka h ) 1 X n h = E n khp h X + n j=1 E n j+1 kh P h W j X(t n ) = E(t n )X + n E(t n s) dw(s)
24 p. 24/32 Error estimates for deterministic parabolic problem Denote F n = E n kh P h E(t n ) We have the following estimates for β 2: F n v C(k β/2 + h β ) v β ( n k F j v 2) 1/2 C(k β/2 + h β ) v β 1 j=1
25 p. 25/32 Strong convergence Theorem. If A (β 1)/2 Q 1/2 HS < for some β [, 2], then, with e n = Xh n X(t n), ) e n L2 (Ω,H) C(k β/2 + h β ) ( X L2(Ω,Ḣ β) + A(β 1)/2 Q 1/2 HS J. Printems (21) (only time-discretization) Y. Yan, BIT (24), SIAM J. Numer. Anal (25)
26 p. 26/32 Weak convergence If A (β 1)/2 Q 1/2 HS < for some β [, 2] and g C 2 b (H,R), then ( ) E g(xh) n g(x(t n )) C(k β + h 2β ) twice the rate of strong convergence Debussche and Printems 27 Geissert, Kovács, and Larsson 27 (with somewhat stronger assumption on g and only spatial discretization) the proof is based on Itô s formula and the Kolmogorov equation in infinite dimension
27 p. 27/32 Weak convergence { } u(t, x) = E g(z(t, t, x)) { ( = E g E(T t)x + T t )} E(T s) dw(s) Kolmogorov equation: u t = u/ t, u x = u/ x, u xx = 2 u/ x ) 2 u t (t, x) (Ax, u x (t, x)) (u Tr xx (t, x)q =, x D(A), t < T, u(t, x) = g(x) ( ) E u(t, X h (T)) = Eg(X h (T)), assume X() = X h () ( ) ( ) ( ) ( ) E u(, X h ()) = E u(, X()) = E g(z(t,, X )) = E g(x(t)) ( Eg(X h (T)) Eg(X(T)) = E u(t, X h (T)) {Itô formula} = E {Kolmogorov eq} = E T T ( ) u(, X h ()) ( ) u t (A h X h, u x ) Tr(u xxq h ) dt ( ) ((A A h )X h, u x ) Tr(u xx(q h Q)) dt
28 p. 28/32 Weak convergence { } u(t, x) = E g(z(t, t, x)) { ( = E g E(T t)x + ( ) u x (t, x) = E E(T t)g (Z(T, t, x)), ( ) u xx (t, x) = E E(T t)g (Z(T, t, x))e(t t) T t )} E(T s) dw(s) Eg(X h (T)) Eg(X(T)) = E T ( ((A Ah )X h (t), u x (t, X h (t)) ) Tr( u xx (t, X h (t))(q h Q) )) dt finite element and semigroup techniques Malliavin calculus??
29 p. 29/32 Ongoing and future work wavelet expansion of noise nonlinear equations dx + AX dt = f(x) dt + g(x) dw adaptive methods based on a posteriori error estimates Cahn-Hilliard-Cook equation dx + A 2 X dt = Af(X) dt + dw Malliavin calculus stochastic wave equation (the semigroup is not analytic) more general noise: H-valued Lévy process
30 p. 3/32 Non-future work We do not study PDE with (smooth) random coefficients u t (a(x, ω) u) = f(x, t, ω)
31 p. 31/32 Implementation Euler s method { U n S h, U = P h u U n U n 1 + ta h U n = P h W n (U n U n 1, χ) + t( U n, χ) = ( W n, χ), χ S h U n (x) = N h k=1 N h k=1 U n k (ϕ k, ϕ j ) + t U n k ϕ k (x), χ = ϕ j, finite element basis functions N h k=1 MU n + tku n = MU n 1 + b n U n k ( ϕ k, ϕ j ) = N h k=1 U n 1 k (ϕ k, ϕ j ) + ( W n, ϕ j )
32 p. 32/32 Implementation How to simulate b n j = ( W n, ϕ j ) = (W(t n ) W(t n 1 ), ϕ j )? Covariance of b n : E(b n i b n j ) = E ( ( W n, ϕ i )( W n, ϕ j ) ) = t(qϕ i, ϕ j ) In other words: E(b n b n ) = tq, Q ij = (Qϕ i, ϕ j ) Cholesky factorization: Q = LL T Take b n = tlβ n, where β n, n = 1, 2,..., are N(, 1) Then E(b n b n ) = E(b n (b n ) T ) = te(lβ n (Lβ n ) T ) = tle(β n (β n ) T )L T = tll T = tq
The semigroup approach to stochastic PDEs and their finite element approximation. Bielefeld IGK, November 2009
The semigroup approach to stochastic PDEs and their finite element approximation Mihály Kovács Department of Mathematics and Statistics University of Otago Dunedin, New Zealand Bielefeld IGK, November
More informationPREPRINT 2009:20. Finite Element Approximation of the Linearized Cahn-Hilliard-Cook Equation STIG LARSSON ALI MESFORUSH
PREPRINT 9: Finite Element Approximation of the Linearized Cahn-Hilliard-Coo Equation STIG LARSSON ALI MESFORUSH Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY
More informationNumerical discretisations of stochastic wave equations
Numerical discretisations of stochastic wave equations David Cohen Matematik och matematisk statistik \ UMIT Research Lab, Umeå universitet Institut für Mathematik, Universität Innsbruck Joint works with
More informationSELECTED TOPICS CLASS FINAL REPORT
SELECTED TOPICS CLASS FINAL REPORT WENJU ZHAO In this report, I have tested the stochastic Navier-Stokes equations with different kind of noise, when i am doing this project, I have encountered several
More informationWEAK CONVERGENCE OF FINITE ELEMENT APPROXIMATIONS OF LINEAR STOCHASTIC EVOLUTION EQUATIONS WITH ADDITIVE NOISE
WEAK CONVERGENCE OF FINITE ELEMENT APPROXIMATIONS OF LINEAR STOCHASTIC EVOLUTION EQUATIONS WITH ADDITIVE NOISE MIHÁLY KOVÁCS, STIG LARSSON,, AND FREDRIK LINDGREN Abstract. A unified approac is given for
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 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 informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationAn Introduction to Stochastic Partial Dierential Equations
An Introduction to Stochastic Partial Dierential Equations Herry Pribawanto Suryawan Dept. of Mathematics, Sanata Dharma University, Yogyakarta 29. August 2014 Herry Pribawanto Suryawan (Math USD) SNAMA
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
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 informationDynamics and stochastic limit for phase transition problems. problems with noise. Dimitra Antonopoulou IACM-FORTH
Dynamics and stochastic limit for phase transition problems with noise IACM-FORTH International Conference on Applied Mathematics, September 16-20, 2013, ACMAC, Heraklion, Greece Table of contents 1 1.1
More informationWHITE NOISE APPROACH TO THE ITÔ FORMULA FOR THE STOCHASTIC HEAT EQUATION
Communications on Stochastic Analysis Vol. 1, No. 2 (27) 311-32 WHITE NOISE APPROACH TO THE ITÔ FORMULA FOR THE STOCHASTIC HEAT EQUATION ALBERTO LANCONELLI Abstract. We derive an Itô s-type formula for
More informationHEAT KERNEL ANALYSIS ON INFINITE-DIMENSIONAL GROUPS. Masha Gordina. University of Connecticut.
HEAT KERNEL ANALYSIS ON INFINITE-DIMENSIONAL GROUPS Masha Gordina University of Connecticut http://www.math.uconn.edu/~gordina 6th Cornell Probability Summer School July 2010 SEGAL-BARGMANN TRANSFORM AND
More informationRough paths methods 4: Application to fbm
Rough paths methods 4: Application to fbm Samy Tindel Purdue University University of Aarhus 2016 Samy T. (Purdue) Rough Paths 4 Aarhus 2016 1 / 67 Outline 1 Main result 2 Construction of the Levy area:
More informationTopics in fractional Brownian motion
Topics in fractional Brownian motion Esko Valkeila Spring School, Jena 25.3. 2011 We plan to discuss the following items during these lectures: Fractional Brownian motion and its properties. Topics in
More informationA NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES
A NOTE ON STOCHASTIC INTEGRALS AS L 2 -CURVES STEFAN TAPPE Abstract. In a work of van Gaans (25a) stochastic integrals are regarded as L 2 -curves. In Filipović and Tappe (28) we have shown the connection
More informationWeak Feller Property and Invariant Measures
Weak Feller Property and Invariant Measures Martin Ondreját, J. Seidler, Z. Brzeźniak Institute of Information Theory and Automation Academy of Sciences Prague September 11, 2012 Outline 2010: Stochastic
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. 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 informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationStochastic Partial Differential Equations with Levy Noise
Stochastic Partial Differential Equations with Levy Noise An Evolution Equation Approach S..PESZAT and J. ZABCZYK Institute of Mathematics, Polish Academy of Sciences' CAMBRIDGE UNIVERSITY PRESS Contents
More informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
More informationSolving Stochastic Partial Differential Equations as Stochastic Differential Equations in Infinite Dimensions - a Review
Solving Stochastic Partial Differential Equations as Stochastic Differential Equations in Infinite Dimensions - a Review L. Gawarecki Kettering University NSF/CBMS Conference Analysis of Stochastic Partial
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationMalliavin calculus and central limit theorems
Malliavin calculus and central limit theorems David Nualart Department of Mathematics Kansas University Seminar on Stochastic Processes 2017 University of Virginia March 8-11 2017 David Nualart (Kansas
More informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationThe stochastic heat equation with a fractional-colored noise: existence of solution
The stochastic heat equation with a fractional-colored noise: existence of solution Raluca Balan (Ottawa) Ciprian Tudor (Paris 1) June 11-12, 27 aluca Balan (Ottawa), Ciprian Tudor (Paris 1) Stochastic
More informationA POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS
A POSTERIORI ERROR ESTIMATES BY RECOVERED GRADIENTS IN PARABOLIC FINITE ELEMENT EQUATIONS D. LEYKEKHMAN AND L. B. WAHLBIN Abstract. This paper considers a posteriori error estimates by averaged gradients
More informationAsymptotic Analysis 88 (2014) DOI /ASY IOS Press
Asymptotic Analysis 88 4) 5 DOI.333/ASY-4 IOS Press Smoluchowski Kramers approximation and large deviations for infinite dimensional gradient systems Sandra Cerrai and Michael Salins Department of Mathematics,
More informationItô formula for the infinite-dimensional fractional Brownian motion
Itô formula for the infinite-dimensional fractional Brownian motion Ciprian A. Tudor SAMOS-MATISSE Université de Panthéon-Sorbonne Paris 1 9, rue de Tolbiac, 75634 Paris France. January 13, 25 Abstract
More informationThe multidimensional Ito Integral and the multidimensional Ito Formula. Eric Mu ller June 1, 2015 Seminar on Stochastic Geometry and its applications
The multidimensional Ito Integral and the multidimensional Ito Formula Eric Mu ller June 1, 215 Seminar on Stochastic Geometry and its applications page 2 Seminar on Stochastic Geometry and its applications
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
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 informationMalliavin Calculus: Analysis on Gaussian spaces
Malliavin Calculus: Analysis on Gaussian spaces Josef Teichmann ETH Zürich Oxford 2011 Isonormal Gaussian process A Gaussian space is a (complete) probability space together with a Hilbert space of centered
More informationStochastic integrals for spde s: a comparison
Stochastic integrals for spde s: a comparison arxiv:11.856v1 [math.pr] 6 Jan 21 Robert C. Dalang Ecole Polytechnique Fédérale de Lausanne Lluís Quer-Sardanyons Universitat Autònoma de Barcelona Abstract
More informationOn the stochastic nonlinear Schrödinger equation
On the stochastic nonlinear Schrödinger equation Annie Millet collaboration with Z. Brzezniak SAMM, Paris 1 and PMA Workshop Women in Applied Mathematics, Heraklion - May 3 211 Outline 1 The NL Shrödinger
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
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 informationNumerical methods for solving stochastic differential equations
Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus
More informationStochastic Integral with respect to Cylindrical Wiener Process
arxiv:math/51151v1 [math.pr] 1 Nov 5 Stochastic Integral with respect to Cylindrical Wiener Process Anna Karczewska Institute of Mathematics, Maria Curie Sk lodowska University pl. M. Curie Sk lodowskiej
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationLévy Processes and Infinitely Divisible Measures in the Dual of afebruary Nuclear2017 Space 1 / 32
Lévy Processes and Infinitely Divisible Measures in the Dual of a Nuclear Space David Applebaum School of Mathematics and Statistics, University of Sheffield, UK Talk at "Workshop on Infinite Dimensional
More informationAnalysis of SPDEs Arising in Path Sampling Part I: The Gaussian Case
Analysis of SPDEs Arising in Path Sampling Part I: The Gaussian Case September 2, 2005 M. Hairer 1, A. M. Stuart 1, J. Voss 1,2, and P. Wiberg 1,3 1 Mathematics Institute, The University of Warwick, Coventry
More informationWhite noise generalization of the Clark-Ocone formula under change of measure
White noise generalization of the Clark-Ocone formula under change of measure Yeliz Yolcu Okur Supervisor: Prof. Bernt Øksendal Co-Advisor: Ass. Prof. Giulia Di Nunno Centre of Mathematics for Applications
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationNonlinear representation, backward SDEs, and application to the Principal-Agent problem
Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationRegularity of Stochastic Partial Differential Equations in Besov Spaces Related to Adaptive Schemes
Regularity of Stochastic Partial Differential Equations in Besov Spaces Related to Adaptive Schemes Petru A. Cioica Philipps-Universität Marburg Workshop on New Discretization Methods for the Numerical
More informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More informationModelling energy forward prices Representation of ambit fields
Modelling energy forward prices Representation of ambit fields Fred Espen Benth Department of Mathematics University of Oslo, Norway In collaboration with Heidar Eyjolfsson, Barbara Rüdiger and Andre Süss
More informationStochastic Differential Equations
CHAPTER 1 Stochastic Differential Equations Consider a stochastic process X t satisfying dx t = bt, X t,w t dt + σt, X t,w t dw t. 1.1 Question. 1 Can we obtain the existence and uniqueness theorem for
More informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More informationInfinite-dimensional calculus under weak spatial regularity of the processes
Infinite-dimensional calculus under weak spatial regularity of the processes Franco Flandoli1, Francesco Russo2, Giovanni Zanco3 November 8, 215 Abstract Two generalizations of Itô formula to infinite-dimensional
More informationContents. 1 Preliminaries 3. Martingales
Table of Preface PART I THE FUNDAMENTAL PRINCIPLES page xv 1 Preliminaries 3 2 Martingales 9 2.1 Martingales and examples 9 2.2 Stopping times 12 2.3 The maximum inequality 13 2.4 Doob s inequality 14
More informationarxiv: v1 [math.pr] 1 May 2014
Submitted to the Brazilian Journal of Probability and Statistics A note on space-time Hölder regularity of mild solutions to stochastic Cauchy problems in L p -spaces arxiv:145.75v1 [math.pr] 1 May 214
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 informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More informationEstimates for the density of functionals of SDE s with irregular drift
Estimates for the density of functionals of SDE s with irregular drift Arturo KOHATSU-HIGA a, Azmi MAKHLOUF a, a Ritsumeikan University and Japan Science and Technology Agency, Japan Abstract We obtain
More informationWeak convergence rates of spectral Galerkin approximations for SPDEs with nonlinear diffusion coefficients
Weak convergence rates of spectral Galerkin approximations for SPDEs with nonlinear diffusion coefficients Daniel Conus, Arnulf Jentzen, and Ryan Kurniawan Lehigh University, USA and EH Zurich, Switzerland
More informationSPDEs, criticality, and renormalisation
SPDEs, criticality, and renormalisation Hendrik Weber Mathematics Institute University of Warwick Potsdam, 06.11.2013 An interesting model from Physics I Ising model Spin configurations: Energy: Inverse
More informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationHölder continuity of the solution to the 3-dimensional stochastic wave equation
Hölder continuity of the solution to the 3-dimensional stochastic wave equation (joint work with Yaozhong Hu and Jingyu Huang) Department of Mathematics Kansas University CBMS Conference: Analysis of Stochastic
More informationON THE STRUCTURE OF GAUSSIAN RANDOM VARIABLES
ON THE STRUCTURE OF GAUSSIAN RANDOM VARIABLES CIPRIAN A. TUDOR We study when a given Gaussian random variable on a given probability space Ω, F,P) is equal almost surely to β 1 where β is a Brownian motion
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationThe Cameron-Martin-Girsanov (CMG) Theorem
The Cameron-Martin-Girsanov (CMG) Theorem There are many versions of the CMG Theorem. In some sense, there are many CMG Theorems. The first version appeared in ] in 944. Here we present a standard version,
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationConvergence analysis of a discontinuous Galerkin method for a sub-diffusion equation
Numer Algor (29 52:69 88 DOI 1.17/s1175-8-9258-8 ORIGINAL PAPER Convergence analysis of a discontinuous Galerkin method for a sub-diffusion equation William McLean Kassem Mustapha Received: 16 September
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationMean-field SDE driven by a fractional BM. A related stochastic control problem
Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,
More informationSTOCHASTIC CONTROLLABILITY OF LINEAR SYSTEMS WITH STATE DELAYS
Int. J. Appl. Math. Comput. Sci., 27, Vol. 17, No. 1, 5 13 DOI: 1.2478/v16-7-1-8 SOCHASIC CONROLLABILIY OF LINEAR SYSEMS WIH SAE DELAYS JERZY KLAMKA Institute of Control Engineering, Silesian University
More informationMoments for the parabolic Anderson model: on a result by Hu and Nualart.
Moments for the parabolic Anderson model: on a result by Hu and Nualart. Daniel Conus November 2 Abstract We consider the Parabolic Anderson Model tu Lu + uẇ, where L is the generator of a Lévy Process
More informationExistence of Invariant Measures of Stochastic Systems with Delay in the Highest Order Partial Derivatives
arxiv:142.2157v1 [math.pr] 1 Feb 214 Existence of Invariant Measures of Stochastic Systems with Delay in the Highest Order Partial Derivatives Kai Liu Department of Mathematical Sciences, School of Physical
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationAn Introduction to Malliavin Calculus. Denis Bell University of North Florida
An Introduction to Malliavin Calculus Denis Bell University of North Florida Motivation - the hypoellipticity problem Definition. A differential operator G is hypoelliptic if, whenever the equation Gu
More informationFINITE-TIME BLOW-UP IN THE ADDITIVE SUPERCRITICAL STOCHASTIC NONLINEAR SCHRÖDINGER EQUATION : THE REAL NOISE CASE
FINITE-TIME BLOW-UP IN THE ADDITIVE SUPERCRITICAL STOCHASTIC NONLINEAR SCHRÖDINGER EQUATION : THE REAL NOISE CASE A. DE BOUARD AND A. DEBUSSCHE Abstract. We review some results concerning the apparition
More informationOptimal series representations of continuous Gaussian random fields
Optimal series representations of continuous Gaussian random fields Antoine AYACHE Université Lille 1 - Laboratoire Paul Painlevé A. Ayache (Lille 1) Optimality of continuous Gaussian series 04/25/2012
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationFrom Fractional Brownian Motion to Multifractional Brownian Motion
From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December
More information(B(t i+1 ) B(t i )) 2
ltcc5.tex Week 5 29 October 213 Ch. V. ITÔ (STOCHASTIC) CALCULUS. WEAK CONVERGENCE. 1. Quadratic Variation. A partition π n of [, t] is a finite set of points t ni such that = t n < t n1
More informationWhat s more chaotic than chaos itself? Brownian Motion - before, after, and beyond.
Include Only If Paper Has a Subtitle Department of Mathematics and Statistics What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond. Math Graduate Seminar March 2, 2011 Outline
More informationStochastic Shear Thickening Fluids: Strong Convergence of the Galerkin Approximation and the Energy Equality 1
Stochastic Shear Thickening Fluids: Strong Convergence of the Galerkin Approximation and the Energy Equality Nobuo Yoshida Contents The stochastic power law fluids. Terminology from hydrodynamics....................................
More informationOn rough PDEs. Samy Tindel. University of Nancy. Rough Paths Analysis and Related Topics - Nagoya 2012
On rough PDEs Samy Tindel University of Nancy Rough Paths Analysis and Related Topics - Nagoya 2012 Joint work with: Aurélien Deya and Massimiliano Gubinelli Samy T. (Nancy) On rough PDEs Nagoya 2012 1
More informationTHE SUBSTITUTION THEOREM FOR SEMILINEAR STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS. Salah-Eldin A. Mohammed and Tusheng Zhang
THE SUBSTITUTION THEOREM FOR SEMILINEAR STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS Salah-Eldin A. Mohammed and Tusheng Zhang Abstract. In this article we establish a substitution theorem for semilinear
More informationContinuous Time Finance
Continuous Time Finance Lisbon 2013 Tomas Björk Stockholm School of Economics Tomas Björk, 2013 Contents Stochastic Calculus (Ch 4-5). Black-Scholes (Ch 6-7. Completeness and hedging (Ch 8-9. The martingale
More informationItô s formula. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Itô s formula Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Itô s formula Probability Theory
More informationBV functions in a Gelfand triple and the stochastic reflection problem on a convex set
BV functions in a Gelfand triple and the stochastic reflection problem on a convex set Xiangchan Zhu Joint work with Prof. Michael Röckner and Rongchan Zhu Xiangchan Zhu ( Joint work with Prof. Michael
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationDiscrete approximation of stochastic integrals with respect to fractional Brownian motion of Hurst index H > 1 2
Discrete approximation of stochastic integrals with respect to fractional Brownian motion of urst index > 1 2 Francesca Biagini 1), Massimo Campanino 2), Serena Fuschini 2) 11th March 28 1) 2) Department
More informationwith deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied
More information