Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Advances in Nonlinear Analysis

Editor-in-Chief: Radulescu, Vicentiu / Squassina, Marco

IMPACT FACTOR 2018: 6.636

CiteScore 2018: 5.03

SCImago Journal Rank (SJR) 2018: 3.215
Source Normalized Impact per Paper (SNIP) 2018: 3.225

Mathematical Citation Quotient (MCQ) 2018: 3.18

Open Access
See all formats and pricing
More options …

Regularity Criteria for Navier-Stokes Equations with Slip Boundary Conditions on Non-flat Boundaries via Two Velocity Components

Hugo Beirão da Veiga / Jiaqi Yang
  • Key Laboratory for Mechanics in Fluid Solid Coupling Systems, Institute of Mechanics, Chinese Academy of Sciences, Beijing, 100190, China
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2019-06-29 | DOI: https://doi.org/10.1515/anona-2020-0017


H.-O. Bae and H.J. Choe, in a 1997 paper, established a regularity criteria for the incompressible Navier-Stokes equations in the whole space ℝ3 based on two velocity components. Recently, one of the present authors extended this result to the half-space case R+3. Further, this author in collaboration with J. Bemelmans and J. Brand extended the result to cylindrical domains under physical slip boundary conditions. In this note we obtain a similar result in the case of smooth arbitrary boundaries, but under a distinct, apparently very similar, slip boundary condition. They coincide just on flat portions of the boundary. Otherwise, a reciprocal reduction between the two results looks not obvious, as shown in the last section below.

Keywords: Navier-Stokes equations; Slip boundary conditions; No flat boundaries; Two components regularity criterium

MSC 2010: 35Q30; 35B65; 76D05

1 Introduction

The starting point of the present paper is the well known Prodi-Serrin (P-S) sufficient condition for regularity of the solutions to the incompressible Navier-Stokes equations


where u = (u1, u2, u3) denotes the unknown velocity of the fluid and p the pressure. To immediately set limits to the circle of our interests, assume for now on that Ω ⊂ ℝ3 is a bounded, smooth domain, even if many results quoted below hold for larger space dimensions. For the time being, assume that suitable boundary conditions are imposed to the velocity u.

The global existence of the so called weak solutions to system (1.1) goes back to J. Leray [1] and E. Hopf [2] classical references. See also A.A. Kiselev and O.A. Ladyzhenskaya [3], and J.L. Lions [4]. Below, solutions of (1.1) are intended in this sense.

A main classical open mathematical problem is to prove, or disprove, that weak solutions are necessarily strong under reasonable but general assumptions, where strong means that


In this context, a remarkable and classical sufficient condition for uniqueness and regularity is the so-called Prodi-Serrin condition, P-S in the sequel, namely


Concerning condition (1.3), we transcribe from [5], Section 1, the following considerations: Assumption (1.3) was firstly considered by G. Prodi in his paper [6] of 1959. He proved uniqueness under this last assumption. See also C. Foias, [7]. Furthermore, J. Serrin, see [8, 9], particularly proved interior spatial regularity under the stronger (non-strict) assumption


Concerning the above problems, see also O.A. Ladyzhenskaya’s contributions [10, 11]. The above setup led to the nomenclature Prodi-Serrin condition.

Complete proofs of the strict regularity result (i.e. under assumption (1.3)) were given by H. Sohr in [12], W. von Wahl in [13], and Y. Giga in [14]. A simplified version of the proof was given in reference [15], to which we refer also for bibliography. For a quite complete overview on the main points, and references, on the initial-boundary value problem for Navier-Stokes equations we strongly recommend Galdi’s contribution [16]. Further, we refer to [9, 17], as sources for information on the historical context of the P-S condition by the initiators themselves.

Finally, we recall that L. Escauriaza, G. Seregin, and V. Šverák, see [18], extended the regularity result to the case (q, p) = (∞, 3).

A significant improvement of the P-S condition was obtained by H.-O. Bae and H.J. Choe [19], see also [20]. They proved, in the whole space case, that it is sufficient for regularity of solutions that two components of the velocity satisfy the above condition (1.3). For convenience we call here this situation as being the restricted P-S condition. In 2017, one of the authors, see [21], extended this result to the half-space R+3 under slip boundary conditions. In this case, the truncated 2-dimensional vector field ū cannot be chosen arbitrarily. The omitted component has to be the normal to the boundary.

