The Two-Dimensional Euler Equations on Singular Domains

We establish the existence of global weak solutions of the two-dimensional incompressible Euler equations for a large class of non-smooth open sets. Loosely, these open sets are the complements (in a simply connected domain) of a finite number of obstacles with positive Sobolev capacity. Existence of weak solutions with L ^p vorticity is deduced from a property of domain continuity for the Euler equations that relates to the so-called γ-convergence of open sets. Our results complete those obtained for convex domains in Taylor (Progress in Nonlinear Differential Equations and their Applications, Vol. 42, 2000), or for domains with asymptotically small holes (Iftimie et al. in Commun Partial Differ Equ 28(1–2), 349–379, 2003; Lopes Filho in SIAM J Math Anal 39(2), 422–436, 2007).     

On the Energy Dissipation Rate of Solutions to the Compressible Isentropic Euler System

In this paper we extend and complement the results in Chiodaroli et al. (Global ill-posedness of the isentropic system of gas dynamics, 2014) on the well-posedness issue for weak solutions of the compressible isentropic Euler system in 2 space dimensions with pressure law p(ρ) = ρ ^γ , γ ≥ 1. First we show that every Riemann problem whose one-dimensional self-similar solution consists of two shocks admits also infinitely many two-dimensional admissible bounded weak solutions (not containing vacuum) generated by the method of De Lellis and Székelyhidi (Ann Math 170:1417–1436, 2009), (Arch Ration Mech Anal 195:225–260, 2010). Moreover we prove that for some of these Riemann problems and for 1 ≤ γ < 3 such solutions have a greater energy dissipation rate than the self-similar solution emanating from the same Riemann data. We therefore show that the maximal dissipation criterion proposed by Dafermos in (J Diff Equ 14:202–212, 1973) does not favour the classical self-similar solutions.     

Weak Solutions of the Navier–Stokes Equations for Compressible Flows in a Half-Space with No-Slip Boundary Conditions

We consider the Navier–Stokes equations for the motion of compressible, viscous flows in a half-space ${\mathbb{R}^n_+,}$ n =  2,  3, with the no-slip boundary conditions. We prove the existence of a global weak solution when the initial data are close to a static equilibrium. The density of the weak solution is uniformly bounded and does not contain a vacuum, the velocity is Hölder continuous in (x, t) and the material acceleration is weakly differentiable. The weak solutions of this type were introduced by D. Hoff in Arch Ration Mech Anal 114(1):15–46, (1991), Commun Pure and Appl Math 55(11):1365–1407, (2002) for the initial-boundary value problem in ${\Omega = \mathbb{R}^n}$ and for the problem in ${\Omega = \mathbb{R}^n_+}$ with the Navier boundary conditions.     

Comparison of dynamical cores for NWP models: comparison of COSMO and Dune

We present a range of numerical tests comparing the dynamical cores of the operationally used numerical weather prediction (NWP) model COSMO and the university code Dune, focusing on their efficiency and accuracy for solving benchmark test cases for NWP. The dynamical core of COSMO is based on a finite difference method whereas the Dune core is based on a Discontinuous Galerkin method. Both dynamical cores are briefly introduced stating possible advantages and pitfalls of the different approaches. Their efficiency and effectiveness is investigated, based on three numerical test cases, which require solving the compressible viscous and non-viscous Euler equations. The test cases include the density current (Straka et al. in Int J Numer Methods Fluids 17:1–22, 1993), the inertia gravity (Skamarock and Klemp in Mon Weather Rev 122:2623–2630, 1994), and the linear hydrostatic mountain waves of (Bonaventura in J Comput Phys 158:186–213, 2000).     

Isentropic Gas Flow for the Compressible Euler Equation in a Nozzle

