The Erpenbeck High Frequency Instability Theorem for Zeldovitch–von Neumann–D?ring Detonations

The rigorous study of spectral stability for strong detonations was begun by Erpenbeck (Phys. Fluids 5:604–614 1962). Working with the Zeldovitch–von Neumann–D?ring (ZND) model (more precisely, Erpenbeck worked with an extension of ZND to general chemistry and thermodynamics), which assumes a finite reaction rate but ignores effects such as viscosity corresponding to second order derivatives, he used a normal mode analysis to define a stability function V(t,e){V(\tau,\epsilon)} whose zeros in ${\mathfrak{R}\tau > 0}${\mathfrak{R}\tau > 0} correspond to multidimensional perturbations of a steady detonation profile that grow exponentially in time. Later in a remarkable paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966; Stability of detonations for disturbances of small transverse wavelength, 1965) he provided strong evidence, by a combination of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of V exist for perturbations of sufficiently large transverse wavenumber e{\epsilon} , even when the von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly 20 years later) by Majda. In spite of a great deal of later numerical work devoted to computing the zeros of V(t,e){V(\tau,\epsilon)} , the paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966) remains one of the few works we know of [another is Erpenbeck (Phys. Fluids 7:684–696, 1964), which considers perturbations for which the ratio of longitudinal over transverse components approaches ∞] that presents a detailed and convincing theoretical argument for detecting them. The analysis in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) points the way toward, but does not constitute, a mathematical proof that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) and provide proofs, together with certain simplifications and extensions, of the main conclusions about stability and instability of detonations contained in that paper. The main mathematical problem, and our principal focus here, is to determine the precise asymptotic behavior as e?￥{\epsilon\to\infty} of solutions to a linear system of ODEs in x, depending on e{\epsilon} and a complex frequency τ as parameters, with turning points x _* on the half-line [0,∞).     

Self-consistent modeling of entangled network strands and linear dangling structures in a single-strand mean-field slip-link model

Linear viscoelastic (LVE) measurements as well as non-linear elongation measurements have been performed on stoichiometrically imbalanced polymeric networks to gain insight into the structural influence on the rheological response (Jensen et al., Rheol Acta 49(1):1–13, 2010). In particular, we seek knowledge about the effect of dangling ends and soluble structures. To interpret our recent experimental results, we exploit a molecular model that can predict LVE data and non-linear stress–strain data. The slip-link model has proven to be a robust tool for both LVE and non-linear stress–strain predictions for linear chains (Khaliullin and Schieber, Phys Rev Lett 100(18):188302–188304, 2008, Macromolecules 42(19):7504–7517, 2009; Schieber, J Chem Phys 118(11):5162–5166, 2003), and it is thus used to analyze the experimental results. Initially, we consider a stoichiometrically balanced network, i.e., all strands in the ensemble are attached to the network in both ends. Next we add dangling strands to the network representing the stoichiometric imbalance, or imperfections during curing. By considering monodisperse network strands without dangling ends, we find that the relative low-frequency plateau, G₀/G_N⁰G_0/G_N^0, decreases linearly with the average number of entanglements. The decrease from G_N⁰G_N^0 to G ₀ is a result of monomer fluctuations between entanglements, which is similar to “longitudinal modes” in tube theory. It is found that the slope of G′ is dependent on the fraction of network strands and the structural distribution of the network. The power-law behavior of G ^″ is not yet captured quantitatively by the model, but our results suggest that it is a result of polydisperse dangling and soluble structures.     

Reynolds number dependence of the inertial range scaling of energy dissipation rate and enstrophy in a cylinder wake

The full energy dissipation rate and enstrophy are measured simultaneously using a probe consisting of four X-wires in the intermediate region of a cylinder wake for Taylor microscale Reynolds number in the range of 120–320. Longitudinal and transverse velocity increments are also obtained temporally using Taylor’s hypothesis. The inertial range scaling exponents indicate that the full enstrophy field has a stronger intermittency than does the full dissipation field for all the Reynolds numbers considered. The approximations of the energy dissipation rate and enstrophy based on isotropy are more intermittent than their corresponding true values. While the scaling exponents of the full energy dissipation rate remain approximately constant for different Reynolds numbers, those of the enstrophy decrease slightly and consistently with the increase of Reynolds number. It is conjectured that the scaling of the energy dissipation rate and the enstrophy may be the same when Reynolds number is extremely high, a trend that is consistent with that suggested by Nelkin (Phys Fluids 11:2202–2204, 1999; Am J Phys 68:310–318, 2000).     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 2: single-crystal plasticity

