This paper develops exponential moment bounds, regularity, and convergence rates for a class of numerical schemes approximating stochastic convolutions, tailored for SPDEs with non-globally Lipschitz nonlinearities.
Contribution
It introduces taming and truncation techniques in exponential Euler-type schemes, enabling strong convergence analysis for challenging SPDEs.
Findings
01
Established exponential moment bounds for the schemes
02
Proved uniform Hölder continuity in time
03
Demonstrated strong convergence rates
Abstract
In this article we establish exponential moment bounds, moment bounds in fractional order smoothness spaces, a uniform H\"older continuity in time, and strong convergence rates for a class of fully discrete exponential Euler-type numerical approximations of infinite dimensional stochastic convolution processes. The considered approximations involve specific taming and truncation terms and are therefore well suited to be used in the context of SPDEs with non-globally Lipschitz continuous nonlinearities.
|\theta|_{T}=\max\!\Big{\{}x\in(0,\infty)\colon\big{(}\exists\,a,b\in\theta\colon\big{[}x=b-a\text{ and }\theta\cap(a,b)=\emptyset\big{]}\big{)}\Big{\}}\in(0,T],
|\theta|_{T}=\max\!\Big{\{}x\in(0,\infty)\colon\big{(}\exists\,a,b\in\theta\colon\big{[}x=b-a\text{ and }\theta\cap(a,b)=\emptyset\big{]}\big{)}\Big{\}}\in(0,T],
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Exponential moment bounds and strong
convergence
rates for
tamed-truncated
numerical
approximations of stochastic convolutions
Arnulf Jentzen1, Felix Lindner2, and Primož Pušnik3
1 Seminar for Applied Mathematics, Department of Mathematics,
In this article we establish
exponential moment bounds,
moment bounds in
fractional order smoothness spaces,
a uniform Hölder continuity in time,
and
strong convergence rates
for a class of fully discrete exponential Euler-type numerical approximations of infinite dimensional stochastic convolution processes.
The considered approximations involve specific taming and truncation terms and are therefore well suited to be used in the context of SPDEs with
non-globally
Lipschitz continuous
nonlinearities.
Stochastic partial differential equations
(SPDEs) of evolutionary type
are important modeling tools
in economics and the natural sciences
(see, e.g., Harms et al. [12, Theorem 3.5],
Filipović et al. [10, Equation (1.2)],
Blömker & Romito [5, Equation (1)],
Hairer [11, Equation (3)],
Mourrat & Weber [21, Equation (1.1)],
Birnir [3, Equation (7)],
and
Birnir [4, Equation (1.5)]).
However, exact solutions to SPDEs are
usually not known explicitly.
Therefore, it has been and still is a very active research area to develop and analyze numerical approximation methods which approximate the exact solutions of SPDEs with a reasonable approximation accuracy in a reasonable computational time.
It is known that in order to
approximate the exact solutions of SPDEs
appropriately, the
numerical methods employed
should enjoy similar statistical properties,
such as finite uniform moment bounds
(see, e.g.,
Hutzenthaler & Jentzen [14] and the references therein).
Unfortunately,
moments of the easily realizable Euler-Maruyama
and exponential Euler approximation methods
are known to diverge
for some
stochastic differential equations (SDEs) and SPDEs with superlinearly growing nonlinearties
(see, e.g., Hutzenthaler et al. [15, 17]).
This poses the challenge to develop new efficient approximation
methods which preserve finite moments
(see, e.g.,
Hutzenthaler et al. [16, Theorem 1.1 and Lemma 3.9],
Gan & Wang [26, Theorem 3.2 and Lemma 3.4],
Hutzenthaler & Jentzen [14, Corollary 2.21 and Theorem 3.15],
Tretyakov & Zhang [25, Theorem 2.1],
Sabanis [22, Theorem 2.2, Corollary 2.3, and Lemmas 3.2, 3.3],
and
Sabanis [23, Theorems 1–3, Lemmas 1,2]
for finite dimensional
stochastic evolution equations
and, e.g.,
[19, Proposition 7.3, Theorem 7.6],
Becker & Jentzen [2, Corollaries 5.1, 6.15, and 6.17],
Becker & Jentzen [1, Lemma 5.4, Theorem 1.1]
for infinite dimensional SPDEs).
In this context, it has been revealed recently in, e.g.,
Hutzenthaler & Jentzen [13, Theorem 1.3]
(cf., e.g.,
Dörsek [9, Proposition 3.1],
Hutzenthaler et al. [18, Corollary 2.10],
and [20, Corollary 3.4])
that finite exponential moments of numerical schemes are crucial for deriving strong convergence with rates in the case of SDEs and SPDEs
with
non-globally Lipschitz continuous
and
non-globally monotone
nonlinearities.
In this article we derive finite uniform exponential moment bounds for a class of fully discrete exponential Euler-type numerical approximations of infinite dimensional stochastic convolution processes
(see Corollary 3.4
in Section 3 below). The considered numerical approximations involve specific taming and truncation terms and are therefore well suited to be applied in the context of semi-linear SPDEs with non-globally Lipschitz continuous and non-globally monotone nonlinearities. In addition to deriving exponential moment bounds we also establish finite uniform moment bounds in fractional order smoothness spaces (see Corollary 3.1
in Section 3 below), a uniform Hölder continuity in time
(see Corollary 3.2 in Section 3 below) as well as strong convergences rates for the considered numerical approximations
(see Corollary 3.3 in Section 3 below).
The application of our results to semi-linear SPDEs such as stochastic Burgers equation will be the subject of a future research article.
In
Theorem 1.1
below
we
illustrate the
results established
in
Corollary 3.3
and
Corollary 3.4.
The stochastic convolution process and its numerical approximations are denoted by
O:[0,T]×Ω→D((−A)γ)
and
OM,N:[0,T]×Ω→PN(H),
M,N∈\mathbbmN, respectively.
Theorem 1.1**.**
Let
(H,⟨⋅,⋅⟩H,∥⋅∥H)
and
(U,⟨⋅,⋅⟩U,∥⋅∥U)
be
separable
\mathbbmR-Hilbert spaces,
let
(hn)n∈\mathbbmN⊆H
be an
orthonormal basis of H,
let
(vn)n∈\mathbbmN⊆\mathbbmR
satisfy
supn∈\mathbbmNvn<0,
let
A:D(A)⊆H→H
be the linear operator
which satisfies
D(A)={v∈H:∑n=1∞∣vn⟨hn,v⟩H∣2<∞}
and
∀v∈D(A):Av=∑n=1∞vn⟨hn,v⟩Hhn,
let
p,T∈(0,∞),
β∈[0,∞),
γ∈[0,\nicefrac12+β),
η∈[0,\nicefrac12+β−γ),
ρ∈[0,\nicefrac12+β−γ)∩[0,\nicefrac12),
B∈HS(U,D((−A)β)),
let
(Ω,F,P,(Ft)t∈[0,T])
be a filtered probability space which fulfills the usual conditions,
let (Wt)t∈[0,T]
be an IdU-cylindrical
(Ft)t∈[0,T]-Wiener process,
let
O:[0,T]×Ω→D((−A)γ)
be a stochastic process which satisfies for every
t∈[0,T] that
P(Ot=∫0te(t−s)ABdWs)=1,
let
(PN)N∈\mathbbmN⊆L(H)
satisfy
for every N∈\mathbbmN,
x∈H
that
PN(x)=∑n=1N⟨hn,x⟩Hhn,
let
WN:[0,T]×Ω→PN(H),
N∈\mathbbmN,
be
stochastic processes
which satisfy for every
N∈\mathbbmN,
t∈[0,T]
that
P(WtN=∫0tPNBdWs)=1,
let
χM,N:[0,T]×Ω→[0,1],
M,N∈\mathbbmN,
be
(Ft)t∈[0,T]-adapted stochastic processes
which satisfy
supM,N∈\mathbbmNsups∈[0,T](E[∣χsM,N−1∣max{p,2}]Mmax{p,2}ρ)<∞,
and
let
OM,N:[0,T]×Ω→PN(H),
M,N∈\mathbbmN,
be
stochastic processes
which satisfy for every
M,N∈\mathbbmN,
m∈{0,1,…,M−1},
t∈[\nicefracmTM,\nicefrac(m+1)TM]
that
O0M,N=0
and
[TABLE]
Then
(i)
it holds that
\inf_{\varepsilon\in(0,\infty)}\sup_{M,N\in{\mathbbm{N}}}\sup_{t\in[0,T]}{\mathbb{E}}\big{[}\exp\!\big{(}\varepsilon\|{\bf O}_{t}^{M,N}\|_{H}^{2}\big{)}\big{]}<\infty and
2. (ii)
there exists a real number
C∈\mathbbmR
such that for every
M,N∈\mathbbmN
it holds that
[TABLE]
Observe that
item (i)
in Theorem 1.1
is a direct consequence of
Corollary 3.4
(with
H=H,
U=U,
H=(hn)n∈\mathbbmN,
v=v,
A=A,
β=β,
T=T,
(Ω,F,P)=(Ω,F,P),
(Ft)t∈[0,T]=(Ft)t∈[0,T],
(Wt)t∈[0,T]=(Wt)t∈[0,T],
B=B,
P{h1,…,hN}=PN,
P^U=IdU,
χ{0,T/M,…,T},{h1,…,hN},U=χM,N,
O{0,T/M,…,T},{h1,…,hN},U=OM,N,
ε=ε
for
M,N∈\mathbbmN,
ε∈[0,\nicefrac1(8[max{∥B∥HS(U,H),1}]2max{T,1})2)
in the notation of
Corollary 3.4)
and Hölder’s inequality.
Moreover, note that
item (ii)
in Theorem 1.1
follows from
Corollary 3.3
(with
H=H,
U=U,
H=(hn)n∈\mathbbmN,
v=v,
A=A,
β=β,
T=T,
(Ω,F,P)=(Ω,F,P),
(Ft)t∈[0,T]=(Ft)t∈[0,T],
(Wt)t∈[0,T]=(Wt)t∈[0,T],
B=B,
P{h1,…,hN}=PN,
P^U=IdU,
χ{0,T/M,…,T},{h1,…,hN},U=χM,N,
O{0,T/M,…,T},{h1,…,hN},U=OM,N,
p=max{p,1},
C=C,
γ=γ,
η=η,
ρ=ρ,
O=O
for
M,N∈\mathbbmN
in the notation of
Corollary 3.3).
We would like to point out
that the exponential moment bound in
Corollary 3.4 below
(see also item (i) above)
is not a direct consequence of the
one
established
in [20].
The difference between the numerical
method in [20, (1) in Section 1]
and (1) above lies, roughly speaking, in the
less restrictive choice for the truncation functions
χM,N:[0,T]×Ω→[0,1],
M,N∈\mathbbmN,
in (1)
compared to [20, (1) in Section 1].
This extended class of discrete approximations allows
to truncate the numerical method
independently of the current value of the approximation process itself.
The flexibility in the choice of the truncation function
is, in turn, important
for applying
Corollaries 3.1–3.4
to establish
strong convergence rates
for fully discrete numerical approximations
in the case of SPDEs
with non-globally Lipschitz continuous
and non-globally monotone nonlinearities.
The remainder of this article is structured as follows.
In Section 2 we analyze a temporally semi-discrete version of our approximation scheme
(see (4) in
Setting 2.1
below).
In Subsection 2.1
the considered numerical approximations are rewritten as Itô processes and finite moment
bounds in fractional order smoothness spaces are derived.
The Itô representation enables us to
establish
Hölder regularity in time in Subsection 2.2.
Furthermore, under additional assumptions on the truncation
function we establish
temporal strong convergence rates of the proposed numerical methods
in Subsection 2.3.
In Subsection 2.4
we further improve our moment estimates from
Subsection 2.1
in order to derive finite exponential moment bounds
in Lemma 2.8.
The results from Section 2
are combined in Section 3 to establish
uniform moment
bounds in fractional order smoothness spaces,
a uniform Hölder regularity in time,
strong convergence rates,
and uniform exponential moment bounds
for fully discrete
tamed-truncated numerical approximations
in Corollaries 3.1–3.4,
respectively.
1.1 General setting
Throughout this article the following setting is frequently used.
Setting 1.1**.**
For every set X let P(X) be the power set of X,
for every set X let
P0(X) be the set given by
P0(X)={θ∈P(X):θhas finitely many elements},
for every
T∈(0,∞)
let ϖT
be the set given by
ϖT={θ∈P0([0,T]):{0,T}⊆θ},
for every T∈(0,∞)
let
∣⋅∣T:ϖT→[0,T]
be the function which satisfies for every
θ∈ϖT that
[TABLE]
for every
θ∈(∪T∈(0,∞)ϖT)
let
└⋅┘θ:[0,∞)→[0,∞)
be the function which satisfies for every
t∈(0,∞)
that
└t┘θ=max([0,t)∩θ)
and
└0┘θ=0,
for every measure space
(Ω,F,μ),
every measurable space (S,S),
every set R,
and every function
f:Ω→R
let
[f]μ,S
be the set given by
[f]μ,S={g:Ω→S:(∃A∈F:[μ(A)=0 and {ω∈Ω:f(ω)=g(ω)}⊆A]) and (∀A∈S:g−1(A)∈F)},
let
(H,⟨⋅,⋅⟩H,∥⋅∥H)
and
(U,⟨⋅,⋅⟩U,∥⋅∥U)
be separable \mathbbmR-Hilbert spaces,
let
H⊆H
be a non-empty orthonormal basis of H,
let
v:H→\mathbbmR
be a function which satisfies
suph∈Hvh<0,
let
A:D(A)⊆H→H
be the linear operator which satisfies
D(A)={v∈H:∑h∈H∣vh⟨h,v⟩H∣2<∞}
and
∀v∈D(A):Av=∑h∈Hvh⟨h,v⟩Hh,
and
let
(Hr,⟨⋅,⋅⟩Hr,∥⋅∥Hr), r∈\mathbbmR, be a family of interpolation spaces associated to −A
(cf., e.g., [24, Section 3.7]).
2 Regularity properties of
temporally
semi-discrete
tamed-truncated
approximations
of stochastic convolutions
Setting 2.1**.**
Assume Setting 1.1,
let
β∈[0,∞),
γ∈[0,\nicefrac12+β),
T∈(0,∞),
θ∈ϖT,
B∈HS(U,Hβ),
let
B∈L(H,U)
be the bounded linear operator
which satisfies for every
u∈U, h∈H that
⟨Bu,h⟩H=⟨u,Bh⟩U,
let
(Ω,F,P,(Ft)t∈[0,T])
be a filtered probability space which fulfills the usual conditions,
let (Wt)t∈[0,T]
be an IdU-cylindrical
(Ft)t∈[0,T]-Wiener process,
let
χ:[0,T]×Ω→[0,1]
be an (Ft)t∈[0,T]-adapted
stochastic process,
and
let
O:[0,T]×Ω→Hγ
be a
stochastic process
which satisfies for every t∈[0,T] that
O0=0
and
[TABLE]
2.1 Moment bounds
for temporally semi-discrete approximations of stochastic convolutions
In this subsection we
provide
in
Lemma 2.1
a representation
of
the approximation process
O:[0,T]×Ω→Hγ
from Setting 2.1
as a mild Itô process
(cf. Da Prato et al. [7, Definition 1]).
This enables us
to obtain certain
moment bounds
for
O:[0,T]×Ω→Hγ
in
Lemmas 2.2
and 2.3.
Lemma 2.1**.**
Assume Setting 2.1
and for every s∈[0,T] let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be an (Ft)t∈[s,T]-adapted stochastic process
with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Then it holds for every t∈[0,T] that
Throughout this proof
let
ψ:H→H be the function
which satisfies for every v∈H that
ψ(v)=1+∥v∥H2v
and let
U⊆U be an orthonormal basis of U.
Note that for every
u,v,z∈H it holds that
[TABLE]
and
[TABLE]
Itô’s formula
(see, e.g.,
Brzeźniak et al. [6, Theorem 2.4]
with
H=U,
E=H,
F=H,
f=([0,T]×H∋(t,x)↦ψ(x)∈H),
Φ=([0,T]×Ω∋(t,ω)↦B∈HS(U,H))
in the notation of
Brzeźniak et al. [6, Theorem 2.4])
hence proves
that
for every
s∈[0,T], t∈[s,T]
it holds that
Assume Setting 2.1,
let
p∈[2,∞),
ρ∈[β,\nicefrac12+β),
η∈[β,1+β),
and
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be an (Ft)t∈[s,T]-adapted
stochastic process
with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Then
it holds
for every
s∈[0,T], t∈[s,T]
that
Note that
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
shows for every s∈[0,T], t∈[s,T]
that
[TABLE]
In addition, observe for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Combining this and (LABEL:eq:some111)
establishes (LABEL:eq:Estimate1).
Furthermore, note that
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
proves that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
The Hölder inequality, the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8],
and
the fact that
∥B∥HS(U,H)≤∣suph∈Hvh∣−β∥B∥HS(U,Hβ)
therefore
ensure that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
Combining this and (LABEL:eq:some112)
assures that (LABEL:eq:Estimate2) holds.
In the next step we observe that
for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
This and the fact that
∥B∥L(Hβ,U)≤∣suph∈Hvh∣−β∥B∥L(H,U)
demonstrate that
for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
The
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
and
the fact that
∥B∥HS(U,H)≤∣suph∈Hvh∣−β∥B∥HS(U,Hβ)
hence
establish for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Moreover, note that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
Combining this and (LABEL:eq:intermediate)
establishes (LABEL:eq:Estimate3).
The proof of Lemma 2.2
is thus completed.
∎
Lemma 2.3**.**
Assume Setting 2.1,
assume that
β≤γ,
and
let
p∈[2,∞),
t∈[0,T].
Then it holds that
Throughout this proof
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,u(ω))(u,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be
an
(Fu)u∈[s,T]-adapted
stochastic process with continuous sample paths
which satisfies for every
u∈[s,T]
that
[Xs,u]P,B(Hβ)=∫suBdWτ.
Lemma 2.1
(with
Xs,u=Xs,u
for
u∈[s,T], s∈[0,T]
in the notation of
Lemma 2.1)
implies that
[TABLE]
Moreover, observe that the
triangle inequality
proves that
[TABLE]
Next note that
Lemma 2.2
(with
p=p,
ρ=γ,
η=γ,
Xs,τ=Xs,τ
for
τ∈[s,T],
s∈[0,T]
in the notation of Lemma 2.2)
shows that
[TABLE]
[TABLE]
and
[TABLE]
Next we combine (LABEL:eq:split_triangle)–(LABEL:eq:helping_lemma)
to obtain that
[TABLE]
This, (LABEL:eq:triangle_inequality),
(LABEL:eq:sec_est), the fact that
2+21≤3,
the
fact that
1+2β−2γ≤2(1+β−γ),
and the fact that
∀x∈(0,1]:\nicefrac1x≤\nicefrac1x
demonstrate that
of temporally semi-discrete approximations of stochastic convolutions
In this subsection we
combine
Lemma 2.1
and
Lemma 2.2
to establish in
Lemma 2.4
a temporal regularity property
for the approximation
process
O:[0,T]×Ω→Hγ
from Setting 2.1.
Lemma 2.4**.**
Assume Setting 2.1,
assume that
β≤γ,
and let
p∈[2,∞),
ρ∈[0,\nicefrac12+β−γ).
Then it holds for every
s∈[0,T], t∈[s,T]
that
Throughout this proof
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be
an
(Ft)t∈[s,T]-adapted
stochastic process with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Lemma 2.1
(with
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T]
in the notation of
Lemma 2.1)
and the triangle inequality
prove for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Furthermore, observe that the triangle inequality implies for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Next note that
Lemma 2.2
(with
p=p,
ρ=γ,
η=γ,
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T]
in the notation of Lemma 2.2)
shows for every
s∈[0,T], t∈[s,T]
that
[TABLE]
[TABLE]
and
[TABLE]
Moreover,
observe that
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
assures
for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Furthermore, note that the fact that
γ+ρ−β∈[0,\nicefrac12)
and the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
hence imply for every
s∈[0,T], t∈[s,T]
that
[TABLE]
In addition, observe that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
Combining this and (LABEL:eq:BDGEstimated)
demonstrates that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
In the next step we observe that for every
s∈[0,T], t∈[s,T]
it holds that
[TABLE]
This, the fact that
∥B∥L(Hβ,U)≤∣suph∈Hvh∣−β∥B∥L(H,U),
and the fact that
∥B∥HS(U,H)≤∣suph∈Hvh∣−β∥B∥HS(U,Hβ)
show for every
s∈[0,T], t∈[s,T]
that
[TABLE]
The
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
therefore proves
for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Moreover, note that for every
s∈[0,T]
it holds that
[TABLE]
Estimate (LABEL:eq:Needed2)
hence establishes
for every
s∈[0,T], t∈[s,T]
that
[TABLE]
Combining (LABEL:eq:first_diff)–(LABEL:eq:est3),
the fact that
∥B∥HS(U,H)≤∣suph∈Hvh∣−β∥B∥HS(U,Hβ),
(LABEL:eq:Diff1),
(LABEL:eq:Diff2),
the fact that
1+2β−2γ≤1+β−γ,
the fact that
1+2(β−γ−ρ)≤1+β−γ−ρ,
and the fact that
21+2p∥B∥HS(U,Hβ)+22∥B∥HS(U,Hβ)2≤2(1+∥B∥HS(U,Hβ)2)p
implies that for every
s∈[0,T], t∈[s,T]
it holds that
2.3 Error estimates for temporally semi-discrete approximations of stochastic convolutions
In this subsection we
combine
Lemma 2.1
and
Lemma 2.2
above to establish in Lemma 2.5 below for every p∈[2,∞) an upper bound for the strong Lp-distance between the approximation process
O:[0,T]×Ω→Hγ
from
Setting 2.1
and a suitable stochastic convolution process related to
O:[0,T]×Ω→Hγ.
Lemma 2.5**.**
Assume Setting 2.1,
let
C∈[1,∞),
B~∈HS(U,Hβ),
p∈[2,∞),
η∈[0,\nicefrac12+β−γ),
ρ∈[0,\nicefrac12+β−γ)∩[0,\nicefrac12),
assume for every
s∈[0,T]
that
∥χ└s┘θ−1∥Lp(P;\mathbbmR)≤C[∣θ∣T]ρ,
and
let
O:[0,T]×Ω→Hγ
be a stochastic process which satisfies for every
t∈[0,T] that
[Ot]P,B(Hγ)=∫0te(t−s)AB~dWs.
Then
Throughout this proof
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be an
(Ft)t∈[s,T]-adapted
stochastic process with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Observe that
Lemma 2.1
(with
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T]
in the notation of
Lemma 2.1)
and the triangle inequality prove
for every t∈[0,T] that
[TABLE]
In the next step we
note that
the fact that
η<\nicefrac12+β−γ
ensures that
max{γ+η−β,0}<\nicefrac12.
The
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
hence shows that
for every t∈[0,T] it holds that
[TABLE]
Furthermore, observe that the triangle inequality
implies for every t∈[0,T] that
[TABLE]
In addition, note that the triangle inequality
assures for every
t∈[0,T]
that
[TABLE]
Moreover, observe that
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
demonstrates that for every
t∈[0,T]
it holds that
[TABLE]
The
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
and the fact that
∥B∥HS(U,H)≤∣suph∈Hvh∣−β∥B∥HS(U,Hβ)
hence ensure that for every t∈[0,T]
it holds that
[TABLE]
Furthermore, observe that
the fact that
1−2max{0,γ−β}=max{1,1−2γ+2β}>0
ensures that
for every
t∈[0,T] it holds that
[TABLE]
Combining this and (LABEL:eq:estimate_first0)
implies that for every
t∈[0,T] it holds that
[TABLE]
In the next step we note that
the fact that
1−2ρ−2max{0,γ−β}=1−2ρ+2min{0,β−γ}>0
and
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
show that for every
t∈[0,T] it holds that
[TABLE]
Moreover, observe that the assumption that
∀u∈[0,T],θ∈ϖT:∥χ└u┘θ−1∥Lp(P;\mathbbmR)≤C[∣θ∣T]ρ
and
the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
establish for every
t∈[0,T]
that
[TABLE]
Combining this,
(LABEL:eq:triangle+burkholder_estimate_error_stochastic),
(LABEL:eq:estimate_first),
(LABEL:eq:estimate_second),
and Hölder’s inequality
demonstrates that for every
t∈[0,T]
it holds that
[TABLE]
Furthermore, note that
the fact that
γ=max{β,γ}+min{0,γ−β}
and
Lemma 2.2
(with
p=p,
ρ=max{β,γ},
η=max{β,γ},
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T],
h∈(0,T]
in the notation of
Lemma 2.2)
prove that for every t∈[0,T] it holds that
[TABLE]
and
[TABLE]
Moreover, observe that
[TABLE]
Combining (LABEL:eq:basic_triangle),
(LABEL:eq:B_defect),
(LABEL:eq:basic_triangle2),
(LABEL:eq:first_part_estimate),
(LABEL:eq:apply_useful_estimate1),
(LABEL:eq:apply_useful_estimate2),
the fact that
1+2β−2max{β,γ}≤1+β−max{β,γ},
and the fact that
23+22+22≤8
hence ensures that
for every t∈[0,T] it holds that
2.4 Exponential moments of temporally semi-discrete approximations of stochastic convolutions
In this subsection
we
first derive two auxiliary lemmas
(see
Lemma 2.6
and
Lemma 2.7
below)
which we then
combine to establish
in Lemma 2.8 below
appropriate
exponential moment bounds
for the
approximation
process
O:[0,T]×Ω→Hγ
from
Setting 2.1.
Lemma 2.6**.**
Assume Setting 2.1,
let
p∈[1,∞),
and
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be an (Ft)t∈[s,T]-adapted
stochastic process
with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Then
it holds
for every
t∈[0,T]
that
Assume Setting 2.1
and
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be an (Ft)t∈[s,T]-adapted
stochastic process
with continuous sample paths
which satisfies for every
t∈[s,T]
that
[Xs,t]P,B(Hβ)=∫stBdWu.
Then
it holds
for every
n∈\mathbbmN,
t∈[0,T]
that
Throughout this proof
let U⊆U be an orthonormal basis of U,
let
Z:[0,T]×Ω→H
be an
(Ft)t∈[0,T]-adapted
stochastic process
which satisfies for every
t∈[0,T] that
[TABLE]
and let
Z:[0,T]×Ω→H
be an
(Ft)t∈[0,T]-adapted
stochastic process
with left-continuous sample paths
and finite right limits
which satisfies for every
t∈[0,T] that
[TABLE]
Note that
Itô’s formula
proves for every
p∈[2,∞),
t∈[0,T] that
[TABLE]
In addition, observe that the
Burkholder-Davis-Gundy type inequality in Lemma 7.7 in Da Prato & Zabczyk [8]
ensures
for every
p∈[2,∞),
t∈[0,T] that
[TABLE]
This implies for every
p∈[2,∞),
t∈[0,T]
that
[TABLE]
Moreover, note that for every
t∈[0,T] it holds that
[TABLE]
Combining this,
(LABEL:eq:ApplyIto),
(LABEL:eq:EstablishMartingale1),
and Tonelli’s theorem
establishes that for every
p∈[2,∞),
t∈[0,T] it holds that
[TABLE]
The Cauchy-Schwarz inequality hence assures for every
p∈[2,∞),
t∈[0,T]
that
[TABLE]
This
demonstrates that for every
n∈\mathbbmN,
t∈[0,T] it holds that
[TABLE]
Therefore, we obtain that for every
n∈\mathbbmN, t0∈[0,T] it holds that
Throughout this proof
for every s∈[0,T]
let
Xs,(⋅)(⋅)=(Xs,t(ω))(t,ω)∈[s,T]×Ω:[s,T]×Ω→Hβ
be
an
(Fu)u∈[s,T]-adapted
stochastic process with continuous sample paths
which satisfies for every
u∈[s,T]
that
[Xs,u]P,B(Hβ)=∫suBdWτ.
Lemma 2.1
(with
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T]
in the notation of
Lemma 2.1),
Lemma 2.6
(with
p=2n,
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T],
n∈\mathbbmN
in the notation of
Lemma 2.6),
Lemma 2.7
(with
Xs,t=Xs,t
for
t∈[s,T],
s∈[0,T],
n∈\mathbbmN
in the notation of
Lemma 2.7),
and the triangle inequality
ensure that for every
n∈\mathbbmN,
t∈[0,T]
it holds that
[TABLE]
This, the fact that for every
x∈[0,∞) it holds that
ex≤2[∑m=0∞(2m)!x2m]
(see, e.g., Lemma 2.4 in Hutzenthaler et al. [18]),
the dominated convergence theorem,
and the fact that
∀m∈\mathbbmN:(4m)!≤24m[(2m)!]2
imply that for every
t∈[0,T] it holds that
tamed-truncated space-time
approximations of stochastic convolutions
Setting 3.1**.**
Assume Setting 1.1,
let
β∈[0,∞),
T∈(0,∞),
let
(Ω,F,P,(Ft)t∈[0,T])
be a filtered probability space which fulfills the usual conditions,
let (Wt)t∈[0,T]
be an IdU-cylindrical
(Ft)t∈[0,T]-Wiener process,
let
B∈HS(U,Hβ),
let
U⊆U
be an orthonormal basis of U,
let
PI:H→H,
I∈P(H),
and
P^J:U→U,
J∈P(U),
be the linear operators
which satisfy for every
x∈H,
y∈U,
I∈P(H),
J∈P(U)
that
PI(x)=∑h∈I⟨h,x⟩Hh
and
P^J(y)=∑u∈J⟨u,y⟩Uu,
let
χθ,I,J:[0,T]×Ω→[0,1],
θ∈ϖT,
I∈P0(H),
J∈P(U),
be
(Ft)t∈[0,T]-adapted
stochastic
processes,
and
let
Oθ,I,J:[0,T]×Ω→PI(H),
θ∈ϖT,
I∈P0(H),
J∈P(U),
be
stochastic
processes
which satisfy for every
θ∈ϖT,
I∈P0(H),
J∈P(U),
t∈[0,T]
that
O0θ,I,J=0
and
[TABLE]
3.1 Main results
Here we apply the results from Section 2 in order to obtain our main results concerning
tamed-truncated
space-time
approximations (see (77) above) of stochastic convolutions.
A uniform boundedness of moments in fractional order smoothness spaces is presented in Corollary 3.1, a uniform Hölder continuity in time is shown in Corollary 3.2, strong convergence rates are established in Corollary 3.3, and Corollary 3.4 concerns a uniform boundedness of exponential moments.
Corollary 3.1**.**
Assume Setting 3.1
and
let
p∈[1,∞),
γ∈[0,\nicefrac12+β).
Then it holds that
Observe that
Lemma 2.3
(with
β=β,
γ=max{γ,β},
T=T,
θ=θ,
(Ω,F,P)=(Ω,F,P),
(Fu)u∈[0,T]=(Fu)u∈[0,T],
(Wu)u∈[0,T]=(Wu)u∈[0,T],
B=(U∋u↦PIBP^J(u)∈Hβ),
χ=χθ,I,J,
O=([0,T]×Ω∋(t,ω)↦Otθ,I,J(ω)∈Hmax{γ,β}),
p=max{p,2}
for
θ∈ϖT,
I∈P0(H),
J∈P(U)
in the notation of Lemma 2.3)
and Hölder’s inequality
show that (78) holds.
The proof of Corollary 3.1
is thus completed.
∎
Corollary 3.2**.**
Assume Setting 3.1
and
let
p∈[1,∞),
γ∈[0,\nicefrac12+β),
ρ∈[0,\nicefrac12+β−γ)∩[0,\nicefrac12).
Then it holds for every
s∈[0,T], t∈[s,T]
that
Note that
Lemma 2.4
(with
H=H,
β=β,
γ=max{γ,β},
T=T,
θ=θ,
(Ω,F,P)=(Ω,F,P),
(Ft)t∈[0,T]=(Ft)t∈[0,T],
(Wt)t∈[0,T]=(Wt)t∈[0,T],
B=(U∋u↦PIBP^J(u)∈Hβ),
χ=χθ,I,J,
O=([0,T]×Ω∋(t,ω)↦Otθ,I,J(ω)∈Hmax{γ,β}),
p=max{p,2},
ρ=ρ
for
θ∈ϖT,
I∈P0(H),
J∈P(U)
in the notation of Lemma 2.4)
establishes (LABEL:eq:DownToGridEstimate).
The proof of
Corollary 3.2
is thus completed.
∎
Corollary 3.3**.**
Assume Setting 3.1,
let
p,C∈[1,∞),
γ∈[0,\nicefrac12+β),
η∈[0,\nicefrac12+β−γ),
ρ∈[0,\nicefrac12+β−γ)∩[0,\nicefrac12),
assume for every
s∈[0,T],
θ∈ϖT,
I∈P0(H),
J∈P(U)
that
∥χ└s┘θθ,I,J−1∥Lmax{p,2}(P;\mathbbmR)≤C[∣θ∣T]ρ,
and
let O:[0,T]×Ω→Hγ
be a stochastic process
which satisfies for every t∈[0,T] that
[Ot]P,B(Hγ)=∫0te(t−s)ABdWs.
Then it holds
for every
I∈P0(H),
J∈P(U),
K∈P(H),
θ∈ϖT
with I⊆K
that
Note that
Lemma 2.5
(with
H=H,
β=β,
γ=γ,
T=T,
(Ω,F,P)=(Ω,F,P),
(Ft)t∈[0,T]=(Ft)t∈[0,T],
(Wt)t∈[0,T]=(Wt)t∈[0,T],
B=(U∋u↦PIBP^J(u)∈Hβ),
χ=χθ,I,J,
O=([0,T]×Ω∋(t,ω)↦Otθ,I,J(ω)∈Hγ),
C=C,
B~=(U∋u↦PKB(u)∈Hβ),
p=max{p,2},
η=η,
ρ=ρ,
O=([0,T]×Ω∋(t,ω)↦PKOt(ω)∈Hγ)
for
θ∈ϖT,
I∈P0(H),
J∈P(U),
K∈P(H)
with
I⊆K
in the notation of Lemma 2.5)
ensures
that for every
θ∈ϖT,
I∈P0(H),
J∈P(U),
K∈P(H)
with
I⊆K
it holds that
Note that
Lemma 2.8
(with
H=H,
β=β,
γ=0,
T=T,
θ=θ,
(Ω,F,P)=(Ω,F,P),
(Ft)t∈[0,T]=(Ft)t∈[0,T],
B=(U∋x↦PIBP^J(x)∈Hβ),
χ=χθ,I,J,
O=([0,T]×Ω∋(t,ω)↦Otθ,I,J(ω)∈H),
ε=ε
for
I∈P0(H),
J∈P(U),
θ∈ϖT
in the notation of
Lemma 2.8)
assures that (LABEL:eq:Exp2) holds.
The proof of Corollary 3.4
is thus completed.
∎
Acknowledgements
This project has been partially supported through the SNSF-Research project 200021_156603 ”Numerical
approximations of nonlinear stochastic ordinary and partial differential equations”.
Bibliography26
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] Becker, S., Gess, B., Jentzen, A., and Kloeden, P. E. Strong convergence rates for explicit space-time discrete numerical approximations of stochastic Allen-Cahn equations. ar Xiv:1711.02423 (2017), 104 pages.
2[2] Becker, S., and Jentzen, A. Strong convergence rates for nonlinearity-truncated Euler-type approximations of stochastic Ginzburg-Landau equations. ar Xiv:1601.05756 (2016), 58 pages.
3[3] Birnir, B. The Kolmogorov-Obukhov statistical theory of turbulence. J. Nonlinear Sci. 23 , 4 (2013), 657–688.
4[4] Birnir, B. The Kolmogorov-Obukhov theory of turbulence . Springer Briefs in Mathematics. Springer, New York, 2013. A mathematical theory of turbulence.
5[5] Blömker, D., and Romito, M. Stochastic PD Es and lack of regularity: a surface growth equation with noise: existence, uniqueness, and blow-up. Jahresber. Dtsch. Math.-Ver. 117 , 4 (2015), 233–286.
6[6] Brzeźniak, Z., van Neerven, J. M. A. M., Veraar, M. C., and Weis, L. Itô’s formula in UMD Banach spaces and regularity of solutions of the Zakai equation. J. Differential Equations 245 , 1 (2008), 30–58.
7[7] Da Prato, G., Jentzen, A., and Röckner, M. A mild Itô formula for SPD Es. ar Xiv:1009.3526 (2012), 39 pages.
8[8] Da Prato, G., and Zabczyk, J. Stochastic equations in infinite dimensions , vol. 44 of Encyclopedia of Mathematics and its Applications . Cambridge University Press, Cambridge, 1992.