Next Article in Journal
AI-Enabled Condition Monitoring Framework for Outdoor Mobile Robots Using 3D LiDAR Sensor
Previous Article in Journal
Studying the Relationship between the Traffic Flow Structure, the Traffic Capacity of Intersections, and Vehicle-Related Emissions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Domain Decomposition Method for Approximation the Spectrum of the Diffusion Operator in a Domain Containing Thin Tubes

by
Andrey Amosov
1,*,
Delfina Gómez
2,
Grigory Panasenko
3 and
Maria-Eugenia Pérez-Martinez
4
1
Department of Mathematical and Computer Modelling, National Research University “Moscow Power Engineering Institute”, Krasnokazarmennay St. 14, 111250 Moscow, Russia
2
Departamento de Matemáticas, Estadística y Computación, Universidad de Cantabria, Av. de los Castros s/n., 39005 Santander, Spain
3
Institute of Applied Mathematics, Vilnius University, Naugarduko Str., 24, 03225 Vilnius, Lithuania
4
Departamento Mathematica Aplicada y Ciencias de la Computación, Universidad de Cantabria, Av. de los Castros s/n., 39005 Santander, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(16), 3592; https://doi.org/10.3390/math11163592
Submission received: 20 June 2023 / Revised: 4 August 2023 / Accepted: 16 August 2023 / Published: 19 August 2023
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
The spectral problem for the diffusion operator is considered in a domain containing thin tubes. A new version of the method of partial asymptotic decomposition of the domain is introduced to reduce the dimension inside the tubes. It truncates the tubes at some small distance from the ends of the tubes and replaces the tubes with segments. At the interface of the three-dimensional and one-dimensional subdomains, special junction conditions are set: the pointwise continuity of the flux and the continuity of the average over a cross-section of the eigenfunctions. The existence of the discrete spectrum is proved for this partially reduced problem of the hybrid dimension. The conditions of the closeness of two spectra, i.e., of the diffusion operator in the full-dimensional domain and the partially reduced one, are obtained.

1. Introduction

The method of asymptotic partial decomposition of the domain (MAPDD) was introduced for partial derivative equations set in thin tube structures in [1] (cf. also [2,3,4]). The thin tube structures (also called thin rod structures) are some unions of thin cylinders. In this case, the method gives an important gain in computational resources reducing the dimension to one everywhere except for small two- or three-dimensional parts of the domain. This method was applied to the spectral problems in [5], where the asymptotic of the spectrum of the Laplacian in two joint thin rectangles with Neumann’s boundary conditions on the lateral boundary were considered. The MAPDD was justified for such a structure. The asymptotic behavior of the spectrum of the Laplacian with Dirichlet’s or Neumann’s boundary conditions in thin domains was considered in a vast body of literature (see [6,7,8,9,10,11,12,13], and the references therein). Furthermore, let us mention some works on the eigenvalues for the Laplacian for several special cases of domains containing thin tubes, described with the help of a small parameter (cf. [14,15,16,17] among others); these spectral problems were treated using asymptotic analysis. As a matter of fact, Refs. [6,7,8,13] address spectral problems in thin planar domains, while Refs. [10,11,12] deal with three-dimensional domains in very different structures from those considered here. This topic is of great interest due to its multiple applications in scattering theory, wave-guides, etc.
In the present paper, we introduce a new version of the MAPDD with special interface conditions between the full-dimensional and one-dimensional parts, that is, the pointwise continuity of the flux and the continuity of the average over a cross-section of the eigenfunctions. These conditions were used for MAPDD approximation of the heat equation in [18]. In particular, some numerical methods also use junction conditions with jumps of the unknown function (see, e.g., [19,20]). Here, we deal with a spectral problem for the Neumann diffusion operator.
The results obtained in former works on the asymptotic analysis of spectral problems in thin domains mostly use the presence in the model of a small parameter, that is, the ratio of the characteristic sizes in the transverse and longitudinal (or in-plane) directions. The justification for these results was provided via the proof of theorems concerning the convergence of the spectrum of the original problem to that of the spectrum of the reduced problem as the small parameter tends to zero. In the present paper, we use another approach related to the method of justification for the MAPDD in [21], but we introduce different junction conditions. There is no explicit small parameter in the description of the domain, but, implicitly, it is introduced via the assumption that the first positive eigenvalues of the Neumann diffusion operator on the cross-sections of the tubes are sufficiently large. The closeness of the spectra of the original and reduced models is proved under the condition that the eigenvalues of both the original and the reduced models are smaller than these first eigenvalues on the cross-sections of the tubes. In addition, here, we do not assume that the domain is thin everywhere. It may have some thick parts connected by thin tubes. The boundary is assumed to be Lipschitz, and we do not require the regularity of the coefficient in the diffusion operator out of the tubes. The proof of the aforementioned closeness of the spectra also relies on the distance from the ends of the tubes to the one-dimensional domain in the reduced model, a distance that needs to be adjusted to achieve the desired accuracy for the approximation between the spectra.
The main results are as follows. The spectral problem for the diffusion operator is considered in a domain containing thin tubes. The reduced model is obtained from the original one by the truncation of the three-dimensional tubes at some small distance from the ends of the tubes replacing the truncated parts of the tubes by the segments. At the interface of the three-dimensional and one-dimensional subdomains, special junction conditions are set: the pointwise continuity of the flux and continuity of the average of the eigenfunctions over a cross-section. The existence of the discrete spectrum is proved for this partially reduced problem of the hybrid dimension. After prescribing an accurate precision ε , to obtain it, the conditions of the closeness of two spectra (i.e., of the diffusion operator in the full-dimensional domain and the partially reduced one) are obtained.
The structure of the paper is as follows: Section 2 contains some preliminary results, the definition of the domain, and the setting of the spectral problem. Moreover, the required properties on the smoothness of the eigenfunctions are obtained (cf. Section 2.3). In Section 3, the reduced approximate spectral problem is formulated. Because of the hybrid dimension, notations, weak formulation, and properties of spaces become more complicated. The required smoothness results for the eigenfunctions are also proved (cf. Section 3.3). In Section 4, we state the main results on the closeness of the spectra of the original and reduced problems (cf. Theorems 2 and 3). As a consequence, we claim the closeness of the spectra in the case where the tubes are cylinders (cf. Theorems 4 and 5). Section 5 and Section 6 contain the proofs of Theorems 2 and 3, respectively. Finally, some concluding remarks are outlined in Section 7.

2. Preliminary Results and Setting of the Problem

In this section, we introduce some notations and preliminary results that will be used throughout the whole paper. Section 2.2 contains the setting of the spectral problem under consideration, and, in Section 2.3, we provide certain smoothness properties for the eigenfunctions.

2.1. Preliminaries

By a tube, we mean a set T R 3 of the form T = Ω × ( 0 , ) in the local coordinates ( y , z ) = ( y 1 , y 2 , z ) , where the base Ω of the tube is a bounded domain in R 2 with a Lipschitz boundary Ω . The set T = Ω × ( 0 , ) is called the lateral surface of the tube T , and stands for its length. Define
D 1 φ = φ y 1 , D 2 φ = φ y 2 , D z φ = φ z , y φ = ( D 1 φ , D 2 φ ) , Δ y φ = D 1 2 φ + D 2 2 φ .
Let μ 1 ( Ω ) be the first positive eigenvalue of the following problem
Δ y φ = μ φ , y Ω , y φ · ν = 0 , y Ω ,
where ν is the outward normal to Ω .
Recall that for a disk Ω with a diameter D, μ 1 ( Ω ) = 4 p 2 / D 2 where p = 1.8142 is the smallest positive root of the derivative J 1 of the Bessel function J 1 . In the general case, we have [22]
μ 1 ( Ω ) μ 1 ( Ω * ) = π p 2 meas ( Ω ) ,
where Ω * is a disk of the same measure as Ω . Inequality (1) becomes an equality only if Ω is the disk.
In the case where Ω is a convex domain of diameter D, we have the lower estimate [23]:
π 2 2 D 2 μ 1 ( Ω ) .
Note that ρ ( Ω ) = 1 / μ 1 ( Ω ) is an optimal constant in the Poincaré inequality
φ L 2 ( Ω ) ρ ( Ω ) y φ L 2 ( Ω ) φ H ˜ 1 ( Ω ) .
Here, H ˜ 1 ( Ω ) = { φ H 1 ( Ω ) | Ω φ d x = 0 } is a Hilbert space equipped with the inner product and norm, respectively
( φ , ψ ) H ˜ 1 ( Ω ) = ( y φ , y ψ ) L 2 ( Ω ) , φ H ˜ 1 ( Ω ) = y φ L 2 ( Ω ) .
Obviously, μ 1 ( ϵ Ω ) = ϵ 2 μ 1 ( Ω ) for the contracted domain, and ϵ Ω = { ϵ y | y Ω } for all 0 < ϵ .
We see that for a wide class of domains, the value ρ ( Ω ) characterizes the size and geometry of Ω . So, we refer to a thin tube T when ρ ( Ω ) . For an Ω that is a disk with a radius of ϵ , the tube is thin if ϵ ; in this case, ρ ( Ω ) = ϵ / p .

2.2. Spectral Problem

Let G be a bounded domain in R 3 with a Lipschitz boundary G . Consider the following spectral problem: find a couple ( λ , u ) R × H 1 ( G ) such that u L 2 ( G ) = 1 , satisfying the equation
div ( K u ) = λ u , x G
and the boundary condition
K u · n = 0 , x G .
Here, K L ( G ) , ess inf x G K ( x ) > 0 , n is the outward normal to G .
The weak formulation is given by the following identity:
( K u , φ ) L 2 ( G ) = λ ( u , φ ) L 2 ( G ) φ H 1 ( G ) .
It is well known that this problem has a countable set 0 = λ 0 < λ 1 λ 2 of eigenvalues and that the corresponding set { φ k } k = 0 of eigenfunctions can be chosen to form an orthogonal base in H 1 ( G ) and an orthonormal base in L 2 ( G ) .
Assume that G contains a set of disjoint thin tubes T j , 1 j N such that their lateral surfaces T j belong to the boundary G . Let each tube T j in the local coordinate system ( y , z ) have the form T j = Ω j × ( 0 , j ) and let T j = Ω j × ( 0 , j ) . Consider the subtube T j , h j = Ω j × ( h j , j h j ) , where h j ( 0 , j / 2 ) , and denote T h = N j = 1 T j , h j , G h = G T ¯ h .
Define the average of the function u L 2 ( Ω j ) over the cross-section of the tube T j
[ u ] j = 1 meas ( Ω j ) Ω j u ( y ) d y .
Assume that K satisfies the following condition:
K ( x ) = K j = const in T j , 1 j N .

