A failure of a class of K-BKZ equations based on principal stretches

The invariants in the K-BKZ constitutive equation for an incompressible viscoelastic fluid are usually taken to be the trace of the Finger strain tensor and its inverse. The basis for this choice of invariants is not derived from the K-BKZ theory, but rather is due to the perception that this is the most natural choice. Research into using other sets of invariants in the K-BKZ equation, such as the principal stretches or the eigenvalues of the Finger strain tensor (i.e., the squares of the principal stretches) is relatively new. We attempt here to derive a K-BKZ equation based on the squares of the principal stretches that models the behavior of a low-density polyethylene melt in simple shear and uniaxial elongational deformation. In doing so, two assumptions are made as to the form of the strain-dependent energy function: first, that there is a function f(q) such that the energy function can be written as the sum of f(q _i ),i = 1, 2, 3, where the q _i 'sare the squares of the principal stretches, and second that f is a power law. We find that the K-BKZ equation resulting from these two assumptions is inadequate to describe both the shear and elongational behavior of our material and we conclude that the second of the above assumptions is not valid. Further investigation, including predictions of the second normal stress difference and some finite element calculations reveals that the first assumption is also invalid for our material.     

Remarks on the links between low‐order DG methods and some finite‐difference schemes for the Stokes problem

In this paper we demonstrate that some well‐known finite‐difference schemes can be interpreted within the framework of the local discontinuous Galerkin (LDG) methods using the low‐order piecewise solenoidal discrete spaces introduced in (SIAM J. Numer. Anal. 1990; 27 (6): 1466–1485). In particular, it appears that it is possible to derive the well‐known MAC scheme using a first‐order Nédélec approximation on rectangular cells. It has been recently interpreted within the framework of the Raviart–Thomas approximation by Kanschat (Int. J. Numer. Meth. Fluids 2007; published online). The two approximations are algebraically equivalent to the MAC scheme, however, they have to be applied on grids that are staggered on a distance h/2 in each direction. This paper also demonstrates that both discretizations allow for the construction of a divergence‐free basis, which yields a linear system with a ‘biharmonic’ conditioning. Both this paper and Kanschat (Int. J. Numer. Meth. Fluids 2007; published online) demonstrate that the LDG framework can be used to generalize some popular finite‐difference schemes to grids that are not parallel to the coordinate axes or that are unstructured. Copyright © 2008 John Wiley & Sons, Ltd.     

On the ��principle of virtual power�� for arbitrary parts of a body

Recently, a “principle of virtual power” has been adopted to model the behavior of materials that involve multiple length scales. In these works, the “principle” is stated for arbitrary parts of a body and this arbitrariness is used, but not to its fullest extent, to draw conclusions concerning the structure of the theory that results. Here, a theorem and an example application are given to illustrate the restrictive nature of the requirement that it hold for arbitrary parts of a body and to draw attention to the full consequences that result from this requirement. Several key results that have been reported in the recent literature are incomplete, and this incompleteness has lead to superficial conclusions.     

Quasipatterns versus superlattices resulting from the superposition of two hexagonal patterns

We present a short review of the experimental observations and mechanisms related to the generation of quasipatterns and superlattices by the Faraday instability with two-frequency forcing. We show how two-frequency forcing makes possible triad interactions that generate hexagonal patterns, twelvefold quasipatterns or superlattices that consist of two hexagonal patterns rotated by an angle α relative to each other. We then consider which patterns could be observed when α does not belong to the set of prescribed values that give rise to periodic superlattices. Using the Swift–Hohenberg equation as a model, we find that quasipattern solutions exist for nearly all values of α. However, these quasipatterns have not been observed in experiments with the Faraday instability for $\alpha \ne \mathrm{\pi }/6$. We discuss possible reasons and mention a simpler framework that could give some hint about this problem.     

Instability of small-amplitude convective flows in a rotating layer with free boundaries

The stability of steady convective flows in a horizontal layer with free boundaries, heated from below and rotating about a vertical axis, is studied in the Boussinesq approximation (Rayleigh-Bénard convection). The flows considered are convective rolls or square cells that are sums of two perpendicular rolls with equal wave numbers k. It is assumed that the Rayleigh number is almost critical in order for convective flows with a wave number k: R = R _c (k) + ε² to arise, the amplitude of the supercritical states being of the order of ε. It is shown that the flows are always unstable relative to perturbations that are the sum of one long-and two short-wave modes corresponding to linear rolls turned through small angles in opposite directions.     

