In this work, we are not trying to consider the extensive references on the Navier-Stokes system, since there are fundamental works in this area [2,3,4] and others. Some special of the above problems were also investigaqted in the works [5,6,7,8] etc.

It is known that the methods of integral transformations in the theory of partial differential equations made it possible to find solutions to many problems [11,12] and clarify the physical meaning of some basic laws and phenomena in fluid mechanics. Therefore, this paper presents one of the developed transformations of the said species.

In this regard, in the present work, we study the Navier-Stokes system, which describe the flow of a viscous incompressible fluid filling all of $R^3$, i.e.: with initial conditions where $R^3∍\psi (x)$is the known initial velocity vector, $R^3∍f(x,t)$is external applied force (e.g. gravity), $0<\mu$is kinematic viscosity, $\rho$is density,$\Delta$is Laplace operator,$\nabla$is Hamilton operator. These equations are to be solved for an unknown velocity vector $\mathsf{\nu }\in {\mathrm{R}}^3$and pressure$P(x,t)$, and equation (1.2) just says that the fluid is incompressible.

Aim of research. The main object of this work is to establish the existence and uniqueness of a solution to the Navier-Stokes system for an incompressible fluid and at the same time, it is proved that:

a) the solutions of the transformed equations are regular with respect to the viscosity coefficient µ, and they simplify the analysis of the original problem,

b) the found pressure distribution law, which is described by a Poisson type equation and plays a fundamental role in the theory of Navier-Stokes systems in constructing analytic conditionally smooth (smooth) solutions.

c) at the same time, the obtained results meet the requirements of the "Navier-Stokes Millennium problem-NSMP " (2000y.).

In the introduced space $G_{3,h}^1(D_0)$, the norm is defined as: where $k=(k_1,k_2,k_3)$ is the multi-index,

A) In this regard, we note that in our early works, for example, in [9], we proposed method for constructing smooth (conditionally smooth) solutions of 3D Navier-Stokes equations in $G_3^1(D_1=R^3\times (0,T_0))$: with the condition:where $\phi (x)$is known scalar function, $R^3∍\lambda$is given vector with positive constant components: $0<{\lambda }_i,(i=\overline{1,3})$. Since it takes place then, we seek the solution of the Navier-Stokes problem in the form:here $R^3\in J$ is known vector-valued function of the form:and $\theta (x,t)$is a new unknown scalar function with the condition:

And since then, on the basis of (1.7), (1.10) and (1.11), the inertial terms of equation (1.1) are equivalently converted to the form:So this means that under condition (1.2), the inertial terms of equation (1.1), taking into account (1.7), are linearized with respect to the newly introduced function $\theta (x,t)$and its derivatives with respect to$x\in R^3$, and the nonlinearity goes over to the known vector of the function $J(x,t)$ and partial derivatives with respect to$x\in R^3$. Which was required to show.

Further, substituting (1.7) into equation (1.1), we obtain a linear inhomogeneous differential equation of the type of heat conduction with variable coefficients:in this case, (1.14) is the condition of unequivocal compatibility for case (1.13), since $\theta$is a scalar function.

a) In transformation (1.7) it is assumeds that: $divf=0$. But this transformation can also be introduced in the case when: $divf\ne 0$. For this purpose, let us introduce the vector function$f^0(x,t),(x\in R^3)$:Then in this case, the vector function $J(x,t)$is represented as:Next, we obtain similar results as in the case of lemma 1.

b) To understand (1.14) we give a concrete example from the field of the system of algebraic equations. For this purpose, consider a system with one unknown quantity z, i.e.:From here we see that if there is then z is uniquely defined in the form But z, can be defined differently, i.e.:or with respect to (0.6), performing some mathematical transformation we obtain This means that the first and second paths are equivalent. So, under condition (0.4), z is indeed uniquely determined from (0.3). Which was reguired to show.

Note that (1.14) has conditions of the type (0.4) for the system with scalar unknown. Therefore, when we solve the system (1.13) we choose the second path as shown in the case (0.6).