Very recently, in reference [5], the result was extended to a cylindrical type three-dimensional domain, consisting on the complement set between two co-axial circular cylinders, with radius ρ0 and ρ1, 0 < ρ0 < ρ1, periodic in the axial direction, under the physical slip boundary condition


where D(u)=u+(u)T2 is the shear stress. The above exclusion of an interior cylinder was done to avoid the radial coordinate singularities on the symmetry axis, which consideration is out of interest in our context. Below we obtain a similar result, extended to domains with general non-flat boundaries, but under the slip boundary condition (2.1). The two boundary conditions coincide just on flat portions of the boundary. Otherwise, a reciprocal reduction between the two results looks not obvious. This claim is shown in the last section.

Again by following [5] we recall that after the contribution by H.-O. Bae and H.J. Choe, related papers appeared that particularly concerned assumptions on two components of velocity or vorticity, see [21, 22, 23, 24, 25]. There are also many papers dedicated to sufficient conditions for regularity which depend merely on one component, see, for instance, [26, 27, 28, 29, 30, 31].

Before going on we want to motivate the particular choice of the domain made below. It takes into account that the real significance of the result has essentially a local character. First of all, a global regular (i.e., without singularities) system of coordinates, two of them parallel and the third orthogonal to the boundary, does not exist in general, even in an arbitrarily thin neighbourhood of the full boundary, as in the case of a sphere and even in the case of a spherical corona. In fact, singularities typically appear, like on the above two cases, and even in full cylinders (due to the symmetry axis). The cylindrical case considered in reference [5] is an exception (see below) due to the removal of a neighbourhood of the symmetry axis.

Luckily, the above type of coordinates’ system exists in sufficient small neighbourhoods of any regular boundary point. Hence, to illustrate the full significance of our thesis in a simple, but still convincing way, it looks sufficient to prove it near any “small” piece of smooth boundary with an arbitrary geometrical shape. This is our aim below. The restrictions on the domain Ω below are made in accordance with these lines, a choice which covers the very basic situation, in the simplest way.

2 Main Results

In the sequel we assume the slip boundary condition


where ω = ∇ × u is the vorticity, n is the outward normal of ∂Ω, and Ω ∈ ℝ3 is a smooth domain satisfying the following condition:

Assumption 2.1

There exists a curvilinear orthogonal system of coordinates


such that Ω can be transited into


where the axis q3 direct to the outward normal on the boundary ∂Ω̂1 := {(q1, q2, q3) : q3 = ρ1} (the inward normal on the boundary ∂Ω̂0 := {(q1, q2, q3) : q3 = ρ0}, respectively), and q1, q2 are periodic.

Remark 2.1

The above “small” piece of a generical smooth boundary is here represented by q3 = ρ0, and q3 = ρ1.

Remark 2.2

It is worth noting that the slip boundary condition (2.1) is equivalent to


where τ stands for any arbitrary unit tangential vector on ∂Ω, and κτ is the principal curvature in the τ direction, positive if the center of curvature lies inside Ω.

The above claim follows immediately by appealing to equation (5.2) in [32], namely


For a mathematical treatment of some aspects related to slip boundary conditions imposed on smooth, but generic, boundaries see also [33], and the pioneering paper [34].

Next we recall some facts on curvilinear coordinates. The Lamé coefficients (scale factors) of the transition to the system of coordinates q are denoted by the letters Hi


Let e^i=1Hixqi, i = 1, 2, 3. Note that |êi| = 1 and êi=1Hiqi. One can write




It is well known (see for example [35, 36]) that




We state our main result as follows.

Theorem 2.2

Let Ω satisfy Assumption 2.1, and suppose that there exist two positive constants c and C such that


for any i, j, k = 1, 2, 3. Let u be a weak solution of the system (1.1) under the boundary condition (2.1), and set ū = û1ê1 + û2ê2. If ū satisfies


then the solution u is strong, namely,


Note that assumption (2.6) implies |iHj|, |ijHk| ≤ C.

It is worth noting that our proof applies to a more general set of geometrical situations. let’s just give some hint in this direction.

Remark 2.3