A study of long separation bubble on thick airfoils and its consequent effects

A parametric study has been performed to analyse the flow around the thick-symmetric NACA 0021 airfoil in order to better understand the characteristics and effects of long separation bubbles (LoSBs) that exist on such airfoils at low Reynolds numbers and turbulence intensities. In the article, the prediction capabilities of two recently-developed transition models, the correlation-based γ–Re_θ model and the laminar-kinetic-energy-based κ–κ_L–ω model are assessed. Two-dimensional steady-state simulations indicated that the κ–κ_L–ω model predicted the separation and reattachment process accurately when compared with published experimental work. The model was then used to study the attributes and the effects of LoSBs as a function of the angle of attack, freestream turbulence intensity and Reynolds number. It was observed that LoSBs considerably degrade the aerodynamic performance of airfoils and lead to abrupt stall behaviour. It is, furthermore, illustrated that the presence of the LoSB leads to an induced camber effect on the airfoil that increases as the airfoil angle of attack increases due to the upstream migration of the bubble. An increase in the Reynolds number or turbulence levels leads to a reduction in the bubble extent, considerably improving the airfoil performance and leading to a progressive trailing-edge stall.     

Some New Implications of Certain Constitutive Inequalities in Plane Isotropic Elasticity

In plane isotropic elasticity a strengthened form of the Ordered–Forces inequality is shown to imply that the restriction of the strain-energy function to the class of deformation gradients which share the same average of the principal stretches is bounded from below by the strain energy corresponding to the conformal deformations in this class. For boundary conditions of place, this property (together with a certain version of the Pressure–Compression inequality) is then used (i) to show that the plane radial conformal deformations are stable with respect to all radial variations of class C ¹ and (ii) to obtain explicit lower bounds for the total energy associated with arbitrary plane radial deformations. For the same type of boundary conditions and together with a different version of the Pressure–Compression inequality, an analogous property in plane isotropic elasticity (established in [3] under the assumption that the material satisfies a strengthened form of the Baker–Ericksen inequality and according to which the restriction of the strain-energy function to the class of deformation gradients which share the same determinant is bounded from below by the strain energy corresponding to the conformal deformations in that class) is used (i) to show that the plane radial conformal deformations are stable with respect to all variations of class C ¹ and (ii) to obtain explicit lower bounds for the total energy associated with any plane deformation.     

In General There Are No Unsheared Tetrads

Suppose the principal stretches are all different at a point P in a deformed body. In this case, it has been shown [1] that generally there is an infinity of non coplanar infinitesimal material line elements at P which remain unsheared following the deformation – that is, the angle between the arms of each pair of material line elements forming the triad remains unchanged. Here it is shown that in this case when all three principal stretches at P are different, there is no set of four infinitesimal material line elements, no three of which are coplanar, and such that the angle between each pair of the six pairs of material line elements is unchanged following the deformation. It is only when all three principal stretches at P are equal to each other, that there are unsheared tetrads at P, and in that case all tetrads are unsheared. This revised version was published online in June 2006 with corrections to the Cover Date.     

Radon Measure-Valued Solutions for a Class of Quasilinear Parabolic Equations

Initial value problems for quasilinear parabolic equations having Radon measures as initial data have been widely investigated, looking for solutions which for positive times take values in some function space. In contrast, it is the purpose of this paper to define and investigate solutions that for positive times take values in the space of the Radon measures of the initial data. We call such solutions measure-valued, in contrast to function-valued solutionspreviously considered in the literature. We first show that there is a natural notion of measure-valued solution of problem (P) below, in spite of its nonlinear character. A major consequence of our definition is that, if the space dimension is greater than one, the concentrated part of the solution with respect to the Newtonian capacity is constant in time. Subsequently, we prove that there exists exactly one solution of the problem, such that the diffuse part with respect to the Newtonian capacity of the singular part of the solution (with respect to the Lebesgue measure) is concentrated for almost every positive time on the set where “the regular part (with respect to the Lebesgue measure) is large”. Moreover, using a family of entropy inequalities we demonstrate that the singular part of the solution is nonincreasing in time. Finally, the regularity problem is addressed, as we give conditions (depending on the space dimension, the initial data and the rate of convergence at infinity of the nonlinearity ψ) to ensure that the measure-valued solution of problem (P) is, in fact, function-valued.     

