Source author record

Carlo Morosi

Carlo Morosi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
6topics
2close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2015arXiv

Large order Reynolds expansions for the Navier-Stokes equations

We consider the Cauchy problem for the incompressible homogeneous Navier-Stokes (NS) equations on a d-dimensional torus, in the C^infinity formulation described, e.g., in [25]. In [22][25] it was shown how to obtain quantitative estimates on the exact solution of the NS Cauchy problem via the "a posteriori" analysis of an approximate solution; such estimates concern the interval of existence of the exact solution and its distance from the approximate solution. In the present paper we consider an approximate solutions of the NS Cauchy problem having the form u^N(t) = sum_{j=0}^N R^j u_j(t), where R is the "mathematical" Reynolds number (the reciprocal of the kinematic viscosity) and the coefficients u_j(t) are determined stipulating that the NS equations be satisfied up to an error O(R^{N+1}). This subject was already treated in [24], where, as an application, the Reynolds expansion of order N=5 in dimension d=3 was considered for the initial datum of Behr-Necas-Wu (BNW). In the present paper, these results are enriched regarding both the theoretical analysis and the applications. Concerning the theoretical aspect, we refine the approach of [24] following [25] and use the symmetries of the initial datum in building up the expansion. Concerning the applicative aspect we consider two more (d=3) initial data, namely, the vortices of Taylor-Green (TG) and Kida-Murakami (KM); the Reynolds expansions for the BNW, TG and KM data are performed via a Python program, attaining orders between N=12 and N=20. Our a posteriori analysis proves, amongst else, that the solution of the NS equations with anyone of the above three data is global if R is below an explicitly computed critical value. Our critical Reynolds numbers are below the ones characterizing the turbulent regime; however these bounds have a sound theoretical support, are fully quantitative and improve previous results of global existence.

preprint2014arXiv

Smooth solutions of the Euler and Navier-Stokes equations from the a posteriori analysis of approximate solutions

The main result of [C. Morosi and L. Pizzocchero, Nonlinear Analysis, 2012] is presented in a variant, based on a C^infinity formulation of the Cauchy problem; in this approach, the a posteriori analysis of an approximate solution gives a bound on the Sobolev distance of any order between the exact and the approximate solution.

preprint2013arXiv

A posteriori estimates for Euler and Navier-Stokes equations

The first two sections of this work review the framework of [6] for approximate solutions of the incompressible Euler or Navier-Stokes (NS) equations on a torus T^d, in a Sobolev setting. This approach starts from an approximate solution u_a of the Euler/NS Cauchy problem and, analyzing it a posteriori, produces estimates on the interval of existence of the exact solution u and on the distance between u and u_a. The next two sections present an application to the Euler Cauchy problem, where u_a is a Taylor polynomial in the time variable t; a special attention is devoted to the case d=3, with an initial datum for which Behr, Necas and Wu have conjectured a finite time blowup [1]. These sections combine the general approach of [6] with the computer algebra methods developed in [9]; choosing the Behr-Necas-Wu datum, and using for u_a a Taylor polynomial of order 52, a rigorous lower bound is derived on the interval of existence of the exact solution u, and an estimate is obtained for the H^3 Sobolev distance between u(t) and u_a(t).

preprint2013arXiv

On the Reynolds number expansion for the Navier-Stokes equations

In a previous paper of ours [Nonlinear Anal. 2012] we have considered the incompressible Navier-Stokes (NS) equations on a d-dimensional torus T^d, in the functional setting of the Sobolev spaces H^n(T^d) of divergence free, zero mean vector fields (n > d/2+1). In the cited work we have presented a general setting for the a posteriori analysis of approximate solutions of the NS Cauchy problem; given any approximate solution u_a, this allows to infer a lower bound T_c on the time of existence of the exact solution u and to construct a function R_n such that || u(t) - u_a(t) ||_n <= R_n(t) for t in [0,T_c). In certain cases it is T_c = + infinity, so global existence is granted for u. In the present paper the framework of [Nonlinear Anal., 2012] is applied using as an approximate solution an expansion u^N(t) = Sum_{j=0}^N R^j u_j(t), where R is the Reynolds number. This allows, amongst else, to derive the global existence of u when R is below some critical value R_{*} (increasing with N in the examples that we analyze). After a general discussion about the Reynolds expansion and its a posteriori analysis, we consider the expansions of orders N=1,2,5 in dimension d=3, with the initial datum of Behr, Necas and Wu [M2AN, 2001]. Computations of order N=5 yield a quantitative improvement of the results previously obtained for this initial datum in [Nonlinear Anal. 2012], where a Galerkin approximate solution was employed in place of the Reynolds expansion.