The above statement does not contain the result proved in reference [5], due to the distinct boundary conditions, see Section 4. On the other hand, we may replace the two circular, vertical, cylinders by more general vertical cylinders where the external circle q3 = ρ1 is replaced by a smooth Jordan curve y1, and the internal circle q3 = ρ0 by a parallel Jordan curve y0, at a sufficient small distance δ > 0 from y1. The coordinate θ is now an arc length coordinate on y1. All points in the same orthogonal segment to y0 and y1 enjoy the same θ coordinate. The coordinate r ∈ (0, δ) is given by the distance to y1. The “vertical” coordinate z preserves his periodic character. Clearly, the role played by the above Jordan curve may be immediately extended to much more general situations.

Another significant application is obtained by replacing the above two cylindrical boundaries by two torus of revolution, generated by revolving two concentric circles y0 and y1 about an axis coplanar with the circles, which does not touch the circles (roughly, we obtain the complement set between two closed tubes). Now z ∈ [0, 2 π) is an angular periodic coordinate, the toroidal coordinate. The result still applies by replacing the two circles by two parallel Jordan curves.

Let’s propose the following benchmark problem:

Problem 2.3

Consider two concentric spheres ΩR and Ωρ, of radius respectively ρ and R, 0 < ρ < R. Let u be a weak solution in ΩR × (0, T] of (1.1) under one of the above slip boundary conditions. Further, assume that the restricted P-S condition holds in (ΩR Ωρ) × (0, T] with respect to the tangential components of the velocity, and holds in Ωρ × (0, T] with respect to two arbitrary components of the velocity. Problem : To prove that u is a strong solution in ΩR × (0, T].

3 Proof of Theorem 2.2


We start by reducing the system (1.1) under the boundary condition (2.1) into the classical vorticity form


Then we take the scalar product with ω, and integrate by parts. One easily gets


Next, we focus on the estimates of I1 and I2.

Control of I1: First, it follows from (2.1) that


Let ∂Ωl = ∂Ω̂l := {(q1, q2, q3) ∈ Ω̂: q3 = ρl}, where l = 0, 1. One can deduce from (3.2) that


Since ∇ ⋅ ω = 0, from (2.4) one gets


which gives


since ω̂1 = ω̂2 = q1 ω̂1 = q2ω̂2 = 0 on ∂Ω̂l. Hence, one obtains


By appealing to (2.6) one shows that


where we have used Gagliardo’s trace theorem, see [37]. See also [38], Theorem 4.2 (for an English recent text see, for example, the Theorem III.2.21 in [39]). It follows that


for all 0 < ϵ < 1.

Control of I2: First, one has


For I21, from (2.6), one has


For I22, we consider separately the three cases j ≠ 3; j = 3 and i ≠ 3; i = j = 3.

  • Case I

    j ≠ 3. By integration by parts, one has


  • Case II

    j = 3 and i ≠ 3. From (2.5) one has


    Hence, by integration by parts, it follows that


  • Case III

    i = j = 3. Note that, due to ∇ ⋅ u = 0, it follows


One has


where the first equality is an identity, the second is obtained by appealing to (3.8), and the third one follows by integration by parts. From (3.6), (3.7), (3.9) and the assumption (2.6), one can obtain


It is easy to get that


and similarly to the proof of (3.5)


Hence, one has


By Hölder’s inequality, interpolation, and a Sobolev’s embedding theorem, one can easily show that


Collecting (3.1) and the estimates (3.4), (3.5), (3.10) and (3.11), one obtains


On the other hand, the following well known estimates (see for instance Theorem IV.4.8 and Theorem IV.4.9 in [39]), hold:


Therefore, from equation (3.12), by letting ϵ be sufficiently small, one has


Finally (1.2) follows by taking into account equations (2.7) (q2pp3 > 2) and (3.12), and by appealing to a well known argument, which is based on Gronwall’s inequality. Recall that weak solutions verify ∥uL2(Ω)L(0, T) and uL6(Ω)2L1(0, T). Hence we have proved that u is a strong solution. □

4 On related slip boundary conditions

In this section we present a first attempt to prove the statement of Theorem 2.2 with the slip boundary condition (2.1) replaced by the slip boundary condition (1.5) (assumed in reference [5]) by means of a simple modification of our proof. This attempt fails. Hence this significant problem remains open to further investigation. This leads us to briefly show our calculations.