On Extension and Torsion of Strain-Stiffening Rubber-Like Elastic Circular Cylinders

This paper is concerned with investigation of the effects of strain-stiffening on the response of solid circular cylinders in the combined deformation of torsion superimposed on axial extension. The cylinders are composed of incompressible isotropic nonlinearly elastic materials. Our primary focus is on materials that undergo severe strain-stiffening in the stress-stretch response. In particular, we consider two particular phenomenological constitutive models for such materials that reflect limiting chain extensibility at the molecular level. The axial stretch γ and twist that can be sustained in cylinders composed of such materials are shown to be constrained in a coupled fashion. It is shown that, in the absence of an additional axial force, a transition value γ=γ _t of the axial stretch exists such that for γ<γ _t , the stretched cylinder tends to elongate on twisting whereas for γ>γ _t , the stretched cylinder tends to shorten on twisting. These results are in sharp contrast with those for classical models such as the Mooney-Rivlin (and neo-Hookean) models that predict that the stretched circular cylinder always tends to further elongate on twisting. We also obtain results for materials modeled by the well-known exponential strain-energy widely used in biomechanics applications. This model reflects a strain-stiffening that is less abrupt than that for the limiting chain extensibility models. Surprisingly, it turns out that the results in this case are somewhat more complicated. For a fixed stiffening parameter, provided that the stretch is sufficiently small, the stretched bar always tends to elongate on twisting in the absence of an additional axial force. However, for sufficiently large stretch, the cylinder tends to shorten on undergoing sufficiently small twist but then tends to elongate on further twisting. These results are of interest in view of the widespread use of exponential models in the context of the mechanics of soft biological tissues. The special case of pure torsion is also briefly considered. In this case, the resultant axial force required to maintain pure torsion is compressive for all the models discussed here. In the absence of such a force, the bar would elongate on twisting reflecting the celebrated Poynting effect.      

Gas and Particle Dynamics of a Contoured Shock Tube for Pre-clinical Microparticle Drug Delivery

We investigate the gas-particle dynamics of a device designed for biological pre-clinical experiments. The device uses transonic/supersonic gas flow to accelerate microparticles such that they penetrate the outer skin layers. By using a shock tube coupled to a correctly expanded nozzle, a quasi-one-dimensional, quasi-steady flow (QSF) is produced to uniformly accelerate the microparticles. The system utilises a microparticle “cassette” (a diaphragm sealed container) that incorporates a jet mixing mechanism to stir the particles prior to diaphragm rupture. Pressure measurements reveal that a QSF exit period – suitable for uniformly accelerating microparticles – exists between 155 and 220 mus after diaphragm rupture. Immediately preceding the QSF period, a starting process secondary shock was shown to form with its (x,t) trajectory comparing well to theoretical estimates. To characterise the microparticle, flow particle image velocimetry experiments were conducted at the nozzle exit, using particle payloads with varying diameter (2.7–48 μm), density (600–16,800 kg/m³) and mass (0.25–10 mg). The resultant microparticle velocities were temporally uniform. The experiments also show that the starting process does not significantly influence the microparticle nozzle exit velocities. The velocity distribution across the nozzle exit was also uniform for the majority of microparticle types tested. For payload masses typically used in pre-clinical drug and vaccine applications (≤ 1 mg), it was demonstrated that payload scaling does not affect the microparticle exit velocities. These characteristics show that the microparticle exit conditions are well controlled and are in agreement with ideal theory. These features combined with an attention to the practical requirements of a pre-clinical system make the device suitable for investigating microparticle penetration into the skin for drug delivery.     

Potential flow of viscous fluids: Historical notes

In this essay I will attempt to identify the main events in the history of thought about irrotational flow of viscous fluids. I am of the opinion that when considering irrotational solutions of the Navier–Stokes equations it is never necessary and typically not useful to put the viscosity to zero. This observation runs counter to the idea frequently expressed that potential flow is a topic which is useful only for inviscid fluids; many people think that the notion of a viscous potential flow is an oxymoron. Incorrect statements like “… irrotational flow implies inviscid flow but not the other way around” can be found in popular textbooks.     

A numerical analysis of flat internal cracks under mixed mode loading