We study the motion of isentropic gas in a nozzle. Nozzles are used to increase the thrust of engines or to accelerate a flow from subsonic to supersonic. Nozzles are essential parts for jet engines, rocket engines and supersonicwind tunnels. In the present paper, we consider unsteady flow, which is governed by the compressible Euler equation, and prove the existence of global solutions for the Cauchy problem. For this problem, the existence theorem has already been obtained for initial data away from the sonic state, (Liu in Commun Math Phys 68:141–172, 1979). Here, we are interested in the transonic flow, which is essential for engineering and physics. Although the transonic flow has recently been studied (Tsuge in J Math Kyoto Univ 46:457–524, 2006; Lu in Nonlinear Anal Real World Appl 12:2802–2810, 2011), these papers assume monotonicity of the cross section area. Here, we consider the transonic flow in a nozzle with a general cross section area. When we prove global existence, the most difficult point is obtaining a bounded estimate for approximate solutions. To overcome this, we employ a new invariant region that depends on the space variable. Moreover, we introduce a modified Godunov scheme. The corresponding approximate solutions consist of piecewise steady-state solutions of an auxiliary equation, which yield a desired bounded estimate. In order to prove their convergence, we use the compensated compactness framework.     

Local Controllability to Trajectories of the Magnetohydrodynamic Equations

In this paper we prove the local controllability to trajectories of the three dimensional magnetohydrodynamic equations by means of two internal controls, one in the velocity equations and the other in the magnetic field equations and both localized in an arbitrary small subset with not empty interior of the domain. This paper improves the previous results (Barbu et al. in Comm Pure Appl Math 56:732–783, 2003; Barbu et al. in Adv Differ Equ 10:481–504, 2005; Havârneanu et al. in Adv Differ Equ 11:893–929, 2006; Havârneanu, in SIAM J Control Optim 46:1802–1830, 2007) where the second control is not localized and it allows to deduce the local controllability to trajectories with boundary controls. The proof relies on the Carleman inequality for the Stokes system of Imanuvilov et al. (Carleman estimates for second order nonhomogeneous parabolic equations, preprint) to deal with the velocity equations and on a new Carleman inequality for a Dynamo-type equation to deal with the magnetic field equations.     

On the Formulation of Continuum Thermodynamic Models for Solids as General Equations for Non-equilibrium Reversible-Irreversible Coupling

The purpose of this work is the comparison of some aspects of the formulation of material models in the context of continuum thermodynamics (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997) with their formulation in the form of a General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC: e.g., Grmela and Öttinger in Phys. Rev. E 56: 6620–6632, 1997; Öttinger and Grmela in Phys. Rev. E 56: 6633–6655, 1997; Öttinger in Beyond equilibrium thermodynamics, Wiley, New York, 2005; Grmela in J. Non-Newton. Fluid Mech. 165: 980–998, 2010). A GENERIC represents a generalization of the Ginzburg-Landau model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to formulate non-equilibrium thermodynamic models for complex fluids, it has recently been applied to anisotropic inelastic solids in a Eulerian setting (Hütter and Tervoort in J. Non-Newton. Fluid Mech. 152: 45–52, 2008; 53–65, 2008; Adv. Appl. Mech. 42: 254–317, 2009) as well as to damage mechanics (Hütter and Tervoort in Acta Mech. 201: 297–312, 2008). In the current work, attention is focused for simplicity on the case of thermoelastic solids with heat conduction and viscosity in a Lagrangian setting (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997, Chaps. 9–12). In the process, the relation of the two formulations to each other is investigated in detail. A particular point in this regard is the concept of dissipation and its model representation in both contexts.     

Singular Limiting Induced from Continuum Solutions and the Problem of Dynamic Cavitation

In the works of Pericak-Spector and Spector (Arch Rational Mech Anal. 101:293–317, 1988, Proc. Royal Soc. Edinburgh Sect A 127:837–857, 1997) a class of self-similar solutions are constructed for the equations of radial isotropic elastodynamics that describe cavitating solutions. Cavitating solutions decrease the total mechanical energy and provide a striking example of non-uniqueness of entropy weak solutions (for polyconvex energies) due to point-singularities at the cavity. To resolve this paradox, we introduce the concept of singular limiting induced from continuum solution (or slic-solution), according to which a discontinuous motion is a slic-solution if its averages form a family of smooth approximate solutions to the problem. It turns out that there is an energetic cost for creating the cavity, which is captured by the notion of slic-solution but neglected by the usual entropic weak solutions. Once this cost is accounted for, the total mechanical energy of the cavitating solution is in fact larger than that of the homogeneously deformed state. We also apply the notion of slic-solutions to a one-dimensional example describing the onset of fracture, and to gas dynamics in Langrangean coordinates with Riemann data inducing vacuum in the wave fan.     