Variational formulations are constructed for rate-independent problems in small-deformation single-crystal strain-gradient plasticity. The framework, based on that of Gurtin (J Mech Phys Solids 50: 5–32, 2002), makes use of the flow rule expressed in terms of the dissipation function. Provision is made for energetic and dissipative microstresses. Both recoverable and non-recoverable defect energies are incorporated into the variational framework. The recoverable energies include those that depend smoothly on the slip gradients, the Burgers tensor, or on the dislocation densities (Gurtin et al. J Mech Phys Solids 55:1853–1878, 2007), as well as an energy proposed by Ohno and Okumura (J Mech Phys Solids 55:1879–1898, 2007), which leads to excellent agreement with experimental results, and which is positively homogeneous and therefore not differentiable at zero slip gradient. Furthermore, the variational formulation accommodates a non-recoverable energy due to Ohno et al. (Int J Mod Phys B 22:5937–5942, 2008), which is also positively homogeneous, and a function of the accumulated dislocation density. Conditions for the existence and uniqueness of solutions are established for the various examples of defect energy, with or without the presence of hardening or slip resistance.     

Long-Time Convergence of an Adaptive Biasing Force Method: The Bi-Channel Case

We present convergence results for an adaptive algorithm to compute free energies, namely the adaptive biasing force (ABF) method (Darve and Pohorille in J Chem Phys 115(20):9169–9183, 2001; Hénin and Chipot in J Chem Phys 121:2904, 2004). The free energy is the effective potential associated to a so-called reaction coordinate ξ(q), where q = (q ₁, … , q _3N ) is the position vector of an N-particle system. Computing free energy differences remains an important challenge in molecular dynamics due to the presence of metastable regions in the potential energy surface. The ABF method uses an on-the-fly estimate of the free energy to bias dynamics and overcome metastability. Using entropy arguments and logarithmic Sobolev inequalities, previous results have shown that the rate of convergence of the ABF method is limited by the metastable features of the canonical measures conditioned to being at fixed values of ξ (Lelièvre et al. in Nonlinearity 21(6):1155–1181, 2008). In this paper, we present an improvement on the existing results in the presence of such metastabilities, which is a generic case encountered in practice. More precisely, we study the so-called bi-channel case, where two channels along the reaction coordinate direction exist between an initial and final state, the channels being separated from each other by a region of very low probability. With hypotheses made on ‘channel-dependent’ conditional measures, we show on a bi-channel model, which we introduce, that the convergence of the ABF method is, in fact, not limited by metastabilities in directions orthogonal to ξ under two crucial assumptions: (i) exchange between the two channels is possible for some values of ξ and (ii) the free energy is a good bias in each channel. This theoretical result supports recent numerical experiments (Minoukadeh et al. in J Chem Theory Comput 6:1008–1017, 2010), where the efficiency of the ABF approach is demonstrated for such a multiple-channel situation.     

A constrained theory for single crystal shape memory wires with application to restrained recovery

The theory of thin wires developed in Dret and Meunier (Comptes Rendus de l’Académie des Sciences. Série I. Mathématique 337:143–147, 2003) is adapted to phase-transforming materials with large elastic moduli in the sense discussed in James and Rizzoni (J Elast 59:399–436, 2000). The result is a one-dimensional constitutive model for shape memory wires, characterized by a small number of material constants. The model is used to analyze self-accommodated and detwinned microstructures and to study superelasticity. It also turns out that the model successfully reproduces the behavior of shape memory wires in experiments of restrained recovery (Tsoi et al. in Mater Sci Eng A 368:299–310, 2004; Tsoi in 50:3535–3544, 2002; S̆ittner et al. in Mater Sci Eng A 286:298–311, 2000; vokoun in Smart Mater Struct 12:680–685, 2003; Zheng and Cui in Intermetallics 12:1305–1309, 2004; Zheng et al. in J Mater Sci Technol 20(4):390–394, 2004). In particular, the model is able to predict the shift to higher transformation temperatures on heating. The model also captures the effect of prestraining on the evolution of the recovery stress and of the martensite volume fraction.     