A numerical/analytical approach is proposed to determine the stress intensity factors K_I, K_II, and K_III of a 3D internal crack. The main point of this approach is the meshing technique that can model very sharp crack fronts. The meshing technique is based on an elliptical coordinate transformation that starts from a circular crack. It allows the obtainment of a curved crack front with elements normal to the crack front. Remarkable accuracy can be obtained for elliptical crack fronts with axes ratio smaller that 0.01. Accuracy demonstration is provided for cylindrical element with an inclined internal crack subjected to uni-axial tension. This case corresponds to crack propagation for all three modes of loading, the solution of which can checked with references’ results.     

The Remarkable Nature of Radially Symmetric Deformation of Spherically Uniform Linear Anisotropic Elastic Solids

Antman and Negron-Marrero [1] have shown the remarkable nature of a sphere of nonlinear elastic material subjected to a uniform pressure at the surface of the sphere. When the applied pressure exceeds a critical value the stress at the center r=0 of the sphere is infinite. Instead of nonlinear elastic material, we consider in this paper a spherically uniform linear anisotropic elastic material. It means that the stress-strain law referred to a spherical coordinate system is the same for any material point. We show that the same remarkable nature appears here. What distinguishes the present case from that considered in [1] is that the existence of the infinite stress at r=0 is independent of the magnitude of the applied traction σ0 at the surface of the sphere. It depends only on one nondimensional material parameter κ. For a certain range of κ a cavitation (if σ0>0) or a blackhole (if σ0<0) occurs at the center of the sphere. What is more remarkable is that, even though the deformation is radially symmetric, the material at any point need not be transversely isotropic with the radial direction being the axis of symmetry as assumed in [1]. We show that the material can be triclinic, i.e., it need not possess a plane of material symmetry. Triclinic materials that have as few as two independent elastic constants are presented. Also presented are conditions for the materials that are capable of a radially symmetric deformation to possess one or more symmetry planes. This revised version was published online in August 2006 with corrections to the Cover Date.     

Interfacial instability in turbulent flow over a liquid film in a channel

We revisit the stability of a deformable interface that separates a fully-developed turbulent gas flow from a thin layer of laminar liquid. Although this problem has received considerable attention previously, a model that requires no fitting parameters and that uses a base-state profile that has been validated against experiments is, as yet, unavailable. Furthermore, the significance of wave-induced perturbations in turbulent stresses remains unclear. To address these outstanding issues, we investigate this problem and introduce a turbulent base-state velocity that requires specification of a flow rate or a pressure drop only; no adjustable parameters are necessary. This base state is validated extensively against available experimental data as well as the results of direct numerical simulations. In addition, the effect of perturbations in the turbulent stress distributions is investigated, and demonstrated to be small for cases wherein the liquid layer is thin. The detailed modelling of the liquid layer also elicits two unstable modes, ‘interfacial’ and ‘internal’, with the former being the more dominant of the two. We show that it is possible for interfacial roughness to reduce the growth rate of the interfacial mode in relation to that of the internal one, promoting the latter, to the status of most dangerous mode. Additionally, we introduce an approximate measure to distinguish between ‘slow’ and ‘fast’ waves, the latter being the case for ‘critical-layer’-induced instabilities; we demonstrate that for the parameter ranges studied, the large majority of the waves are ‘slow’. Finally, comparisons of our linear stability predictions are made with experimental data in terms of critical parameters for onset of wave-formation, wave speeds and wavelengths; these yield agreement within the bounds of experimental error.     

One, no one, and one hundred thousand crack propagation laws: A generalized Barenblatt and Botvina dimensional analysis approach to fatigue crack growth

Barenblatt and Botvina with elegant dimensional analysis arguments have elucidated that Paris’ power-law is a weak form of scaling, so that the Paris’ parameters C and m should not be taken as material constants. On the contrary, they are expected to depend on all the dimensionless parameters of the problem, and are really “constants” only within some specific ranges of all these. In the present paper, the dimensional analysis approach by Barenblatt and Botvina is generalized to explore the functional dependencies of m and C on more dimensionless parameters than the original Barenblatt and Botvina, and experimental results are interpreted for a wider range of materials including both metals and concrete. In particular, we find that the size-scale dependencies of m and C and the resulting correlation between C and m are quite different for metals and for quasi-brittle materials, as it is already suggested from the fact the fatigue crack propagation processes lead to m=2-5 in metals and m=10-50 in quasi-brittle materials. Therefore, according to the concepts of complete and incomplete self-similarities, the experimentally observed breakdowns of the classical Paris’ law are discussed and interpreted within a unified theoretical framework. Finally, we show that most attempts to address the deviations from the Paris’ law or the empirical correlations between the constants can be explained with this approach. We also suggest that “incomplete similarity” corresponds to the difficulties encountered so far by the “damage tolerant” approach which, after nearly 50 years since the introduction of Paris’ law, is still not a reliable calculation of damage, as Paris himself admits in a recent review.     