2.3. Some Properties of Eigenfunctions

Let u be an eigenfunction corresponding to the eigenvalue λ . Denote by u j the restriction of the function u to T j , 1 j N . Let r j = u j [ u j ] j . We will consider u j and r j as functions on ( 0 , j ) with values in H 1 ( Ω j ) . Note that u j , r j L 2 ( 0 , j ; H 1 ( Ω j ) ) , D z u j , D z r j L 2 ( 0 , j ; L 2 ( Ω j ) ) . We set λ ˜ j = λ / K j .
Lemma 1.
The function [ u j ] j belongs to the space H 2 ( 0 , j ) , and
D z 2 [ u j ] j = λ ˜ j [ u j ] j .
Proof. 
It is clear that [ u j ] j H 1 ( 0 , j ) and D z [ u j ] j = [ D z u j ] j .
Let η C 0 ( 0 , j ) . Substituting into (2) the function φ , equal to K j 1 η ( z ) on T j and equal to zero on G T j , we obtain the identity
0 j D z [ u j ] j ( z ) η ( z ) d z = 0 j λ ˜ j [ u j ] j ( z ) η ( z ) d z η C 0 ( 0 , j ) ,
which implies the existence of the weak derivative D z 2 [ u j ] j by definition and equality (3).
The lemma is proved. □
Lemma 2.
The functions u j and r j have the derivatives D z 2 u j and D z 2 r j , which belong to the space L l o c 2 ( 0 , j ; L 2 ( Ω j ) ) .
Proof. 
We introduce the finite-difference analogues of the derivative D z Φ
h Φ ( z ) = Φ ( z + h ) Φ ( z ) h , ¯ h Φ ( z ) = Φ ( z ) Φ ( z h ) h ,
where h > 0 .
Assume that ( a , b ) ( 0 , j ) , and let ζ C 0 ( 0 , j ) be a cut-off function, equal to one on ( a , b ) . By substituting into identity (2) the test function φ equal to K j 1 ¯ h ( ζ 2 h u j ) on T j and zero on G T j , we obtain
0 j D z u j , ¯ h ( ζ 2 h D z u j ) L 2 ( Ω j ) d z 0 j D z u j , ¯ h ( 2 ζ ζ h u j ) L 2 ( Ω j ) d z 0 j y u j , ¯ h ( ζ 2 h y u j ) L 2 ( Ω j ) d z = λ ˜ j 0 j u j , ¯ h ( ζ 2 h u j ) L 2 ( Ω j ) d z .
Here, h is small enough.
Using the finite-difference analogue of the formula for integration by parts, we obtain
0 j ζ h D z u j L 2 ( Ω j ) 2 d z + 0 j ζ h y u j L 2 ( Ω j ) 2 d z = 0 j ( ζ h D z u j , 2 ζ h u j ) L 2 ( Ω j ) d z + λ ˜ j 0 j ζ h u j L 2 ( Ω j ) 2 d z 1 2 0 j ζ h D z u j L 2 ( Ω j ) 2 d z + 1 2 0 j 2 ζ h u j L 2 ( Ω j ) 2 d z + λ ˜ j 0 j ζ h u j L 2 ( Ω j ) 2 d z .
From here,
0 j ζ h D z u j L 2 ( Ω j ) 2 d z c ζ D z u j L 2 ( 0 , j ; L 2 ( Ω j ) ) 2 .
Thus, we have the following estimate uniform in h
h D z u j L 2 ( a , b ; L 2 ( Ω j ) ) 2 c ζ D z u j L 2 ( 0 , j ; L 2 ( Ω j ) ) 2 ,
where c ζ is a constant depending on ζ .
It implies the existence of the derivative D z 2 u j L 2 ( a , b ; L 2 ( Ω j ) ) . Thus, D z 2 u j L l o c 2 ( 0 , j ; L 2 ( Ω j ) ) . As a consequence, there exists a derivative D z 2 r j L l o c 2 ( 0 , j ; L 2 ( Ω j ) ) .
The lemma is proved. □
Lemma 3.
For almost all z ( 0 , j ) the following identities hold
( D z 2 u j ( z ) , φ ) L 2 ( Ω j ) = ( y u j ( z ) , y φ ) L 2 ( Ω j ) λ ˜ j ( u j ( z ) , φ ) L 2 ( Ω j ) φ H 1 ( Ω j ) .
( D z 2 r j ( z ) , φ ) L 2 ( Ω j ) = ( y r j ( z ) , y φ ) L 2 ( Ω j ) λ ˜ j ( r j ( z ) , φ ) L 2 ( Ω j ) φ H 1 ( Ω j ) .
Proof. 
Let η C 0 ( 0 , j ) and φ H 1 ( Ω j ) . Since D z 2 u j L l o c 2 ( 0 , j ; L 2 ( Ω j ) ) , then the substitution into (2) of the test function equal to K j 1 η ( z ) φ ( y ) on T j and zero on G T j and an integration by parts gives
0 j ( D z 2 u j ( z ) , φ ) L 2 ( Ω j ) η ( z ) d z + 0 j ( u j ( z ) , φ ) L 2 ( Ω j ) η ( z ) d z = λ ˜ j 0 j ( u ( z ) , φ ) L 2 ( Ω j ) η ( z ) d z .
This implies identity (5). As a consequence (cf. (4)), identity (6) holds.
The lemma is proved. □

3. Approximate Spectral Problem

This section is devoted to the setting of the problem referred to by us as the approximate spectral problem of hybrid dimensions. Section 3.1. contains the notations and preliminary results convenient for the formulation of the problem. In Section 3.2, we formulate the spectral problem along with its variational formulation in the suitable Hilbert spaces and show the discreteness of the spectrum. Moreover, we prove the required smoothness results for the eigenfunctions which, due to the junctions conditions (cf. (15) and (16)), can be weaker than those obtained for the eigenfunctions of the original problem in Section 2.3, but which somehow justify the junction conditions (cf. Corollaries 2 and 3).

3.1. Spaces L h 2 ( G ) , L ˜ h 2 ( G ) , H h 1 ( G ) , and H ˜ h 1 ( G )

