![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
A class of nonlinear variational problems describing incompressible fluids and solids by stationary Stokes equations given in a planar domain with a crack (infinitely thin flat plate in fluids) is considered. Based on the Fourier asymptotic analysis, general analytical solutions are obtained in polar coordinates as the power series with respect to the distance to the crack tip. The logarithm terms and angular functions are accounted in the asymptotic expansion using recurrence relations. Then boundary conditions imposed between the opposite crack faces in the sector of angle determine admissible exponents and parameters in the power series. For the specific conditions of Dirichlet, Neumann, impermeability, non-penetration and shear crack, the principal asymptotic terms are derived, which verify the singular behaviour. In particular, the analytical solution answers the questions of a square-root singularity at the crack tip and the presence of log-oscillations of variational solutions for the Stokes problems.
1. Introduction
The Stokes system can describe stationary flow of incompressible fluids [Citation1] as well as incompressible linearly elastic solids that are subject to the divergence-free condition [Citation2]. We refer to suitable qualitative properties of solutions given in Ref. [Citation3] and related regularity results in Ref. [Citation4,Citation5] for dual strategies solving the Stokes problem. Motivated by applications to fracture mechanics, in our study we are interested in the presence of a crack in a reference domain (see the variational theory of crack problems in the monograph by Khludnev and Kovtunenko [Citation6]). For fluids, the discontinuity may be suitable when simplifying aeroplane wings to infinitely thin flat plates. This research has further important extensions to nonlinear continuum mechanics [Citation7,Citation8], shape optimization [Citation9–11], overdetermined [Citation12] and inverse problems [Citation13].
The main question of crack problems consists in finding a singularity at the crack tip. The generic singularity is of power type with respect to the distance to the crack tip and may contain logarithmic terms, whose presence or absence is of especial importance for engineers. If we consider linear elastic boundary value problems for cracked configurations, then we have the Williams series solution [Citation14]. The boundary value problems analyzed in the paper have long history. Highlighting questions related to the singularity in the vicinity of the crack we cite [Citation15], the respective approaches use regular and singular perturbations [Citation16–18], integral Fourier and Mellin transforms [Citation19]. Here we focus on asymptotic series solution for incompressible continua given in a planar domain corresponding to the sector of angle .
The crack singularities for general elliptic systems in planar domains were analyzed in Ref. [Citation20] in dependence of boundary conditions. Thus, the Lamé system describing isotropic bodies with a stress-free crack (the homogeneous Neumann boundary condition) obeys the classic square-root singularity without logarithms. It is worth noting that Stokes equations correspond to the limit case of the Lamé equations when the Poisson ratio (hence the Lamé parameter
). Setting
in asymptotic formulas was suggested in Ref. [Citation21] to get the displacement appropriate for the special case of an incompressible body. However, in the limit state of Stokes equations under non-penetration conditions, log-oscillations may occur for high-order asymptotic terms. In Ref. [Citation22, Chapters 5 and 6], Stokes problems in conical domains were investigated with respect to its eigenvalues and generalized eigenvalues generated by Dirichlet and Neumann conditions.
Commonly, boundary conditions are incorporated into the operator pencil for differential equations in order to determine admissible exponents and parameters in the power series after solving this eigenvalue problem [Citation3,Citation4,Citation20,Citation22]. This approach is highly dependent on the specific choice of conditions imposed on the boundary. In contrast, we obtain first a general asymptotic solution for the Stokes equations using Fourier asymptotic expansion and recurrence relations for logarithmic terms and angular functions. Then we can specify coefficients in the power series for every kind of boundary conditions imposed between the opposite crack faces, even for nonlinear ones. For engineers, it is important to note that we establish: if log-oscillations occur in the power series, or not. In this way, we provide the principal asymptotic terms for the condition of Dirichlet, Neumann, impermeability, non-penetration and shear crack presented next.
In the domain with crack , let unknown displacement (velocity for fluids) vector
and scalar pressure p satisfy the homogeneous Stokes equations:
(1)
(1) where the shear modulus
and Young's modulus E>0. This constitutes the Cauchy stress
in terms of the linearized strain
and pressure p as
(2)
(2) where
is the identity transformation. At the crack
with a normal vector
, we can distinguish traces of the discontinuous functions
,
, and p across the opposite crack faces
, thus the jump:
(3)
(3) and the mean:
(4)
(4) Using (Equation2
(2)
(2) )–(Equation4
(4)
(4) ), we set the following boundary conditions at
for the stick (homogeneous Dirichlet):
(5)
(5) impermeability (mixed homogeneous Dirichlet–Neumann):
(6)
(6) stress-free crack (homogeneous Neumann):
(7)
(7) non-penetration (complementarity):
(8)
(8) and shear crack (transmission with slip):
(9)
(9) The power series solution for the Stokes equation (Equation1
(1)
(1) ) when expressed in polar
-coordinates is found in the general form of sum
(10)
(10) where the angular functions are determined by the recurrence relations for
:
(11)
(11)
(12)
(12) with free coefficients
. The four eigenvectors
and
(where
) of the operator pencil are found from the condition that for every
functions
and
solve the homogeneous Stokes equation (Equation1
(1)
(1) ).
Inserting ansatz (Equation10(10)
(10) )–(Equation12
(12)
(12) ) into the boundary conditions (Equation5
(5)
(5) )–(Equation9
(9)
(9) ) as
and
, for every condition we determine the admissible eigenvalues γ, the presence of generalized eigenvectors by means of
, and coefficients
for
. The important observation is that, if only trivial
and
are possible in (Equation11
(11)
(11) ) and (Equation12
(12)
(12) ), then
and
too, thus
, and there is no logarithms. Otherwise, if non-trivial
exist, then natural number
may be arbitrary.
The energy -solution
excepting possible constant
restricts to
. In particular, we found the crack-tip singular of energy solutions (the principal asymptotic terms) according to the following expansions as
for the boundary conditions of Dirichlet (Equation5
(5)
(5) ), mixed Dirichlet–Neumann (Equation6
(6)
(6) ), Neumann (Equation7
(7)
(7) ) and non-penetration (Equation8
(8)
(8) ):
(13)
(13) using the Landau notations, whereas in the case (Equation9
(9)
(9) ) of shear crack:
(14)
(14) The structure of the paper is as follows. In Section 2, we reformulate in polar coordinates the boundary value problems for the Stokes model (Equation1
(1)
(1) ), (Equation2
(2)
(2) ) under the boundary conditions (Equation5
(5)
(5) )–(Equation9
(9)
(9) ) at the crack. Introducing proper feasible sets, we prove weak solutions to the corresponding variational Stokes problems in bounded domains with cracks. Based on the Fourier asymptotic expansion, in Section 3, we justify rigorously the power series solution (Equation10
(10)
(10) )–(Equation12
(12)
(12) ) and derive explicitly the eigenvectors
. In Section 4, the asymptotic expansion is specified in a sector of angle
around the crack tip for every boundary conditions (Equation5
(5)
(5) )–(Equation9
(9)
(9) ) by means of the eigenvalues γ, eigenvectors
, and generalized eigenvectors
. This allows us to write the principal asymptotic terms in (Equation13
(13)
(13) ) and (Equation14
(14)
(14) ) and validate singular properties of the solution. Some concluding remarks and further perspectives are outlined in Section 5.
2. Variational Stokes problems
We begin with the geometric description of a planar sectorial domain around a crack tip. Let be a star-shaped domain with respect to the origin
. Let it have a Lipschitz continuous boundary
, which consists of mutually disjoint parts
and
, and the normal vector
outward from Ω. We assume that the origin
coincides with the tip of a semi-infinite direct crack, whose intersection with Ω forms a line segment
. The cracked domain
implies a finite part of sector of angle
bounded by
. We introduce the polar coordinates r>0 and
such that the upper and lower crack faces correspond to
, respectively, as illustrated in Figure .
In the sectorial domain , we look for unknown vector
and scalar
functions. They constitute the linearized strain
and the Cauchy stress
defined according to (Equation2
(2)
(2) ), which in the polar coordinates are represented by symmetric tensors in
as
(15)
(15) where the strain components are
(16)
(16) using the convention for partial derivatives
and
.
Applying the constitutive relations (Equation15(15)
(15) ), the Stokes equation (Equation1
(1)
(1) ) turn into the equilibrium equations
(17)
(17) where the identity for the trace
was used, such that in polar coordinates
(18)
(18) At the outer boundary
, we endow (Equation17
(17)
(17) ) with the mixed, homogeneous Dirichlet–inhomogeneous Neumann boundary conditions
(19)
(19) for a prescribed force
, where
stands for boundary traction.
Next we set admissible function spaces for the Stokes problems (Equation17(17)
(17) ), (Equation19
(19)
(19) ) under boundary conditions (Equation5
(5)
(5) )–(Equation9
(9)
(9) ) at the crack
. Taking into account the no-slip condition in (Equation19
(19)
(19) ), we introduce the Sobolev space
(20)
(20) Admissible at crack functions are restricted in a feasible set
meeting the Dirichlet case (Equation5
(5)
(5) ) as
(21)
(21) the case of impermeability (Equation6
(6)
(6) ) as
(22)
(22) respectively in the Neumann case (Equation7
(7)
(7) )
(23)
(23) the case of non-penetration (Equation8
(8)
(8) ) implies
(24)
(24) and the shear crack from (Equation9
(9)
(9) ) reads
(25)
(25) In all cases (Equation21
(21)
(21) )–(Equation25
(25)
(25) ), the topological set K yields a convex closed cone.
For stress tensors such that
, the following Green's formula takes place using (Equation18
(18)
(18) ), and the notation from (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ):
(26)
(26) for smooth functions
such that
on
, and the strain operator
defined according to (Equation16
(16)
(16) ). After integration of Equation (Equation17
(17)
(17) ) with the help of constitutive relations (Equation15
(15)
(15) ), boundary conditions (Equation19
(19)
(19) ) and Green's formula (Equation26
(26)
(26) ), we arrive at the variational inequality: find
satisfying
(27)
(27) for all test functions
. This system is completed with the weak form of incompressibility as
(28)
(28) for all test functions
. Conversely, for the smooth solution
, the boundary conditions on the stress in (Equation5
(5)
(5) )–(Equation9
(9)
(9) ) follow from the variational inequality (Equation27
(27)
(27) ) after applying Green's formula (Equation26
(26)
(26) ).
Proposition 2.1
(Well-posedness of the weak Stokes problems). Let be a generic feasible set, which is convex and closed (e.g. (Equation21
(21)
(21) )–(Equation25
(25)
(25) )). Then there exists a pair
, the stress and strain
determined from (Equation15
(15)
(15) ) and (Equation16
(16)
(16) ), which solve uniquely the variational relations (Equation27
(27)
(27) ) and (Equation28
(28)
(28) ).
Proof.
In the proof we recall the well-known facts from the theory of mixed variational problems for the Stokes equations (see, e.g. Ref. [Citation5]). The weak formulation (Equation27(27)
(27) ) and (Equation28
(28)
(28) ) gives rise to the Lagrangian function
(the free energy), which is well-defined by
(29)
(29) The mapping
in (Equation29
(29)
(29) ) is quadratic, convex and coercive due to Poincare inequality provided by the Dirichlet condition in (Equation20
(20)
(20) ) such that
(30)
(30) The operator
in (Equation18
(18)
(18) ) is surjective. Therefore, the inf-sup condition (see Ref. [Citation1]) for the linear mapping
in (Equation29
(29)
(29) ) follows directly from the definition of dual norm in
:
(31)
(31) As the consequence of Ladyzhenskaya–Babuška–Brezzi–Nečas minimax theorem, there exists a unique saddle point
such that
(32)
(32) for all test functions
. After differentiation with respect to
and p, the optimality conditions for (Equation32
(32)
(32) ) lead to the variational inequality (Equation27
(27)
(27) ) and Equation (Equation28
(28)
(28) ). The proof is completed.
Evidently, the assertion of Proposition 2.1 holds for much more general geometries than assumed at the beginning of this section. In the subsequent sections, we use these geometric assumptions for construction of a semi-analytic solution for the variational Stokes problems locally in a neighbourhood of the crack tip.
3. Power series solution in the general form
To write the equilibrium equations in terms of , we substitute the stress
from (Equation15
(15)
(15) ) and strain
from (Equation16
(16)
(16) ) into the first two equations in (Equation18
(18)
(18) ), and use the incompressibility
such that
Due to the symmetry of mixed derivatives, differentiating (Equation16
(16)
(16) ) yields two compatibility conditions
Inserting these conditions justifies the equilibrium equations as
(33)
(33) In its turn, the incompressibility condition in terms of
reads
(34)
(34) We look for the solution of Stokes Equations (Equation33
(33)
(33) ) and (Equation34
(34)
(34) ) by the Fourier series (Equation10
(10)
(10) ) with respect to powers γ (see Ref. [Citation22]). For every fixed
, we consider the asymptotic terms
(35)
(35) with unknown angular vector
and scalar
functions. Here, the
-power for the pressure is due to compatibility of the terms in the stress
entering the constitutive equation (Equation15
(15)
(15) ). Our task needs few auxiliary lemmas.
Lemma 3.1
(Recurrence equations for angular functions). If the angular functions satisfy the following recurrence equations:
(36)
(36)
(37)
(37)
(38)
(38) then ansatz (Equation35
(35)
(35) ) fulfils the Stokes equations (Equation33
(33)
(33) ) and (Equation34
(34)
(34) ).
Proof.
For differentiation of the power series with respect to r, we exploit the following calculus:
(39)
(39) where shifting of the summation index was used. Inserting ansatz (Equation35
(35)
(35) ) into the incompressibility equation (Equation34
(34)
(34) ), after the differentiation of
with respect to θ and r as in (Equation39
(39)
(39) ), gathering like terms gives
(40)
(40) and necessitates relations (Equation36
(36)
(36) ).
Akin to (Equation40(40)
(40) ), we calculate the expression of
as
(41)
(41) and, similarly to (Equation39
(39)
(39) ), its derivative with respect to r as
(42)
(42) Inserting into equilibrium equation (Equation33
(33)
(33) ), the series for
from (Equation35
(35)
(35) ) and its derivative with respect to r as
the series (Equation41
(41)
(41) ) and (Equation42
(42)
(42) ) for
, this provides respectively
which leads to (Equation37
(37)
(37) ) and (Equation38
(38)
(38) ) and completes the proof.
With the help of differentiation with respect to θ, we decouple the relations in Lemma 3.1 as follows.
Lemma 3.2
(Decoupled equations for angular functions). The recurrence equations (Equation36(36)
(36) )–(Equation38
(38)
(38) ) are equivalent to the following subsequent system of equations for
,
, and
:
(43)
(43)
(44)
(44)
(45)
(45) excluding non-zero constant values of the left-hand side in Equation (Equation38
(38)
(38) ).
Proof.
Indeed, after differentiation of (Equation38(38)
(38) ) with respect to θ and using (Equation37
(37)
(37) ), it follows
for all
supposing
, hence (Equation43
(43)
(43) ). Inserting
and
from (Equation36
(36)
(36) ) into (Equation37
(37)
(37) ) and setting
yields (Equation44
(44)
(44) ). The relations (Equation37
(37)
(37) ) and (Equation45
(45)
(45) ) coincide, constant in the left-hand side of (Equation38
(38)
(38) ) due to the differentiation is excluded.
Since Equations (Equation43(43)
(43) ) and (Equation44
(44)
(44) ) imply an inhomogeneous Sturm–Liouville problem for two eigenvalues
and
, its solutions can be constructed as a linear span of four eigenvectors
(46)
(46) in the form of (Equation11
(11)
(11) ) and (Equation12
(12)
(12) ). For further use, we remind the properties of derivatives for (Equation46
(46)
(46) ):
(47)
(47)
(48)
(48) formulas of the Sturm–Liouville operator applied to these functions:
(49)
(49) and after differentiation with respect to γ:
(50)
(50)
Lemma 3.3
Recurrence formula of angular functions
The angular functions solving (Equation43(43)
(43) )–(Equation45
(45)
(45) ) for all
are determined by the recurrence sequence
(51)
(51)
(52)
(52)
(53)
(53) with arbitrary 10 coefficients
,
,
satisfying the six relations
(54)
(54)
(55)
(55)
The proof of Lemma 3.3 is leaded by induction over indexes and given in Appendix A.
Based on auxiliary Lemmas 3.1–3.3, we formulate the main result of this section.
Theorem 3.1
(Power series solution to Stokes equations). The power series solution for the Stokes equations in polar coordinates (Equation18(18)
(18) ) has the general form of the sum over powers of r:
(56)
(56) where the angular functions are determined as m = 0, and for
by the recurrence relations:
(57)
(57)
(58)
(58) If the eigenvalues
, the coefficients
are free, and
satisfy two relations
(59)
(59) the eigenvectors
are given according to
in (Equation46
(46)
(46) ) by
(60)
(60) If
, the coefficients
and
are free, such that all
in (Equation57
(57)
(57) ) depend only on the eigenvectors
and
from (Equation60
(60)
(60) ).
Proof.
Indeed, setting the ansatz in the form of power series (Equation56(56)
(56) ), formulas (Equation57
(57)
(57) )–(Equation60
(60)
(60) ) follow straightforwardly from the relations (Equation51
(51)
(51) )–(Equation55
(55)
(55) ), after solving Equation (Equation55
(55)
(55) ) with respect to
by division over
. If
, we get
, and coefficients
,
in (Equation55
(55)
(55) ) can be arbitrary. Since
and
as
, the corresponding eigenvector components
and
are trivial, whereas constant
and
are excluded in Lemma 3.2 since do not satisfy Equations (Equation37
(37)
(37) ) and (Equation38
(38)
(38) ). If
, then the pressure
due to (Equation59
(59)
(59) ) is skipped from the series for p in (Equation56
(56)
(56) ). The proof is completed.
It can be observed that if for fixed γ, then no logarithm terms occur in the expansion (Equation56
(56)
(56) ), and the system (Equation57
(57)
(57) )–(Equation59
(59)
(59) ) for angular functions consists of
and
only. Otherwise, if
, we have the log-oscillations in (Equation56
(56)
(56) ). In this respect, the important criterion of log-oscillations is given next.
Corollary 3.1
(Criterion of log-oscillations). For fixed in the system (Equation57
(57)
(57) ) and (Equation58
(58)
(58) ) as m = 0, 1:
(61)
(61) if only trivial functions
and
solve the inhomogeneous system from (Equation43
(43)
(43) ) to (Equation45
(45)
(45) ) as M = 1:
(62)
(62) then
and
too, thus continuing m>1 we conclude that
. Otherwise, if non-trivial solutions
to (Equation61
(61)
(61) ) and (Equation62
(62)
(62) ) exist, then
may be arbitrary natural number.
In the following section, we apply the general solution obtained in Theorem 3.1 to the variational Stokes problems under linear and nonlinear boundary conditions at the crack from (Equation21(21)
(21) ) to (Equation25
(25)
(25) ).
4. Specific boundary conditions and singularity at crack
We begin this section with deriving stresses from Theorem 3.1.
Corollary 4.1
(Power series for stresses). The stress components ,
in (Equation15
(15)
(15) ) are represented by the following series:
(63)
(63) The angular functions are determined as m = 0, and for
by the recurrence relations
(64)
(64)
(65)
(65) where
,
, and
are from (Equation57
(57)
(57) ) and (Equation58
(58)
(58) ).
Proof.
The representation (Equation63(63)
(63) )–(Equation65
(65)
(65) ) follows directly from Equations (Equation56
(56)
(56) )–(Equation58
(58)
(58) ) using the incompressibility condition in (Equation18
(18)
(18) ), after calculus of (Equation15
(15)
(15) ) for the components of stress, and (Equation16
(16)
(16) ) for strain. The calculation employs the differentiation of
with respect to r and θ akin to formulas (Equation39
(39)
(39) ) and (Equation40
(40)
(40) ).
Based on Corollaries 3.1 and 4.1, for fixed powers γ, we compute the first two terms in the general solution (Equation57(57)
(57) ), (Equation64
(64)
(64) ) and (Equation65
(65)
(65) ) as
explicitly:
(66)
(66) and
(67)
(67) for
. Here we have used
, and
and
from (Equation58
(58)
(58) ) due to (Equation59
(59)
(59) ) as
(68)
(68) If
, then we put
in (Equation66
(66)
(66) )–(Equation68
(68)
(68) ), such that
agrees (Equation38
(38)
(38) ).
At the crack faces corresponding to the angles
, the functions in (Equation46
(46)
(46) ) are
which jump and mean according to (Equation3
(3)
(3) ) and (Equation4
(4)
(4) ) are
(69)
(69) Therefore, inserting (Equation69
(69)
(69) ) into (Equation66
(66)
(66) )–(Equation68
(68)
(68) ) we find at the crack
the values
(70)
(70)
(71)
(71)
(72)
(72)
(73)
(73) and
(74)
(74)
(75)
(75)
(76)
(76)
(77)
(77) Based on (Equation70
(70)
(70) )–(Equation77
(77)
(77) ), we distinguish each case of the boundary conditions in (Equation21
(21)
(21) )–(Equation25
(25)
(25) ).
4.1. Crack under Dirichlet conditions of stick
According to the Dirichlet conditions in (Equation5(5)
(5) ), two pairs of equations
and
in (Equation70
(70)
(70) ) and (Equation71
(71)
(71) ) compose two homogeneous matrix equations with respect to free coefficients
:
(78)
(78) The solvability of (Equation78
(78)
(78) ) needs matrix determinant
equals to zero, i.e.
(79)
(79) which holds for the half-integer powers
(80)
(80) These eigenvalues distinguish three cases of possible solutions to (Equation78
(78)
(78) ):
(81)
(81) Here,
as
due to (Equation59
(59)
(59) ) and follows only trivial
in (Equation78
(78)
(78) ), and the corresponding pressure
is constant.
As m = 1, equations and
from (Equation74
(74)
(74) ) and (Equation75
(75)
(75) ) now build two inhomogeneous
-matrix equations with respect to coefficients
as follows:
(82)
(82) with the same system matrix as in (Equation78
(78)
(78) ). Then (Equation82
(82)
(82) ) is solvable when
(83)
(83) The conditions (Equation81
(81)
(81) ) and (Equation83
(83)
(83) ) together lead to trivial coefficients
for all γ. On the basis of Corollary 3.1, using (Equation80
(80)
(80) ) and (Equation81
(81)
(81) ) we conclude with the following result.
Theorem 4.1
(Power series for the Dirichlet crack problem). The Stokes equation (Equation18(18)
(18) ) with the crack under Dirichlet boundary conditions (Equation5
(5)
(5) ) have the power series solution without logarithms:
(84)
(84) since
and
, where the angular functions for every n are given by
(85)
(85) with the eigenvectors
and
from (Equation60
(60)
(60) ) and (Equation46
(46)
(46) ) as
. If
, six coefficients
satisfy the four relations in dependence of
:
(86)
(86) If n = 2, then
, and
is constant.
As a consequence of Theorem 4.1, for the variational solution from Proposition 2.1 requiring that n>0 (up to constant when
), we derive the principal asymptotic term in (Equation84
(84)
(84) ) as n = 1.
Corollary 4.2
(Singular solution for the Dirichlet crack problem). The variational solution to the Dirichlet problem has the following singularity at the crack tip as :
(87)
(87) the pressure
and stresses
Proof.
As n = k = 1, i.e. , we calculate
, then
,
in (Equation86
(86)
(86) ), and from (Equation84
(84)
(84) ) it follows
which results in the explicit expression (Equation87
(87)
(87) ). The formula for pressure and stresses follows respectively from (Equation85
(85)
(85) ), (Equation86
(86)
(86) ) and (Equation63
(63)
(63) )–(Equation65
(65)
(65) ).
The asymptotic stresses corresponding to the dominate singular term in (Equation87
(87)
(87) ) are depicted around the crack tip in Figure for
in plots (a)–(c), and for
in plots (d)–(f).
4.2. Crack under mixed Dirichlet–Neumann impermeability condition
Because in the polar coordinates the normal displacement at the crack is , and the tangential stress is
, the mixed Dirichlet–Neumann boundary conditions in (Equation6
(6)
(6) ) imply the following pairs of equations
and
. According to their representation in (Equation71
(71)
(71) ) and (Equation72
(72)
(72) ), we get two matrix equations:
(88)
(88) which solvability necessitates again the condition (Equation79
(79)
(79) ), half-integer powers (Equation80
(80)
(80) ), and admissible solutions
(89)
(89) Here, if
, then
holds together with
in (Equation88
(88)
(88) ).
As m = 1, equations and
according to (Equation75
(75)
(75) ) and (Equation76
(76)
(76) ) read
(90)
(90) and they are solvable when
(91)
(91) From (Equation90
(90)
(90) ) and (Equation91
(91)
(91) ) we infer that no logarithms occur except for
.
Theorem 4.2
(Power series for the impermeability crack problem). The Stokes equation (Equation18(18)
(18) ) with the crack under mixed Dirichlet–Neumann impermeability conditions (Equation6
(6)
(6) ) have the power series solution:
(92)
(92) since
, and
is determined according to (Equation56
(56)
(56) ) and (Equation57
(57)
(57) ) at
. For
, the functions
and
are given in (Equation85
(85)
(85) ) involving the six coefficients
as follows:
(93)
(93) For the variational solution n>0, the crack singularity has the explicit form as n = k = 1 in (Equation93
(93)
(93) ):
(94)
(94) the pressure and stresses
which are depicted in Figure for
in plots (a)–(c), and for
in plots (d)–(f).
4.3. Stress-free crack under Neumann conditions
The four equations constituting homogeneous Neumann boundary conditions (Equation7(7)
(7) ) are decoupled into two pairs, using expressions (Equation72
(72)
(72) ) and (Equation73
(73)
(73) ) as m = 0 for
and
:
(95)
(95) Zero matrix determinant in (Equation95
(95)
(95) ) leads to γ in (Equation79
(79)
(79) ) and (Equation80
(80)
(80) ), such that admissible coefficients are
(96)
(96) where
as
was used from Theorem 3.1.
For m = 1, equations and
in (Equation76
(76)
(76) ) and (Equation77
(77)
(77) ) read:
(97)
(97) The inhomogeneous system (Equation97
(97)
(97) ) is solvable when
(98)
(98) Together with (Equation96
(96)
(96) ), conditions (Equation98
(98)
(98) ) imply
for
and the next theorem.
Theorem 4.3
(Power series for the Neumann crack problem). The Stokes equation (Equation18(18)
(18) ) under the stress-free crack conditions (Equation7
(7)
(7) ) have the power series solution:
(99)
(99) where
is determined according to (Equation56
(56)
(56) ) and (Equation57
(57)
(57) ). The angular functions for
are defined in (Equation85
(85)
(85) ) and (Equation60
(60)
(60) ) as
, with the coefficients determined according to formula (Equation96
(96)
(96) ) for
:
(100)
(100) For the variational solution as n>0, the crack-tip singularity as
is expressed explicitly:
(101)
(101) when
and
are set for n = k = 1 in (Equation100
(100)
(100) ), and
in (Equation60
(60)
(60) ). The pressure
and stresses are
as depicted around the crack tip in Figure for
in plots (a)–(c), and for
in plots (d)–(f).
4.4. Crack under non-penetration conditions
The boundary conditions describing non-penetration between the crack faces (Equation8(8)
(8) ) in the polar coordinates at m = 0 yield as in the Neumann case
, hence the second system in (Equation95
(95)
(95) ) providing the determinant in (Equation79
(79)
(79) ) and the half-integers powers in (Equation80
(80)
(80) ). The other conditions consisting of
and complementarity relations
,
,
, using (Equation71
(71)
(71) )–(Equation73
(73)
(73) ) yield
(102)
(102) The second system in (Equation95
(95)
(95) ) for
and nonlinear relations (Equation102
(102)
(102) ) for
have admissible solutions
(103)
(103) where
at
follows
in (Equation102
(102)
(102) ), and
from the second system in (Equation95
(95)
(95) ).
At m = 1, linear conditions imply the second system in (Equation97
(97)
(97) ). According to (Equation98
(98)
(98) ) we find
for
or
, and
for
and
. Together with (Equation103
(103)
(103) ) this means that logarithms for the eigenvectors
and
may appear only if
.
The nonlinear conditions and
,
,
in (Equation75
(75)
(75) ) and (Equation77
(77)
(77) ), after division by 2 and
are
(104)
(104) If
in (Equation104
(104)
(104) ) for
, using
from the first string in (Equation103
(103)
(103) ) shortens it to
which leads to the complementary cases: either
(following
due to (Equation103
(103)
(103) )) or
(105)
(105) But the last inequality in (Equation105
(105)
(105) ) contradicts to the opposite inequality in (Equation103
(103)
(103) ) when
. Thus, (Equation105
(105)
(105) ) is possible and may lead to the logarithms only for
.
If for
and
, from conditions (Equation104
(104)
(104) ) we have
(106)
(106) The last inequality in (Equation106
(106)
(106) ) contradicts to the opposite inequality in the second string of (Equation103
(103)
(103) ) when
. For
, this complementarity holds when either
and for
when
and
due to (Equation59
(59)
(59) ), thus excluding logarithms. For
, the relations (Equation104
(104)
(104) ) hold when
. Thus, logarithms for the eigenvectors
and
may appear for the powers
and
. We summarize the results in the next theorem.
Theorem 4.4
(Power series for the non-penetration crack problem). The Stokes equation (Equation18(18)
(18) ) describing the non-penetrating crack under unilateral conditions (Equation8
(8)
(8) ) possess the power series solution, for
:
(107)
(107)
(108)
(108) where reduced functions
and
for the powers
and
are defined by the recurrence:
(109)
(109) For the variational solution as k>0, the singularity of the non-penetrating crack is the same as in (Equation101
(101)
(101) ) for the stress-free crack, the next asymptotic term is of order
, and the coefficients
are such that
(110)
(110)
4.5. Shear crack under transmission with slip conditions
The transmission with slip boundary conditions at the crack (Equation9(9)
(9) ) expressed in the polar coordinates as
and
, using formulas (Equation71
(71)
(71) )–(Equation73
(73)
(73) ), are
(111)
(111) The solvability of (Equation111
(111)
(111) ) requires that
, that holds for
,
, and admissible coefficients
(112)
(112) For m = 1, equations
and
from (Equation75
(75)
(75) )–(Equation77
(77)
(77) ) imply
(113)
(113) Since
, the inhomogeneous system (Equation97
(97)
(97) ) is solvable when
(114)
(114) Conditions (Equation112
(112)
(112) ) and (Equation114
(114)
(114) ) together imply that
and no log-oscillations for
.
Theorem 4.5
(Power series for the share crack problem). The Stokes equation (Equation18(18)
(18) ) under the transmission boundary conditions with slip at the crack (Equation9
(9)
(9) ) possess the power series solution:
(115)
(115) where
is given in (Equation56
(56)
(56) ) and (Equation57
(57)
(57) ) as n = 0. The functions
and
are expressed by the sum
(116)
(116) of the eigenvectors from (Equation60
(60)
(60) ) as
with the coefficients satisfying relations (Equation112
(112)
(112) ) and (Equation59
(59)
(59) ).
For the variational solution as n>0, the crack-tip singularity is described by the asymptotic formula
(117)
(117) the pressure
and stresses
,
,
.
5. Concluding remarks
The main purpose of the analytical derivation of asymptotic solutions near the crack tip is to describe singular behaviour for the underlying crack-like problems. Thus, the principal asymptotic term obtained from the power series answers the questions of square-root and log-oscillations. The important point for analysis is that the derived singular term restricts maximal regularity of variational solutions for the Stokes problem, which can be expressed in weighted Sobolev spaces, see Ref. [Citation22, Section 5.8].
The power series solutions for the Stokes equation obtained in our work are validated theoretically by the rigorous asymptotic analysis. The main advantage concerns nonlinear and non-standard boundary conditions which can be treated within our approach. In comparison with the classical Dirichlet and Neumann problems, the present solutions coincide with those ones known from the literature. For numerical techniques which are suitable for the singular field modelling, we cite the methods of potential [Citation23], locking-free elements [Citation24], extended finite element [Citation25] and boundary element methods [Citation26,Citation27].
The example solutions are presented for the crack that is a void in a solid body. A crack-like object in an incompressible fluid can describe flow around a thin plate subject to boundary conditions imposed on the plate faces, e.g. the fluid adhesion [Citation28]. In the general case, a thin plate in a body can be considered as an inclusion made from another phase of material (called anti-crack, stiffener, defect). The problem definition was extended to inclusions in the works by Khludnev and Popova [Citation29] with the help of Robin-type boundary conditions:
where
is the rigidity parameter. The limit case
describes a rigid inclusion corresponding to the Dirichlet conditions (Equation5
(5)
(5) ), and the crack-free state is recovered as
. The asymptotic solution was extended to a poroelastic body with a fluid-filled crack in Ref. [Citation30]. For the other possible nonlinear boundary conditions due to the stick–slip, we refer to Ref. [Citation5].
Acknowledgments
V.A.K. thanks the Austrian Science Fund (FWF) project P26147-N26 (PION); K.O. is supported by the Research Institute for Mathematical Sciences.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Correction Statement
This article has been republished with minor changes. These changes do not impact the academic content of the article.
References
- Ladyzhenskaya OA. The mathematical theory of viscous incompressible flow. New York: Science Publishers; 1969.
- Kikuchi N, Oden JT. Contact problems in elasticity: a study of variational inequalities and finite element methods. Philadelphia, PA: SIAM; 1988.
- Beneš M, Kučera P. Solutions to the Navier–Stokes equations with mixed boundary conditions in two-dimensional bounded domains. Math Nachrichten. 2016;289(2–3):194–212.
- Dauge M. Stationary Stokes and Navier–Stokes systems on two- or three-dimensional domains with corners. Part I. Linearized equations. SIAM J Math Anal. 1989;20(1):74–97.
- Haslinger J, Kučera R, Sassi T, et al. Dual strategies for solving the Stokes problem with stick–slip boundary conditions in 3D. Math Comput Simul. 2021;189:191–206.
- Khludnev AM, Kovtunenko VA. Analysis of cracks in solids. Southampton: WIT-Press; 2000.
- González Granada JR, Kovtunenko VA. A shape derivative for optimal control of the nonlinear Brinkman–Forchheimer equation. J Appl Numer Optim. 2021;3(2):243–261.
- Itou H, Kovtunenko VA, Rajagopal KR. On an implicit model linear in both stress and strain to describe the response of porous solids. J Elasticity. 2021;144(1):107–118.
- González Granada JR, Gwinner J, Kovtunenko VA. On the shape differentiability of objectives: a Lagrangian approach and the Brinkman problem. Axioms. 2018;7(4):76.
- Kovtunenko VA. Sensitivity of cracks in 2D-Lamé problem via material derivatives. Z Angew Math Phys. 2001;52(6):1071–1087.
- Kovtunenko VA, Ohtsuka K. Shape differentiability of Lagrangians and application to Stokes problem. SIAM J Control Optim. 2018;56:3668–3684.
- Kovtunenko VA, Ohtsuka K. Shape differentiability of Lagrangians and application to overdetermined problems. In: Itou H, Hirano S, Kimura M, Kovtunenko VA, Khludnev AM, editors. Mathematical analysis of continuum mechanics and industrial applications III (Proc. CoMFoS18), Mathematics for Industry Vol. 34. Singapore: Springer; 2020. p. 97–110.
- Kovtunenko VA, Ohtsuka K. Inverse problem of shape identification from boundary measurement for Stokes equations: shape differentiability of Lagrangian. J Inv Ill-Posed Prob. 2022;30. https://doi.org/10.1515/jiip-2020-0081.
- Williams ML. On the stress distribution at the base of a stationary crack. J Appl Mech. 1956;24(1):109–114.
- Kondrat'ev VA, Oleinik OA. Boundary-value problems for partial differential equations in non-smooth domains. Russ Math Surv. 1983;38(2):1–86.
- Itou H, Khludnev AM, Rudoy EM, et al. Asymptotic behaviour at a tip of a rigid line inclusion in linearized elasticity. Z Angew Math Mech. 2012;92(9):716–730.
- Itou H, Kovtunenko VA, Tani A. The interface crack with Coulomb friction between two bonded dissimilar elastic media. Appl Math. 2011;56(1):69–97.
- Khludnev AM, Kozlov VA. Asymptotics of solutions near crack tips for Poisson equation with inequality type boundary conditions. Z Angew Math Phys. 2008;59(2):264–280.
- Duduchava R, Wendland WL. The Wiener–Hopf method for systems of pseudodifferential equations with an application to crack problems. Integr Equ Oper Theory. 1995;23:294–335.
- Costabel M, Dauge M. Crack singularities for general elliptic systems. Math Nachrichten. 2002;235(1):29–49.
- Stephenson RA. The equilibrium field near the tip of a crack for finite plane strain of incompressible elastic materials. J Elasticity. 1982;12(1):65–99.
- Kozlov VA, Maz'ya VG, Rossmann J. Spectral problems associated with corner singularities of solutions to elliptic equations. Providence, UT: AMS; 2001.
- Hintermüller M, Kovtunenko VA, Kunisch K. A Papkovich–Neuber-based numerical approach to cracks with contact in 3D. IMA J Appl Math. 2009;74(3):325–343.
- Belhachmi Z, Sac-Epée J-M, Tahir S. Locking-free finite elements for unilateral crack problems in elasticity. Math Model Nat Phenom. 2009;4(1):1–20.
- Wagner GJ, Moës N, Liu WK, et al. The extended finite element method for rigid particles in Stokes flow. Int J Numer Meth Eng. 2001;51(3):293–313.
- Feng W-Z, Gao L-F, Dai Y-W, et al. DBEM computation of T-stress and mixed-mode SIFs using interaction integral technique. Theor Appl Fract Mech. 2020;110:102795.
- Feng W-Z, Gao L-F, Qu M, et al. An improved singular curved boundary integral evaluation method and its application in dual BEM analysis of two- and three-dimensional crack problems. Eur J Mech A/Solids. 2020;84:104071.
- Berdnik Y, Beskopylny A. The approximation method in the problem on a flow of viscous fluid around a thin plate. Aircr Eng Aerosp Technol. 2019;91(6):807–813.
- Khludnev AM, Popova TS. Equilibrium problem for elastic body with delaminated T-shape inclusion. J Comput Appl Math. 2020;376:112870.
- Itou H, Kovtunenko VA, Lazarev NP. Asymptotic series solution for plane poroelastic model with non-penetrating crack driven by hydraulic fracture. Appl Eng Sci. 2022;10:100089.
Appendix 1
Proof of Lemma 3.3
The proof of (Equation51(51)
(51) )–(Equation55
(55)
(55) ) is leaded by induction over indexes
. As i = 0, the general solution
is the direct consequence of the second equation in (Equation49
(49)
(49) ):
which leads to representation of
by the sum of general and particular solutions from (Equation49
(49)
(49) ):
and to the similar expression of
:
according to (Equation47
(47)
(47) ). As i = 1, the solution
in the form of a linear span
follows from the second equations in (Equation49
(49)
(49) ) and (Equation50
(50)
(50) ). This again leads to the combination of general and particular solutions satisfying (Equation49
(49)
(49) ) and (Equation50
(50)
(50) ) such that
Let us suppose that equations in (Equation43
(43)
(43) ) hold for all
, that is
(A1)
(A1) Inserting the ansatz (Equation51
(51)
(51) ) into the left-hand side of Equation (Equation43
(43)
(43) ) for
, using the second equation in (Equation49
(49)
(49) ) and Equation (EquationA1
(A1)
(A1) ), after shifting the summation indexes, gives the expression
(A2)
(A2) With the help of calculus using the triangle numbers, for
:
(A3)
(A3)
(A4)
(A4) we proceed the expression (EquationA2
(A2)
(A2) ) as
which proves formula (Equation43
(43)
(43) ) for the series (Equation51
(51)
(51) ) by the induction argument.
Now let equations in (Equation44(44)
(44) ) hold for all
, such that setting
we have
(A5)
(A5) Inserting the ansatz (Equation51
(51)
(51) ) and (Equation52
(52)
(52) ) into the left-hand side of Equation (Equation44
(44)
(44) ) for
and applying Equations (Equation49
(49)
(49) ), (EquationA5
(A5)
(A5) ) and shift of the summation indexes akin to (EquationA2
(A2)
(A2) ) result in
Using similar to (EquationA3
(A3)
(A3) ) and (EquationA4
(A4)
(A4) ) identities
we continue the calculation
This equality leads to the relation between coefficients (Equation54
(54)
(54) ) and proves Equation (Equation44
(44)
(44) ).
Finally, using (Equation47(47)
(47) ) we calculate the derivative of
in (Equation53
(53)
(53) ) as
(A6)
(A6) Based on the assumption that the terms in the series (Equation52
(52)
(52) ) and (Equation53
(53)
(53) ) satisfy the equations in (Equation45
(45)
(45) ):
the sum in (EquationA6
(A6)
(A6) ) can be represented in such a way that
(A7)
(A7) Consequently, inserting (EquationA7
(A7)
(A7) ) into (EquationA6
(A6)
(A6) ) we derive the relation between coefficients (Equation55
(55)
(55) ) and Equation (Equation45
(45)
(45) ). This completes the proof.