A particular case of initial data for the two-dimensional Euler equations is studied numerically. The results show that the Godunov method does not always converge to the physical solution, at least not on feasible grids. Moreover, they suggest that entropy solutions (in the weak entropy inequality sense) are not well posed.
(TODO: mention that @cite_28 fig 30 p. 296 is our example if the wedge is replaced by stagnation air; see p. 345 Fig 69. Quote: All these and other mathematically possible flow patterns with a singular center Z are at our disposal for interpreting experimental evidence. Which, if any, of these possibilities occurs under given circumstances is a question that cannot possibly be decided within the framework of a theory with such a high degree of indeterminacy. Here we have a typical instance of a theory incomplete and oversimplified in its basic assumptions; only by going more deeply into the physical basis of our theory, i.e. by accounting for heat conduction and viscosity, can we hope to clarify completely the phenomena at a three-shock singularity. It may well be that the boundary layer which develops along the constant discontinuity line modifies the flow pattern sufficiently to account for the observed deviation; [...quote Liepmann paper]''
@cite_28
2261902652 1982359839 2055457008 1993711495
In this article we consider large energy wave maps in dimension 2+1, as in the resolution of the threshold conjecture by Sterbenz and Tataru (Commun. Math. Phys. 298(1):139–230, 2010; Commun. Math. Phys. 298(1):231–264, 2010), but more specifically into the unit Euclidean sphere ( S ^ n-1 R ^ n ) with ( n ), and study further the dynamics of the sequence of wave maps that are obtained in Sterbenz and Tataru (Commun. Math. Phys. 298(1):231–264, 2010) at the final rescaling for a first, finite or infinite, time singularity. We prove that, on a suitably chosen sequence of time slices at this scaling, there is a decomposition of the map, up to an error with asymptotically vanishing energy, into a decoupled sum of rescaled solitons concentrating in the interior of the light cone and a term having asymptotically vanishing energy dispersion norm, concentrating on the null boundary and converging to a constant locally in the interior of the cone, in the energy space. Similar and stronger results have been recently obtained in the equivariant setting by several authors (Cote, Commun. Pure Appl. Math. 68(11):1946–2004, 2015; Cote, Commun. Pure Appl. Math. 69(4):609–612, 2016; Cote, Am. J. Math. 137(1):139–207, 2015; , Am. J. Math. 137(1):209–250, 2015; Krieger, Commun. Math. Phys. 250(3):507–580, 2004), where better control on the dispersive term concentrating on the null boundary of the cone is provided, and in some cases the asymptotic decomposition is shown to hold for all time. Here, however, we do not impose any symmetry condition on the map itself and our strategy follows the one from bubbling analysis of harmonic maps into spheres in the supercritical regime due to Lin and Riviere (Ann. Math. 149(2):785–829, 1999; Duke Math. J. 111:177–193, 2002), which we make work here in the hyperbolic context of Sterbenz and Tataru (Commun. Math. Phys. 298(1), 231–264, 2010). We propose and analyze a discontinuous Galerkin approximation for the Stokes problem. The finite element triangulation employed is not required to be conforming and we use discontinuous pressures and velocities. No additional unknown fields need to be introduced, but only suitable bilinear forms defined on the interfaces between the elements, involving the jumps of the velocity and the average of the pressure. We consider hp approximations using ℚk′–ℚk velocity-pressure pairs with k′ = k + 2, k + 1, k. Our methods show better stability properties than the corresponding conforming ones. We prove that our first two choices of velocity spaces ensure uniform divergence stability with respect to the mesh size h. Numerical results show that they are uniformly stable with respect to the local polynomial degree k, a property that has no analog in the conforming case. An explicit bound in k which is not sharp is also proven. Numerical results show that if equal order approximation is chosen for the velocity and pressure, no spurious pressure modes are present but the method is not uniformly stable either with respect to h or k. We derive a priori error estimates generalizing the abstract theory of mixed methods. Optimal error estimates in h are proven. As for discontinuous Galerkin methods for scalar diffusive problems, half of the power of k is lost for p and hp pproximations independently of the divergence stability. The viscously dominated, low Reynolds' number dynamics of multi-phase, compacting media can lead to nonlinear, dissipationless dispersive behavior when viewed appropriately. In these systems, nonlinear self-steepening competes with wave dispersion, giving rise to dispersive shock waves (DSWs). Example systems considered here include magma migration through the mantle as well as the buoyant ascent of a low density fluid through a viscously deformable conduit. These flows are modeled by a third-order, degenerate, dispersive, nonlinear wave equation for the porosity (magma volume fraction) or cross-sectional area, respectively. Whitham averaging theory for step initial conditions is used to compute analytical, closed form predictions for the DSW speeds and the leading edge amplitude in terms of the constitutive parameters and initial jump height. Novel physical behaviors are identified including backflow and DSW implosion for initial jumps sufficient to cause gradient catastrophe in the Whitham modulation equations. Theoretical predictions are shown to be in excellent agreement with long-time numerical simulations for the case of small to moderate amplitude DSWs. Verifiable criteria identifying the breakdown of this modulation theory in the large jump regime, applicable to a wide class of DSW problems, are presented. We study the phase separation phenomenon in the Ising model in dimensions d ≥ 3. To this end we work in a large box with plus boundary conditions and we condition the system to have an excess amount of negative spins so that the empirical magnetization is smaller than the spontaneous magnetization m*. We confirm the prediction of the phenomenological theory by proving that with high probability a single droplet of the minus phase emerges surrounded by the plus phase. Moreover, the rescaled droplet is asymptotically close to a definite deterministic shape, the Wulff crystal, which minimizes the surface free energy. In the course of the proof we establish a surface order large deviation principle for the magnetization. Our results are valid for temperatures T below a limit of slab-thresholds T c conjectured to agree with the critical point T c . Moreover, T should be such that there exist only two extremal translation invariant Gibbs states at that temperature, a property which can fail for at most countably many values and which is conjectured to be true for every T. The proofs are based on the Fortuin-Kasteleyn representation of the Ising model along with coarse-graining techniques. To handle the emerging macroscopic objects we employ tools from geometric measure theory which provide an adequate framework for the large deviation analysis. Finally, we propose a heuristic picture that for subcritical temperatures close enough to T c , the dominant minus spin cluster of the Wulff droplet permeates the entire box and has a strictly positive local density everywhere.