h-Principles for the Incompressible Euler Equations

In Dissipative Euler Flows and Onsager’s Conjecture. arxiv.1205.3626, preprint, 2012, De Lellis and Székelyhidi construct Hölder continuous, dissipative (weak) solutions to the incompressible Euler equations in the torus ${{\mathbb T}^3}$ . The construction consists of adding fast oscillations to the trivial solution. We extend this result by establishing optimal h-principles in two and three space dimensions. Specifically, we identify all subsolutions (defined in a suitable sense) which can be approximated in the H ^?1-norm by exact solutions. Furthermore, we prove that the flows thus constructed on ${{\mathbb T}^3}$ are genuinely three-dimensional and are not trivially obtained from solutions on ${{\mathbb T}^2}$ .     

Concerning the Existence of Classical Solutions to the Stokes System. On the Minimal Assumptions Problem

In this notes we consider the stationary Stokes system in a bounded, connected, three-dimensional smooth domain, with homogeneous Dirichlet boundary condition. Proofs also apply to the n-dimensional case, and to other boundary conditions, like Navier-slip ones. We say here that a solution is classical if all derivatives appearing in the equations are continuous up to the boundary. It is well known, for long time, that solutions of the Stokes system are classical if the external forces belong to the H?lder space \({C^{0,\; \lambda}(\bar{\Omega})}\) . It is also well known that, in general, solutions are not classical in the presence of continuous external forces. Hence, a very challenging problem is to find Banach spaces, strictly containing the H?lder spaces \({C^{0,\; \lambda}(\bar{\Omega})}\) such that solutions to the Stokes problem corresponding to forces in the above space are classical. We prove this result for external forces in a suitable functional space, denoted \({{\rm C}_*(\bar{\Omega})}\) , introduced in references Beirão da Veiga (On the solutions in the large of the two-dimensional flow of a non-viscous incompressible fluid, 1982) and Beirão da Veiga (J Differ Equ 54(3):373–389, 1984) in connection with the Euler equations.     

Analytical Validation of a Continuum Model for Epitaxial Growth with Elasticity on Vicinal Surfaces

Within the context of heteroepitaxial growth of a film onto a substrate, terraces and steps self-organize according to misfit elasticity forces. Discrete models of this behavior were developed by Duport et al. (J Phys I 5:1317–1350, 1995) and Tersoff et al. (Phys Rev Lett 75:2730–2733, 1995). A continuum limit of these was in turn derived by Xiang (SIAM J Appl Math 63:241–258, 2002) (see also the work of Xiang and Weinan Phys Rev B 69:035409-1–035409-16, 2004; Xu and Xiang SIAM J Appl Math 69:1393–1414, 2009). In this paper we formulate a notion of weak solution to Xiang’s continuum model in terms of a variational inequality that is satisfied by strong solutions. Then we prove the existence of a weak solution.     

On the Inviscid Limit for the Compressible Navier–Stokes System in an Impermeable Bounded Domain

In this paper we investigate the issue of the inviscid limit for the compressible Navier–Stokes system in an impermeable fixed bounded domain. We consider two kinds of boundary conditions. The first one is the no-slip condition. In this case we extend the famous conditional result (Kato, T.: Remarks on zero viscosity limit for nonstationary Navier–Stokes flows with boundary. In: Seminar on nonlinear partial differential equations, vol. 2, pp. 85–98. Math. Sci. Res. Inst. Publ., Berkeley 1984) obtained by Kato in the homogeneous incompressible case. Kato proved that if the energy dissipation rate of the viscous flow in a boundary layer of width proportional to the viscosity vanishes then the solutions of the incompressible Navier–Stokes equations converge to some solutions of the incompressible Euler equations in the energy space. We provide here a natural extension of this result to the compressible case. The other case is the Navier condition which encodes that the fluid slips with some friction on the boundary. In this case we show that the convergence to the Euler equations holds true in the energy space, as least when the friction is not too large. In both cases we use in a crucial way some relative energy estimates proved recently by Feireisl et al. in J. Math. Fluid Mech. 14:717–730 (2012).     

Traveling Wave Solutions for a Predator?CPrey System With Sigmoidal Response Function

