MSC 35R30, 35Q35, 65M60, 65M32 DOI: 10.14529/mmp190107
ON SOME INVERSE COEFFICIENT PROBLEMS WITH THE POINTWISE OVERDETERMINATION FOR MATHEMATICAL MODELS OF FILTRATION
S.N. Shergin1, E.I. Safonov1, S.G. Pyatkov12
:Ugra State University, Khanty-Mansyisk, Russian Federation 2South Ural State University, Chelyabinsk, Russian Federation E-mails: S_pyatkov@ugrasu.ru, Ssn@ugrasu.ru, dc.gerz.hd@gmail.com
We examine inverse problems of recovering coefficients in a linear pseudoparabolic equation arising in the filtration theory. Boundary conditions of the Neumann type are supplemented with the overtermination conditions which are the values of the solution at some interior points of a domain. We expose existence and uniqueness theorems in the Sobolev spaces. The solution is regular, i. e., it possesses all generalized derivatives occurring in the equation containing in some Lebesgue space. The method of the proof is constructive. The problem is reduced to a nonlinear operator equation with a contraction operator whenever the time interval is sufficiently small. Involving the method of the proof, we construct a numerical algorithm, the corresponding software bundle, and describe the results of numerical experiments in the two-dimensional case in the space variables. The unknowns are a solution to the equation and the piezo-conductivity coefficient of a fissured rock. The main method of numerical solving the problem is the finite element method together with a difference scheme for solving of the corresponding system of ordinary differential equations. Finally, the problem is reduced to a system of nonlinear algebraic equations which solution is found by the iteration procedure. The results show a good convergence of the algorithms.
Keywords: inverse problem; pseudoparabolic equation; filtration; fissured rock; numerical solution.
Introduction
In 1960, G.I. Barenblatt, Iu.P. Zheltov and I.N. Kochina (see [1]) proposed the basic concept in the theory of seepage (filtration) of homogeneous liquids in fissured rocks. A fissured rock is considered as a material consisting of pores and permeable blocks which are generally separated from each other by a system of fissures. In contrast to to the conventional arguments of filtration in a porous medium, the significant feature lies in the fact that two liquid pressures, both in the pores and in the fissures, are introduced at any point in a space and the transfer of liquids between the fissures and the pores is taken into consideration. The corresponding model is written as
vA«1 + a(u2 — u1) = 0, m0d1 + d2)u2t + a(u2 — u1) = 0, (1)
where u1(t, x), u2(t, x) are the pressures of the liquid in the fissures and pores, respectively; A is the Laplacian, d1 and d2 are the coefficients of compressibility of the liquid and the blocks, m0 stands for the magnitude of the porosity of the blocks at standard pressure; ^ is the viscosity of the liquid, and v represents the permeability of fissures. The dimensionless coefficient a characterizes the intensity of the liquid transfer between the blocks and fissures. More general models can include the nonlinearities arising from fluid type (liquid or gas), concentration (porosity, absorption or saturation) and the exchange rate [2].
Eliminating u2 from (1) we obtain for the pressure of the liquid in the fissures the so-called fissured medium equation of pseudoparabolic type
u1t — nAu1t — kAu1 = 0, k = v/(m0d1 + d2), n = v/a. (2)
The parameter k corresponds to the piezo-conductivity of fissured rock. The pressure of the liquid in the pores u2 satisfies a similar equation. Since the natural stratum is involved, the parameters of fissured rock in (2) should be determined on the basis of the investigation of their behavior under the natural nonsteady-state conditions. This leads to the interest in studying the inverse problems for equation (2) and its analogues. More general equations can be written in form
L(t, x, D)ut — M(t, x, D)u = f, (x,t) G Q = G x (0, T), (3)
where L, M are second order operators and G is a bounded domain in Rn. The equation (3) is furnished with initial and boundary conditions of form
u(0,x) = uo(x), Ru|S = g (t,x) (4)
n
with Ru = u or Ru = Y%(t, x)uXi+&(t, x)u (other boundary conditions are also possible).
i=i i
Pseudoparabolic equations of form (3) with various differential operators L1 and L2 of the even order in spacial variables also arise in the mathematical models of the heat conduction, wave processes quasistationary processes in semiconductors and magnetics, in the models for filtration of the two-phase flow in porous media with the dynamic capillary pressure (see [3,4, Sect. 0.1.4] and the bibliography therein). Detailed bibliography and the results concerning the solvability of direct problems for pseudoparabolic equations can be found, for instance, in [5,6]. The first results devoted to inverse problems for pseudoparabolic equations were obtained in [7], where an inverse problems of recovering an unknown source f of a special form in (3) is considered. Large number of results is exposed in the monographs [8-10]. We mention also the article [11] devoted to some coefficients inverse problems with lower order coefficients depending on spatial variables being the unknowns. The problems of recovering coefficients, in particular, the coefficients k(t) and n are studied in [12,13], where integral overdetermination conditions are used. Closed results with pointwise overdetermination conditions are presented in [14]. Exposition of numerical methods for solving inverse problem can be found, for instance, in [15,16]. We can refer also to articles [17-24] devoted to different numerical methods of solving boundary value problems for pseudo-parabolic equations. At the same time, the number of articles devoted to a numerical solution of inverse problems for pseudo-parabolic equations is rather limited (see, for instance, [24,25]). Most of the articles are devoted to different model problems. In the present article we generalize the results in [14] to the case of the oblique derivative problem, describe numerical methods applicable to a wide class of inverse problems with pointwise ovedetermination, and present the results of numerical experiments in the case of recovering the piezo-conductivity coefficient of fissured rock in equation (2).
1. Preliminaries
We consider a general inverse problem on recovering functions occurring into the right-hand side and left-hand side of the equation. Let the right-hand in (3) be of form
ro
f = ^ Ci(t)fi(x,t) + f°(x,t), fi E L^(0,T; LP(G)), (5)
i=1
where functions fi are given. We assume that
n n
L aij(t, x)3XiXj + ai(t, x)3Xi + ao(t, x)
i,j=1 i=1
and operator M is representable as
r
Mu = Mr0u + Y^ ck(t)Mku,
k=ro+1 nn
Mku = Y, bkj (x, t)uxiXj + E bk (x, t)uxi + b° (x, t)u, i,j=1 i=1
where functions ci(t) are unknown. Our problem is stated as follows: find functions [ci(t)}Tr=1 and solution u to the problem (3) - (4) such that
u(xi,t) = ai(t), (i = 1,2,..,r), (6)
where Xi are arbitrary points lying in G and «¿(t) are given functions. Put (u,v) = jG u(x)v(x) dx.
We employ the Sobolev spaces W^{G) and Holder spaces Ca(G) (see the definitions in [26]). Symbol Lp(0, T; H) (H is a Banach space) stands for the space of strongly measurable functions defined on [0,T] with values in H. Given an interval J = (0,T) and domain G E Rn, put Q = (0,T) x G and W;,s(Q) = W;(J; Lp(G)) n Lp(J; Ws(G)). Respectively, W*'r(S) = W;(J]LP(T)) n LP(J-W;(T)) (S = (0,T) x DC). Similarly, we can define the Holder spaces Cr's(Q).
Next, we describe the condition on the data of the problem. We assume that operator L is elliptic, i. e., there exists constant 8 > 0 such that
m
V£GiT, W(x,t)eQ. (7)
i,j=1
Fix parameter p > n and assume that
bkj E L^ (Q),bk,bk E L^ (0,T; LP(G)) (k = ro,...,r), (8)
al3eC(Q), at,aoeC{[0,T];Lp{G)) {i,j = l,2,..,n). (9)
We also suppose that
a0(x,t) < 0 a.e. (almost everywhere) in Q in the case of the Dirichlet boundary conditions and a° < 0 a.e. in Q and a0 < 0 a.e. in some neighborhood about S in the case of the oblique derivative problem.
(A)
We also assume that
G C1/2'1 (S), i = l,2(10)
ai(0) = uo(xi,t) (i = 1, 2,..,r), B(0,x,D)uo\r = g(0,x). (11)
Let s0 = 2 — 1/p in the case of the Dirichlet boundary condtions and so = 1 — 1/p otherwise. We can determine function $ G C([0,T]; Wp(G)) (p > n) such that $t G Lp(0,T; W£(G)), $\t=o = uo(x), and R$\S = g. Construct matrix B with the rows
L-Ifl(xj ,t),L-1fro (xj ,t), —L-lMro+1 $(xj ,t),.., —L-1Mr $(xj ,t),
where j = 1, 2, ...,r and assume that the there exists constant So > 0 such that
\ det B\ > So Vt G [0,T]. (12)
Here L-1 fi is solution Ui to problem LUi = fi, Ui\t=0 = 0, RU0\S = 0. In this case, it is easy to justify that locally in time conditions (12) does not depend on the choice of function $.
The following theorem is proven in [14] for the case of the Dirichlet boundary conditions. Here we expose the claim in the general case.
Theorem 1. Let the conditions (7) - (12), (A) be fulfilled. Then there exists a constant Y0 > 0 such that on interval [0,y0] problem (3) - (6) has unique solution (u,c1,... ,cr) such that
u G W1(0, Yo; W2(G)), Ci (t) G Lp(0, Yo)(i = 1, 2,...,r).
Proof. In the case of the oblique derivative problem the proof almost repeats that in the case of the Dirichlet boundary conditions. So we present the proof omitting some details. Let function $ G C([0,T]; W2(G)) such that $t G Lp([0,T]; W2(G)) be a solution to the direct problem (see the existence theorem in [14])
L$t + Mr0$ = fo, $\t=o = uo, R$\s = g. In this case function v = u — $ is a solution to equation
ro r
Lvt + Mv = £ Ci(t)fi(x,t) — E Ci(t)Mi$ = F(c), (13)
i=1 i=ro+1
satisfying homogeneous conditions (4). There exists yo > 0 such that condition (12), with So/2 rather than So, holds on [0,Yo]. In this case
v(t,x3) = aj (t) — $(t,xj) = otj (t). (14)
Function $ G Wp(G; W^(0,T)) after a possible modification on a set of zero measure possesses property $ G Ca(G; W£(0,T)) for a < 2 - ^ (see (5.4) in [28]). In particular, ,t) G W^(0,T), and thus
aj (t) G W1(0,T) (j = 1, 2,..,r). (15)
Inverting L in (13), we have
vt + L-1 Mv — vo (v) = L-1 F (c), (16)
where vo(v) = 0 in the case of the Dirichlet boundary conditions and is a solution to problem Lvo = 0, Rvo\S = —Rtv otherwise.
Note that L-1 Mi$(t,xj) and L-1Mi$(t,xj) G L^(0,T) for almost all t. Actually, L-1 Mi$(x,t) G W2(G; Lp(0,T)) and thereby L-1Mi$(t,xj) G Lp(0,T) (see [28, (5.4)]). Let x = xj in (16). We have
vt(t,xj) + L-1Mv(t,xj) — vo (v)(t,xj) = L-1F (c)(t,xj).
The overdetermination conditions yield
oijt + R(v)(t,xj) = L-1F (c)(t,xj),
r
R(v) = L-1Mr0 v — vo(v) + Y^ ci (t)L-1 Mi v.
i=ro+1
We can rewrite the last equations in form
( a1t + R(v)(t,x1) \
Bc
(17)
As a result, we infer
V art + R(v)(t,xr)) ( au + R(v)(t,x1) \
B
-1
S(c). (18)
\ art + R(v)(t,xr)) /
The right-hand side can be viewed as an operator S(c) taking vector c(t) in solution v(x, t) to equation (13) satisfying homogeneous conditions (4) and then taking v(x,t) into the right-hand side of (18). System (18) is the desired system. The properties of this operator S can be described as follows. Fix
a
Ro = 2\\B la||lp(o,t), a = (au, ...,aH)T,
where norm \\B-1 a\\lp(o,t) is just the sum of the norms of the coordinates of vector B-1 in space Lp(0,T). We assume that
c(t) G broy = {c G Lp(0,y) : \\c1\lp(o;7) < Ro}, Y < T.
We now can state that there exists constant c > 0 and Yo > 0 such that operator S : br0y ^ BRoY, y < Yo takes ball BRoY into itself and is a contraction. More exactly, we have the estimate
\\S(c1) — S(c2)\\lp(Oyo) < cy 1/q\\c1 — C21lp(o7o) Vc1 ,c2 G bro
and thus this operator is a contraction for all y < Yo with cyo/q = 1/2. The method of successive approximations cn = S(cn-1) (n = 1, 2,...), c° = 0 converges to a solution to equation (18). The proof of these facts is in line with that in [14]. So we omit it.
□
2. Description of the Algorithm
To simplify the presentation, we describe the idea of the algorithm in the model case which in particular includes the filtration problems described in the introduction. We rely on some integral identities. Consider problem
Lo ut + k(t)Li u = f, (19)
du
— |г = g, и(0) = щ(х), (20)
dn
u(xo,t) = фо (t), (21)
where
L0u = -div(a0(x, t)Vut) + b0(x, t) ■ Vu + c0(x, t)u,
L1 u = -div(al(x, t)Vu) + bl(x, t) ■ Vu + cl(x, t)u
and ao,al,co,cl are scalar functions and b0,bl are vector-function of length n. Functions u and k(t) are unknown. We assume that all conditions of Theorem 1 for the data are fulfilled. Function Ф is a solution to problem (19), (20). Since the existence theorem is local in time, we can replace condition (12) with the following condition: there exists t0 > 0 such that
lL-lLl^(xo,t)l > So > 0 yt e [0,to]
or in a simpler form
lL-lLluo(xo)l > So > 0 yt e [0,to]- (22)
Let ç e Lq(0,T; W^(G)) (1/q + 1/p =1) be a test function and let a function u be a solution to problem (19), (20) from the class pointed out in Theorem 1. Integrating by parts in identity
(Lout, ç) + k(t)(Llu, ç) = (f,ç), ç e Lq (0, T ; Wql(G)), (23)
we arrive at equality
a(ut, ç) + k(t)b(u, ç) = l(ç) + k(t)ll(ç), (24)
where a(ut, ç) = (aoVut, Vç) + (bo■Vut, + cout, ç), b(u, ç) = (alVu, Vç) + (bl■Vu + clu, ç),
l(ç) = (f, ç) + lo(ç), lo(ç) = / aogtçdT, ll(ç) = / algçdt.
г г
Next, we look for solution ço(x,t) to problem
L*ip0 = S(x-x0), a0ff+ 6-^o|r = 0, (25)
where Lo is a formally adjoint to Lo and S is the Dirac delta-function. Inserting ço in (24), we obtain that
фг + k(t)(b(u, ço) - ll(ç)) = l(ço). (26)
Hence, we conclude that
k(t)(b(u, ço) - ll(ço)) = l(ço) - фь- (27)
Note that definition of p0 implies that b(u,p0) — l1(p0) = L-1 L1u(x0,t). Since IL-1 L1 u0(x0)| > 60 > 0 on some segment [0,t0] we can say that I(b(u0, p0) —l^p0))| > 0 on some segment of time. Integral identities (27) allows us to construct the iteration procedure realized in the proof of Theorem 1. To avoid an excessive calculations, we can omit the determination of $ in our case. Let k0 = (l0(p0) — ^t)/(b(u0,p0) — l1(^0))|t=0. Given function ki, we can construct ui+1 as a solution to problem (19), (20) with k(t) = ki(t) and to determine next iteration ki+1 from one of the equalities
ki+1 (t) = (l(<po) — 4H)/(b(ui+1, <po) — h(<po)), (28)
ki+1 = (l(po) — ^t — kib(ui+1 — uo, <po))/(b(uo, po) — h(po)). (29)
The latter formula almost corresponds to the iteration procedure in the proof of the fixed point theorem for operator S constructed in Theorem 1. The peculiarity is that we have replaced function $ with function u0. It can be easily clarified that all arguments remain valid in this case. Denominator in (29) is different from zero on some segment [0, t0]. Former formula (28) corresponds an another operator S but the same arguments can be used to prove that the iteration procedure corresponding to formula (28) converges as well on some small time segment [0, t0].
3. Numerical Algorithm
The algorithm is iterative and relies on the finite element method. We define triangulation of G, mesh nodes, X1, X2, Xn , and corresponding piecewise linear functions [pi(x)} (thus, pi(xj) = 6ij, where 5ij is the Kronecker symbol. Without loss of generality, we can assume that observation point x0 is mesh node xj0. Approximate solution to (19),
N
(20) is sought in form um+1 = ^ ci(t)pi(x). To determine the functions ci we employ
i=1
integral identity (24). Vector-function C(t) = (c1 (t), c2(t),..., cN(t))T is a solution to system of ordinary differential equations
AC\ + k(t)BC = F, C(0) = (uo(x1 ),uo(x2),... ,uo(xn))T, (30)
where A,B are matrices with entries aij = a(pj,Pi), bij = b(pj,pi), and F = (F1,..., Fn)t, Fi(t) = l0(pi) + k(t)l1(pi). To solve (30), we involve the finite difference method (FDM) (the implicit scheme) and replace (30) with finite difference equation
AnCn~Cn~l + knBnCn = F0n + knF1>n, Co = C{0), (31)
T
where n = 1, 2,..., M, t = T/M, and F0n, F1n, An, Bn are the values of the right-hand side in (30), and matrices A,B at nT. We assume here that approximation k of k is a piecewise constant function taking value kn on ((n — 1)t,vt]. Respectively, a piecewise constant approximation of solution C(t) to (30) is a piecewise constant function equal to vector Cn on set ((n — 1)t, nr]. An analog of the overdetermination condition is as follows:
Cn~Cn~l =ijot((n-l)r). (32)
T
Since our results are local in time, we employ the predictor-corrector arguments. Given vector C, denote its j-th coordinate by (C)j. Find quantity k0 = ((A-1F0,0)j0 —
^t(0))/((Ao 1B0C0)j0 — (A0 1 F1,0)j0). Put k\ = k0. Let the approximation k\ of quantity k1 be known. Find solution C[ to equation (31) with n = 1 and k1 = k\, i.e.
A1Cl~Co + k\BnC\ = Fq i + k\Fltl. (33)
T
The next approximation is defined as
k\+1 = ((A-1Fo,i)j0 — fr (O))/((A01BiC\ )j0 — (A-1 F^j). (34)
We repeat the arguments until Ik^ — k\l < e with e a prescribed small number. Next we take k1 := k]11. We have find quantity k1 and corresponding vector C1 := C1+1. Assume that we have found k1,... ,km (m < M) and vectors C1,..., Cm. Take k]m+1 := km and calculate quantities klm+1 for all i as follows. Assume that quantity klm+1 is known. Find solution Clm+1 to equation (31) with n = m +1 and kn = klm+1) i.e.
C — Cm
(35)
Next, we calculate
K+h = ((Am\1Fo,m+1)jo — MmT))/rm+1, (36)
where llm+1 = ((Am1+1Bm+1 Clm+1)jo — (A~m+1F\,m+1 )jo). We repeat the arguments until Wm^ — ktm+1I < e. In this case we put km+1 := klr+1^1 and Cm+1 := Clr+1^1. The arguments are repeated until m = M. It is possible of course that the process stops at some iteration (for instance, if the denominator in (36) becomes equal to zero). It is possible due to the nonlinearity of the problem. Formula (36) corresponds to equality (28).
4. The Results of Numerical Experiments
In this section we analyze the results of numerical experiments. The characteristics of the computer are as follows: processor Intel(R) Core(TM) i7-3517U CPU @ 1.90GHz 2.40GHz, 10.00 GB RAM, 64-digit operating system Windows 7 Enterprise.
As a result of calculations, we obtain approximate values of solution (u(x,y,t), k(t)) of problem (19) - (20) at points (t1 ,t2,... tN). Here point (x,y) belongs to the unit circle centered at (0, 0). We present the results of calculations only for function k.
To solve the problem numerically, we use three meshes for this domain with the number of nodes N1 = 197, N2 = 751 and N3 = 2933 (Fig. 1). Additional information (21) is given at observation point Xj0 = (x0, y0) = (0, 3, —0, 3).
All numerical experiments are divided into two groups in dependence on unknown functions u, k, the boundary conditions, noise percentage 5, error between iterations e, coefficients a0,a1 ,b0>1,b0>2,bi>1,b1>2,c0,c1, and right-hand sides f.
We will use the Neumann boundary conditions from (24) which are represented as
du dn
|r = (uxx + Uy y) = g.
The data of the first group are as follows:
- solution: u(x, y, t) = (x2 + y2 + 1) ■ (1 + t);
- initial data: uIt=0 = x2 + y2 + 1;
Graph with points xl yl ( Graph with points x1 y 1, Graph with pointe x1 y1
-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8
XXX
a) b) c)
Fig. 1. Meshes: a) N1 = 197; b) N2 = 751; c) N3 = 2933
- Neumann boundary conditions: g = 2(x2 + y2)(1 + t);
- additional information: ^(xj0, t) = (xl + y02 + 1) • (1 + t);
- unknown function: k(t) = t2;
-coefficients: a0 = 0, 5, a1 = 4/(1+1), b0^ = b0^2 = b1^1 = b1^2 = 0, c0 = 1/(x1+x<2 + 1), c1 = 1/(1 + t);
- right-hand side: f = t2 • (x2 + y2 — 15) — 1.
We compare the results of calculations for three meshes for the data without noise, i.e. S = 0. Take e = 10-3 (the error defined by the user). Denote by r the time of the work of the algorithm in seconds. Introduce one more error of the algorithm by equality e0 = max Ikm — k (mAt)I, where m = 1, 2,... ,N. Consider segment [0,T], T = 1, and take
m
At = T/N, N = 100.
1.2 i-,-1-1-T-
0.8 / -
0.6 //
0.4
0.2
ol -—-1-1-1-
0 0.2 0.4 0.6 0.8 1
t
Fig. 2. Results of calculations for the meshes: : т = 9, 81, ^o = 0, 0058; N2: т = 131, 62, so = 0, 0015; N3: т = 3820, 1, ^o = 0, 00035
The obtained plots of original functions k(t) and their approximation for three different meshes are almost identical, so we present only one of them for the first mesh. For the following experiments, add 10, 20 and 40 percent random noise to the overdetermination data: ijjt (x0, t) = фt (x0, t)(1 + S(2a — 1)), where a is random function normally distributed on interval [0,1]. Due to the almost complete coincidence of the graphs of the original and the approximated function, we give only the graphics with noise 40 percent and s = 10-3.
a) b) c)
Fig. 3. Results of calculations for the meshes: a) N1: t = 8, 41, e0 = 0, 0286; b) N2: t = 110, 5, £q = 0, 0231; c) N3: t = 3983,1, e0 = 0, 0232
Fig. 3 shows that, the obtained plots of original functions k(t) and their approximation are almost identical and regardless of the level of input noise 5, the results of the calculations repeat the ones sought or are located next to them. Table 1 shows the results for different input data.
Table 1
Results of numerical experiments (first group)
No ехр. Grid 5 T
1 N1 0 0,0058 8,41
2 N1 0,1 0,0106 8,8
3 N1 0,2 0,0166 8,51
4 N, 0,3 0,0187 8,3
5 n2 0 0,0015 110,5
6 n2 0,1 0,0066 109,4
7 n2 0,2 0,0114 114,32
8 n2 0,3 0,018 118,57
9 N3 0 0,00035 3983,1
10 N3 0,1 0,0058 3773,3
11 N3 0,2 0,0117 3730,4
12 N3 0,3 0,0174 3821,9
Based on the results of numerical experiments for the first group of data, we can conclude that the increase in the number of nodes in 4 times leads to the increase in the calculation time by an average of « 12 times for grids N1 and N2 and « 42 times for grids N2 and N3, but does not lead to a significant increase in accuracy and with the increase in the noise level.
Let us describe the data of the second group:
- solution: u(x, y, t) = (x2 + 1) • (y2 + 1) • (1 + t);
- initial data: uIt=0 = (x2 + 1) • (y2 + 1);
- Neumann boundary conditions: g = 2(t + 1)(y(x2 + 1) + x(y2 + 1));
- additional information: ^ (xj0, t) = (xI + 1) • (y02 + 1) • (1 + t);
- unknown function: k(t) = (t — 0, 5)3;
-coefficients: a0 = (t+1)(x2 + 1), a1 = (t2+y+4)/(x2+1), b= x3*t, b02 = (x+y)*t, bM = y2/(t + 1), b12 = xt/(y + 1), Co = ty/(x2 + 1), C1 = (x + t)/(2 + y);
- right-hand side: f = (2(x2 + 1)2)/(t + 1) + (t — 0.5)3(((x2 + 1)(y2 + 1))/(t + 1)2 — (2y + 2t2 + 8)/(t + 1) — (2y )/(t +1) — (2(y2 + 1)(y +12 + 4))/((x2 + 1)(t + 1)) + (2xy2(y2 + 1))/(t + 1)2 + (4x2(y2 + 1)(y + t2 + 4))/((x2 + 1)2(t + 1)) + (2txy (x2 + 1))/((t + 1)(y + 1))) + (2(x2 + 1)(y2 + 1))/(t + 1) + (4x2(y2 + 1))/(t +1) — (2tx4 (y2 + 1))/(t + 1)2 — ((x2 + 1)(y2 + 1))/((t + 1)2 (x2 + y2 + 1)) — (2ty (x2 + 1)(x + y ))/(t + 1)2.
Since the results for different grids are similar, we present them for the grid N1. As in the case of the first group, we vary parameters e and S. Table 2
Results of numerical experiments for the second group of data
No exp. Grid 8 £ r
1 N, 0 10"4 0,0035 7,1
2 N! 0 10"5 0,0035 6,35
3 N! 0 10"6 0,0035 6,59
4 N! 0 10"7 0,0035 7,02
5 N! 0,2 10"4 0,0156 6,2
6 N! 0,2 10"5 0,0126 5,85
7 N! 0,2 10"6 0,0151 6
8 N! 0,2 10"7 0,0137 6,6
9 N! 0,4 10"4 0,0291 6,38
10 N! 0,4 10"5 0,0246 6,94
11 N! 0,4 10"6 0,0262 8,08
12 N! 0,4 10"7 0,0273 8,2
As we can see the decrease in variable e does not lead to a significant increase in the
accuracy and decreasing the time of calculations. In the next two experiments we reduce
the time step in 5 times. , , „
Table 3
Results of numerical experiments for the second group of data,
time step 0,002
No exp. Grid 5 £ So r
1 Ni 0 10"4 0,0031 36,41
2 N! 0 10"5 0,0031 32,69
3 N! 0 10"6 0,0031 32,41
4 N! 0 10"7 0,0031 35,56
5 N! 0,2 10"4 0,0147 47,44
6 N! 0,2 10"5 0,0146 37,77
7 N! 0,2 10"6 0,0142 36,28
8 N! 0,2 10"7 0,0153 42,46
9 N! 0,4 10"4 0,0263 45,27
10 N! 0,4 10"5 0,0269 56,34
11 N! 0,4 10"6 0,0273 56,35
12 N! 0,4 10"7 0,0272 54,81
Summing up, we can say that the use of a grid with a large number of nodes shows better accuracy, but the time of the calculation increases by an amount equal to the ratio
of the number of nodes of the three grids. Decreasing the variable e leads to an increase in the computation time T, but does not lead to a significant increase in accuracy. We also see that the dependence of the time of a calculation on the time step is inversely proportional.
Conclusions
We establish the existence and uniqueness theorems in inverse problems of recovering the coefficients of a pseudoparabolic equation with pointwise overdetermination conditions. Actually, we have proven that the problems in questions are well-posed in the usual sense at least locally in time. The problem is reduced to an operator equation of Volterra type. This allows to construct a numerical algorithms based on the conventional methods (in our case FEM and difference schemes) with a sufficiently good convergence to a solution. The results of numerical experiments are presented.
Acknowledgements. This work was supported by the science foundation of Yugra State University (Grant no. 13-01-20/16) and by the Act 211 of the Government of the Russian Federation, contract no. 02.A03.21.0011.
References
1. Barenblatt G.I., Zheltov Iu.P., Kochina I.N. Basic Concepts in the Theory of Seepage of Homogeneous Liquids in Fissured Rocks. Journal of Applied Mathematics and Mechanics, 1960, vol. 24, no. 5, pp. 1286-1303. DOI: 10.1016/0021-8928(60)90107-6
2. Bohm M., Showalter R.E. Diffusion in Fissured Media. SIAM Journal on Mathematical Analysis, 1985, vol. 16, no. 3, pp. 500-509. DOI: 10.1137/0516036
3. Lyubanova A.Sh., Tani A. On Inverse Problems for Pseudoparabolic and Parabolic Equations of Filtration. Inverse Problems in Science and Engineering, 2011, vol. 19, no. 7, pp. 1023-1042. DOI: 10.1080/17415977.2011.569712
4. Al'shin A.B., Korpusov M.O., Sveshnikov A.G. Blow-up in Nonlinear Sobolev Type Equations. Berlin, N.Y., De Gruyter, 2011. DOI: 10.1515/9783110255294
5. Sviridyuk G.A., Fedorov V.E. Linear Sobolev Type Equations and Degenerate Semigroups of Operators. Utrecht, VSP, 2003. DOI: 10.1515/9783110915501
6. Favini A., Yagi A. Degenerate Differential Equations in Banach Spaces. N.Y., Marcel Dekker, 1999.
7. Rundell W. Determination of an Unknown Nonhomogeneous Term in a Linear Partial Differential Equation from Overspecified Boundary Data. Journal of Applied Analysis, 1980, vol. 10, no. 3, pp. 231-242. DOI: 10.1080/00036818008839304
8. Kozhanov A.I. Composite Type Equations and Inverse Problems. Utrecht, VSP, 1999. DOI: 10.1515/9783110943276
9. Asanov A., Atamanov E.R. Nonclassical and Inverse Problems for Pseudoparabolic Equations. Berlin, De Gruyter, 2014.
10. Favini A., Lorenzi A. Differential Equations. Inverse an Direct Problems. Abingdon, Tylor and Francis Group, 2006. DOI: 10.1201/9781420011135
11. Mamayusupov M.Sh. The Problem of Determining Coefficients of a Pseudoparabolic Equation. Studies in Integro-Differential Equations, 1983, vol. 16, pp. 290-297.
12. Lyubanova A.Sh. Identification of a Coefficient in the Leading Term of a Pseudoparabolic Equation of Filtration. Siberian Mathematical Journal, 2013, vol. 54, no. 6, pp. 1046-1058. DOI: 10.1134/S0037446613060116_
13. Lyubanova A.Sh. The Inverse Problem for the Nonlinear Pseudoparabolic Equation of Filtration Type. Journal of Siberian Federal University. Series: Mathematics and Physics, 2017, vol. 10, no. 1, pp. 4-15.
14. Pyatkov S.G., Shergin S.N. On Some Mathematical Models of Filtration Theory. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2015, vol. 8, no. 2, pp. 105-116.
15. Kabanikhin S.I. Inverse and Ill-Posed Problems. Berlin, Boston, De Gruyter, 2012.
16. Samarskii A.A., Vabishchevich P.N. Numerical Methods for Solving Inverse Problems of Mathematical Physics. Berlin, Boston, De Gruyter, 2007. DOI: 10.1515/9783110205794
17. Ewing R.E. Numerical Solution of Sobolev Partial Differential Equations. SIAM Journal on Numerical Analysis, 1975, vol. 12, no. 3, pp. 345-363. DOI: 10.1137/0712028
18. Cuesta C.M., Pop I.S. Numerical Schemes for a Pseudo-Parabolic Burgers Equation, Discontinuous Data and Long-Time Behaviour. Journal of Computational and Applied Mathematics, 2009, vol. 224, no. 1, pp. 269-283. DOI: 10.1016/j.cam.2008.05.001
19. Beshtokov M.Kh. Differential and Difference Boundary Value Problem for Loaded Third-Order Pseudo-Parabolic Differential Equations and Difference Methods for Their Numerical Solution. Computational Mathematics and Mathematical Physics, 2017, vol. 57, no. 12, pp. 1973-1993. DOI: 10.1134/S0965542517120089
20. Vabishchevich P.N., Grigor'ev A.V. Splitting Schemes for Pseudoparabolic Equations. Differential Equations, 2013, vol. 49, no. 7, pp. 807-814. DOI: 10.1134/S0012266113070033
21. Guezane-Lakoud A., Belakroum D. Time-Discretization Schema for an Integrodifferential Sobolev Type Equation with Integral Conditions. Applied Mathematics and Computation, 2012, vol. 218, no. 9, pp. 4695-4702. DOI: 10.1016/j.amc.2011.11.077
22. Luoa Z.D. Teng F. A Reduced-Order Extrapolated Finite Difference Iterative Scheme Based on POD Method for 2D Sobolev Equation. Applied Mathematics and Computation, 2018, vol. 329, pp. 374-383. DOI: 10.1016/j.amc.2018.02.022
23. Xia H., Luo Z. An Optimized Finite Difference Crank-Nicolson Iterative Scheme for the 2D Sobolev Equation. Advances in Difference Equations, 2017, vol. 196, 12 p. DOI: 10.1186/s13662-017-1253-8.
24. Keller A.V., Zagrebina S.A. Some Generalizations of the Schowalter-Sidorov Problem for Sobolev-Type Models. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2015, vol. 8, no. 2, pp. 5-23.
25. Hasan A., Aamo A.M., Foss B. Boundary Control for a Class of Pseudo-Parabolic Differential Equations. Systems and Control Letters, 2013, vol. 62, no. 1, pp. 63-69. DOI: 10.1016/j.sysconle.2012.10.009
26. Triebel H. Interpolation Theory, Function Spaces, Differential Operators. North-Holland Mathematical Library. Amsterdam, North-Holland Publishing, 1978.
27. Ladyzhenskaya O.A., Solonnikov V.A., Ural'tseva N.N. Linear and Quasi-Linear Equations of Parabolic Type. Providence, American Mathematical Society, 1968. DOI: 10.1090/mmono/023
28. Amann H. Operator-Valued Fourier Multipliers, Vector-Valued Besov Spaces and Applications. Mathematische Nachrichten, 1997, vol. 186, no. 1, pp. 5-56. DOI: 10.1002/mana.3211860102
Received December 7, 2018
УДК 517.956
БСТ: 10.14529/mmp190107
0 НЕКОТОРЫХ КОЭФФИЦИЕНТНЫХ ОБРАТНЫХ ЗАДАЧАХ
С ТОЧЕЧНЫМ ПЕРЕОПРЕДЕЛЕНИЕМ ДЛЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ФИЛЬТРАЦИИ
С.Н. Шергин1, Е.И. Сафонов1, С.Г. Пятков1,2,
1 Югорский государственный университет, г. Ханты-Мансийск, Российская Федерация
2Южно-Уральский государственный университет, г. Челябинск, Российская Федерация
Рассматриваются обратные задачи восстановления коэффициентов линейного псевдопараболического уравнения, возникающие в теории фильтрации. Граничные условия типа Неймана дополняются условиями переопределения, которые есть значения решения в некотором наборе внутренних точек области. Мы приводим теоремы существования и единственности решений в пространствах Соболева. Полученное решение является регулярным, то есть обладает всеми обобщенными производными, входящими в уравнение, принадлежащими некоторому пространству Лебега. Метод доказательства является конструктивным. Задача сводится к нелинейному операторному уравнению с сжимающим оператором, если временной промежуток достаточно мал. Используя метод доказательства, мы строим численный алгоритм определения решения, соответствующий программный комплекс и описываем результаты численных экспериментов в двухмерном случае по пространственным переменным. Определению подлежат само решение уравнения и коэффициент пьезопроводимости трещиноватой среды. Основной метод для численного определения решения - метод конечных элементов, который дополняется разностной схемой для решения соответствующей системы обыкновенных дифференциальных уравнений. В конечном счете задача сводится к решению нелинейной алгебраической системы, решение которой находится при помощи итерационной процедуры. Результаты показывают очень хорошую сходимость численного алгоритма.
Ключевые слова: обратная задача; псевдопараболическое уравнение; фильтрация; численное 'решение.
Сергей Николаевич Шергин, аспирант, ведущий программист, отдел по набору и профориентационной работе, Югорский государственный университет (г. Ханты-Мансийск, Российская Федерация), ssn@ugrasu.ru.
Егор Иванович Сафонов, кандидат физико-математических наук, доцент, Высшая цифровая школа, Югорский государственный университет (г. Ханты-Мансийск, Российская Федерация), dc.gerz.hd@gmail.com.
Сергей Григорьевич Пятков, доктор физико-математических наук, профессор, заведующий кафедрой «Высшая математика:», Югорский государственный университет (г. Ханты-Мансийск, Российская Федерация); научно-исследовательская лаборатория «Неклассические уравнений математической физики», Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), S_pyatkov@ugrasu.ru
Поступила в редакцию 7 декабря 2018 г.