the objective of this investigation is to study three - dimensional ( 3d ) flows of viscous incompressible fluids which are constructed to exhibit extreme growth of enstrophy .
it is motivated by the question whether the solutions to the 3d incompressible navier - stokes system on unbounded or periodic domains corresponding to smooth initial data may develop a singularity in finite time @xcite . by formation of a `` singularity '' we mean
the situation when some norms of the solution corresponding to smooth initial data have become unbounded after a finite time .
this so - called `` blow - up problem '' is one of the key open questions in mathematical fluid mechanics and , in fact , its importance for mathematics in general has been recognized by the clay mathematics institute as one of its `` millennium problems '' @xcite .
questions concerning global - in - time existence of smooth solutions remain open also for a number of other flow models including the 3d euler equations @xcite and some of the `` active scalar '' equations @xcite .    while the blow - up problem is fundamentally a question in mathematical analysis , a lot of computational studies
have been carried out since the mid-90s in order to shed light on the hydrodynamic mechanisms which might lead to singularity formation in finite time .
given that such flows evolving near the edge of regularity involve formation of very small flow structures , these computations typically require the use of state - of - the - art computational resources available at a given time .
the computational studies focused on the possibility of finite - time blow - up in the 3d navier - stokes and/or euler system include @xcite , all of which considered problems defined on domains periodic in all three dimensions .
recent investigations by @xcite focused on the time evolution of vorticity moments and compared it with the predictions derived from analysis based on rigorous bounds .
we also mention the studies by @xcite and @xcite , along with the references found therein , in which various complexified forms of the euler equation were investigated .
the idea of this approach is that , since the solutions to complexified equations have singularities in the complex plane , singularity formation in the real - valued problem is manifested by the collapse of the complex - plane singularities onto the real axis .
overall , the outcome of these investigations is rather inconclusive : while for the navier - stokes flows most of recent computations do not offer support for finite - time blow - up , the evidence appears split in the case of the euler system . in particular , the recent studies by @xcite and @xcite hinted at the possibility of singularity formation in finite time . in this connection we also mention the recent investigations by @xcite in which blow - up was observed in axisymmetric euler flows in a bounded ( tubular ) domain .
a common feature of all of the aforementioned investigations was that the initial data for the navier - stokes or euler system was chosen in an ad - hoc manner , based on some heuristic arguments .
on the other hand , in the present study we pursue a fundamentally different approach , proposed originally by @xcite and employed also by @xcite for a range of related problems , in which the initial data leading to the most singular behaviour is sought systematically via solution of a suitable variational optimization problem .
we carefully analyze the time evolution induced by the extreme vortex states first identified by @xcite and compare it to the time evolution corresponding to a number of other candidate initial conditions considered in the literature @xcite .
we demonstrate that the taylor - green vortex , studied in the context of the blow - up problem by @xcite , is in fact a particular member of the family of extreme vortex states maximizing the instantaneous rate of enstrophy production in the limit of vanishing enstrophy .
in addition , based on these findings , we identify the set of initial data , parameterized by its energy and enstrophy , for which one can a priori guarantee global - in - time existence of smooth solutions .
this result therefore offers a physically appealing interpretation of an `` abstract '' mathematical theorem concerning global existence of classical solutions corresponding to `` small '' initial data @xcite .
we also emphasize that , in order to establish a direct link with the results of the mathematical analysis discussed below , in our investigation we follow a rather different strategy than in most of the studies referenced above . while these earlier studies relied on data from a relatively small number of simulations performed at a high ( at the given time ) resolution , in the present investigation we explore a broad range of cases , each of which is however computed at a more moderate resolution ( or , equivalently , reynolds number ) . with such an approach to the use of available computational resources , we are able to reveal trends resulting from the variation of parameters which otherwise would be hard to detect .
systematic computations conducted in this way thus allow us to probe the sharpness of the mathematical analysis relevant to the problem .
the question of regularity of solution to the navier - stokes system is usually addressed using `` energy '' methods which rely on finding upper bounds ( with respect to time ) on certain quantities of interest , typically taken as suitable sobolev norms of the solution .
a key intermediate step is obtaining bounds on the rate of growth of the quantity of interest , a problem which can be studied with ode methods . while for the navier - stokes system
different norms of the velocity gradient or vorticity can be used to study the regularity of solutions , the use of enstrophy @xmath5 ( see equation ) is privileged by the well - known result of @xcite , where it was established that if the uniform bound @xmath6 holds , then the regularity of the solution @xmath7 is guaranteed up to time @xmath8 ( to be precise , the solution remains in a suitable gevrey class ) . from the computational point of view , the enstrophy @xmath9 is thus a convenient indicator of the regularity of solutions , because in the light of , singularity formation must manifest itself by the enstrophy becoming infinite .
while characterization of the maximum possible finite - time growth of enstrophy in the 3d navier - stokes flows is the ultimate objective of this research program , analogous questions can also be posed in the context of more tractable problems involving the one - dimensional ( 1d ) burgers equation and the two - dimensional ( 2d ) navier - stokes equation .
although global - in - time existence of the classical ( smooth ) solutions is well known for both these problems @xcite , questions concerning the sharpness of the corresponding estimates for the instantaneous and finite - time growth of various quantities are relevant , because these estimates are obtained using essentially the same methods as employed to derive their 3d counterparts .
since in 2d flows on unbounded or periodic domains the enstrophy may not increase ( @xmath10 ) , the relevant quantity in this case is the palinstrophy @xmath11 , where @xmath12 is the vorticity ( which reduces to a pseudo - scalar in 2d ) .
different questions concerning sharpness of estimates addressed in our research program are summarized together with the results obtained to date in table [ tab : estimates ] .
we remark that the best finite - time estimate for the 1d burgers equation was found _ not _ to be sharp using the initial data obtained from both the instantaneous and the finite - time variational optimization problems @xcite . on the other hand , in 2d the bounds on both the _ instantaneous _ and _ finite - time _ growth of palinstrophy were found to be sharp and , somewhat surprisingly , both estimates were realized by the same family of incompressible vector fields parameterized by energy @xmath13 and palinstrophy @xmath14 , obtained as the solution of an _ instantaneous _ optimization problem @xcite .
it is worth mentioning that while the estimate for the instantaneous rate of growth of palinstrophy @xmath15 ( see table [ tab : estimates ] ) was found to be sharp with respect to variations in palinstrophy , the estimate is in fact not sharp with respect to the prefactor @xmath16 @xcite , with the correct prefactor being of the form @xmath17 .
we add that what distinguishes the 2d problem , in regard to both the instantaneous and finite - time bounds , is that the rhs of these bounds are expressed in terms of two quantities , namely , energy @xmath13 and enstrophy @xmath5 , in contrast to the enstrophy alone appearing in the 1d and 3d estimates . as a result ,
the 2d instantaneous optimization problem had to be solved subject to _ two _ constraints .    in the present investigation we advance the research program summarized in table [ tab : estimates ] by assessing to what extent the finite - time growth of enstrophy predicted by the analytic estimates and can be actually realized by flow evolution starting from different initial conditions , including the extreme vortex states found by @xcite to saturate the instantaneous estimate .
the key finding is that , at least for the range of modest enstrophy values we considered , the growth of enstrophy corresponding to this initial data , which has the form of two colliding axisymmetric vortex rings , is rapidly depleted and there is no indication of singularity formation in finite time .
thus , should finite - time singularity be possible in the navier - stokes system , it is unlikely to result from initial conditions instantaneously maximizing the rate of growth of enstrophy .
we also provide a comprehensive characterization of the extreme vortex states which realize estimate together with the resulting flow evolutions .
the structure of the paper is as follows : in the next section we present analytic estimates on the instantaneous and finite - time growth of enstrophy in 3d flows . in  [ sec:3d_instopt ]
we formulate the variational optimization problems which will be solved to find the vortex states with the largest rate of enstrophy production and in
[ sec:3d_instopt_e0to0 ] we provide an asymptotic representation for these optimal states in the limit of vanishing enstrophy . in  [ sec:3d_instopt_e ] we present numerically computed extreme vortex states corresponding to intermediate and large enstrophy values , while in  [ sec : timeevolution ] we analyze the temporal evolution corresponding to different initial data in order to compare it with the predictions of estimates and .
our findings are discussed in  [ sec : discuss ] , whereas conclusions and outlook are deferred to
[ sec : final ] .
l|c|c &    estimate +    &    realizability     + & & + & & + & @xmath18 + @xmath19 & +    2d navier - stokes + finite - time    & @xmath20 + @xmath21 & + & & +    3d navier - stokes + finite - time    & @xmath22 + @xmath23 $ ] & +
we consider the incompressible navier - stokes system defined on the 3d unit cube @xmath24 ^ 3 $ ] with periodic boundary conditions    [ eq : nse3d ] @xmath25    where the vector @xmath26 $ ] is the velocity field , @xmath27 is the pressure and @xmath28 is the coefficient of kinematic viscosity ( hereafter we will set @xmath29 which is the same value as used in the seminal study by @xcite ) . the velocity gradient @xmath30 is the tensor with components @xmath31_{ij } =    \partial_j u_i$ ] , @xmath32 .
the fluid density @xmath33 is assumed to be constant and equal to unity ( @xmath34 ) .
the relevant properties of solutions to system can be studied using energy methods , with the energy @xmath35 and its rate of growth given by @xmath36 where `` @xmath37 '' means `` equal to by definition '' .
the enstrophy @xmath38 and its rate of growth are given by @xmath39 for incompressible flows with periodic boundary conditions we also have the following identity @xcite @xmath40 hence , combining  , the energy and enstrophy satisfy the system of ordinary differential equations    @xmath41    a standard approach at this point is to try to upper - bound @xmath2 and using standard techniques of functional analysis it is possible to obtain the following well - known estimate in terms of @xmath13 and @xmath5 @xcite @xmath42 for @xmath43 an absolute constant .
a related estimate expressed entirely in terms of the enstrophy @xmath5 is given by @xmath44 by simply integrating the differential inequality in with respect to time we obtain the finite - time bound @xmath45 which clearly becomes infinite at time @xmath46 $ ] .
thus , based on estimate , it is not possible to establish the boundedness of the enstrophy @xmath47 globally in time and hence the regularity of solutions .
therefore , the question about the finite - time singularity formation can be recast in terms of whether or not estimate can be saturated . by this
we mean the existence of initial data with enstrophy @xmath48 such that the resulting time evolution realizes the largest growth of enstrophy @xmath47 allowed by the right - hand side ( rhs ) of estimate .
a systematic search for such most singular initial data using variational optimization methods is the key theme of this study .
although different notions of sharpness of an estimate can be defined , e.g. , sharpness with respect to constants or exponents in the case of estimates in the form of power laws , the precise notion of sharpness considered in this study is the following    [ def : notionsharpness ] given a parameter @xmath49 and maps @xmath50 , the estimate @xmath51 is declared sharp in the limit @xmath52 if and only if @xmath53    from this definition , the sharpness of estimates in the form @xmath54 for some @xmath55 and @xmath56 can be addressed in the limit @xmath57 by studying the adequacy of the exponent @xmath58 .    the question of sharpness of estimate was addressed in the seminal study by @xcite , see also @xcite , who constructed a family of divergence - free velocity fields saturating this estimate .
more precisely , these vector fields were parameterized by their enstrophy and for sufficiently large values of @xmath5 the corresponding rate of growth @xmath1 was found to be proportional to @xmath59 .
therefore , in agreement with definition [ def : notionsharpness ] , estimate was declared sharp up to a numerical prefactor .
however , the sharpness of the instantaneous estimate alone does not allow us to conclude about the possibility of singularity formation , because for this situation to occur , a sufficiently large enstrophy growth rate would need to be sustained over a _ finite _ time window @xmath60 .
in fact , assuming the instantaneous rate of growth of enstrophy in the form @xmath61 for some @xmath62 , any exponent @xmath63 will produce blow - up of @xmath47 in finite time if the rate of growth is sustained .
the fact that there is no blow - up for @xmath64 follows from grnwall s lemma and the fact that one factor of @xmath5 in can be bounded in terms of the initial energy using as follows @xmath65 \leq \frac{1}{2\nu } \k(0 ) .
\label{eq : kt}\ ] ] this relation also leads to an alternative form of the estimate for the finite - time growth of enstrophy , namely @xmath66 \label{eq : evs_t_fixe}\end{aligned}\ ] ] which is more convenient than from the computational point of view and will be used in the present study .
we note , however , that since the rhs of this inequality can not be expressed entirely in terms of properties of the initial data , this is _ not _ in fact an a priori estimate .
estimate also allows us to obtain a condition on the size of the initial data , given in terms of its energy @xmath67 and enstrophy @xmath68 , which guarantees that smooth solutions will exist globally in time , namely , @xmath69 from which it follows that @xmath70 thus , flows with energy and enstrophy satisfying inequality are guaranteed to be smooth for all time , in agreement with the regularity results available under the assumption of small initial data @xcite .
sharpness of instantaneous estimate , in the sense of definition [ def : notionsharpness ] , can be probed by constructing a family of `` extreme vortex states '' @xmath71 which , for each @xmath72 , have prescribed enstrophy @xmath73 and produce the largest possible rate of growth of enstrophy @xmath74 . given the form of , the fields @xmath71 can be expected to exhibit ( at least piecewise ) smooth dependence on @xmath0 and we will refer to the mapping @xmath75 as a `` maximizing branch '' .
thus , information about the sharpness of estimate can be deduced by analyzing the relation @xmath0 versus @xmath74 obtained for a possibly broad range of enstrophy values .
a maximizing branch is constructed by finding , for different values of @xmath0 , the extreme vortex states @xmath71 as solutions of a variational optimization problem defined below .
hereafter , @xmath76 will denote the sobolev space of functions with square - integrable second derivatives endowed with the inner product @xcite @xmath77 where @xmath78 are parameters with the meaning of length scales ( the reasons for introducing these parameters in the definition of the inner product will become clear below ) .
the inner product in the space @xmath79 is obtained from by setting @xmath80 .
the notation @xmath81 will refer to the sobolev space @xmath76 of functions with zero mean . for every fixed value @xmath0 of enstrophy we will look for a divergence - free vector field @xmath71 maximizing the objective function @xmath82 defined in .
we thus have the following    [ pb : maxdedt_e ] given @xmath83 and the objective functional @xmath84 from equation , find @xmath85    which will be solved for enstrophy @xmath0 spanning a broad range of values .
this approach was originally proposed and investigated by @xcite . in the present study
we extend and generalize these results by first showing how other fields considered in the context of the blow - up problem for both the euler and navier - stokes system , namely the taylor - green vortex , also arise from variational problem [ pb : maxdedt_e ] .
we then thoroughly analyze the time evolution corresponding to our extreme vortex states and compare it with the predictions of the finite - time estimates and . as discussed at the end of this section , some important aspects of our approach to solving problem [ pb : maxdedt_e ] are also quite different from the method adopted by @xcite .
the smoothness requirement in the statement of problem [ pb : maxdedt_e ] ( @xmath86 ) follows from the definition of the objective functional @xmath84 in equation , where both the viscous term @xmath87 contain derivatives of order up to two .
the constraint manifold @xmath88 can be interpreted as an intersection of the manifold ( a subspace ) @xmath89 of divergence - free fields and the manifold @xmath90 of fields with prescribed enstrophy @xmath0 .
the structure of these constraint manifolds is reflected in the definition of the corresponding projections @xmath91 ( without a subscript , @xmath92 refers to a generic manifold ) which is given for each of the two constraints as follows :    * ( _ div - free_)-constraint : the projection of a field @xmath93 onto the subspace of solenoidal fields @xmath94 is performed using the helmholtz decomposition ; accordingly , every zero - mean vector field @xmath86 can be decomposed uniquely as @xmath95 where @xmath96 and @xmath97 are scalar and vector potentials , respectively ; it follows from the identity @xmath98 , valid for any sufficiently smooth vector field @xmath97 , that the projection @xmath99 is given simply by @xmath100 and is therefore calculated as @xmath101 ,   \label{eq : phodge}\ ] ] where @xmath102 is the inverse laplacian associated with the periodic boundary conditions ; the operator @xmath103 is also known as the leray - helmholtz projector . * @xmath104-constraint : the projection onto the manifold @xmath105 is calculated by the normalization @xmath106    thus , composing with , the projection onto the manifold @xmath88 defined in problem [ pb : maxdedt_e ] is constructed as @xmath107 this approach , which was already successfully employed by @xcite , allows one to enforce the enstrophy constraint essentially with the machine precision .    for a given value of @xmath0
, the maximizer @xmath71 can be found as @xmath108 using the following iterative procedure representing a discretization of a gradient flow projected on @xmath109 @xmath110 where @xmath111 is an approximation of the maximizer obtained at the @xmath112-th iteration , @xmath113 is the initial guess and @xmath114 is the length of the step in the direction of the gradient .
it is ensured that the maximizers @xmath71 obtained for different values of @xmath0 lie on the same maximizing branch by using the continuation approach , where the maximizer @xmath71 is employed as the initial guess @xmath113 to compute @xmath115 at the next enstrophy level for some sufficiently small @xmath116 . as will be demonstrated in  [ sec:3d_instopt_e0to0 ] , in the limit @xmath117 optimization problem [ pb : maxdedt_e ] admits a discrete family of closed - form solutions and each of these vortex states is the limiting ( initial ) member @xmath118 of the corresponding maximizing branch . as such
, these limiting extreme vortex states are used as the initial guesses @xmath113 for the calculation of @xmath119 , i.e. , they serve as `` seeds '' for the calculation of an entire maximizing branch ( as discussed in  [ sec : discuss ] , while there exist alternatives to the continuation approach , this technique in fact results in the fastest convergence of iterations and also ensures that all computed extreme vortex states lie on a single branch ) .
the procedure outlined above is summarized as algorithm [ alg : optimalg ] , whereas all details are presented below .
set @xmath120 set @xmath121 @xmath122 @xmath123 @xmath124 compute @xmath125    compute the @xmath126 gradient @xmath127 , see equation    compute the sobolev gradient @xmath128 , see equation    compute the step size @xmath114 , see equation    set @xmath129 , see equations     set @xmath130    compute the ` relative error ` @xmath131    set @xmath132    set @xmath133    a key step of algorithm [ alg : optimalg ] is the evaluation of the gradient @xmath134 of the objective functional @xmath135 , cf .
, representing its ( infinite - dimensional ) sensitivity to perturbations of the velocity field @xmath93 , and it is essential that the gradient be characterized by the required regularity , namely , @xmath136 .
this is , in fact , guaranteed by the riesz representation theorem @xcite applicable because the gteaux differential @xmath137 , defined as @xmath138 $ ] for some perturbation @xmath139 , is a bounded linear functional on @xmath140 .
the gteaux differential can be computed directly to give @xmath141\,d{\mathbf{x}}-2\nu\int_{\omega}{\delta}^2{\mathbf{u}}\cdot{\mathbf{u}}'\,d{\mathbf{x}}\label{eq : dr}\ ] ] from which , by the riesz representation theorem , we obtain @xmath142 with the riesz representers @xmath134 and @xmath143 being the gradients computed with respect to the @xmath144 and @xmath126 topology , respectively , and the inner products defined in .
we remark that , while the @xmath144 gradient is used exclusively in the actual computations , cf . , the @xmath126 gradient is computed first as an intermediate step . identifying the gteaux differential with the @xmath126 inner product and performing integration by parts yields @xmath145 similarly , identifying the gteaux differential with the @xmath144 inner product , integrating by parts and using ,
we obtain the required @xmath144 gradient @xmath146 as a solution of the following elliptic boundary - value problem @xmath147 \nabla\r = \nabla^{l_2 } \r   \qquad \text{in } \
\omega , \\ & \text{periodic boundary conditions}. \end{aligned } \label{eq : gradrh2}\ ] ] the gradient fields @xmath143 and @xmath134 can be interpreted as infinite - dimensional sensitivities of the objective function @xmath135 , cf .  , with respect to perturbations of the field @xmath93 .
while these two gradients may point towards the same local maximizer , they represent distinct `` directions '' , since they are defined with respect to different topologies ( @xmath126 vs.  @xmath144 ) .
as shown by @xcite , extraction of gradients in spaces of smoother functions such as @xmath76 can be interpreted as low - pass filtering of the @xmath126 gradients with parameters @xmath148 and @xmath149 acting as cut - off length - scales and the choice of their numerical values will be discussed in  [ sec:3d_instopt_e ] .
the step size @xmath114 in algorithm is computed as @xmath150 \right\}\ ] ] which is done using a suitable derivative - free line - search algorithm @xcite .
equation can be interpreted as a modification of a standard line search method where the optimization is performed following an arc ( a geodesic ) lying on the constraint manifold @xmath109 , rather than a straight line .
this approach was already successfully employed to solve similar problems in @xcite .
it ought to be emphasized here that the approach presented above in which the projections  and gradients  are obtained based on the infinite - dimensional ( continuous ) formulation to be discretized only at the final stage is fundamentally different from the method employed in the original study by @xcite in which the optimization problem was solved in a fully discrete setting ( the two approaches are referred to as `` optimize - then - discretize '' and `` discretize - then - optimize '' , respectively , cf .
a practical advantage of the continuous ( `` optimize - then - discretize '' ) formulation used in the present work is that the expressions representing the sensitivity of the objective functional @xmath84 , i.e. the gradients @xmath151 and @xmath146 , are independent of the specific discretization approach chosen to evaluate them .
this should be contrasted with the discrete ( `` discretize - then - optimize '' ) formulation , where a change of the discretization method would require rederivation of the gradient expressions .
in addition , the continuous formulation allows us to strictly enforce the regularity of maximizers required in problem [ pb : maxdedt_e ] .
finally and perhaps most importantly , the continuous formulation of the maximization problem makes it possible to obtain elegant closed - form solutions of the problem in the limit @xmath117 , which is done in  [ sec:3d_instopt_e0to0 ] below .
these analytical solutions will then be used in  [ sec:3d_instopt_e ] to guide the computation of maximizing branches by numerically solving problem [ pb : maxdedt_e ] for a broad range of @xmath0 , as outlined in algorithm [ alg : optimalg ] .
it is possible to find analytic solutions to problem [ pb : maxdedt_e ] in the limit @xmath152 using perturbation methods . to simplify the notation , in this section
we will drop the subscript @xmath0 when referring to the optimal field .
the euler - lagrange system representing the first - order optimality conditions in optimization problem [ pb : maxdedt_e ] is given by @xcite    [ eq : kkt_e ] @xmath153    where @xmath154 and @xmath155 are the lagrange multipliers associated with the constraints defining the manifold @xmath88 , and @xmath156 , given by @xmath157 is the bilinear form from equation . using the formal series expansions with @xmath158    [ eq : series3d ] @xmath159    in and collecting terms proportional to different powers of @xmath160 , it follows from that , at every order @xmath161 in @xmath160 , we have @xmath162 similarly , equation leads to @xmath163 at every order @xmath164 in @xmath160 .
it then follows from equation that @xmath165\e_0^{2\alpha } + \ldots \\                & = & \e_0,\end{aligned}\ ] ] which , for @xmath166 , forces @xmath167 .
hence , @xmath168 , @xmath169 and @xmath170 .
the systems at orders @xmath171 and @xmath172 are given by :    [ eq : maxdedt_asympt_1 ] @xmath173    [ eq : maxdedt_asympt_2 ] @xmath174    where the fact that @xmath175 for all @xmath176 has been used . while continuing this process to larger values of @xmath164 may lead to some interesting insights , for the purpose of this investigation
it is sufficient to truncate expansions at the order @xmath177 .
the corresponding approximation of the objective functional then becomes @xmath178 it is worth noting that , in the light of relation , the maximum rate of growth of enstrophy in the limit of small @xmath0 is in fact negative , meaning that , for sufficiently small @xmath0 , the enstrophy itself is a decreasing function for all times .
this observation is consistent with the small - data regularity result discussed in introduction .    as regards problem defining the triplet @xmath179 , taking the divergence of equation and using the condition @xmath180 leads to the laplace equation @xmath181 in @xmath182 . since for zero - mean functions defined on @xmath182 , @xmath183
, it follows that @xmath184 and equation is reduced to the eigenvalue problem @xmath185 with @xmath186 satisfying the incompressibility condition .
direct calculation using equation and condition leads to an asymptotic expression for the objective functional in the limit of small enstrophy @xmath187 solutions to the eigenvalue problem in equation can be found using the fourier expansion of @xmath186 given as ( with hats denoting fourier coefficients ) @xmath188 where @xmath189 is a set of wavevectors @xmath190 for which @xmath191 .
the eigenvalue problem then becomes @xmath192\widehat{{\mathbf{u}}}_1({\mathbf{k } } ) & = & 0 \qquad\forall\,{\mathbf{k}}\in{\mathcal{w } } , \\
\widehat{{\mathbf{u}}}_1({\mathbf{k}})\cdot{\mathbf{k } } & = & 0 \qquad\forall\,{\mathbf{k}}\in{\mathcal{w}},\end{aligned}\ ] ] with solutions obtained by choosing , for any @xmath193 , a set of wavevectors with the following structure @xmath194 and @xmath195 with an appropriate form satisfying the incompressibility condition @xmath196 . for the solutions to equation constructed in such manner
it then follows that @xmath197 and the optimal asymptotic value of @xmath84 is given by @xmath198    since the fields @xmath186 are real - valued , their fourier modes must satisfy @xmath199 , where @xmath200 denotes the complex conjugate ( c.c . ) of @xmath201 . depending on the choice of @xmath202 ,
a number of different solutions of can be constructed and below we focus on the following three most relevant cases characterized by the largest values of @xmath203 :    1 .
@xmath204 , where @xmath205 , @xmath206 , is the @xmath207 unit vector of the canonical basis of @xmath208 ; the most general solution can then be constructed as @xmath209 with the complex - valued constant vectors @xmath210 $ ] , @xmath211 $ ] and @xmath212 $ ] suitably chosen so that @xmath170 ; hereafter we will use the values @xmath213 ; it follows that @xmath214 @xmath215 , and the optimal asymptotic value of @xmath84 obtained from equation is given by @xmath216 [ c1 ] 2 .
@xmath217 , where @xmath218 denotes the set whose elements are the negatives of the elements of set @xmath219 , for @xmath220 ; the most general solution can be then constructed as @xmath221\cdot{\mathbf{x } } } } +                        \mathbf{b}{\textrm{e}^{2\pi i[1,-1,0]\cdot{\mathbf{x } } } } +                       \mathbf{c}{\textrm{e}^{2\pi i[1,0,1]\cdot{\mathbf{x } } } } + \nonumber \\                 &    & \mathbf{d}{\textrm{e}^{2\pi i[1,0,-1]\cdot{\mathbf{x } } } } +                       \mathbf{e}{\textrm{e}^{2\pi i[0,1,1]\cdot{\mathbf{x } } } } +                        \mathbf{f}{\textrm{e}^{2\pi i[0,1,-1]\cdot{\mathbf{x } } } } + \textrm{c.c . }
\label{eq : uvec_3d_k2}\end{aligned}\ ] ] with the constants @xmath222 suitably chosen so that @xmath223 = 0 $ ] , @xmath224 = 0 , \ldots,\mathbf{f}\cdot[0,1,-1 ] = 0 $ ] , which ensures that incompressibility condition is satisfied , and that @xmath170 ; in this case , @xmath225 , @xmath226 , and the optimal asymptotic value of @xmath84 is @xmath227 [ c2 ] 3 .   @xmath228 for @xmath229 ; the most general solution can then be constructed as @xmath230\cdot{\mathbf{x } } } } +                        \mathbf{b}{\textrm{e}^{2\pi i[-1,1,1]\cdot{\mathbf{x } } } } + \nonumber \\                 &    & \mathbf{c}{\textrm{e}^{2\pi i[1,-1,1]\cdot{\mathbf{x } } } } +                        \mathbf{d}{\textrm{e}^{2\pi i[1,1,-1]\cdot{\mathbf{x } } } } + \textrm{c.c . }
\label{eq : uvec_3d_k3}\end{aligned}\ ] ] with the constants @xmath231 suitably chosen so that @xmath232 = 0 $ ] , @xmath233 = 0 $ ] , @xmath234 = 0 $ ] and @xmath235 = 0 $ ] , which ensures that incompressibility condition is satisfied , and that @xmath170 ; in this case , @xmath236 , @xmath237 , and the optimal asymptotic value of @xmath84 is @xmath238 [ c3 ]    the three constructions of the extremal field @xmath186 given in , and are all defined up to arbitrary shifts in all three directions , reflections with respect to different planes and rotations by angles which are multiples of @xmath239 about the different axes . as a result of this nonuniqueness , there is some freedom in choosing the constants @xmath240 .
given that the optimal asymptotic value of @xmath84 depends exclusively on the wavevector magnitude @xmath241 , cf .
, any combination of constants @xmath240 will produce the _ same _ optimal rate of growth of enstrophy .
thus , to fix attention , in our analysis we will set @xmath242 in case ( i ) , @xmath243 in case ( [ c2 ] ) and @xmath244 in case ( [ c3 ] ) . with these choices ,
the contribution from each component of the field @xmath186 to the total enstrophy is the same .
the maximum ( i.e. , least negative ) value of @xmath84 can be thus obtained by choosing the smallest possible @xmath245 .
this maximum is achieved in case ( [ c1 ] ) with the wavevectors @xmath246 $ ] , @xmath247 $ ] , @xmath248 $ ] , and @xmath249 , @xmath250 and @xmath251 , for which @xmath252 .
because of this maximization property , this is the field we will focus on in our analysis in  [ sec:3d_instopt_e ] and  [ sec : timeevolution ] .
the three fields constructed in , and are visualized in figure [ fig : maxdedt_vortexcells ] .
this analysis is performed using the level sets @xmath253 defined as @xmath254 for a suitable function @xmath255 . in figures [ fig : maxdedt_vortexcells](a  c ) we choose @xmath256 with @xmath257 . to complement this information , in figures
[ fig : maxdedt_vortexcells](d  f ) we also plot the isosurfaces and cross - sectional distributions of the @xmath258 component of the field @xmath186 .
the fields shown in figure [ fig : maxdedt_vortexcells ] reveal interesting patterns involving well - defined `` vortex cells '' .
more specifically , we see that in case ( [ c1 ] ) , given by equation and shown in figures [ fig : maxdedt_vortexcells](a , d ) , the vortex cells are staggered with respect to the orientation of the cubic domain @xmath182 in all three planes , whereas in case ( [ c3 ] ) , given by equation and shown in figures [ fig : maxdedt_vortexcells](c , f ) , the vortex cells are aligned with the domain @xmath182 in all three planes . on the other hand , in case ( [ c2 ] ) , given by equation and shown in figures [ fig : maxdedt_vortexcells](b , e ) , the vortex cells are staggered in one plane and aligned in another with the arrangement in the third plane resulting from the arrangement in the first two .
these geometric properties are also reflected in the @xmath258-component of the field @xmath186 which is independent of @xmath258 in cases ( [ c1 ] ) and ( [ c2 ] ) , but exhibits , respectively , a staggered and aligned arrangement of the cells in the @xmath259 plane in these two cases . in case ( [ c3 ] )
the cells exhibit an aligned arrangement in all three planes .
the geometric properties of the extreme vortex states obtained in the limit @xmath117 are summarized in table [ tab : e0 ] .
we remark that an analogous structure of the optimal fields , featuring aligned and staggered arrangements of vortex cells in the limiting case , was also discovered by @xcite in their study of the maximum palinstrophy growth in 2d .
while due to a smaller spatial dimension only two optimal solutions were found in that study , the one characterized by the staggered arrangement also lead to a larger ( less negative ) rate of palinstrophy production .
l|c|c|c|c    case +    &    formula for the velocity field    &    arrangement of cells in @xmath259 plane    &    dependence of @xmath258 component of @xmath186 on @xmath258    &    remarks     +    ( [ c1 ] ) +    &    &    staggered    &    uniform    &    staggered abc flow     +    ( [ c2 ] ) +    &    &    aligned    &    uniform    &    aligned abc flow     +    ( [ c3 ] ) +    &    &    aligned    &    cell - like    &    taylor - green vortex     +    it is also worth mentioning that the initial data for two well - known flows , namely , the arnold - beltrami - childress ( abc ) flow @xcite and the taylor - green flow @xcite , are in fact particular instances of the optimal field @xmath186 corresponding to , respectively , cases ( [ c1 ] ) and ( [ c3 ] ) . following the notation of @xcite ,
general abc flows are characterized by the following velocity field @xmath260 where @xmath261 , @xmath262 and @xmath263 real constants .
the vector field in equation can be obtained from equation by choosing @xmath264 $ ] , @xmath265 $ ] and @xmath266 $ ] . by analogy
, we will refer to the state described by as the `` aligned abc flow '' .
likewise , the well - known taylor - green vortex can be obtained as a particular instance of the field @xmath186 from equation again using a suitable choice of the constants @xmath267 .
traditionally , the velocity field @xmath268 $ ] characterizing the taylor - green vortex is defined as @xcite @xmath269 for @xmath270 . for given values of @xmath271,@xmath272 and @xmath273 in , the corresponding constants @xmath267 in
can be found by separating them into their real and imaginary parts denoted , respectively , @xmath274 and @xmath275 .
then , after choosing @xmath276,\ ] ] the imaginary parts can be determined by solving the following system of linear equations @xmath277 \left [ \begin{array}{c } { \mathbf{a}_{\textrm{im } } } \\ { \mathbf{b}_{\textrm{im } } } \\ { \mathbf{c}_{\textrm{im } } } \\
{ \mathbf{d}_{\textrm{im } } } \end{array } \right ] =   \left [ \begin{array}{c } \mathbf{0 } \\ \gamma_1\mathbf{e}_1 \\ \gamma_2\mathbf{e}_2 \\ \gamma_3\mathbf{e}_3 \end{array } \right],\ ] ] where @xmath278 is the @xmath279 identity matrix .
the values of @xmath280 are thus given by @xmath281,\ ,   { \mathbf{b}_{\textrm{im } } } = -\frac{1}{8}\left [ \begin{array}{c}-\gamma_1 \\ \gamma_2 \\
\gamma_3 \end{array } \right],\ , { \mathbf{c}_{\textrm{im } } } = -\frac{1}{8}\left [ \begin{array}{c}\gamma_1 \\ -\gamma_2 \\
\gamma_3 \end{array } \right],\ , { \mathbf{d}_{\textrm{im } } } = -\frac{1}{8}\left [ \begin{array}{c}\gamma_1 \\ \gamma_2 \\ -\gamma_3 \end{array } \right].\ ] ] a typical choice of the parameters used in the numerical studies performed by @xcite and @xcite is @xmath282 and @xmath283 .
we remark that the taylor - green vortex has been employed as the initial data in a number of studies aimed at triggering singular behaviour in both the euler and navier - stokes systems @xcite .
it is therefore interesting to note that it arises in the limit @xmath117 as one of the extreme vortex states in the variational formulation considered in the present study .
it should be emphasized , however , that out of the three optimal states identified above ( see table [ tab : e0 ] ) , the taylor - green vortex is characterized by the smallest ( i.e. , the most negative ) instantaneous rate of enstrophy production @xmath1 . on the other hand ,
we are not aware of any prior studies involving abc flows in the context of extreme behaviour and potential singularity formation .
the time evolution corresponding to these states and some other initial data will be analyzed in detail in
[ sec : timeevolution ] .
in this section we analyze the optimal vortex states @xmath71 obtained for finite values of the enstrophy in which we extend the results obtained in the seminal study by @xcite . as was also the case in the analogous study in 2d @xcite
, there is a distinct branch of extreme states @xmath71 parameterized by the enstrophy @xmath0 and corresponding to each of the three limiting states discussed in  [ sec:3d_instopt_e0to0 ] ( cf .
figure [ fig : maxdedt_vortexcells ] and table [ tab : e0 ] ) .
each of these branches is computed using the continuation approach outlined in algorithm [ alg : optimalg ] . as a key element of the gradient - based maximization technique , the gradient expressions
are approximated pseudo - spectrally using standard dealiasing of the nonlinear terms and with resolutions varying from @xmath284 in the low - enstrophy cases to @xmath285 in the high - enstrophy cases , which necessitated a massively parallel implementation using the message passing interface ( mpi ) . as regards the computation of the sobolev @xmath144 gradients , cf .  , we set @xmath286 , whereas the second parameter @xmath149 was adjusted during the optimization iterations and was chosen so that @xmath287 $ ] , where @xmath288 is the length scale associated with the spatial resolution @xmath289 used for computations and @xmath290 is the characteristic length scale of the domain @xmath182 , that is , @xmath291 and @xmath292 .
we remark that , given the equivalence of the inner products corresponding to different values of @xmath148 and @xmath149 ( as long as @xmath293 ) , these choices do not affect the maximizers found , but only how rapidly they are approached by iterations . for further details concerning the computational approach
we refer the reader to the dissertation by @xcite .
as was the case in the analogous 2d problem studied by @xcite , the largest instantaneous growth of enstrophy is produced by the states with vortex cells staggered in all planes , cf .
case ( [ c1 ] ) in table [ tab : e0 ] .
therefore , in our analysis we will focus exclusively on this branch of extreme vortex states which has been computed for @xmath294 $ ] .    the optimal instantaneous rate of growth of enstrophy @xmath295 and the energy of the optimal states @xmath296 are shown as functions of @xmath0 for small @xmath0 in figures [ fig : rvse0_fixe_small](a ) and [ fig : rvse0_fixe_small](b ) , respectively .
as indicated by the asymptotic form of @xmath84 in and the poincar limit @xmath297 , both of which are marked in these figures , the behaviour of @xmath298 and @xmath296 as @xmath117 is correctly captured by the numerically computed optimal states . in particular , we note that @xmath298 is negative for @xmath299 and exhibits the same trend as predicted in for @xmath117 . for larger values of @xmath0
the rate of growth of enstrophy becomes positive .
likewise , the asymptotic behaviour of the energy of the optimal fields does not come as a surprise since , as discussed in
[ sec:3d_instopt_e0to0 ] , in the limit @xmath152 the maximizers of @xmath84 are eigenfunctions of the laplacian , which also happen to saturate poincar s inequality .
the structure of the optimal vortex states @xmath71 is analyzed next .
they are visualized using in which the vortex cores are identified as regions @xmath300 for @xmath301 defined as @xcite @xmath302,\ ] ] where @xmath303 and @xmath304 are the symmetric and anti - symmetric parts of the velocity gradient tensor @xmath305 , that is , @xmath306_{ij } = \frac{1}{2}(\partial_j u_i + \partial_i u_j)$ ] and @xmath307_{ij } = \frac{1}{2}(\partial_j u_i - \partial_i u_j ) $ ] , @xmath32 .
the quantity @xmath301 can be interpreted as the local balance between the strain rate and the vorticity magnitude .
the isosurfaces @xmath308 representing the optimal states @xmath71 with selected values of @xmath0 are shown in figures [ fig : rvse0_fixe_small](c)-(e ) .
for the smallest values of @xmath0 , the optimal fields exhibit a cellular structure already observed in figure [ fig : maxdedt_vortexcells](a ) . for increasing values of @xmath0
this cellular structure transforms into a vortex ring , as seen in figure [ fig : rvse0_fixe_small](e ) .
the component of vorticity normal to the plane @xmath309 for @xmath310 $ ] and @xmath311 $ ] is shown in figures [ fig : rvse0_fixe_small](f)-(h ) , where the transition from cellular structures to a localized vortex structure as enstrophy increases is evident .
the results corresponding to large values of @xmath0 are shown in figure [ fig : rvse0_fixe_large ] with the maximum rate of growth of enstrophy @xmath298 plotted as a function of @xmath0 in figure [ fig : rvse0_fixe_large](a ) .
we observe that , as @xmath0 increases , this relation approaches a power law of the form @xmath312 . in order to determine the prefactor @xmath313 and the exponent @xmath314 we perform a local least - squares fit of the power law to the actual relation @xmath298 versus @xmath0 for increasing values of @xmath0 starting with @xmath315 ( this particular choice the starting value is justified below ) .
then , the exponent @xmath314 is computed as the average of the exponents obtained from the local fits with their standard deviation providing the error bars , so that we obtain @xmath316 ( the same approach is also used to determine the exponents in other power - law relations detected in this study ) .
we note that the exponent @xmath314 obtained in is in fact very close to 3 which is the exponent in estimate . for the value of the viscosity coefficient used in the computations ( @xmath29 ) ,
the constant factor @xmath317 in estimate has the value @xmath318 which is approximately seven orders of magnitude larger than @xmath319 given in . to shed more light at the source of this discrepancy
, the objective functional @xmath84 from equation can be separated into a negative - definite viscous part @xmath320 and a cubic part @xmath321 defined as    @xmath322    so that @xmath323 .
the values of @xmath324 are also plotted in figure [ fig : rvse0_fixe_large](a ) and it is observed that this quantity exhibits a power - law behaviour of the form @xmath325 while the value of @xmath326 is slightly larger than the value of @xmath319 in , it is still some six orders of magnitude smaller than the constant factor @xmath317 from estimate .
these differences notwithstanding , we may conclude that estimate is sharp in the sense of definition [ def : notionsharpness ] . the power laws from equations and are consistent with the results first presented by @xcite , where the authors reported a power - law with exponent @xmath327 and a constant of proportionality @xmath328 .
the energy of the optimal fields @xmath296 for large values of @xmath0 is shown in figure [ fig : rvse0_fixe_large](b ) in which we observe that the energy stops to increase at about @xmath329 .
this transition justifies using @xmath330 as the lower bound on the range of @xmath0 where the power laws are determined via least - square fits .
figures [ fig : rvse0_fixe_large](c)-(e ) show the isosurfaces @xmath308 representing the optimal fields @xmath71 for selected large values of @xmath0 .
the formation of these localized vortex structures featuring two rings as @xmath0 increases is evident in these figures .
the formation process of localized vortex structures is also visible in figures [ fig : rvse0_fixe_large](f)(h ) , where the component of vorticity normal to the plane @xmath331 for @xmath332 $ ] and @xmath333 $ ] is shown ( we note that the planes used in figures [ fig : rvse0_fixe_small](c)(e ) and [ fig : rvse0_fixe_large](c)(e ) have different orientations ) .
+   +     +   +    next we examine the variation of different diagnostics applied to the extreme states @xmath71 as enstrophy @xmath0 increases .
the maximum velocity @xmath334 and maximum vorticity @xmath335 of the optimal fields are shown , respectively , in figures [ fig : scalinglaws_fixe](a ) and [ fig : scalinglaws_fixe](b ) as functions of @xmath0 .
for each quantity , two distinct power laws are observed in the forms    @xmath336    and    @xmath337    in order to quantify the variation of the relative size of the vortex structures , we will introduce two characteristic length scales . the first one is based on the energy and enstrophy , and was defined by @xcite as @xmath338^{1/2}.\ ] ] it is therefore equivalent to the taylor microscale @xmath339 used in turbulence research @xcite .
another length scale , better suited to the ring - like vortex structures shown in figures [ fig : rvse0_fixe_large](c)-(e ) , is the average radius @xmath340 of one of the vortex rings calculated as @xmath341 and @xmath342 is the characteristic function of the set @xmath343 , \ { \mathbf{x}}_0 = [ 1/2,1/2,1/2 ] \}.\end{aligned}\ ] ] in the above definition of the set @xmath344 , the intersection of the two regions is necessary to restrict the set @xmath340 to only one of the two ring structures visible in figures [ fig : rvse0_fixe_large](c)(e ) . the quantity @xmath345 can be therefore interpreted as the geometric centre of one of the vortex rings .
the dependence of @xmath346 and @xmath347 on @xmath0 is shown in figures [ fig : scalinglaws_fixe](c , d ) in which the following power laws can be observed    @xmath348    by comparing the error bars in the key power laws and with the error bars in power - law relations , , and , we observe that there is less uncertainty in the first case , indicating that the quantities @xmath334 , @xmath335 , @xmath346 and @xmath347 tend to be more sensitive to approximation errors than @xmath349 .
non - negligible error bars may also indicate that , due to modest enstrophy values attained in our computations , the ultimate asymptotic regime corresponding to @xmath350 has not been reached in some power laws .
a useful aspect of employing the average ring radius @xmath340 as the characteristic length scale is that its observed scaling with respect to @xmath0 can be used as an approximate indicator of the resolution @xmath351 required to numerically solve problem [ pb : maxdedt_e ] for large values of enstrophy . from the scaling in relation , it is evident that a two - fold increase in the value of @xmath0 will be accompanied by a similar reduction in @xmath347 , thus requiring an eight - fold increase in the resolution ( a two - fold increase in each dimension ) .
this is one of the reasons why computation of extreme vortex states @xmath71 for large enstrophy values is a very challenging computational task . in particular
, this relation puts a limit on the largest value of @xmath0 for which problem [ pb : maxdedt_e ] can be in principle solved computationally at the present moment : a value of @xmath352 , a mere order of magnitude above the largest value of @xmath0 reported here , would require a resolution of @xmath353 used by some of the largest navier - stokes simulations to date .
to summarize , as the enstrophy increases from @xmath354 to @xmath355 , the optimal vortex states change their structure from cellular to ring - like . while with the exception of @xmath74 and @xmath296 , all of the diagnostic quantities behave in a monotonous manner , the corresponding power laws change at about @xmath356 , which approximately marks the transition from the cellular to the ring - like structure ( cf .
figure [ fig : rvse0_fixe_small](e ) vs.  [ fig : rvse0_fixe_large](c ) ) .
this is also the value of the enstrophy beyond which the energy @xmath296 starts to decrease ( figure [ fig : rvse0_fixe_large](b ) ) .
this transition also coincides with a change of the symmetry properties of the extreme vortex states @xmath71  while in the limit @xmath117 these fields feature reflection and discrete rotation symmetries ( cf .
[ sec:3d_instopt_e0to0 ] ) , for @xmath357 the optimal states are characterized by axial symmetry .
the asymptotic ( as @xmath358 ) extreme vortex states on locally maximizing branches corresponding to the aligned abc flow and the taylor - green vortex ( cf .
table [ tab : e0 ] ) are similar to the fields shown in figures [ fig : rvse0_fixe_large](c)(h ) , except for a different orientation of their symmetry axes with respect to the periodic domain @xmath182 ( these results are not shown here for brevity ) .
the different power laws found here are compared to the corresponding results obtained in 2d in  [ sec : discuss ] .
it is also worth mentioning that , as shown by @xcite , all power laws discussed in this section , cf .
, , , and , can be deduced rigorously using arguments based on dimensional analysis under the assumption of axisymmetry for the optimal fields @xmath71 .
+    finally , the findings of this section allow us to shed some light on the `` small data '' result which provides the conditions on the size of the initial data @xmath359 , given in terms of its energy @xmath67 and enstrophy @xmath68 , in the navier - stokes system guaranteeing that smooth solutions exist globally in time .
the power - law fits and allow us to sharpen condition be replacing the constant on the rhs with either @xmath360 or @xmath361 , so that we obtain @xmath362 the region of the `` phase space '' @xmath363 described by condition is shown in white in figure [ fig : k0e0 ] . the gray region represents the values of @xmath67 and @xmath68 for which long - time existence of smooth solutions can not be a priori guaranteed ( the two shades of gray correspond to the two constants which can be used in ) .
solid circles represent the different extreme states found in this section , whereas the thin curves mark the time - dependent trajectories which will be analyzed in  [ sec : timeevolution ] .
we conclude from figure [ fig : k0e0 ] that the change of the properties of the optimal states @xmath71 discussed above occurs in fact for the values of enstrophy ( @xmath364 ) for which the states @xmath71 are on the boundary of the region of guaranteed long - time regularity .
the goal of this section is to analyze the time evolution , governed by the navier - stokes system , with extreme vortex states identified in  [ sec:3d_instopt_e ] used as the initial data @xmath359 . in particular , we are interested in the finite - time growth of enstrophy @xmath47 and how it relates to estimates , and .
we will compare these results with the growth of enstrophy obtained using other types of initial data which have also been studied in the context of the blow - up problem for both the euler and navier - stokes systems , namely , the taylor - green vortex @xcite , the kida - pelz flow @xcite , colliding lamb - chaplygin dipoles @xcite and perturbed antiparallel vortex tubes @xcite .
precise characterization of these different initial conditions is provided in table [ tab : initialconditions ] and , for the sake of completeness , the last three states are also visualized in figure [ fig : summaryic ] .
we comment that , with the exception of the taylor - green vortex which was shown in  [ sec:3d_instopt_e0to0 ] to be a local maximizer of problem [ pb : maxdedt_e ] in the limit @xmath117 , all these initial conditions were postulated based on rather ad - hoc physical arguments .
we also add that , in order to ensure a fair comparison , the different initial conditions listed in table [ tab : initialconditions ] are rescaled to have the same enstrophy @xmath0 , which is different from the enstrophy values used in the original studies where these initial conditions were investigated @xcite .
as regards our choices of the initial enstrophy @xmath0 , to illustrate different possible behaviours , we will consider initial data located in the two distinct regions of the phase space @xmath363 shown in figure [ fig : k0e0 ] , corresponding to values of @xmath365 and @xmath0 for which global regularity may or may not be a priori guaranteed according to estimates .    .characterization of the different initial data used in time evolution studies in  [ sec : timeevolution ] . [ cols="^,^,^",options="header " , ]     system is solved numerically with an approach combining a pseudo - spectral approximation of spatial derivatives with a third - order semi - implicit runge - kutta method @xcite used to discretize the problem in time . in the evaluation of the nonlinear term dealiasing was used based on the @xmath366 rule together with the gaussian filtering proposed by @xcite .
massively parallel implementation based on mpi and using the fftw routines @xcite to perform fourier transforms allowed us to use resolutions varying from @xmath367 to @xmath368 in the low - enstrophy and high - enstrophy cases , respectively .
a number of different diagnostics were checked to ensure that all flows discussed below are well resolved .
we refer the reader to the dissertation by @xcite for additional details and a validation of this approach .
the time - dependent results will be shown with respect to a normalized time defined as @xmath369 with @xmath370 and @xmath371 ( cf .
equation ) playing the roles of the characteristic velocity and length scale .
we begin by showing the time evolution of the enstrophy @xmath372 corresponding to the five different initial conditions listed in table [ tab : initialconditions ] with @xmath373 and @xmath374 in figures [ fig : fvs_t_fixe](a ) and [ fig : fvs_t_fixe](b ) , respectively ( because of the faster time - scale , the time axis in the latter figure is scaled logarithmically ) .
we see that the maximizers @xmath71 of problem [ pb : maxdedt_e ] are the only initial data which triggers growth of enstrophy for these values of the initial enstrophy and , as expected , this growth is larger for @xmath374 than for @xmath373 .
the other initial condition which exhibits some tendency for growth when @xmath374 is the taylor - green vortex . in all cases
the enstrophy eventually decays to zero for large times .
next we examine whether the flow evolutions starting from the instantaneous maximizers @xmath71 as the initial data saturate the finite - time estimate .
we do this by defining functions    [ eq : fg ] @xmath375 \label{eq : fgb}\end{aligned}\ ] ]    representing , respectively , the left- and right - hand side of the estimate and then plotting them with respect to the normalized time @xmath376 , which is done in figures [ fig : fg](a ) and [ fig : fg](b ) for @xmath373 and @xmath374 , respectively . the constant @xmath62 in the definition of @xmath377 is numerically computed from the power - law fit in .
it follows from estimate that @xmath378 pointwise in time .
the hypothetical extreme event of a finite - time blow - up can be represented graphically by an intersection of the graph of @xmath379 with the horizontal line @xmath380 , which is also shown in figures [ fig : fg](a)(b ) .
the behaviour of @xmath377 , representing the upper bound in estimate , is quite distinct in figures [ fig : fg](a ) and [ fig :
fg](b ) reflecting the fact that the initial data @xmath71 in the two cases comes from different regions of the phase diagram in figure [ fig : k0e0 ] . in figure
[ fig : fg](a ) , corresponding to @xmath373 , the upper bound @xmath377 never reaches @xmath381 , in agreement with the fact that the finite - time blow - up is a priori ruled out in this case . on the other hand , in figure [ fig : fg](b ) , corresponding to @xmath374 , the upper bound @xmath377 does intersect @xmath381 implying that , in principle , finite - time blow - up might be possible in this case .
the sharpness of estimate can be assessed by analyzing how closely the behaviour of @xmath379 matches that of @xmath377 . in both figures [ fig : fg](a ) and
[ fig : fg](b ) we observe that for a short period of time @xmath379 exhibits a very similar growth to the upper bound @xmath377 , but then this growth slows down and @xmath379 eventually starts to decrease short of ever approaching the limit @xmath381 .
we further characterize the time evolution by showing the maximum enstrophy increase @xmath382 and the time when the maximum is achieved @xmath383 as functions of @xmath0 in figures [ fig : emax_vse0_fixe](a ) and [ fig : emax_vse0_fixe](b ) , respectively . in both cases approximate
power laws in the form @xmath384 are detected in the limit @xmath358 ( as regards the second result , we remark that @xmath385 is not equivalent to the time until which the enstrophy grows at the sustained rate proportional to @xmath386 , cf .
figure [ fig : fg ] ) . to complete presentation of the results ,
the dependence of the quantities @xmath387\ ] ] on the initial enstrophy @xmath0 is shown in figures [ fig : emax_vse0_fixe](c ) and [ fig : emax_vse0_fixe](d ) , respectively .
it is observed that both quantities approximately exhibit a power - law behaviour of the form @xmath388 .
discussion of these results in the context of the estimates recalled in
[ sec : intro ] is presented in the next section .
in this section we provide some comments about the results reported in   [ sec:3d_instopt_e0to0 ] , [ sec:3d_instopt_e ] and [ sec : timeevolution ] .
first , we need to mention that our gradient - based approach to the solution of optimization problem [ pb : maxdedt_e ] can only yield local maximizers and , due to nonconvexity of the problem , it is not possible to guarantee a priori that the maximizers found are global . to test for the possible presence of branches other than the ones found using the continuation approach described in  [ sec:3d_instopt ] , cf .
algorithm [ alg : optimalg ] , we tried to find new maximizers by initializing the gradient iterations with different initial guesses @xmath113 .
they were constructed as solenoidal vector fields with prescribed regularity and random structure , which was achieved by defining the fourier coefficients of @xmath113 as @xmath389 with the amplitude @xmath390 and the phases @xmath391 chosen as random numbers uniformly distributed in @xmath392 $ ] .
however , in all such tests conducted for @xmath393 the gradient optimization algorithm would always converge to maximizers @xmath71 belonging to one of the branches discussed in
[ sec:3d_instopt_e ] ( modulo possible translations in the physical domain ) . while far from settling this issue definitely , these observations lend some credence to the conjecture that the branch identified in  [ sec:3d_instopt_e ]
corresponds in fact to the global maximizers .
these states appear identical to the maximizers found by @xcite and our search has also yielded two additional branches of locally maximizing fields , although we did not capture the lower branch reported by @xcite . however , since that branch does not appear connected to any state in the limit @xmath117 , we speculate that it might be an artifact of the `` discretize - then - optimize '' formulation used by @xcite , in contrast to the `` optimize - then - discretize '' approach employed in our study which provides a more direct control over the analytic properties of the maximizers .
we add that the structure of the maximizing branches found here is in fact quite similar to what was discovered by @xcite in an analogous problem in 2d .
since the 2d problem is more tractable from the computational point of view , in that case we were able to undertake a much more thorough search for other maximizers which did not however yield any solutions not associated with the main branches .
the results reported in  [ sec:3d_instopt_e ] and  [ sec : timeevolution ] clearly exhibit two distinct behaviours , depending on whether or not global - in - time regularity can be guaranteed a priori based on estimates .
these differences are manifested , for example , in the power laws evident in figures [ fig : scalinglaws_fixe ] and [ fig : emax_vse0_fixe ] , as well as in the different behaviours of the rhs of estimate with respect to time in figures [ fig : fg](a ) and [ fig : fg](b ) .
however , for the initial data for which global - in - time regularity can not be ensured a priori there is no evidence of sustained growth of enstrophy strong enough to signal formation of singularity in finite time .
indeed , in figure [ fig : emax_vse0_fixe](c ) one sees that the quantity @xmath394 behaves as @xmath395 , where @xmath396 , when @xmath0 increases , revealing no tendency to approach @xmath381 which is a necessary precursor of a singular behaviour ( cf .
discussion in
[ sec : timeevolution ] ) . to further illustrate how the rate of growth of enstrophy achieved initially by the maximizers @xmath71
is depleted in time , in figure [ fig : dedte ] we show the flow evolution corresponding to @xmath71 with @xmath374 as a trajectory in the coordinates @xmath397 . from the discussion in introduction we know that in order for the singularity to occur in finite time , the enstrophy must grow at least at a sustained rate @xmath398 for some @xmath63 .
in other words , a `` blow - up trajectory '' will be realized only if the trajectory of the flow , expressed in @xmath399 coordinates , is contained in the region @xmath400 , for some positive constants @xmath401 and @xmath402 .
for the flow corresponding to the instantaneous optimizers @xmath71 , the initial direction of a trajectory in @xmath399 coordinates is determined by the vector @xmath403 $ ] and , for initial conditions @xmath359 satisfying @xmath404 , it follows that @xmath405 since the optimal rate of growth is sustained only over a short interval of time , the trajectory of the flow in the @xmath406 coordinates approaches the region @xmath407 only tangentially following the direction of the lower bound @xmath408 , and remains outside @xmath407 for all subsequent times .
this behaviour is clearly seen in the inset of figure [ fig : dedte ] .
an interpretation of this behaviour can be proposed based on equation from which it is clear that the evolution of the flow energy is closely related to the growth of enstrophy . in particular ,
if the initial energy @xmath67 is not sufficiently large , then its depletion due to the initial growth of enstrophy may render the flow incapable of sustaining this growth over a longer period of time .
this is in fact what seems to be happening in the present problem as evidenced by the data shown in figure [ fig : k0e0 ] .
we remark that , for a prescribed enstrophy @xmath68 , the flow energy can not be increased arbitrarily as it is upper - bounded by poincar s inequality @xmath409 .
this behaviour can also be understood in terms of the geometry of the extreme vortex states @xmath71 .
figure [ fig : ring ] shows a magnification of the pair of vortex rings corresponding to the optimal field @xmath71 with @xmath410 .
it is observed that the vorticity field @xmath411 inside the vortex core has an azimuthal component only which exhibits no variation in the azimuthal direction .
thus , in the limit @xmath412 the vortex ring shrinks with respect to the domain @xmath182 ( cf .
figure [ fig : scalinglaws_fixe](d ) ) and the field @xmath71 ultimately becomes axisymmetric ( i.e. , in this limit boundary effects vanish ) . at the same time , it is known that the 3d navier - stokes problem on an unbounded domain and with axisymmetric initial data is globally well posed @xcite , a results which is a consequence of the celebrated theorem due to @xcite .
with @xmath374 in the coordinates @xmath406 . for comparison , in the inset
the thin line represents the borderline growth at the rate @xmath413,scaledwidth=90.0% ]    .
[ fig : dedte ]     with @xmath374 .
the colour coding of the vortex lines is for identification purposes only.,scaledwidth=60.0% ]    we close this section by comparing the different power laws characterizing the maximizers @xmath71 and the corresponding flow evolutions with the results obtained in analogous studies of extreme behaviour in 1d and 2d ( see also table [ tab : estimates ] ) .
first , we note that the finite - time growth of enstrophy @xmath414 in 3d , cf .
figure [ fig : emax_vse0_fixe](a ) , exhibits the same dependence on the enstrophy @xmath0 of the instantaneously optimal initial data as in 1d , i.e. , is directly proportional to @xmath0 in both cases @xcite .
this is also analogous to the maximum growth of palinstrophy @xmath14 in 2d which was found by @xcite to scale with the palinstrophy @xmath415 of the initial data , when the instantaneously optimal initial condition was computed subject to _ one _ constraint only ( on @xmath415 ) .
when the instantaneously optimal initial data was determined subject to _ two _ constraints , on @xmath365 and @xmath415 , then the maximum finite - time growth of palinstrophy was found to scale with @xmath416 @xcite . on the other hand , the time @xmath385 when the maximum enstrophy is attained , cf .
figure [ fig : emax_vse0_fixe](b ) , scales as @xmath417 , which should be contrasted with the scalings @xmath418 and @xmath419 found in the 1d and 2d cases , respectively .
this implies that the time interval during which the extremal growth of enstrophy is sustained in 3d is shorter than the corresponding intervals in 1d and 2d .
by constructing the initial data to exhibit the most extreme behaviour allowed for by the mathematically rigorous estimates , this study offers a fundamentally different perspective on the problem of searching for potentially singular solutions from most earlier investigations . indeed , while the corresponding flow evolutions did not reveal any evidence for finite - time singularity formation , the initial data obtained by maximizing @xmath1 produced a significantly larger growth of enstrophy in finite time than any other candidate initial conditions ( cf .
table [ tab : initialconditions ] and figure [ fig : fvs_t_fixe ] ) . admittedly , this observation is limited to the initial data with @xmath420 , which corresponds to reynolds numbers @xmath421 lower than the reynolds numbers achieved in other studies concerned with the extreme behaviour in the navier - stokes flows @xcite .
given that the definitions of the reynolds numbers applicable to the various flow configurations considered in these studies were not equivalent , it is rather difficult to make a precise comparison in terms of specific numerical values , but it is clear that the largest reynolds numbers attained in these investigations were at least an order of magnitude higher than used in the present study ; for euler flows such a comparison is obviously not possible at all .
however , from the mathematical point of view , based on estimates  , there is no clear indication that a very large initial enstrophy @xmath0 ( or , equivalently , a high reynolds number ) should be a necessary condition for singularity formation in finite time .
in fact , blow - up can not be a priori ruled out as soon as condition is violated , which happens for all initial data lying on the gray region of the phase space in figure [ fig : k0e0 ] .
we remark that additional results were obtained ( not reported in this paper ) by studying the time evolution corresponding to the optimal initial data @xmath71 , but using smaller values of the viscosity coefficient @xmath422 , thereby artificially increasing the reynolds number at the price of making the initial data suboptimal .
although these attempts did increase the amplification of enstrophy as compared to what was observed in figures [ fig : fg ] and [ fig : emax_vse0_fixe ] , no signature of finite - time singularity formation could be detected either .
our study confirmed the earlier findings of @xcite about the sharpness of the instantaneous estimate .
we also demonstrated that the finite - time estimate is saturated by the flow evolution corresponding to the optimal initial data @xmath71 , but only for short times , cf .
figure [ fig : fg ] , which are not long enough to trigger a singular behaviour .    in  [ sec : discuss ] we speculated that a relatively small initial energy @xmath67 , cf .
figure [ fig : rvse0_fixe_large](b ) , might be the property of the extreme vortex states @xmath71 preventing the resulting flow evolutions from sustaining a significant growth of enstrophy over long times .
on the other hand , in introduction we showed that estimate need not be saturated for blow - up to occur in finite time and , in fact , sustained growth at the rate @xmath423 with any @xmath63 will also produce singularity in finite time .
thus , another strategy to construct initial data which could lead to a more sustained growth of enstrophy in finite time might be to increase its kinetic energy by allowing for a smaller instantaneous rate of growth ( i.e. , with an exponent @xmath424 instead of @xmath425 ) .
this can be achieved by prescribing an additional constraint in the formulation of the variational optimization problem , resulting in    [ pb : maxdedt_ke ] @xmath426    it differs from problem [ pb : maxdedt_e ] in that the maximizers are sought at the intersection of the original constraint manifold @xmath88 and the manifold defined by the condition @xmath427 , where @xmath428 is the prescribed energy . while computation of such maximizers is more complicated , robust techniques for the solution of optimization problems of this type have been developed and were successfully used in the 2d setting by @xcite .
preliminary results obtained in the present setting by solving problem [ pb : maxdedt_ke ] for @xmath429 are indicated in figure [ fig : k0e0 ] , where we see that the flow evolutions starting from @xmath430 do not in fact produce a significant growth of enstrophy either .
an alternative , and arguably more flexible approach , is to formulate this problem in terms of multiobjective optimization @xcite in which the objective function @xmath135 in problem [ pb : maxdedt_e ] would be replaced with @xmath431 where @xmath432 $ ] .
solution of such a multiobjective optimization problem has the form of a `` pareto front '' parameterized by @xmath433 .
clearly , the limits @xmath434 and @xmath435 correspond , respectively , to the extreme vortex states already found in  [ sec:3d_instopt_e0to0 ] and  [ sec:3d_instopt_e ] , and to the poincar limit .
another interesting possibility is to replace the energy @xmath35 with the helicity @xmath436 in the multiobjective formulation , as this might allow one to obtain extreme vortex states with a more complicated topology ( i.e. , a certain degree of `` knottedness '' ) .
we note that all the extreme vortex states found in the present study were `` unknotted '' , i.e. , were characterized by @xmath437 , as the vortex rings were in all cases disjoint ( cf .
figure [ fig : ring ] ) .    finally , another promising possibility to find initial data producing a larger growth of enstrophy is to solve a _
finite - time _ optimization problem of the type already studied by @xcite in the context of the 1d burgers equation , namely    [ pb : maxde ] @xmath438    where @xmath439 is the length of the time interval of interest and @xmath359 the initial data for the navier - stokes system .
in contrast to problems [ pb : maxdedt_e ] and [ pb : maxdedt_ke ] , solution of problem [ pb : maxde ] is more complicated as it involves flow evolution .
it represents therefore a formidable computational task for the 3d navier - stokes system .
however , it does appear within reach given the currently available computational resources and will be studied in the near future .
the authors are indebted to charles doering for many enlightening discussions concerning the research problems studied in this work .
the authors are grateful to nicholas kevlahan for making his parallel navier - stokes solver available , which was used to obtain the results reported in  [ sec : timeevolution ] .
anonymous referees provided many insightful comments which helped us improve this work .
this research was funded through an early researcher award ( era ) and an nserc discovery grant , whereas the computational time was made available by sharcnet under its dedicated resource allocation program .
diego ayala was funded in part by nsf award dms-1515161 at the university of michigan and by the institute for pure and applied mathematics at ucla .