preprint2012arXiv

On power series solutions for the Euler equation, and the Behr-Necas-Wu initial datum

We consider the Euler equation for an incompressible fluid on a three dimensional torus, and the construction of its solution as a power series in time. We point out some general facts on this subject, from convergence issues for the power series to the role of symmetries of the initial datum. We then turn the attention to a paper by Behr, Necas and Wu in ESAIM: M2AN 35 (2001) 229-238; here, the authors chose a very simple Fourier polynomial as an initial datum for the Euler equation and analyzed the power series in time for the solution, determining the first 35 terms by computer algebra. Their calculations suggested for the series a finite convergence radius τ_3 in the H^3 Sobolev space, with 0.32 < τ_3 < 0.35; they regarded this as an indication that the solution of the Euler equation blows up. We have repeated the calculations of Behr, Necas and Wu, using again computer algebra; the order has been increased from 35 to 52, using the symmetries of the initial datum to speed up computations. As for τ_3, our results agree with the original computations of Behr, Necas and Wu (yielding in fact to conjecture that 0.32 < τ_3 < 0.33). Moreover, our analysis supports the following conclusions: (a) The finiteness of τ_3 is not at all an indication of a possible blow-up. (b) There is a strong indication that the solution of the Euler equation does not blow up at a time close to τ_3. In fact, the solution is likely to exist, at least, up to a time θ_3 > 0.47. (c) Pade' analysis gives a rather weak indication that the solution might blow up at a later time.

preprint2012arXiv

On the constants in a basic inequality for the Euler and Navier-Stokes equations

We consider the incompressible Euler or Navier-Stokes (NS) equations on a d-dimensional torus T^d; the quadratic term in these equations arises from the bilinear map sending two velocity fields v, w : T^d -> R^d into v . D w, and also involves the Leray projection L onto the space of divergence free vector fields. We derive upper and lower bounds for the constants in some inequalities related to the above quadratic term; these bounds hold, in particular, for the sharp constants K_{n d} = K_n in the basic inequality || L(v . D w)||_n <= K_n || v ||_n || w ||_{n+1}, where n in (d/2, + infinity) and v, w are in the Sobolev spaces H^n, H^{n+1} of zero mean, divergence free vector fields of orders n and n+1, respectively. As examples, the numerical values of our upper and lower bounds are reported for d=3 and some values of n. Some practical motivations are indicated for an accurate analysis of the constants K_n.

preprint2011arXiv

On approximate solutions of the incompressible Euler and Navier-Stokes equations

We consider the incompressible Euler or Navier-Stokes (NS) equations on a torus T^d in the functional setting of the Sobolev spaces H^n(T^d) of divergence free, zero mean vector fields on T^d, for n > d/2+1. We present a general theory of approximate solutions for the Euler/NS Cauchy problem; this allows to infer a lower bound T_c on the time of existence of the exact solution u analyzing a posteriori any approximate solution u_a, and also to construct a function R_n such that || u(t) - u_a(t) ||_n <= R_n(t) for all t in [0,T_c). Both T_c and R_n are determined solving suitable "control inequalities", depending on the error of u_a; the fully quantitative implementation of this scheme depends on some previous estimates of ours on the Euler/NS quadratic nonlinearity [15][16]. To keep in touch with the existing literature on the subject, our results are compared with a setting for approximate Euler/NS solutions proposed in [3]. As a first application of the present framework, we consider the Galerkin approximate solutions of the Euler/NS Cauchy problem, with a specific initial datum considered in [2]: in this case our methods allow, amongst else, to prove global existence for the NS Cauchy problem when the viscosity is above an explicitly given bound.

preprint2010arXiv

On the constants in a Kato inequality for the Euler and Navier-Stokes equations

