The qualitative analysis of solution of the Stokes and Navier-Stokes system in non-smooth domains with weighted Sobolev spaces

: This study typically emphasizes analyzing the geometrical singularities of weak solutions of the mixed boundary value problem for the stationary Stokes and Navier-Stokes system in two-dimensional non-smooth domains with corner points and points at which the type of boundary conditions change. The existence of these points on the boundary generally generates local singularities in the solution. We will see the impact of the geometrical singularities of the boundary or the mixed boundary conditions on the qualitative properties of the solution including its regularity. The solvability of the underlying boundary value problem is analyzed in weighted Sobolev spaces and the regularity theorems are formulated in the context of these spaces. To compute the singular terms for various boundary conditions, the generalized form of the boundary eigenvalue problem for the stationary Stokes system is derived. The emerging eigenvalues and eigenfunctions produce singular terms, which permits us to evaluate the optimal regularity of the corresponding weak solution of the Stokes system. Additionally, the obtained results for the Stokes system are further extended for the non-linear Navier-Stokes system.


Introduction
Let D ⊂ R 2 be a 2-dimensional bounded domain, whose boundary ∂D comprises the corner points and points at which the type of boundary conditions change.Note that a point P ∈ ∂D is said to be a corner point if there exists a neighborhood η(P) of P such that D ∩ η(P) is diffeomorphic to an angle κ intersected with the unit circle.For simplicity, we are considering a bounded plane polygonal domain (see Figure 1) with corner points (ω π) and points (ω = π) where the boundary conditions change.
The boundary points where the boundary conditions change are also referred to as corner points or vertices.The obtained results for a polygonal domain can be extended to a 2-dimensional bounded domain, i.e. (Lipschitz continuous) C 0, 1 with corner points.We considered one point as a special case of interest of corner points with an angle ω = π on one side of the domain D, where the Neumann boundary condition, the Dirichlet boundary condition, respectively, is prescribed.For the polygonal domain D with the vertices P 1 , ..., P N , we introduce the following notations.Let P N+1 = P 1 , J = 1, ..., N , Γ i (i ∈ J) is the open edge connecting the vertices P i+1 and P i , Γ 0 = Γ N , and ω i (i ∈ J) is the interior angle made by Γ i−1 , Γ i .Let J D = i ∈ J : on Γ i the Dirichlet boundary conditions are prescribed and J N = i ∈ J : on Γ i the Neumann boundary conditions are prescribed .Let N denote the set of the boundary points of D, i.e., N = P i , i ∈ J.
We assume that J D , J N are non-empty disjoint sets and J = J D ∪ J N .Moreover, let Γ 0 , Γ 1 be given by Γ 0 = i∈J D Γ i , Γ 1 = i∈J N Γ i .We have We consider the Navier-Stokes equations for an incompressible, viscous fluid, i.e.
where u = (u 1 , u 2 ) is the velocity vector field with the cartesian components u 1 , u 2 , ν is the viscosity parameter of the fluid flow, i.e., (ν > 0), q is the hydrostatic pressure and f is a given volume force density.
If the viscosity coefficient ν is sufficiently large, then the flow can be described by the following stationary Stokes system on a domain D: 2) The following mixed boundary conditions are considered on the boundary ∂D: S u, q n = h 2 on Γ 1 , ( where n = (n 1 , n 2 ) is the unit outward normal vector to the boundary and S u, q is the hydrostatic stress tensor with the cartesian components S u, q = −qδ i j + ν ∂u i ∂x j + ∂u j ∂x i . (1.5) Here, δ i j is the Kronecker symbol.Further note that if the second equation of (1.2) becoming −div u = g for a given function g satisfying the property D g dx = 0, then a particular regularity of g is required for proving the regularity of the pressure function or for handling the non-zero boundary data.Generally, for incompressible flows the function g is set equal to zero to satisfy the incompressibility condition.For simplicity, we are considering g equal to zero.Therefore, for a smooth boundary, smooth given data and boundary conditions, the system (1.2) has a smooth solution [41].The system (1.2) with the boundary conditions (1.3)- (1.4) is known as the stationary Stokes system with mixed boundary conditions [30,36].The Navier-Stokes equations or even the Stokes equations are solved for Dirichlet boundary conditions [7,10,11,15,17] but this is not common in some situations like finite channel flow models [17,28].Usually, these boundary conditions are used in the upstream of the channel and on the fixed walls but not downstream of the channel, because the downstream velocity depends on the flow in a channel which is unknown.The situation becomes more intricate when the boundary of the domain has corners or edges and the Neumann boundary conditions are applied on parts of the boundary [27,35].Therefore, the second equation of (1.2) helps to characterize the different types of Neumann boundary conditions with the Green theorem.In numerical methods, the condition (1.4) is used on the downstream boundary [13].
Presently, the corner singularity theory has been constructed for compressible viscous Stokes and Navier-Stokes systems for polygonal domains in 2-dimensions and polyhedral domains in 3dimensions.The mathematical techniques to analyze the singular structure of the solutions near corners, edges, and cusps have been discussed in [12,16,21,22,26,27,37].The key point of the corner singularity theory is the decomposition of the solution of the given problem into a regular part and a locally acting singular part which is a linear combination of explicit model singular solutions s m with unknown coefficients c m .The special singular functions s m rely on the geometry of the model problem, the differential operator, and the characteristic boundary conditions.
In the singularity expansion method for the Stokes problem, the spectral problems related to the corner singularities of solutions to elliptic equations were discussed in [7,8].Kweon [31] has considered zero Dirichlet boundary conditions to investigate the regularity results of the incompressible Navier-Stokes equations in a non-convex polygonal domain.[32,34] have extended these results for a non-convex polyhedral cylinder in R 3 with inflow boundary conditions for compressible Navier-Stokes equations.The Helmholtz decomposition to obtain regularity results of the compressible Stokes system in a non-convex polygonal domain with no-slip boundary conditions is used in [33].Recently, Anjam [4] has comprehensively discussed the singularities and regularity results of the stationary Stokes and Navier-Stokes equations on polygonal domains with convex or non-convex corners.
This study typically emphasizes analyzing the boundary singularities and regularity results of the stationary Stokes and Navier-Stokes system in two-dimensional non-smooth domains with corner points and points at which the type of boundary conditions change.We use the theory developed by Kondrat íev [24,25] and further extended by [39] for scalar problems in the context of weighted Sobolev spaces.The solvability of the considered boundary value problem is analyzed in the context of these weighted Sobolev spaces and the regularity theorems are formulated.To compute the singular terms for various boundary conditions, the Fourier transform is used to obtain the generalized form of the boundary eigenvalue problem for the stationary Stokes system.The emerging eigenvalues and eigenfunctions produce singular terms, which permits us to evaluate the optimal regularity of the corresponding weak solution of the stationary Stokes system.
The main result for the stationary Stokes system is presented in Theorem 3, and the regularity results that are direct consequences of Theorem 3 are given in Section 4.Moreover, it is proved that the weak solution (u, q) of the underlying boundary value problem belongs to W 2−γ, 2 (D) 2 × W 1−γ, 2 (D), where γ is an arbitrarily small positive real number that depends on the apex angle ω 0 .Additionally, the obtained results for the Stokes system are further extended for the non-linear Navier-Stokes system.It is proved by using the local diffeomorphism theorem that the solution of the Navier-Stokes equations has similar regularity results as the solution of the generalized Stokes problem near the corner points if the norm of the body force is sufficiently small.So far the problem has been ignored with such type of boundary conditions and domain.
The rest of this paper is organized as follows: Section 2 is devoted to present the weak formulation of the Stokes problem and introduce some function spaces.In Section 3, we determine a parametric boundary eigenvalue problem with a complex parameter ξ, the stationary Stokes system is being considered for various combinations of Dirichlet, Neumann, and mixed boundary conditions.The main regularity and expansion theorem for the stationary Stokes system is given in Theorem 3. The transcendental equations for various conditions whose zeros are the eigenvalues of the operator Û(ξ) are derived.Further, the distribution of eigenvalues and eigenfunctions are discussed.In Section 4, some regularity results for the stationary Stokes system are investigated.The results of Section 4 are further extended for non-linear Navier-Stokes system in Section 5. Section 6 is devoted to conclusions.

The weak formulation of the Stokes problem
In this section, we consider the weak formulation of the stationary Stokes problem (1.2)- (1.4).The variational formulation, solvability, and the uniqueness of the solution are offered.For a weak solution, we have restricted ourself to the homogenous boundary conditions.Denote 0 , and suppu ⊂ D. Let V m, p be a closure of E(D) in the norm of W m, p (D) 2 , 1 ≤ p < ∞ and m ≥ 0 (m need not be an integer), it is a Banach space with the norm of W m, p (D) 2 .For simplicity, we represent V 0, 2 and V 1, 2 , respectively, as H and V. Usually, these spaces are used to find the solution of the Navier-Stokes equations with homogeneous Dirichlet boundary conditions, and they are closed subspaces of the spaces L 2 (D) 2 and W 1, 2 (D) 2 .They are Hilbert spaces with the scalar products where The pressure q is a scalar function, such that, the pair (u, q) satisfy the Eqs (1.2)-(1.4) in domain D in the sense of distributions.The bilinear form a(., .) is elliptic and continuous, whereas the bilinear form b(., .) is continuous and verifies the inf-sup condition (see [15,20,41]).It is proved in [6,38] that for every f ∈ H, there exists a unique weak solution (u, q) of the boundary value problem (1.2)-(1.4)and for it the following estimate holds: where c = c(D).Hence, we have to analyze the smoothness of the weak solution (u, q) and see how it depends on the sizes of the angles ω i , i = 1, ..., N.
Remarks 1.If the given data on the right-hand sides of (1.2)-(1.4)are smoother, for example, , if the domain is sufficiently smooth and the boundary conditions do not change their types, then it is proved (see [41]) that the weak solution (u, q) ∈ [H 2 (D)] 2 × [H 1 (D)].On the other hand, the same regularity result does not hold, if the domain has corner points or points at which the type of boundary conditions changes (see [16,23]).As a matter of fact, in these cases, the regularity can be described by a decomposition of the twodimensional solution field into singular and regular parts of the form (2.5) Here, the regular part u reg belongs to [H 2 (D)] 2 × [H 1 (D)], the corner points are shown by k with the equivalent polar coordinates (r k , θ k ), the exponents ξ j, k are the eigenvalues of a Sturm-Liouville problem, and Φ j, k are the corresponding generalized eigenvector fields.

Weighted Sobolev spaces
To investigate the regularity results of the weak solution (u, q) of the corresponding boundary value problem, firstly, we introduce some function spaces in line with [1,14,24,39].
Let N be the set of the corner points and of points where the type of boundary conditions change (shortly called the singular boundary points), i.e., N ⊂ ∂D.Denote where the supp v is bounded.Let α = α 1 , ..., α N be an N−tuple of real numbers which satisfying 0 < α i < 1 for 1 ≤ i ≤ N. Therefore, the weight function is characterized by where m is an any integer and r i (x) = dist (x, P i ).We denote by D β v be the multi-index notation for higher-order derivatives and in cartesian coordinates is defined as Let W m, p α (D) be the weighted Sobolev spaces and is the closure of C ∞ N (D) equipped with the norm Analogously to the factor spaces, the trace spaces are also defined as where W m, p 0, α (D) is the closure of C ∞ 0 (D) with respect to the norm of W m, p α (D).This approach is classical for domains with conical points.Principally, for 1 ≤ p < ∞, we denote L p α (D) = W 0, p α (D).For negative integers m, i.e., m ∈ Z, m < 0, we describe the spaces W m, p α (D) as the closure of the set C ∞ N (D) equipped with the subsequent norm where q = p p−1 is known as the inverse of p.The dual space of W m, p α (D) is given as W −m, q −α (D).Similarly for the trace spaces, defining the space W m+ 1 q , p α (∂D) for m < 0, m ∈ Z, as the dual space to W −m− 1 q , q −α (∂D).Therefore, the consequent continuous imbeddings are considered directly from the definition of the above spaces (2.10) For a bounded plane domain D, we have the subsequent continuous imbeddings and provided that p 2 ≥ p 1 and the infinite strip with positive width ω 0 .For any real h > 0 and for an integer m ≥ 0, the spaces are defined as

Singularities of the Stokes problem in an infinite cone
In this section, we will see the occurrence of the singular terms of the solution of the mixed boundary value problem for the stationary Stokes problem near the corners and the structure which they have.Further, the distribution of the eigenvalues and eigenfunctions are given.

Localization and the model problem
Assume that D is a polygonal domain.To show that the weak solution (u, q) of the underlying boundary value problem is regular, we have to investigate its behaviour near the corner points P i (i ∈ J).Let us consider the point P N as origin and denote ω N = ω 0 ∈ (0, 2π).An appropriate infinite differentiable cut-off function χ(|x|) = χ(r) depending on the distance r from the point P N is defined as The number is so small that P N is the only corner point of the domain D that lies inside the circle {x : |x| ≤ 2 }.We multiply the both sides of (1.2) and (1.3)-(1.4)by the smooth cut-off function χ, then substitute (v, p) = (χu, χq) in (1.2) and likewise in (1.3)-(1.4).The derivatives are considered in the distribution sense.Thus, the boundary value problem is set into an infinite cone and coincides with the original problem near the P N .The Stokes system (1.2) becomes where The behavior of (v, p) near the corner point P N illustrates the regularity of the solution (u, q) in the neighborhood of the point P N .If we suppose  2 and G ∈ H 1 (S ).Besides, the following boundary conditions are prescribed on the subsequent edges Γ S , 0 (θ = 0) and Γ S , ω 0 (θ = ω 0 ) of the cone (see Figure 2).Just one condition is considered per edge to differentiate between the mixed boundary conditions.Therefore, the obtained boundary conditions are: Dirichlet boundary conditions: where Neumann boundary conditions: where and the notation (•) denotes the vector direct product between two vectors.Mixed boundary conditions: It is observed that the right-hand sides of the obtained boundary conditions have similar smoothness as the original problem in the domain D. To analyze the regularity results of the boundary value problem (3.1)-(3.4),we rewrite the operators in polar coordinates.Hence, the transformed form is where (v r , v θ ) are the polar components of the velocity vector v, (F r , F θ ) are the polar components of F and are given by Correspondingly, the boundary conditions (3.2)-(3.4)emerge as and , where m = 1 for Dirichlet and m = 2 for Neumann boundary conditions.They hold in the infinite cone where v(r, θ) Now, the variable τ is introduced by the relation r = e τ .Accordingly, the system (3.5) is set on the infinite strip with width ω 0 as Here, S = (τ, θ) : −∞ < τ < ∞, 0 < θ < ω 0 and ṽ = v(e τ , θ), p = e τ p(e τ , θ), F = e 2τ F(e τ , θ) and G = e τ G(e τ , θ).The Dirichlet, Neumann and mixed boundary conditions also yield the transformed form with the boundary data H l+1 = e lτ H l+1 (e τ , θ), l = 0, 1 as (3.12) To obtain the boundary eigenvalue value problem, the complex Fourier transform with respect to τ is introduced as and the inverse Fourier transform is It defines an isomorphic mapping, i.e., ) We have Moreover, it is noted that if h 1 < h 2 and the following properties are satisfied then v(ξ) is holomorphic in the strip h 1 < Im ξ < h 2 .Therefore, the relationship between the Fourier transform and the Mellin transform for any α ∈ C is given by Now, by applying (3.13) to (3.9)-(3.12)with respect to τ, the two-point boundary value problem for the unknown functions (v τ , vθ , p) is obtained.It depends on the complex parameter ξ and holds on the interval I = (0, ω 0 ).The transformed form of the problem (3.9) is given by For complex parameter ξ, we have v ∈ W 2, 2 (I) 2 , p ∈ W 1, 2 (I), F ∈ L 2 (I) 2 and Ĝ ∈ W 1, 2 (I).Let L(ξ) denote the matrix differential operator analogous to the system (3.19) and maps where The operator L(ξ) is considered for all parameter ξ ∈ C with various combinations of the boundary conditions to analyze the qualitative properties of the solution of the underlying problem near the corner points.Additionally, the Fourier transformed form of the boundary conditions is also expressed as follows: (3.23) (3.24) Additionally, the matrix boundary operators for different kinds of boundary conditions can be written as: For Dirichlet boundary conditions For Neumann boundary conditions For mixed boundary conditions Analogously, the operator B[. .](ξ) is used below to define the general transformed form of the matrix boundary operators for different kinds of boundary conditions Accordingly, the generalized form of the operator pencil Û(ξ) for the two-point boundary value problem can be written as Note that Û(ξ) can be defined for every boundary point in the sense of [2,3].So, Û(ξ)(θ, ξ) = 0 is used to describe a generalized eigenvalue problem and the solvability of these type of problems is discussed in [26].The operator Û(ξ) realizes an isomorphism for all ξ ∈ C apart from some isolated points (known as the eigenvalues of Û(ξ)).So, the resolvent R(ξ) = Û(ξ) −1 is an operator-valued, meromorphic function of ξ has poles of finite multiplicity.The eigenvalues of Û(ξ) are obtaining with the determinant method, which means that the nontrivial solution of the generalized eigenvalue problem leads a transcendental equation whose zeros are the eigenvalues of Û(ξ).
Besides, the properties of the operator Û(ξ) in the neighborhood of the corner points can be obtained by the properties of the special operator Û0 (ξ) which is explicitly given by the principal parts of the matrix differential operator L(ξ) and the matrix boundary operators B[. .](ξ).Hence, we have and Likewise, we can write for mixed boundary conditions.To compute the eigenvalues ξ µ (generally referred for multiple eigenvalues) and the corresponding eigenfunctions, we proceed as.
Remarks 2. It is noted [24][25][26] that if the complex number ξ is not an eigenvalue of the operator Û(ξ), then Û(ξ) is an isomorphism between the spaces W 2, 2 (I , and it holds for all ξ ∈ l h : with the constant c is independent of ξ.
Therefore, the Lemma 1 provides us an opportunity to prove the following theorem of the solvability of the problem (3.1)- (3.4).
Lemma 2. [29].Let Û−1 (ξ) be the inverse operator of Û(ξ).Û−1 (ξ) is a meromorphic operatorvalued function with poles which are the eigenvalues of Û(ξ).The order m of a pole ξ 0 is the largest of the lengths of the Jordan chains corresponding to ξ 0 .Moreover, the operator Û−1 (ξ) has the following expansion in the neighborhood of ξ 0 : where q i : i = 1, ..., m are the finite-dimensional operators which do not depend on ξ and Γ(ξ) is holomorphic.The operator q m behaves into the space of eigenfunctions of Û(ξ) corresponding to ξ 0 , while the operators q m−1 , ..., q 1 behave into the subspaces of the corresponding associate functions.
The next theorem describes the expansion and the regularity of the problem (3.1)-(3.4).
It is recognized for elliptic boundary value problems that the eigenvalues of the operator Û(ξ) which lies on the strip have a significant role in the regularity results.The assertions 1 and 2 of Theorem 3 represents the regularity and expansion of the solutions of the system (1.2)-(1.4)near the corner points.
Remarks 3. The technique of Mellin transform, the method of special ansatzes, and spherical coordinates are used in [26,27] to obtain the generalized form of the boundary eigenvalue problem for the stationary Stokes system with Dirichlet and mixed boundary conditions.The existence of the generalized eigenvalues is discussed in a strip Re ξ ∈ (0, 1).Here, we use the Fourier transform technique to obtain the generalized form of the boundary eigenvalue problem for the stationary Stokes system with mixed boundary conditions.Moreover, the existence of the generalized eigenvalues in a strip Im ξ ∈ (−1, 0) with various combinations of the boundary conditions that depend on the apex angle ω 0 are studied.

General solutions
Let (v τ , vθ , p) be denoting the general solution of the system (3.19) by considering the right-hand side functions equal to zero, i.e., L(ξ)(v, p) = 0, where (v τ , vθ ) stands the components of the velocity vector v and p for the pressure function.On the other hand, If the system (3.19) has the non-zero right-hand side functions, then the general solution can be written as vτ , vθ , p T = I hom + I p (ξ, θ).
Here, I p denote the particular solution corresponding to the non-zero right-hand side functions and I hom is the homogenous or general solution for zero right-hand sides.Here, our interest is to find the I hom solution.
For simplicity, we substitute ξ = −iz from [18] into (3.19),then a system of linear ordinary differential equations is obtained that depends on the complex parameter z.

The fundamental system to the ordinary differential equations
A system of linear ordinary differential equations in θ with complex parameter z is considered (3.48) Furthermore, the system (3.48)provides a linear homogenous fourth-order ordinary differential equation with constant complex coefficients, i.e., It is noted from the general theory of ordinary differential equations that (3.49) gives four independent solutions, and the general form of the fundamental solution for (z 0, ±1) or (ξ 0, ±i) can be written as Therefore, it is necessary to consider the other cases for various values of z or ξ, and the general forms of their fundamental solutions are described as follows: Case 1. (For z = 0; or ξ = 0): The coefficients B = (B 1 , B 2 , B 3 , B 4 ) T would be determined according to the types of boundary conditions.In line with the above cases, we can substitute ξ = −iz into Dirichlet, Neumann, and mixed boundary conditions.As well, their transformed forms also depend on the complex parameter z.

The calculation of the eigenvalues
To evaluate the eigenvalues and corresponding eigenfunctions of the stationary Stokes system, the solution (3.50) with various combinations of the boundary conditions is considered to obtain a system of four linear homogeneous equations with our unknowns B 1 , B 2 , B 3 and B 4 .The resulting matrix of coefficients of these equations depends on the complex parameter z, and a nontrivial solution exists if the determinant of the resulting matrix of coefficients vanishes (see [5]).Further, it produces the transcendental equations whose roots are the eigenvalues, namely, ξ µ wherein (µ is used for multiple eigenvalues, i.e., µ = 1, ..., M).To compute these results, we proceed as follows.

Dirichlet boundary conditions (DD)
It means that the Dirichlet boundary conditions are given on both sides of the corner point.The determinant method is used from [5] to obtain a system of linear homogeneous equations.A non-trivial solution exists if the determinant det D DD (z) of the corresponding system of the matrix of coefficients vanishes.So, the computation leads to the transcendental equation ) that the zeros of these equations are symmetric with respect to the origin and the real axis lies in the complex plane.Therefore, the eigenvalues of the operator Û(ξ) are positioned in the complex plane symmetrically with respect to the origin and the imaginary axis.
Proof.For Dirichlet boundary conditions The Eq (3.57) is studied in (3.54) with Dirichlet boundary conditions and is satisfied for z 0, ±1.Furthermore, the various cases of z are considered with Dirichlet boundary conditions.
For z = 0, the general solution is taken from (3.51) with Dirichlet boundary conditions and a system of linear homogeneous equations Σ B = 0 is obtained, wherein the symbol Σ denotes the matrix of coefficients.It follows from the existence of the non-trivial solution which means that the determinant of the corresponding matrix of coefficients for the linear homogeneous equations is zero.We have det(Σ) = − sin 2 (ω 0 ) + ω 2 0 > 0, (3.60) which implies that B = 0 and zero is not an eigenvalue.
For z = ±1, the general solutions of the system (3.48) is given in (3.52) and (3.53).By using Dirichlet boundary conditions, a system of linear homogeneous equations is obtained.Moreover, the determinant of the matrix of coefficients is zero and the non-trivial solution exists.However, we are not interested in the null space analogous to the eigenvalues ξ = ±i.
For Neumann boundary conditions The Eq (3.58) is calculated for Neumann boundary conditions in (3.55) and is satisfied for z 0, ±1.
For z = 0, we consider the general solution given in (3.51) with the Neumann boundary conditions, and a system of linear homogeneous equations is obtained.By the use of the method of the determinant, we get zero determinant of the matrix of coefficients.Further, the rank of the matrix of coefficients is two and (0, 0, 0, 1) T and (0, 0, 1, 0) T are two linearly independent solutions.Thus, the corresponding eigenfunctions are e 1 = (cos θ, − sin θ) T and e 2 = (sin θ, cos θ) T which represent the translation in x and y directions.
For z = 1, the general solution is given in (3.53).The determinant of the matrix of coefficients for this case is zero and the non-trivial solution exists.
For z = −1, the general solution is given in (3.52).The use of the Neumann boundary conditions gives the following determinant det(Σ) = 2 sin 2 (ω 0 ) of the corresponding matrix of coefficients.
For z = 0, the general solution is given in (3.51) and the mixed boundary conditions are used to get a system of linear homogenous equations.The determinant of the corresponding system is det(Σ) = 1 > 0, and is not an eigenvalue of Û(ξ).
For z = 1, the general solution is given in (3.53) and the determinant of the matrix of coefficients for this case is det(Σ) = − cos 2(ω 0 ).
Consequently, z = ±1, are the eigenvalues of the corresponding problem if the corresponding determinants are equal to zero.

