\documentclass[reqno]{amsart} \usepackage{hyperref} \AtBeginDocument{{\noindent\small \emph{Electronic Journal of Differential Equations}, Vol. 2013 (2013), No. 83, pp. 1--17.\newline ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu \newline ftp ejde.math.txstate.edu} \thanks{\copyright 2013 Texas State University - San Marcos.} \vspace{9mm}} \begin{document} \title[\hfilneg EJDE-2013/83\hfil 3D Navier-Stokes equations] {Existence, uniqueness and smoothness of a solution for 3D Navier-Stokes equations with any smooth initial velocity} \author[A. Tsionskiy, M. Tsionskiy \hfil EJDE-2013/83\hfilneg] {Arkadiy Tsionskiy, Mikhail Tsionskiy} % in alphabetical order \address{Arkadiy Tsionskiy, Mikhail Tsionskiy\newline Department of Civil, Environmental and Ocean Engineering, Stevens Institute of Technology, Hoboken, NJ 07030, USA} \email{amtsionsk@gmail.com} \thanks{Submitted December 10, 2012. Published April 5, 2013.} \subjclass[2000]{35Q30, 76D05} \keywords{3D Navier-Stokes equations; Fourier transform; Laplace transform; \hfill\break\indent fixed point principle} \begin{abstract} Solutions of the Navier-Stokes and Euler equations with initial conditions for 2D and 3D cases were obtained in the form of converging series, by an analytical iterative method using Fourier and Laplace transforms in \cite{TT10,TT11}. There the solutions are infinitely differentiable functions, and for several combinations of parameters numerical results are presented. This article provides a detailed proof of the existence, uniqueness and smoothness of the solution of the Cauchy problem for the 3D Navier-Stokes equations with any smooth initial velocity. When the viscosity tends to zero, this proof applies also to the Euler equations. \end{abstract} \maketitle \numberwithin{equation}{section} \newtheorem{theorem}{Theorem}[section] \newtheorem{remark}[theorem]{Remark} \allowdisplaybreaks \section{Introduction} Many authors have obtained results regarding the Euler and Navier-Stokes equations. Existence and smoothness of solution for the Navier-Stokes equations in two dimensions have been known for a long time. Leray (1934) showed that the Navier-Stokes equations in three dimensional space have a weak solution. Scheffer (1976, 1993) and Shnirelman (1997) obtained weak solution of the Euler equations with compact support in space-time. Caffarelli, Kohn and Nirenberg (1982) improved Scheffer's results, and Lin (1998) simplified the proof of the results by Leray. Many problems and conjectures about behavior of weak solutions of the Euler and Navier-Stokes equations are described in the books by Ladyzhenskaya (1969), Temam (1977), Constantin (2001), Bertozzi and Majda (2002) or Lemari\'e-Rieusset (2002). The solution of the Cauchy problem for the 3D Navier-Stokes equations is described in this article. We will consider an initial velocity that is infinitely differentiable and decreasing rapidly to zero in infinity. The applied force is assumed to be identically zero. A solution of the problem will be presented in the following stages: \textbf{First stage} (sections 2, 3). We move the non-linear parts of equations to the right side. Then we solve the system of linear partial differential equations with constant coefficients. We have obtained the solution of this system using Fourier transforms for the space coordinates and Laplace transform for time. From theorems about applications of Fourier and Laplace transforms, for system of linear partial differential equations with constant coefficients, we see that in this case if initial velocity and applied force are smooth enough functions decreasing in infinity, then the solution of such system is also a smooth function. Corresponding theorems are presented in Bochner \cite{SB59}, Palamodov \cite{VP70}, Shilov \cite{gS01}, Hormander \cite{LH83}, Mizohata \cite{SM73}, Treves \cite{JFT61}. The result of this stage is an integral equation for the vector-function of velocity. \textbf{Second stage} (sections 4, 5). We introduce perfect spaces of functions and vector-functions (Gel'fand, Chilov \cite{GC68}), in which we look for the solution of the problem. We demonstrate the equivalence of solving the Cauchy problem in differential form and in the form of an integral equation. \textbf{Third stage} (section 6). We divide all parts of the integral equation by an appropriate constant depending on value of initial fluid velocity, and obtain the equivalent integral equation. We also replace the corresponding integration variables in the integral operators. This newly received equivalent integral equation allowed us to analyze the Cauchy problem for the 3D Navier-Stokes equations for any value of initial fluid velocity. According to a priori estimate of the solution of the Cauchy problem for the 3D Navier-Stokes equations \cite{oL69,LK63}, the described constant is proportional to max of the norms of the initial velocity in the spaces $C^2$ and L$_2$. \textbf{Fourth stage} (section 6). We use the newly obtained equivalent integral equation to prove the existence and uniqueness of the solution of the Cauchy problem in the time range $[0,\infty)$ based on the Caccioppoli-Banach fixed point theorem (Kantorovich, Akilov \cite{KA64}, Trenogin \cite{VT80}, Rudin \cite{WR73}, Kirk and Sims \cite{KS01}, Granas and Dugundji \cite{GD03}, Ayerbe Toledano, Dominguez Benavides, Lopez Acedo \cite{ADL97}). For this purpose the following three theorems are proven in this article: Theorem \ref{thm1}: the integral operator of the problem is a contraction operator; Theorem \ref{thm2}: the existence and uniqueness of the solution of the problem is valid for any $t\in [0,\infty)$; Theorem \ref{thm3}: the solution of the problem depends continuously on $t$. \textbf{Fifth stage} (section 6). By using a priori estimate of the solution of the Cauchy problem for the 3D Navier-Stokes equations \cite{oL69,LK63}, we show that the energy of the whole process has a finite value for any $t$ in $[0,\infty)$. \section{Mathematical setup} The Navier-Stokes equations describe the motion of a fluid in $\mathbb{R}^N$ ($N = 3$). We look for a viscous incompressible fluid filling all of $\mathbb{R}^N$ here. The Navier-Stokes equations are then given by \begin{gather}\label{eqn1} \frac{\partial u_k}{\partial t} + \sum_{n=1}^N u_n \frac{\partial u_k}{\partial x_n} = \nu\Delta u_k - \frac{\partial p}{\partial x_k} +f_k(x,t)\quad x\in \mathbb{R}^N,\;t\geq 0,\; 1\leq k \leq N\,, \\ \label{eqn2} \operatorname{div}\vec{u}= \sum_{n=1}^N \frac{\partial u_n}{\partial x_n} = 0\quad x\in \mathbb{R}^N,\;t\geq 0\,, \end{gather} with initial conditions $$\label{eqn3} \vec{u}(x,0)= \vec{u}^0(x)\quad x\in \mathbb{R}^N\,.$$ Here $\vec{u}(x,t)=(u_k(x,t)) \in \mathbb{R}^N$ $(1\leq k \leq N)$ is an unknown velocity vector, $N = 3$; $p(x,t)$ is an unknown pressure; $\vec{u}^0(x)$ is a given $C^{\infty}$ divergence-free vector field; $f_k(x,t)$ are components of a given, externally applied force $\vec{f}(x,t)$; $\nu$ is a positive coefficient of the viscosity (if $\nu = 0$ then \eqref{eqn1}--\eqref{eqn3} are the Euler equations); and $\Delta= \sum_{n=1}^N \frac{\partial^2}{\partial x_n^2}$ is the Laplacian in the space variables. Equation \eqref{eqn1} is Newton's law for a fluid element. Equation \eqref{eqn2} says that the fluid is incompressible. For physically reasonable solutions, we accept $$\label{eqn4} u_k(x,t) \to 0, \quad \frac{\partial u_k}{\partial x_n} \to 0\quad \text{as }| x | \to \infty\quad 1\leq k \leq N ,\; 1\leq n \leq N\,.$$ Hence, we will restrict attention to initial conditions $\vec{u}^0$ and force $\vec{f}$ that satisfy $$\label{eqn5} |\partial_{x}^{\alpha}\vec{u}^0(x)|\leq C_{\alpha K}(1+| x |)^{-K} \quad \text{on \mathbb{R}^N for any \alpha and any K}.$$ and $$\label{eqn6} |\partial_{x}^{\alpha}\partial_{t}^{\beta}\vec{f}(x,t)| \leq C_{\alpha \beta K}(1+| x | +t)^{-K} \quad \text{on \mathbb{R}^N\times[0,\infty) for any \alpha,\beta and any K}.$$ To start the process of solution let us add $- \sum_{n=1}^N u_n\frac{\partial u_k}{\partial x_n}$ to both sides of the equations \eqref{eqn1}. Then we have \begin{gather}\label{eqn7} \frac{\partial u_k}{\partial t}=\nu\,\Delta\,u_k -\frac{\partial p}{\partial x_k}+f_k(x,t)-\sum_{n=1}^N u_n\frac{\partial u_k}{\partial x_n}\quad x\in \mathbb{R}^N,\;t\geq 0,\; 1\leq k \leq N,\\ \label{eqn8} \operatorname{div}\vec{u}= \sum_{n=1}^N \frac{\partial u_n}{\partial x_n} = 0 \quad x\in \mathbb{R}^N,\;t\geq 0,\\ \label{eqn9} \vec{u}(x,0)= \vec{u}^0(x)\quad x\in \mathbb{R}^N,\\ \label{eqn10} u_k(x,t) \to 0\quad \frac{\partial u_k}{\partial x_n}\to 0\quad\text{as } | x | \to \infty\quad 1\leq k \leq N,\; 1\leq n \leq N,\\ \label{eqn11} |\partial_{x}^{\alpha}\vec{u}^0(x)|\leq C_{\alpha K}(1+| x |)^{-K} \quad \text{on $\mathbb{R}^N$ for any $\alpha$ and any $K$},\\ \label{eqn12} |\partial_{x}^{\alpha}\partial_{t}^{\beta}\vec{f}(x,t)| \leq C_{\alpha \beta K}(1+| x | +t)^{-K} \quad \text{on $\mathbb{R}^N\times[0,\infty)$ for any $\alpha,\beta$ and any $K$}. \end{gather} Let us denote $$\label{eqn13} \tilde{f}_k(x,t)= f_k(x,t) - \sum_{n=1}^N u_n\frac{\partial u_k}{\partial x_n}\quad 1\leq k \leq N\,.$$ We can present it in the vector form as $$\label{eqn14} \vec{\tilde{f}}(x,t)= \vec{f}(x,t) -(\vec{u}\cdot \nabla)\vec{u}\,.$$ \section{Solution of the system \eqref{eqn7}--\eqref{eqn14}} Let us assume that all operations below are valid. The validity of these operations will be proved in the next sections. Taking into account our substitution \eqref{eqn13} we see that \eqref{eqn7}--\eqref{eqn9} are in fact system of linear partial differential equations with constant coefficients. The solution of this system will be presented by the following steps: \textbf{First step.} We use Fourier transform \eqref{A3} to solve equations \eqref{eqn7}--\eqref{eqn14}. We obtain: \begin{gather*} U_k( \gamma_1, \gamma_2, \gamma_3, t)=F[u_k ( x_1, x_2, x_3, t)],\\ -\gamma^2_{s}U_k( \gamma_1, \gamma_2, \gamma_3, t) =F[ \frac{\partial^2u_k( x_1, x_2, x_3, t)}{\partial x^2_{s}}] \quad \text{(use \eqref{eqn10})},\\ U_k^0( \gamma_1 ,\gamma_2 ,\gamma_3)=F[u_k^0 ( x_1, x_2, x_3)],\\ P( \gamma_1, \gamma_2, \gamma_3, t)=F[p ( x_1, x_2, x_3, t)],\\ \tilde{F}_k( \gamma_1, \gamma_2, \gamma_3, t) =F[\tilde{f}_k ( x_1, x_2, x_3, t)], \end{gather*} for $k,s=1,2,3$. Then \begin{gather}\label{eqn134} \begin{aligned} \frac{d U_1( \gamma_1, \gamma_2, \gamma_3, t )}{d t} & =-\nu ( \gamma_1^2 +\gamma_2^2 +\gamma_3^2) U_1( \gamma_1, \gamma_2, \gamma_3, t )+i\gamma_1 P( \gamma_1, \gamma_2, \gamma_3, t )\\ &\quad + \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3,t ), \end{aligned} \\ \label{eqn135} \begin{aligned} \frac{d U_2( \gamma_1, \gamma_2, \gamma_3, t )}{d t} &=-\nu ( \gamma_1^2 +\gamma_2^2 +\gamma_3^2) U_2( \gamma_1, \gamma_2, \gamma_3, t )+i\gamma_2 P( \gamma_1, \gamma_2, \gamma_3, t )\\ &\quad + \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3,t ), \end{aligned} \\ \label{eqn136} \begin{aligned} \frac{d U_3( \gamma_1, \gamma_2, \gamma_3, t )}{d t} &=-\nu ( \gamma_1^2 +\gamma_2^2 +\gamma_3^2) U_3( \gamma_1, \gamma_2, \gamma_3, t ) +i\gamma_3 P( \gamma_1, \gamma_2, \gamma_3, t )\\ &\quad + \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3,t ), \end{aligned}\\ \label{eqn137} \gamma_1 U_1( \gamma_1, \gamma_2, \gamma_3, t ) + \gamma_2\, U_2( \gamma_1, \gamma_2, \gamma_3, t ) + \gamma_3\, U_3( \gamma_1, \gamma_2, \gamma_3, t ) =0,\\ \label{eqn138} U_1(\gamma_1, \gamma_2, \gamma_3, 0)= U_1^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn139} U_2(\gamma_1, \gamma_2, \gamma_3, 0)= U_2^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn140} U_3(\gamma_1, \gamma_2, \gamma_3, 0)= U_3^0(\gamma_1 ,\gamma_2 ,\gamma_3)\,. \end{gather} Hence, we have received a system of linear ordinary differential equations with constant coefficients \eqref{eqn134}-\eqref{eqn140} according to Fourier transforms. At the same time the initial conditions are set only for Fourier transforms of velocity components $U_1( \gamma_1, \gamma_2, \gamma_3, t )$, $U_2( \gamma_1, \gamma_2, \gamma_3, t )$, $U_3( \gamma_1, \gamma_2, \gamma_3, t )$. Because of that we can eliminate Fourier transform for pressure $P( \gamma_1, \gamma_2, \gamma_3, t )$ from equations \eqref{eqn134}--\eqref{eqn136} on the next step of the solution process. \textbf{Second step.} From here assuming that $\gamma_1 \neq 0, \gamma_2 \neq 0,\gamma_3 \neq 0$, we eliminate $P( \gamma_1, \gamma_2, \gamma_3, t )$ from equations $\eqref{eqn134}- \eqref{eqn136}$ and find \begin{gather}\label{eqn141} \begin{aligned} &\frac{d}{dt} [ U_2( \gamma_1, \gamma_2, \gamma_3, t ) - \frac{\gamma_2}{\gamma_1} \,U_1( \gamma_1, \gamma_2, \gamma_3, t) ]\\ &= -\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2) [ U_2( \gamma_1, \gamma_2, \gamma_3, t ) - \frac{\gamma_2}{\gamma_1} U_1( \gamma_1, \gamma_2, \gamma_3, t) ]\\ &\quad + [ \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3, t ) -\frac{\gamma_2}{\gamma_1} \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, t) ]\,, \end{aligned}\\ \label{eqn142} \begin{aligned} &\frac{d}{dt} [ U_3( \gamma_1, \gamma_2, \gamma_3, t ) -\frac{\gamma_3}{\gamma_1}\, U_1( \gamma_1, \gamma_2, \gamma_3, t) ] \\ &= -\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2) [ U_3( \gamma_1, \gamma_2, \gamma_3, t ) - \frac{\gamma_3}{\gamma_1}\, U_1( \gamma_1, \gamma_2, \gamma_3, t) ]\\ &\quad + [ \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, t ) -\frac{\gamma_3}{\gamma_1} \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, t) ]\,, \end{aligned}\\ \label{eqn143} \gamma_1 U_1( \gamma_1, \gamma_2, \gamma_3, t ) + \gamma_2\, U_2( \gamma_1, \gamma_2, \gamma_3, t ) + \gamma_3\, U_3( \gamma_1, \gamma_2, \gamma_3, t ) =0,\\ \label{eqn144} U_1(\gamma_1, \gamma_2, \gamma_3, 0)= U_1^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn145} U_2(\gamma_1, \gamma_2, \gamma_3, 0)= U_2^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn146} U_3(\gamma_1, \gamma_2, \gamma_3, 0)= U_3^0(\gamma_1 ,\gamma_2 ,\gamma_3)\,. \end{gather} \textbf{Third step.} We use Laplace transform \eqref{A4}, \eqref{A5} for a system of linear ordinary differential equations with constant coefficients \eqref{eqn141}--\eqref{eqn143} and have as a result the system of linear algebraic equations with constant coefficients: \begin{gather} U_k^{\otimes} (\gamma_1, \gamma_2, \gamma_3, \eta) =L[U_k(\gamma_1, \gamma_2, \gamma_3, t)] \quad k=1,2,3;\\ \tilde{F}_k^{\otimes} (\gamma_1, \gamma_2, \gamma_3, \eta) =L[\tilde{F}_k(\gamma_1, \gamma_2, \gamma_3, t)] \quad k=1,2,3; \\ \label{eqn147} \begin{aligned} &\eta [ U_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_2}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) ]\\ &\quad - [ U_2( \gamma_1, \gamma_2, \gamma_3, 0 ) -\frac{\gamma_2}{\gamma_1} U_1( \gamma_1, \gamma_2, \gamma_3, 0) ] \\ &= -\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)[ U_2^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_2}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) ]\\ &\quad + [ \tilde{F}_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_2}{\gamma_1} \tilde{F}_1^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) ], \end{aligned}\\ \label{eqn148} \begin{aligned} &\eta [ U_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_3}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) ]\\ &\quad - [ U_3( \gamma_1, \gamma_2, \gamma_3, 0 ) -\frac{\gamma_3}{\gamma_1} U_1( \gamma_1, \gamma_2, \gamma_3, 0) ] \\ &= -\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)[ U_3^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_3}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) ]\\ &\quad + [ \tilde{F}_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) -\frac{\gamma_3}{\gamma_1} \tilde{F}_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) ], \end{aligned} \\ \label{eqn149} \gamma_1 U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) + \gamma_2\, U_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) + \gamma_3\, U_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) =0,\\ \label{eqn150} U_1(\gamma_1, \gamma_2, \gamma_3, 0)= U_1^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn151} U_2(\gamma_1, \gamma_2, \gamma_3, 0)= U_2^0(\gamma_1 ,\gamma_2 ,\gamma_3),\\ \label{eqn152} U_3(\gamma_1, \gamma_2, \gamma_3, 0)= U_3^0(\gamma_1 ,\gamma_2 ,\gamma_3)\,. \end{gather} Let us rewrite system of equations \eqref{eqn147}--\eqref{eqn149} in the form \begin{gather}\label{eqn147a} \begin{aligned} &[\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] \frac{\gamma_2}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta)\\ &\quad - [\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] U_2^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) \\ &= [ \frac{\gamma_2}{\gamma_1} \tilde{F}_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) -\tilde{F}_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) ]\\ &\quad + [ \frac{\gamma_2}{\gamma_1} U_1( \gamma_1, \gamma_2, \gamma_3, 0) -U_2( \gamma_1, \gamma_2, \gamma_3, 0 ) ], \end{aligned} \\ \label{eqn148a} \begin{aligned} &[\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] \frac{\gamma_3}{\gamma_1} U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta)\\ &\quad - [\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] U_3^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) \\ &=[ \frac{\gamma_3}{\gamma_1} \tilde{F}_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) -\tilde{F}_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) ]\\ &\quad + [ \frac{\gamma_3}{\gamma_1} U_1( \gamma_1, \gamma_2, \gamma_3, 0) -U_3( \gamma_1, \gamma_2, \gamma_3, 0 ) ], \end{aligned}\\ \label{eqn149a} \gamma_1 U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) + \gamma_2\, U_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) + \gamma_3\, U_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta ) =0 \end{gather} The determinant of this system is \label{eqn1481b} \begin{aligned} \Delta &= \begin{vmatrix} [\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)]\frac{\gamma_2}{\gamma_1} & -[\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] & 0 \\ [\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)]\frac{\gamma_3}{\gamma_1} & 0 & -[\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)] \\ \gamma_1 & \gamma_2 & \gamma_3 \end{vmatrix} \\ & = \frac{[\eta +\nu( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)]^2( \gamma_1^2 +\gamma_2^2 +\gamma_3^2)}{\gamma_1} \neq 0\,. \end{aligned} Consequently the system of equations \eqref{eqn147}--\eqref{eqn149} and/or \eqref{eqn147a}--\eqref{eqn149a} has a unique solution. Taking into account formulas \eqref{eqn150}--\eqref{eqn152} we can write this solution in the form \begin{gather}\label{eqn153} \begin{aligned} &U_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta )\\ &=\frac{[( \gamma_2^2 +\gamma_3^2) \tilde{F}_1^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_1\gamma_2 \tilde{F}_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_1\gamma_3 \tilde{F}_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta)]} { (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) [\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\\ &\quad + \frac{ U_1^0(\gamma_1 , \gamma_2 , \gamma_3)} {[\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\,, \end{aligned}\\ \label{eqn154} \begin{aligned} &U_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta )\\ &=\frac{[( \gamma_3^2 +\gamma_1^2) \tilde{F}_2^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_2\gamma_3 \tilde{F}_3^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_2\gamma_1 \tilde{F}_1^{\otimes} ( \gamma_1, \gamma_2, \gamma_3, \eta)]}{ (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) [\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\\ &\quad + \frac{ U_2^0(\gamma_1 , \gamma_2 , \gamma_3)}{[\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\,, \end{aligned} \\ \label{eqn155} \begin{aligned} &U_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta )\\ &=\frac{[( \gamma_1^2 +\gamma_2^2) \tilde{F}_3^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_3\gamma_1 \tilde{F}_1^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta) - \gamma_3\gamma_2 \tilde{F}_2^{\otimes}( \gamma_1, \gamma_2, \gamma_3, \eta)]} { (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) [\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\\ &\quad + \frac{ U_3^0(\gamma_1 , \gamma_2 , \gamma_3)} {[\eta+\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2)] }\,. \end{aligned} \end{gather} Then we use the convolution theorem with the convolution formula \eqref{A6} and integral \eqref{A7} for \eqref{eqn153}--\eqref{eqn155} to obtain \begin{gather}\label{eqn156} \begin{aligned} &U_1(\gamma_1, \gamma_2, \gamma_3, t)\\ &=\int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)}\\ &\quad\times \frac{[( \gamma_2^2 +\gamma_3^2) \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau) -\gamma_1\gamma_2 \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3, \tau ) -\gamma_1\gamma_3 \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, \tau ) ]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \\ &\quad +e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_1^0(\gamma_1 ,\gamma_2 ,\gamma_3), \end{aligned}\\ \label{eqn157} \begin{aligned} &U_2(\gamma_1, \gamma_2, \gamma_3, t)\\ &=\int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \\ &\quad\times \frac{[( \gamma_3^2 +\gamma_1^2) \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3, \tau) -\gamma_2\gamma_3 \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, \tau ) -\gamma_2\gamma_1 \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau ) ]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \\ &\quad +e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_2^0(\gamma_1 ,\gamma_2 ,\gamma_3), \end{aligned} \\ \label{eqn158} \begin{aligned} &U_3(\gamma_1, \gamma_2, \gamma_3, t)\\ &=\int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \\ &\quad\times \frac{[( \gamma_1^2 +\gamma_2^2) \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, \tau) -\gamma_3\gamma_1 \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau ) -\gamma_3\gamma_2 \tilde{F}_2 ( \gamma_1, \gamma_2, \gamma_3, \tau ) ]}{ (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \\ &\quad +e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_3^0(\gamma_1 ,\gamma_2 ,\gamma_3)\,. \end{aligned} \end{gather} Using the Fourier inversion formula \eqref{A3} we obtain \begin{align} &u_1(x_1, x_2, x_3, t) \nonumber \\ &=\frac{1}{(2\pi)^{3/2}} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \frac{[ ( \gamma_2^2 +\gamma_3^2) \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau)]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \,d\tau \nonumber\\ &\quad - \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)}\frac{ [\gamma_1\gamma_2 \tilde{F}_2 ( \gamma_1, \gamma_2, \gamma_3, \tau ) + \gamma_1\gamma_3 \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, \tau )]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \nonumber\\ &\quad + e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_1^0(\gamma_1 ,\gamma_2 ,\gamma_3)\Big] e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &=\frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{( \gamma_2^2 +\gamma_3^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2+\tilde x_3 \gamma_3)} \tilde{f}_1(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2 d\gamma_3 \nonumber \\ &\quad - \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ \gamma_1\gamma_2} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber\\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty} e^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2+\tilde x_3\gamma_3)} \tilde{f}_2(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2 d\gamma_3 \nonumber \\ &\quad- \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty}\frac{ \gamma_1\gamma_3} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1 +\tilde x_2\gamma_2+\tilde x_3\gamma_3)} \tilde{f}_3(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big]\nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)} \,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &\quad + \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t}\nonumber \\ &\quad\times \Big[ \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} u_1^0(\tilde x_1,\tilde x_2, \tilde x_3)\,d\tilde x_1d\tilde x_2 d\tilde x_3\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &= S_{11}(\tilde{f}_1)+ S_{12}(\tilde{f}_2)+ S_{13}(\tilde{f}_3)+B(u_1^0) \label{eqn160} \end{align} (see Remark \ref{rmkA8}); \begin{align} &u_2(x_1, x_2, x_3, t) \nonumber \\ &= \frac{1}{(2\pi)^{3/2}} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \frac{[( \gamma_3^2 +\gamma_1^2) \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3, \tau)]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \,d\tau \nonumber\\ &\quad - \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)}\frac{[\gamma_2\gamma_3 \tilde{F}_3 ( \gamma_1, \gamma_2, \gamma_3, \tau ) +\gamma_2\gamma_1 \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau )]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \nonumber \\ &\quad +e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_2^0(\gamma_1 ,\gamma_2 ,\gamma_3)\Big] e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &= -\frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty}\frac{ \gamma_2\gamma_1} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} \tilde{f}_1(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1 d\gamma_2d\gamma_3 \nonumber \\ &\quad + \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ (\gamma_3^2 + \gamma_1^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} \tilde{f}_2(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber\\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)} \,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &\quad - \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ \gamma_2\gamma_3} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber\\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} \tilde{f}_3(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber\\ &\quad \times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1 d\gamma_2d\gamma_3 \nonumber \\ &\quad + \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t}\nonumber \\ &\quad\times \Big[ \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} u_2^0(\tilde x_1,\tilde x_2, \tilde x_3)\,d\tilde x_1d\tilde x_2 d\tilde x_3\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1 d\gamma_2d\gamma_3 \nonumber \\ &= S_{21}(\tilde{f}_1)+ S_{22}(\tilde{f}_2)+ S_{23}(\tilde{f}_3)+B(u_2^0), \label{eqn161} \end{align} (see Remark \ref{rmkA8}); % \begin{align} &u_3(x_1, x_2, x_3, t) \nonumber \\ &= \frac{1}{(2\pi)^{3/2}} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \frac{[( \gamma_1^2 +\gamma_2^2) \tilde{F}_3( \gamma_1, \gamma_2, \gamma_3, \tau)]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \,d\tau \nonumber \\ &\quad - \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \frac{[\gamma_3\gamma_1 \tilde{F}_1( \gamma_1, \gamma_2, \gamma_3, \tau ) +\gamma_3\gamma_2 \tilde{F}_2( \gamma_1, \gamma_2, \gamma_3, \tau )]} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }\,d\tau \nonumber \\ &\quad + e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} U_3^0 (\gamma_1 ,\gamma_2 ,\gamma_3)\Big] e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)} \,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &=-\frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ \gamma_3\gamma_1} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} \tilde{f}_1(\tilde x_1,\tilde x_2, \tilde x_3,\tau) \,d\tilde x_1d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1d\gamma_2 d\gamma_3 \nonumber \\ &\quad -\frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ \gamma_3\gamma_2} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2 +\tilde x_3\gamma_3)} \tilde{f}_2(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad \times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)} \,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &\quad +\frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{ (\gamma_1^2 + \gamma_2^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \Big[ \int_0^t e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) (t-\tau)} \nonumber \\ &\quad\times \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty} e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2+\tilde x_3\gamma_3)} \tilde{f}_3(\tilde x_1,\tilde x_2, \tilde x_3,\tau)\,d\tilde x_1 d\tilde x_2 d\tilde x_3d\tau\Big] \nonumber \\ &\quad \times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)}\,d\gamma_1 d\gamma_2d\gamma_3 \nonumber \\ &\quad + \frac{1}{8\pi^3} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e ^{-\nu (\gamma_1^2 +\gamma_2^2 +\gamma_3^2) t} \nonumber \\ &\quad\times \Big[ \int_{-\infty}^{\infty} \int_{-\infty}^{\infty}\int_{-\infty}^{\infty} e ^{i(\tilde x_1\gamma_1+\tilde x_2\gamma_2+\tilde x_3\gamma_3)} u_3^0(\tilde x_1,\tilde x_2, \tilde x_3)\,d\tilde x_1d\tilde x_2 d\tilde x_3\Big] \nonumber \\ &\quad\times e ^{-i(x_1\gamma_1+x_2\gamma_2+x_3\gamma_3)} \,d\gamma_1d\gamma_2d\gamma_3 \nonumber \\ &= S_{31}(\tilde{f}_1)+ S_{32}(\tilde{f}_2)+ S_{33}(\tilde{f}_3)+B(u_3^0) , \label{eqn162} \end{align} (see Remark \ref{rmkA8}). Here $S_{11}()$, $S_{12}()$, $S_{13}()$, $S_{21}()$, $S_{22}()$, $S_{23}()$, $S_{31}()$, $S_{32}()$, $S_{33}()$, $B()$ are integral operators, and satisfy $S_{12}()=S_{21}(),\quad S_{13}()=S_{31}() ,\quad S_{23}()=S_{32}()\,.$ From the three expressions above for $u_1,u_2,u_3$, it follows that the vector $\vec{u}$ can be represented as: $$\label{eqn164} \vec{u}=\bar{\bar{S}}\cdot\vec{\tilde{f}}+\bar{\bar{B}}\cdot\vec{u}^0 =\bar{\bar{S}}\cdot\vec{f}- \bar{\bar{S}}\cdot(\vec{u}\cdot\nabla)\vec{u} +\bar{\bar{B}}\cdot\vec{u}^0\,,$$ where $\vec{\tilde{f}}$ is determined by formula \eqref{eqn14}. Here $\bar{\bar{S}}$ and $\bar{\bar{B}}$ are the matrix integral operators: $\begin{pmatrix} S_{11} & S_{12} & S_{13} \\ S_{21} & S_{22} & S_{23} \\ S_{31} & S_{32} & S_{33} \end{pmatrix}, \quad \begin{pmatrix} B & 0 & 0 \\ 0 & B & 0 \\ 0 & 0 & B \end{pmatrix}.$ \section{Spaces $S$ and $\vec{TS}$} As in \cite{GC68, RR64}, we consider the space $S$ of all infinitely differentiable functions $\varphi$(x) defined in $N$-dimensional space $\mathbb{R}^N$ ($N = 3$), such that these functions tend to 0 as $| x | \to \infty$, as well as their derivatives of any order, more rapidly than any power of $1/|x|$. To define a topology in the space $S$ let us introduce countable system of norms $$\label{eqn200} \|\varphi\|_{p}= \sup_{x}\big\{ | x^{k}D^{q} \varphi(x)| ,\, |k|\leq p, | q |\leq p\big\} \quad p = 0, 1, 2,\dots,$$ where \begin{gather*} | x^{k}D^{q} \varphi(x)|= | x_1^{k_1}\dots x_{N}^{k_{N}} \frac{\partial^{q_1+\dots + q_{N}}\varphi(x)}{\partial {x_1^{q_1}} \dots \partial {x_{N}^{q_{N}}} }|,\\ k = (k_1, \dots , k_{N}),\quad q = (q_1, \dots , q_{N}),\quad x^{k} = x_1^{k_1}\dots x_{N}^{k_{N}} ,\\ D^{q} = \frac{\partial^{q_1+\dots + q_{N}}}{\partial {x_1^{q_1}}\dots \partial {x_{N}^{q_{N}}} },\quad q_1, \dots , q_{N} = 0, 1, 2, \dots\,. \end{gather*} Note that $S$ is a perfect space (complete countably normed space, in which the bounded sets are compact). The space $\vec{TS}$ of vector-functions $\vec{\varphi}$ is a direct sum of $N$ perfect spaces $S$ ($N = 3$) \cite{VT80}: $\vec{TS} = S \oplus S \oplus S.$ To define a topology in the space $\vec{TS}$ let us introduce countable system of norms $$\label{eqn201} \|\vec{\varphi}\|_{p}= \sum_{i = 1}^N\|\varphi_{i}\|_{p} = \sum_{i = 1}^N\sup_{x}\big\{ | x^{k}D^{q} \varphi_{i}(x)|,\, | k |\leq p, |q| \leq p\big\},$$ $N = 3$, $p = 0, 1, 2,\dots$. The Fourier transform maps the space $S$ onto the whole space $S$, and maps the space $\vec{TS}$ onto the whole space $\vec{TS}$ \cite{gS01,GC68}. \section{Equivalence of the Cauchy problem in differential form \eqref{eqn1}--\eqref{eqn3} and in integral form} Let us denote solution of \eqref{eqn1}--\eqref{eqn3} as \{$\vec{u}(x_1, x_2, x_3, t)$, p($x_1, x_2, x_3,$ t)\}; in other words let us consider the infinitely differentiable by $t \in [0,\infty)$ vector-function $\vec{u}(x_1, x_2, x_3, t) \in \vec{TS}$, and infinitely differentiable function $p(x_1, x_2, x_3, t) \in S$, that turn equations \eqref{eqn1} and \eqref{eqn2} into identities. Vector-function $\vec{u}(x_1, x_2, x_3, t)$ also satisfies the initial condition \eqref{eqn3} $(\vec{u}^0 (x_1, x_2, x_3)\in \vec{TS})$: $$\label{eqn202} \vec{u}(x_1, x_2, x_3, t)|_{t = 0}=\vec{u}^0(x_1, x_2, x_3)$$ Let us put $\{\vec{u}(x_1, x_2, x_3, t), p(x_1, x_2, x_3, t)\}$ into equations \eqref{eqn1}, \eqref{eqn2} and apply Fourier and Laplace transforms to the result identities considering initial condition \eqref{eqn3}. After all required operations (as in sections 2 and 3) we receive that vector-function $\vec{u}(x_1, x_2, x_3, t)$ satisfies integral equation $$\label{eqn203} \vec{u}=\bar{\bar{S}}\cdot \vec{f}- \bar{\bar{S}}\cdot(\vec{u}\cdot\nabla)\vec{u} +\bar{\bar{B}}\cdot\vec{u}^0\; = \bar{\bar{S}}^{\nabla}\cdot\vec{u}$$ Then the vector-function grad $p\in \vec{TS}$ is defined by equations \eqref{eqn1} where vector-function $\vec{u}$ is defined by \eqref{eqn203}. Here $\vec{f} \in \vec{TS}$, $\vec{u}^0 \in \vec{TS}$ and $\bar{\bar{S}},\bar{\bar{B}},\bar{\bar{S}}^{\nabla}$ are matrix integral operators. Vector-functions $\bar{\bar{S}}\cdot\vec{f}$, $\bar{\bar{B}}\cdot\vec{u}^0$, $\bar{\bar{S}}\cdot(\vec{u}\cdot\nabla)\vec{u}$ also belong $\vec{TS}$ since Fourier transform maps perfect space $\vec{TS}$ onto $\vec{TS}$. Going from the other side, let us assume that $\vec{u}(x_1, x_2, x_3, t)\in \vec{TS}$ is continuous in $t \in [0,\infty)$ solution of integral equation \eqref{eqn203}. Integral-operators $S_{ij}\cdot(\vec{u}\cdot\nabla)\vec{u}$ are continuous in $t \in [0,\infty)$ [see \eqref{eqn160}--\eqref{eqn162}]. From here we obtain that according to \eqref{eqn203}, $$\vec{u}(x_1, x_2, x_3, 0) = \vec{u}^0(x_1, x_2, x_3)$$ also that $\vec{u}(x_1, x_2, x_3, t)$ is differentiable by $t \in [0,\infty)$. As described before, the Fourier transform maps the perfect space $\vec{TS}$ on itself. Hence, $\{\vec{u}(x_1, x_2, x_3, t) \text{ and}\\ p(x_1, x_2, x_3, t)\}$ is the solution of the Cauchy problem \eqref{eqn1}--\eqref{eqn3}. From here we see that solving the Cauchy problem \eqref{eqn1}--\eqref{eqn3} is equivalent to finding continuous in $t \in [0,\infty)$ solution of integral equation \eqref{eqn203}. \section{The Caccioppoli-Banach fixed point principle} See \cite{KA64,VT80,WR73,KS01,GD03,ADL97}. Further we have $\vec{f}\equiv 0$. Let us rewrite integral equation \eqref{eqn203} with this condition as $$\label{eqn203a} \vec{u}=- \bar{\bar{S}}\cdot(\vec{u}\cdot\nabla)\vec{u}+\bar{\bar{B}} \cdot\vec{u}^0$$ Let us divide all parts of the integral equation \eqref{eqn203a} by some constant $V$, that we will define appropriately below. Then we receive modified integral equation equivalent to equation \eqref{eqn203a}: $$\label{eqn203ab} \vec{u}_{\scriptscriptstyle V}=- \bar{\bar{S}}_{\scriptscriptstyle V}\cdot(\vec{u}_{\scriptscriptstyle V}\cdot\nabla_{\scriptscriptstyle V})\vec{u}_{\scriptscriptstyle V} +\bar{\bar{B}}_{\scriptscriptstyle V}\cdot\vec{u}_{\scriptscriptstyle V}^0 = \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V})$$ here $$\label{eqn203z} \begin{gathered} \vec{u}_{\scriptscriptstyle V}=\frac{\vec{u}}{\scriptscriptstyle V},\quad \vec{u}_{\scriptscriptstyle V}^0=\frac{\vec{u}^0}{\scriptscriptstyle V},\quad \nabla_{\scriptscriptstyle V}=V\cdot\nabla,\\ x_{k{\scriptscriptstyle V}}=\frac{x_k}{\scriptstyle V},\quad \gamma_{k{\scriptscriptstyle V}}=V\cdot\gamma_k,\quad \nu_{\scriptscriptstyle V}=\frac{\nu}{{\scriptstyle V}^2},\quad 1\leq k \leq N\,. \end{gathered}$$ The selection of the constant $V$ is based on a priori estimate of the solution of the Cauchy problem for the Navier-Stokes equations \cite{oL69,LK63}: $$\label{eqn203c} \|{\vec{u}}\|_{L_2}\leq \|{\vec{u^0}}\|_{L_2}$$ where $\|\vec{u}\|_{L_2}=\Big(\int_{\mathbb{R}^3}|\vec{u}|^2dx \Big)^{1/2}.$ The space of vector-functions $\vec{TS}$ is divided into two subsets. First subset consists of the vector-functions rapidly rapidly''decreasing to zero for $| x | \to \infty$. This subset is defined by the inequality $$\label{eqn203c1} \|{\vec{u}}\|_{L_2}\leq B \|{\vec{u}}\|_{C^2}\,,$$ where $B$ is a fixed constant. Second subset consists of the vector-functions rapidly'' decreasing to zero for $| x | \to \infty$. This subset is defined by the inequality $$\label{eqn203h1} \|{\vec{u}}\|_{C^2}\leq B' \|{\vec{u}}\|_{L_2},$$ where $B'$ is a fixed constant. $\|\vec{u}\|_{C^2}=\sum_{i = 1}^N\|u_{i}\|_{C^2} = \sum_{i = 1}^N\sum_{q = 0}^2\max_{x} | D^{q} u_{i}(x)|,\quad N = 3.$ For the first subset by using appropriate sequence of inequalities \cite{GC68}, \eqref{eqn203c}, \eqref{eqn203c1}, $$\label{eqn203c2} \|{\vec{u}}\|_{L_2}\leq B \|{\vec{u}}\|_{C^2} \leq K \|{\vec{u}}\|_{L_2}\leq K \|{\vec{u^0}}\|_{L_2} \leq K B \|{\vec{u^0}}\|_{C^2}$$ we obtain $$\label{eqn203d} \|{\vec{u}}\|_{C^2}\leq K\|{\vec{u^0}}\|_{C^2},$$ where $K$ is a constant. For the second subset by using appropriate sequence of inequalities \cite{GC68}, \eqref{eqn203c}, \eqref{eqn203h1}, $$\label{eqn203c2m} \|{\vec{u}}\|_{C^2}\leq B' \|{\vec{u}}\|_{L_2}\leq B' \|{\vec{u^0}}\|_{L_2}$$ we obtain $$\label{eqn203ds} \|{\vec{u}}\|_{C^2}\leq B'\cdot\|{\vec{u^0}}\|_{L_2}\,.$$ Let us choose the constant $V$ as: $$\label{eqn203e} \max \{K\|{\vec{u^0}}\|_{C^2}, \; B'\|{\vec{u^0}}\|_{L_2}\} < V$$ Then we have $$\label{eqn203f} \|{\vec{u}_{\scriptscriptstyle V}}\|_{C^2} < 1$$ We can use the fixed point principle to prove existence and uniqueness of the solution of integral equation \eqref{eqn203ab}. For this purpose we will operate with the following properties of matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$: 1. The matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ depends continuously on its parameter $t \in [0,\infty)$ (based on formulas \eqref{eqn160}--\eqref{eqn162}). 2. The matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ maps vector-functions $\vec{u}_{\scriptscriptstyle V}$ from perfect space $\vec{TS}$ onto perfect space $\vec{TS}$. This property directly follows from the properties of Fourier transform \cite{GC68}, and the form of integrands of integral operators $S_{ij}, B$ (based on formulas \eqref{eqn160}--\eqref{eqn162}). Hence in this case we can consider the convergence of functions from $\vec{TS}$ not only in countable system of norms $\eqref{eqn201}$, but also in norm $C^2$. 3. Matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}$ is quadratic''. 4. $\|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}\cdot\vec{u}_{\scriptscriptstyle V}-\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla} \cdot\vec{u}_{\scriptscriptstyle V}'\|_{C^2}\;<\;\|\vec{u}_{\scriptscriptstyle V}-\vec{u}_{\scriptscriptstyle V}' \|_{C^2}$ for any $\vec{u}_{\scriptscriptstyle V},\;\vec{u}_{\scriptscriptstyle V}' \in \vec{TS}$ (${\vec{u}_{\scriptscriptstyle V}\;\neq \;\vec{u}_{\scriptscriptstyle V}'})$ and any $t \in [0,\infty)$ (based on properties 1, 2, 3 and formulas \eqref{eqn160}--\eqref{eqn162}, \eqref{eqn203f}). The properties mentioned above allow us to prove that matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ is a contraction operator. \begin{theorem}[\cite{KA64}] \label{thm1} The matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ maps the perfect space $\vec{TS}$ onto the perfect space $\vec{TS}$, and for any $\vec{u}_{\scriptscriptstyle V},\vec{u}_{\scriptscriptstyle V}'\in \vec{TS}$ (${\vec{u}_{\scriptscriptstyle V}\neq \vec{u}_{\scriptscriptstyle V}'}$) the condition 4 is valid. Then the matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ is a contraction operator; i.e., the following condition is true $$\label{eqn206} \|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}) -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla} (\vec{u}_{\scriptscriptstyle V}')\|_{C^2}\leq \alpha\cdot\|\vec{u}_{\scriptscriptstyle V}-\vec{u}_{\scriptscriptstyle V}' \|_{C^2}\;$$ where $\alpha<1$ and is independent from $\vec{u}_{\scriptscriptstyle V},\;\vec{u}_{\scriptscriptstyle V}' \in \vec{TS}$ for any $t \in [0,\infty)$. \end{theorem} \begin{proof} By contradiction, let us assume that the opposite is true. Then there exist such $\vec{u}_{{\scriptscriptstyle V}n},\vec{u}'_{{\scriptscriptstyle V}n} \in \vec{TS}$ $(n=1,2,\dots)$ and $\lim_{ n\to \infty}\vec{u}_{{\scriptscriptstyle V}n}, \lim_{ n\to \infty}\vec{u}'_{{\scriptscriptstyle V}n} \in \vec{TS}$ that $$\label{eqn207} \|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}n}) -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}'_{{\scriptscriptstyle V}n})\|_{C^2} =\alpha_n\cdot\|\vec{u}_{{\scriptscriptstyle V}n}-\vec{u}'_{{\scriptscriptstyle V}n} \|_{C^2}\quad n=1,2,\dots;\, \alpha_n\to 1$$ Then the limiting result in \eqref{eqn207} would lead to equality $\|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V})-\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}')\|_{C^2} =\|\vec{u}_{\scriptscriptstyle V}-\vec{u}_{\scriptscriptstyle V}' \|_{C^2},$ which contradicts condition 4. Hence, $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ is a contraction operator. \end{proof} Next, we have the existence and uniqueness of a solution \cite{KA64}. \begin{theorem} \label{thm2} Let us consider a contraction operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$. Then there exists the unique solution $\vec{u}_{\scriptscriptstyle V}^{*}$ of equation \eqref{eqn203ab} in space $\vec{TS}$ for any $t \in [0,\infty)$. Also in this case it is possible to obtain $\vec{u}_{\scriptscriptstyle V}^{*}$ as a limit of sequence $\{\vec{u}_{{\scriptscriptstyle V}n}\}$, where $\vec{u}_{{\scriptscriptstyle V},n+1} = \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla} (\vec{u}_{{\scriptscriptstyle V}n}) \quad n=0,1,\dots,$ and $\vec{u}_{\scriptscriptstyle V0}= 0$. The rate of conversion of the sequence $\{\vec{u}_{{\scriptscriptstyle V}n}\}$ to the solution can be defined from the inequality $$\label{eqn208} \|\vec{u}_{{\scriptscriptstyle V}n} -\vec{u}_{\scriptscriptstyle V}^{*}\|_{C^2} \leq \frac{\alpha^{n}}{(1 - \alpha)}\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2} \quad n=0,1,\dots\,.$$ \end{theorem} \begin{proof} It is clear that $\vec{u}_{{\scriptscriptstyle V},n+1} = \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}n}) ,\quad \vec{u}_{{\scriptscriptstyle V}n} = \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{{\scriptscriptstyle V},n-1}).$ It follows from \eqref{eqn206} that $\|\vec{u}_{{\scriptscriptstyle V},n+1} -\vec{u}_{{\scriptscriptstyle V}n}\|_{C^2} \leq \alpha\cdot\|\vec{u}_{{\scriptscriptstyle V}n}-\vec{u}_{{\scriptscriptstyle V},n-1} \|_{C^2}.$ Using similar inequalities one after another while decreasing $n$ we will obtain $\|\vec{u}_{{\scriptscriptstyle V},n+1} -\vec{u}_{{\scriptscriptstyle V}n}\|_{C^2} \leq \alpha^{n}\cdot\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2}.$ From this result it follows that \label{eqn209} \begin{aligned} \|\vec{u}_{{\scriptscriptstyle V},n+l} -\vec{u}_{{\scriptscriptstyle V}n}\|_{C^2} &\leq \|\vec{u}_{{\scriptscriptstyle V},n+l} -\vec{u}_{{\scriptscriptstyle V},n+l-1}\|_{C^2}+\dots +\|\vec{u}_{{\scriptscriptstyle V},n+1} -\vec{u}_{{\scriptscriptstyle V}n}\|_{C^2}\\ &\leq (\alpha^{n+l-1} +\dots+\;\alpha^{n})\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2}\\ &\leq\frac{\alpha^{n}}{(1 - \alpha)}\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2}. \end{aligned} Because of $\alpha^{n} \to 0$ for $n \to \infty$, the obtained estimate \eqref{eqn209} shows that sequence $\{\vec{u}_{{\scriptscriptstyle V}n}\}$ is a Cauchy sequence. Since the space $\vec{TS}$ is a perfect space, this sequence converges to an element $\vec{u}_{\scriptscriptstyle V}^{*}\in \vec{TS}$, such that $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})$ has sense. We use inequality $\eqref{eqn206}$ again and have: $\|\vec{u}_{{\scriptscriptstyle V},n+1} -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})\|_{C^2} =\|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}n}) -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})\|_{C^2} \leq \alpha \|\vec{u}_{{\scriptscriptstyle V}n}-\vec{u}_{\scriptscriptstyle V}^{*} \|_{C^2}\quad n=0,1,2,\dots$ The right part of the above inequality tends to 0 as $n \to \infty$ and it means that $\vec{u}_{{\scriptscriptstyle V},n+1} \to \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})$ and $\vec{u}_{\scriptscriptstyle V}^{*} = \bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})$. In other words, $\vec{u}_{\scriptscriptstyle V}^{*}$ is the solution of equation \eqref{eqn203ab}. Uniqueness of the solution also follows from \eqref{eqn206}. In fact, if there would exist another solution $\widetilde{\vec{u}_{\scriptscriptstyle V}} \in \vec{TS}$, then $\|\widetilde{\vec{u}_{\scriptscriptstyle V}} -\vec{u}_{\scriptscriptstyle V}^{*}\|_{C^2} =\|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\widetilde{\vec{u}_{\scriptscriptstyle V}}) -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}(\vec{u}_{\scriptscriptstyle V}^{*})\|_{C^2} \leq \alpha\cdot\|\widetilde{\vec{u}_{\scriptscriptstyle V}}-\vec{u}_{\scriptscriptstyle V}^{*} \|_{C^2}.$ Such situation could happen only if $\|\widetilde{\vec{u}_{\scriptscriptstyle V}} -\vec{u}_{\scriptscriptstyle V}^{*}\|_{C^2}=0$, or $\widetilde{\vec{u}_{\scriptscriptstyle V}}=\vec{u}_{\scriptscriptstyle V}^{*}$. We can also obtain the estimate \eqref{eqn208} from estimate \eqref{eqn209} as a limiting result as $l\to \infty$. Now let us show that continuous dependence of operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ on t leads to continuous dependence of the solution of the problem on $t$. We will say that matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ is continuous in $t$ at a point $t_0 \in [0,\infty)$, if for any sequence $\{t_n\} \in [0,\infty)$ with $t_n \to t_0$ for $n\to \infty$, the following is true: $$\label{eqn210} \bar{\bar{S}}_{{\scriptscriptstyle V}t_n}^{\nabla}(\vec{u}_{\scriptscriptstyle V})\to \bar{\bar{S}}_{{\scriptscriptstyle V}t_0}^{\nabla}(\vec{u}_{\scriptscriptstyle V})\quad \text{for any }\vec{u}_{\scriptscriptstyle V} \in \vec{TS}.$$ \end{proof} From Theorem \ref{thm2} it follows that for any $t \in [0,\infty)$, equation \eqref{eqn203ab} has the unique solution, which depends on $t$. Let us denote it as $\vec{u}_{{\scriptscriptstyle V}t}^{*}$. We will say that solution of equation \eqref{eqn203ab} depends continuously on $t$ at $t=t_0$, if for any sequence $\{t_n\} \in [0\infty)$ with $t_n \to t_0$ for $n\to \infty$, the following is true: $\vec{u}_{{\scriptscriptstyle V}t_n}^{*}\to \vec{u}_{{\scriptscriptstyle V}t_0}^{*}.$ Next, we have the continuous dependence of the solution on $t$. \begin{theorem}[\cite{KA64}] \label{thm3} Let us consider operator $\bar{\bar{S}}_{{\scriptscriptstyle V}t}^{\nabla}$ that satisfies condition \eqref{eqn206} for any $t \in [0,\infty)$, where $\alpha$ is independent from $t$ and that operator $\bar{\bar{S}}_{{\scriptscriptstyle V}t}^{\nabla}$ is continuous in $t$ at a point $t_0 \in [0,\infty)$. Then for $t=t_0$ the solution of \eqref{eqn203ab} depends continuously on $t$. \end{theorem} \begin{proof} Let us consider any $t \in [0,\infty)$. We will construct the solution $\vec{u}_{{\scriptscriptstyle V}t}^{*}$ of equation \eqref{eqn203ab} as a limit of sequence $\{\vec{u}_{{\scriptscriptstyle V}n}\}$: $$\label{eqn211} \vec{u}_{{\scriptscriptstyle V},n+1} = \bar{\bar{S}}_{{\scriptscriptstyle V}t}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}n}) \quad n=0,1,\dots;\;\vec{u}_{{\scriptscriptstyle V}0}=\vec{u}_{{\scriptscriptstyle V}t_0}^{*}$$ Let us rewrite inequality \eqref{eqn208} for n = 0: $$\label{eqn212} \|\vec{u}_{\scriptscriptstyle V}^{*} -\vec{u}_{{\scriptscriptstyle V}0}\|_{C^2}\leq \frac{1}{(1 - \alpha)}\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2}$$ Since $\vec{u}_{{\scriptscriptstyle V}t_0}^{*} = \bar{\bar{S}}_{{\scriptscriptstyle V}t_0}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}t_0}^{*})$, because of \eqref{eqn211} and \eqref{eqn212} we have $$\label{eqn213} \|\vec{u}_{{\scriptscriptstyle V}t}^{*} -\vec{u}_{{\scriptscriptstyle V}t_0}^{*}\|_{C^2} \leq \frac{1}{(1 - \alpha)}\|\vec{u}_{{\scriptscriptstyle V}1}-\vec{u}_{{\scriptscriptstyle V}0} \|_{C^2} =\frac{1}{(1 - \alpha)}\|\bar{\bar{S}}_{{\scriptscriptstyle V}t}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}t_0}^{*}) -\bar{\bar{S}}_{{\scriptscriptstyle V}t_0}^{\nabla}(\vec{u}_{{\scriptscriptstyle V}t_0}^{*}) \|_{C^2}$$ Now with the help of \eqref{eqn210} we obtain the required continuity of $\vec{u}_{{\scriptscriptstyle V}t}$ for $t=t_0$. \end{proof} Following \eqref{eqn203z} and \eqref{eqn203e} we obtain the result: $$\label{eqn203z1} \vec{u}=\vec{u}_{\scriptscriptstyle V}\cdot V,\quad \nu=\nu_{\scriptscriptstyle V}\cdot V^2.$$ Then vector-function $\nabla p \in \vec{TS}$ is defined by \eqref{eqn1} where vector-function $\vec{u}$ is received from equation \eqref{eqn203z1}. Function $p$ is defined up to an arbitrary constant. \begin{remark} \rm From the above statements, it follows that there exists the unique set of smooth functions $u_{\infty i}(x, t)$, $p_{\infty}(x, t)$ $(i = 1, 2, 3)$ $\mathbb{R}^3 \times [0,\infty)$ that satisfies \eqref{eqn1}, \eqref{eqn2}, \eqref{eqn3} and $$\label{eqn186b} u_{\infty i},p_{\infty} \in C^{\infty}(\mathbb{R}^3 \times [0,\infty)),$$ Then, using the inequality $\|{\vec{u}}\|_{L_2}\leq \|{\vec{u^0}}\|_{L_2}$ from \eqref{eqn203c}, \cite{oL69}, \cite{LK63}, we have $$\label{eqn186c} \int_{\mathbb{R}^3}|\vec{u}_{\infty}(x, t)|^2dx < C,\quad \forall t\geq 0.$$ Let us consider $\nu \to$ 0 in integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$. Then we see that Theorems \ref{thm1}--\ref{thm3} are correct also in case of Euler equations; i.e., there exists unique smooth solution in all time range for this case. \end{remark} \section{Appendix} The Fourier integral can be stated in the forms: \label{A3} \begin{gathered} \begin{aligned} U( \gamma_1 , \gamma_2 , \gamma_3) &=F[ u(x_1 , x_2 , x_3)]\\ &= \frac{1}{(2\pi)^{3/2}} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} u( x_1 , x_2 , x_3) e ^{ i( \gamma_1 x_1 + \gamma_2 x_2 + \gamma_3 x_3) } dx_1 dx_2 dx_3 \end{aligned}\\ u( x_1 , x_2 , x_3)= \frac{1}{(2\pi)^{3/2}} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} U( \gamma_1 , \gamma_2,\gamma_3 )\, e ^{- i( \gamma_1 x_1 + \gamma_2 x_2 + \gamma_3 x_3) } d\gamma_1 d\gamma_2 d\gamma_3 \end{gathered} The Laplace integral is usually stated in the form $$\label{A4} U^{\otimes}(\eta)=L[u(t)]= \int_0^{\infty}u(t) e ^{-\eta t}dt \quad u(t)=\frac{1}{2\pi i}\int_{c- i \infty }^{c + i \infty} U^{\otimes}(\eta) e ^{\eta t}d\eta \quad c > c_0.$$ Then $$\label{A5} L[u'(t)]=\eta U^{\otimes}(\eta)-u(0).$$ The convolution theorem \cite{DP65,DW46} is stated as: If integrals $U_1^{\otimes}(\eta)= \int_0^{\infty}u_1(t) e ^{-\eta t}d\,t \quad U_2^{\otimes}(\eta)= \int_0^{\infty}u_2(t) e ^{-\eta t}d\,t$ converge absolutely for $\operatorname{Re} \eta > \sigma_{d}$, then $U^{\otimes}(\eta)=U_1^{\otimes}(\eta) U_2^{\otimes}(\eta)$ is Laplace transform of $$\label{A6} u(t)=\int_0^tu_1(t-\tau)\,u_2(\tau)\,d\,\tau$$ A useful Laplace integral is $$\label{A7} L[e ^{\eta_kt}]=\int_0^{\infty}e ^{-(\eta-\eta_k)\,t}d\,t = \frac{1}{(\eta-\eta_k)}\quad \operatorname{Re}\eta >\eta_k$$ \begin{remark} \label{rmkA8} \rm In the calculations of integrals \eqref{eqn160}--\eqref{eqn162} for components of velocity $u_1, u_2, u_3$ for the inverse Fourier transforms, we have each integrand $\tilde{f}(\gamma_1,\gamma_2,\gamma_3$) as a product of functions $\chi(\gamma_1,\gamma_2,\gamma_3)$ and $\varphi(\gamma_1,\gamma_2,\gamma_3)$, $\tilde{f}(\gamma_1,\gamma_2,\gamma_3) = \chi(\gamma_1,\gamma_2,\gamma_3)\cdot \varphi(\gamma_1,\gamma_2,\gamma_3),$ where $\varphi(\gamma_1,\gamma_2,\gamma_3)$ belongs to space $S$ (functions of $\gamma_1,\gamma_2,\gamma_3$) \cite{GC68} and $\chi(\gamma_1,\gamma_2,\gamma_3)$ is one of the fractions: \begin{gather*} \frac{( \gamma_2^2 +\gamma_3^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) },\quad \frac{( \gamma_1\cdot\gamma_2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) },\quad \frac{( \gamma_1\cdot\gamma_3)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) }, \\ \frac{( \gamma_3^2 +\gamma_1^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) },\quad \frac{( \gamma_2\cdot\gamma_3)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) },\quad \frac{( \gamma_1^2 +\gamma_2^2)} { (\gamma_1^2 +\gamma_2^2+\gamma_3^2 ) } \end{gather*} These fractions are infinitely differentiable functions for $\gamma_1\neq 0$, $\gamma_2\neq 0$, $\gamma_3\neq 0$ with one point of discontinuity $\gamma_1=0$, $\gamma_2=0$, $\gamma_3=0$ (The discontinuities have finite values at this point). For these calculations the inverse Fourier transforms are defined as Lebesgue integrals with Cauchy principal values. \end{remark} \begin{theorem} \label{thmA61}: The inverse Fourier transform of $\tilde{f}(\gamma_1,\gamma_2,\gamma_3) = \chi(\gamma_1,\gamma_2,\gamma_3) \cdot \varphi(\gamma_1,\gamma_2,\gamma_3)$, $$F[\tilde{f}]\equiv \psi(\sigma)\equiv \int_{-\infty}^{\infty} e ^{ i( \gamma,\sigma) }\chi(\gamma)\cdot\varphi(\gamma) d\gamma \label{A61}$$ as function of $\sigma$ belongs to space $S$ (functions of $\sigma$); i.e., $\psi(\sigma)$ has two criteria: \begin{itemize} \item[(1)] $\psi(\sigma)$ is infinitely differentiable function, \item[(2)] when $| \sigma | \to \infty$, $\psi(\sigma)$ tends to 0, as well as its derivatives of any order, more rapidly than any power of $1/|\sigma|$. \end{itemize} \end{theorem} Integral in \eqref{A61} admits of differentiation with respect to the parameter $\sigma_{j}$, since the integral obtained after formal differentiation remains absolutely convergent: $$\frac{\partial \psi(\sigma)}{\partial \sigma_{j}} \equiv \int_{-\infty}^{\infty} i \gamma_{j} e ^{ i( \gamma,\sigma) } \chi(\gamma) \varphi(\gamma) d\gamma \label{A62}$$ The properties of function $\varphi(\gamma)$ permit this differentiation to be continued without limit. This means that the function $\psi(\sigma)$ is infinitely differentiable (see criterion 1). To prove criterion 2 we create a function with parameter $n$: $f_n(\gamma_1,\gamma_2,\gamma_3) = n e ^{- \frac{1}{n^2 (\gamma_1^2 +\gamma_2^2+\gamma_3^2 )}} \chi(\gamma_1,\gamma_2,\gamma_3) \cdot \varphi(\gamma_1,\gamma_2,\gamma_3)$ $n = 1,2,3 \dots$. Then $f_n(\gamma_1,\gamma_2,\gamma_3$) belongs to space $S$ (functions of $\gamma_1,\gamma_2,\gamma_3$). Then we estimate the integral in formula \eqref{A61} using the inverse Fourier transform of $f_n(\gamma)$ for $| \sigma |>> 0$, $| \sigma | \to \infty$: $$|\psi(\sigma)| \equiv \big|\int_{-\infty}^{\infty} e ^{ i( \gamma,\sigma) } \chi(\gamma) \varphi(\gamma) d\gamma \big| \leq \big| n \int_{-\infty}^{\infty} e ^{ i( \gamma,\sigma) } e ^{- \frac{1}{n^2 \gamma^2 }}\chi(\gamma) \varphi(\gamma) d\gamma \big|,\label{A63}$$ $n >> 1, \dots <\infty$. Since $f_n(\gamma_1,\gamma_2,\gamma_3$) belongs to space $S$ (functions of $\gamma_1,\gamma_2,\gamma_3$) then the inverse Fourier transform of $f_n(\gamma_1,\gamma_2,\gamma_3$) belongs to space $S$ (functions of $\sigma_1,\sigma_2,\sigma_3$). Hence $F[f_n]$ is the function, such that when $| \sigma | \to \infty$ this function tends to $0$, as well as its derivatives of any order, more rapidly than any power of $1/|\sigma|$. So from formula \eqref{A63}, we have that $\psi(\sigma)$ is a function, such that when $| \sigma | \to \infty$ this function tends to $0$, more rapidly than any power of $1/|\sigma|$. Then we estimate the integral in formula \eqref{A62}. Using formula like \eqref{A63} for $\frac{\partial \psi(\sigma)}{\partial \sigma_{j}}$, we obtain that $\frac{\partial \psi(\sigma)}{\partial \sigma_{j}}$ tends to 0 more rapidly than any power of $1/|\sigma|$. The properties of function $\varphi(\gamma)$ permit this differentiation to be continued without limit and further using of formulas like \eqref{A63} to derivatives leads to the conclusion that all derivatives of function $\psi(\sigma)$ tend to 0 more rapidly than any power of $1/|\sigma|$. We have proved that $\psi(\sigma)$ belongs to space $S$ (functions of $\sigma$). \subsection*{Acknowledgments} We express our sincere gratitude to Professor L. Nirenberg, whose suggestion led to conduction of this research. We are also very thankful to Professor A. B. Gorstko for helpful friendly discussions. \begin{thebibliography}{35} \bibitem{ADL97} J. M. Ayerbe Toledano, T. Dominguez Benavides, G. Lopez Acedo; \emph{Measures of Noncompactness in Metric Fixed Point Theory,} Birkhauser Verlag, Basel - Boston - Berlin, 1997. \bibitem{BM02} A. Bertozzi, A. Majda; \emph{Vorticity and Incompressible Flows,} Cambridge U. Press, Cambridge, 2002. \bibitem{SB59} S. Bochner; \emph{Lectures on Fourier integrals; with an author's supplement on monotonic functions, Stieltjes integrals, and harmonic analysis.} Princeton, N.J., Princeton University Press, 1959., 333. \bibitem{CKN82} L. Caffarelli, R. Kohn, L. Nirenberg; \emph{Partial regularity of suitable weak solutions of the Navier-Stokes equations.} Communications on Pure and Applied Math. 35 (1982), 771-831. \bibitem{pC01} P. Constantin; \emph{Some open problems and research directions in the mathematical study of fluid dynamics,} in Mathematics Unlimited-2001 and Beyond, Springer Verlag, Berlin, 2001, 353-360. \bibitem{DP65} V. A. Ditkin, A.P. Prudnikov; \emph{Integral transforms and operational calculus.} Pergamon Press in Oxford, 1965. \bibitem{GC68} I. M. Gel'fand, G. E. Chilov; \emph{Generalized functions./Volume 2, Spaces of fundamental and generalized functions.} New York; London: Academic Press, 1968. \bibitem{GD03} A. Granas, J. Dugundji; \emph{Fixed Point Theory,} Springer-Verlag, New York, 2003. \bibitem{LH83} L. Hormander; \emph{The Analysis of Linear Partial Differential Operators I - IV.} Berlin, New York, Springer Verlag, 1983 - 1985. \bibitem{KA64} L. V. Kantorovich, G. P. Akilov; \emph{Functional analysis in normed spaces.} Oxford London Edinburgh New York Paris Frankfurt, Pergamon Press., 1964. \bibitem{KS01} W. A. Kirk, B. Sims; \emph{Handbook of Metric Fixed Point Theory,} Kluwer Academic, London, 2001. \bibitem{LK63} O. Ladyzhenskaya, A. Kiselev; \emph{On the existence and uniqueness of the solution of the non-stationary problem for a viscous incompressible fluid,} Izv. Akad.Nauk SSSR Ser. Mat. 21 (1957), 665-680; English transl., Amer. Math. Soc. Transl. (2) 24 (1963), 79-106. \bibitem{oL69} O. Ladyzhenskaya; \emph{The Mathematical Theory of Viscous Incompressible Flows,} (2nd edition), Gordon and Breach, 1969. \bibitem{pL02} P. G. Lemari\'e-Rieusset; \emph{Recent developments in the Navier-Stokes problem,} CRC Press (Boca Raton), Research Notes in Mathematics Series , 2002. \bibitem{jL34} J. Leray; \emph{Sur le Mouvement d'un Liquide Visquex Emplissent l'Espace.} Acta Math. J. 63 (1934), 193-248. \bibitem{fL98} F.-H. Lin; \emph{A new proof of the Caffarelli-Kohn-Nirenberg theorem.} Communications on Pure and Applied Math. 51 (1998), 241-257. \bibitem{SM73} S. Mizohata; \emph{The Theory of Partial Differential Equations.} Cambridge Univ. Press, 1973. \bibitem{VP70} V. P. Palamodov; \emph{Linear Differential Operators with Constant Coefficients.} Berlin, New York, Springer-Verlag, 1970. \bibitem{RR64} A. P. Robertson, W. J. Robertson; \emph{Topological Vector Spaces.} Cambrige University Press, 1964. \bibitem{WR73} W. Rudin; \emph{Functional analysis.} New York St. Louis San Francisco Dusseldorf Johannesburg Kuala Lumpur London Mexico Montreal New Delhi Panama Rio de Janeiro Singapore Sydney Toronto, McGraw-Hill book company, 1973. \bibitem{vS76} V. Scheffer; \emph{Turbulence and Hausdorff dimension, in Turbulence and the Navier-Stokes Equations.} Lecture Notes in Math. No. 565, Springer Verlag, 1976, pp. 94-112. \bibitem{vS93} V. Scheffer; \emph{An inviscid flow with compact support in spacetime.} J. Geom. Analysis 3 No. 4 (1993), 343-401. \bibitem{gS01} G. E. Shilov; \emph{Elementary functional analysis,} Cambridge, Mass.: MIT Press, 1974. \bibitem{aS97} A. Shnirelman; \emph{On the nonuniqueness of weak solutions of the Euler equation.} Communications on Pure and Applied Math. 50 (1997), 1260-1286. \bibitem{RT77} R. Temam; \emph{Navier-Stokes Equation: Theory and Numerical Analysis,} North-Holland Pub. Co, 1977. \bibitem{VT80} V. A. Trenogin; \emph{Functional'nyiy analiz (Functional analysis).} Nauka, Moskva, GRFML, Russian, 1980. \bibitem{JFT61} J. F. Treves; \emph{Lectures on linear partial differential equations with constant coefficients.} Rio de Janeiro: Instituto de Matem\'atica Pura e Aplicada do Conselho Nacional de Pesquisas, 1961. \bibitem{TT10} A. Tsionskiy, M. Tsionskiy; \emph{Solution of the Cauchy problem for the Navier - Stokes and Euler equations,} arXiv:1009.2198v3, 2010. \bibitem{TT11} A. Tsionskiy, M. Tsionskiy; \emph{Research of convergence of the iterative method for solution of the Cauchy problem for the Navier - Stokes equations based on estimated formula,} arXiv:1101.1708v2, 2011. \bibitem{DW46} D. W. Widder; \emph{The Laplace Transform.} Princeton, 1946. \end{thebibliography} \section*{Addendum posted by the editor on December 31, 2013} Two anonymous readers informed us that the results in this article are incorrect. I sent the questions to the authors who answer them, and posted a new version of their article in arXiv in September 2013; see reference [2] below. However, the answers seem to be still erroneous. Here is an extract of our correspondence. \medskip \noindent\textbf{Question 1.} The proof of Point 4 page 12 (which is the core of the paper) does not make sense. I doubt that in the book by Gelfand one find the appropriate sequence of inequalities'' allowing to control a $C^2$ norm by a $L^2$ norm. \noindent\textbf{Answer 1.} To improve the estimate of the constant $V$ we have applied an approach that avoids the sequence of inequalities. Let us choose the constant $V$ as: $\max\{\|{\vec{u}^0}\|_{C^2}, \|{\vec{u}^0}\|_{L_2}, (\|{\vec{u}}\|_{C_2} + \|{\vec{u}}\|_{L_2})\}\ll V$ Then we have $|{\vec{u}_{\scriptscriptstyle V}}|<1,\quad |{\vec{u}^0_{\scriptscriptstyle V}}|<1$ In property 4 of the matrix integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$, and in Theorems 6.1--6.3, where two vector-functions $\vec{u}_{\scriptscriptstyle V}$ and $\vec{u}_{\scriptscriptstyle V}'$ are presented, constant $V$ is selected as maximum between $V(\vec{u}$) and $V'(\vec{u}')$. In other words we do not use appropriate sequence of inequalities'' in the new version, see reference [2] below. \noindent\textbf{Reader's reply:} As $V$ is constant and as $\vec{u}$ is time-varying, does the inequality holds on $[0;1)$? Moreover, the role of $V$ is still more dubious. Indeed, the aim seems to ensure that $\vec{u}$ is bounded by 1. But the role of this boundedness is unclear. It seems that the author want to use the fact that $S$ is a perfect set (bounded subsets are relatively compact); but boundedness in $S$ means being bounded for a whole countable family of semi-norms involving polynomial weights and higher derivatives, so that the mere boundedness of the modulus of $\vec{u}$ does not bear any compactness property in $S$. \medskip \noindent\textbf{Question 2.} The proof of Theorem 6.1 is false, since $u$ could be equal to $u'$, leading to no contradiction with point 4. \noindent\textbf{Answer 2.} Below you can see the condition (${\vec{u}_{\scriptscriptstyle V}\neq \vec{u}_{\scriptscriptstyle V}'}$) in Theorem 6.1. The matrix integral operator $\bar{\bar{S}_{\scriptscriptstyle V}}^{\nabla}$ maps the perfect space $\vec{TS}$ onto the perfect space $\vec{TS}$, and for any $\vec{u}_{\scriptscriptstyle V},\vec{u}_{\scriptscriptstyle V}'\in \vec{TS}$ (${\vec{u}_{\scriptscriptstyle V}\neq \vec{u}_{\scriptscriptstyle V}'}$) the condition 4 is valid. Then the matrix integral operator $\bar{\bar{S}_{\scriptscriptstyle V}}^{\nabla}$ is a contraction operator, i.e. the following condition is true: $\|\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}\cdot\vec{u}_{\scriptscriptstyle V} -\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla} \cdot\vec{u}_{\scriptscriptstyle V}'\|_{p} \leq \alpha\cdot\|\vec{u}_{\scriptscriptstyle V}-\vec{u}_{\scriptscriptstyle V}' \|_{p}\tag{6.13}$ where $\alpha<1$ and is independent from $\vec{u}_{\scriptscriptstyle V},\vec{u}_{\scriptscriptstyle V}' \in \vec{TS}$ for any t $\in [0,\infty)$. \noindent\textbf{Reader's reply:} The condition $\vec{u}\neq \vec{u}'$ is stated in the assumption of Theorem 6.1. However, the proof uses two limit points of sequences $\vec{u}_n$ and $\vec{u}'_n$ with $\vec{u}\neq \vec{u}'$ and those two limit points may satisfy $\vec{u}=\vec{u}'$ and thus the alleged contradiction is not proved. As a matter of fact, there are a lot of examples of contraction on compacts spaces that are not strictly contractive: for instance $t\mapsto \sin t$ on $[0, \pi/2]$: we have $|\sin t- \sin t'| < |t-t'|$ when $t\neq t'$, however there is no $\alpha < 1$ such that $|\sin t- \sin t'| \leq \alpha |t-t'|$ for all $t, t'\in [0,\pi/2]$ (as the derivative of $\sin t$ at $t = 0$ equals $1$). \medskip \noindent\textbf{Question 3.} The proof of Theorem 6.2 is false: convergence in $C^2$ norm of a sequence of Schwartz functions does not imply that the limit is a Schwartz function \noindent\textbf{Answer 3.} We use $\|\cdot\|_{p}$ countable system of norms (4.2) for the property 4 of the matrix integral operator $\bar{\bar{S}_{\scriptscriptstyle V}}^{\nabla}$, and in Theorems 6.1-6.3 now. The limit of a sequence of Schwartz functions is also a Schwartz function by $\|\cdot\|_{p}$-countable system of norms (4.2). \noindent\textbf{Reader's reply:} The convergence in $C^2$ norm has been replaced by convergence in the countable family of norms that define the topology of $S$. This is better than in the published version, except that the convergence in those norms is not proved. The new Theorem 6.1 is false: the assertion that $\|S_{\scriptscriptstyle V}^\nabla \vec{u}_{\scriptscriptstyle V}-S_{\scriptscriptstyle V}^\nabla \vec{u}'_{\scriptscriptstyle V}\|_p < \|\vec{u}_{\scriptscriptstyle V}- \vec{u}'_{\scriptscriptstyle V}\|_p$ (condition 4) follows from formulas (3.30)-(3.31)-(3.32) and (6.5) is purely an act of faith, but not a plausible proof. There is absolutely no scientific reason that this infinite family of inequalities holds for a providential choice of $V$. \medskip \noindent\textbf{Question 4.} Remark 6.4 is strange. How can we deal the case of the vanishing viscosity? \noindent\textbf{Answer 4.} Let us consider $\nu \rightarrow 0$. We see that integral operator $\bar{\bar{S}}_{\scriptscriptstyle V}^{\nabla}$ and Theorems 6.1-6.3 are correct also in case of Euler equations, i.e. there exists unique smooth solution in all time range for this case too. \noindent\textbf{Reader's reply:} The answer that we see that theorems 6.1, 6.2 and 6.3 are correct in case of Euler equations'' without any further explanations on how we can see it is meaningless. \medskip \noindent\textbf{Question 5.} Theorem 7.2 is incorrect, as a discontinuous function cannot belong to the Schwartz class \noindent\textbf{Answer 5.} $\varphi(\gamma)$ as function of $\gamma$ belongs to space $S$ (functions of $\gamma$). I.e., $\varphi(\gamma)$ has two properties: \\ (1) $\varphi(\gamma)$ is infinitely differentiable function, \\ (2) when $| \gamma | \to \infty$, $\varphi(\gamma)$ tends to 0, as well as its derivatives of any order, more rapidly than any power of $1/| \gamma|$. Function $\tilde{f}(\gamma_1,\gamma_2,\gamma_3) = \chi(\gamma_1,\gamma_2,\gamma_3)\cdot \varphi(\gamma_1,\gamma_2,\gamma_3),$ does not belong to space $S$ (functions of $\gamma$) then $\chi(\gamma_1,\gamma_2,\gamma_3)$ is the discontinuous function for $\gamma = 0$. But the inverse Fourier transform of function $\tilde{f}(\gamma_1,\gamma_2,\gamma_3)$ $F[\tilde{f}]\equiv \psi(\sigma) \equiv \int_{-\infty}^{\infty} e ^{ i( \gamma,\sigma) }\chi(\gamma) \cdot\varphi(\gamma) d\gamma$ $\psi(\sigma)$ as function of $\sigma$ belongs to space $S$ (functions of $\sigma$). I.e., $\psi(\sigma)$ has two properties: \\ (1) $\psi(\sigma)$ is infinitely differentiable function, \\ (2) when $| \sigma |\to \infty$ $\psi(\sigma)$ tends to $0$, as well as its derivatives of any order, more rapidly than any power of $1/|\sigma |$. We can see the correct proof of two properties (1) and (2) for $\psi(\sigma)$ in Theorem 7.2. \noindent\textbf{Reader's reply:} Theorem 7.2 is incorrect. It is a basic fact of the theory of distributions that the Fourier transform is a bijection of the Schwartz class onto itself, so that obviously a non-smooth function cannot have an inverse Fourier transform in the Schwartz class. The author still argues that his proof is correct. The proof is based on the inequality $\big|\int e^{i(\gamma,\sigma)} \chi(\gamma)\,d\gamma\big| \leq n\big|\int e^{i(\gamma,\sigma)} e^{-1/(n^2\gamma^2)} \chi(\gamma) \varphi(\gamma)\,d\gamma \big|$ which is not proved (and false). As a matter of fact, the inequality is valid for large $n$ when $\sigma$ is fixed, but is false for large $\sigma$ when $n$ is fixed. \medskip \noindent\textbf{Question 6.} What are the main properties of the solution? \noindent\textbf{Answer 6.} Solution $\vec{u} \in \vec{TS}$ and $p \in S$ In other words there exists the unique set of smooth functions $u_{\infty i}(x, t)$, $p_{\infty}(x, t)$ $(i = 1, 2, 3)$ on $\mathbb{R}^3 \times [0,\infty)$ that satisfies (2.1), (2.2), (2.3), and $u_{\infty i}, p_{\infty} \in C^{\infty}(\mathbb{R}^3 \times [0,\infty)),$ Then, using the inequality $\|\vec{u}\|_{L_2}\leq\|\vec{u}^0\|_{L_2}$, we have $\int_{\mathbb{R}^3}|\vec{u}_{\infty}(x, t)|^2dx < C$ for all $t \geq 0$, see Fefferman [1], below. \medskip \noindent\textbf{Question 7.} I think that (6.7) and (6.9) are wrong ($L^2$ never controls $C^2$ with universal constants). After that, the authors get (6.8). From there using Beale-Kato-Majda they may conclude the global existence. Without (6.7) and (6.9), their definition of $V$ makes no sense, as the new $V$ will depend on time. \noindent\textbf{Answer 7.} Instead inequalities (6.5)--(6.12) we choose the constant $V$ as: $\max \{\|{\vec{u}^0}\|_{C^2}, \|{\vec{u}^0}\|_{L_2}, (\|{\vec{u}}\|_{C_2} + \|{\vec{u}}\|_{L_2})\}\ll V$ Then we have that $|{\vec{u}_{\scriptscriptstyle V}}|<1$, $|{\vec{u}^0_{\scriptscriptstyle V}}|<1$. (Please see Answer 1.) $V$ is a constant, and hence is not dependent on time. \noindent\textbf{Reader's reply:} I do not understand how $V$ can be constant and control the varying values of $\vec{u}$. Thus, the revised version is still totally incorrect. As a matter of fact, it is well known that generically one cannot have a better decay than $O(|x|^{-4})$ for the solutions, so that any attempt to prove a fixed-point theorem in the Schwartz class for the Navier-Stokes equations can be but a disastrous failure. \medskip Additional references:\begin{itemize} \item[{[1]}] C. L. Fefferman; \emph{Existence and smoothness of the Novier-Stokes equation}, The Clay Mathematics Institute, Official Problem Description. \item[{[2]}] A. Tsionskiy, M. Tsionskiy; \emph{Existence, uniqueness and smoothness of solution for 3D Navier-Stokes equations with any smooth initial velocity}, \\ arXiv:1201.1609v8, September 1, 2013. \end{itemize} \end{document}