We study the existence of traveling wave solutions for a diffusive predator?Cprey system. The system considered in this paper is governed by a Sigmoidal response function which in some applications is more realistic than the Holling type I, II responses, and more general than a simplified form of the Holling type III response considered before. Our method is an improvement to the original method introduced in the work of Dunbar (J Math Biol 17:11?C32, 1983; SIAM J Appl Math 46:1057?C1078, 1986). A bounded Wazewski set is used in this work while unbounded Wazewski sets were used in Dunbar (1983, 1986). The existence of traveling wave solutions connecting two equilibria is established by using the original Wazewski??s theorem which is much simpler than the extended version in Dunbar??s work.     

Experimental Measurement of Local Burning Velocity Within a Rotating Flow

The work presented in this paper details the implementation of a new technique for the measurement of local burning velocity using asynchronous particle image velocimetry. This technique uses the local flow velocity ahead of the flame front to measure the movement of the flame by the surrounding fluid. This information is then used to quantify the local burning velocity by taking into account the translation of the flame via convection. In this paper the developed technique is used to study the interaction between a flame front and a single toroidal vortex for the case of premixed stoichiometric methane and air combustion. This data is then used to assess the impact of vortex structure on flame propagation rates. The burning velocity data demonstrates that there is a significant enhancement to the rate of flame propagation where the flame directly interacts with the rotating vortex. The increases found were significantly higher than expected but are supported by burning velocities (Filatyev et al, Combust Flame 141:1?C21, 2005; Kobayashi et al, Proc Combust Inst 29:1793?C1800, 2002; Shepherd et al. 1998) found in turbulent flames of the same mixture composition. Away from this interaction with the main vortex core, the flame exhibits propagation rates around the value recorded in literature for unperturbed laminar combustion (Tahtouh et al, Combust Flame 159:1735?C1743, 2009; Hassan et al, Combust Flame 115:539?C550, 1998); Halter et al, Proc Combust Inst 30:201?C208, 2005; Coppens et al, Exp Therm Fluid Sci 31:437?C444, 2007).     

Synchronization of piecewise continuous systems of fractional order

This paper proves analytically that synchronization of a class of piecewise continuous fractional-order systems can be achieved. Since there are no dedicated numerical methods to integrate differential equations with discontinuous right-hand sides for fractional-order models, Filippov’s regularization (Filippov, Differential Equations with Discontinuous Right-Hand Sides, 1988) is applied, and Cellina’s Theorem (Aubin and Cellina, Differential Inclusions Set-valued Maps and Viability Theory, 1984; Aubin and Frankowska, Set-valued Analysis, 1990) is used. It is proved that the corresponding initial value problem can be converted to a continuous problem of fractional-order systems, to which numerical methods can be applied. In this way, the synchronization problem is transformed into a standard problem for continuous fractional-order systems. Three examples are presented: the Sprott’s system, Chen’s system, and Shimizu–Morioka’s system.     

Existence Result for a Dynamical Equations of Generalized Marguerre-von Kármán Shallow Shells

In a recent work in the static case, Gratie (Appl. Anal. 81:1107–1126, 2002) has generalized the classical Marguerre-von Kármán equations studied by Ciarlet and Paumier in (Comput. Mech. 1:177–202, 1986), where only a portion of the lateral face of the shallow shell is subjected to boundary conditions of von Kármán type, while the remaining portion is subjected to boundary conditions of free edge. Then Ciarlet and Gratie (Math. Mech. Solids 11:83–100, 2006) have established an existence theorem for these equations. In Chacha et al. (Rev. ARIMA 13:63–76, 2010), we extended formally these studies to the dynamical case. More precisely, we considered a three-dimensional dynamical model for a nonlinearly elastic shallow shell with a specific class of boundary conditions of generalized Marguerre-von Kármán type. Using technics from formal asymptotic analysis, we showed that the scaled three-dimensional solution still leads to two-dimensional dynamical boundary value problem called the dynamical equations of generalized Marguerre-von Kármán shallow shells. In this paper, we establish the existence of solutions to these equations using a compactness method of Lions (Quelques Méthodes de Résolution des Problèmes aux Limites non Linéaires, Dunod, Paris, 1969).     