Steady and transient sliding under rate-and-state friction

The physics of dry friction is often modelled by assuming that static and kinetic frictional forces can be represented by a pair of coefficients usually referred to as μ_s and μ_k, respectively. In this paper we re-examine this discontinuous dichotomy and relate it quantitatively to the more general, and smooth, framework of rate-and-state friction. This is important because it enables us to link the ideas behind the widely used static and dynamic coefficients to the more complex concepts that lie behind the rate-and-state framework. Further, we introduce a generic framework for rate-and-state friction that unifies different approaches found in the literature.We consider specific dynamical models for the motion of a rigid block sliding on an inclined surface. In the Coulomb model with constant dynamic friction coefficient, sliding at constant velocity is not possible. In the rate-and-state formalism steady sliding states exist, and analysing their existence and stability enables us to show that the static friction coefficient μ_s should be interpreted as the local maximum at very small slip rates of the steady state rate-and-state friction law.Next, we revisit the often-cited experiments of Rabinowicz (J. Appl. Phys., 22:1373–1379, 1951). Rabinowicz further developed the idea of static and kinetic friction by proposing that the friction coefficient maintains its higher and static value μ_s over a persistence length before dropping to the value μ_k. We show that there is a natural identification of the persistence length with the distance that the block slips as measured along the stable manifold of the saddle point equilibrium in the phase space of the rate-and-state dynamics. This enables us explicitly to define μ_s in terms of the rate-and-state variables and hence link Rabinowicz's ideas to rate-and-state friction laws.This stable manifold naturally separates two basins of attraction in the phase space: initial conditions in the first one lead to the block eventually stopping, while in the second basin of attraction the sliding motion continues indefinitely. We show that a second definition of μ_s is possible, compatible with the first one, as the weighted average of the rate-and-state friction coefficient over the time the block is in motion.     

Computer simulation of turbulence in internal combustion engines

The paper presents two- and three-dimensional computations of the in-cylinder turbulent flow in a diesel engine. The mathematical formulation is presented first, with emphasis on the modifications made to the standard k-ε model of turbulence, to account for rapid compression/expansion, and on the k-w model also used in the computations. Then, the results of two-and three-dimensional transient calculations are presented and compared with experimental data. It is realized that two-dimensional computations may be of little value to real engines, which would probably require three-dimensional analyses. However, two-dimensional studies are still useful in allowing the testing of new ideas easily and economically. It is concluded that the standard k-ε model may lead to poor predictions when used for internal combustion (IC) engine simulations, and that the modified model leads to more reasonable length-scale distributions, and it improves significantly the overall agreement of velocity predictions with experiment. The effect of the k-ε modification is apparent in both the two- and three-dimensional simulations. It is also demonstrated that the k-w model provides better turbulence predictions than the unmodified k-ε model, for the cases considered, and that a similar modification of the k-w model, to account for rapid compression/expansion, might improve its predictions even further.     

A theory of grain boundaries that accounts automatically for grain misorientation and grain-boundary orientation

This work is an attempt to answer the question:

Is there a physically natural method of characterizing the possible interactions between the slip systems of two grains that meet at a grain boundary—a method that could form the basis for the formulation of grain-boundary conditions?

Here we give a positive answer to this question based on the notion of a Burgers vector as described by a tensor field G on the grain boundary [Gurtin, M.E., Needleman, A., 2005. Boundary conditions in small-deformation single-crystal plasticity that account for the Burgers vector. J. Mech. Phys. Solids 53, 1-31]. We show that the magnitude of G can be expressed in terms of two types of moduli: inter-grain moduli that characterize slip-system interactions between the two grains; intra-grain moduli that for each grain characterize interactions between any two slip systems of that grain.We base the theory on microscopic force balances derived using the principle of virtual power, a version of the second law in the form of a free-energy imbalance, and thermodynamically compatible constitutive relations dependent on G and its rate. The resulting microscopic force balances represent flow rules for the grain boundary; and what is most important, these flow rules account automatically—via the intra- and inter-grain moduli—for the relative misorientation of the grains and the orientation of the grain boundary relative to those grains.