We study the spatial discretisation of the Cauchy problem for a multidimensional linear parabolic PDE of second order, with nondivergent operator and unbounded time- and space-dependent coefficients. The equation free term and the initial data are also allowed to grow. Under a nondegeneracy assumption, we consider the PDE solvability in the framework of the variational approach and approximate in space the PDE problem's generalised solution, with the use of finite-difference methods. The rate of convergence is estimated.
1. Introduction
In this paper, we study the discretisation in space of the Cauchy problem
(1)∂u∂t=Lu+fin[0,T]×ℝd,u(0,x)=g(x)onℝd,
where L is the second-order partial differential operator in the nondivergence form
(2)L(t,x)=aij(t,x)∂2∂xi∂xj+bi(t,x)∂∂xi+c(t,x),iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii,j=1,…,d,
with real coefficients (written with the usual summation convention), f and g are given real-valued functions, and T∈(0,∞) is a constant. We assume that operator ∂/∂t-L is uniformly parabolic and allows the growth in the spatial variables of the first- and second-order coefficients in L (linear and quadratic growth, resp.) and of the data f and g (polynomial growth).
Multidimensional partial differential equation (PDE) problems arise in Financial Mathematics and in Mathematical Physics. We are mainly motivated by the application to a large class of stochastic models in Financial Mathematics comprising the non-path-dependent options, with fixed exercise, written on multiple assets (basket options, exchange options, compound options, European options on future contracts and foreign-exchange, and others) and also to a particular type of path-dependent options: the Asian options (see, e.g., [1]).
Let us consider the stochastic modelling of a multiasset option of European type under the framework of a general version of Black-Scholes model, where the vector of asset appreciation rates and the volatility matrix are taken to be time- and space-dependent, and the riskless interest rate is a function of time. Owing to a Feynman-Kač type formula, pricing this option can be reduced to solving the Cauchy problem (with terminal condition) for the degenerate second-order linear parabolic PDE of nondivergent type, with null term and unbounded coefficients (see, e.g., [1]),
(3)∂V∂t+12σij(t,S)SiSj∂2V∂Si∂Sj+r(t)Si∂V∂Si-r(t)V=0in[0,T]×ℝ+d,V(T,S)=ϕ(S)onℝ+d,
where ℝ+d≡{x∈ℝd:xi>0,i=1,…,d}, V is the (unknown) option value, Si the price of the ith underlying asset, (σij) the volatility matrix, r the risk-free interest rate, and ϕ the pay-off function.
Therefore, as an alternative to approximating the option price with probabilistic numerical methods, we can approximate the solution of the corresponding PDE problem (3) with the use of nonprobabilistic techniques.
When problem (1) is considered in connection with option pricing, we see that the growth of the Black-Scholes PDE coefficients is appropriately matched. Also, the general case where the asset appreciation rate vector, the volatility matrix, and the risk-free interest rate are variable is covered. Finally, by imposing weak conditions on the initial data g, we will allow the financial derivative pay-off to be specified in a large class of functions. The free term f is included to further improve generality.
In this paper, we study the approximation in space (for the time approximation, we refer to [2–4], where a general evolution equation problem of parabolic type is discretised) of the second-order parabolic problem (1), in the challenging case where the coefficients are unbounded (as well as the free data f and g). The results are obtained under the strong assumption that the PDE does not degenerate but by imposing weak regularity assumptions. In order to facilitate the approach, we avoid any numerical methods’ sophistication and make use of basic one-step finite-difference schemes. Also, an estimate for the rate of convergence of the discretised problem's generalised solution to the exact problem's generalised solution is given.
The numerical methods and possible approximation results are strongly linked to the theory on the solvability of the PDEs. We make use of the L2 theory of solvability of linear PDEs in weighted Sobolev spaces. In particular, we consider the PDE solvability in a class of weighted Sobolev spaces used by O. G. Purtukhia (the references for Purtukhia's works can be found in [5]) for the treatment of linear stochastic partial differential equations (SPDEs) and further generalised by Gyöngy and Krylov (see [5]), the so-called well-weighted Sobolev spaces. By constructing discrete versions of these spaces, we set a suitable discretised framework and investigate the spatial approximation to the PDE generalised solution with the use of standard variational techniques.
We emphasize some points.
Firstly, we note that many PDE problems related to finance are Cauchy problems: initial-boundary value problems arise mostly after a localisation procedure for the purpose of obtaining implementable numerical schemes. Therefore, we do not find in many of these problems the complex domain geometries which are one important reason to favour other numerical methods (e.g., finite-element methods).
Also, although the finite-difference method for approximating PDEs is a well-developed area, and the theory could be considered reasonably complete since three decades ago, some important research is still currently pursued (see, e.g., the recent works [8–10]). (We refer to [6] for a brief summary of the method's history, and also for the references of the seminal works by R. Courant, K. O. Friedrichs, and H. Lewy, and further major contributions by many others. For the application of the finite-difference method to option pricing, we refer to the review paper [7] for the references of the original publications by M. Brennan and E. S. Schwartz and further major research.)
Secondly, we observe that the usual procedure for obtaining implementable numerical schemes for problem (1) is to localise it to a bounded domain in [0,T]×ℝd and then to approximate the localised problem (see, e.g., [11–13]; see also [14], where the approximation is pursued for more complex financial models but using the same localisation technique). In this case, there is no need to consider weighted functional spaces for the solvability and approximation study, as the PDE coefficients are bounded in the truncated domain.
An alternative procedure is to (semi)discretise problem (1) in [0,T]×Zhd, with Zhd being the h-grid on ℝd, and then localise the discretised problem to a bounded domain in [0,T]×Zhd by imposing a discrete artificial boundary condition (see, e.g., [15–17], where several types of initial-value problems on unbounded domains are approximated; we refer to [16, 17] for the procedure discussion). Our study is meaningful in this latter case, as the coefficient unboundedness remains a problem that must be dealt with.
Finally, we remark that (i) the partial differential operators arising in finance are of nondivergent type and (ii) we do not assume the operator coefficients to be smooth enough to be possible to obtain an equivalent divergent operator. Therefore, although there are definite advantages in considering the operator in the divergent form when the approach is variational, this is not available for the present work.
We outline the paper. In Section 2, we establish some well-known facts on the solvability of linear PDEs under a general framework and introduce the well-weighted Sobolev spaces. In Section 3, we discretise in space problem (1), with the use of a finite-difference scheme. We set a discrete framework and deduce the existence and uniqueness of the discretised problem's generalised solution. In Section 4, we investigate the approximation properties of the scheme and compute a rate of convergence. In Section 5, we make a few final comments.
2. Preliminaries and Classical Results
We establish some facts on the solvability of PDEs under a general framework.
Let V be a reflexive separable Banach space embedded continuously and densely into a Hilbert space H with inner product (,). Then H*, the dual space of H, is also continuously and densely embedded into V*, the dual of V. Let us use the notation 〈,〉 for the duality. Let H* be identified with H in the usual way, by the help of the inner product. Then we have the so-called normal triple V↪H≡H*↪V*, with continuous and dense embeddings.
Let us consider the Cauchy problem for an evolution equation
(4)dudt=A(t)u+f(t)in[0,T],u(0)=g,
with T∈(0,∞) and where for every t∈[0,T]A(t) is a linear operator from V to V*, f(t)∈V*, and g∈H.
We assume that the operator A(t) is continuous and impose a coercivity condition, as well as some regularity on the free data f and g.
Assumption 1.
There exist constants λ>0, K, M, and N such that
〈A(t)v,v〉+λ∥v∥V2≤K∥v∥H2, ∀v∈V, ∀t∈[0,T];
∥A(t)v∥V*≤M∥v∥V, ∀v∈V, ∀t∈[0,T];
∫0T∥f(t)∥V*2dt≤N and ∥g∥H≤N.
We define the generalised solution of problem (4).
Definition 2.
One says that u∈C([0,T];H) is a generalised solution of (4) on [0,T] if
u∈L2([0,T];V);
for every t∈[0,T],
(5)(u(t),v)=(g,v)+∫0t〈A(s)u(s),v〉ds+∫0t〈f(s),v〉ds
holds for all v∈V.
Notation 1.
Let W be a Banach space with norm ∥∥. We denote by C([0,T];W) the space of continuous W-valued functions on [0,T] and by L2([0,T];W) the space of W-valued functions w on [0,T] such that ∥w∥L2([0,T];W):=(∫0T∥w(t)∥2dt)1/2<∞.
The following well-known result states that, under Assumption 1, problem (4) has a unique generalised solution (see, e.g., [5], where the result is stated for the more general case of a linear stochastic evolution equation problem and [18]).
Theorem 3.
Under (1)–(3) in Assumption 1, problem (4) has a unique generalised solution on [0,T]. Moreover
(6)supt∈[0,T]∥u(t)∥H2+∫0T∥u(t)∥V2dt≤N(∥g∥H2+∫0T∥f(t)∥V*2dt),
where N is a constant.
Let us now consider the particular PDE problem
(7)ut=Lu+finQ,u(0,x)=g(x)onℝd,
where L is the second-order operator with real coefficients
(8)L(t,x)=aij(t,x)∂2∂xi∂xj+bi(t,x)∂∂xi+c(t,x),Q=[0,T]×ℝd, with T∈(0,∞), and f and g are given functions. We allow the growth, in the spatial variables, of the coefficients aij(t,x) and bi(t,x), i,j=1,…,d, and of the free data f(t,x) and g(x).
In order to set the framework for problem (7), we introduce a suitable class of weighted Sobolev spaces, the so-called well-weighted Sobolev spaces (we refer to [5] for a complete description of this class of spaces).
Let U be a domain in ℝd, that is, an open subset of ℝd. Let r>0, ρ>0 be smooth functions in U and m≥0 an integer. The weighted Sobolev space Wm,2(U;r,ρ) is the Banach space of locally integrable functions v:U→ℝ such that for each multi-index α, with |α|≤m, Dαv exists in the weak sense, and
(9)∥v∥Wm,2(U;r,ρ):=(∑|α|≤m∫Ur2|ρ|α|Dαv|2dx)1/2<∞.
Endowed with the inner product which generates the above norm
(10)(v,w)Wm,2(U;r,ρ):=∑|α|≤m∫Ur2ρ2|α|DαvDαwdx,
for all v,w∈Wm,2(U;r,ρ), Wm,2(U;r,ρ) is a Hilbert space.
Remark 4.
Setting the weight functions r=ρ=1, for all x∈U, we obtain the particular case of the Sobolev spaces Wm,2(U).
Notation 2.
In the sequel, when U=ℝd we drop the argument in the function space notation. For instance, we denote Wm,2(ℝd;r,ρ)=:Wm,2(r,ρ).
We make some assumptions on the behaviour of the weight functions r and ρ (see [5]).
Assumption 5.
Let m≥0 be an integer and r>0, ρ>0 smooth functions on ℝd. There exists a constant K such that
|Dαρ|≤Kρ1-|α|, for all multi-indexes α such that |α|≤m-1 if m≥2;
|Dαr|≤K(r/ρ|α|), for all multi-indexes α such that |α|≤m.
Remark 6.
In (1) in Assumption 5, if m<2 nothing is required.
Example 7.
The following functions (taken from [5], citing O. G. Purtukhia) satisfy Assumption 5:
ρ(x)=exp(-(1+|x|2)γ), γ≥0; each weight function r(x) in examples (1)–(5).
Now, we switch our point of view and consider the functions w:Q→ℝ as mappings of t into certain spaces of functions of x we make precise below such that, for all t∈[0,T], x∈ℝd, (w(t))(x):=w(t,x).
We impose a coercivity condition and make some assumptions on the growth and regularity of the operator's coefficients and also on the regularity of the free data f and g (see [5], where the assumptions are made for the more general case of an SPDE problem).
Assumption 8.
Let r>0, ρ>0 be smooth functions on ℝd and m≥0 an integer.
There exists a constant λ>0 such that ∑i,j=1daij(t,x)ξiξj≥λρ2(x)|ξ|2, for all t≥0, x∈ℝd, and ξ∈ℝd.
The coefficients of L are measurable functions in [0,T]×ℝd. The derivatives in x of the coefficients aij up to order m∨1 and of the coefficients bi and c up to order m exist for any t∈[0,T]. Moreover, there exists a constant K such that
(11)|Dxαaij|≤Kρ2-|α|∀|α|≤m∨1,|Dxαbi|≤Kρ1-|α|,|Dxαc|≤K∀|α|≤m,
for all t∈[0,T], x∈ℝd, with Dxα denoting the |α|th partial derivative operator with respect to x.
Consider f∈L2([0,T];Wm-1,2(r,ρ)) and g∈Wm,2(r,ρ).
Notation 3.
In the above assumption, p∨q:=max(p,q), with p, q integers. Also, for m=0 we use the notation W-1,2(r,ρ):=(W1,2(r,ρ))*, where (W1,2(r,ρ))* is the dual of W1,2(r,ρ).
We define the generalised solution of problem (7).
Definition 9.
One says that u∈C([0,T];W0,2(r,ρ)) is a generalised solution of (7) on [0,T] if
u∈L2([0,T];W1,2(r,ρ));
for every t∈[0,T],
(12)(u(t),φ)=(g,φ)+∫0t{-(aij(s)Dxiu(s),Dxjφ)+(bi(s)Dxiu(s)-Dxjaij(s)Dxiu(s),φ)+(c(s)u(s),φ)+〈f(s),φ〉(aij(s)Dxiu(s),Dxjφ)}ds
holds for all φ∈C0∞.
Notation 4.
The notation (,) in the above definition stands for the inner product in W0,2(r,ρ). Also, we denote C0∞:=C0∞(ℝd).
Remark 10.
Note that as an alternative to the infinite differentiability of φ in (2) it can be required that φ∈W1,2(r,ρ).
Finally, we state a result on the existence and uniqueness of the solution of problem (7). This result can be obtained from the general result in abstract spaces by using the suitable triple of spaces (see, e.g., [5], where the result is proved for the more general case of an SPDE problem).
Theorem 11.
Under (1)-(2) in Assumption 5, with m+1 in place of m, with m≥0 being an integer, and (1)–(3) in Assumption 8, problem (7) admits a unique generalised solution u on [0,T]. Moreover
(13)u∈C([0,T];Wm,2(r,ρ))∩L2([0,T];Wm+1,2(r,ρ)),sup0≤t≤T∥u(t)∥Wm,2(r,ρ)2+∫0T∥u(t)∥Wm+1,2(r,ρ)2dt≤N(∥g∥Wm,2(r,ρ)2+∫0T∥f(t)∥Wm-1,2(r,ρ)2dt),
with N being a constant.
3. The Discrete Framework
We now proceed to the discretisation of problem (7) in the spatial variable. We set a suitable discrete framework with the use of a finite-difference scheme and, by showing that discretised problem can be cast into the general problem (4), we prove an existence and uniqueness result for the discretised problem's generalised solution.
We emphasize that this study mirrors the study of problem (7), in the sense that the framework we now set is a discrete version of the framework set for problem (7), and the techniques used for proving the existence and uniqueness results are the same for both problems.
Define the h-grid on ℝd, with h∈(0,1],
(14)Zhd={x∈ℝd:x=h∑i=1deini,ni=0,±1,±2,…},
where ei, for i=1,2,…,d, is the unit vector in ℝd whose ith entry is 1.
For every x∈Zhd, denote
(15)∂i+u=∂i+u(t,x)=h-1(u(t,x+hei)-u(t,x)),∂i-u=∂i-u(t,x)=h-1(u(t,x)-u(t,x-hei))
the forward and backward difference quotients in space, respectively. Define the discrete operator
(16)Lh(t,x)=aij(t,x)∂j-∂i++bi(t,x)∂i++c(t,x).
We consider the discrete problem
(17)ut=Lhu+fhinQ(h),u(0,x)=gh(x)onZhd,
where Q(h)=[0,T]×Zhd, with T∈(0,∞), and fh and gh are functions such that fh:Q(h)→ℝ and gh:Zhd→ℝ.
Consider functions v:Zhd→ℝ. We introduce the space l0,2(r), the discrete version of the weighted Sobolev space W0,2(r,ρ),
(18)l0,2(r):={v:∥v∥l0,2(r)<∞},
where the norm ∥v∥l0,2(r) is given by
(19)∥v∥l0,2(r)=(∑x∈Zhdr2(x)|v(x)|2hd)1/2.
We define, for any v,w∈l0,2(r), the inner product
(20)(v,w)l0,2(r)=∑x∈Zhdr2(x)v(x)w(x)hd,
which induces the norm.
The inner product space l0,2(r) has a good structure: it can be easily shown that it is complete, therefore a Hilbert space.
For functions w:Zhd→ℝ, we introduce also the discrete version of the weighted Sobolev space W1,2(r,ρ), the space l1,2(r,ρ) defined by
(21)l1,2(r,ρ)={w:∥w∥l1,2(r,ρ)<∞},
with norm
(22)∥w∥l1,2(r,ρ)2:=∥w∥l0,2(r)2+∑i=1d∥ρ∂i+w∥l0,2(r)2.
We endow l1,2(r,ρ) with the inner product, inducing the norm,
(23)(w,z)l1,2(r,ρ)=(w,z)l0,2(r)+∑i=1d(ρ∂i+w,ρ∂i+z)l0,2(r),
for any functions w,z∈l1,2(r,ρ).
We want to show that the discrete framework we have set is a particular case of the general framework considered in Section 2.
It can be easily checked that l1,2(r,ρ) is a reflexive and separable Banach space, continuously and densely embedded into the Hilbert space l0,2(r).
As l1,2(r,ρ), endowed with the inner product (,)l1,2(r,ρ), is clearly a Hilbert space therefore it is reflexive, and the proof for the separability is trivial. The continuity of the embedding follows immediately from ∥v∥l0,2(r)≤∥v∥l1,2(r,ρ), for all v∈l1,2(r,ρ). Finally, the denseness can be checked by noticing that, for an arbitrary function w∈l0,2(r), and B a ball in Zhd, the function z defined by
(24)z(x)={w(x),x∈B,0,otherwise,
belongs obviously to l1,2(r,ρ) and that, for any given ε>0, ∥w-z∥l0,2(r)<ε if the diameter of B is chosen sufficiently large.
As in the previous section, we switch our viewpoint and consider the functions z:Q(h)→ℝ as mappings of t into certain spaces of functions of x, defined by (z(t))(x):=z(t,x), for all t∈[0,T] and for all x∈Zhd. For these functions, we consider the space C([0,T];l0,2(r)) of continuous l0,2(r)-valued functions on [0,T] and the spaces
(25)L2([0,T];lm,2(r,ρ))={z:[0,T]⟶lm,2(r,ρ):∫0T∥z(t)∥lm,2(r,ρ)2dt<∞},
with m=0,1.
Notation 5.
We identify l0,2(r,ρ) with l0,2(r).
Remark 12.
Clearly, if u∈C([0,T];l0,2(r)) then supt∈[0,T]∥u(t)∥l0,2(r)<∞.
We make some assumptions on the regularity of the data fh and gh in (17).
Assumption 13.
Let r>0 be a smooth function on ℝd. We assume
fh∈L2([0,T];l0,2(r));
gh∈l0,2(r).
Remark 14.
In Assumption 13, (1) can be replaced with the weaker assumption fh∈L2([0,T];(l1,2(r,ρ))*), where (l1,2(r,ρ))* denotes the dual space of l1,2(r,ρ).
Remark 15.
We note that |∂i+aij|≤Kρ can be obtained from (2) in Assumption 8. In fact, by the mean-value theorem,
(26)|∂i+aij(t,x)|=|h-1(aij(t,x+hei)-aij(t,x))|=|∂∂xiaij(t,x+τei)|,
for some τ such that 0<τ<h. Thus |(∂/∂xi)aij|≤Kρ implies |∂i+aij|≤Kρ.
We define the generalised solution of problem (17).
Definition 16.
One says that u∈C([0,T];l0,2(r))∩L2([0,T];l1,2(r,ρ)) is a generalised solution of (17) if, for every t∈[0,T],
(27)(u(t),φ)=(gh,φ)+∫0t{-(aij(s)∂i+u(s),∂j+φ)+(bi(s)∂i+u(s)-∂j+aij(s)∂i+u(s),φ)+(c(s)u(s),φ)+〈fh(s),φ〉(aij(s)∂i+u(s),∂j+φ)}ds
holds for all φ∈l1,2(r,ρ).
Notation 6.
In the above definition, (,) denotes the inner product in l0,2(r). We keep this convention for the remainder of present section.
Finally, we prove an existence and uniqueness result for the solution of the discrete problem (17). This result shows that the scheme is stable, that is, that, informally, the discrete problem's solution remains bounded independent of the space-step h. The result is obtained as a consequence of Theorem 3, remaining only to prove that, under the discrete framework we constructed, (1)-(2) in Assumption 1 hold.
Theorem 17.
Under (1)-(2) in Assumption 8 and (1)-(2) in Assumption 13, problem (17) has a unique generalised solution u in [0,T]. Moreover
(28)sup0≤t≤T∥u(t)∥l0,2(r)2+∫0T∥u(t)∥l1,2(r,ρ)2dt≤N(∥gh∥l0,2(r)2+∫0T∥fh(t)∥l0,2(r)2dt),
with N a constant independent of h.
Proof.
Let Lh(s):l1,2(r,ρ)→(l1,2(r,ρ))*, for every s∈[0,T]. We define
(29)〈Lh(s)ψ,φ〉≔-(aij(s)∂i+ψ,∂j+φ)+(bi(s)∂i+ψ-∂j+aij(s)∂i+ψ,φ)+(c(s)ψ,φ),
for all s∈[0,T], φ, ψ∈l1,2(r,ρ).
It suffices to prove that the following properties hold:
For the first property, owing to (1) and (2) in Assumption 8, we have
(30)〈Lh(s)ψ,ψ〉=-∑i,j∑xr2aij(s)∂i+ψ∂j+ψhd+∑i∑xr2(bi(s)-∂j+aij(s))∂i+ψψhd+∑xr2c(s)ψψhd≤-λ∑i∑xr2|ρ∂i+ψ|2hd+2K∑i∑xr2ρ|∂i+ψψ|hd+K∑xr2|ψ|2hd=-λ∑i∥ρ∂i+ψ∥l0,2(r)2+2K∑i∑xr2ρ|∂i+ψψ|hd+K∥ψ∥l0,2(r)2,
where the variable x∈Zhd is omitted, ∑x denotes the summation over Zhd, and ∑i, ∑j the summation over {1,2,…,d}. We use Cauchy’s inequality on the second term in estimate (30) and obtain
(31)〈Lh(s)ψ,ψ〉≤-λ∑i∥ρ∂i+ψ∥l0,2(r)2+εK∑i∑xr2|ρ∂i+ψ|2hd+Kε∑i∑xr2|ψ|2hd+K∥ψ∥l0,2(r)2=-λ∑i∥ρ∂i+ψ∥l0,2(r)2-λ∥ψ∥l0,2(r)2+εK∑i∥ρ∂i+ψ∥l0,2(r)2+Kε∥ψ∥l0,2(r)2+(K+λ)∥ψ∥l0,2(r)2≤-λ∥ψ∥l1,2(r,ρ)2+K∥ψ∥l0,2(r)2,
with λ>0,K constants, by taking ε sufficiently small. The first property is proved.
The second property follows from (2) in Assumption 8 and Cauchy-Schwarz inequality
(32)|〈Lh(s)ψ,φ〉|=|-∑i,j∑xr2aij(s)∂i+ψ∂j+φhd+∑i∑xr2bi(s)∂i+ψφhd-∑i,j∑xr2∂j+aij(s)∂i+ψφhd+∑xr2c(s)ψφhd|≤K∑i,j∑xr2|ρ2∂i+ψ∂j+φ|hd+K∑i∑xr2|ρ∂i+ψφ|hd+K∑xr2|ψφ|hd≤K∑i∥ρ∂i+ψ∥l0,2(r)∑j∥ρ∂j+φ∥l0,2(r)+K∑i∥ρ∂i+ψ∥l0,2(r)∥φ∥l0,2(r)+K∥ψ∥l0,2(r)∥φ∥l0,2(r)≤K∥ψ∥l1,2(r,ρ)·∥φ∥l1,2(r,ρ),
where the same writing conventions are kept.
Owing to Theorem 3 the result follows.
4. Approximation Results
In this section, we study the approximation properties of scheme (17). We begin by investigating the consistency of the scheme and prove that the difference quotients approximate the partial derivatives (with accuracy of order 1). In addition, we estimate the rate of convergence of the difference quotients to the partial derivatives.
The result is obtained by using a Sobolev inequality, under stronger regularity assumptions, and imposing that the weights ρ are bounded from below by a positive constant. In practice, the latter restriction amounts to assuming that the weights ρ are increasing functions of |x|, which is precisely the case we are interested in.
Also, we note that the way we set our discrete framework, in strong connection with the framework for problem (7), plays a crucial role in obtaining the convergence rate.
Theorem 18.
Let r>0 and ρ>0 be functions on ℝd and m an integer strictly greater than d/2. Assume that (1)-(2) in Assumption 5 are satisfied and also that ρ(x)≥C on ℝd, with C>0 being a constant. Let u(t)∈Wm+2,2(r,ρ), v(t)∈Wm+3,2(r,ρ), for all t∈[0,T]. Then there exists a constant N independent of h such that
The following remarks will be used in the proof of the theorem.
Under the conditions of the theorem, function u(t) (function v(t)) has a modification in x which is continuously differentiable in x up to order 2 (up to order 3), and the derivatives equal the weak derivatives, for every t∈[0,T]. This can be proved by Sobolev’s embedding of Wm,2(B) into Cn(B¯), for balls B in ℝd, if m>d/2+n (see, e.g., [18, 19]). We consider these modifications in the theorem’s proof.
Note that if U, V are open subsets of ℝd with V⊂U and w∈Wm,2(U) then w∈Wm,2(V). Also, if w∈Wm,2(U) and ζ∈C0∞(U) then ζ∈Wm,2(U) and ζw∈Wm,2(U) (see, e.g., [18, 19]).
Proof of Theorem 18.
Let us prove (1). We define a suitable geometric setting and then obtain an estimate for
(33)r2(x)|uxi(t,x)-∂i+u(t,x)|2ρ2(x),
with x∈Zhd, using Sobolev’s inequality on a fixed ball.
Let us consider d-cells
(34)Rh={(x1,x2,…,xd)∈ℝd:xhi<xi<xhi+h,i=1,2,…,d(x1,x2,…,xd)∈ℝd},
with xh=(xh1,xh2,…,xhd)∈Zhd fixed. Consider the particular d-cell, where h=1 and x1=(0,0,…,0), and denote it by R10. Now, take open balls Bh such that Bh⊃Rh, with the vertices {xhi,xhi+h,i=1,2,…,d} lying on the limiting sphere. Denote by B10 the ball containing R10.
For every xh∈Zhd, taking in mind (1) in Remark 19, we have, by the mean-value theorem,
(35)∂i+u(t,xh)=h-1(u(t,xh+hei)-u(t,xh))=uxi(t,xh+θhei),(36)|uxi(t,xh)-∂i+u(t,xh)|=|uxi(t,xh)-uxi(t,xh+θhei)|≤h|uxixi(t,xh+θ′hei)|,
for some 0<θ′<θ<1. Clearly,
(37)|uxixi(t,xh+θ′hei)|≤supx∈Rh|uxixi(t,x)|,
and then, from (36) and (37),
(38)|uxi(t,xh)-∂i+u(t,xh)|2≤h2supx∈Rh|uxixi(t,x)|2.
We change variable in order to have the supremum in (38) calculated over the fixed d-cell R10:
(39)supx∈Rh|uxixi(t,x)|=supx∈R10|uxixi(t,xh+hx)|.
As
(40)supx∈R10|uxixi(t,xh+hx)|2≤supx∈B10|uxixi(t,xh+hx)|2,
from (38)–(40) we immediately obtain
(41)r2(xh)|uxi(t,xh)-∂i+u(t,xh)|2ρ2(xh)≤h2supx∈R10(r2(xh+hx)|uxixi(t,xh+hx)|2ρ2(xh+hx))≤h2supx∈B10(r2(xh+hx)|uxixi(t,xh+hx)|2ρ2(xh+hx)).
Now, taking in mind (2) in Remark 19, we have for m>d/2 by using Sobolev’s inequality
(42)supx∈B10|r(xh+hx)uxixi(t,xh+hx)ρ(xh+hx)|2≤N∑|α|≤m∫B10|Dxα(r(xh+hx)uxixi(t,xh+hx)ρ×(xh+hx))|2dx,
with N being a constant independent of h. Observe that the Leibniz formula
(43)|Dxα(ruxixiρ)|=|∑β≤α(αβ)Dβ(rρ)Dxα-βuxixi|=|∑β≤α(αβ)(∑γ≤β(βγ)DγrDβ-γρ)Dxα-βuxixi|
holds (the arguments of r, ρ and uxixi are omitted). Also, keeping the same convention, owing to Assumption 5(44)|Dγr|≤Krρ-|γ|,|Dβ-γρ|≤Kρ1-(|β|-|γ|),
with K a constant, and then
(45)|∑γ≤β(βγ)DγrDβ-γρ|≤N∑γ≤β(βγ)rρ-|γ|ρ1-(|β|-|γ|)≤Nrρ1-|β|,
with N a constant. From (42)–(45), we get
(46)supx∈B10|r(xh+hx)uxixi(t,xh+hx)ρ(xh+hx)|2≤N∑|α|≤m∑β≤α∫B10r2(xh+hx)|ρ1-|β|(xh+hx)|2×|Dxα-βuxixi(t,xh+hx)|2dx.
Note also that, owing to Hölder inequality and to the hypotheses on function ρ, the integral in (46) can be estimated by
(47)∫B10r2(xh+hx)×|ρ1-|β|(xh+hx)Dxα-βuxixi(t,xh+hx)|2dx≤N∫B10r2(xh+hx)|Dxα-βρ2+(|α|-|β|)(xh+hx)×Dxα-βuxixi(t,xh+hx)|2dx·supx∈B10|ρ-1-|α|(xh+hx)|2≤N∫B10r2(xh+hx)|Dxα-βρ2+(|α|-|β|)(xh+hx)×Dxα-βuxixi(t,xh+hx)|2dx.
Thus, from (46) and (47),
(48)supx∈B10|r(xh+hx)uxixi(t,xh+hx)ρ(xh+hx)|2≤N∑|α|≤m∑β≤α∫B10r2(xh+hx)×|Dxα-βρ2+(|α|-|β|)(xh+hx)·Dxα-βuxixi(t,xh+hx)|2dx≤N∑|α|≤m∫B10r2(xh+hx)×|ρ2+|α|(xh+hx)Dxαuxixi(t,xh+hx)|2dx≤N∑|α|≤m+2∫B10r2(xh+hx)×|ρ|α|(xh+hx)Dxαu(t,xh+hx)|2dx=N∑|α|≤m+2∫Bhr2(x)|ρ|α|(x)Dxαu(t,x)|2h-dh2|α|dx≤N∑|α|≤m+2∫Bhr2(x)|ρ|α|(x)Dxαu(t,x)|2h-ddx.
Finally, owing to the particular geometry of the framework we have set, from (41) and (48) we obtain
(49)∑x∈Zhdr2(x)|uxi(t,x)-∂i+u(t,x)|2ρ2(x)hd≤Nh2∑|α|≤m+2∑xh∈Zhd∫Bh(xh)r2(x)|ρ|α|(x)Dxαu(t,x)|2dx≤Nh2∑|α|≤m+2∑xh∈Zhd∫Rh(xh)r2(x)|ρ|α|(x)Dxαu(t,x)|2dx≤h2N∥u(t)∥Wm+2,2(r,ρ)2,
where Bh(xh):=Bh, Rh(xh):=Rh, and the proof for (1) is complete.
The proof for (2) follows the same steps.
Finally, owing to the stability and consistency properties of the scheme (Theorems 17 and 18, resp.), we prove the convergence of the discrete problem’s solution to the PDE problem's solution and compute a rate of convergence. The accuracy obtained is of order 1.
Theorem 20.
Let r>0 and ρ>0 be functions on ℝd and m an integer strictly greater than d/2. Assume that the hypotheses of Theorems 11 and 17 are satisfied and that ρ(x)≥C on ℝd, with C>0 being a constant. Denote by u the solution of problem (7) in Theorem 11 and by uh the solution of problem (17) in Theorem 17. Assume additionally that u∈L2([0,T];Wm+3,2(r,ρ)). Then
(50)sup0≤t≤T∥u(t)-uh(t)∥l0,2(r)2+∫0T∥u(t)-uh(t)∥l1,2(r,ρ)2dt≤h2N∫0T∥u(t)∥Wm+3,2(r,ρ)2dt+N(∥g-gh∥l0,2(r)2+∫0T∥f(t)-fh(t)∥l0,2(r)2dt),
with N being a constant independent of h.
Remark 21.
Under the conditions of the above theorem, there are modifications in x such that the data f(t) and g are continuous in x, for every t∈[0,T] (see Remark 19). We will consider these modifications in the proof of the theorem.
Proof of Theorem 20.
From (7) and (17), we have that u-uh satisfies the problem
(51)(u-uh)t=Lh(u-uh)+(L-Lh)u+(f-fh)inQ(h)(u-uh)(0,x)=(g-gh)(x)onZhd.
Taking in mind Remark 21, we see clearly that f-fh∈L2([0,T];l0,2(r)) and g-gh∈l0,2(r).
With respect to the term (L-Lh)u, note that if u(t)∈Wm+3,2(r,ρ), for all t∈[0,T],
(52)∑x∈Zhdr2(x)|(L-Lh)(t)u(t)|2hd=∑x∈Zhdr2(x)|aij(t,x)(∂2∂xi∂xj-∂j-∂i+)u(t,x)+bi(t,x)(∂∂xi-∂i+)u(t,x)(∂2∂xi∂xj-∂j-∂i+)|2hd≤h2N∥u(t)∥Wm+3,2(r,ρ)2<∞,
owing to (2) in Assumption 8 and to Theorem 18. Thus (L-Lh)(t)u(t)∈l0,2(r), for every t∈[0,T]. Moreover, as by assumption u∈L2([0,T];Wm+3,2(r,ρ)), we obtain immediately (L-Lh)u∈L2([0,T];l0,2(r)).
We have shown that problem (51) satisfies the hypotheses of Theorem 17, therefore holding the estimate
(53)sup0≤t≤T∥u(t)-uh(t)∥l0,2(r)2+∫0T∥u(t)-uh(t)∥l1,2(r,ρ)2dt≤N(∥g-gh∥l0,2(r)2+∫0T∥f(t)-fh(t)∥l0,2(r)2dt+∫0T∥(L-Lh)(t)u(t)∥l0,2(r)2dt).
Owing again to (2) in Assumption 8 and to Theorem 18, the result follows.
The following result is an immediate consequence of Theorem 20.
Corollary 22.
Let the hypotheses of Theorem 20 be satisfied, and denote by u the solution of (7) in Theorem 11 and by uh the solution of (17) in Theorem 17. If there is a constant N independent of h such that
(54)∥g-gh∥l0,2(r)2+∫0T∥f(t)-fh(t)∥l0,2(r)2dt≤h2N(∥g∥Wm,2(r,ρ)2+∫0T∥f(t)∥Wm-1,2(r,ρ)2dt)
then
(55)sup0≤t≤T∥u(t)-uh(t)∥l0,2(r)2+∫0T∥u(t)-uh(t)∥l1,2(r,ρ)2dt≤h2N(∫0T∥u(t)∥Wm+3,2(r,ρ)2dt+∥g∥Wm,2(r,ρ)2+∫0T∥f(t)∥Wm-1,2(r,ρ)2dt).
5. Conclusions
In this paper, we investigated the finite-difference spatial approximation of the Cauchy problem for a second-order linear parabolic PDE, in the framework of the variational approach.
By considering a suitable class of weighted Sobolev spaces, and its zero and first-order discrete versions, we could deal with the growth in space of the PDE coefficients (as well as with the spatial growth of the free data f and g). Moreover, as the framework and techniques used to study the discrete problem mirror the framework and techniques for the corresponding continuous problem we could estimate a rate of convergence.
The approximation was studied under the strong assumption that the PDE does not degenerate. But the framework we used is broadly the appropriate framework for a future investigation of the related degenerate case.
Other possible further research directions include the use of splitting-up methods (see [10]), following Richardson's idea to accelerate numerical schemes, and also the use of techniques reducing the volume of computational work (e.g., sparse grid techniques), in order to deal with the computational challenge posed by the possible high dimensionality of the problem.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by FCT, Fundação para a Ciência e Tecnologia (PEst-OE/EGE/UI0491/2011, Scientific Cooperation Agreement between Portugal and Slovakia—Analysis of Nonlinear Partial Differential Equations in Mathematical Finance), and by the European Commission, Marie Curie Multi-ITN STRIKE—Novel Methods in Computational Finance (Grant no. 304617).
LambertonD.LapeyreB.1996London, UKChapman & HallMR1422250GonçalvesF. F.2007Edinburgh, UKUniversity of EdinburghGonçalvesF. F.GrossinhoM. D. R.MoraisE.Discretisation of abstract linear evolution equations of parabolic type20122012article 142910.1186/1687-1847-2012-14MR2944124ZBL06241416MoraisE.2013Lisbon, PortugalUniversity of LisbonGyöngyI.KrylovN. V.Stochastic partial differential equations with unbounded coefficients and applications. I1990321-2539110.1080/17442509008833653ThoméeV.CiarletP. G.LionsJ. L.Finite difference methods for linear parabolic equations19901Amsterdam, The NetherlandsNorth-Holland3196Handbook of Numerical AnalysisMR1039324ZBL0875.65080BroadieM.DetempleJ. B.Option pricing: valuation models and applications2004509114511772-s2.0-494424040310.1287/mnsc.1040.0275GyöngyI.KrylovN. V.First derivatives estimates for finite-difference schemes2009782682019204610.1090/S0025-5718-09-02229-7MR2521277ZBL1198.65165GyöngyI.KrylovN. V.Higher order derivative estimates for finite-difference schemes for linear elliptic and parabolic equations200916218721610.4310/MAA.2009.v16.n2.a3MR2563747ZBL1183.65107GyöngyI.KrylovN. V.Accelerated finite difference schemes for second order degenerate elliptic and parabolic problems in the whole space2011802751431145810.1090/S0025-5718-2011-02478-6MR2785464ZBL1223.65068BarlesG.DaherC.RomanoM.Convergence of numerical schemes for parabolic equations arising in finance theory19955112514310.1142/S0218202595000085MR1315000ZBL0822.65056LötstedtP.PerssonJ.von SydowL.TyskJ.Space-time adaptive finite difference method for European multi-asset options20075381159118010.1016/j.camwa.2006.09.014MR2327671ZBL1154.91462TangmanD. Y.GopaulA.BhuruthM.Numerical pricing of options using high-order compact finite difference schemes2008218227028010.1016/j.cam.2007.01.035MR2427213ZBL1146.91338PooleyD. M.ForsythP. A.VetzalK. R.Numerical convergence properties of option pricing PDEs with uncertain volatility200323224126710.1093/imanum/23.2.241MR1974225ZBL1040.91053DouglasJ.Jr.RussellT. F.Numerical methods for convection-dominated diffusion problems based on combining the method of characteristics with finite element or finite difference procedures198219587188510.1137/0719063MR672564ZBL0492.65051EhrhardtM.MickensR. E.Finite difference schemes on unbounded domains20052River Edge, NJ, USAWorld Scientific34338410.1142/9789812703316_0008MR2212203ZBL1085.65076EhrhardtM.Discrete transparent boundary conditions for Schrödinger-type equations for non-compactly supported initial data200858566067310.1016/j.apnum.2007.02.002MR2407739ZBL1138.65069LionsJ.-L.MagenesE.19681Paris, FranceDunod, Gauthier-VillarsMR0247243EvansL. C.199819Providence, RI, USAAmerican Mathematical SocietyGraduate Studies in MathematicsMR1625845