Henceforth, we will use the notation u j for the restriction of a function u L 2 ( G ) to T j , 1 j N .
Let L h 2 ( G ) be the close subspace of L 2 ( G ) of functions u such that u j ( y , z ) = [ u j ] j ( z ) for x = ( y , z ) T j , h j , 1 j N . Introduce also the space
L ˜ h 2 ( G ) = { u L h 2 ( G ) | G u d x = 0 } .
As noted in the introduction, the approximate spectral problem deals with the domain G with truncated tube-like parts. The following lemma pertains to the topological structure of the set G h = G N j = 1 T j , h j ¯ .
Lemma 4.
The set G h has the following structure:
G h = m k = 1 G h , k ,
where m N + 1 and G h , k are domains with Lipschitz boundaries such that G h , k ¯ G h , ¯ = for k .
Proof. 
Let G ˜ be an open subset of the set G. We will write that points x , y G ˜ are connected in G ˜ if there exists a continuous curve L G ˜ starting at x and ending at y or y = x .
We set G h 0 = G and G h j = G h j 1 T j , h j ¯ , 1 j N . Note that the boundary G h j of each of the sets G h j satisfies the Lipschitz condition because G h j = G h j 1 ( Ω j × ( 0 , j ) ) ( Ω j × { 0 , j } ) .
Let us show that for all 1 j N the following equality holds:
G h j = m j k = 1 G h , k ,
where m j j + 1 and G h , k are domains with Lipschitz boundaries such that G h , k ¯ G ¯ h , = for k .
Introduce the sets T j , h j and T j , h j r , 1 j N , which in the local coordinate system associated with T j have the following form:
T j , h j = Ω j × ( 0 , h j ) , T j , h j r = Ω j × ( j h j , j ) .
Step 1. Let us show that equality (8) is valid for j = 1 . If the set G h 1 = G T 1 , h 1 ¯ is connected, then (8) is true with m 1 = 1 and G h , 1 = G h 1 .
Assume that G h 1 is not connected. Fix points x T 1 , h 1 and x r T 1 , h 1 r . Denote by G h 1 , the set of points y G h 1 , connected with x in G h 1 . Similarly, denote by G h 1 , r the set of points y G h 1 , connected with x r in G h 1 . It is clear that G h 1 , and G h 1 , r are disjoint domains and T 1 , h 1 G h 1 , , T 1 , h 1 r G h 1 , r .
Let us prove that
G h 1 = G h 1 , G h 1 , r
Assume that y G h 1 G h 1 , . Since set G is connected, there is a curve L = { x = ζ ( t ) , 0 t 1 } G , where the function ζ : [ 0 , 1 ] G is continuous, and ζ ( 0 ) = y , ζ ( 1 ) = x .
Since y G h 1 , , then L T 1 , h 1 ¯ , and there exists t * ( 0 , 1 ) such that ζ ( t * ) T ¯ 1 , h 1 and ζ ( t ) G h 1 for 0 t < t * . Moreover, there exists t ( 0 , t * ) such that ζ ( t ) T 1 , h 1 T 1 , h 1 r . If ζ ( t ) T 1 , h 1 then points y and x are connected in G h 1 , which contradicts the assumption y G h 1 , . Hence, ζ ( t ) T 1 , h 1 r , and so y G h 1 , r . Thus, equality (9) is true, i.e., (8) is true with m 1 = 2 and G h , 1 = G h 1 , , G h , 2 = G h 1 , r .
Note also that G h 1 = G h 1 , G h 1 , r , where G h 1 , G h 1 , r = . Since G h 1 satisfies the Lipschitz condition, then G h 1 , and G h 1 , r satisfy the Lipschitz condition also.
Step n. Assume that (8) is true for j = n 1 . Since the set T j , h j G h j 1 is connected, then T j , h j G h , k for some 1 k m n 1 . If G ˜ h , k = G h , k T j , h j ¯ is connected then
G h n = 1 i m n 1 , i k G h , i G ˜ h , k .
If the set G ˜ h , k is not connected, then repeating the argument at step 1 (with G h 1 replaced by G ˜ h , k and T 1 , h 1 , T 1 , h 1 r replaced by T j , h j , T j , h j r ), we see that G ˜ h , k is represented as the union of two non-intersecting regions G ˜ h , k and G ˜ h , k r . In this case
G h n = 1 i m n 1 , i k G h , i G ˜ h , k G ˜ h , k r .
Thus, representation (8) also holds for j = n .
Having done N steps, we claim equality (7).
The lemma is proved. □
Introduce the space H h 1 ( G ) consisting of functions u L h 2 ( G ) such that u H 1 ( G h , k ) for all 1 k m , u j H 1 ( T j , h ) for all 1 j N , and the following junction conditions
[ u ] j | z = h j 0 = u j | z = h j + 0 , u j | z = j h j 0 = [ u ] j | z = j h j + 0
hold.
Naturally, u j = [ u j ] j H 1 ( 0 , j ) and u j ( x ) = ( 0 , 0 , D z u j ( z ) ) for x = ( y , z ) T j , h , 1 j N .
Note that H h 1 ( G ) is a Hilbert space with the inner product and the norm
( u , v ) H h 1 ( G ) = ( K u , v ) L 2 ( G ) + ( u , v ) L 2 ( G ) , u H h 1 ( G ) = K u L 2 ( G ) 2 + u L 2 ( G ) 2 1 / 2 .
Lemma 5.
The embedding of H h 1 ( G ) into L h 2 ( G ) is compact.
Proof. 
Assume that { u n } n = 1 H h 1 ( G ) is a sequence that converges weakly in H h 1 ( G ) . As a result, it weakly converges in H 1 ( G h , k ) for all 1 k m and in H 1 ( T j , h ) for all 1 j N . So, { u n } n = 1 converges strongly in L 2 ( G h , k ) for all 1 k m and in L 2 ( T j , h ) for all 1 j N . Thus, it converges strongly in L h 2 ( G ) .
The lemma is proved. □
Lemma 6.
Suppose that u H h 1 ( G ) and u = 0 . Then, u = const .
Proof. 
It follows from u = 0 that u ( x ) = c h , k = const in G h , k for all 1 k m and u ( x ) = c j = const in T j , h j for all 1 j m . Hence, due to the connectedness of G and condition (10), it implies that u ( x ) = const in G.
The lemma is proved. □
The following analogue of the Poincaré inequality holds in H h 1 ( G ) (see analogous inequalities in some spaces of discontinuous functions in [19,24]).
Lemma 7.
There exists a constant C > 0 such that
u H h 1 ( G ) 2 C [ ( K u , u ) L 2 ( G ) + ( G u d x ) 2 ] u H h 1 ( G ) .
Proof. 
Assume the opposite. Then, there is a sequence { u n } n = 1 H h 1 ( G ) such that
u n H h 1 ( G ) 2 > n [ ( K u n , u n ) L 2 ( G ) + ( G u n d x ) 2 ] n 1 .
Putting v n = 1 u n H h 1 ( G ) u n , we have v n H h 1 ( G ) = 1 and
1 n > ( K v n , v n ) L 2 ( G ) + ( G v n d x ) 2 n 1 .
Using the compactness of embedding of H h 1 ( G ) into L h 2 ( G ) , we choose a subsequence { v n k } k = 1 such that v n k v weakly in H h 1 ( G ) and strongly in L h 2 ( G ) .
It follows from (12) that v n k L 2 ( G ) 0 and G v n k d x 0 . So, v = 0 , v = c = const and G v d x = c · meas G = 0 . Thus, v = 0 and v n k 0 in H h 1 ( G ) which contradicts the equality v n k H h 1 ( G ) = 1 .
The lemma is proved. □
Let us introduce in H h 1 ( G ) the closed subspace
H ˜ h 1 ( G ) = H h 1 ( G ) L ˜ h 2 ( G ) = { u H h 1 ( G ) | G u d x = 0 } .
It follows from (11) that H ˜ h 1 ( G ) is a Hilbert space with the following inner product and the norm:
( u , v ) H ˜ h 1 ( G ) = ( K u , v ) L 2 ( G ) , u H ˜ h 1 ( G ) = K u L 2 ( G ) .
Moreover, it follows from Lemma 5 that the embedding of H ˜ h 1 ( G ) into L ˜ h 2 ( G ) is compact.

3.2. Setting of the Approximate Spectral Problem

Consider the following spectral problem: find a couple ( λ ^ , U ) × H h 1 ( G ) , such that U L 2 ( G ) = 1 , satisfying the equation
div ( K Δ U ) = λ ^ U , x G h T h ,
the boundary condition
K U · n = 0 , x G
and the following junction conditions:
[ U j ] | z = h j 0 = U j | z = h j + 0 , U j | z = j h j 0 = [ U j ] | z = j h j + 0 , 1 j N ,
D z U j | z = h j 0 = D z U j | z = h j + 0 , D z U j | z = j h j 0 = D z U j | z = j h j + 0 , 1 j N .
where U j is the restriction of U on T j .
The weak formulation of this problem is given by the following identity:
( K U , φ ) L 2 ( G ) = λ ^ ( U , φ ) L 2 ( G ) φ H h 1 ( G ) .
Note that conditions (15) are satisfied since U H h 1 ( G ) . We will show also that (17) implies a validity of junction conditions (16) (see Corollaries 2 and 3).
It is clear that each eigenvalue λ ^ is non-negative, the minimal eigenvalue is λ ^ 0 = 0 , and the corresponding eigenfunction U ^ 0 is a constant.
Since the embedding of H ˜ h 1 ( G ) into L ˜ h 2 ( G ) is dense and compact, we are in a classical abstract framework of bilinear, continuous, coercive forms on a couple of Hilbert spaces H ˜ h 1 ( G ) and L ˜ h 2 ( G ) , cf., for instance, Theorem 5.5 of Chapter I in [25]. So, the following theorem is true.
Theorem 1.
There exists a system { φ ^ k } k = 1 of eigenfunctions to problem (17) corresponding to the eigenvalues 0 < λ ^ 1 λ ^ 2 λ ^ k forming an orthonormal base in H ˜ h 1 ( G ) and an orthogonal base in L ˜ h 2 ( G ) .
Corollary 1.
There is a system { φ ^ k } k = 0 of eigenfunctions to problem (17) corresponding to the eigenvalues 0 = λ ^ 0 < λ ^ 1 λ ^ 2 λ ^ k , forming an orthonormal base in L h 2 ( G ) and an orthogonal base in H h 1 ( G ) .

3.3. Some Properties of Eigenfunctions