Let’s start by explaining our guess. As still shown in Remark 2.2 condition (1.5) is equivalent to


We may replace the arbitrary tangent vector τ simply by a couple of independent vectors like, for instance, the principal direction’s vectors τ1 and τ2. In this case κ1 = κτ1 and κ2 = κτ2 are the maximum and the minimum principal curvatures.

A more natural choice here is to consider the couple of tangent, orthogonal, vectors ê1 and ê2. In this case κ1 and κ2 are the related curvatures. This second choice easily leads to the couple of linear equations


Hence to replace the slip boundary condition (2.1) by [D(u)n] ⋅ τ = 0 means to replace assumption (3.2) by


To prove our main statement with the boundary condition (2.1) replaced by the boundary condition (4.1) we have to control some new boundary integrals, which no longer vanish since now ω̂1 and ω̂2 do not vanish. However, by (4.2), ω̂1 and ω̂2 can be expressed in terms of the (lower order) velocity components û2 and û1. Well known inverse trace theorems allow us to control boundary-norms of these two components by suitable internal norms. Since our P-S assumption guarantees additional regularity just for these two velocity components, one could expect that the above internal norms could be estimated in a convenient way. Unfortunately this device seems not sufficient to prove our goal. So this interesting problem remains open.

Next we pass to showing our calculations. Let’s turn back to equation (3.3), by taking into account that now we can not apply to ω̂1 = ω̂2 = 0. One has


We will drop terms which could be easily manipulated, called here "lower order terms". Dropping lower order terms and also cancelling non significant multiplication coefficients, lead us to introduce the symbols ”≃” and ”⪯”, which have a clear meaning here.

One has


Since ∇ ⋅ ω = 0, from (2.4) one gets


which gives, on ∂Ω̂,


Under the new boundary conditions we can not apply to ω̂1 = ω̂2 = q1 ω̂1 = q2 ω̂2 = 0 on ∂Ω̂l to claim the cancellation of the above right hand side. By noting that the two terms on the right hand side are symmetric, with respect to the index 1 and 2, we may consider just the first one.

One has


Note that the smooth coefficients Hj, as their derivatives, are not significant on our estimate below. Further, since q1 is a tangential derivative, we may apply to the second equality (4.2) to assume that q1 ω̂1 ≃ –q1 û2 on ∂Ω̂. Hence


By appealing to Gagliardo’s theorem we show that the above right hand side is bounded by C(ϵ) ∥∇ u2 + ϵ ∥∇2 u2, which is sufficient to our purposes.

Let’s now consider in equation (4.4) the terms q3 (H1 H2 ω̂j) H31 ω̂j, for j = 1, 2. Assume, for instance, j = 1. One has q3(H1 H2 ω̂1) H31 ω̂1 ≃ (q3 ω̂1) û2. Hence we need to control the integral


Roughly speaking the above integrand has the same order as that on the right hand side of (4.5). However in (4.6) the derivation symbol q3 appears now in the “bad position". A suitable control of the above integral turns out to be not obvious.