The role of dissipation and defect energy in variational formulations of problems in strain-gradient plasticity. Part 1: polycrystalline plasticity

A general set of flow laws and associated variational formulations are constructed for small-deformation rate-independent problems in strain-gradient plasticity. The framework is based on the thermodynamically consistent theory due to Gurtin and Anand (J Mech Phys Solids 53:1624–1649, 2005), and includes as variables a set of microstresses which have both energetic and dissipative components. The flow law is of associative type. It is expressed as a normality law with respect to a convex but otherwise arbitrary yield function, or equivalently in terms of the corresponding dissipation function. Two cases studied are, first, an extension of the classical Hill-Mises or J ₂ flow law and second, a form written as a linear sum of the magnitudes of the plastic strain and strain gradient. This latter form is motivated by work of Evans and Hutchinson (Acta Mater 57:1675–1688, 2009) and Nix and Gao (J Mech Phys Solids 46:411–425, 1998), who show that it leads to superior correspondence with experimental results, at least for particular classes of problems. The corresponding yield function is obtained by a duality argument. The variational problem is based on the flow rule expressed in terms of the dissipation function, and the problem is formulated as a variational inequality in the displacement, plastic strain, and hardening parameter. Dissipative components of the microstresses, which are indeterminate, are absent from the formulation. Existence and uniqueness of solutions are investigated for the generalized Hill-Mises and linear-sum dissipation functions, and for various combinations of defect energy. The conditions for well-posedness of the problem depend critically on the choice of dissipation function, and on the presence or otherwise of a defect energy in the plastic strain or plastic strain gradient, and of internal-variable hardening.     

Convergence Rates in <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript> for Elliptic Homogenization Problems

We study rates of convergence of solutions in L ² and H ^1/2 for a family of elliptic systems {L_e}{\{\mathcal{L}_\varepsilon\}} with rapidly oscillating coefficients in Lipschitz domains with Dirichlet or Neumann boundary conditions. As a consequence, we obtain convergence rates for Dirichlet, Neumann, and Steklov eigenvalues of {L_e}{\{\mathcal{L}_\varepsilon\}} . Most of our results, which rely on the recently established uniform estimates for the L ² Dirichlet and Neumann problems in Kenig and Shen (Math Ann 350:867–917, 2011; Commun Pure Appl Math 64:1–44, 2011) are new even for smooth domains.     

High-speed fiber spinning process with spinline flow-induced crystallization and neck-like deformation

The dynamics and stability of the high-speed fiber spinning process with spinline flow-induced crystallization and neck-like deformation have been studied using a simulation model equipped with governing equations of continuity, motion, energy, and crystallinity, along with the Phan-Thien–Tanner constitutive equation. Despite the fact that a simple one-phase model was incorporated into the governing equations to describe the spinline crystallinity, as opposed to the best-known two-phase model [Doufas et al. J Non-Newton Fluid Mech, 92:27–66, 2000a]; [Kohler et al. J Macromol Sci Phys, 44:185–202, 2005] that treats amorphous and crystalline phases separately in computing the spinline stress, the simulation has successfully portrayed the typical nonlinear characteristic of the high-speed spinning process called neck-like spinline deformation. It has been found that the criterion for the neck-like deformation to occur on the spinline is for the extensional viscosity to decrease on the spinline, so that the spinning is stabilized by the formation of the spinline neck-like deformation. The accompanying linear stability analysis explains this stabilizing effect of the spinline neck-like deformation, corroborating a recent experimental finding [Takarada et al. Int Polym Process, 19:380–387, 2004].This paper was presented at the 2nd Annual European Rheology Conference 2005 on April 21–23, 2005, in Grenoble, France.     

Group classification of one-dimensional nonisentropic equations of fluids with internal inertia