Next, the equation for the pressure is derived:We are taking into account the operation div with respect to (1.13), (that's tantamount to applying the operation div with respect to equation (1.1), since (1.1) is equivalently converted to the form (1.13) based on (1.7)), since takes places:Here, formula (1.15) modifies the Landau – Lifshitz formula (see [2]: (15.11)) and is an equation of Poisson type. Then it follows from (1.15), at that where (1.16) is called the Newtonian potential [11]. On the other hand, a solution to the Poisson equation (1.15) tending to zero at infinity will be unique if the function${\theta }_{x_i},(i=1,2,3)$is unique, since the function$\Omega (x,t)$ contains these functions.

To prove the above, we note that the obtained pressure distribution law allows us to express the velocity in integral form when $\mathsf{\nu }\in R^3$. In fact, substituting (1.17) into equation (1.13) with allowance for (1.14), we obtain an inhomogeneous linear integro-differential heat conduction equation with the Cauchy condition:here the known functions contained in (1.18) are introduced on the basis of the notation:

As a result, problem (1.18) is transformed to a system of Volterra and Volterra-Abel equations of second kind, where the solvability of this problem in $G_{}^1(D_1)$ follows from the solvability of this system. Therefore, we obtain similar conclusions for problem (1.1)-(1.3) in $G_3^1(D_1)$.

B) Similar issues were investigated in [10], i.e., equation (1.1), (1.2) with the condition:at that$f_i(x,t)$ is the component of a given external force $f$ admits the conditions:and this means that $R^3∍\lambda$ is a vector with constant components: $0<{\lambda }_i,(i=1,2,3)$, therefore, it becomes possible to use modification of the method (1.7) of the previous sector, i.e.:where (1.8) is taken into account. At that $\theta (x,t)$is a new unknown scalar function with the condition:Here (1.22) differs from (1.7), since $t\in R_{+}$, therefore as a multiplier function, we introduced: ${\Omega }_0(t)\equiv {(1+t)}^{-q}$. Hence, follow (1.10) and (1.11) we have And Then, taking into account (1.20), (1.24) and (1.25), the inertial terms of equation (1.1) are equivalently converted to the form:Constraints on external force of the form (1.21) make it possible to simplify the Navier-Stokes problem and transform it into a system of integral equations of the second kind. In fact, on the basis of (1.22) and (1.26), from (1.1) follows the equation:since $\theta$ is a scalar function, then the condition:is a univocal compatibility condition for (1.28). From where the equation for pressure is derived:where On the other hand, we note that it follows from (1.29):here (1.30) tends to zero at infinity, and there are second-order partial continuous derivatives, and for the first-order partial derivatives it takes place:Therefore, excluding pressure from (1.27), we obtain a linear differential equation with variable coefficients and with the Cauchy condition of the form:where In such an approach the solution of problem is reduced to finding two functions, $\theta (x,t)$and $\zeta (x,t).$ The latter can be usually determined without difficulties and we will solve the equation with respect to $\zeta (x,t)$ by the Picard’s method [11]. Besides, since the function $\theta$ has continuous partial derivatives up to the second order inclusive with respect to spatial coordinates, and a first order time derivative, then problem (1.32) with sufficiently smooth data is solvable in $W^0(D)$. Therefore, based on (1.22) we have similar conclusions for problem (1.1)-(1.3) we obtain in $W_3^0(D)$:

Methods of analysis of physical phenomena are based on statements of corresponding mathematical problems formulated by means of various kinds of functional equations and certain additional conditions. The solutions of these problems may be considered the main aim of the theoretical investigation.

Let is the velocity vector satisfies conditions (1.2), (1.3) and takes place where $R^3∍\lambda$ is a known vector with positive constant components: $0<{\lambda }_i,(i=\overline{1,3})$, and then applying the transformation:where $\theta (x,t)$ is the new unknown scalar function, and $0<{\delta }_0$ is the introduced constant, which ensures the application of the Banach principle and the Picard's method for the system of integral equations of Volterra-Abel type of the second kind, into which the original problem is transformed, $R^3∍J(x,t)$ is the given vector:where $\theta (x,t)$ is a new unknown scalar function with the condition:at that Then, taking into account (1.11) and (2.5), the inertial terms of equation (1.1) are equivalently converted to the form:The conditions of the form (2.6) make it possible to simplify the Navier-Stokes problem and transform it into a system of Volterian type integral equations of the second kind. In fact, on the basis of (2.2) and (2.6), from (1.1) follows the equation:since $\theta$ is a scalar function, then the condition:is a univocal compatibility condition for (2.8). In addition, the Poisson equation for pressure is derived in the form:and it is obtained on the basis of (2.2) by applying the operation $\mathrm{div}$ to equation (1.1), (or (2.7)), since On the other hand, we note that it follows from (2.9):here (2.10) tends to zero at infinity, and there are second-order partial continuous derivatives, and for the first-order partial derivatives it takes place:Therefore, excluding pressure from (3.3), we obtain a linear differential equation with variable coefficients and with the Cauchy condition of the form:where Further, the solution of the problem under study is reduced to the determination of functions from the equations: So, taking into account (2.12), (2.14) and $\zeta (x,t)$ we have at that Similar assessments can be made regarding ${\gamma }_1$ in $G_h^1(D_0)$, therefore we have ${\gamma }_1+{\gamma }_2=\gamma \in G_{}^1(D_0).$Since, operator ${\Gamma }_0$ of system (2.15) contains small viscositie $\delta$and $\sqrt{{\delta }_0}$, then the proof of solvability and the construction of the approximate solution can be realized on the basis of the Banach principle and the Picard's method.