Processing math: 100%
World Scientific
Skip main navigation

Cookies Notification

We use cookies on this site to enhance your user experience. By continuing to browse the site, you consent to the use of our cookies. Learn More
×

System Upgrade on Tue, May 28th, 2024 at 2am (EDT)

Existing users will be able to log into the site and access content. However, E-commerce and registration of new users may not be available for up to 12 hours.
For online purchase, please visit us again. Contact us at customercare@wspc.com for any enquiries.

Stochastic complex integrals in a two-dimensional flow

    https://doi.org/10.1142/S2661335221500027Cited by:1 (Source: Crossref)

    Abstract

    Two-dimensional flow is considered in the xy plane. A flat plate is placed parallel to the x-axis. The circulation of the flow is investigated and it is the real part of the stochastic complex integral. In that study, analogs of Karhunen–Loève expansions for stochastic complex integrals are studied.

    1. Introduction

    We consider two-dimensional flow of a liquid in the xy plane. Suppose every fluid particle moves with the constant speed U parallel to the x-axis. Place a flat plate on the xy plane with its leading edge coinciding with x=0, as shown in Fig. 1. For example, see the example on p. 80 in Chorin and Marsden’s book1 or Sec. 23.20 in Milne-Thomson’s book.2

    Fig. 1.

    Fig. 1. Velocity profile near the plate.

    The fluid velocity on the plate surface is the same as that of the plate, that is, it is equal to zero. On the other hand, the fluid velocity in the regions far from the plate is equal to U. The fluid velocity drastically changes in the neighborhood of the plate surface. This region is called a boundary layer. In Fig. 1, d denotes the thickness of the boundary layer. When the flow rate is small, the flow is laminar. Then we denote by u(x,y)=(u(x,y),v(x,y)) the velocity of the fluid particle that is moving through (x,y). Here, the flow is steady. When the flow rate increases, we suppose that a random force acts on the flow in the y-axis direction. Under this assumption, we investigate the circulation in the boundary layer. We adopt the following stochastic process {Zy} as a random force :

    Zy={Xy,for y0,˜X(y),for y<0.
    Here, {Xy:y0} and {˜Xy:y0} are independent Lévy processes on such that
    E[eiξXy]=eyΨ(ξ)andE[eiξ˜Xy]=eyΨ(ξ),
    where
    Ψ(ξ)=a02ξ2+iγξ+(eiξs1iξs1{|s|1}(s))Π(ds),
    a00, γ, min{1,|s|2}Π(ds)<, and Π({0})=0. Then we have the decomposition
    Zy=k=0(Zky+γky),
    where each {Zky:y} is a second-order stochastic process with mean zero and
    Zky={Xky,for y0,˜Xk(y),for y<0.
    Here, {Xky:y0} and {˜Xky:y0} are independent Lévy processes on such that
    E[eiξXky]=eyΨk(ξ)andE[eiξ˜Xky]=eyΨk(ξ),
    where
    Ψ0(ξ)=a02ξ2,
    and, for k1,
    Ψk(ξ)=(eiξs1iξs)Πk(ds),
    and
    Πk(ds)=1{k1<|s|k}(s)Π(ds).
    In addition, we have γ0=γ, γ1=0 and, for k2
    γk=k1<|s|ksΠ(ds).
    Here, we note that {Zky}, k=0,1,2,, are independent. Let 0<A< and let C be a closed curve in the boundary layer such that
    C(0,A)×(0,A).
    The circulation Γ around C is defined to be the line integral
    Γ=Cudt,
    where dt=(dx,dZy). Similarly, we define the flux Θ across C by the line integral
    Θ=Cudn,
    where dn=(dZy,dx). Then the circulation is represented as the real part of the stochastic complex integral
    Cw(z)dZ=Γ+iΘ,
    where dZ=dx+idZy and w=uiv. In particular, w is called a complex velocity. Let β>0. In addition, we consider a stochastic process {Wky:y} such that
    Wk(y)Wk(y0)=βyy0Wk(s)ds+Zk(y)Zk(y0),(1)
    for <y0<y<. This process is called an Ornstein–Uhlenbeck type process (OU type process). The solution of (1) is represented as
    Wky=yeβ(ys)dZks.(2)
    Each {Wky:y} is stationary in the sense that, for every s
    {Wky+s:y}d={Wky:y}.
    See Refs. 3 and 4 for more details. We define random variables
    Φkn=AAWk(y)ϕn(y)dy,
    for n and k{0}, where the symbol ϕn(y) is found in Proposition 3.3. By virtue of Theorem 6.3 in Sec. 6, we can obtain the circulation in the boundary layer.

    Theorem 1.1. Let the length of  C  be finite. If  u  is of class  C1, the circulation  Γ  is given by

    Γ=Cu(x,y)dx+k=0(γkCv(x,y)dy+Cv(x,y)dWky+βCv(x,y)Wk(y)dy)=Cu(x,y)dx+k=0(γkCv(x,y)dy+l.i.m.NNn=NΦkn{Cϕn(y)v(x,y)xdx+Cϕn(y)(βv(x,y)v(x,y)y)dy}).
    Here, l.i.m. stands for limit in quadratic mean.

    Remark 1.1. Let ν>0. Suppose the pressure gradient dpdx is equal to zero, and consider the boundary layer equation

    uux+vuy=ν2uy2,(3)
    and the continuity equation
    ux+vy=0,(4)
    where u(x,0)=v(x,0)=0 for x0 and u=U as y. We see from (4) that there is a stream function ψ such that
    u=ψyandv=ψx.
    With the change of variables
    η=yUνx,
    the stream function can be expressed as
    ψ(x,η)=νUxf(η).
    Hence, (3) implies the Blasius equation
    ff+2f=0.
    The boundary conditions are that f=f=0 for η=0 and f=1 as η. Then it follows that
    u(x,y)=Uf(η),v(x,y)=12νUx(ηf(η)f(η)).
    See Ref. 5 for the Blasius problem.

    Stochastic integrals of nonrandom functions with respect to additive processes on were studied in Refs. 6,7, and 8. See also Refs. 911. Sato developed this theory in Refs. 12 and 13. These days, a lot of papers related to this theory exist, for example, Refs. 1418. Stochastic complex integrals are constructed on basis of this theory in Ref. 19. Our study extends to the application to Blasius’ formula of fluid mechanics in Ref. 20. The paper is written to be as self-contained as possible. The terminologies follow Ref. 21.

    2. Covariance Functions

    We find the covariance functions, before we define the stochastic integrals with respect to {Wky}.

    Lemma 2.1. The law of  Wky  is represented as 

    Eexp(iξWky)=exp[yΨk(ξeβ(ys))ds].

    Lemma 2.1 is derived from Proposition 2.6 of Ref. 7. Lemma 2.1 tells us that the mean of Wky is equal to 0, that is

    E[Wky]=0.
    Define the covariance function of Wky by
    Rk(s,y)=E[(WksE[Wks])(WkyE[Wky])]=E[WksWky].

    Lemma 2.2. The second moment of  Wky  is represented as 

    E[(Wky)2]=12βΨk(0).

    Proof. From Lemma 2.1, we obtain that

    E[(Wky)2exp(iξWky)]=(ye2β(ys)Ψk(ξeβ(ys))ds+(yeβ(ys)Ψk(ξeβ(ys))ds)2)×exp[yΨk(ξeβ(ys))ds].
    Substituting 0 for ξ, we get the lemma. □

    Lemma 2.3. For  h>0, it holds that 

    E[h0eβ(hs)dZks]=0
      and 
    E[(h0eβ(hs)dZks)2]=1e2βh2βΨk(0).

    Proof. Since we have

    Eexp(iξh0eβ(hs)dZks)=exp[h0Ψk(ξeβ(hs))ds],
    the lemma is obtained. □

    Lemma 2.4. For  h>0, it holds that 

    E[(Wk0+Wkh)2]=1+eβhβΨk(0).

    Proof. We have

    Wk0+Wkh=0eβsdZks+heβ(hs)dZks=(1+eβh)0eβsdZks+h0eβ(hs)dZks.
    From Lemmas 2.2 and 2.3, we obtain that
    E[(Wk0+Wkh)2]=(1+eβh)2E[(0eβsdZks)2]+2(1+eβh)E[0eβsdZks]×E[h0eβ(hs)dZks]+E[(h0eβ(hs)dZks)2]=(1+eβh)2(12βΨk(0))1e2βh2βΨk(0)=1+eβhβΨk(0).
    The lemma has been proved. □

    We give the covariance function of Wkt.

    Theorem 2.1. The covariance function of  Wky  is represented as 

    Rk(s,y)=Υkeβ|sy|,
      where 
    Υk=Ψk(0)2β.

    Proof. Let h=|sy|. From Lemmas 2.2 and 2.4, we obtain that

    E[Wk0Wkh]=21(E[(Wk0+Wkh)2]E[(Wk0)2]E[(Wkh)2])=21(E[(Wk0+Wkh)2]2E[(Wk0)2])=21(1+eβhβΨk(0)+1βΨk(0))=eβh2βΨk(0).
    Hence, we have
    Rk(s,y)=Rk(0,h)=eβh2βΨk(0).
    The theorem has been proved. □

    3. Karhunen–Loève Expansions

    Let <Aa<bA<. We first discuss the existence of the integral

    bag(y)Wk(y)dy.
    Subdivide [a,b] as
    a=yn0<yn1<yn2<<ynkn=b,
    and denote this partition by
    Δn={yn0,yn1,,ynkn},
    where we take Δn such that
    ΔnΔn+1.
    We define the size of Δn by
    |Δn|=max{yniyni1:i=1,2,,kn}.
    Then for the partition Δn, we define
    Q(g,Δn)=kni=1g(ni)Wk(ni)(yniyni1),
    where ni[yni1,yni] for each i.

    Definition 3.1. If the quadratic mean limit of Q(g,Δn) exists as |Δn|0 and is independent of the choice of {Δn} and, for each Δn, is independent of the choice of {ni}, then the limit is denoted by

    bag(y)Wk(y)dy.

    Theorem 3.1. If  g(y)  is continuous on  [a,b], then the integral  bag(y)Wk(y)dy  exists as a quadratic mean limit.

    Proof. We have

    E[Q(g,Δn)Q(g,Δm)]=kni=1kmj=1g(ni)g(mj)×Rk(ni,mj)(yniyni1)(ymjymj1).
    Since Rk(s,y) is continuous, this limit exists as n,m and is independent of the choice of {Δn}. Hence, we obtain that
    E[(Q(g,Δn)Q(g,Δm))2]=E[Q(g,Δn)2]+E[Q(g,Δm)2]2E[Q(g,Δn)Q(g,Δm)]0,
    as n,m. The theorem has been proved. □

    Assume that Υk0. We next consider the integral equation

    AARk(s,y)η(s)ds=λη(y),AyA,(5)
    where η(y) satisfies
    AA|η(y)|2dy<.
    Then λ and η(y) are called an eigenvalue and the corresponding eigenfunction. We see from (5) that η(y) is continuous. Hence, we have
    λAAη(y)2dy=AAAARk(s,y)η(s)η(y)dsdy=E(AAWk(y)η(y)dy)20.
    This implies λ0. Theorem 2.1 tells us that
    ΥkAAeβ|sy|η(s)ds=λη(y).
    This integral equation is equal to
    TTe|uv|f(u)du=μf(v),(6)
    where T=βA, u=βs, v=βy, μ=Υ1kβλ, and f(u)=η(uβ1). Equation (6) can be written
    μf(v)=vTeuvf(u)du+Tveu+vf(u)du.
    Differentiate both sides with respect to v twice
    μf(v)=vTeuvf(u)du+Tvevuf(u)du,μf(v)=TTe|uv|f(u)du2f(v)=μf(v)2f(v).
    Hence
    μf(v)+(2μ)f(v)=0,(7)
    where μ0. If μ=0 or if μ=2, then f(v)=0.

    Lemma 3.1. The integral equation (6is not satisfied for any  μ>2.

    Proof. Equation (7) has the general solution

    f(v)=c1eαv+c2eαv,
    where c1,c2 and
    0<α=μ2μ<1.
    Substituting this expression in (6), we obtain that
    eαv[c1α+1+c11α]+eαv[c2α+1+c21α]+ev[c1e(α+1)Tα+1c2e(α1)T1α]+ev[c2e(α+1)Tα+1c1e(α1)T1α]=μ(c1eαv+c2eαv).
    As the coefficients of ev and ev are equal to zero, it follows that
    c1e(α+1)Tα+1+c2e(α1)T1α=0,c2e(α+1)Tα+1+c1e(α1)T1α=0.(8)
    Adding these equations gives
    (c1+c2)(e(α+1)Tα+1+e(α1)T1α)=0.
    Hence, c1+c2=0. It follows from (8) that if c10 and c20
    1α1+α=e2αT.
    The left-hand side is smaller than 1, but the right-hand side is bigger than 1. This is a contradict. Hence, in the case μ>2, (6) does not have a nontrivial solution. □

    Hence, we see from Lemma 3.1 that if (6) has a nontrivial solution, then 0<μ<2. We consider the equation

    f(v)+α2f(v)=0,(9)
    where
    α=2μμ>0.

    Proposition 3.1. Eigenvalues for (5are 

    λkn=2Υkβ(1+α2n),
      and the corresponding eigenfunctions are 
    ηn(y)=cos(αnβy),
      if  n0, and 
    ηn(y)=sin(αnβy),
      if  n1. Here, if  n0, then  αn  satisfies 
    αntan(αnT)=1and0<αn<αn+1,
      and if  n1, then  αn  satisfies 
    tan(αnT)=αnand0<αn<αn1.

    Proof. Now we have

    μf(T)=TTeuTf(u)du,μf(T)=TTeuTf(u)du,μf(T)=TTeTuf(u)du,μf(T)=TTeTuf(u)du.
    Hence, the boundary condition is that
    f(T)+f(T)=0,(10)
    f(T)f(T)=0.(11)
    Equation (9) has the general solution
    f(v)=c3sin(αv)+c4cos(αv),
    where c3,c4. The boundary condition yields that
    (f(T)+f(T)f(T)f(T))=(c3c4αc3α+c4c3c4αc3α+c4)(sin(αT)cos(αT))=(00).
    Then we have
    |c3c4αc3α+c4c3c4αc3α+c4|=2c3c4(1+α2).
    If c3c40, then
    sin(αT)=0,cos(αT)=0.
    This is a contradiction. If c3=0 and c40, then
    c4αsin(αT)+c4cos(αT)=0,
    that is
    αtan(αT)=1.
    If c30 and c4=0, then
    c3sin(αT)+c3αcos(αT)=0,
    that is
    tan(αT)=α.
    Recall that
    μ=21+α2.
    Hence, in the case where αtan(αT)=1, we have
    λ=2Υkβ(1+α2)andη(y)=cos(αβy).
    In the case where tan(αT)=α, we have
    λ=2Υkβ(1+α2)andη(y)=sin(αβy).
    The proposition has been proved. □

    Here, we show the range of each eigenvalue.

    Proposition 3.2. If  n0, then 

    2T2Υkβ(T2+(n+21)2π2)<λkn<2T2Υkβ(T2+n2π2).
      If  n1, then 
    2T2Υkβ(T2+n2π2)<λkn<2T2Υkβ(T2+(n+21)2π2).

    Proof. If n0, then

    tan(αnT)=TαnT.
    This implies
    nπ<αnT<(n+21)π.
    If n1, then
    tan(αnT)=1TαnT.
    This implies
    (n+21)π<αnT<nπ.
    It follows from Proposition 3.1 that the proposition holds. □

    Now we find the orthonormal eigenfunctions required for the Karhunen–Loève expansion.

    Proposition 3.3. Let 

    ϕn(y)=2αnβ2αnβA+ϵnsin(2αnβA)ηn(y),
      where  ϵn=1  if  n0  and  ϵn=1  if  n<0. The eigenfunctions  {ϕn(y):n}  are orthonormal.

    Proof. Let fn(u)=ηn(β1u). Let nm. We see from (9) that

    fn(u)+α2nfn(u)=0,fm(u)+α2mfm(u)=0.
    These equations yield
    (α2mα2n)fn(u)fm(u)=fm(u)fn(u)fn(u)fm(u).
    Hence, we obtain that
    AAηn(y)ηm(y)dy=β1TTfn(u)fm(u)du=1(α2mα2n)βTT(fm(u)fn(u)fn(u)fm(u))du=1(α2mα2n)β(fm(T)fn(T)fn(T)fm(T)fm(T)fn(T)+fn(T)fm(T)).
    Using (10) and (11), we see that
    (fm(u)+fm(u))(fn(u)fn(u))=0,(fn(u)+fn(u))(fm(u)fm(u))=0
    for u=T,T. Subtracting both sides, we see that
    fn(u)fm(u)fm(u)fn(u)=0,
    for u=T,T. Hence
    AAηn(y)ηm(y)dy=0,
    and thereby {ηn(y):n} are orthogonal. Furthermore, we have
    AAηn(y)2dy=2αnβA+ϵnsin(2αnβA)2αnβ.
    The proposition has been proved. □

    Let δn,m=1 if n=m and δn,m=0 if nm.

    Proposition 3.4. The random variables  {Φkn}  are orthogonal, that is

    E[ΦknΦlm]=λknδk,lδn,m.

    Proof. We have

    E[ΦknΦlm]=AAAAE[Wk(y)Wl(s)]ϕn(y)ϕm(s)dyds.
    If kl, then Wk(y) and Wl(s) are independent. Hence
    E[Wk(y)Wl(s)]=E[Wk(y)]E[Wl(s)]=0.
    This implies E[ΦknΦlm]=0.

    Let k=l. We have

    E[ΦknΦkm]=AA(AARk(s,y)ϕn(y)dy)ϕm(s)ds=λknAAϕn(s)ϕm(s)ds=λknδn,m.
    The proposition has been proved. □

    The Karhunen–Loève expansion is represented as follows: see Refs. 22 and 23 for the Karhunen–Loève expansion.

    Theorem 3.2. We have 

    Wk(y,ω)=l.i.m.NNn=NΦkn(ω)ϕn(y),
      uniformly for  AyA, and 
    E[ΦknΦkm]=λknδn,m.

    Theorem 3.2 tells us a proposition.

    Proposition 3.5. For any  m, we have 

    E[WkyΦkm]=λkmϕm(y).

    Proof. As N, we see that

    |E[Nn=Nϕn(y)ΦknΦkmWkyΦkm]|2E[(Nn=Nϕn(y)ΦknWky)2]E|Φkm|20.
    Theorem 3.2 yields that
    E[WkyΦkm]=limNE[Nn=Nϕn(y)ΦknΦkm]=ϕm(y)E[ΦkmΦkm]=ϕm(y)λkm.
    The proposition has been proved. □

    4. Stochastic Integrals with Respect to Momentum

    Let <Aa<bA<. Subdivide [a,b] as

    a=yn0<yn1<yn2<<ynkn=b,
    and denote this partition by
    Δn={yn0,yn1,,ynkn},
    where we take Δn such that
    ΔnΔn+1.
    For the partition Δn, we define
    S(g,Δn)=kni=1g(ni)(Wk(yni)Wk(yni1)),
    where ni[yni1,yni] for each i.

    Definition 4.1. If the quadratic mean limit of S(g,Δn) exists as |Δn|0 and is independent of the choice of {Δn} and, for each Δn, is independent of the choice of {ni}, then the limit is denoted by

    bag(y)dWky,
    and g is said to be {Wky}-integrable.

    We introduce two half planes

    H+={(x,y):xy}andH={(x,y):xy}.
    Define
    Dk(I)=Rk(t,v)Rk(s,v)Rk(t,u)+Rk(s,u)=Υk(eβ|tv|eβ|sv|eβ|tu|+eβ|su|),
    for any rectangle I=[s,t]×[u,v]. Let
    Ini,j=[yni1,yni]×[ynj1,ynj],
    for 1i,jkn. Denote by |Ini,j| the length of the diagonal of Ini,j.

    Lemma 4.1. (i) If  ij, then we have 

    Dk(Ini,j)=Ψk(0)β2(yniyni1)(ynjynj1)eβ|yni1ynj1|+(yniyni1)(ynjynj1)l(Ini,j).
      Here, l(Ini,j)uniformly converges to 0 as|Ini,j|0.

    (ii) If  i=j, then we have 

    Dk(Ini,i)=Ψk(0)(yniyni1)eβ(yniyni1)+(yniyni1)2L(Ini,i).
      Here, L(Ini,i)is uniformly bounded as|Ini,i|0.

    Proof. We first prove (i). We consider the case Ini,jH+. Then

    Dk(Ini,j)=Υk(eβ(yniynj)eβ(yni1ynj)eβ(yniynj1)+eβ(yni1ynj1)).
    Let yniyni1ynjynj1. There exists
    ξ1(yni1ynj,yniynj),
    such that
    eβ(yniynj)eβ(yni1ynj)=β(yniyni1)eβ(yni1ynj)+β2(yniyni1)22eβξ1,
    and exists
    ξ2(yni1ynj1,yniynj1),
    such that
    eβ(yniynj1)eβ(yni1ynj1)=β(yniyni1)eβ(yni1ynj1)+β2(yniyni1)22eβξ2.
    In addition, we have
    eβ(yni1ynj)eβ(yni1ynj1)=β(ynjynj1)eβ(yni1ynj1)+β2(ynjynj1)22eβξ3,
    for some ξ3(yni1ynj,yni1ynj1). Hence, we obtain that
    Dk(Ini,j)Υk=β(yniyni1)(eβ(yni1ynj)eβ(yni1ynj1))+β2(yniyni1)22(eβξ1eβξ2)=β2(yniyni1)(ynjynj1)eβ(yni1ynj1)+β2(yniyni1)22(eβξ1eβξ2)β32(yniyni1)(ynjynj1)2eβξ3.
    Next, let yniyni1>ynjynj1. There exists
    ξ4(yniynj,yniynj1),
    such that
    eβ(yniynj)eβ(yniynj1)=β(ynjynj1)eβ(yniynj1)+β2(ynjynj1)22eβξ4,
    and exists
    ξ5(yni1ynj,yni1ynj1),
    such that
    eβ(yni1ynj)eβ(yni1ynj1)=β(ynjynj1)eβ(yni1ynj1)+β2(ynjynj1)22eβξ5.
    In addition, we have
    eβ(yniynj1)eβ(yni1ynj1)=β(yniyni1)eβ(yni1ynj1)+β2(yniyni1)22eβξ6,
    for some ξ6(yni1ynj1,yniynj1). Hence, we obtain that
    Dk(Ini,j)Υk=β(ynjynj1)(eβ(yniynj1)eβ(yni1ynj1))+β2(ynjynj1)22(eβξ4eβξ5)=β2(yniyni1)(ynjynj1)eβ(yni1ynj1)+β2(ynjynj1)22(eβξ4eβξ5)+β32(ynjynj1)(yniyni1)2eβξ6.
    This implies that (i) holds.

    We next consider the case Ini,jH. Then we have

    Dk(Ini,j)=Υk(eβ(ynjyni)eβ(ynjyni1)eβ(ynj1yni)+eβ(ynj1yni1)).
    Hence, exchanging i for j in the case Ini,jH+, we have
    Dk(Ini,j)Υk=β2(yniyni1)(ynjynj1)eβ(ynj1yni1)+(yniyni1)(ynjynj1)l(Ini,j).
    Assertion (i) remains true.

    Finally, we prove (ii). We consider the case i=j. Then we have

    Dk(Ini,i)Υk=2(eβ0eβ(yniyni1))=2β(yniyni1)eβ(yniyni1)+β2(yniyni1)2eβξ7,
    for some ξ7(0,yniyni1). The lemma has been proved. □

    Lemma 4.2. Then  1i,jkn|Dk(Ini,j)|  is uniformly bounded regardless of how  [a,b]×[a,b]  is divided.

    Proof. Let ij. From Lemma 4.1(i), we see that there is M1>0 such that

    |Dk(Ini,j)||Ψk(0)|β2(yniyni1)(ynjynj1)+M1(yniyni1)(ynjynj1).
    Let i=j. From Lemma 4.1(ii), we see that there is M2>0 such that
    |Dk(Ini,i)||Ψk(0)|(yniyni1)+M2(yniyni1)2(|Ψk(0)|+M2(ba))(yniyni1).
    Consequently, we have
    1i,jkn|Dk(Ini,j)|(|Ψk(0)|β2+M1)ij(yniyni1)(ynjynj1)+(|Ψk(0)|+M2(ba))kni=1(yniyni1)(|Ψk(0)|β2+M1)(ba)2+(|Ψk(0)|+M2(ba))(ba).
    Hence, the lemma is true. □

    Through this paper, we use “piecewise continuous” as the following meaning:

    Definition 4.2. Let a1,,an be discontinuous points of g(y) on (a,b). For each interval [ai1,ai], we define a function ği(y) by

    ği(y)={limtai1g(y),(y=ai1),g(y),(ai1<y<ai),limtaig(y),(y=ai).
    Here, a0=a and an+1=b. If ği(y) is continuous on [ai1,ai] for i=1,2,,n+1, g(y) is said to be piecewise continuous.

    Theorem 4.1. If  g(y)  is piecewise continuous on  [a,b], then the integral  bag(y)dWky  exists as a quadratic mean limit.

    Proof. Let n>m. Since ΔnΔm, we have

    S(g,Δn)S(g,Δm)=kni=1gn,m(i)(Wk(yni)Wk(yni1)),
    where gn,m(i)=g(ni)g(mj) if [yni1,yni][ymj1,ymj]. Let a1,,an be discontinuous points of g(y) on (a,b) and in addition, a0=a and an+1=b. Since ği(y) is uniformly continuous on [ai1,ai], it follows from Lemma 4.2 that
    limn,mE[(S(g,Δn)S(g,Δm))2]=limn,mkni=1knj=1gn,m(i)gn,m(j)(Rk(yni,ynj)Rk(yni,ynj1)Rk(yni1,ynj)+Rk(yni1,ynj1))=limn,mkni=1knj=1gn,m(i)gn,m(j)Dk(Ini,j)=0.
    Here, we used the fact that the number of discontinuous points of g is finite. The theorem has been proved. □

    The following theorem is obvious. The proof is omitted.

    Theorem 4.2. If  f  and  g  are  {Wky}-integrable, then  c1f+c2g  is  {Wky}-integrable and 

    ba(c1f(y)+c2g(y))dWky=c1baf(y)dWky+c2bag(y)dWkya.s.,
      for any  c1,c2.

    We show one lemma to prove Proposition 4.1. The proof is omitted.

    Lemma 4.3. Let  Xn  and  X  be random variables and suppose  supnEX2n<. If  Xn  converges to  X  in quadratic mean, then 

    limnEX2n=EX2.

    Proposition 4.1. Let  [a,b][aj,bj]  for  j=1,2. If  f(y)  and  g(y)  are piecewise continuous on  [a1,b1]  and on  [a2,b2], respectively, then we have 

    E[b2a2g(y)dWky]=0,
      and 
    E[b1a1f(y)dWkyb2a2g(y)dWky]=Ψk(0)β2b1a1b2a2f(s)g(y)eβ|sy|dydsΨk(0)[a1,b1][a2,b2]f(y)g(y)dy.

    Proof. We have

    E[S(g,Δn)]=kni=1g(ni)(E[Wk(yni)]E[Wk(yni1)])=0.
    Hence, it follows that
    E[bag(y)dWky]=lim|Δn|0E[S(g,Δn)]=0.
    Let h(y) be piecewise continuous on [a,b]. From Lemma 4.1, we obtain that
    E[S(h,Δn)2]=kni=1knj=1h(ni)h(nj)Dk(Ini,j)=Ψk(0)β2kni=1knj=1h(ni)h(nj)eβ|yni1ynj1|×(yniyni1)(ynjynj1)Ψk(0)β2kni=1h(ni)2(yniyni1)2+ijh(ni)h(nj)l(Ini,j)(yniyni1)(ynjynj1)Ψk(0)kni=1h(ni)2(yniyni1)eβ(yniyni1)+kni=1h(ni)2L(Ini,i)(yniyni1)2.
    Hence, it follows from Lemmas 4.2 and 4.3 that
    E[(bah(y)dWky)2]=lim|Δn|0E[S(h,Δn)2]=Ψk(0)β2babah(s)h(y)eβ|sy|dydsΨk(0)bah(y)2dy.
    Here, we define f0 and g0 by
    f0(y)={f(y),y[a1,b1],0,y[a,b][a1,b1],
    and
    g0(y)={g(y),y[a2,b2],0,y[a,b][a2,b2],
    respectively. Theorem 4.2 tells us that
    E[(baf0(y)dWky+bag0(y)dWky)2]=E[(ba(f0(y)+g0(y))dWky)2]=Ψk(0)β2baba(f0(s)+g0(s))(f0(y)+g0(y))×eβ|sy|dydsΨk(0)ba(f0(y)+g0(y))2dy.
    Hence, we obtain that
    E[baf0(y)dWkybag0(y)dWky]=12(E[(baf0(y)dWky+bag0(y)dWky)2]E[(baf0(y)dWky)2]E[(bag0(y)dWky)2])=Ψk(0)β2babaf0(s)g0(y)eβ|sy|dydsΨk(0)baf0(y)g0(y)dy.
    The proposition has been proved. □

    Theorem 4.3. If  g(y)  is of class  C1, then 

    bag(y)dWky=g(b)Wk(b)g(a)Wk(a)bag(y)Wk(y)dy,
      almost surely.

    Proof. From Theorem 4.1, we see that g is {Wky}-integrable. For the partition Δn, we take

    S(g,Δn)=kni=1g(yni)(Wk(yni)Wk(yni1)).
    Then we see that
    S(g,Δn)=g(b)Wk(b)g(yn1)Wk(a)kni=2(g(yni)g(yni1))Wk(yni1)=g(b)Wk(b)g(yn1)Wk(a)kni=2(g(yni1)+ϵi)(yniyni1)Wk(yni1).
    Here, ϵi0 uniformly in i as |Δn|0. Hence, the theorem holds. □

    5. Representation Through OU Type Processes

    Stochastic integrals based on {Zy} are represented in terms of Wky.

    Theorem 5.1. Let  <a<b<. If  g(y)  is a function of class  C1, then 

    bag(y)dZy=k=0(bag(y)dWky+βbag(y)Wk(y)dy+γkbag(y)dy)a.s..

    In order to prove Theorem 5.1, we prepare a lemma.

    Lemma 5.1. For any function  g(y)  of class  C1, we have 

    l.i.m.nkni=1g(yni1)yniyni1(Wk(y)Wk(yni1))dy=0.

    Proof. We have

    E(kni=1g(yni1)yniyni1(Wk(y)Wk(yni1))dy)2=i,jg(yni1)g(ynj1)×E[yniyni1(Wk(y)Wk(yni1))dy×ynjynj1(Wk(s)Wk(ynj1))ds]=i,jg(yni1)g(ynj1)yniyni1ynjynj1(Rk(y,s)Rk(yni1,s)Rk(y,ynj1)+Rk(yni1,ynj1))dsdyK.
    Let ϵ>0. As we take |Ini,j| small enough, we have
    |Rk(y,s)Rk(yni1,s)Rk(y,ynj1)+Rk(yni1,ynj1)|<ϵ
    in (y,s)Ini,j=[yni1,yni]×[ynj1,ynj]. Then it follows that
    Kϵi,j|g(yni1)g(ynj1)|(yniyni1)(ynjynj1)ϵ(ba|g(y)|dy)2as |Δn|0.
    Since ϵ can be taken arbitrarily, K goes to zero as n. □

    Now we prove Theorem 5.1.

    Proof. Proof of Theorem 5.1

    Since {Zky}, k=0,1,2,, are independent, the integrals

    bag(y)dZky,k=0,1,2,,
    are independent. By virtue of Lévy’s theorem, we have
    bag(y)dZy=k=0(bag(y)dZky+γkbag(y)dy),
    almost surely. Now we have
    bag(y)dZky=p-limnkni=1g(yni1)(ZkyniZkyni1).
    Here, p-lim stands for limit in probability. On the other hand, we see from (1) that
    kni=1g(yni1)(ZkyniZkyni1)=kni=1g(yni1)(Wk(yni)Wk(yni1)+βyniyni1Wk(y)dy)=kni=1g(yni1)(Wk(yni)Wk(yni1))+βkni=1g(yni1)Wk(yni1)(yniyni1)+βkni=1g(yni1)yniyni1(Wk(y)Wk(yni1))dy.
    Theorems 3.1 and 4.1 and Lemma 5.1 tell us that almost surely
    bag(y)dZky=bag(y)dWky+βbag(y)Wk(y)dy.
    The theorem has been proved. □

    Theorem 5.1 is also represented as follows:

    Theorem 5.2. Let  <Aa<bA<. If  g(y)  is a function of class  C1, then 

    bag(y)dZy=k=0(l.i.m.NNn=NΦkn{g(b)ϕn(b)g(a)ϕn(a)+ba(βg(y)g(y))ϕn(y)dy}+γkbag(y)dy),
      almost surely.

    Proof. From Theorems 3.2 and 4.3, we obtain that

    bag(y)dWky+βbag(y)Wk(y)dy=g(b)Wk(b)g(a)Wk(a)+ba(βg(y)g(y))Wk(y)dy=l.i.m.NNn=NΦkn(g(b)ϕn(b)g(a)ϕn(a)+ba(βg(y)g(y))ϕn(y)dy),
    almost surely. The theorem has been proved. □

    6. Stochastic Complex Integrals

    In this section, z is any complex number and represented as z=x+iy, where x and y are real numbers. Furthermore, z=x+iy is identified with (x,y)2. Hence, we often use the representation g(x,y) instead of any complex function g(z). Only in this section, g is a complex-valued function. A curve Λ is represented as a function

    r:[0,1]2,
    which is of class C1. Let Λ be a curve with an initial point z0=(a,c) and a terminal point z1=(b,d) such that
    Λ˜D[A,A]×[A,A],
    where ˜D is some open set. From now on, we write
    r(t)=z(t)=(x(t),y(t)),
    so we have z0=z(0)=(a,c) and z1=z(1)=(b,d). We subdivide [0,1] as
    0=tm0<tm1<tm2<<tmkm=1
    and denote this partition by
    Δm={tm0,tm1,,tmkm},
    where we suppose Δm+1Δm. We define
    Δy(tmi)=y(tmi)y(tmi1)
    and
    Q(g,Δm,Λ)=kmi=1g(z(tmi1))Wk(y(tmi1))Δy(tmi).

    Definition 6.1. If Q(g,Δm,Λ) converges in quadratic mean as |Δm|0, and if the limit does not depend on the choice of the sequence {Δm}, then the limit is denoted by

    Λg(x,y)Wk(y)dy.

    To show quadratic mean convergence of stochastic integrals, the following criterion is useful. See Ref. 22.

    Lemma 6.1. The Loève criterion

    The random variable sequence  Xn  converges in quadratic mean if and only if  E[XmˉXn]  has a finite limit, when  mand  n  tend to infinity independently of each other.

    Remark 6.1. For any z, we denote by ˉz the complex-conjugate of z. Hence, ˉXn means the complex-conjugate of Xn.

    The following lemma is known as the Mercer’s theorem. See Refs. 24 and 25.

    Lemma 6.2. In our setting, we have 

    limNsupAs,yA|Rk(s,y)Nn=Nλknϕn(s)ϕn(y)|=0.

    For any random variable X, the norm of X is defined by

    X=E[XˉX].

    Lemma 6.3. Suppose the length of  Λ  is finite. If  g  is continuous on  ˜D, then 

    limmQ(g,Δm,Λ)Nn=NΦknΛg(z)ϕn(y)dyϵ(N),
      where  ϵ(N)  converges to 0 as  N.

    Proof. We obtain from Proposition 3.5 that

    Q(g,Δm,Λ)Nn=NΦknΛg(z)ϕn(y)dy2=kmi=1kmj=1g(z(tmi1))¯g(z(tmj1))Rk(y(tmi1),y(tmj1))×Δy(tmi)Δy(tmj)kmi=1g(z(tmi1))Nn=NE[Wk(y(tmi1))Φkn]×Λ¯g(z)ϕn(y)dyΔy(tmi)kmj=1¯g(z(tmj1))Nl=NE[Wk(y(tmj1))Φkl]×Λg(z)ϕl(y)dyΔy(tmj)+Nn=NNl=NE[ΦknΦkl]×Λ¯g(z)ϕn(y)dyΛg(z)ϕl(y)dy=kmi=1kmj=1g(z(tmi1))¯g(z(tmj1))Rk(y(tmi1),y(tmj1))×Δy(tmi)Δy(tmj)kmi=1g(z(tmi1))Nn=Nλknϕn(y(tmi1))×Λ¯g(z)ϕn(y)dyΔy(tmi)kmj=1¯g(z(tmj1))Nl=Nλklϕl(y(tmj1))×Λg(z)ϕl(y)dyΔy(tmj)+Nn=NλknΛ¯g(z)ϕn(y)dyΛg(z)ϕn(y)dy=kmi=1kmj=1g(z(tmi1))¯g(z(tmj1))Δy(tmi)Δy(tmj)×Nn=Nλknϕn(y(tmi1))ϕn(y(tmj1))+ϵ(Δm,N)Nn=Nλknkmi=1g(z(tmi1))ϕn(y(tmi1))Δy(tmi)×Λ¯g(z)ϕn(y)dyNl=Nλklkmj=1¯g(z(tmj1))ϕl(y(tmj1))Δy(tmj)×Λg(z)ϕl(y)dy+Nn=NλknΛ¯g(z)ϕn(y)dyΛg(z)ϕn(y)dy=Nn=Nλkn|kmi=1g(z(tmi1))ϕn(y(tmi1))Δy(tmi)Λg(z)ϕn(y)dy|2+ϵ(Δm,N).
    Here
    ϵ(Δm,N)=kmi=1kmj=1g(z(tmi1))¯g(z(tmj1))×(Rk(y(tmi1),y(tmj1))Nn=Nλknϕn(y(tmi1))ϕn(y(tmj1)))×Δy(tmi)Δy(tmj).
    Since the length of Λ is finite, Lemma 6.2 tells us that ϵ(Δm,N) goes to 0 uniformly in m. The lemma has been proved. □

    Lemma 6.4. Suppose the length of  Λ  is finite. If  g  is continuous on  ˜D, then 

    n=λkn|Λg(z)ϕn(y)dy|210|g(z(t))|2|dydt|dt10Rk(y,y)|dydt|dt<.

    Proof. We have

    |Λg(z)ϕn(y)dy|210|g(z(t))|2|dydt|dt10ϕn(y)2|dydt|dt.
    Hence, we obtain that
    n=λkn|Λg(z)ϕn(y)dy|210|g(z(t))|2|dydt|dt10n=λknϕn(y)2|dydt|dt10|g(z(t))|2|dydt|dt10Rk(y,y)|dydt|dt.
    The lemma is true. □

    Theorem 6.1. Suppose the length of  Λ  is finite. If  g  is continuous on  ˜D, then the integral  Λg(z)Wk(y)dy  exists as a quadratic mean limit, which is represented as 

    Λg(z)Wk(y)dy=l.i.m.NNn=NΦknΛg(z)ϕn(y)dy.

    Proof. From Lemma 6.2, we see that

    E[Q(g,Δm,Λ)Q(,Δ˜m,Λ)]=kmi=1k˜mj=1g(z(tmi1))¯g(z(t˜mj1))×Rk(y(tmi1),y(t˜mj1))Δy(tmi)Δy(t˜mj)=limNNn=Nλknkmi=1g(z(tmi1))ϕn(y(tmi1))Δy(tmi)×k˜mj=1¯g(z(t˜mj1))ϕn(y(t˜mj1))Δy(t˜mj).
    By using uniform convergence in m and ˜m, it follows from Lemma 6.4 that the limit is equal to
    n=λkn|Λg(z)ϕn(y)dy|2,
    as |Δm|0 and |Δ˜m|0. We obtain from Lemma 6.1 that
    Λg(z)Wk(y)dy=l.i.m.|Δm|0Q(g,Δm,Λ).
    Now we have
    Λg(z)Wk(y)dyNn=NΦknΛg(z)ϕn(y)dyΛg(z)Wk(y)dyQ(g,Δm,Λ)+Q(g,Δm,Λ)Nn=NΦknΛg(z)ϕn(y)dy.
    Hence, Lemma 6.3 tells us that the theorem holds. □

    Definition 6.2. We define

    S(g,Δm,Λ)=kmi=1g(z(tmi))(Wk(y(tmi))Wk(y(tmi1))).
    If S(g,Δm,Λ) converges in quadratic mean as |Δm|0, and if the limit does not depend on the choice of the sequence {Δm}, then the limit is denoted by
    Λg(z)dWky.

    We prepare two lemmas.

    Lemma 6.5. Suppose the length of  Λ  is finite. Let  Hm(t)  be a complex-valued random variable sequence. If 

    limmsup0t1E|Hm(t)|2=0,
      then 
    limmkmi=1tmitmi1Hm(t)dxdtdt=0,limmkmi=1tmitmi1Hm(t)dydtdt=0.

    Proof. There is δ>0 such that

    sup0t1E|Hm(t)|2<δ,
    for sufficiently large m. If m is sufficiently large, then we have
    kmi=1tmitmi1Hm(t)dxdtdt2kmi=1kmj=1tmitmi1tmjtmj1E|Hm(t)|2E|Hm(s)|2×|dx(s)ds||dx(t)dt|dsdtδ(kmi=1tmitmi1|dx(t)dt|dt)2.
    Since the length of Λ is finite, the last expression is bounded. As δ0, we get the first limit of the lemma. The second limit can be obtained as well. □

    Lemma 6.6. Suppose that a random function  Hm(t)  on  [0,1]  is defined as follows:

    Hm(0)=0,Hm(t)=gx(z(t))Wk(y(t))gx(z(tmi1))Wk(y(tmi1)),
    ift(tmi1,tmi]. Ifgis of classC1on˜D, then
    limmsup0t1E|Hm(t)|2=0.

    Remark 6.2. In the case where

    Hm(t)=gy(z(t))Wk(y(t))gy(z(tmi1))Wk(y(tmi1)),
    the assertion of Lemma 6.6 remains true.

    Proof. Let t(tmi1,tmi]. We have

    E|Hm(t)|2=|gx(z(t))|2Rk(y(t),y(t))gx(z(tmi1))¯gx(z(t))Rk(y(tmi1),y(t))¯gx(z(tmi1))gx(z(t))Rk(y(tmi1),y(t))+|gx(z(tmi1))|2Rk(y(tmi1),y(tmi1)).
    Hence, the lemma is true. □

    Theorem 6.2. Suppose the length of  Λ  is finite. If  g  is of class  C1  on  ˜D, then the integral  Λg(z)dWky  exists as a quadratic mean limit, which is represented as 

    Λg(z)dWky=g(z1)Wk(d)g(z0)Wk(c)ΛWk(y)(gxdx+gydy)=l.i.m.NNn=NΦkn{g(z1)ϕn(d)g(z0)ϕn(c)Λϕn(y)(gxdx+gydy)}.

    Remark 6.3. If g is continuous on ˜D, the integral Λg(z)Wk(y)dx exists as a quadratic mean limit, which is represented as

    Λg(z)Wk(y)dx=l.i.m.NNn=NΦknΛg(z)ϕn(y)dx.

    Proof. Let Δg(z(tmi))=g(z(tmi))g(z(tmi1)). We have

    S(g,Δm,Λ)=g(z1)Wk(d)g(z(tm1))Wk(c)kmi=2Δg(z(tmi))Wk(y(tmi1))=g(z1)Wk(d)g(z(tm1))Wk(c)kmi=2(gx(z(tmi1))Δx(tmi)+gy(z(tmi1))Δy(tmi)+ϵ1iΔx(tmi)+ϵ2iΔy(tmi))Wk(y(tmi1)).
    Here, ϵji0 uniformly in i as |Δm|0 for j=1,2. This implies that
    S(g,Δm,Λ)g(z1)Wk(d)+g(z0)Wk(c)+ΛWk(y)(gxdx+gydy)|g(z0)g(z(tm1))|Wk(c)+|gx(z0)Δx(tm1)|Wk(c)+kmi=1tmitmi1(gx(z(t))Wk(y(t))gx(z(tmi1))Wk(y(tmi1)))dxdtdt+|gy(z0)Δy(tm1)|Wk(c)+kmi=1tmitmi1(gy(z(t))Wk(y(t))gy(z(tmi1))Wk(y(tmi1)))dydtdt+kmi=2(ϵ1iΔx(tmi)+ϵ2iΔy(tmi))Wk(y(tmi1)).
    We obtain from Lemmas 6.5 and 6.6 that the above expression goes to 0 as m. The theorem has been proved. □

    Finally, we introduce line integrals with respect to Zky, which is not assumed that Λ is closed.

    Definition 6.3. We define

    S0(g,Δm,Λ)=kmi=1g(z(tmi))(Zk(y(tmi))Zk(y(tmi1))).
    If S0(g,Δm,Λ) converges in quadratic mean as |Δm|0, and if the limit does not depend on the choice of the sequence {Δm}, then the limit is denoted by
    Λg(z)dZky.

    Proposition 6.1. Suppose the length of  Λ  is finite. If  g(z)  is of class  C1  on  ˜D, then the integral  Λg(z)dZky  exists as a quadratic mean limit and 

    Λg(z)dZky=Λg(z)dWky+βΛg(z)Wk(y)dy,
      almost surely.

    Proof. We see from (1) that

    S0(g,Δm,Λ)=kmi=1g(z(tmi))(Zk(y(tmi))Zk(y(tmi1)))=kmi=1g(z(tmi))(Wk(y(tmi))Wk(y(tmi1))+βy(tmi)y(tmi1)Wk(s)ds)=kmi=1g(z(tmi))(Wk(y(tmi))Wk(y(tmi1)))+βkmi=1g(z(tmi1))Wk(y(tmi1))Δy(tmi)+βkmi=1(g(z(tmi))g(z(tmi1)))×Wk(y(tmi1))Δy(tmi)+βkmi=1g(z(tmi))×y(tmi)y(tmi1)(Wk(s)Wk(y(tmi1)))ds.
    Since a lemma similar to Lemma 5.1 holds, we can get the proposition. □

    We do not assume that Λ is closed. Let g be a function of class C1 on ˜D and let g1 and g2 be the real part and the imaginary part of g, respectively. By virtue of Proposition 6.1, we have

    p-lim|Δm|0S0(gj,Δm,Λ)=Λgj(z)dZky,
    for j=1,2. This implies that
    E[eiξΛgj(z)dZky]=exp[ΛΨk(gj(z)ξ)dy].
    As a result, we obtain that
    limnE[eiξnk=0(γkΛgj(z)dy+Λgj(z)dZky)]=exp[ΛΨ(gj(z)ξ)dy].
    We here note that γkΛgj(z)dy+Λgj(z)dZky,k=0,1,2,, are independent. Since nk=0(γkΛgj(z)dy+Λgj(z)dZky) converges in probability, we can define line integrals along Λ with respect to Z as follows:
    Λg(z)dZ=Λg(z)dx+ik=0(γkΛg(z)dy+Λg(z)dZky).

    Theorem 6.3. Suppose the length of  Λ  is finite. If  g  is of class  C1  on  ˜D, then 

    Λg(z)dZ=Λg(z)dx+ik=0(γkΛg(z)dy+Λg(z)dWky+βΛg(z)Wk(y)dy)=Λg(z)dx+ik=0(γkΛg(z)dy+l.i.m.NNn=NΦkn×{g(z1)ϕn(d)g(z0)ϕn(c)Λϕn(y)(gxdx+gydy)+βΛg(z)ϕn(y)dy}a.s..

    Proof. Proposition 6.1 tells us that almost surely

    Λg(z)dZ=Λg(z)dx+ik=0(γkΛg(z)dy+Λg(z)dWky+βΛg(z)Wk(y)dy).
    Hence, it follows from Theorems 6.1 and 6.2 that the theorem holds. □

    Acknowledgments

    The author appreciates the referee’s comments for improving the paper.

    Competing Interests

    The author declares that he has no competing interests.

    Contribution

    This work has been done alone. The author read and approved the final paper.