A systematic application of the group analysis method for modeling fluids with internal inertia is presented. The equations studied include models such as the nonlinear one-velocity model of a bubbly fluid (with incompressible liquid phase) at small volume concentration of gas bubbles (Iordanski Zhurnal Prikladnoj Mekhaniki i Tekhnitheskoj Fiziki 3, 102–111, 1960; Kogarko Dokl. AS USSR 137, 1331–1333, 1961; Wijngaarden J. Fluid Mech. 33, 465–474, 1968), and the dispersive shallow water model (Green and Naghdi J. Fluid Mech. 78, 237–246, 1976; Salmon 1988). These models are obtained for special types of the potential function W(r,[(r)\dot],S){W(\rho,\dot \rho,S)} (Gavrilyuk and Teshukov Continuum Mech. Thermodyn. 13, 365–382, 2001). The main feature of the present paper is the study of the potential functions with W _S  ≠ 0. The group classification separates these models into 73 different classes.     

Full 3D correlation tensor computed from double field stereoscopic PIV in a high Reynolds number turbulent boundary layer

The turbulence structure near a wall is a very active subject of research and a key to the understanding and modeling of this flow. Many researchers have worked on this subject since the fifties Hama et al. (J Appl Phys 28:388–394, 1957). One way to study this organization consists of computing the spatial two-point correlations. Stanislas et al. (C R Acad Sci Paris 327(2b):55–61, 1999) and Kahler (Exp Fluids 36:114–130, 2004) showed that double spatial correlations can be computed from stereoscopic particle image velocimetry (SPIV) fields and can lead to a better understanding of the turbulent flow organization. The limitation is that the correlation is only computed in the PIV plane. The idea of the present paper is to propose a new method based on a specific stereoscopic PIV experiment that allows the computation of the full 3D spatial correlation tensor. The results obtained are validated by comparison with 2D computation from SPIV. They are in very good agreement with the results of Ganapthisubramani et al. (J Fluid Mech 524:57–80, 2005a).     

From Rate-Dependent to Rate-Independent Brittle Crack Propagation

On the base of many experimental results, e.g., Ravi-Chandar and Knauss (Int. J. Fract. 26:65–80, 1984), Sharon et al. (Phys. Rev. Lett. 76(12):2117–2120, 1996), Hauch and Marder (Int. J. Fract. 90:133–151, 1998), the object of our analysis is a rate-dependent model for the propagation of a crack in brittle materials. Restricting ourselves to the quasi-static framework, our goal is a mathematical study of the evolution equation in the geometries of the ‘Single Edge Notch Tension’ and of the ‘Compact Tension’. Besides existence and uniqueness, emphasis is placed on the regularity of the evolution making reference also to the ‘velocity gap’. The transition to the rate-independent model of Griffith is obtained by time rescaling, proving convergence of the rescaled evolutions and of their energies. Further, the discontinuities of the rate-independent evolution are characterized in terms of unstable points of the free energy. Results are illustrated by a couple of numerical examples in the above mentioned geometries.     

Convergence of Ginzburg–Landau Functionals in Three-Dimensional Superconductivity

In this paper we consider the asymptotic behavior of the Ginzburg–Landau model for superconductivity in three dimensions, in various energy regimes. Through an analysis via Γ-convergence, we rigorously derive a reduced model for the vortex density and deduce a curvature equation for the vortex lines. In the companion paper (Baldo et al. Commun. Math. Phys. 2012, to appear) we describe further applications to superconductivity and superfluidity, such as general expressions for the first critical magnetic field H _c1, and the critical angular velocity of rotating Bose–Einstein condensates.     

Use of particle image velocimetry to study heterogeneous drag reduction