We continue an analysis, started in [10], of some issues related to the incompressible Euler or Navier-Stokes (NS) equations on a d-dimensional torus T^d. More specifically, we consider the quadratic term in these equations; this arises from the bilinear map (v, w) -> v . D w, where v, w : T^d -> R^d are two velocity fields. We derive upper and lower bounds for the constants in some inequalities related to the above bilinear map; these bounds hold, in particular, for the sharp constants G_{n d} = G_n in the Kato inequality | < v . D w | w >_n | <= G_n || v ||_n || w ||^2_n, where n in (d/2 + 1, + infinity) and v, w are in the Sobolev spaces H^n, H^(n+1) of zero mean, divergence free vector fields of orders n and n+1, respectively. As examples, the numerical values of our upper and lower bounds are reported for d=3 and some values of n. When combined with the results of [10] on another inequality, the results of the present paper can be employed to set up fully quantitative error estimates for the approximate solutions of the Euler/NS equations, or to derive quantitative bounds on the time of existence of the exact solutions with specified initial data; a sketch of this program is given.

preprint2009arXiv

An H^1 setting for the Navier-Stokes equations: quantitative estimates

We consider the incompressible Navier-Stokes (NS) equations on a torus, in the setting of the spaces L^2 and H^1; our approach is based on a general framework for semi- or quasi-linear parabolic equations proposed in the previous work [9]. We present some estimates on the linear semigroup generated by the Laplacian and on the quadratic NS nonlinearity; these are fully quantitative, i.e., all the constants appearing therein are given explicitly. As an application we show that, on a three dimensional torus T^3, the (mild) solution of the NS Cauchy problem is global for each H^1 initial datum u_0 with zero mean, such that || curl u_0 ||_{L^2} <= 0.407; this improves the bound for global existence || curl u_0 ||_{L^2} <= 0.00724, derived recently by Robinson and Sadowski [10]. We announce some future applications, based again on the H^1 framework and on the general scheme of [9].

preprint2009arXiv

New results on multiplication in Sobolev spaces

We consider the Sobolev (Bessel potential) spaces H^ell(R^d, C), and their standard norms || ||_ell (with ell integer or noninteger). We are interested in the unknown sharp constant K_{ell m n d} in the inequality || f g ||_{ell} \leqs K_{ell m n d} || f ||_{m} || g ||_n (f in H^m(R^d, C), g in H^n(R^d, C); 0 <= ell <= m <= n, m + n - ell > d/2); we derive upper and lower bounds K^{+}_{ell m n d}, K^{-}_{ell m n d} for this constant. As examples, we give a table of these bounds for d=1, d=3 and many values of (ell, m, n); here the ratio K^{-}_{ell m n d}/K^{+}_{ell m n d} ranges between 0.75 and 1 (being often near 0.90, or larger), a fact indicating that the bounds are close to the sharp constant. Finally, we discuss the asymptotic behavior of the upper and lower bounds for K_{ell, b ell, c ell, d} when 1 <= b <= c and ell -> + Infinity. As an example, from this analysis we obtain the ell -> + Infinity limiting behavior of the sharp constant K_{ell, 2 ell, 2 ell, d}; a second example concerns the ell -> + Infinity limit for K_{ell, 2 ell, 3 ell, d}. The present work generalizes our previous paper [16], entirely devoted to the constant K_{ell m n d} in the special case ell = m = n; many results given therein can be recovered here for this special case.

preprint2009arXiv

On the averaging principle for one-frequency systems. An application to satellite motions

This paper is related to our previous works [1][2] on the error estimate of the averaging technique, for systems with one fast angular variable. In the cited references, a general method (of mixed analytical and numerical type) has been introduced to obtain precise, fully quantitative estimates on the averaging error. Here, this procedure is applied to the motion of a satellite in a polar orbit around an oblate planet, retaining only the J_2 term in the multipole expansion of the gravitational potential. To exemplify the method, the averaging errors are estimated for the data corresponding to two Earth satellites; for a very large number of orbits, computation of our estimators is much less expensive than the direct numerical solution of the equations of motion.

preprint2008arXiv

On the averaging principle for one-frequency systems. Seminorm estimates for the error

We extend some previous results of our work [1] on the error of the averaging method, in the one-frequency case. The new error estimates apply to any separating family of seminorms on the space of the actions; they generalize our previous estimates in terms of the Euclidean norm. For example, one can use the new approach to get separate error estimates for each action coordinate. An application to rigid body under damping is presented. In a companion paper [2], the same method will be applied to the motion of a satellite around an oblate planet.