The first author is partially supported by FCT (Portugal) under grant UID/MAT/04561/2019.


  • [1]

    J. Leray, Sur le mouvement d’un liquide visqueux emplissant l’espace, Acta Math. 63, (1934), 193–248. CrossrefGoogle Scholar

  • [2]

    E. Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen, Math. Nachr. 4, (1951), 213–231. Google Scholar

  • [3]

    A.A. Kiselev, O.A. Ladyzhenskaya, On the existence and uniqueness of the solution of the nonstationary problem for a viscous, incompressible fluid, Izv. Akad. Nauk SSSR Ser. Mat. 21, (1957), 655–680. Google Scholar

  • [4]

    J.L. Lions, Sur l’existence de solutions des équations de Navier-Stokes, C. R. Acad. Sci. Paris 248, (1959), 2847–2849. Google Scholar

  • [5]

    H. Beirão da Veiga, J. Bemelmans, J. Brand, On a two components condition for regularity of the 3D Navier-Stokes equations under physical slip boundary conditions on non-flat boundaries, Mathematische Annalen (2018), 1–38. Web of ScienceGoogle Scholar

  • [6]

    G. Prodi, Un teorema di unicità per le equazioni di Navier-Stokes, Ann. Mat. Pura Appl. 48, (1959), 173–182. CrossrefGoogle Scholar

  • [7]

    C. Foias, Une remarque sur l’unicité des solutions des équations de Navier-Stokes en dimension n, Bull. Soc. Math. Fr. 89, (1961), 1–8. Google Scholar

  • [8]

    J. Serrin, On the interior regularity of weak solutions of the Navier-Stokes equations, Arch. Ration. Mech. Anal. 9, (1962), 187–195. CrossrefGoogle Scholar

  • [9]

    J. Serrin, The initial value problem for the Navier-Stokes equations, In: Langer, R.E. (ed.) Nonlinear Problems, 69–98. University of Wisconsin Press, Madison, 1963. Google Scholar

  • [10]

    O.A. Ladyzhenskaya, On uniqueness and smoothness of generalized solutions to the Navier-Stokes equations, Zap. Nauchn. Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) 5, (1967), 169–185. Google Scholar

  • [11]

    O.A. Ladyzhenskaya, La théorie mathématique des fluides visqueux incompressibles, Moscou, 1961. English edition. 2nd edn. Gordon & Breach, New York, 1969. Google Scholar

  • [12]

    H. Sohr, Zur Regularitätstheorie der instationären Gleichungen von Navier-Stokes, Math. Z. 184, (1983), 359–375. CrossrefGoogle Scholar

  • [13]

    W. von Wahl, Regularity of weak solutions of the Navier-Stokes equations, Proc. Sympos. Pure Math. 45, (1986), 497–503. Google Scholar

  • [14]

    Y. Giga, Solutions for Semilinear Parabolic Equations in Lp and Regularity of Weak Solutions of the Navier–Stokes System, J. Differential Equations 62, (1986), 186–212. CrossrefGoogle Scholar

  • [15]

    G.P. Galdi, P. Maremonti, Sulla regolarità delle soluzioni deboli al sistema di Navier-Stokes in domini arbitrari, Ann. Univ. Ferrara Sez. VII Sci. Mat. 34, (1988), 59–73. Google Scholar

  • [16]

    G.P. Galdi, An Introduction to the Navier-Stokes Initial-Boundary Value Problems, In: G.P. Galdi, M.I. Heywood, R. Rannacher, (eds.) Fundamental Directions in Mathematical Fluid Mechanics. Advances in Mathematical Fluid Mechanics, pp. 1–70, Birkhäuser, Basel, 2000. Google Scholar

  • [17]

    G. Prodi, Résultats récents et problèmes anciens dans la théorie des équations de Navier-Stokes, In: Les Équations aux Dérivées Partielles. Colloques Intern. du CNRS 117, pp. 181–196, Paris, 1962. Google Scholar

  • [18]

    L. Escauriaza, G. Seregin, V. Šverák, L3,∞-Solutions to the Navier-Stokes Equations and Backward Uniqueness, Russian Math. Surveys 58, (2003), 211–250. CrossrefGoogle Scholar

  • [19]

    H.-O. Bae, H.J. Choe, L-bound of weak solutions to Navier-Stokes equations, In: Proceedings of the Korea-Japan Partial Differential Equations Conference (Taejon, 1996). Lecture Notes Ser. 39. Seoul Nat. Univ., Seoul, p. 13, 1997. Google Scholar

  • [20]

    H.-O. Bae, H.J. Choe, A regularity criterion for the Navier-Stokes equations, Commun. Partial Differ. Equations 32, (2007), 1173–1187. CrossrefWeb of ScienceGoogle Scholar

  • [21]

    H. Beirão da Veiga, On the extension to slip boundary conditions of a Bae and Choe regularity criterion for the Navier-Stokes equations. The half space case, J. Math. Anal. Appl. 453, (2017), 212–220. CrossrefWeb of ScienceGoogle Scholar

  • [22]

    H.-O. Bae, J. Wolf, A local regularity condition involving two velocity components of Serrin-type for the Navier–Stokes equations, C. R. Acad. Sci. Paris, Ser. I. 354, (2016), 167–174. CrossrefGoogle Scholar

  • [23]

    H. Beirão da Veiga, On the Smoothness of a Class of Weak Solutions to the Navier–Stokes equations, J. Math. Fluid Mech. 2, (2000), 315–323. CrossrefGoogle Scholar

  • [24]

    L.C. Berselli, A note on regularity of weak solutions of the Navier-Stokes equations in ℝn, Jpn. J. Math. 28, (2002), 51–60. CrossrefGoogle Scholar

  • [25]

    D. Chae, H.-J. Choe, Regularity of Solutions to the Navier-Stokes Equation, Electron. J. Differential Equations 05, (1999), 1–7. Google Scholar

  • [26]

    C. Cao, E.S. Titi, Regularity Criteria for the Three-dimensional Navier–Stokes Equations, Indiana Univ. Math. J. 57, (2008), 2643–2661. CrossrefWeb of ScienceGoogle Scholar

  • [27]

    C. He, Regularity for solutions to the Navier-Stokes equations with one velocity component regular, Electron. J. Differential Equations 29, (2002), 1–13. Google Scholar

  • [28]

    I. Kukavica, M. Ziane, Navier-Stokes equations with regularity in one direction, J. Math. Phys. 48, (2007), 2643–2661. Web of ScienceGoogle Scholar

  • [29]

    J. Neustupa, P. Penel, Anisotropic and Geometric Criteria for Interior Regularity of Weak Solutions to the 3D Navier–Stokes Equations, In: J. Neustupa, P. Penel, (eds.) Mathematical Fluid Mechanics. Advances in Mathematical Fluid Mechanics, pp. 237–265. Birkhäuser, Basel, 2001. Google Scholar

  • [30]

    Z. Zhang, D. Zhong, X. Huang, A refined regularity criterion for the Navier-Stokes equations involving one non-diagonal entry of the velocity gradient, J. Math. Anal. Appl. 453, (2017), 1145–1150. CrossrefWeb of ScienceGoogle Scholar

  • [31]

    Y. Zhou, M. Pokorný, On the regularity of the solutions of the Navier–Stokes equations via one velocity component, Nonlinearity 23, (2010), 1097–1107. Web of ScienceCrossrefGoogle Scholar

  • [32]

    H. Beirão da Veiga, F. Crispo, Concerning the Wk,p-inviscid limit for 3-D flows under a slip boundary condition, Journal of Mathematical Fluid Mechanics 13, (2011), no. 1, 117–135. CrossrefWeb of ScienceGoogle Scholar

  • [33]

    H. Beirão da Veiga, Regularity for Stokes and generalized Stokes systems under nonhomogeneous slip-type boundary conditions, Adv. Differential Equations 9, (2004), 1079–1114. Google Scholar

  • [34]

    V.A. Solonnikov, V.E. Ščadilov, On a boundary value problem for a stationary system of Navier-Stokes equations, Proc. Steklov Inst. Math. 125, (1973), 186–199. Google Scholar

  • [35]

    G.K. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, Cambridge, 1967. Google Scholar

  • [36]

    N.E. Kochin, A.K. Il’ja, N.V. Roze, Theoretical hydromechanics, Interscience, 1964. Google Scholar

  • [37]

    E. Gagliardo, Caratterizzazioni delle tracce sulla frontiera relative ad alcune classi di funzioni in n variabili, Rend. Sem. Mat. Univ. Padova 27, (1957), 284–305. Google Scholar

  • [38]

    J. Nečas, Les Méthodes Directes en Théorie des Équations Elliptiques, Academia, Prague, 1967. Google Scholar

  • [39]

    F. Boyer, P.Fabrie, Mathematical tools for the study of the incompressible Navier-Stokes equations and related models, Springer Science & Business Media, 2012. Google Scholar

About the article

Received: 2018-12-24

Accepted: 2019-01-18

Published Online: 2019-06-29

Published in Print: 2019-03-01

Citation Information: Advances in Nonlinear Analysis, Volume 9, Issue 1, Pages 633–643, ISSN (Online) 2191-950X, ISSN (Print) 2191-9496, DOI: https://doi.org/10.1515/anona-2020-0017.

Export Citation

© 2020 Hugo Beirão da Veiga and Jiaqi Yang, published by De Gruyter. This work is licensed under the Creative Commons Attribution 4.0 Public License. BY 4.0

Comments (0)

Please log in or register to comment.
Log in