Large polymer filaments can form when drag reducing polymers are injected through wall slots. The presence of these structures enhances the performance of the drag reducing function by mechanisms which are not understood. This paper shows how particle image velocimetry (PIV) techniques can be used to study changes in the configuration of the injected polymer and in the structure of the velocity field with increasing drag reduction. The filaments are found to behave as solid bodies which break up in high shear regions close to a boundary. The breakup process provides an explanation of why the filaments are not observed close to a wall and offers the possibility of providing a heterogeneous distribution of small aggregates of polymers which could be more effective than uniformly distributed molecules as suggested by Hoyer and Gyr (J Non-Newton Fluid Mech 65:221–240, 1996; J Fluids Eng 120:818–823, 1998), Dunlop and Cox (Phys Fluids 20:203–213, 1977) and Vlachogiannis et al. (Phys Fluid 15:3786–3794, 2004). PIV measurements show dramatic qualitative changes in the velocity patterns at maximum drag reduction.     

On Entropy Flux of Transversely Isotropic Elastic Bodies

Recently (Liu in J. Elast. 90:259–270, 2008) thermodynamic theory of elastic (and viscoelastic) material bodies has been analyzed based on the general entropy inequality. It is proved that for isotropic elastic materials, the results are identical to the classical results based on the Clausius-Duhem inequality (Coleman and Noll in Arch. Ration. Mech. Anal. 13:167–178, 1963), for which one of the basic assumptions is that the entropy flux is defined as the heat flux divided by the absolute temperature. For anisotropic elastic materials in general, this classical entropy flux relation has not been proved in the new thermodynamic theory. In this note, as a supplement of the theory presented in (Liu in J. Elast. 90:259–270, 2008), it will be proved that the classical entropy flux relation need not be valid in general, by considering a transversely isotropic elastic material body.      

Fractional order theory in thermoelastic solid with three-phase lag heat transfer

In this work, the field equations of the linear theory of thermoelasticity have been constructed in the context of a new consideration of Fourier law of heat conduction with time-fractional order and three-phase lag. A uniqueness and reciprocity theorems are proved. One-dimensional application for a half-space of elastic material in the presence of heat sources has been solved using Laplace transform and state space techniques Ezzat (Canad J Phys Rev 86:1241–1250, 2008). According to the numerical results and its graphs, conclusion about the new theory has been established.     

Reliable treatment of a new analytical method for solving MHD boundary-layer equations

The purpose of this study is to implement a new analytical method which is a combination of the homotopy analysis method (HAM) and the Padé approximant for solving magnetohydrodynamic boundary-layer flow. The solution is compared with the numerical solution. Comparisons between the HAM–Padé and the numerical solution reveal that the new technique is a promising tool for solving MHD boundary-layer equations. The effects of the various parameters on the velocity and temperature profiles are presented graphically form. Favorable comparisons with previously published works (Crane, J. Appl. Math. Phys. 21:645–647, 1970, and Vajravelu and Hadjinicolaou, Int. J. Eng. Sci. 35:1237–1244, 1997) are obtained. It is predicted that HAM–Padé can have wide application in engineering problems (especially for boundary-layer and natural convection problems).     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

Rice Local Phase Angle Study for a Delamination Problem Between a Shape Memory Alloy and an Elastic Material

The present study is an extension of a recent paper of Freed et al. (J Mech Phys Solids 56:3003–3020, 2008). The final aim is to describe the transformation toughening behavior of a static crack along an interface between a shape memory alloy (SMA) and a linear elastic isotropic material. With an SMA as an equivalent Huber–Von Mises stress model (hypothesis of symmetric behavior between tension and compression), Freed et al. determine the initiation (ending) phase transformation yield surfaces in terms of the local phase angle introduced by Rice et al. (Metal ceramic interfaces, Pergamon Press, New York, pp 269–294, 1990). In this paper we give the general framework to determine this angle for a model integrating the asymmetry between tension and compression (experimentally measured: Vacher and Lexcellent in Proc ICM 6:231–236, 1991; Orgéas and Favier in Acta Mater 46(15):5579–5591, 2000), the Huber–Von Mises model being only a particular case. We demonstrate the local phase angle existence in an appropriate framing domain and give a sufficient hypothesis for its uniqueness and an algorithm to obtain it. Estimates are obtained in terms of physical quantities such as the Young modulus ratio, the bimaterial Poisson modulus values and also the choice of the yield loading functions. Finally, we illustrate this theoretical study by an application linking the asymmetry intensity on the width and the shape on predicted phase transformation surfaces and by a comparison with the symmetric case.