The regularity results
Let (u, q) ∈ W 1, 2 (D) 2 × L 2 (D) be the unique weak solution of the stationary Stokes problem.It follows from the theory of Kondrat íev [24] that the pair (u, q) ∈ W 2, 2 γ+1 (D) 2 × W 1, 2 γ+1 (D), where γ is a small positive real number.To obtain further qualitative regularity results, the theory proposed by [39] is employed and we will analyze that the weak solution (u, q) ∈ W 2, 2 γ+1 (D) 2 × W 1, 2 γ+1 (D).Firstly, the case of L 2 -data is considered for the direct consequences of Theorem 3 and the observations of Section 3. Let ω 0DN denote the maximal angle from the set of all angles of such corner points wherein types of boundary conditions change.Analogously, ω 0DD and ω 0NN are the maximal angle of those corner points which have similar types of boundary conditions, i.e., Dirichlet-Dirichlet and Neumann-Neumann, on both adjacent sides of the corner point.If no such type of points exists or the angles ω 0DD , ω 0NN < π, then strongly we can set ω 0DD , ω 0NN = π as the minimum value.
The following propositions hold to formulate the regularity results of the weak solution (u, q) of the stationary Stokes system with various combinations of the boundary conditions.Lemma 3. Suppose that if the strip α − 1 ≤ Im ξ < , α > 0, is free of the zeros of the equations and ÛNN (ξ) = sin 2 (iξω 0 ) − (iξ) 2 sin 2 (ω 0 ) = 0, ξ 0, ( for any arbitrary small > 0, then the solution (u, q) ∈ W 2, 2 α (D) 2 ×W 1, 2 α (D) and satisfies the following estimate Proof.To show that no eigenvalue of (4.1) and (4.2) lies in the strip α − 1 ≤ Im ξ < for an angle ω 0 , where α is a small positive real number.Firstly, we are considering the case of Dirichlet boundary conditions ÛDD (ξ) for an apex angle ω 0 = π 2 and ω 0 = π, respectively, or for an any arbitrary angle ω 0 ∈ (0, π].We note that no eigenvalue of (4.1) is found that lie on the line h = Im ξ = t for t > −1.Consequently, (ṽ, p) ∈ [W 2, 2 α (S )] 2 × [W 1, 2 α (S )] for a small positive real number α.In addition, the singularities appear for this case at corners with an apex angle ω 0 is greater than π.(See Figure 3).In all the following graphs, the red lines indicate the pure imaginary values, while the black lines indicate the complex parts of the complex ones.Analogously, the case of Neumann boundary conditions ÛNN (ξ) is considered.There is no eigenvalue of (4.2) that lie on the line h = Im ξ = t for t > −1, for an angle ω 0 = π 2 and ω 0 = π, respectively, or for any arbitrary angle ω 0 ∈ (0, π].It produces singularities when the apex angle ω 0 is greater than π.Hence, it follows that the solution (u, q) ∈ W 2, 2 α (D) 2 × W 1, 2 α (D).Additionally, we have a bounded domain D and the corresponding continuous imbeddings (2.9) and (2.12).Let we have It is well-known that for every f ∈ V * , where V * is the dual space of V, a unique weak solution (u, q) ∈ V × L 2 (D) of some generalized steady Stokes problem exists.
Therefore, the following lemma describes the regularity results for L p -data.
Proof.The statement can be followed directly from Theorem 3 by considering p = 2 2−µ , α 2 = 0 and the similar result is obtained by applying the same process used in Lemmas 3-4.

The Navier-Stokes equations
We consider the steady Navier-Stokes equations with the homogenous mixed boundary conditions (1.3)-(1.4).If the given right-hand sides have a sufficiently small norm, then we prove by using the local diffeomorphism theorem that the Navier-Stokes equations have similar regularity results as the solution of the generalized Stokes problem near the corner points.Denote Additionally, we denote H and V are the closures of E(D) equipped with the norms of L 2 (D) 2 and W 1, 2 (D).Recall that V and H are the Hilbert spaces and their scalar products are given in (2.1).Theorem 6.Let P 1 be a continuous one-to-one operator from X onto Y, and P 2 be a compact linear operator from X into Y, where X and Y are Banach spaces.Therefore, the subsequent statements are equivalent: 1. P 1 + P 2 is one-to-one; 2. P 1 + P 2 is onto Y.
Theorem 7. Let f ∈ H, the norm of f is sufficiently small.Then a uniquely determined u ∈ M exists, such that ν (u, v) + b u, u, v = f, v ∀v ∈ V. (5.12) Proof.By (5.8), B u is the sum of operators Q and Z u .It is noted from the above results that Q : M → H is the one-to-one operator and onto H, and Z u is a compact operator.Furthermore, from Lemma 7, B u is one-to-one.These facts and Theorem 6 produces that B u is onto H. Hence, the Theorem 5 yields a unique solution.

Conclusions
In this article, we have studied the boundary singularities and regularity of the weak solution of the mixed boundary value problem for the stationary Stokes and Navier-Stokes system in a two-dimensional non-smooth domain with corner points and points at which the type of boundary conditions change.The solvability of the considered boundary value problem has been analyzed in the context of the weighted Sobolev spaces with Kondrat íev type weights and the regularity theorems are formulated.To compute the singular terms for various boundary conditions, the complex Fourier transform has been used to obtain the generalized form of the boundary eigenvalue problem for the stationary Stokes system.The emerging eigenvalues and eigenfunctions produce singular terms, which permits us to evaluate the optimal regularity of the corresponding weak solution of the stationary Stokes system.
The main regularity and expansion theorem for the stationary Stokes system is presented in Theorem 3. We have discussed the regularity results of the corresponding boundary value problem for the case of L 2 and L p -data which are the direct consequences of Theorem 3. It is seen for the case of Dirichlet and Neumann boundary conditions that if the domain D has reentrant corners (ω i > π : i = 1, 2, ...N), then the weak solution (u, q) of the considered problem produces singularities.On the other hand, for the case of mixed conditions, the singularities appear at corners with (ω i > π 4 : i = 1, 2, ...N).Moreover, it is observed that if singularities exist, then splitting the solution into a singular part which defines a linear combination of explicit model singularity functions s m for the Stokes operator with corresponding unknown coefficients C m and a regular part that belongs to H 2 × H 1 .Finally, it is proved that the weak solution (u, q) of the underlying boundary value problem belongs to W 2−γ, 2 (D) 2 × W 1−γ, 2 (D), where γ is an arbitrarily small positive real number that depends on the apex angle ω 0 .
Additionally, we have extended the obtained results for the Stokes system for the non-linear Navier-Stokes system.We have proved this by using the local diffeomorphism theorem that the solution of the Navier-Stokes equations has similar regularity results as the solution of the generalized Stokes problem near the corner points if the given body force has a sufficiently small norm.To prove this, an operator R relating to the Navier-Stokes equations is defined and has shown that it is Fréchet differentiable at the point u = 0. Furthermore, the Fréchet derivative of R at the point u is agreed with the Stokes problem.
Presently, the Stokes and Navier-Stokes equations with the Navier-slip boundary conditions and the free-boundary problems in domains with corners have very interesting phenomena.The issues regarding their existence and regularity are considered for smooth domains but theoretical results for the corner singularity decomposition are still not obtained.Therefore, these issues are numerically interesting.In future works, it is important to show the unique existence of the approximations for the regular parts and coefficients, and to derive their error estimates.On the other hand, it is also observed that the non-stationary compressible Stokes and Navier-Stokes equations on polygonal domains could be considered.

( 3 .
45)Now, substituting r = e τ , applyingLemma 2 and ( [27], Theorem 1.1.5),we obtain(3.40).Now, we derive our fundamental regularity and expansion theorem for the mixed boundary value problem for the stationary Stokes system in a two-dimensional bounded domain with corner points.By considering the substitution Re α = −Im ξ − 2, it improves the theorems ([26], Theorem 8.2.1 and Theorem 8.2.2) which are based on the Mellin transform and used for the solvability of the elliptic systems.Theorem 3. (Regularity and expansion theorem): Let α 1 and α 2 be real numbers and satisfying be a solution of the stationary Stokes system (1.2) with the homogenous Dirichlet, Neumann, and mixed boundary conditions (1.3)-(1.4)and f

Figure 3 .
Figure 3. Distribution of eigenvalues for D-D and N-N boundary conditions.

Definition 2 .Definition 3 .
Let f ∈ V * .u is called the weak solution of the problem (5.1) with the homogenous mixed boundary conditions(1.3)-(1.4),if u ∈ V and satisfies ν (u, v) + b u, u, v = f, v ∀v ∈ V. (5.2) Further, b u, v, w describes the trilinear continuous form for every u, v, w ∈ V by b u, v, w = D Let H represent the closure of E(D) in the norm L 2 (D) 2 and describe the Banach space M = u; there exists f ∈ H such that ν (u, v) = f, v for all v ∈ V .