Third order nonoscillatory central scheme for hyperbolic conservation laws

Summary. A third-order accurate Godunov-type scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: 1. A non-oscillatory piecewise-quadratic reconstruction of pointvalues from their given cell averages; and 2. A central differencing based on staggered evolution of the reconstructed cell averages. This results in a third-order central scheme, an extension along the lines of the second-order central scheme of Nessyahu and Tadmor \cite{NT}. The scalar scheme is non-oscillatory (and hence – convergent), in the sense that it does not increase the number of initial extrema (– as does the exact entropy solution operator). Extension to systems is carried out by componentwise application of the scalar framework. In particular, we have the advantage that, unlike upwind schemes, no (approximate) Riemann solvers, field-by-field characteristic decompositions, etc., are required. Numerical experiments confirm the high-resolution content of the proposed scheme. Thus, a considerable amount of simplicity and robustness is gained while retaining the expected third-order resolution. Received April 10, 1996 / Revised version received January 20, 1997     

The numerical radius and specttural matrices

In this paper we investigate spectral matrices, i.e., matrices with equal spectral and numerical radii. Various characterizations and properties of these matrices are given.     

Finite-temperature quasicontinuum: molecular dynamics without all the atoms

Using a combination of statistical mechanics and finite-element interpolation, we develop a coarse-grained (CG) alternative to molecular dynamics (MD) for crystalline solids at constant temperature. The new approach is significantly more efficient than MD and generalizes earlier work on the quasicontinuum method. The method is validated by recovering equilibrium properties of single crystal Ni as a function of temperature. CG dynamical simulations of nano-indentation reveal a strong dependence on temperature of the critical stress to nucleate dislocations under the indenter.     

The numerical radius and specttural matrices

In this paper we investigate spectral matrices, i.e., matrices with equal spectral and numerical radii. Various characterizations and properties of these matrices are given.     

Construction of Approximate Entropy Measure-Valued Solutions for Hyperbolic Systems of Conservation Laws

Entropy solutions have been widely accepted as the suitable solution framework for systems of conservation laws in several space dimensions. However, recent results in De Lellis and Székelyhidi Jr (Ann Math 170(3):1417–1436, 2009) and Chiodaroli et al. (2013) have demonstrated that entropy solutions may not be unique. In this paper, we present numerical evidence that state-of-the-art numerical schemes need not converge to an entropy solution of systems of conservation laws as the mesh is refined. Combining these two facts, we argue that entropy solutions may not be suitable as a solution framework for systems of conservation laws, particularly in several space dimensions. We advocate entropy measure-valued solutions, first proposed by DiPerna, as the appropriate solution paradigm for systems of conservation laws. To this end, we present a detailed numerical procedure which constructs stable approximations to entropy measure-valued solutions, and provide sufficient conditions that guarantee that these approximations converge to an entropy measure-valued solution as the mesh is refined, thus providing a viable numerical framework for systems of conservation laws in several space dimensions. A large number of numerical experiments that illustrate the proposed paradigm are presented and are utilized to examine several interesting properties of the computed entropy measure-valued solutions.     

Numerical radius of positive matrices

We study some properties of the numerical radius of matrices with non-negative entries, and explicit ways to compute it. We also characterize positive matrices with equal spectral and numerical radii, i.e., positive spectral matrices.     

Flocking With Short-Range Interactions

Journal of Statistical Physics - We study the large-time behavior of continuum alignment dynamics based on Cucker–Smale (CS)-type interactions which involve short-range kernels, that is,...     

On the role of the three-phase contact line in surface deformation

Viscoelastic braking theories developed by Shanahan and de Gennes and by others predict deformation of a solid surface at the solid-liquid-air contact line. This phenomenon has only been observed for soft smooth surfaces and results in a protrusion of the solid surface at the three-phase contact line, in agreement with the theoretical predictions. Despite the large (enough to break chemical bonds) forces associated with it, this deformation was not confirmed experimentally for hard surfaces, especially for hydrophobic ones. In this study we use superhydrophobic surfaces composed of an array of silicon nanostructures whose Young modulus is 4 orders of magnitude higher than that of surfaces in earlier recorded viscoelastic braking experiments. We distinguish between two cases: when a water drop forms an adhesive contact, albeit small, with the apparent contact angle θ < 180° and when the drop-surface adhesion is such that the conditions for placing a resting drop on the surface cannot be reached (i.e., θ = 180°). In the first case we show that there is a surface deformation at the three-phase contact line which is associated with a reduction in the hydrophobicity of the surface. For the second case, however, there cannot be a three-phase contact line associated with a drop in contact with the surface, and indeed, if we force-place a drop on the surface by holding it with a needle, no deformation is detected, nor is there a reduction in the hydrophobic properties of the surface. Yet, if we create a long horizontal three-phase contact line by partially immersing the superhydrophobic substrate in a water bath, we see a localized reduction in the hydrophobic properties of the surface in the region where the three-phase contact line used to be. The SEM scan of that region shows a narrow horizontal stripe where the nanorods are no longer there, and instead there is only a shallow structure that is lower than the nanorods height and resembles fused or removed nanorods. Away from that region, either on the part of the surface which was exposed to bulk water or the part which was exposed to air, no change in the hydrophobic properties of the surface is observed, and the SEM scan confirms that the nanorods seem intact in both regions.     

Shechtman’s three question marks: possibility, impossibility, and quasicrystals

The revolutionary discovery of actual quasicrystals, thanks to Dan Shechtman’s stamina, is a golden opportunity to analyze once again the role that pure (“theoretical”) possibilities and saving them plays in scientific progress. Some theoreticians, primarily Alan Mackay, contributed to saving pure possibilities of quasicrystalline structures and to opening materials science for them. My analysis rests upon an original modal metaphysics—panenmentalism—which I introduced and have been developing since 1999, quite independently of any familiarity with modern crystallography, and which deals with saving pure possibilities as indispensably contributing, inter alia, to our knowledge and sciences.