Loading [MathJax]/jax/output/SVG/jax.js
Research article

A Fourier cosine expansion method for pricing FX-TARN under Lévy processes

  • Received: 22 April 2023 Revised: 14 May 2023 Accepted: 22 May 2023 Published: 29 May 2023
  • JEL Codes: C63, C69, G13

  • In this paper, we extend the Fourier cosine expansion (COS) method to the pricing of {foreign exchange} target redemption note (FX-TARN), a popular exotic currency option. We take the FX spot rate and the cumulated positive cash flow as two state variables and factor the joint distribution by two marginals that can be approximated by Fourier cosine expansions. To recover the Fourier coefficients recursively, we approximate the two-dimensional integration by higher-order quadratures such as Gauss-Legendre or Clenshaw-Curtis quadrature for the integration over the spot rate. We derive the analytical formulas for the price under different knock-out types. We demonstrate that fast Fourier transform (FFT) can be employed to obtain the Fourier coefficients efficiently. We also evaluate the performance and accuracy of the method through a number of numerical experiments.

    Citation: Kevin Z. Tong. A Fourier cosine expansion method for pricing FX-TARN under Lévy processes[J]. Quantitative Finance and Economics, 2023, 7(2): 261-286. doi: 10.3934/QFE.2023014

    Related Papers:

    [1] Jens Lorenz, Wilberclay G. Melo, Suelen C. P. de Souza . Regularity criteria for weak solutions of the Magneto-micropolar equations. Electronic Research Archive, 2021, 29(1): 1625-1639. doi: 10.3934/era.2020083
    [2] Yazhou Wang, Yuzhu Wang . Regularity criterion of three dimensional magneto-micropolar fluid equations with fractional dissipation. Electronic Research Archive, 2024, 32(7): 4416-4432. doi: 10.3934/era.2024199
    [3] Xiuli Xu, Lian Yang . Global well-posedness of the 3D nonlinearly damped Boussinesq magneto-micropolar system without heat diffusion. Electronic Research Archive, 2025, 33(4): 2285-2294. doi: 10.3934/era.2025100
    [4] Noelia Bazarra, José R. Fernández, Ramón Quintanilla . Numerical analysis of a problem in micropolar thermoviscoelasticity. Electronic Research Archive, 2022, 30(2): 683-700. doi: 10.3934/era.2022036
    [5] J. Bravo-Olivares, E. Fernández-Cara, E. Notte-Cuello, M.A. Rojas-Medar . Regularity criteria for 3D MHD flows in terms of spectral components. Electronic Research Archive, 2022, 30(9): 3238-3248. doi: 10.3934/era.2022164
    [6] Jie Zhang, Gaoli Huang, Fan Wu . Energy equality in the isentropic compressible Navier-Stokes-Maxwell equations. Electronic Research Archive, 2023, 31(10): 6412-6424. doi: 10.3934/era.2023324
    [7] Anatoliy Martynyuk, Gani Stamov, Ivanka Stamova, Yulya Martynyuk–Chernienko . Regularization scheme for uncertain fuzzy differential equations: Analysis of solutions. Electronic Research Archive, 2023, 31(7): 3832-3847. doi: 10.3934/era.2023195
    [8] Cheng He, Changzheng Qu . Global weak solutions for the two-component Novikov equation. Electronic Research Archive, 2020, 28(4): 1545-1562. doi: 10.3934/era.2020081
    [9] Xun Wang, Qunyi Bie . Energy equality for the multi-dimensional nonhomogeneous incompressible Hall-MHD equations in a bounded domain. Electronic Research Archive, 2023, 31(1): 17-36. doi: 10.3934/era.2023002
    [10] Lin Shen, Shu Wang, Yongxin Wang . The well-posedness and regularity of a rotating blades equation. Electronic Research Archive, 2020, 28(2): 691-719. doi: 10.3934/era.2020036
  • In this paper, we extend the Fourier cosine expansion (COS) method to the pricing of {foreign exchange} target redemption note (FX-TARN), a popular exotic currency option. We take the FX spot rate and the cumulated positive cash flow as two state variables and factor the joint distribution by two marginals that can be approximated by Fourier cosine expansions. To recover the Fourier coefficients recursively, we approximate the two-dimensional integration by higher-order quadratures such as Gauss-Legendre or Clenshaw-Curtis quadrature for the integration over the spot rate. We derive the analytical formulas for the price under different knock-out types. We demonstrate that fast Fourier transform (FFT) can be employed to obtain the Fourier coefficients efficiently. We also evaluate the performance and accuracy of the method through a number of numerical experiments.



    In this paper we present some regularity criteria for weak solutions of the following magneto-micropolar fluid system in three space dimensions:

    {ut+uu+(p+12|b|2)=(μ+χ)Δu+bb+χ×w,wt+uw=γΔw+κ(w)+χ×u2χw,bt+ub=νΔb+bu,u=b=0,u(,0)=u0(),w(,0)=w0(),b(,0)=b0(), (1)

    where u(x,t)=(u1(x,t),u2(x,t),u3(x,t))R3 denotes the incompressible velocity field, w(x,t)=(w1(x,t),w2(x,t), w3(x,t))R3 denotes the micro-rotational velocity, b(x,t) =(b1(x,t),b2(x,t), b3(x,t)) R3 the divergence free magnetic field and p(x,t)R the hydrostatic pressure. As usual, xR3 denotes the space variable and 0t<T the time variable. The positive constants μ,χ,ν,κ, and γ are associated with specific properties of the fluid; precisely, μ is the kinematic viscosity, χ is the vortex viscosity, κ and γ are spin viscosities and, lastly, ν1 is the magnetic Reynolds number. The initial data for the velocity and magnetic fields, given by u0 and b0 in (1), are divergence free, i.e., u0=b0=0.

    Let us list some recent papers which discuss regularity of weak solutions of the magneto-micropolar equations (1) and systems that are particular cases of these equations as, for example, the classical Navier-Stokes equations.

    In 2010, Y. Baoquan [18] established the following regularity criteria: a solution (u,w,b)(t) of (1) can be extended smoothly beyond t=T if

    T0u(t)βαdt<,where 3α+2β1,3<α,

    or

    T0u(t)βαdt<,where 3α+2β2,32<α,

    provided that the initial data (u0,w0,b0) belong to H1(R3) and u0=b0=0. Here (u,w,b)C([0,T);H1(R3))C((0,T);H2(R3)) is assumed. Notice that Y. Baoquan [18] used only the velocity field or its gradient in the regularity criteria described above.

    In 2013, Y. Wang [14] showed that a weak solution (u,w,b)(t) of (1), defined on the interval [0,T), can be extended smoothly beyond t=T if one assumes

    T03u(t)βαdt<where 3α+2β1,α3, (2)

    provided that (u0,w0,b0)H1(R3) and u0=b0=0. Let us point out that Y. Wang [14] established the regularity criterion given above by considering only one component of the gradient of the velocity field (see (2)). The paper [14] deals with an extension of a regularity criterion obtained for weak solutions of the magneto-hydrodynamic equations (see [2]). This latter system is the special case of(1) obtained for w=0 and χ=0. Further results related to the current paper are in [2,21].

    The papers [14] and [18] (see also [3,8,12,13,15,17,19,20,22,23,24,28,29]), raised our interest to obtain regularity criteria for weak solutions of the magneto-micropolar system (1), which involve only one component of the velocity field u.

    Note that the magneto-micropolar system (1) reduces to the classical Navier-Stokes equations for the velocity field u and the pressure p if one assumes that w=b=0 and χ=0. Regularity criteria for weak solutions of Navier-Stokes equations, involving only one component of the velocity field, have been published recently. We want to comment on two results.

    First, Z. Zhang and X. Yang [22] present a regularity criterion for the Navier-Stokes equations involving the gradient of one component of the velocity field. Precisely, if u(t) is a weak solution of the Navier-Stokes equations in (0,T) and

    T0u3(t)3272dt<, (3)

    then u is C in R3×(0,T) and the solution can be extended as a C function beyond T. Here u0L2(R3) with u0=0 is assumed. In the current paper we prove that the criterion (3) can be extended from the Navier-Stokes equations to the magneto-micropolar system (1) under appropriate conditions on the other fields w and b. See Theorem 1.1 below.

    The second paper we cite is [23]. Z. Zhang and X. Yang [23] deal with one component of the gradient of one component of the velocity field. More precisely, regularity of a weak solution u(t) of the Navier-Stokes equations is obtained under the assumption

    3u3(t)L(0,T;L2(R3)). (4)

    Theorem 1.2 below establishes an extension of criterion (4) from the Navier-Stokes equations to the magneto-micropolar system (1).

    Further regularity results for weak solutions of the Navier-Stokes equations are established in [1,4,5,6,7,10,25,26].

    The main results of the current paper are:

    Theorem 1.1. Let (u0,w0,b0)L2(R3) with u0,b0=0. Let T>0 and let

    (u,w,b)C([0,T);H1(R3))C((0,T);H2(R3)) (5)

    denote a weak solution of the magneto-micropolar equations (1) in (0,T) satisfying the initial condition (u,w,b)(0)=(u0,w0,b0). If

    u3,hw,hbL327(0,T;L2(R3)), (6)

    then (u,w,b) is C in R3×(0,T) and the solution can be extended as a C function beyond T.

    Theorem 1.2. Let (u0,w0,b0)L2(R3) with u0,b0=0. Let T>0 and let

    (u,w,b)C([0,T);H1(R3))C((0,T);H2(R3)) (7)

    denote a weak solution of the magneto-micropolar equations (1) in (0,T) satisfying the initial condition (u,w,b)(0)=(u0,w0,b0). If

    3u3,3w,3bL(0,T;L2(R3)), (8)

    then (u,w,b) is C in R3×(0,T) and the solution can be extended as a C function beyond T.

    An outline of the paper follows: There are two sections after the Introduction. In Section 2, we list definitions and notations used throughout the paper and recall results that play an important role in our proofs of the main results. Section 3 presents the proofs of Theorems 1.1 and 1.2.

    We introduce notations and definitions used in the paper.

    ● Boldface letters denote vector fields; for example,

    a=a(x,t)=(a1(x,t),a2(x,t),a3(x,t)),xR3,t0.

    ● The Euclidean norm of any vector a=(a1,...,an)Rn is denoted and given by |a|=ni=1a2i.

    ● The notation Lα is used for the standard Lebesgue space equipped with the norm α, where 1α; more specifically,

    aα:=(R3|a(x)|αdx)1α,1α<, 

    and

    a:=ess supxR3{|a(x)|}, 

    where a:R3Rn (nN) is a mensurable function. We define the L2-inner product of two vector functions by

    (a,b)2:=R3a(x)b(x)dx, 

    where cd:=ni=1cidi for c=(c1,...,cn),d=(d1,...,dn)Rn; and a,b:R3Rn (nN) are mensurable functions.

    ● Let a=(a1,...,an) denote the gradient of a=(a1,...,an)Rn, where aj=(1aj,2aj,3aj), with i=/xi for all i=1,2,3 and j=1,...,n.

    ● The horizontal gradient is denoted by ha=(ha1,...,han), where a=(a1,...,an)Rn and haj=(1aj,2aj), with j=1,...,n.

    ● Here a:=3i=1aii, where a=(a1,a2,a3)R3.

    ×a denotes the curl of the vector field a=(a1,a2,a3)R3.

    ● Denote a=3i=1iai, where a=(a1,a2,a3)R3.

    Δ represents the standard Laplacian operator.

    ● The horizontal Laplacian is denoted by Δha=(Δha1,...,Δhan), where a=(a1,...,an)Rn and Δhaj=2i=12iaj, with j=1,...,n.

    ● Let (X,) be a Banach space and assume that 1β, c,dR,c<d. We denote by Lβ(c,d;X) the space of all measurable functions f:(c,d)X with f()Lβ(c,d) endowed with the norm fLβ(c,d;X):=(dcf(t)βdt)1β, where β<, and also fL(c,d;X)=ess supc<t<d{f(t)}.

    ● We define a weak solution of (1) as follows: Let T>0 and (u0,w0,b0)L2(R3), with u0=b0=0. A mensurable function (u,w,b)(t) is called a weak solution for (1) on [0,T) if the following conditions hold

    1. (u,w,b)(t)L(0,T;L2(R3))L2(0,T;H1(R3));

    2. the system (1) is satisfied in the sense of distributions;

    3. the energy inequality holds, i.e.,

    (u,w,b)(t)22+2(μ+χ)t0u(τ)22dτ+2γt0w(τ)22dτ+2νt0b(τ)22dτ+2κt0w(τ)22dτ+2χt0w(τ)22dτ(u0,w0,b0)22, (9)

    for all 0t<T.

    ● For brevity, dependencies on the variables x and t are often suppressed in our notation. For example, the function tu(,t)2 may also be written as u(t)2 or u2. Furthermore, as usual, the value of constants may change line by line in the paper.

    Now, we enunciate the lemmas that will be applied in the proofs of our main results. The first one is proved in [16].

    Lemma 2.1 (see [16]). Let i,j,k be a permutation of 1,2,3. Assume that

    f,g,ig,jg,h,jh,khL2(R3).

    Then,

    R3fghdxCf2g142ig122jg142h142kh122jh142.

    The second one was established in [23].

    Lemma 2.2 (see [23]). Let fL6(R3) and gL2(R3) with 3fL2(R3) and hgL2(R3). Then,

    R3f2g2dx22f3263f122g2hg2.

    The third one was written in [2].

    Lemma 2.3 (see [2]). Assume that θ,λ,ϑR satisfy

    1θ,λ,ϑ<,1θ+2λ>1,1+3ϑ=1θ+2λ.

    Consider that fH1(R3),1f,2fLλ(R3) and 3fLθ(R3). Then, there exists a positive constant C such that

    fϑC1f13λ2f13λ3f13θ.

    In particular, if λ=2 and fH1(R3), 3fLθ(R3) (with 1θ<), then there is a positive constant C such that

    f3θC1f1322f1323f13θ.

    In this section we prove Theorems 1.1 and 1.2. In both results, it is necessary to consider ϵ(0,T) arbitrary in order to obtain δ(0,ϵ) such that (u,w,b)(δ)L2(R3). It is known that there is a unique strong solution (u,w,b)C([δ,T);H1(R3))L2(δ,T;H2(R3)) for the system (1) (see [9,11] and references therein) that satisfies (u,w,b)(δ)=(u,w,b)(δ), where t=T is the maximum time of existence for this solution, and (u,w,b)C(R3×(0,T)) (since ϵ>0 is arbitrary). Therefore, if T<T, one concludes that (u,w,b)(t) is regular in R3×(0,T). On the other hand, assuming TT, we will prove below that (u,w,b)(t)2 is uniformly bounded for t[δ,T) (in Theorem 1.2) and bounded as tT (in Theorem 1.1). However, this is not possible. In fact, this boundedness would imply that (u,w,b)(t) could be extended beyond t=T, which contradicts the definition of T.

    First, notice that by applying the product (,Δhu)2 to the first equation of the magneto-micropolar system (1) and using the fact that u=0, one gets

    12ddthu22+(μ+χ)hu22=(uu,Δhu)2(bb,Δhu)2χ(×w,Δhu)2. (10)

    Similarly, from the second and third equations in (1), we obtain

    12ddthw22+γhw22+κh(w)22+2χhw22=(uw,Δhw)2χ(×u,Δhw)2 (11)

    and also

    12ddthb22+νhb22=(ub,Δhb)2(bu,Δhb)2. (12)

    By adding the results (10), (11) and (12) one infers

    12ddt(hu,hw,hb)22+(μ+χ)hu22+γhw22+νhb22+κh(w)22+2χhw22=(uu,Δhu)2(bb,Δhu)2χ(×w,Δhu)2+(uw,Δhw)2χ(×u,Δhw)2+(ub,Δhb)2(bu,Δhb)2. (13)

    Let us examine the terms on the right hand side of the above equality. We have

    (bb,Δhu)2=3i,j=12k=1R3biibj2kujdx=3i,j=12k=1R3kbiibjkujdx+3i,j=12k=1R3bikibjkujdx.

    Similarly, we get

    (bu,Δhb)2=3i,j=12k=1R3biiuj2kbjdx=3i,j=12k=1R3kbiiujkbjdx3i,j=12k=1R3bikibjkujdx,

    where we have used that b is divergence free. Hence,

    (bb,Δhu)2(bu,Δhb)2=3i,j=12k=1R3kbiibjkujdx+3i,j=12k=1R3kbiiujkbjdxC(R3|hb||b||hu|dx+R3|hb||u||hb|dx).

    Furthermore,

    (uw,Δhw)2=3i,j=12k=1R3uiiwj2kwjdx=3i,j=12k=1R3kuiiwjkwjdx3i,j=12k=1R3uikiwjkwjdx.

    On the other hand, by analysing the last term above it is easy to prove that it is actually null. In fact, since u=0, one has

    3i,j=12k=1R3uikiwjkwjdx=3i,j=12k=1R3uikwjikwjdx. 

    Therefore,

    (uw,Δhw)2CR3|hw||w||hu|dx. 

    Similarly, we obtain

    (ub,Δhb)2CR3|hb||b||hu|dx. 

    Also notice that, by applying Cauchy-Schwarz's inequality, one has

    χ(×u,Δhw)2χ(×w,Δhu)2χhu22+χhw22. (14)

    By [22], the following estimate holds:

    (uu,Δhu)2CR3|u3||u||hu|dx. 

    Consequently, from (13), we infer

    12ddt(hu,hw,hb)22+μhu22+γhw22+νhb22+κh(w)22+χhw22CR3|(u3,hw,hb)||(u,w,b)||(hu,hw,hb)|dx.

    By applying Lemma 2.1, one obtains

    12ddt(hu,hw,hb)22+μhu22+γhw22+νhb22+κh(w)22+χhw22C(u3,hw,hb)2(hu,hw,hb)142(u,w,b)142×(hu,hw,hb)322.

    By using Young's inequality, it follows that

    ddt(hu,hw,hb)22+α(hu,hw,hb)22+2κh(w)22+2χhw22C(u3,hw,hb)42(hu,hw,hb)2×(u,w,b)2,

    where α=min{μ,ν,γ}. Now, by integrating over [Tτ,t] the inequality above (τ will be chosen a posteriori), we obtain

    (hu,hw,hb)(t)22+αtTτ(hu,hw,hb)(s)22dsC+CtTτ(u3,hw,hb)(s)42(hu,hw,hb)(s)2×(u,w,b)(s)2ds. (15)

    In order to estimate the term (u,w,b)(t)2 for all t[Tτ,T), we define

    I(t):=sups[Tτ,t]{(hu,hw,hb)(s)2}+(tTτ(hu,hw,hb)(s)22ds)12 (16)

    and

    J(t):=sups[Tτ,t]{(u,w,b)(s)2}+(tTτ(Δu,Δw,Δb)(s)22ds)12, (17)

    where Tτt<T. First of all, let us establish a relationship between I and J. By (15), one has

    I2(t)2sups[Tτ,t]{(hu,hw,hb)(s)22}
    +2tTτ(hu,hw,hb)(s)22dsC+CI(t)J(t)34tTτ(u3,hw,hb)(s)42(u,w,b)(s)142dsC+CI(t)J(t)34T0(u3,hw,hb)(s)3272ds+CI(t)J(t)34T0(u,w,b)(s)22dsC+CI(t)J(t)34,

    where we have applied Young's inequality, (6) and (9). By using Young's inequality again, we get

    I2(t)C+CJ32(t)+12I2(t), 

    or equivalently,

    I(t)C+CJ34(t),t(Tτ,T). (18)

    The inequality (18) is useful to prove that (u,w,b)(t)2 is bounded in the interval [Tτ,T). In order to establish this last statement, we start noting that the system (1) lets us conclude the following:

    12ddtu22+(μ+χ)Δu22=(uu,Δu)2(bb,Δu)2χ(×w,Δu)2, 
    12ddtw22+γΔw22+κ(w)22+2χw22=(uw,Δw)2χ(×u,Δw)2

    and also

    12ddtb22+νΔb22=(ub,Δb)2(bu,Δb)2, 

    where we used the fact that u=0. Hence, by adding the three equalities above, one obtains

    12ddt(u,w,b)22+(μ+χ)Δu22+γΔw22+νΔb22+κ(w)22+2χw22=(uu,Δu)2(bb,Δu)2χ(×w,Δu)2+(uw,Δw)2χ(×u,Δw)2+(ub,Δb)2(bu,Δb)2. (19)

    Let us examine all the terms on the right hand side of the equality above. We have

    (uw,Δw)2=3i,j,k=1R3uiiwj2kwjdx=3j=12k=1R3u33wj2kwjdx+3j,k=12i=1R3uiiwj2kwjdx+3j=1R3u33wj23wjdx=:I1(t)+I2(t)+I3(t). (20)

    Here

    I1(t)=3j=12k=1R3u33wj2kwjdx=3j=12k=1R3ku33wjkwjdx3j=12k=1R3u3k3wjkwjdx=3j=12k=1R3ku33wjkwjdx+123j=12k=1R33u3(kwj)2dxCR3|u||w||hw|dx+CR3|u||hw|2dx.

    Similarly,

    I2(t)=3j,k=12i=1R3uiiwj2kwjdx=3j,k=12i=1R3kuiiwjkwjdx3j,k=12i=1R3uikiwjkwjdx=3j,k=12i=1R3kuiiwjkwjdx+123j,k=12i=1R3iui(kwj)2dxCR3|u||hw||w|dx+CR3|hu||w|2dx.

    By using that u is divergence free, one also has

    I3(t)=3j=1R3u33wj23wjdx=123j=1R33u3(3wj)2dx=123j=12k=1R3kuk(3wj)2dxCR3|hu||w|2dx.

    Therefore, using the above estimates, the equality (20) yields

    (uw,Δw)2CR3|u||hw|2dx+CR3|u||hw||w|dx+CR3|hu||w|2dx.

    Following the same process, we conclude that

    (ub,Δb)2CR3|u||hb|2dx+CR3|u||hb||b|dx+CR3|hu||b|2dx.

    It is important to point out that the technique applied to (uw,Δw)2 may not be useful when we consider the terms (bb,Δu)2 and (bu,Δb)2 (such expressions were obtained in (19)). However, we can argue as follows.

    (bb,Δu)2(bu,Δb)2=3i,j,k=1R3biibj2kujdx3i,j,k=1R3biiuj2kbjdx=3i,j,k=1R3kbiibjkujdx+3i,j,k=1R3bikibjkujdx+3i,j,k=1R3kbiiujkbjdx+3i,j,k=1R3bikiujkbjdx.

    Consequently,

    (bb,Δu)2(bu,Δb)2=3i,j,k=1R3kbiibjkujdx+3i,j,k=1R3bikibjkujdx+3i,j,k=1R3kbiiujkbjdx3i,j,k=1R3bikibjkujdx.

    By using b=0, we have

    (bb,Δu)2(bu,Δb)2=3i,j,k=1[R3kbiibjkujdx+R3kbiiujkbjdx] (21)
    =3j=12k=1[R3kb33bjkujdx+R3kb33ujkbjdx]+3j,k=12i=1[R3kbiibjkujdx+R3kbiiujkbjdx]+3j=1[R33b33bj3ujdx+R33b33uj3bjdx]=:J1(t)+J2(t)+J3(t). (22)

    Let us estimate each term Ji(t) for i=1,2,3. Hence,

    J1(t)=3j=12k=1[R3kb33bjkujdx+R3kb33ujkbjdx]CR3|hb||b||u|dx.

    Similarly, one obtains

    J2(t)=3j,k=12i=1[R3kbiibjkujdx+R3kbiiujkbjdx]CR3|b||hb||u|dx+CR3|hu||b|2dx

    and, by applying b=0, we have

    J3(t)=3j=1[R33b33bj3ujdx+R33b33uj3bjdx]=3j=12k=1[R3kbk3bj3ujdx+R3kbk3uj3bjdx]CR3|hb||b||u|dx.

    Replacing, in (21), the estimates obtained for Ji(t), i=1,2,3, we get

    (bb,Δu)2(bu,Δb)2CR3|hu||b|2dx+CR3|hb||b||u|dx. 

    Furthermore, notice that

    χ(×w,Δu)2χ(×u,Δw)2χw22+χΔu22, 

    where we have applied Cauchy-Schwarz's inequality. At last, Y. Zhou and M. Pokorný [29] proved that

    (uu,Δu)2CR3|hu||u|2dx. 

    Therefore, (19) reads

    12ddt(u,w,b)22+μΔu22+γΔw22+νΔb22+κ(w)22+χw22CR3|(hu,hw,hb)||(u,w,b)|2dx. (23)

    By using Lemma 2.1, one gets

    ddt(u,w,b)22+2α(Δu,Δw,Δb)22+2κ(w)22+2χw22C(hu,hw,hb)2(u,w,b)122(hu,hw,hb)2×(Δu,Δw,Δb)122, 

    where α=min{μ,γ,ν}. Now, by integrating over [Tτ,s], st, the inequality above yields

    (u,w,b)(s)22+2αsTτ(Δu,Δw,Δb)(τ)22dτC+CI(t)sTτ(u,w,b)(τ)122(hu,hw,hb)(τ)2×(Δu,Δw,Δb)(τ)122dτ, 

    where we applied the definition of I(t) given in (16). By Hölder's inequality, we have

    (u,w,b)(s)22+2αsTτ(Δu,Δw,Δb)(τ)22dτC+CI2(t)J12(t)(sTτ(u,w,b)(τ)22dτ)14, 

    for all st. In the inequality above the definitions of I(t) and J(t) were applied (see (16) and (17)). It follows, by using (17), that

    J2(t)C+CI2(t)J12(t)(tTτ(u,w,b)(τ)22dτ)14. 

    By using Young's inequality, we infer

    J2(t)C+CI83(t)(tTτ(u,w,b)(τ)22dτ)13+12J2(t). 

    Consequently, by applying (18), we obtain

    J(t)C+[C+CJ(t)](tTτ(u,w,b)(τ)22dτ)16. (24)

    From the energy inequality (9), one concludes that there exists 0<τ1 such that

    TTτ(u,w,b)(τ)22dτ1(2C)6.

    Now, we can obtain the desired estimate for (u,w,b)(t) in [Tτ,T). In fact, by replacing the bound above in (24), we get

    J(t)C, t[Tτ,T). 

    The definition (17) establishes the proof of Theorem 1.1.

    In order to prove Theorem 1.2 let us examine all the terms on the right hand side of (13) in an alternative way. We have

    (uw,Δhw)2=3i,j=12k=1R3uiiwj2kwjdx=3i,j=12k=1R3kuiiwjkwjdx=3i,j=12k=1R3wjkuikiwjdxCR3|w||u||hw|dx,

    since u=0. By arguing in the same way, one gets

    (ub,Δhb)2CR3|b||u||hb|dx.

    Notice that

    (bb,Δhu)2=3i,j=12k=1R3biibj2kujdxCR3|b||b||hu|dx

    and also

    (bu,Δhb)2=3i,j=12k=1R3biiuj2kbjdxCR3|b||u||hb|dx.

    The reader might check that (14) assures the following estimate:

    χ(×u,Δhw)2χ(×w,Δhu)2χhu22+χhw22. 

    At last, Y. Zhou and M. Pokorný [29] proved that

    (uu,Δhu)2CR3|u3||u||hu|dx. 

    By replacing all these last results obtained above in (13) and by using Young's inequality, one has

    12ddt(hu,hw,hb)22+α(hu,hw,hb)22+κh(w)22+χhw22CR3|(u3,w,b)|2|(u,w,b)|2dx+α2R3|(hu,hw,hb)|2dx, 

    where α=min{μ,γ,ν}. Hence,

    ddt(hu,hw,hb)22+α(hu,hw,hb)22+2κh(w)22+2χhw22CR3|(u3,w,b)|2|(u,w,b)|2dx. 

    By Lemmas 2.2 and 2.3, and also by (8), we obtain

    ddt(hu,hw,hb)22+α(hu,hw,hb)22+2κh(w)22+2χhw22C(hu,hw,hb)2(u,w,b)2×(hu,hw,hb)2. 

    By Young's inequality, one concludes

    ddt(hu,hw,hb)22+α2(hu,hw,hb)22+2κh(w)22+2χhw22C(hu,hw,hb)22(u,w,b)22. 

    By applying Gronwall's inequality, we get

    (hu,hw,hb)(t)2(hu,hw,hb)(δ)2×exp{CTδ(u,w,b)(s)22ds}, 

    for all t[δ,T). By energy inequality (9), we can guarantee the following estimate:

    (hu,hw,hb)(t)2C, t[δ,T). (25)

    In order to prove that the term (u,w,b)(t)2 is bounded in [δ,T), we recall that (23) and Hölder's inequality imply

    ddt(u,w,b)22+2α(Δu,Δw,Δb)22+2κ(w)22+2χw22C(hu,hw,hb)2(u,w,b)24, 

    where α=min{μ,γ,ν}. By using (25) and Gagliardo-Nirenberg's inequality, one has

    ddt(u,w,b)22+2α(Δu,Δw,Δb)22+2κ(w)22+2χw22C(u,w,b)122(Δu,Δw,Δb)322, 

    for all t[δ,T). By Young's inequality, we infer

    ddt(u,w,b)22+α(Δu,Δw,Δb)22+2κ(w)22+2χw22C(u,w,b)22. 

    By Gronwall's inequality,

    (u,w,b)(t)2C(u,w,b)(δ)2, t[δ,T).

    This completes the proof of Theorem 1.2.

    The authors would like to thank reviewers for their precious suggestions.



    [1] Arias LAS, Cirillo P, Oosterlee CW (2022) A new self-exciting jump-diffusion process for option pricing. arXiv preprint arXiv: 2205.13321. https://doi.org/10.48550/arXiv.2205.13321
    [2] Arregui I, Ráfales J (2020) A stochastic local volatility technique for TARN options. Int J Comput Math 97: 1133–1149. https://doi.org/10.1080/00207160.2019.1608357 doi: 10.1080/00207160.2019.1608357
    [3] Bandelier V (2017) Pricing FX-TARN under Lévy processes using numerical methods. Thesis, EPFL.
    [4] Carr P, Wu L (2004) Time-changed Lévy processes and option pricing. J Financ Econ 71: 113–141. https://doi.org/10.1016/S0304-405X(03)00171-5 doi: 10.1016/S0304-405X(03)00171-5
    [5] Caspers P (2015) Fast approximate pricing for FX target redemption forwards. SSRN. http://dx.doi.org/10.2139/ssrn.2606957 doi: 10.2139/ssrn.2606957
    [6] Fang F, Jönsson H, Oosterlee CW. et al. (2010) Fast valuation and calibration of credit default swaps under Lévy dynamics. J Comput Financ 14: 57–86. https://doi.org/10.21314/JCF.2010.209 doi: 10.21314/JCF.2010.209
    [7] Fang F, Oosterlee CW (2008) A novel pricing method for European options based on Fourier-cosine series expansions. SIAM J Sci Comput 31: 826–848. https://doi.org/10.1137/080718061 doi: 10.1137/080718061
    [8] Fang F, Oosterlee CW (2009) Pricing early-exercise and discrete barrier options by fourier-cosine series expansions. Numer Math 114: 27–62. https://doi.org/10.1007/s00211-009-0252-4 doi: 10.1007/s00211-009-0252-4
    [9] Fang F, Oosterlee CW (2011) A Fourier-based valuation method for Bermudan and barrier options under Heston's model. SIAM J Financ Math 2: 439–463. https://doi.org/10.1137/100794158 doi: 10.1137/100794158
    [10] Gohberg I, Olshevsky V (1994) Fast algorithms with preprocessing for matrix-vector multiplication. J Complex 10: 411–427. https://doi.org/10.1006/jcom.1994.1021 doi: 10.1006/jcom.1994.1021
    [11] Junike G, Pankrashkin K (2022) Precise option pricing by the COS method-How to choose the truncation range. Appl Math Comput 421: 1–14. https://doi.org/10.1016/j.amc.2022.126935 doi: 10.1016/j.amc.2022.126935
    [12] Lord R, Fang, F, Bervoets F, et al. (2008) A fast and accurate FFT-based method for pricing early-exercise options under Lévy processes. SIAM J Sci Comput 30: 1678–1705. https://doi.org/10.1137/070683878 doi: 10.1137/070683878
    [13] Luo X, Shevchenko P (2015) Pricing TARNs using a finite difference method. J Deriv 23: 62–72. https://doi.org/10.3905/jod.2015.23.1.062 doi: 10.3905/jod.2015.23.1.062
    [14] Ruijter MJ, Oosterlee CW (2012) Two-dimensional Fourier cosine series expansion method for pricing financial options. SIAM J Sci Comput 34:B642–B671. https://doi.org/10.1137/120862053 doi: 10.1137/120862053
    [15] Ruijter MJ, Versteegh M, Oosterlee CW (2015) On the application of spectral filters in a Fourier option pricing technique. J Comput Financ 19: 75–106. https://doi.org/10.21314/JCF.2015.306 doi: 10.21314/JCF.2015.306
    [16] Sato K (1999) Lévy processes and infinitely divisible distribution, Cambridge University Press, Cambridge.
    [17] Schoutens W (2003) Lévy Processes in Finance, Wiley, Chichester. https://doi.org/10.1002/0470870230
    [18] Zhang B, Oosterlee CW (2014) Pricing early-exercise Asian options under Lévy processes based on Fourier cosine expansions. Appl Numer Math 78: 14–30. https://doi.org/10.1007/s00211-009-0252-4 doi: 10.1007/s00211-009-0252-4
  • This article has been cited by:

    1. Diego Chamorro, David Llerena, Interior ϵ-regularity theory for the solutions of the magneto-micropolar equations with a perturbation term, 2022, 8, 2296-9020, 555, 10.1007/s41808-022-00163-y
    2. Muhammad Naqeeb, Amjad Hussain, Ahmad Mohammed Alghamdi, Blow-up criteria for different fluid models in anisotropic Lorentz spaces, 2023, 8, 2473-6988, 4700, 10.3934/math.2023232
    3. Diego Chamorro, David Llerena, A crypto-regularity result for the micropolar fluids equations, 2023, 520, 0022247X, 126922, 10.1016/j.jmaa.2022.126922
    4. Xinliang Li, Dandan Ding, Global well-posedness and large-time behavior of solutions to the 3D inviscid magneto-micropolar equations with damping, 2024, 103, 0003-6811, 1963, 10.1080/00036811.2023.2271946
    5. Diego Chamorro, David Llerena, Gastón Vergara-Hermosilla, Some remarks about the stationary micropolar fluid equations: Existence, regularity and uniqueness, 2024, 536, 0022247X, 128201, 10.1016/j.jmaa.2024.128201
    6. Takashi Kurihara, Possible involvement of FFAR1 signaling in mouse emotional behaviors through the regulation of brain monoamine releases, 2023, 158, 0015-5691, 454, 10.1254/fpj.23054
    7. Diego Chamorro, David Llerena, Partial suitable solutions for the micropolar equations and regularity properties, 2025, 31, 2118-7436, 137, 10.5802/ambp.428
    8. Diego Chamorro, David Llerena, Partial regularity and L3-norm concentration effects around possible blow-up points for the micropolar fluid equations, 2025, 32, 1021-9722, 10.1007/s00030-025-01059-1
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1722) PDF downloads(159) Cited by(0)

Figures and Tables

Figures(7)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog