Next Article in Journal
On a Boundary Value Problem for the Biharmonic Equation with Multiple Involutions
Next Article in Special Issue
Positive Solutions for a Singular Elliptic Equation Arising in a Theory of Thermal Explosion
Previous Article in Journal
An Intelligent Web Service Composition and Resource-Optimization Method Using K-Means Clustering and Knapsack Algorithms
Previous Article in Special Issue
Nonstationary Radiative–Conductive Heat Transfer Problem in a Semitransparent Body with Absolutely Black Inclusions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nonhomogeneous Boundary Value Problem for Steady State Navier-Stokes Equations in a Multiply-Connected Cusp Domain

by
Kristina Kaulakytė
and
Konstantinas Pileckas
*
Institute of Applied Mathematics, Vilnius University, Naugarduko Str., 24, 03225 Vilnius, Lithuania
*
Author to whom correspondence should be addressed.
Submission received: 27 July 2021 / Revised: 13 August 2021 / Accepted: 17 August 2021 / Published: 24 August 2021
(This article belongs to the Special Issue Nonlinear Boundary Value Problems and Their Applications)

Abstract

:
The boundary value problem for the steady Navier–Stokes system is considered in a 2 D multiply-connected bounded domain with the boundary having a power cusp singularity at the point O. The case of a boundary value with nonzero flow rates over connected components of the boundary is studied. It is also supposed that there is a source/sink in O. In this case the solution necessarily has an infinite Dirichlet integral. The existence of a solution to this problem is proved assuming that the flow rates are “sufficiently small”. This condition does not require the norm of the boundary data to be small. The solution is constructed as the sum of a function with the finite Dirichlet integral and a singular part coinciding with the asymptotic decomposition near the cusp point.

1. Introduction

In the paper we study the nonhomogeneous stationary boundary value problem for the Navier-Stokes equations
ν Δ u + ( u · ) u + p = f , x Ω , div u = 0 , x Ω , u = a , x Ω ,
in a multiply-connected domain Ω R 2 with a cusp point O = ( 0 , 0 ) on the boundary. We assume that Ω = j = 1 N Γ j Γ consists of N + 1 disjoint components Γ , Γ j , j = 1 , . . . , N ,
Ω = Ω 0 G = D 0 \ j = 1 N D ¯ j G , D ¯ j D 0 , j = 1 , , N ,
where Γ j = D j , j = 1 , , N , Γ = Ω \ j = 1 N Γ j and G = { x R 2 : | x 1 | φ ( x 2 ) , x 2 ( 0 , H ] } , φ ( x 2 ) = γ 0 x 2 λ , γ 0 = const , λ > 1 . Moreover, we suppose that Ω Ω 0 is C 2 (see Figure 1). In (1) the velocity vector u and the pressure function p are the unknowns while the boundary value a W 1 / 2 , 2 ( Ω ) and the external force f L 2 ( Ω ) are given; ν > 0 denotes a constant coefficient of the kinematic viscosity.
We assume that the support of a is separated from the cusp point O, i.e.,
supp a Γ Ω 0 j = 1 N Γ j .
Let
F 0 = Γ Ω 0 a · n d S , F j = Γ j a · n d S , j = 1 , , N ,
be the flow rates of the boundary value a over the outer boundary Γ Ω 0 and the inner boundaries Γ j , j = 1 , , N , where n denotes the unit vector of the outward normal to Ω . By the incompressibility of the fluid it follows that
σ ( R ) u · n d S = F 0 + j = 1 N F j , 0 < R < H ,
where σ ( R ) = ( φ ( R ) , φ ( R ) ) is a cross section of G by the straight line x 2 = R parallel to the x 1 -axis. We assume that the total flux may be nonzero, i.e., F 0 + j = 1 N F j 0 . This nonzero condition means that there is a source or sink in the cusp point O. Then, due to the geometry of the domain, the velocity vector field u necessarily has infinite Dirichlet integral Ω | u ( x ) | 2 d x = (see, e.g., [1]).
The point source/sink approach is widely used in physics, astronomy and in fluid and aerodynamics. The behaviour of solutions to the Stokes and Navier–Stokes equations in singularly perturbed domains became of growing interest during the last fifty years. There is an extensive literature concerning these issues for various elliptic problems, e.g., [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. In particular, the steady Navier–Stokes equations are studied in a punctured domain Ω = Ω 0 \ { O } with O Ω 0 assuming that the point O is a sink or source of the fluid [19,20,21] (see also [22] for the review of these results). We also mention the papers [23,24,25] where the existence of a solution (with an infinite Dirichlet integral) to the Navier–Stokes problem with a sink or source in the cusp point O was proved for arbitrary data and the papers [26,27,28] where the asymptotics of a solution to the nonstationary Stokes problem is studied in domains with conical points and conical outlets to infinity.
The existence of singular solutions to the time-periodic and initial boundary value problems for the linear Stokes and the nonlinear Navier–Stokes equations in domains with a cusp point on the boundary were studied in recent papers [29,30,31,32], where the case with a sink/source in the cusp point O was considered. In [23], the existence of a generic stationary solution with infinite Dirichlet integral was proved. However, the behaviour of the solution near the cusp point was not found. The asymptotic decomposition near the cusp point of the solution u to problem (1) was constructed and the existence of a unique solution which is represented as a sum of this decomposition and a vector field belonging to a suitable second order weighted Sobolev space is proved in [1]. In [1], it is assumed that a W 3 / 2 , 2 ( Ω ) and the results are obtained under the condition that the norm a W 3 / 2 , 2 ( Ω ) is sufficiently small.
In this paper we extend the results of [1] in two directions: first, we study the case of domains with multiply-connected boundaries and, second, we prove the existence of the solution coinciding near the cusp point with the formal asymptotic decomposition assuming only that the flow rates F 0 , F 1 , , F N of the boundary value a are sufficiently small. The proof is based on the construction of an extension of the boundary value which coincides near the cusp point with the asymptotic decomposition and allows to obtain needed a priori estimates assuming only that flow rates are sufficiently small. Note that in this case the norm of a is not obliged to be small. It is worth to mention the papers [33,34,35] where the nonhomogeneous boundary value problem for the stationary Navier–Stokes equations was studied in bounded domains with multiply-connected boundaries having C 2 -regularity.

2. Notation and Auxiliary Results

2.1. Function Spaces

We will use the letter “c” for a generic constant which numerical value or dependence on parameters is unessential to our considerations; “c” may have different values in a single computation. Vector valued functions are denoted by bold letters while function spaces for scalar and vector valued functions are denoted in the same way.
Let D be a bounded domain in R n with Lipschitz boundary. C ( D ) denotes the set of all infinitely differentiable in D functions and C 0 ( D ) is the subset of all functions from C ( D ) with compact supports in D. For given non-negative integers k and q > 1 , L q ( D ) and W k , q ( D ) denote the usual Lebesgue and Sobolev spaces; W k 1 / q , q ( D ) is the trace space on D of functions from W k , q ( D ) . W ˚ 1 , 2 ( D ) is the closure of C 0 ( D ) in W 1 , 2 -norm. J 0 ( D ) is the set of all solenoidal ( div u = 0 ) vector fields u from C 0 ( D ) and H ( D ) is the closure of J 0 ( D ) in the gradient norm · L 2 ( D ) .
Lemma 1
([36,37], Chapter 1, Lemma 1). Let D R 2 be a bounded domain. If u W 1 , 2 ( D ) , then the following estimate
u L 4 ( D ) 4 c u W 1 , 2 ( D ) 2 u L 2 ( D ) 2 c u W 1 , 2 ( D ) 4
holds. Moreover, if u W ˚ 1 , 2 ( D ) then
u L 4 ( D ) 4 2 u L 2 ( D ) 2 u L 2 ( D ) 2 c u L 2 ( D ) 4 .
Consider the domain Ω with a cusp point. We introduce a family of subdomains Ω k Ω with Lipschitz boundaries:
Ω k = Ω k 1 { x R 2 : | x 1 | φ ( x 2 ) , x 2 ( h k , h k 1 ) } = Ω k 1 ω k ,
where
h 0 = H , h k = h k 1 φ ( h k 1 ) 2 L , k = 1 , 2 , ,
and L is the Lipschitz constant for the function φ .
We write u W l o c l , q ( Ω ¯ \ { O } ) if u W l , q ( Ω k ) for k .
Lemma 2.
Let u W l o c 1 , 2 ( Ω ¯ \ { O } ) , u = 0 on Ω \ { O } . If G | φ ( x 2 ) | 2 α | u ( x ) | 2 d x < , then the integral G | φ ( x 2 ) | 2 α 2 | u ( x ) | 2 d x is finite and the following inequality
H 1 H 2 φ ( x 2 ) φ ( x 2 ) | φ ( x 2 ) | α 2 | u ( x ) | 2 d x 1 d x 2 4 π 2 H 1 H 2 φ ( x 2 ) φ ( x 2 ) | φ ( x 2 ) | α | u ( x ) | 2 d x 1 d x 2
holds, where H 1 , H 2 are any numbers from the interval [ 0 , H ]
The proof of this lemma can be found in [32] (see Lemma 2.1).

2.2. Formal Asymptotic Decomposition

The formal asymptotic decomposition U [ J ] , P [ J ] of the solution ( u , p ) of problem (1) near the cusp point O was constructed in [1]. It has the form
U 1 [ J ] ( y 1 , y 2 ) = y 2 1 U 1 , 0 ( y 1 ) + k = 1 J y 2 1 + k ( λ 1 ) U 1 , k ( y 1 ) , U 2 [ J ] ( y 1 , y 2 ) = F κ 0 y 2 λ Φ ( y 1 ) + k = 1 J y 2 λ + k ( λ 1 ) U 2 , k ( y 1 ) , P [ J ] ( y 1 , y 2 ) = F κ 0 ( 1 3 λ ) y 2 1 3 λ + y 2 1 λ Q 0 ( y 1 ) + k = 1 J y 2 1 3 λ + k ( λ 1 ) C k + y 2 1 λ + k ( λ 1 ) Q k ( y 1 ) ,
where y 1 = x 1 x 2 λ , y 2 = x 2 , Φ ( y 1 ) = 1 2 ν | y 1 | 2 γ 0 2 , the functions U 1 , k , U 2 , k , Q k are regular, and
σ ( R ) U [ J ] ( x 1 x 2 λ , x 2 ) · n ( x ) d S = F , F = F 0 + j = 1 N F j .
It was proved in [1] that U j [ J ] ( y 1 , y 2 ) , j = 1 , 2 , P [ J ] ( y 1 , y 2 ) satisfy the estimates
| l U 1 [ J ] ( y 1 , y 2 ) y 1 l | c | F | y 2 1 + l , | l U 2 [ J ] ( y 1 , y 2 ) y 1 l | c | F | y 2 λ + l , l = 0 , 1 , , | P [ J ] ( y 1 , y 2 ) | c | F | y 2 3 λ 1 .
The asymptotic decomposition U ^ [ J ] ( x ) , P ^ [ J ] ( x ) = U [ J ] ( x 1 x 2 λ , x 2 ) , P [ J ] ( x 1 x 2 λ , x 2 ) is defined in G and, by construction, div x U ^ [ J ] = 0 , U ^ [ J ] | G Ω = 0 . Moreover, it was proved in [1] that for a sufficiently large J = J ( λ ) ( J > 4 λ + 1 2 ( λ 1 ) ) holds the relation
ν Δ U ^ [ J ] + ( U ^ [ J ] · ) U ^ [ J ] + ^ P [ J ] = H ^ [ J ]
with H ^ [ J ] L 2 ( G ) . Moreover, the discrepancy H ^ [ J ] satisfies the estimate
H ^ [ J ] L 2 ( G ) c | F | + | F | 2 .

3. Extension of Boundary Value

3.1. Flux Carrier from Inner Boundaries

In this subsection we construct a solenoidal vector function having the flow rates F j on inner components of the boundaries Γ j , j = 1 , , N . We call such a function the flux carrier. The construction used below is based on ideas proposed by H. Fujita in [38] for the case of symmetric domains. In [38] such functions are called virtual drains.
First we define some auxiliary functions. Let κ ( 0 , 1 / 2 ) be a parameter. We introduce non-negative, even functions β κ ( t ) C 0 ( ; + ) such that
β κ ( t ) 1 t for 0 < t < + and β κ ( t ) = 0 , | t | 1 , 1 / t , κ | t | 1 / 2 .
Define y κ = β κ ( t ) d t . Then
y κ = + β κ ( t ) d t = 1 1 β κ ( t ) d t 2 κ 1 2 1 t d t + as κ + 0 .
Define a smooth non-negative functions s κ ( t ) = s κ ( t , δ ) = 1 y κ δ β κ ( t δ ) such that s κ C 0 ( ; + ) and supp s κ [ δ ; δ ] , where δ is a small positive number. Then
+ s κ ( t ) d t = δ δ s κ ( t ) d t = 1 .
Choose one of the domains D j , j = 1 , , N , and take two points X j D j and X j 0 Γ Ω 0 such that the line X j X j 0 intersects Γ j and Γ only at one point and, if X j X j 0 intersects other boundaries, say, Γ k m , k m = k 1 , k 2 , , K j , 0 K j N 1 , then-at even number of points (if K j = 0 , then X j X j 0 does not intersect any of Γ k m ). Let us introduce in Ω the local coordinates z ( j ) = z 1 ( j ) , z 2 ( j ) such that the origin of this coordinate system coincides with the point X j and z 2 ( j ) axis is directed over the vector X j X j 0 .
The points X j 0 and X j in the local coordinates z ( j ) have the form X j 0 = ( 0 , Z 0 ( j ) ) , Z 0 ( j ) > 0 and X j = ( 0 , 0 ) . Let us take a small number μ 0 > 0 and define the strip:
Υ ( j ) = [ δ j , δ j ] × [ 0 , Z 0 ( j ) + μ 0 ] ,
where we choose a small number δ j so that the segments { z ( j ) : z 1 ( j ) = ± δ j } Υ ( j ) intersect Γ j and Γ only at one point and if intersect other boundaries, then - at even number of points.
In Υ ( j ) Ω we define a vector field:
b j ( z ( j ) ) = 0 , F j s κ ( z 1 ( j ) ) .
Notice that b j ( z ( j ) ) defined on Υ ( j ) Ω can be extended by zero into the whole domain Ω , because the bottom of Υ ( j ) is outside the domain Ω . For the sake of simplicity we keep the same notation for this extension, i.e., in the whole domain Ω we have:
b j ( z ( j ) ) = 0 , F j s κ ( z 1 ( j ) ) in Υ ( j ) Ω , ( 0 , 0 ) in Ω ¯ \ Υ ( j ) .
We shall show that
Γ i b j ( z ( j ) ) · n d S = F i , i = j , 0 , i j .
Let us introduce the domain Υ ˜ ( j ) Υ ( j ) Ω with the boundary Υ ˜ ( j ) which is the union of: Γ j Υ ( j ) , { z ( j ) : δ j z 1 ( j ) δ j , z 2 ( j ) = μ j * } and the lines z 1 ( j ) = ± δ j , where μ j * > 0 is a such small number that Υ ˜ ( j ) is a simple connected set (see Figure 2). Since, due to the construction, b j ( z ( j ) ) is solenoidal and b j ( z ( j ) ) | z 1 ( j ) = ± δ j = 0 , we get 0 = Υ ˜ ( j ) div b j ( z ( j ) ) d z ( j ) = Υ ˜ ( j ) b j ( z ( j ) ) · n d S = Υ ( j ) Γ j b j ( z ( j ) ) · n d S + { z ( j ) : δ j z 1 ( j ) δ j , z 2 ( j ) = μ j * } b j ( z ( j ) ) · e 2 d S = Γ j b j ( z ( j ) ) · n d S + δ j δ j 0 , F j s κ ( z 1 ( j ) ) · ( 0 , 1 ) d z 1 ( j ) = Γ j b j k ( z ( j ) ) · n d S F j δ j δ j s κ ( z 1 ( j ) ) d z 1 ( j ) , where the vector field n denotes the unit outward normal to Ω on Γ j , while the vector e 2 denotes the unit normal to Υ ˜ ( j ) on { z ( j ) : δ j z 1 ( j ) δ j , z 2 ( j ) = μ j * } . Due to (10), from the last equality we get (11). Notice that for the case i j , when Υ ( j ) does not intersect or touch Γ i , the vector field b j vanishes on Γ i (by construction). Otherwise, if Υ ( j ) intersects Γ i at even number of points, then flow rates of b j across Γ i are equal to zero: the flow rates of b j over not intersecting parts of Υ ( j ) Γ i cancel each other.
In order to rewrite vector field b j ( z ( j ) ) in global coordinates let us take the orthogonal matrix A j with det A j = 1 such that z ( j ) = A j ( x x 0 . ) Then it is easy to verify that
b j ( i n n ) ( x ) = A j T b j ( z ( j ) ) | z ( j ) = A j ( x x 0 ) .
Therefore, the flux carrier from the inner boundaries has the form:
b ( i n n ) ( x ) = j = 1 N b j ( i n n ) ( x ) .
Lemma 3.
The vector field b ( i n n ) is smooth and solenoidal. Moreover, supp b ( i n n ) Ω 0 ,
Γ j b ( i n n ) · n d S = F j , j = 1 , , N , Γ Ω 0 b ( i n n ) · n d S = j = 1 N F j
and the following estimate
| b ( i n n ) ( x ) | + | b ( i n n ) ( x ) | c j = 1 N | F j | , x Ω ¯ 0
holds.

3.2. Flux Carrier from the Outer Boundary

The boundary condition u = a is prescribed on Γ Ω 0 j = 1 N Γ j . After subtracting the constructed flux carrier b ( i n n ) , which “removes” the fluxes F j from the inner boundaries Γ j , j = 1 , , N , we get a modified boundary value a 1 = a b ( i n n ) | Ω such that supp a 1 Γ Ω 0 j = 1 N Γ j and the flow rates of a 1 over the inner boundaries Γ j , j = 1 , , N , are equal to zero:
Γ j a 1 · n d S = 0 ,
and the flow rate of a 1 over the outer boundary Γ Ω 0 is equal to F 0 + j = 1 N F j :
Γ Ω 0 a 1 · n d S = F 0 + j = 1 N F j = F .
Now we remove the nonzero flux from the outer boundary Γ Ω 0 . For this we will need the notion of Stein’s regularised distance. Let M be a closed set in R 2 . Stein’s regularised distance Δ M ( x ) from the point x to the set M is an infinitely differentiable function in R 2 \ M and the following inequalities
a 1 d M ( x ) Δ M ( x ) a 2 d M ( x ) , | D α Δ M ( x ) | a 3 d M 1 | α | ( x ) ,
hold, where d M ( x ) = d i s t ( x , M ) is the distance from x to M . The positive constants a 1 , a 2 and a 3 are independent of M (see [39], Chapter VI, Sections 1 and 2, 167–171, Theorem 2).
Let γ be a smooth simple curve, which intersects the outer boundary at some point x ( o u t ) Γ Ω 0 , does not intersect or touch any inner boundary Γ j , j = 1 , , N , and coincides with the straight line x 1 = 0 in G (see Figure 3).
Let us introduce a function
ξ ( x ) = Ψ ln ρ ( Δ γ ( x ) ) Δ Ω \ ( Γ Ω 0 ) ( x ) ,
where Ψ and ρ are infinitely differentiable monotonic functions such that Ψ ( t ) = 0 for t 0 , Ψ ( t ) = 1 for t 1 , and ρ ( τ ) = a 1 2 d 0 for τ a 2 2 , ρ ( τ ) = τ for τ a 2 d 0 , d 0 is the distance between the curve γ and Ω \ ( Γ Ω 0 ) . The functions Δ γ ( x ) and Δ Ω \ ( Γ Ω 0 ) are regularised distances from x to γ and Ω \ ( Γ Ω 0 ) , respectively.
Lemma 4.
The function ξ ( x ) vanishes at those points x Ω ¯ \ { O } , where ρ ( Δ γ ( x ) ) Δ Ω \ ( Γ Ω 0 ) ( x ) , and ξ ( x ) = 1 at points x Ω ¯ \ { O } where Δ Ω \ ( Γ Ω 0 ) ( x ) e 1 ρ ( Δ γ ( x ) ) . Moreover, the following inequalities
ξ ( x ) x k c Δ Ω \ ( Γ Ω 0 ) ( x ) c 1 , 2 ξ ( x ) x k x l c Δ Ω \ ( Γ Ω 0 ) 2 ( x ) c 1 , x Ω ¯ 4 ,
hold with the constant c 1 dependent only on a 1 , a 2 and d 0 .
The proof of this lemma can be found in [40] (see Lemma 2).
Let us define a vector field
b ( o u t ) ( x ) = F ξ ˜ ( x ) x 2 , ξ ˜ ( x ) x 1 ,
where ξ ˜ ( x ) coincides with ξ ( x ) on the right side of the curve γ and ξ ˜ ( x ) = 0 on the left of γ .
By construction, the vector field b ( o u t ) ( x ) is smooth, solenoidal and b ( o u t ) ( x ) | Ω \ ( Γ Ω 0 ) = 0 .
Lemma 5.
There hold the relation
Γ Ω 4 b ( o u t ) · n d S = F
and the estimate
| b ( o u t ) ( x ) | + | b ( o u t ) ( x ) | c | F | x Ω ¯ 4 .
Proof. 
Since div b ( o u t ) = 0 , we have
Γ Ω 4 b ( o u t ) · n d S = σ ( h 4 ) b ( o u t ) · n d S = F φ ( h 4 ) φ ( h 4 ) ξ ˜ ( x ) x 2 , ξ ˜ ( x ) x 1 · ( 0 , 1 ) d x = F φ ( h 4 ) φ ( h 4 ) ξ ˜ ( x ) x 1 d x 1 = F ξ ˜ ( φ ( h 4 ) , h 4 ) ξ ˜ ( φ ( h 4 ) , h 4 ) = F .
Estimate (16) follows from Lemma 4 and properties of the regularised distance. □
The modified boundary value a 2 = a 1 b ( o u t ) | Ω has a support on Γ Ω 0 j = 1 N Γ j and the flow rates of a 2 on Γ Ω 0 and Γ j , j = 1 , , N , are equal to zero:
Γ Ω 0 a 2 · n d S = 0 , Γ j a 2 · n d S = 0 , j = 1 , , N .

3.3. Extension of a 2

The extension of the boundary value function having zero flux over the boundary was constructed by O.A. Ladyzhenskaya (see [37], Chapter V, Section 4, 127–128). To be more precise, in [37] was proved the following result
Lemma 6.
Let D R 2 be a bounded domain with Lipschitz boundary D , L D , meas ( L ) > 0 . Assume that the vector field h W 1 / 2 , 2 ( D ) satisfies the conditions L h · n d S = 0 , supp h L . Then h can be extended inside D in the form
H ( x , ε ) = ( χ ( x , ε ) E ( x ) ) x 2 , ( χ ( x , ε ) E ( x ) ) x 1 ,
where E W 2 , 2 ( D ) , ( E ( x ) x 2 , E ( x ) x 1 ) | D = h and χ = χ ( x , ε ) is Hopf’s type cut-off function, i.e., χ is smooth, χ ( x , ε ) = 1 on L , supp χ is contained in a small neighborhood of L and
| χ ( x , ε ) | ε c dist ( x , L ) .
The constant c in (19) is independent of ε > 0 .
The vector field H W 1 , 2 ( D ) is solenoidal, H | D = h , supp H is contained in a small neighbourhood of L and there holds the estimate
H W 1 , 2 ( D ) c ( ε ) h W 1 / 2 , 2 ( D ) .
Moreover, for any ε > 0 the vector field H satisfies the Leray-Hopf inequality
D w · H · w d x c ε D | w | 2 d x w H ( D )
with the constant c independent of ε.
Because of the condition (17) we can apply Lemma 6 to a 2 and we obtain the following result.
Lemma 7.
There exists a vector field b 0 W 1 , 2 ( Ω ) such that b 0 | Ω = a 2 , div b 0 = 0 , supp b 0 Ω ¯ 1 ,
b 0 W 1 , 2 ( Ω ) c a 2 W 1 / 2 , 2 ( Ω ) .
Moreover,
Ω k w · b 0 · w d x c ε Ω k | w | 2 d x w H ( Ω k ) .
The constant c in (23) is independent of ε > 0 and k.

3.4. Construction of Extension Coinciding with Asymptotic Decomposition near Cusp Point

Now we “glue” the above constructed vector field B = b 0 + b ( i n n ) + b ( o u t ) with the asymptotic decomposition U [ J ] .
Let ζ be a smooth cut-off function such that ζ ( x 2 ) = 1 for x 2 h 3 , ζ ( x 2 ) = 0 for x 2 h 4 , 0 ζ ( x 2 ) 1 . We put
A = b 0 + b ( i n n ) + ζ b ( o u t ) + 1 ζ U [ J ] + V [ J ] ,
where V [ J ] is the solution of the following problem
div V [ J ] = ζ · b ( o u t ) + ζ · U [ J ] , x ω 4 , V [ J ] = 0 , x ω 4 .
Notice that ω 4 ζ · b ( o u t ) ζ · U [ J ] d x = ω 4 div b ( o u t ) + 1 ζ U [ J ] d x = 0 . Indeed,
ω 4 div ζ b ( o u t ) + 1 ζ U [ J ] d x = ω 4 ζ b ( o u t ) · n d S + ω 4 1 ζ U [ J ] · n d S = σ ( h 3 ) b ( o u t ) · n d S σ ( h 4 ) U [ J ] · n d S = F 0 + j = 1 N F j φ ( h 4 ) φ ( h 4 ) U 2 [ J ] d x 1 = 0 ,
where we used the fact that b ( i n n ) = 0 in ω 4 . Therefore, there exists a solution V [ J ] W ˚ 1 , 2 ( ω 4 ) of problem (25) satisfying the estimate
V [ J ] L 2 ( ω 4 ) c ζ · b ( o u t ) + U [ J ] L 2 ( ω 2 ) c | F | ,
see [41].
Since div b 0 + b ( i n n ) = 0 and supp b 0 + b ( i n n ) Ω ¯ 4 , from the construction we conclude the following result.
Lemma 8.
The vector field A W l o c 1 , 2 ( Ω ¯ \ { O } ) satisfies the boundary condition A | Ω = a , A is solenoidal and A ( x ) = U [ J ] ( x ) for x 2 h 4 .

4. Existence and Uniqueness of Weak Solution

In this section we prove the existence of the weak solution of problem (1).
First assume that u , p is a classical solution of (1). Multiplying (1)1 by the test function η C 0 ( Ω ) and integrating by parts, we obtain
ν Ω u · η d x + Ω u · u · η d x Ω p div η d x = Ω f · η d x η C 0 ( Ω ) .
We look for the solution u , p in the form
u = A + v , p = ( 1 ζ ) P [ j ] + p ˜ ,
where A is the extension of the boundary value a constructed in the previous section, P [ j ] is defined by (6)3 and v H ( Ω ) . Substituting (28) into (27) we obtain
ν Ω v · η d x + Ω ( v + A ) · v · η d x + Ω v · A · η d x = f ^ , η Ω η J 0 ( Ω ) ,
where
f ^ , η Ω = Ω f · η d x ν Ω A · η d x Ω A · A · η d x Ω ( 1 ζ ) P [ J ] · η d x .
The vector field v will be found as a limit of the sequence { v k } , where v k are weak solutions in the domains Ω k , that is, the vector fields v k H ( Ω k ) , satisfy the integral identities
ν Ω k v k · η d x + Ω k ( v k + A ) · v k · η d x + Ω k v k · A · η d x = f ^ , η Ω k η H ( Ω k ) .
Theorem 1.
Let f L 2 ( Ω ) , a W 1 / 2 , 2 ( Ω ) and supp a Γ Ω 0 j = 1 N Γ j . There exists a number F 0 > 0 such that if
j = 0 N | F j | F 0 ,
then problem (31) admits at least one solution v k H ( Ω k ) . There holds the estimate
v k L 2 ( Ω ) 2 c f L 2 ( Ω k ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4
with the constant c independent of k.
Proof. 
It is well known (see [37]) that integral identity (31) is equivalent to the operator equation
v k = B v k
with a completely continuous operator B : H ( Ω k ) H ( Ω k ) , defined by the relation
[ B v k , η ] Ω k = ν 1 Ω k ( v k + A ) · v k · η d x ν 1 Ω k v k · A · η d x + ν 1 f ^ , η Ω k ,
where [ w , η ] Ω k = Ω k w · η d x is the scalar product in H ( Ω k ) .
So, the solvability of Equation (34) will follow from the Leray–Schauder theorem provided we prove that the norms of all possible solutions of the operator equations
v k ( λ ) = λ B v k ( λ ) , λ [ 0 ; 1 ] ,
are bounded by a constant independent of λ .
Operator Equation (35) is equivalent to the identity
ν Ω k v k ( λ ) · η d x + λ Ω k ( v k ( λ ) + A ) · v k ( λ ) · η d x + λ Ω k v k ( λ ) · A · η d x = λ f ^ , η Ω k η H ( Ω k ) .
Taking in (36) η = v k ( λ ) we obtain
ν Ω k | v k ( λ ) | 2 d x = λ Ω k v k ( λ ) · A · v k ( λ ) d x + λ f ^ , v k ( λ ) Ω k .
To estimate the term λ f ^ , v k ( λ ) Ω k in the right hand side of (37), we use the representation (24) for the vector field A . We denote A 1 = b 0 + b ( i n n ) + ζ b ( o u t ) + V [ J ] and A 2 = ( 1 ζ ) U [ J ] , so that A = A 1 + A 2 . Since supp A 1 Ω 4 , using estimates (13), (16), (22) and (26), the embedding W 1 , 2 ( Ω 4 ) L 4 ( Ω 4 ) and the definition of a 2 = a b ( i n n ) + b ( o u t ) | Ω , we obtain the following inequality
λ ν | Ω 4 A 1 · v k ( λ ) d x | + λ | Ω 4 A 1 · v k ( λ ) · A 1 d x | c A 1 L 2 ( Ω k ) 2 + A 1 L 4 ( Ω k ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c j = 0 N | F j | 2 + j = 0 N | F j | 4 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 .
Since supp ( b 0 + b ( i n n ) ) Ω 3 , we have ( b 0 + b ( i n n ) ) · A 2 + A 2 · ( b 0 + b ( i n n ) ) = 0 . Thus, by (7), (16) and (26),
λ | Ω 4 A 1 · v k ( λ ) · A 2 d x | + λ | Ω 4 A 2 · v k ( λ ) · A 1 d x | | ω 4 ( V [ J ] + b ( o u t ) ) · v k ( λ ) · U [ J ] d x | + | ω 4 U [ J ] · v k ( λ ) · ( V [ J ] + b ( o u t ) ) d x | c V [ J ] L 4 ( ω 4 ) 4 + b ( o u t ) L 4 ( ω 4 ) 4 + U [ J ] L 4 ( ω 4 ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c V [ J ] W 1 , 2 ( ω 4 ) 4 + | F | 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c j = 0 N | F j | 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c a W 1 / 2 , 2 ( Ω ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 .
Further, the straightforward calculations give the equality
λ ν Ω k \ Ω 3 A 2 · v k ( λ ) d x + λ Ω k \ Ω 3 A 2 · A 2 · v k ( λ ) d x + λ Ω k \ Ω 3 ( 1 ζ ) P [ J ] · v k ( λ ) d x = λ ν Ω k \ Ω 3 ( 1 ζ ) U [ J ] · v k ( λ ) d x + λ Ω k \ Ω 3 ( 1 ζ ) U [ J ] · ( 1 ζ ) U [ J ] · v k ( λ ) d x + λ Ω k \ Ω 3 ( 1 ζ ) P [ J ] · v k ( λ ) d x = λ Ω k \ Ω 4 ν Δ U [ J ] + ( U [ J ] · U [ J ] + P [ J ] ) · v k ( λ ) d x + λ ν ω 3 ( 1 ζ ) Δ U [ J ] + 2 ζ · U [ J ] + ζ U [ J ] · v k ( λ ) d x + λ ω 3 P [ J ] ζ + ( 1 ζ ) P [ J ] · v k ( λ ) d x + λ ω 3 ( 1 ζ ) 2 U [ J ] · U [ J ] · v k ( λ ) ( 1 ζ ) U [ J ] · ζ U [ J ] · v k ( λ ) d x = λ k = 1 4 J k .
Because the asymptotic decomposition ( U [ J ] , P [ J ] ) satisfies Equation (8), by (9), we obtain
λ J 1 = λ Ω k \ Ω 4 ν Δ U [ J ] + ( U [ J ] · U [ J ] + P [ J ] ) · v k ( λ ) d x = λ Ω k \ Ω 4 H [ J ] · v k ( λ ) d x c Ω k \ Ω 4 | H [ J ] | 2 d x + ε v k ( λ ) L 2 ( Ω k ) 2 c | F | 2 + | F | 4 + ε v k ( λ ) L 2 ( Ω k ) 2 c a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 .
The integrals J k , k = 2 , 3 , 4 , can be estimated using (7), and we get
λ | J 2 | + | J 3 | + | J 4 c a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 + ε v k ( λ ) L 2 ( Ω k ) 2 .
From (39)–(41) it follows that
| λ f ^ , v k ( λ ) Ω k | c f L 2 ( Ω k ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 + c ε v k ( λ ) L 2 ( Ω k ) 2 .
Substituting this estimate into (37) and choosing ε sufficiently small we obtain
ν 2 Ω k | v k ( λ ) | 2 d x λ Ω k v k ( λ ) · A · v k ( λ ) d x + c f L 2 ( Ω k ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 .
The constant c in (42) is independent of k.
Consider now the integral λ Ω k v k ( λ ) · A · v k ( λ ) d x = λ Ω k v k ( λ ) · v k ( λ ) · A d x . In virtue of Lemmas 3 and 5, we have
λ | Ω k v k ( λ ) · v k ( λ ) · b ( i n n ) + ζ b ( o u t ) d x | v k ( λ ) L 4 ( Ω k ) v k ( λ ) L 2 ( Ω k ) b ( i n n ) L 4 ( Ω 4 ) + b ( o u t ) L 4 ( Ω 4 ) c i = 0 N | F i | v k ( λ ) L 2 ( Ω k ) 2 .
By (26),
λ | Ω k v k ( λ ) · V [ J ] · v k ( λ ) d x | v k ( λ ) L 4 ( Ω k ) 2 V [ J ] L 2 ( ω 4 ) c i = 0 N | F i | v k ( λ ) L 2 ( Ω k ) 2 .
For the integral λ Ω k v k ( λ ) · v k ( λ ) · U ( [ J ] ( 1 ζ ) d x inequalities (7) and (5) yield
λ | Ω k \ Ω 3 v k ( λ ) · v k ( λ ) · U ( [ J ] ( 1 ζ ) d x | Ω k \ Ω 3 | v k ( λ ) | 2 | U ( [ J ] | 2 d x 1 / 2 v k ( λ ) L 2 ( Ω k ) c | F | h k h 3 φ ( x 2 ) φ ( x 2 ) | v k ( λ ) | 2 φ 2 ( x 2 ) d x 1 d x 2 1 / 2 v k ( λ ) L 2 ( Ω k ) c i = 0 N | F i | v k ( λ ) L 2 ( Ω k ) 2 .
Finally, using Leray–Hopf’s inequality (23), we estimate the integral Ω k v k ( λ ) · b 0 · v k ( λ ) d x :
Ω k v k ( λ ) · b 0 · v k ( λ ) d x c 0 ε Ω k | v k ( λ ) | 2 d x .
Estimates (43)–(46) yield the inequality
λ Ω k v k ( λ ) · A · v k ( λ ) d x c 0 ε + c * i = 0 N | F i | Ω k | v k ( λ ) | 2 d x λ [ 0 , 1 ] ,
where the constants c 0 and c * are independent of k and λ . Thus, estimate (42) takes the form
ν 2 Ω k | v k ( λ ) | 2 d x c 0 ε + c * F 0 Ω k | v k ( λ ) | 2 d x + c f L 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 .
Choosing ε sufficiently small, say ε = ν 8 c 0 and assuming that F 0 = ν 8 c * from the last inequality we derive
ν 2 Ω k | v k ( λ ) | 2 d x c f L 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 .
So, the norms v k ( λ ) L 2 ( Ω k ) 2 of all possible solutions v k ( λ ) of operator Equation (35) are bounded by a constant independent of λ [ 0 , 1 ] and, according to the Leray–Schauder theorem, operator Equation (34) has at least one solution v k H ( Ω k ) . Moreover, for v k holds estimate (33). □
Theorem 2.
Suppose that the conditions of Theorem 1 are fulfilled. Then problem (29) admits a solution v H ( Ω ) satisfying the following estimate
v L 2 ( Ω ) 2 c f L 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 .
Proof. 
Let us take the sequence of solutions v k constructed in Theorem 1. Extending v k by zero into Ω \ Ω k we get vector fields v ˜ k H ( Ω ) . Notice that v ˜ k satisfy integral identity (31) in which we can integrate over the domain Ω instead of Ω k . Taking an arbitrary function η J 0 ( Ω ) we can find a number k such that supp η Ω k . Since the sequence { v ˜ k } is bounded in H ( Ω ) , there exists a subsequence { v ˜ k m } which converges weakly in the space H ( Ω ) and converges strongly in L 4 ( Ω k ) for any k, as the embedding H ( Ω k ) L 4 ( Ω k ) is compact. Such subsequence can be constructed using Cantor’s diagonal argument. Then we can pass to the limit as k m in integral identity (31) taking any test function η J 0 ( Ω ) . For the limit function v H ( Ω ) we obtain the integral identity (29). Obviously, the limit function v obeys estimate (48). □
Remark 1.
Since the space J 0 ( Ω ) is dense in H ( Ω ) , integral identity (29) remains valid for every test function η H ( Ω ) .
Let us prove now the uniqueness of the solution to problem (1) having representation (28).
Theorem 3.
Let a W 1 / 2 , 2 ( Ω ) , f L 2 ( Ω ) and let the boundary value a satisfies the conditions of Theorem 1, in particular, the fluxes F i , i = 0 , 1 , , N , satisfy condition (32). Let u 1 and u 2 be two solutions of problem (1) admitting representation (28): u i = A + v i , i = 1 , 2 , with v i H ( Ω ) . There exists a number a 0 > 0 such that if
v 2 L 2 ( Ω ) a 0 ,
then the solutions u 1 and u 2 coincide.
Proof. 
Suppose problem (1) has two solutions u 1 and u 2 admitting representation (28), i.e., u 1 = A + v 1 , u 2 = A + v 2 , where v 1 , v 2 H ( Ω ) and satisfy integral identity (29). Denote v = u 1 u 2 = v 1 v 2 H ( Ω ) . Subtracting integral identity (29) for v 2 from the one for v 1 we obtain
ν Ω v · η d x + Ω v 1 · v · η d x + Ω v · v 2 · η d x + Ω A · v · η d x + Ω v · A · η d x = 0 η H ( Ω ) .
Taking in (50) η = v yields
ν Ω | v | 2 d x = Ω v · v 2 · v d x Ω v · A · v d x .
The integral Ω v · A · v d x admits the estimate
| Ω v · A · v d x | c 0 ε + c * i = 0 N | F i | Ω | v | 2 d x c 0 ε + c * F 0 Ω | v | 2 d x
see (47), (32), and
| Ω v · v 2 · v d x | v L 4 ( Ω ) 2 v 2 L 2 ( Ω ) c 1 a 0 Ω | v | 2 d x .
Then from (51) it follows
ν Ω | v | 2 d x c 0 ε + c * F 0 + c 1 a 0 Ω | v | 2 d x .
Remind that F 0 is equal to ν 8 c * (see the proof of Theorem 1). Taking ε = ν 8 c 0 and assuming that a 0 = ν 4 c 1 we get
ν 2 Ω | v | 2 d x 0 .
Thus, v 1 = v 2 . □
Remark 2.
If v 2 satisfies estimate (48), that is
v 2 L 2 ( Ω ) 2 c f L 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 ,
then (49) follows from the condition
c f L 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 2 + a W 1 / 2 , 2 ( Ω ) 4 a 0 2 .

Author Contributions

Both authors contributed equally in this article. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pileckas, K.; Raciene, A. On singular solutions of the stationary Navier–Stokes system in power cusp domains. Math. Model. Anal. 2021, 26, 982–1010. [Google Scholar]
  2. Albişoru, A.-F. On transmission-type problems for the generalized Darcy-Forchheimer-Brinkman and Stokes systems in complementary Lipschitz domains in R3. Filomat 2019, 33, 3361–3373. [Google Scholar] [CrossRef] [Green Version]
  3. Cardone, G.; Nazarov, S.A.; Sokolowski, J. Asymptotics of solutions of the Neumann problem in a domain with closely posed components of the boundary. Asymptot. Anal. 2009, 62, 41–88. [Google Scholar] [CrossRef]
  4. Kamotski, I.V.; Maz’ya, V.G. On the third boundary value problem in domains with cusps. J. Math. Sci. 2011, 173, 609–631. [Google Scholar] [CrossRef]
  5. Liu, Q. Regularity for the 3D MHD Equations via One Directional Derivative of the Pressure. Bull. Braz. Math. Soc. New Ser. 2020, 51, 157–167. [Google Scholar] [CrossRef]
  6. Maz’ya, V.G.; Plamenevskii, B.A. Estimates in Lp and Hölder classes and the Miranda–Agmon maximum principle for solutions of elliptic boundary value problems in domains with singular points on the boundary. Math. Nachr. 1978, 81, 25. [Google Scholar]
  7. Maz’ya, V.G.; Nazarov, S.A.; Plamenevskij, B.A. Asymptotic Theory of Elliptic Boundary Value Problems in Singularly Perturbed Domains; Birkhäuser-Verlag: Basel, Switzerland; Boston, MA, USA; Berlin, Germany, 2000; Volume 2. [Google Scholar]
  8. Movchan, A.B.; Nazarov, S.A. Asymptotics of the stressed-deformed state near a spatial peak-type inclusion. Mekh. Kompozit. Mater. 1985, 5, 792–800. [Google Scholar]
  9. Movchan, A.B.; Nazarov, S.A. The stressed-deformed state near a vertex of a three-dimensional absolutely rigid peak imbedded in an elastic body. Prikl. Mech. 1989, 25, 10–19. [Google Scholar]
  10. Nazarov, S.A. The structure of solutions of elliptic boundary value problems in thin domains. Vestn. Leningrad. Univ. Ser. Mat. Mekh. Astr. 1982, 7, 65–68. [Google Scholar]
  11. Nazarov, S.A. Asymptotic solution of the Navier–Stokes problem on the flow of a thin layer of fluid. Sib. Math. J. 1990, 31, 131–144. [Google Scholar] [CrossRef]
  12. Nazarov, S.A.; Pileckas, K. The Reynolds flow of a fluid in a thin three dimensional channel. Litov. Matem. Sb. 1990, 30, 772–783. [Google Scholar] [CrossRef]
  13. Nazarov, S.A.; Polyakova, O.R. Asymptotic behaviour of the stress–strain state near a spatial singularity of the boundary of the “beak ti” type. Prikl. Mat. Mekh. 1993, 57, 130. [Google Scholar]
  14. Nazarov, S.A. On the flow of water under a still stone. Math. Sb. 1995, 11, 75. [Google Scholar] [CrossRef]
  15. Nazarov, S.A.; Pileckas, K. Asymptotics of solutions to Stokes and Navier–Stokes equations in domains with paraboloidal outlets to infinity. Rend. Sem. Math. Univ. Padova 1998, 99, 1–43. [Google Scholar]
  16. Nazarov, S.A. On the essential spectrum of boundary value problems for systems of differential equations in a bounded domain with a peak. Funkt. Anal. I Prilozhen. 2009, 43, 55–67. [Google Scholar]
  17. Nazarov, S.A.; Sokolowski, J.; Taskinen, J. Neumann Laplacian on a domain with tangential components in the boundary. Ann. Acad. Sci. Fenn. Math. 2009, 34, 131–143. [Google Scholar]
  18. Ragusa, M.A.; Wu, F. Global regularity and stability of solutions to the 3D double-diffusive convection system with Navier boundary conditions. Adv. Differ. Equ. 2021, 26, 341–362. [Google Scholar]
  19. Kim, H.; Kozono, H. A removable isolated singularity theorem for the stationary Navier–Stokes equations. J. Diff. Equ. 2006, 220, 68–84. [Google Scholar] [CrossRef] [Green Version]
  20. Pukhnachev, V.V. Singular solutions of Navier–Stokes equations. Advances in Mathematical Analysis of PDEs. AMS Transl. 2014, 232, 193–218. [Google Scholar]
  21. Russo, A.; Tartaglione, A. On the existence of singular solutions of the stationary Navier–Stokes problem. Lith. Math. J. 2013, 53, 423–437. [Google Scholar] [CrossRef]
  22. Korobkov, M.B.; Pileckas, K.; Pukhnachev, V.V.; Russo, R. The flux problem for the Navier–Stokes equations. Russ. Math. Surv. 2014, 69, 115–176. [Google Scholar] [CrossRef]
  23. Kaulakyte, K.; Kloviene, N.; Pileckas, K. Nonhomogeneous boundary value problem for the stationary Navier–Stokes equations in a domain with a cusp. Z. Angew. Math. Phys. 2019, 70, 36. [Google Scholar] [CrossRef]
  24. Kaulakyte, K.; Pileckas, K. Nonhomogeneous boundary value problem for the time periodic linearized Navier–Stokes system in a domain with outlet to infinity. J. Math. Anal. Appl. 2020, 489, 1241262020. [Google Scholar] [CrossRef]
  25. Kaulakyte, K.; Kloviene, N. On nonhomogeneous boundary value problem for the stationary Navier–Stokes equations in a symmetric cusp domain. Math. Model. Anal. 2021, 26, 55–71. [Google Scholar] [CrossRef]
  26. Kozlov, V.; Rossmann, J. On the nonstationary Stokes system in a cone. J. Diff. Equ. 2016, 260, 8277–8315. [Google Scholar] [CrossRef] [Green Version]
  27. Kozlov, V.; Rossmann, J. On the nonstationary Stokes system in a cone: Asymptotics of solutions at infinity. J. Math. Anal. Appl. 2020, 486, 123821. [Google Scholar] [CrossRef] [Green Version]
  28. Kozlov, V.; Rossmann, J. On the nonstationary Stokes system in a cone (Lp-theory). J. Math. Fluid Mech. 2020, 22, 42. [Google Scholar] [CrossRef]
  29. Eismontaite, A.; Pileckas, K. On singular solutions of time-periodic and steady Stokes problems in a power cusp domain. Appl. Anal. 2018, 97, 415–437. [Google Scholar] [CrossRef]
  30. Eismontaite, A.; Pileckas, K. On singular solutions of the initial boundary value problem for the Stokes system in a power cusp domain. Appl. Anal. 2019, 98, 2400–2422. [Google Scholar] [CrossRef]
  31. Pileckas, K.; Raciene, A. Non-stationary Navier–Stokes equation in 2D power cusp domain. I. Construction of the formal asymptotic decomposition. Adv. Nonlinear Anal. 2021, 10, 982–1010. [Google Scholar] [CrossRef]
  32. Pileckas, K.; Raciene, A. Non-stationary Navier-Stokes equation in 2D power cusp domain. II. Existence of the solution. Adv. Nonlinear Anal. 2021, 10, 1011–1038. [Google Scholar] [CrossRef]
  33. Korobkov, M.B.; Pileckas, K.; Russo, R. On the flux problem in the theory of steady Navier–Stokes equations with nonhomogeneous boundary conditions. Arch. Ration. Mech. Anal. 2013, 207, 185–213. [Google Scholar] [CrossRef] [Green Version]
  34. Korobkov, M.B.; Pileckas, K.; Russo, R. An existence theorem for steady Navier–Stokes equations in the axially symmetric case. Ann. Sc. Norm. Super. Pisa Cl. Sci. 2015, 15, 233–262. [Google Scholar]
  35. Korobkov, M.B.; Pileckas, K.; Russo, R. Solution of Leray’s problem for stationary Navier–Stokes equations in plane and axially symmetric spatial domains. Ann. Math. 2015, 181, 769–807. [Google Scholar] [CrossRef] [Green Version]
  36. Adams, R.A. Sobolev Spaces; Academic Press: New York, NY, USA; San Francisco, CA, USA; London, UK, 1975. [Google Scholar]
  37. Ladyzhenskaya, O.A. The Mathematical Theory of Viscous Incompressible Flow; Gordon and Breach: New York, NY, USA, 1969. [Google Scholar]
  38. Fujita, H. On stationary solutions to Navier-Stokes equation in symmetric plane domain under general outflow condition. Pitman research notes in mathematics. In Proceedings of the International Conference on Navier-Stokes Equations, Theory and Numerical Methods, Varenna, Italy, 2–6 June 1997; pp. 16–30. [Google Scholar]
  39. Stein, E.M. Singular Integrals and Differentiability Properties of Functions; Princeton University Press: Princeton, NJ, USA, 1970. [Google Scholar]
  40. Solonnikov, V.A.; Pileckas, K. Certain spaces of solenoidal vectors and the solvability of the boundary problem for the Navier–Stokes system of equations in domains with noncompact boundaries. J. Math. Sci. 1986, 34, 2101–2111. [Google Scholar] [CrossRef]
  41. Ladyzhenskaya, O.A.; Solonnikov, V.A. On some problems of vector analysis and generalized formulations of boundary value problems for the Navier–Stokes equations. Zap. Nauchn. Semin. Leningr. Otd. Mat. Inst. Steklova 1976, 59, 81–116. [Google Scholar] [CrossRef]
Figure 1. Domain Ω = Ω 0 G .
Figure 1. Domain Ω = Ω 0 G .
Mathematics 09 02022 g001
Figure 2. The strip Υ ( 1 ) Dashed area is Υ ˜ ( 1 ) .
Figure 2. The strip Υ ( 1 ) Dashed area is Υ ˜ ( 1 ) .
Mathematics 09 02022 g002
Figure 3. Curve γ .
Figure 3. Curve γ .
Mathematics 09 02022 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kaulakytė, K.; Pileckas, K. A Nonhomogeneous Boundary Value Problem for Steady State Navier-Stokes Equations in a Multiply-Connected Cusp Domain. Mathematics 2021, 9, 2022. https://0-doi-org.brum.beds.ac.uk/10.3390/math9172022

AMA Style

Kaulakytė K, Pileckas K. A Nonhomogeneous Boundary Value Problem for Steady State Navier-Stokes Equations in a Multiply-Connected Cusp Domain. Mathematics. 2021; 9(17):2022. https://0-doi-org.brum.beds.ac.uk/10.3390/math9172022

Chicago/Turabian Style

Kaulakytė, Kristina, and Konstantinas Pileckas. 2021. "A Nonhomogeneous Boundary Value Problem for Steady State Navier-Stokes Equations in a Multiply-Connected Cusp Domain" Mathematics 9, no. 17: 2022. https://0-doi-org.brum.beds.ac.uk/10.3390/math9172022

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