Henceforth, in this section, λ ^ is an eigenvalue and U H h 1 ( G ) is corresponding eigenfunction such that U L 2 ( G ) = 1 . Moreover, λ ¯ j = λ ^ / K j , 1 j N .
Lemma 8.
The function U satisfies the identity
( K U , φ ) L 2 ( G ) = λ ^ ( U , φ ) L 2 ( G ) φ H 1 ( G ) .
Proof. 
Let φ H 1 ( G ) . Note that the function
φ ˜ ( x ) = φ ( x ) , x G h , φ ] j ( x ) , x T j , h j , 1 j N
belongs to H h 1 ( G ) and
( K U , φ ) L 2 ( G ) λ ^ ( U , φ ) L 2 ( G ) = ( K U , φ ˜ ) L 2 ( G ) λ ^ ( U , φ ˜ ) L 2 ( G )
So, (18) holds.
The lemma is proved. □
Remember that U j is the restriction of U on T j , 1 j m . We will consider U j as a function on ( 0 , j ) with values in H 1 ( Ω j ) . We put W j ( z ) = D z U j ( z ) for z ( 0 , h j ) ( j h j , j ) and W j ( z ) = D z [ U j ] j ( z ) for z ( h j , j h j ) . It is clear that U j L 2 ( 0 , j ; H 1 ( Ω j ) ) and W j L 2 ( 0 , j ; L 2 ( Ω j ) ) . Moreover, [ U j ] j L 2 ( 0 , j ) and D z [ U j ] j L 2 ( 0 , j ) .
Lemma 9.
The function W j has a derivative D z W j L 2 ( 0 , j ; ( H 1 ( Ω j ) ) ) .
Proof. 
It follows from (17) that
0 j ( W j ( z ) , φ ) L 2 ( Ω j ) η ( z ) d z + 0 j ( y U j ( z ) , y φ ) L 2 ( Ω j ) λ ¯ j ( U j ( z ) , φ ) L 2 ( Ω j ) η ( z ) d z = 0
for all φ H 1 ( Ω j ) and η C 0 ( 0 , j ) . Remember that λ ¯ j = λ ^ / K j
Therefore, the function W j has a derivative D z W j : ( 0 , j ) ( H 1 ( Ω j ) ) such that
D z W j ( z ) , φ ( H 1 ( Ω j ) ) × H 1 ( Ω j ) = ( y U j ( z ) , y φ ) L 2 ( Ω j ) λ ¯ j ( U j ( z ) , φ ) L 2 ( Ω j ) φ H 1 ( Ω j ) .
for almost all z ( 0 , j ) .
It follows from the estimate
D z W j ( z ) , φ ( H 1 ( Ω j ) ) × H 1 ( Ω j ) y U j ( z ) L 2 ( Ω j ) y φ L 2 ( Ω j ) + λ ¯ j U j ( z ) L 2 ( Ω j ) φ L 2 ( Ω j )                                                 max { 1 , λ ¯ j } U j ( z ) H 1 ( Ω j ) φ H 1 ( Ω j )
that
D z W j L 2 ( 0 , j ; ( H 1 ( Ω j ) ) ) max { 1 , λ ¯ j } U j L 2 ( 0 , j ; H 1 ( Ω j ) ) .
Thus, D z W j L 2 ( 0 , j ; ( H 1 ( Ω j ) ) ) .
The lemma is proved. □
Corollary 2.
W j C ( [ 0 , j ] ; ( H 1 ( Ω j ) ) ) , and so, the eigenfunction U satisfies the junction conditions (16) in ( H 1 ( Ω j ) ) .
Lemma 10.
The function [ U j ] j belongs to the space H 2 ( 0 , j ) , and
D z 2 [ U j ] j = λ ¯ j [ U j ] j .
Proof. 
The substitution of φ = K j 1 η with η C 0 ( 0 , j ) into (17) gives
0 j D z [ U j ] j ( z ) η ( z ) d z = λ ¯ j 0 j [ U j ] j ( z ) η ( z ) d z = 0 η C 0 ( 0 , j ) .
So, D z 2 [ U j ] j = λ ¯ j [ U ] j .
The lemma is proved. □
We put R j = U j [ U j ] j , 1 j N and will consider R j as a function on ( 0 , j ) with values in H 1 ( Ω j ) . It is clear that R j = 0 , D z R j = 0 on ( h j , j h j ) and D z R j = W j D z [ U j ] j on ( 0 , h j ) ( j h j , j ) .
It follows from Lemmas 9 and 10 that D z R j C ( [ 0 , j ] ; ( H 1 ( Ω j ) ) ) . So D z R j ( h j ) = 0 and D z R j ( j h j ) = 0 in ( H 1 ( Ω j ) ) .
Lemma 11.
For each 0 < a < h j the function R j has derivatives D z 2 R j L 2 ( a , h j ; L 2 ( Ω j ) ) and D z 2 R j L 2 ( j h j , j a ; L 2 ( Ω j ) ) .
Proof. 
We put R ˜ j ( z ) = R j ( z ) for 0 < z h j and R ˜ j ( z ) = R j ( 2 h j z ) for h j z < 2 h j . Note that R ˜ j C ( 0 , 2 h j ; L 2 ( Ω j ) ) and D z R ˜ j L 2 ( 0 , 2 h j ; L 2 ( Ω j ) ) .
Let φ L 2 ( 0 , 2 h j ; H 1 ( Ω j ) ) be an arbitrary function identically equal to zero in some neighborhoods of the points z = 0 and z = 2 h j and such that D z φ L 2 ( 0 , 2 h j ; L 2 ( Ω j ) ) .
As R j ( z ) = 0 for z ( h j , j h j ) then (17) and (20) imply that
0 h j ( D z R j , D z φ ) L 2 ( Ω j ) d z + 0 h j ( y R j , y φ ) L 2 ( Ω j ) d z = λ ¯ j 0 h j ( R j , φ ) L 2 ( Ω j ) d z .
Moreover,
h j 2 h j ( D z R ˜ j , D z φ ˜ ) L 2 ( Ω j ) d z + h j 2 h j ( y R ˜ j , y φ ˜ ) L 2 ( Ω j ) d z = λ ¯ j h j 2 h j ( R ˜ j , φ ˜ ) L 2 ( Ω j ) d z ,
where φ ˜ ( z ) = φ ( 2 h j z ) . Consequently,
0 2 h j ( D z R ˜ j , D z φ ) L 2 ( Ω j ) d z + 0 2 h j ( y R ˜ j , φ ) L 2 ( Ω j ) d z = λ ¯ j 0 2 h j ( R ˜ j , φ ) L 2 ( Ω j ) d z .
Assume that 0 < a < h j < b < 2 h j and ζ C 0 ( 0 , 2 h j ) is a cutoff function, equal to one on ( a , b ) . Substituting into (21) the test function φ = ¯ h ( ζ 2 h R ˜ j ) , where h is small enough, we have
0 2 h j ( D z R ˜ j , ¯ h ( ζ 2 h D z R ˜ j ) ) L 2 ( Ω j ) d z 0 2 h j ( D z R ˜ j , ¯ h ( 2 ζ ζ h R ˜ j ) ) L 2 ( Ω j ) d z 0 2 h j ( y R ˜ j , y ¯ h ( ζ 2 h R ˜ j ) ) L 2 ( Ω j ) d z = λ ¯ j 0 2 h j ( R ˜ j , ¯ h ( ζ 2 h R ˜ j ) ) L 2 ( Ω j ) d z .
Using the difference analogue of the formula for integration by parts, we obtain
0 2 h j ζ h D z R ˜ j L 2 ( Ω j ) 2 d z + 0 2 h j ζ h y R ˜ j L 2 ( Ω j ) 2 d z = 0 2 h j ( ζ h D z R ˜ j , 2 ζ h R ˜ j ) L 2 ( Ω j ) d z + λ ˜ j 0 2 h j ζ h R ˜ j L 2 ( Ω j ) 2 d z 1 2 0 2 h j ζ h D z R ˜ j L 2 ( Ω j ) 2 d z + 1 2 0 2 h j 2 ζ h R ˜ j L 2 ( Ω j ) 2 d z + λ ¯ j 0 2 h j ζ h R ˜ j L 2 ( Ω j ) 2 d z .
From here, we have the following estimate uniform in h
h D z R ˜ j L 2 ( a , b ; L 2 ( Ω j ) ) 2 c ζ D z R ˜ j L 2 ( 0 , 2 h j ; L 2 ( Ω j ) ) 2 = 2 c ζ D z R j L 2 ( 0 , h j ; L 2 ( Ω j ) ) 2 .
This implies the existence of the derivative D z 2 R ˜ j L 2 ( a , b ; L 2 ( Ω j ) ) . Thus, D z 2 R j L 2 ( a , h j ; L 2 ( Ω j ) ) .
The existence of the derivative D z 2 R j L 2 ( j h j , j a ; L 2 ( Ω j ) ) is proved in a similar way.
The lemma is proved. □
Corollary 3.
The following properties hold:
(1) 
D z R j C ( ( 0 , h j ] ; L 2 ( Ω j ) ) and D z R j ( h j 0 ) = 0 in L 2 ( Ω j ) .
(2) 
D z R j C ( [ j h j , j ) ; L 2 ( Ω j ) ) and D z R j ( j h j + 0 ) = 0 in L 2 ( Ω j ) .
(3) 
The eigenfunction U satisfies the junction conditions in (16) in L 2 ( Ω j ) .
Proof. 
To prove properties (1) and (2), it should be taken into account that, due to Corollary 2, D z R j ( h j 0 ) = 0 , D z R j ( j h j + 0 ) = 0 in ( H 1 ( Ω j ) ) , and H 1 ( Ω j ) is dense in L 2 ( Ω j ) .
Property (3) holds because D z U j = D z [ U j ] j + D z R j , where D z R j ( z ) = 0 for z ( h j , j h j ) and D z [ U j ] j C [ 0 , j ] .
The corollary is proved. □
Lemma 12.
For almost all z ( 0 , h j ) ( j h j , j ) , the following identities hold
( D z 2 U j ( z ) , φ ) L 2 ( Ω j ) = ( y U j ( z ) , y φ ) L 2 ( Ω j ) λ ¯ j ( U j ( z ) , φ ) L 2 ( Ω j ) φ H 1 ( Ω j ) .
( D z 2 R j ( z ) , φ ) L 2 ( Ω j ) = ( y R j ( z ) , y φ ) L 2 ( Ω j ) λ ¯ j ( R j ( z ) , φ ) L 2 ( Ω j ) φ H 1 ( Ω j ) .
The proof is similar to the proof of Lemma 3.

4. Formulation of the Main Results

Let { λ k } k = 0 be the set of eigenvalues to problem (2), and { λ ^ k } k = 0 be the set of eigenvalues to problem (17).
Let ε be the desired accuracy, 0 < ε < 1 and ε 1 = ε 1 + ε .
Let 0 = μ 0 , j < μ 1 , j μ 2 , j be the set of eigenvalues to the Neumann problem
Δ y ψ = μ ψ , y Ω j , y ψ · ν = 0 , y Ω j ,
and { ψ k , j } k = 0 be the set of corresponding eigenfunctions, forming an orthogonal base in H 1 ( Ω j ) and an orthonormal base in L 2 ( Ω j ) .
Theorem 2.
Let λ be the eigenvalue of problem (2) and λ ˜ j = λ / K j < μ 1 , j for all 1 j N . Let the numbers h j be such that
1 ν 1 , j max ln 2 ε 1 2 , ln 96 2 K j 2 ν 1 , j 4 ε 1 2 , ln 24 2 K j λ ν 1 , j 2 ε 1 2 h j < j / 2 , 1 j N ,
where ν 1 , j = μ 1 , j λ ˜ j . Then, the following estimate holds:
min k 0 | λ ^ k λ | ε .
Theore 3.
Let λ ^ be the eigenvalue of problem (17) and λ ¯ j = λ ^ / K j < μ 1 , j for all 1 j N . Let the numbers h j be such that
1 ν ^ 1 , j max ln 4 ε 1 2 , ln 96 2 K j 2 ν ^ 1 , j 4 ε 1 2 , ln 12 2 K j λ ^ ν ^ 1 , j 2 ε 1 2 h j < j / 2 , 1 j N ,
where ν ^ 1 , j = μ 1 , j λ ¯ j . Then, the following estimate holds:
min k 0 | λ k λ ^ | ε .
Consider a special case where, for each of the tubes T j , its cross section is a disk with a radius of ϵ j . Remember that, in this case, μ 1 , j = p 2 / ϵ j 2 , where p = 1.8142 is the smallest positive root of the derivative J 1 of the Bessel function J 1 . So, Theorems 2 and 3 can be reformulated as follows.
Theorem 4.
Let λ be the eigenvalue of problem (2) and λ < K j p 2 / ϵ j 2 for all 1 j N . Let the numbers h j be such that
K j ρ j max ln 2 ε 1 2 , ln 96 2 ρ j 2 ε 1 2 , ln 24 2 λ ρ j ε 1 2 h j < j / 2 , 1 j N ,
where ρ j = K j p 2 / ϵ j 2 λ . Then, the following estimate holds:
min k 0 | λ ^ k λ | ε .
Theorem 5.
Let λ ^ be the eigenvalue of problem (17) and λ ^ < K j p 2 / ϵ j 2 for all 1 j N . Let the numbers h j be such that
K j ρ ^ j max ln 4 ε 1 2 , ln 96 2 ρ ^ j 2 ε 1 2 , ln 12 2 λ ^ ρ ^ j ε 1 2 h j < j / 2 , 1 j N ,
where ρ ^ j = K j p 2 / ϵ j 2 λ ^ . Then, the following estimate holds:
min k 0 | λ k λ ^ | ε .

5. Proof of Theorem 2

