Some Indentities for the Gradient of the Principal Invariants, Traces and Determinants via Grassmann Calculus

Let A be a second order tensor in a finite dimensional space. In this work we determine the gradient of the principal invariants of A and obtain some trace and determinant identities using only some standard rigorous statements concerning Grassmann calculus. We recover some of the results of Dui et al. (J. Elast. 75:193–196, 2004) and of Truesdell and Noll (The Non-linear Field Theories of Mechanics, Springer, Berlin, 2002) and solve an old problem proposed in SIAM Review concerning a determinant identity from a new perspective in a concise and simple manner.     

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 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.     

Well-Posedness for the Motion of Physical Vacuum of the Three-dimensional Compressible Euler Equations with or without Self-Gravitation

This paper concerns the well-posedness theory of the motion of a physical vacuum for the compressible Euler equations with or without self-gravitation. First, a general uniqueness theorem of classical solutions is proved for the three dimensional general motion. Second, for the spherically symmetric motions, without imposing the compatibility condition of the first derivative being zero at the center of symmetry, a new local-in-time existence theory is established in a functional space involving less derivatives than those constructed for three-dimensional motions in (Coutand et al., Commun Math Phys 296:559–587, 2010; Coutand and Shkoller, Arch Ration Mech Anal 206:515–616, 2012; Jang and Masmoudi, Well-posedness of compressible Euler equations in a physical vacuum, 2008) by constructing suitable weights and cutoff functions featuring the behavior of solutions near both the center of the symmetry and the moving vacuum boundary.     

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).     

Reciprocity in elastic media with rheology

The classic reciprocity theorem of Betti (Nuovo Cimento 2(VII–VIII):69–97, 1872) for perfectly elasto-static media has been extended to the dynamic fields by Graffi (Ann. Mat. Pura Appl. IV(18):173–200, 1939) for isotropic media with variable density. In this note the theorem is extended to a wide class of media with rheology.     

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.     

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).     

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.     

Microscale Fracture Toughness of Bismuth Doped Copper Bicrystals Using Double Edge Notched Microtensile Tests

The availability of focused ion beam (FIB) milling, nanoindentation, and microelectromechanical systems (MEMS) based test platforms has enabled small-scale mechanical testing to become an increasingly popular approach for measuring material properties. While great emphasis has been placed on measuring plastic properties at the micro- and nanoscale [1, 2], an area that has received significantly less consideration is the measurement of fracture toughness. A technique for performing small-scale, in situ fracture toughness tests using double edge notched tensile (DENT) specimens has been developed and used to measure a nearly 40 % reduction in toughness associated with the addition of Bi to the grain boundary of a Cu bicrystal. That Bi embrittles Cu grain boundaries is well known [3–10], however, as shown herein, the DENT technique offers certain advantages over existing boundary fracture tests, especially when used with ductile materials.     

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).     

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.     

A Quantitative Modulus of Continuity for the Two-Phase Stefan Problem

We derive the quantitative modulus of continuity $$\omega(r)=\left[ p+\ln \left( \frac{r_0}{r}\right)\right]^{-\alpha (n, p)},$$ which we conjecture to be optimal for solutions of the p-degenerate two-phase Stefan problem. Even in the classical case p = 2, this represents a twofold improvement with respect to the early 1980’s state-of-the-art results by Caffarelli– Evans (Arch Rational Mech Anal 81(3):199–220, 1983) and DiBenedetto (Ann Mat Pura Appl 103(4):131–176, 1982), in the sense that we discard one logarithm iteration and obtain an explicit value for the exponent α(n, p).     

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.     

Particle-based numerical modeling of AE statistics in disordered materials

We present some numerical simulations of AE due to damage propagation in disordered materials under compression and bending. To this purpose, the AE cumulative number, the time frequency analysis and the statistical properties of AE time series will be numerically simulated adopting the so-called “particle method strategy” (Cundall in Proceedings ISRM Symp., Nancy, France, vol. 2, pp. 129–136, 1971). The method provides the velocity of particles in a set simulating the behavior of a granular system and, therefore, is suitable to model the compressive wave propagation and acoustic emission (corresponding to cracking) in a solid body. The numerical simulations (Abe et al. in Pure Appl. Geophys. 161:2265–2277, 2004) correctly describe the compression test in terms of mean stress-strain response and crack pattern (Invernizzi et al. in Proceedings of the SEM annual conference, Society for Experimental Mechanics Inc., USA, SEM Annual Conference, Indianapolis, Indiana, USA, June 7–10, 2010). The size effects on the peak compressive strength and on the AE count are correctly reproduced. In addition, the amplitude distribution (b-value) and temporal evolution of AE events due to cracking, crucial for the evaluation of damage and remaining lifetime, were simulated and result in agreement with the experimental evidences.     

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).     

Quasi-linear versus potential-based formulations of force–flux relations and the GENERIC for irreversible processes: comparisons and examples

An essential part in modeling out-of-equilibrium dynamics is the formulation of irreversible dynamics. In the latter, the major task consists in specifying the relations between thermodynamic forces and fluxes. In the literature, mainly two distinct approaches are used for the specification of force–flux relations. On the one hand, quasi-linear relations are employed, which are based on the physics of transport processes and fluctuation–dissipation theorems (de Groot and Mazur in Non-equilibrium thermodynamics, North Holland, Amsterdam, 1962, Lifshitz and Pitaevskii in Physical kinetics. Volume 10, Landau and Lifshitz series on theoretical physics, Pergamon Press, Oxford, 1981). On the other hand, force–flux relations are also often represented in potential form with the help of a dissipation potential (?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997). We address the question of how these two approaches are related. The main result of this presentation states that the class of models formulated by quasi-linear relations is larger than what can be described in a potential-based formulation. While the relation between the two methods is shown in general terms, it is demonstrated also with the help of three examples. The finding that quasi-linear force–flux relations are more general than dissipation-based ones also has ramifications for the general equation for non-equilibrium reversible–irreversible coupling (GENERIC: e.g., Grmela and Öttinger in Phys Rev E 56:6620–6632, 6633–6655, 1997, Öttinger in Beyond equilibrium thermodynamics, Wiley Interscience Publishers, Hoboken, 2005). This framework has been formulated and used in two different forms, namely a quasi-linear (Öttinger and Grmela in Phys Rev E 56:6633–6655, 1997, Öttinger in Beyond equilibrium thermodynamics, Wiley Interscience Publishers, Hoboken, 2005) and a dissipation potential–based (Grmela in Adv Chem Eng 39:75–129, 2010, Grmela in J Non-Newton Fluid Mech 165:980–986, 2010, Mielke in Continuum Mech Therm 23:233–256, 2011) form, respectively, relating the irreversible evolution to the entropy gradient. It is found that also in the case of GENERIC, the quasi-linear representation encompasses a wider class of phenomena as compared to the dissipation-based formulation. Furthermore, it is found that a potential exists for the irreversible part of the GENERIC if and only if one does for the underlying force–flux relations.