Estimation of the largest Lyapunov exponent-like (LLEL) stability measure parameter from the perturbation vector and its derivative dot product (part 2) experiment simulation

Controlling system dynamics with use of the Largest Lyapunov Exponent (LLE) is employed in many different areas of the scientific research. Thus, there is still need to elaborate fast and simple methods of LLE calculation. This article is the second part of the one presented in Dabrowski (Nonlinear Dyn 67:283–291, 2012). It develops method LLEDP of the LLE estimation and shows that from the time series of two identical systems, one can simply extract value of the stability parameter which value can be treated as largest LLE. Unlike the method presented in part, one developed method (LLEDPT) can be applied to the dynamical systems of any type, continuous, with discontinuities, with time delay and others. The theoretical improvement shows simplicity of the method and its obvious physical background. The proofs for the method effectiveness are based on results of the simulations of the experiments for Duffing and Van der Pole oscillators. These results were compared with ones obtained with use of the Stefanski method (Stefanski in Chaos Soliton Fract 11(15):2443–2451, 2000; Chaos Soliton Fract 15:233–244, 2003; Chaos Soliton Fract 23:1651–1659, 2005; J Theor Appl Mech 46(3):665–678, 2008) and LLEDP method. LLEDPT can be used also as the criterion of stability of the control system, where desired behavior of controlled system is explicitly known (Balcerzak et al. in Mech Mech Eng 17(4):325–339, 2013). The next step of development of the method can be considered in direction that allows estimation of LLE from the real time series, systems with discontinuities, with time delay and others.     

On Finite Time Singularity and Global Regularity of an Axisymmetric Model for the 3D Euler Equations

In (Comm Pure Appl Math 62(4):502–564, 2009), Hou and Lei proposed a 3D model for the axisymmetric incompressible Euler and Navier–Stokes equations with swirl. This model shares a number of properties of the 3D incompressible Euler and Navier–Stokes equations. In this paper, we prove that the 3D inviscid model with an appropriate Neumann-Robin or Dirichlet-Robin boundary condition will develop a finite time singularity in an axisymmetric domain. We also provide numerical confirmation for our finite time blowup results. We further demonstrate that the energy of the blowup solution is bounded up to the singularity time, and the blowup mechanism for the mixed Dirichlet-Robin boundary condition is essentially the same as that for the energy conserving homogeneous Dirichlet boundary condition. Finally, we prove that the 3D inviscid model has globally smooth solutions for a class of large smooth initial data with some appropriate boundary condition. Both the analysis and the results we obtain here improve the previous work in a rectangular domain by Hou et al. (Adv Math 230:607–641, 2012) in several respects.     

Flow temporal reconstruction from non-time-resolved data part I: mathematic fundamentals

At least two circumstances point to the need of postprocessing techniques to recover lost time information from non-time-resolved data: the increasing interest in identifying and tracking coherent structures in flows of industrial interest and the high data throughput of global measuring techniques, such as PIV, for the validation of computational fluid dynamics (CFD) codes. This paper offers the mathematic fundamentals of a space--time reconstruction technique from non-time-resolved, statistically independent data. An algorithm has been developed to identify and track traveling coherent structures in periodic flows. Phase-averaged flow fields are reconstructed with a correlation-based method, which uses information from the Proper Orthogonal Decomposition (POD). The theoretical background shows that the snapshot POD coefficients can be used to recover flow phase information. Once this information is recovered, the real snapshots are used to reconstruct the flow history and characteristics, avoiding neither the use of POD modes nor any associated artifact. The proposed time reconstruction algorithm is in agreement with the experimental evidence given by the practical implementation proposed in the second part of this work (Legrand et al. in Exp Fluids, 2011), using the coefficients corresponding to the first three POD modes. It also agrees with the results on similar issues by other authors (Ben Chiekh et al. in 9 Congrès Francophone de Vélocimétrie Laser, Bruxelles, Belgium, 2004; Van Oudheusden et al. in Exp Fluids 39-1:86?C98, 2005; Meyer et al. in 7th International Symposium on Particle Image Velocimetry, Rome, Italy, 2007a; in J Fluid Mech 583:199?C227, 2007b; Perrin et al. in Exp Fluids 43-2:341?C355, 2007). Computer time to perform the reconstruction is relatively short, of the order of minutes with current PC technology.