Let u H 1 ( G ) be the eigenfunction of problem (2) corresponding to the eigenvalue λ such that u L 2 ( G ) = 1 . Recall that according to the hypothesis of the theorem, λ ˜ j = λ / K j < μ 1 , j for all 1 j N and ν k , j = μ k , j λ ˜ j for k 1 .
Henceforth, we will use the following notations:
T j , h j / 2 , h j = Ω j × [ ( h j / 2 , h j ) ( j h j , j h j / 2 ) ] , T j , h j / 4 , h j = Ω j × [ ( h j / 4 , h j ) ( j h j , j h j / 4 ) ] .
We introduce the function u ^ H h 1 ( G ) by the formula
u ^ ( x ) = u ( x ) , x G h / 2 = G N j = 1 T j , h j / 2 ¯ , u j ( y , z ) ζ j ( z ) r j ( y , z ) , x = ( y , z ) T j , h j / 2 , 1 j N ,
where r j ( y , z ) = u j ( y , z ) [ u j ] j ( z ) , ( y , z ) T j , h j / 2 , 1 j N and
ζ j ( z ) = 1 2 + 3 h j ( z 3 h j / 4 ) 16 h j 3 ( z 3 h j / 4 ) 3 , z [ h j / 2 , h j ] , 1 2 3 h j ( z j + 3 h j / 4 ) + 16 h j 3 ( z j + 3 h j / 4 ) 3 , z [ j h j , j h j / 2 ] , 1 , z [ h j , j h j ] .
Note that ζ j H 2 ( h j / 2 , j h j / 2 ) . Moreover,
ζ j ( h j / 2 ) = 0 , ζ j ( h j / 2 ) = 0 , ζ j ( j h j / 2 ) = 0 , ζ j ( j h j / 2 ) = 0 , 0 ζ j 1 , | ζ j | 3 h j , | ζ j | 24 h j 2 .
Let us expand r j ( z ) and D z r j ( z ) for almost all z ( 0 , j ) into a Fourier series converging in L 2 ( Ω j )
r j ( z ) = k = 1 c k , j ( z ) ψ k , j , D z r j ( z , y ) = k = 1 c k , j ( z ) ψ k , j ,
where c k , j ( z ) = ( r j ( z ) , ψ k , j ) L 2 ( Ω j ) , c k , j ( z ) = ( D z r j ( z ) , ψ k , j ) L 2 ( Ω j ) . It is clear that c 0 , j ( z ) = 0 , c 0 , j ( z ) = 0 .
It follows from (2) and (20) that
0 j ( D z r j ( z ) , η ( z ) ψ k , j ) L 2 ( Ω j ) d z + 0 j ( y r j ( z ) , η ( z ) y ψ k , j ) L 2 ( Ω j ) d z = λ ˜ j 0 j ( r j ( z ) , η ( z ) ψ k , j ) L 2 ( Ω j ) d z η C 0 ( 0 , j ) .
Taking into account that
( y r j ( z ) , y ψ k , j ) L 2 ( Ω j ) = μ k , j ( r j ( z ) , ψ k , j ) L 2 ( Ω j ) ,
we have
0 j c k , j ( z ) η ( z ) d z + 0 j ν k , j 2 c k , j ( z ) η ( z ) d z = 0 η C 0 ( 0 , j ) .
Thus, for all k 1 ,
c k , j ( z ) = ν k , j 2 c k , j ( z ) , z ( 0 , j ) .
Lemma 13.
Let h j < j / 2 , 1 j N . Then, for all k 1 , the following estimates hold
c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 ,
c k , j L 2 ( h j / 2 , j h j / 2 ) 2 2 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 ,
c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 , .
Proof. 
Since the coefficients c k , j satisfy Equation (30), then for all 0 < ξ < h j / 2 < z < h j the following formulas hold:
c k , j ( z ) = [ e ν k , j ( z ξ ) e ν k , j ( 2 j 3 ξ z ) ] c k , j ( ξ ) + [ e ν k , j ( j ξ z ) e ν k , j ( j 3 ξ + z ) ] c k , j ( j ξ ) 1 e 2 ν k , j ( j 2 ξ ) = e ν k , j ( z ξ ) [ 1 e 2 ν k , j ( j ξ z ) ] c k , j ( ξ ) + e ν k , j ( j ξ z ) [ 1 e 2 ν k , j ( z ξ ) ] c k , j ( j ξ ) 1 e 2 ν k , j ( j 2 ξ ) , c k , j ( j z ) = [ e ν k , j ( j z ξ ) e ν k , j ( j 3 ξ + z ) ] c k , j ( ξ ) + [ e ν k , j ( z ξ ) e ν k , j ( 2 j 3 ξ z ) ] c k , j ( j ξ ) 1 e 2 ν k , j ( j 2 ξ ) = e ν k , j ( j z ξ ) [ 1 e 2 ν k , j ( z ξ ) ] c k , j ( ξ ) + e ν k , j ( z ξ ) [ 1 e 2 ν k , j ( j z ξ ) ] c k , j ( j ξ ) 1 e 2 ν k , j ( j 2 ξ ) ,
From here,
c k , j 2 ( z ) 2 e 2 ν k , j ( z ξ ) c k , j 2 ( ξ ) + e 2 ν k , j ( j ξ z ) c k 2 ( j ξ ) , c k , j 2 ( j z ) 2 e 2 ν k , j ( j ξ z ) c k , j 2 ( ξ ) + e 2 ν k , j ( z ξ ) c k , j 3 ( j ξ ) .
Summing up these estimates, we have
c k , j 2 ( z ) + c k , j 2 ( j z ) 2 e 2 ν k , j ( z ξ ) + e 2 ν k , j ( j ξ z ) ] [ c k , j 2 ( ξ ) + c k , j 2 ( j ξ ) 2 e 2 ν 1 , j ( z ξ ) + e 2 ν 1 , j ( j ξ z ) ] [ c k , j 2 ( ξ ) + c k , j 2 ( j ξ ) .
From here,
e 2 ν 1 , j ξ c k , j 2 ( z ) + c k , j 2 ( z ) 2 e 2 ν 1 , j z + e 2 ν 1 , j ( j z ) c k , j 2 ( ξ ) + c k , j 2 ( ξ ) .
Integrating this inequality over ξ ( 0 , h j / 2 ) and over z ( h j / 2 , h j ) , we arrive at the estimate
c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 2 e ν 1 , j h j + e ν 1 , j ( 2 j 3 h j ) c k , j L 2 ( 0 , h j / 2 ) 2 + c k , j L 2 ( j h j / 2 , j ) 2 4 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 .
Integrating inequality (34) over ξ ( 0 , h j / 2 ) and over z ( h j / 2 , j / 2 ) , we have
( 1 e ν 1 , j h j ) c k , j L 2 ( h j / 2 , j h j / 2 ) 2 2 ( e ν 1 , j h j e 2 ν 1 , j ( j h j / 2 ) ) c k , j L 2 ( 0 , h j / 2 ) 2 + c k , j L 2 ( j h j / 2 , j ) 2 2 e ν 1 , j h j c k , j L 2 ( 0 , h j / 2 ) 2 + c k , j L 2 ( j h j / 2 , j ) 2 .
From here,
c k , j L 2 ( h j / 2 , j h j / 2 ) 2 2 e ν 1 , j h j c k , j L 2 ( 0 , h j / 2 ) 2 + c k , j L 2 ( j h j / 2 , j ) 2 + e ν 1 , j h j c k , j L 2 ( h j / 2 , j h j / 2 ) 2 2 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 .
It follows from (30) that the coefficients c k , j satisfy the equation ( c k , j ) = ν k , j 2 c k , j . Therefore, estimate (33) is proved in exactly the same way as estimate (31).
The lemma is proved. □
Corollary 4.
Let h j < j / 2 , 1 j N . Then, the following estimates hold:
r j L 2 ( T j , h j / 2 , h j ) 2 4 e ν 1 , j h j u j L 2 ( T j ) 2 ,
D z r j L 2 ( T j , h j / 2 , h j ) 2 4 e ν 1 , j h j D z u j L 2 ( T j ) 2 .
Proof. 
Indeed, it follows from (31) that
r j L 2 ( T j , h j / 2 , h j ) 2 = k = 1 c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν 1 , j h j k = 1 c k , j L 2 ( 0 , j ) 2 = 4 e ν 1 , j h j u j [ u j ] j L 2 ( T j ) 2 4 e ν 1 , j h j u j L 2 ( T j ) 2 .
Similarly, estimate (37) follows from (33).
The corollary is proved. □
Lemma 14.
Let 1 ν 1 , j ln 2 ε 1 2 h j j / 2 , 1 j N . Then, the following estimate holds:
u ^ u L 2 ( G ) ε 1 .
Proof. 
Using estimate (32) and noting that 2 e ν 1 , j h j ε 1 2 , we have
u ^ u L 2 ( G ) 2 j = 1 N r j L 2 ( T j , h j / 2 ) 2 = j = 1 N k = 1 c k , j L 2 ( h j / 2 , j h j / 2 ) 2 j = 1 N k = 1 2 e ν 1 , j h j c k , j L 2 ( 0 , j ) 2 j = 1 N u j [ u j ] j L 2 ( T j ) 2 ε 1 2 j = 1 N u j L 2 ( T j ) 2 ε 1 2 u L 2 ( G ) 2 ε 1 2 = ε 1 2 .
The lemma is proved. □
Lemma 15.
The function u ^ satisfies the identity
( K u ^ , φ ) L 2 ( G ) λ ( u ^ , φ ) L 2 ( G ) = ( Ψ , φ ) L 2 ( G ) φ H h 1 ( G ) ,
where
Ψ ( x ) = 0 , x G N j = 1 T j , h j / 2 , h j , Ψ j ( x ) = K j [ ζ j ( z ) r j ( y , z ) + 2 ζ j ( z ) D z r j ( y , z ) ] , x = ( y , z ) T j , h j / 2 , h j , 1 j N .
Proof. 
Let φ H h 1 ( G ) , and let functions η j C [ 0 , j ] , where 1 j N , be such that η j ( z ) = 0 for z [ 0 , h j / 4 ] [ j h j / 4 , j ] , and η j ( z ) = 1 for z [ h j / 2 , j h j / 2 ] .
Introduce functions
ψ 0 ( x ) = φ ( x ) , x G N j = 1 T j , h j / 4 , ( 1 η j ( z ) ) φ ( y , z ) , x = ( y , z ) T j , h j / 4 , 1 j N . ψ j ( x ) = 0 , x G T j , h j / 4 , η j ( z ) φ ( y , z ) , x = ( y , z ) T j , h j / 4 .
As ϕ = j = 0 N ψ j then
( K u ^ , φ ) L 2 ( G ) λ ( u ^ , φ ) L 2 ( G ) = ( K u ^ , ψ 0 ) L 2 ( G ) λ ( u ^ , ψ 0 ) L 2 ( G ) + j = 1 N K j ( u ^ , ψ j ) L 2 ( T j , h j / 4 ) λ ˜ j ( u ^ , ψ j ) L 2 ( T j , h j / 4 )
Note that ψ 0 H 1 ( G ) , ψ 0 = 0 in G G h / 2 and u ^ = u in G h / 2 . So,
( K u ^ , ψ 0 ) L 2 ( G ) λ ( u ^ , ψ 0 ) L 2 ( G ) = ( K u , ψ 0 ) L 2 ( G ) λ ( u , ψ 0 ) L 2 ( G ) = 0 .
Note also that
( u ^ , ψ j ) L 2 ( T j , h j / 4 ) λ ˜ j ( u ^ , ψ j ) L 2 ( T j , h j / 4 ) = ( D z u ^ , D z ψ j ) L 2 ( T j , h j / 4 , h j ) + ( y u ^ , y ψ j ) L 2 ( T j , h j / 4 , h j ) λ ˜ j ( u ^ , ψ j ) L 2 ( T j , h j / 4 , h j ) + ( D z [ u j ] j , D z ψ j ) L 2 ( T j , h j ) λ ˜ j ( [ u j ] j , ψ j ) L 2 ( T j , h j ) .
Integrating by parts over z and taking into account that ψ j ( h j / 4 ) = 0 , ψ j ( j h j / 4 ) = 0 , and D z 2 [ u j ] j = λ ˜ j [ u j ] j in T j , h j , we have
( u ^ , ψ j ) L 2 ( T j , h j / 4 ) λ ˜ j ( u ^ , ψ j ) L 2 ( T j , h j / 4 ) = ( D z 2 u ^ + λ ˜ j u ^ , ψ j ) L 2 ( T j , h j / 4 , h j ) + ( y u ^ , y ψ j ) L 2 ( T j , h j / 4 , h j ) + ( D z u ^ , ψ j ) L 2 ( Ω j ) | z = h j 0 + ( D z [ u j ] j , ψ j ) L 2 ( Ω j ) | z = j h j 0 ( D z [ u j ] j , ψ j ) L 2 ( Ω j ) | z = h j + 0 ( D z u ^ , ψ j ) L 2 ( Ω j ) | z = j h j 0 .
Noting that D z u ^ ( h j 0 ) = D z [ u j ] j ( h j 0 ) does not depend on y and [ ψ j ] j ( h j 0 ) = ψ j ( h j + 0 ) , we have
( D z u ^ , ψ j ) L 2 ( Ω j ) | z = h j 0 ( D z [ u j ] j , ψ j ) L 2 ( Ω j ) | z = h j + 0 = ( D z [ u j ] j , [ ψ j ] j ) L 2 ( Ω j ) | z = h j 0 ( D z [ u j ] j , ψ j ) L 2 ( Ω j ) | z = h j + 0 = 0 .
Likewise
( D z u ^ , ψ j ) L 2 ( Ω j ) | z = l h j 0 ( D z [ u ] j , ψ j ) L 2 ( Ω j ) | j h j + 0 = 0 .
Thus, using (5) we obtain
( u ^ , ψ j ) L 2 ( T j , h j / 4 ) λ ˜ j ( u ^ , ψ j ) L 2 ( T j , h j / 4 ) = ( D z 2 u ^ + λ ˜ j u ^ , ψ j ) L 2 ( T j , h j / 4 , h j ) + ( y u ^ , y ψ j ) L 2 ( T j , h j / 4 , h j ) = ( D z 2 u + λ ˜ j u , ψ j ) L 2 ( T j , h j / 4 , h j ) + ( y u , y ψ j ) L 2 ( T j , h j / 4 , h j ) + ( D z 2 ( ζ j r j ) + λ ˜ j ( ζ j r j ) , ψ j ) L 2 ( T j , h j / 2 , h j ) ( y ( ζ j r j ) , y ψ j ) L 2 ( T j , h j / 2 , h j ) = ( ζ j r j + 2 ζ j D z r j , ψ j ) L 2 ( T j , h j / 2 , h j ) + + ( D z 2 r j + λ ˜ j r j , ζ j ψ j ) L 2 ( T j , h j / 2 , h j ) ( y r j , ζ j y ψ j ) L 2 ( T j , h j / 2 , h j ) .
By virtue of identity (6), the last two terms vanish. Noting that ψ j = φ on T j , h j / 2 , h j , we obtain
( K u ^ , φ ) L 2 ( G ) λ ( u ^ , φ ) L 2 ( G ) = j = 1 N K j ( ζ j r j + 2 ζ j D z r j , ψ j ) L 2 ( ) = ( Ψ , φ ) L 2 ( G ) .
The lemma is proved. □
Lemma 16.
Assume that the numbers h j satisfy the conditions in (24). Then, the following inequalities hold:
2 K j 2 48 2 h j 4 e ν 1 , j h j 1 2 ε 1 2 , 2 K j 12 2 h j 2 e ν 1 , j h j 1 2 λ ε 1 2 .
Proof. 
Let us set t = ν 1 , j h j and transform the inequalities in (40) to the form
t + 4 ln t α j , t + 2 ln t β j
where α j = ln 96 2 K j 2 ν 1 , j 4 ε 1 2 , β j = ln 24 2 K j λ ν 1 , j 2 ε 1 2 .
The function g ( t ) = t + 4 ln t is increasing. If α j 1 , then g ( t ) g ( 1 ) = 1 α j for t 1 . If α j 1 , then g ( t ) g ( α j ) α j for t α j . Thus, the first inequality in (41) holds for h j 1 ν 1 , j max α j , 1 .
Similarly, the second inequality (41) holds for h j 1 ν 1 , j max β j , 1 .
It should be noted that (24) implies that h j ν 1 , j max α j , β j , 1 . Furthermore, note that ν 1 , j h j > ln ( 2 / ε 1 2 ) already implies that ν 1 , j h j > 1 .
The lemma is proved. □
Lemma 17.
Assume that the numbers h j satisfy the conditions in (24). Then, the following estimate is true:
Ψ L 2 ( G ) ε 1 .
Proof. 
Using estimates (29), (36), (37), we have
Ψ L 2 ( G ) 2 = j = 1 N Ψ j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 ζ j r j L 2 ( T j , h j / 2 , h j ) 2 + 2 ζ j D z r j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 24 2 h j 4 r j L 2 ( T j , h j / 2 , h j ) 2 + 6 2 h j 2 D z r j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 e ν 1 , j h j 48 2 h j 4 u j L 2 ( T j ) 2 + 12 2 h j 2 u j L 2 ( T j ) 2 .
Due to the estimates in (40), we have
Ψ L 2 ( G ) 2 j = 1 N u j L 2 ( T j ) 2 + 1 λ K j u j L 2 ( T j ) 2 ε 1 2 2 u L 2 ( G ) + 1 λ K u L 2 ( G ) 2 ε 1 2 2 .
Taking into account that 1 λ K u L 2 ( G ) 2 = u L 2 ( G ) 2 = 1 , we arrive at inequality (42).
The lemma is proved. □
Let us now show that estimate (25) holds under the conditions of Theorem 2.
Let { λ ^ k } k = 0 be the set of eigenvalues to problem (17) and { φ ^ k } k = 0 be the set of corresponding eigenfunctions, forming an orthonormal base in L h 2 ( G ) .
From (39) and equality
( K u ^ , φ ^ k ) L 2 ( G ) = λ ^ k ( u ^ , φ ^ k ) L 2 ( G ) , k 0
it follows that
( λ ^ k λ ) ( u ^ , φ ^ k ) L 2 ( G ) = ( Ψ , φ ^ k ) L 2 ( G ) , k 0 .
So,
min k 0 ( λ ^ k λ ) 2 u ^ L 2 ( G ) 2 k = 0 ( λ ^ k λ ) 2 ( u ^ , φ ^ k ) L 2 ( G ) 2 k = 0 ( Ψ , φ ^ k ) L 2 ( G ) 2 Ψ L 2 ( G ) 2 .
Hence,
min k 0 | λ ^ k λ | ( u L 2 ( G ) u ^ u L 2 ( G ) ) min k 0 | λ ^ k λ | u ^ L 2 ( G ) Ψ L 2 ( G ) .
As u L 2 ( G ) = 1 , u ^ u L 2 ( G ) ε 1 and Ψ L 2 ( G ) ε 1 , we obtain the estimate
min k 0 | λ ^ k λ | ( 1 ε 1 ) ε 1 ,
which is equivalent to estimate (25).
Theorem 2 is proved.

6. Proof of Theorem 3

Let U H h 1 ( G ) be the eigenfunction of problem (17) corresponding to the eigenvalue λ ^ and such that U L 2 ( G ) = 1 . Recall that according to the hypothesis of the theorem, λ ¯ j = λ ^ / K j < μ 1 , j for all 1 j N and ν ^ k , j = μ k , j λ ¯ j for k 1 .
We introduce the function U ^ H 1 ( G ) by the formula
U ^ ( x ) = U ( x ) , x G N j = 1 T j , h j / 2 , h j U j ( y , z ) ζ j ( z ) R j ( y , z ) , x = ( y , z ) T j , h j / 2 , h j , 1 j N ,
where
R j ( y , z ) = U j ( y , z ) [ U j ] j ( z ) , ( y , z ) T j , h j / 2 , h j , 1 j N ,
and the function ζ j is defined by Formula (28).
Let us expand R j ( z ) and D z R j ( z ) for almost all z ( 0 , h j ) ( j h j , j ) into a Fourier series converging in L 2 ( Ω j ) :
R j ( z ) = k = 1 c k , j ( z ) ψ k , j , D z R j ( z ) = k = 1 c k , j ( z ) ψ k , j .
where c k , j ( z ) = ( R j ( z ) , ψ k , j ) L 2 ( Ω j ) , c k , j ( z ) = ( D z R j ( z ) , ψ k , j ) L 2 ( Ω j ) .
It follows from (17) that
0 h j ( D z R j ( z ) , η ( z ) ψ k , j ) L 2 ( Ω j ) d z + 0 h j ( y R j ( z ) , η ( z ) y ψ k , j ) L 2 ( Ω j ) d z = λ ¯ j 0 h j ( R j ( z ) , η ( z ) ψ k , j ) L 2 ( Ω j ) d z η C 0 ( 0 , h j ) .
Taking into account that
( y R j ( z ) , y ψ k , j ) L 2 ( Ω j ) = μ k , j ( R j ( z ) , ψ k , j ) L 2 ( Ω j ) ,
for almost all z ( 0 , h j ) ( j h j , j ) , we have
0 h j c k , j ( z ) η ( z ) d z + 0 h j ν ^ k , j 2 c k , j ( z ) η ( z ) d z = 0 η C 0 ( 0 , h j ) C 0 ( j h j , j ) .
Thus,
c k , j ( z ) = ν ^ k , j 2 c k , j ( z ) , z ( 0 , h j ) ( j h j , j ) .
Remember that D z R j ( h j ) = 0 , D z R j ( j h j ) = 0 (cf. Corollary 3.5). Hence,
c k , j ( h j ) = 0 , c k , j ( j h j ) = 0 k 1 .
Lemma 18.
Let h j j / 2 , 1 j N . Then, for all k 1 , the following estimates hold:
c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν ^ 1 , j h j c k , j L 2 ( 0 , j ) 2 ,
c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 e ν ^ 1 , j h j c k , j L 2 ( 0 , j ) 2 .
Proof. 
Let 0 < ξ < h j / 2 < z < h j . As the coefficients c k , j satisfy equation (43) and the condition c k , j ( h j ) = 0 , we have
c k , j ( z ) = c k , j ( ξ ) e ν ^ k , j ( z ξ ) 1 + e 2 ν ^ k , j ( h j z ) 1 + e 2 ν ^ k , j ( h j ξ ) .
From here,
c k , j 2 ( z ) 4 c k , j 2 ( ξ ) e 2 ν ^ k , j ( z ξ ) 4 c k , j 2 ( ξ ) e 2 ν ^ 1 , j ( z ξ ) .
So,
c k , j 2 ( z ) e 2 ν ^ 1 , j ξ 4 c k , j 2 ( ξ ) e 2 ν ^ 1 , j z .
Integrating this inequality over ξ ( 0 , h j / 2 ) and z ( h j / 2 , h j ) , we come to an inequality
c k , j L 2 ( h j / 2 , h j ) 2 4 e ν ^ 1 , j h j c k , j L 2 ( 0 , h j / 2 ) 2 .
The following inequality is proved in a similar way
c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν ^ 1 , j h j c k L 2 ( j h j / 2 , j ) 2 .
By adding inequalities (46) and (47) and coarsening the result, we arrive at inequality (44).
From (43), it follows that the coefficients c k , j satisfy the equation ( c k ) = ν ^ k , j 2 c k . Taking into account that c k , j ( h j ) = 0 , we have for all 0 < ξ < h j / 2 < z < h j
c k , j ( z ) = c k , j ( ξ ) e ν ^ k , j ( z ξ ) 1 e 2 ν ^ k , j ( h j z ) 1 e 2 ν ^ k , j ( h j ξ ) .
From here,
| c k ( z ) | 2 | c k ( ξ ) | 2 e 2 ν ^ k , j ( z ξ ) | c k ( ξ ) | 2 e 2 ν ^ 1 , j ( z ξ ) .
Hence,
e 2 ν ^ 1 , j ξ | c k , j ( z ) | 2 | c k , j ( ξ ) | 2 e 2 ν ^ 1 , j z .
Integrating this inequality over ξ ( 0 , h j / 2 ) and z ( h j / 2 , h j ) , we come to an inequality
c k , j L 2 ( h j / 2 , h j ) 2 e ν ^ 1 , j h j c k , j L 2 ( 0 , h j / 2 ) 2 .
In the same way, we prove the inequality
c k , j L 2 ( j h j , j h j / 2 ) 2 e ν ^ 1 , j h j c k , j L 2 ( j h j / 2 , j ) 2 .
By adding inequalities (48) and (49) and coarsening the result, we arrive at (45).
The lemma is proved. □
Corollary 5.
Let h j < j / 2 , 1 j N . Then, the following estimates hold:
R j L 2 ( T j , h j / 2 , h j ) 2 4 e ν ^ 1 , j h j U j L 2 ( T j ) 2 ,
D z R j L 2 ( T j , h j / 2 , h j ) 2 e ν ^ 1 , j h j D z U j L 2 ( T j ) 2 .
Proof. 
Indeed,
R j L 2 ( T j , h j / 2 , h j ) 2 = k = 1 c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 4 e ν ^ 1 , j h j k = 1 c k , j L 2 ( 0 , j ) 2 = 4 e ν ^ 1 , j h j U j [ U j ] j L 2 ( T j ) 2 4 e ν ^ 1 , j h j U j L 2 ( T j ) 2 , D z R j L 2 ( T j , h j / 2 , h j ) 2 = k = 1 c k , j L 2 ( h j / 2 , h j ) 2 + c k , j L 2 ( j h j , j h j / 2 ) 2 e ν ^ 1 , j h j k = 1 c k , j L 2 ( 0 , j ) 2 = e ν ^ 1 , j h j D z U j [ D z U j ] j L 2 ( T j ) 2 e ν ^ 1 , j h j D z U j L 2 ( T j ) 2 .
The lemma is proved. □
Lemma 19.
Let 1 ν ^ 1 , j ln 4 ε 1 2 h j < j / 2 for all 1 j N . Then, the following estimate holds:
U ^ U L 2 ( G ) ε 1 .
Proof. 
Using estimate (50) and taking into account that 4 e ν ^ 1 , j h j ε 1 2 , we have
U ^ U L 2 ( G ) 2 = j = 1 N ζ j R j L 2 ( T j , h j / 2 , h j ) 2 j = 1 N R j L 2 ( T j , h j / 2 , h j ) 2 j = 1 N 4 e ν ^ 1 , j h j U j L 2 ( T j ) 2 j = 1 N U j L 2 ( T j ) 2 ε 1 2 U L 2 ( G ) 2 ε 1 2 = ε 1 2 .
The lemma is proved. □
Lemma 20.
The following identity is true:
( K U ^ , φ ) L 2 ( G ) λ ^ ( U ^ , φ ) L 2 ( G ) = ( Φ , φ ) L 2 ( G ) φ H 1 ( G ) ,
where
Φ ( x ) = 0 , x G N j = 1 T j , h j / 2 , h j , Φ j ( x ) = K j [ ζ j ( z ) R j ( y , z ) + 2 ζ j ( z ) D z R j ( y , z ) ] , x = ( y , z ) T j , h j / 2 , h j , 1 j N .
Proof. 
From (18), it follows that
( K U ^ , φ ) L 2 ( G ) λ ^ ( U ^ , φ ) L 2 ( G ) = ( K ( U ^ U ) , φ ) L 2 ( G ) λ ^ ( U ^ U , φ ) L 2 ( G ) = j = 1 N K j ( D z ( ζ j R j ) , D z φ ) L 2 ( T j , h j / 2 , h j ) + ( ζ j y R j , y φ ) L 2 ( T j , h j / 2 , h j ) λ ¯ j ( ζ j R j , φ ) L 2 ( T j , h j / 2 , h j ) .
Using the fact that D z ( ζ j R j ) ( z ) = 0 in L 2 ( Ω j ) for z = h j / 2 , h j , j h j , j h j / 2 , we perform the transformation
( D z ( ζ j R j ) , D z φ ) L 2 ( T j , h j / 2 , h j ) = ( D z 2 ( ζ j R j ) , φ ) L 2 ( T j , h j / 2 , h j ) = ( ζ j D z 2 R j , φ ) L 2 ( T j , h j / 2 , h j ) ( ζ j R j + 2 ζ j D z R j , φ ) L 2 ( T j , h j / 2 , h j ) .
As a result, we obtain
( K U ^ , φ ) L 2 ( G ) λ ^ ( U ^ , φ ) L 2 ( G ) = ( Φ , φ ) L 2 ( G ) + j = 1 N K j [ ( D z 2 R j + λ ¯ j R j , ζ j φ ) L 2 ( T j , h j / 2 , h j ) ( y R j , ζ j y φ ) L 2 ( T j , h j / 2 , h j ) ] .
Using identity (23), we arrive at (53).
The lemma is proved. □
Lemma 21.
Assume that the numbers h j satisfy the conditions in (26). Then, the following inequalities hold:
2 K j 2 48 2 h j 4 e ν ^ 1 , j h j ε 1 2 2 , 2 K j 6 2 h j 2 e ν ^ 1 , j h j ε 1 2 2 λ ^ .
The proof of this lemma repeats the proof of lemma 16.
Lemma 22.
Assume that the numbers h j satisfy the conditions in (26). Then, the following estimate holds:
Φ L 2 ( G ) ε 1 .
Proof. 
Using estimates (29), (50), and (51), we have
Φ L 2 ( G ) 2 = j = 1 N Φ j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 ζ R j L 2 ( T j , h j / 2 , h j ) 2 + 2 ζ D z R j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 24 2 h j 4 R j L 2 ( T j , h j / 2 , h j ) + 6 2 h j 2 D z R j L 2 ( T j , h j / 2 , h j ) 2 2 j = 1 N K j 2 e ν ^ 1 , j h j 48 2 h j 4 U j L 2 ( T j ) 2 + 6 2 h j 2 D z U j L 2 ( T j ) 2 .
Using estimates (54), we have
Φ L 2 ( G ) 2 j = 1 N U j L 2 ( T j ) 2 + 1 λ ^ K j D z U j L 2 ( T j ) 2 ε 1 2 2 U L 2 ( G ) 2 + 1 λ ^ K U L 2 ( G ) 2 ε 1 2 2 .
Taking into account that 1 λ ^ K U L 2 ( G ) 2 = U L 2 ( G ) 2 = 1 , we arrive at inequality (55).
The lemma is proved. □
Let us now establish the validity of Theorem 3.
Remember that U ^ H 1 ( G ) . So
( K U ^ , φ k ) L 2 ( G ) = λ k ( U ^ , φ k ) L 2 ( G ) , k 0
It follows from (53) and (56) that
( λ k λ ^ ) ( U ^ , φ k ) L 2 ( G ) = ( Φ , φ k ) L 2 ( G ) , k 0 .
So,
min k 0 ( λ k λ ^ ) 2 U ^ L 2 ( G ) 2 k = 0 ( λ k λ ^ ) 2 ( U ^ , φ k ) L 2 ( G ) 2 k = 0 ( Φ , φ k ) L 2 ( G ) 2 = Φ L 2 ( G ) 2 .
Hence,
min k 0 | λ k λ ^ | ( U L 2 ( G ) U ^ U L 2 ( G ) ) min k 0 | λ k λ ^ | U ^ L 2 ( G ) Φ L 2 ( G ) .
As U L 2 ( G ) = 1 , U ^ U L 2 ( G ) ε 1 and Φ L 2 ( G ) ε 1 , we obtain the estimate
min k 0 | λ k λ ^ | ( 1 ε 1 ) ε 1 ,
which is equivalent to estimate (27).
Theorem 3 is proved.

7. Conclusions

A new method reducing computational resources is introduced to find a set of first eigenvalues of the Neumann diffusion operator in a three-dimensional domain containing thin tubes. The method consists of the truncation of the tubes at some small distance from the ends of the tubes. The truncated parts are replaced by one-dimensional segments, and special junction conditions are stated on the interfaces of the three-dimensional parts and one-dimensional segments: pointwise continuity of the fluxes and continuity of an average of the eigenfunction. The method is justified by the theorems pointing out at what distance we can truncate the tubes while keeping the given accuracy for the approximations of eigenvalues. In the case in which the domain consists of thin tubes only, the method significantly reduces the computational time: if the ratio of the thicknesses of tubes to their lengths is 1 / m , the time is reduced m times. This acceleration brings a significant gain in time and allows the computations of the eigenvalues in domains of complex geometry. Using the technique developed here, the approach to the eigenvalues of the original and reduced problems preserving their multiplicities (up to a prescribed accuracy), and the approaches of the corresponding eigenfunctions will be addressed in a forthcoming paper by the authors. This extension shell likely involve the construction of sets of almost orthogonal eigenfunctions of each problem (cf. [26,27] for an abstract framework and [28] for the technique in a singularly perturbed spectral problem). The analysis of the method will be extended to the comparison of the eigenfunction of the original and partially decomposed problems.

Author Contributions

This article was written in equal co-authorship. All authors have read and agreed to the published version of the manuscript.

Funding

The study by the first author was supported by a grant from the Russian Science Foundation (project no. 19-11-00033); the second and fourth authors were supported by the grant Gob. Cantabria-UC, Ref. 20.VP66.64662; and the third author was supported by the European Social Fund (project No 09.3.3-LMT-K-712-17-003) under a grant agreement with the Research Council of Lithuania (LMTLT).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Panasenko, G. Method of asymptotic partial decomposition of domain. Math. Model. Methods Appl. Sci. 1998, 8, 139–156. [Google Scholar] [CrossRef]
  2. Panasenko, G.P. Asymptotic analysis of bar systems. I. Russ. J. of Math. Phys. 1994, 2, 325–352. [Google Scholar]
  3. Panasenko, G. Multi-Scale Modelling for Structures and Composites; Springer: Dordrecht, The Netherlands, 2005; 398p. [Google Scholar]
  4. Panasenko, G. Method of asymptotic partial domain decomposition for non-steady problems: Heat equation on a thin structure. Math. Commun. 2014, 19, 453–468. [Google Scholar]
  5. Panasenko, G.; Pérez, M.-E. Method of asymptotic partial decomposition of domain for spectral problems in rod structures. J. Math. Pures Appl. 2007, 87, 1–36. [Google Scholar] [CrossRef]
  6. Borisov, D.; Freitas, P. Singular asymptotic expansions for Dirichlet eigenvalues and eigenfunctions on thin planar domains. Ann. I.H.Poincaré - AN 2009, 26, 547–560. [Google Scholar] [CrossRef]
  7. Friedlander, L.; Solomyak, M. On the spectrum of the Dirichlet Laplacian in a narrow strip. Israel J. Math. 2009, 170, 337–354. [Google Scholar] [CrossRef]
  8. Gaudiello, A.; Gómez, D.; Pérez-Martínez, M.-E. Asymptotic analysis of the high frequencies for the Laplace operator in a thin T-like shaped structure. J. Math. Pures Appl. 2020, 134, 299–327. [Google Scholar] [CrossRef]
  9. Gaudiello, A.; Gómez, D.; Pérez-Martínez, M.-E. A spectral problem for the Laplacian in joined thin films. Calc. Var. Partial Differ. Equ. 2023, 62, 129. [Google Scholar] [CrossRef]
  10. Kuchment, P.; Zheng, H. Asymptotics of spectra of Neumann Laplacians in thin domains. In Contemporary Mathematics; AMS: Providence, RI, USA, 2003; Volume 327, pp. 199–213. [Google Scholar]
  11. Molchanov, S.; Vainberg, B. Laplace operator in networks of thin fibers: Small diameter asymptotics. Comm. Math. Phys. 2007, 273, 533–559. [Google Scholar] [CrossRef]
  12. Nazarov, S.A.; Pérez, E.; Taskinen, J. Localization effect for Dirichlet eigenfunctions in thin non-smooth domains. Trans. Amer. Math. Soc. 2016, 368, 4787–4829. [Google Scholar] [CrossRef]
  13. Rubinstein, J.; Schatzmann, M. Variational problems on multiply connected thin strips I: Basic estimates and convergence of the Laplacian spectrum. Arch. Rat. Mech. Anal. 2001, 160, 271–308. [Google Scholar] [CrossRef]
  14. Arrieta, J.M. Rates of eigenvalues on a dumbbell domain. Simple eigenvalue case. Trans. Am. Math. Soc. 1995, 347, 3503–3531. [Google Scholar] [CrossRef]
  15. Gadyl’shin, R.R. On the eigenvalues of a “dumbbell with a thin handle”. Izv. Ross. Akad. Nauk Ser. Mat. 2005, 69, 45–110, Translation in Izv. Math. 2005, 69, 265–329. [Google Scholar]
  16. Melnik, T.A.; Nazarov, S.A. Asymptotic behavior of the solution of the Neumann spectral problem in a domain of “tooth comb” type. Tr. Semin. Im. I. G. Petrovsk. 1996, 19, 138–173, Translation in J. Math. Sci. 1997, 85, 2326–2346. [Google Scholar]
  17. Melnik, T.A.; Nazarov, S.A. Asymptotic analysis of the Neumann problem on the junction of a body and thin heavy rods. Algebra i Anal. 2000, 12, 188–238, Translation in St. Petersburg Math. J. 2001, 12, 317–351. [Google Scholar]
  18. Panasenko, G.; Viallon, M.-C. Method of asymptotic partial decomposition with discontinuous junctions. Comput. Math. Appl. 2022, 105, 75–93. [Google Scholar] [CrossRef]
  19. Brenner, S.C.; Scott, L.R. The Mathematical Theory of Finite Element Methods; Springer: New York, NY, USA, 2008. [Google Scholar]
  20. Zhan, Q.; Sun, Q.; Ren, Q.; Fang, Y.; Wang, H.; Liu, Q.H. A discontinuous Galerkin method for simulating the effects of arbitrary discrete fractures on elastic wave propagation. Geophys. J. Int. 2017, 210, 1219–1230. [Google Scholar] [CrossRef]
  21. Amosov, A.A.; Panasenko, G.P. Partial dimension reduction for the heat equation in a domain containing thin tubes. Math. Methods Appl. Sci. 2018, 41, 9529–9545. [Google Scholar] [CrossRef]
  22. Szego, G. Inequalities for certain eigenvalues of a membrane of given area. J. Ration. Mech. Anal. 1954, 3, 343–356. [Google Scholar] [CrossRef]
  23. Schoen, R.; Yau, S.-T. Lectures on Differential Geometry; International Press: Cambridge, MA, USA, 1994. [Google Scholar]
  24. Brenner, S.C. Poincaré-Friedrichs inequalities for piecewise H1 functions. SIAM J. Numer. Anal. 2003, 41, 304–324. [Google Scholar] [CrossRef]
  25. Sanchez-Hubert, J.; Sanchez-Palencia, E. Vibration and Coupling of Continuous Systems; Asymptotic Methods; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  26. Lazutkin, V.F. KAM Theory and Semiclassical Approximations to Eigenfunctions, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  27. Visik, M.I.; Lusternik, L.A. Regular degeneration and boundary layer for linear differential equations with small parameter. Trans. Am. Math. Soc. Ser. 2 1957, 20, 239–364. [Google Scholar]
  28. Lobo, M.; Pérez, E. On the local vibrations for systems with many concentrated masses near the boundary. Comptes Rendus Acad. Sci. Paris 1997, 324, 323–329. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Amosov, A.; Gómez, D.; Panasenko, G.; Pérez-Martinez, M.-E. Asymptotic Domain Decomposition Method for Approximation the Spectrum of the Diffusion Operator in a Domain Containing Thin Tubes. Mathematics 2023, 11, 3592. https://doi.org/10.3390/math11163592

AMA Style

Amosov A, Gómez D, Panasenko G, Pérez-Martinez M-E. Asymptotic Domain Decomposition Method for Approximation the Spectrum of the Diffusion Operator in a Domain Containing Thin Tubes. Mathematics. 2023; 11(16):3592. https://doi.org/10.3390/math11163592

Chicago/Turabian Style

Amosov, Andrey, Delfina Gómez, Grigory Panasenko, and Maria-Eugenia Pérez-Martinez. 2023. "Asymptotic Domain Decomposition Method for Approximation the Spectrum of the Diffusion Operator in a Domain Containing Thin Tubes" Mathematics 11, no. 16: 3592. https://doi.org/10.3390/math11163592

APA Style

Amosov, A., Gómez, D., Panasenko, G., & Pérez-Martinez, M. -E. (2023). Asymptotic Domain Decomposition Method for Approximation the Spectrum of the Diffusion Operator in a Domain Containing Thin Tubes. Mathematics, 11(16), 3592. https://doi.org/10.3390/math11163592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop