Entropy analysis for third grade fluid flow with Vogel model viscosity in annular pipe

The steady-state flow of a third grade fluid between concentric circular cylinders is considered and entropy generation due to fluid friction and heat transfer in the annular pipe is examined. Depending upon the fluid viscosity, entropy generation in the flow system varies. The third grade fluid is employed to account for the non-Newtonian effect while Vogel model is accommodated for temperature-dependent viscosity. The analysis is based on perturbation technique. The closed form solutions for velocity, temperature and entropy fields are presented. Entropy generation due to fluid friction and heat transfer in the flow system is formulated. The influence of viscosity parameters A and B on the entropy generation number is investigated. It is found that entropy generation number reduces with increasing viscosity parameter A, which is more pronounced in the region close to the annular pipe inner wall and opposite is true for increasing viscosity parameter B.     

Finite deformation analysis of crack-tip opening in elastic-plastic materials and implications for fracture

Analyses of the stress and strain fields around smoothly-blunting crack tips in both non-hardening and hardening elastic-plastic materials, under contained plane-strain yielding and subject to mode I opening loads, have been carried out by use of a finite element method suitably formulated to admit large geometry changes. The results include the crack-tip shape and near-tip deformation field, and the crack-tip opening displacement has been related to a parameter of the applied load, the J-integral. The hydrostatic stresses near the crack tip are limited due to the lack of constraint on the blunted tip, limiting achievable stress levels except in a very small region around the crack tip in power-law hardening materials. The J-integral is found to be path-independent except very close to the crack tip in the region affected by the blunted tip. Models for fracture are discussed in the light of these results including one based on the growth of voids. The rate of void-growth near the tip in hardening materials seems to be little different from the rate in non-hardening ones when measured in terms of crack-tip opening displacement, which leads to a prediction of higher toughness in hardening materials. It is suggested that improvement of this model would follow from better understanding of void-void and void-crack coalescence and void nucleation, and some criteria and models for these effects are discussed. The implications of the finite element results for fracture criteria based on critical stress or strain, or both, is discussed with respect to transition of fracture mode and the angle of initial crack-growth. Localization of flow is discussed as a possible fracture model and as a model for void-crack coalescence.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

Étude de la rupture fragile d'un stratifié quasi-isotrope à plis tissés : mise en évidence d'une longueur interne

This Note is an experimental and numerical study of the fracture of a quasi-isotropic laminate. It is shown that a local criteria fail to predict accurately such a fracture in the presence of large stress gradients. Good agreement is obtained with a nonlocal criteria involving a characteristic length which is determined experimentally. To cite this article: N. Lahellec et al., C. R. Mecanique 333 (2005).     

Étude numérique de la topologie du champ résiduel de déplacements d'un milieu ductile fissuré : influence sur la formation des caustiquesNumerical study of the topology of displacements residual field for a cracked ductile medium: Influence on the shape of the caustic curve

The purpose of our work is to illustrate that, by using a simple optical method (phenomenon of caustics) it is possible to locate, near the crack tip, the boundaries of the plastic zone for a cracked body. In this Note, we show that, through numerical simulations, in the case of a ductile SEN specimen, loaded in mode I and following a loading history, the appearance and development of the plastic zone involves irreversible topological changes in the vicinity of the crack tip, made obvious by caustics shape variations. To cite this article: O. Pop et al., C. R. Mecanique 332 (2004).     

MHD flow and heat transfer of a UCM fluid over a stretching surface with variable thermophysical properties

In this paper we investigate the effects of temperature-dependent viscosity, thermal conductivity and internal heat generation/absorption on the MHD flow and heat transfer of a non-Newtonian UCM fluid over a stretching sheet. The governing partial differential equations are first transformed into coupled non-linear ordinary differential equation using a similarity transformation. The resulting intricate coupled non-linear boundary value problem is solved numerically by a second order finite difference scheme known as Keller-Box method for various values of the pertinent parameters. Numerical computations are performed for two different cases namely, zero and non-zero values of the fluid viscosity parameter. That is, 1/?? _r ??0 and 1/?? _r ??0 to get the effects of the magnetic field and the Maxwell parameter on the velocity and temperature fields, for several physical situations. Comparisons with previously published works are presented as special cases. Numerical results for the skin-friction co-efficient and the Nusselt number with changes in the Maxwell parameter and the fluid viscosity parameter are tabulated for different values of the pertinent parameters. The results obtained for the flow characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid phenomena, especially the UCM fluid phenomena. Maxwell fluid reduces the wall-shear stress.     

Oscillatory fracture paths in thin elastic sheets

We report a novel mode of quasi-static oscillatory crack propagation when a cutting tip of moderately large width is driven through a thin brittle polymer film. Experiments show that the amplitude and wavelength of the oscillatory crack paths scale linearly with the width of the cutting tip over a wide range of length scales but are independent of the width of the sheet and of the cutting speed. We propose a mechanism for this instability, based on the coupling between crack propagation and out-of-plane deformations of the film. To cite this article: B. Roman et al., C. R. Mecanique 331 (2003).     

Propagation of a plane-strain hydraulic fracture with a fluid lag: Early-time solution

This paper studies the propagation of a plane-strain fluid-driven fracture with a fluid lag in an elastic solid. The fracture is driven by a constant rate of injection of an incompressible viscous fluid at the fracture inlet. The leak-off of the fracturing fluid into the host solid is considered negligible. The viscous fluid flow is lagging behind an advancing fracture tip, and the resulting tip cavity is assumed to be filled at some specified low pressure with either fluid vapor (impermeable host solid) or pore-fluids infiltrating from the permeable host solid. The scaling analysis allows to reduce problem parametric space to two lumped dimensionless parameters with the meaning of the solid toughness and of the tip underpressure (difference between the specified pressure in the tip cavity and the far field confining stress). A constant lumped toughness parameter uniquely defines solution trajectory in the parametric space, while time-varying lumped tip underpressure parameter describes evolution along the trajectory. Further analysis identifies the early and large time asymptotic states of the fracture evolution as corresponding to the small and large tip underpressure solutions, respectively. The former solution is obtained numerically herein and is characterized by a maximum fluid lag (as a fraction of the crack length), while the latter corresponds to the zero-lag solution of Spence and Sharp [Spence, D.A., Sharp, P.W., 1985. Self-similar solution for elastohydrodynamic cavity flow. Proc. Roy. Soc. London, Ser. A (400), 289–313]. The self-similarity at small/large tip underpressure implies that the solution for crack length, crack opening and net fluid pressure in the fluid-filled part of the crack is a given power-law of time, while the fluid lag is a constant fraction of the increasing fracture length. Evolution of a fluid-driven fracture between the two limit states corresponds to gradual expansion of the fluid-filled region and disappearance of the fluid lag. For small solid toughness and small tip underpressure, the fracture is practically devoid of fluid, which is localized into a narrow region near the fracture inlet. Corresponding asymptotic solution on the fracture lengthscale corresponds to that of a crack loaded by a pair of point forces which magnitude is determined from the coupled hydromechanical solution in the fluid-filled region near the crack inlet. For large solid toughness, the fluid lag is vanishingly small at any underpressure and the solution is adequately approximated by the zero-lag self-similar large toughness solution at any stage of fracture evolution. The small underpressure asymptotic solutions obtained in this work are sought to provide initial condition for the propagation of fractures which are initially under plane-strain conditions.     

Equilibrium of a pressurized plastic fluid in a wellbore crack

This paper investigates equilibrium of a pressurized plastic fluid invading a tensile wellbore crack in a linear elastic, permeable rock. The crack is initially filled by pore fluid at ambient pressure, that is immiscibly displaced by the plastic fluid invading from the wellbore. The plastic fluid comes to rest to form a “plug” within the elastically deformed crack when the limit equilibrium between the shear stresses generated at the fracture walls and the pressure drop between the wellbore wall and the crack tip is reached. The model assumes that the leak-off of the plastic fluid into the rock is negligible, while the displaced pore fluid in the crack tip region is freely exchanged with the surrounding permeable rock to maintain the ambient pressure level. When the crack length ? is small or large compared to the wellbore radius R, the problem reduces to that of a pressurized edge or Griffith’s crack, respectively, subjected to a uniform far-field confining stress. In these two end-member cases, the normalized solution for the net pressure distribution, the plug length, and the stress intensity factor at the crack tip is obtained as a function of two numbers – the normalized net fluid pressure at the crack inlet and at the crack tip (partial plugs only) – that embody the solution’s dependence on the wellbore and the far field loading, the fluid yield strength, and the rock modulus. In the general case of an intermediate crack length (? ～ R), the normalized solution is a function of two additional parameters, the length-to-radius ratio and a normalized measure of the far field stress anisotropy, respectively, which accurate approximation is devised from an end-member solution using a rescaling argument. The equilibrium plug solutions are used to evaluate the breakdown pressure, the critical wellbore pressure at which the crack propagation condition is first met, and to analyze the stability of the ensuing crack propagation.     

Entropy generation for pipe flow of a third grade fluid with Vogel model viscosity

The flow of fluid-solid mixtures in a pipe can be treated as non-Newtonian fluids of third grade. Depending upon the fluid viscosity, entropy generation in the flow system varies. In the present study, flow of third grade fluid in a pipe is considered. The Vogel model is introduced to account for the temperature-dependent viscosity. Entropy generation due to fluid friction and heat transfer in the flow system is formulated. The influence of viscosity parameters A and B on the entropy generation number is investigated. It is found that increasing viscosity parameter A reduces the entropy generation number and opposite is true for increasing viscosity parameter B.     

Theoretical development and experimental validation of a thermally dissipative cohesive zone model for dynamic fracture of amorphous polymers

A thermally dissipative cohesive zone model is developed for predicting the temperature increase at the tip of a crack propagating dynamically in a nominally brittle material exhibiting a cohesive-type failure such as crazing. The model assumes that fracture energy supplied to the crack tip region that is in excess of that needed for the creation of new free surfaces during crack advance is converted to heat within the cohesive zone. Bulk dissipation mechanisms, such as plasticity, are not accounted for. Several cohesive traction laws are examined, and the model is then used to make predictions of crack tip heating at various crack propagation speeds in the nominally brittle amorphous polymer PMMA, observed to fail by a crazing-type mechanism. The heating predictions are compared to experimental data where the temperature field surrounding a high speed crack in PMMA was measured. Measurements are made in real time using a multi-point high speed HgCdTe infrared radiation detector array. At the same time as temperature, simultaneous measurement of fracture energy is made by a strain gauge technique, and crack tip speed is monitored through a resistance ladder method. Material strength can be estimated through uniaxial tension tests, thus minimizing the need for parameter fitting in the stress-opening traction law. Excellent agreement between experiments and theory is found for two of the cohesive traction law temperature predictions, but only for the case where a single craze is active during the dynamic fracture of PMMA, i.e. crack tip speed up to approximately 0.2c_R. For higher speed fracture where subsurface damage becomes prominent, the line dissipation model of a cohesive zone is inadequate, and a distributed damage model is needed.     

Waves in an elastic tube filled with a heterogeneous fluid of variable viscosity

By treating the artery as a prestressed thin elastic tube and the blood as an incompressible heterogeneous fluid with variable viscosity, we studied the propagation of weakly non-linear waves in such a composite medium through the use of reductive perturbation method. By assuming a variable density and a variable viscosity for blood in the radial direction we obtained the perturbed Korteweg-deVries equation as the evolution equation when the viscosity is of order of ε^3/2. We observed that the perturbed character is the combined result of the viscosity and the heterogeneity of the blood. A progressive wave type of solution is presented for the evolution equation and the result is discussed. The numerical results indicate that for a certain value of the density parameter sigma, the wave equation loses its dispersive character and the evolution equation degenerates. It is further shown that, for the perturbed KdV equation both the amplitude and the wave speed decay in the time parameter τ.     

Analysis of the thermodynamic conditions for brittle fractureAnalyse des conditions thermodynamique pour la rupture fragile.

The analysis of the thermodynamic conditions for brittle fracture is given for the two known models of an isolated defect. In the first model the stresses on the external surface of the solid remain the same as before and after formation of a defect. In the second model the displacements on the external surface of the solid at a defect formation remain the same as in the solid without a defect. It is shown that the first model in the isothermal case of deformation leads to Griffith condition and the second model leads to the other proposed energy condition of fracture which, unlike the Griffith condition, contains an increment of the entropy component of the internal energy which is not zero in a general case. To cite this article: I.M. Dunaev, V.I. Dunaev, C. R. Mecanique 332 (2004).     

Linear instability of the electrified free interface between two cylindrical shells of viscoelastic fluids through porous media

In this paper, we have discussed the linear stability analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interfacial insoluble surfactant and surface charge through porous media. The case of long waves interfacial stability has been studied. The dispersion relation is solved numerically and hence the effects of various parameters are illustrated graphically. Our results reveal that the influence of the physicochemical parameter β is to shrink the instability region of the surface and reduce the growth rate of the unstable normal modes. Such important effects of the surfactant on the shape of interfacial structures are more sensitive to the variation of the β corresponding to non-Newtonian fluids-model compared with the Newtonian fluids model. In the case of long wave limit, it is demonstrated that increasing β, has a dual role in-fluence (de-stabilizing effects) depending on the viscosity of the core fluid. It has a destabilizing effect at the large values of the core fluid viscosity coefficient, while this role is exchanged to a regularly stabilizing influence at small values of such coefficient.     

Approche analytique de la convection forcée des fluides de Bingham dans un tubeAnalytical approach for forced convection of Bingham fluids in a tube

The asymptotic behaviour of laminar forced convection for Bingham fluid in a circular tube subjected to an axially varying wall heat flux is studied analytically. The effect of viscous dissipation is taken into account while the axial heat conduction in the fluid is considered as negligible. The asymptotic temperature profile and the asymptotic Nusselt number are determined for various axial wall heat flux distributions which yield a thermally developed region. The results obtained show a diminution in asymptotic Nusselt number when the Brinkman number and the dimensionless radius of the plug flow region increase. Comparisons with the results in the literature for Newtonian fluids show the validity of the present analysis. To cite this article: R. Khatyr et al., C. R. Mecanique 330 (2002) 69–75.     

A new approach to the pseudo-orthogonal properties of eigenfunction expansion form of the crack-tip complex potential function in anisotropic and piezoelectric fracture mechanics

Over the past twenty years, the well-known weight function theory based on the Bueckner work conjugate integral has been widely used to calculate crack tip fracture dominant parameter such as the stress intensity factor, the energy release rate (or the J-integral) and the T-stress in various kinds of cracked materials (e.g. isotropic materials, anisotropic materials and piezoelectric materials). Meanwhile, the pseudo-orthogonal property of the eigenfunction expansion form of the crack tip stress complex potential function has been proved to play a very important role in the theory. In this paper, we provide a new approach to establish the pseudo-orthogonal properties for crack problems in anisotropic and/or piezoelectric materials. In the latter case associated with mechanical-electric coupling, the electrical boundary conditions under both impermeable and permeable crack models are considered. The approach developed is much simpler than the classical complex variable separation technique proposed by previous researchers and hence the cumbersome and lengthy manipulations are avoided. Moreover, it is shown that, unlike previous works, the orthogonal properties of the material characteristic matrices A and B induced by the Stroh theory are no longer necessary in establishing the pseudo-orthogonal properties of eigenfunction expansion form in cracked piezoelectric materials. The approach can be easily extended to treat many other different crack problems concerning the Bueckner integral involving several complex arguments.     

Weakly-imposed Dirichlet boundary conditions for non-Newtonian fluid flow

We propose a new formulation for weakly imposing Dirichlet boundary conditions in non-Newtonian fluid flow. It is based on the Gerstenberger–Wall formulation for Newtonian fluids [1], but extended to non-Newtonian fluids. It uses a stabilization term in the weak form that is independent from the actual fluid model used, except for an adjustable parameter κ, having the physical dimension of a viscosity. The new formulation is tested, combined with an extended finite element method, for the flow past a cylinder between two walls using both a generalized Newtonian and a viscoelastic fluid. It is shown that the convergence is optimal for the generalized Newtonian fluid by comparing with a converged boundary-fitted solution using traditional strong boundary conditions. Also the solution of the viscoelastic fluid compares very well with a traditional solution using a boundary-fitted mesh and strong Dirichlet boundary conditions. For both fluid models we also test various values of the κ parameter and it turns out that a value equal to the zero-shear-viscosity gives good results. But, it is also shown that a wide range of κ values can be chosen without sacrificing accuracy.     

A new criterion for mixed mode fracture initiation based on the crack tip plastic core region

In this work, we propose a new criterion for mixed mode I-II crack initiation angles based on the characteristics of the plastic core region surrounding the crack tip. The shape and size of the plastic core region are thoroughly analyzed under different loading conditions and a new formulation for the non-dimensional variable radius of the core region is presented for mixed mode (K_I−K_II) fracture. The proposed criterion states that the crack extends in the direction of the local or global minimum of the plastic core region boundary depending on the resultant stress state at the crack tip. The results show a well-defined correlation between the plastic core region characteristics and crack extension angles predicted by other criteria. The proposed criterion is formulated for various loading conditions and is compared with other available criteria against the limited available experimental data. It is shown that the proposed criterion provides a better agreement with the experimental data.     

Lattice Boltzmann simulations of viscoplastic fluid flows through complex flow channels

We present the results of lattice Boltzmann (LB) simulations for the planar-flow of viscoplastic fluids through complex flow channels. In this study, the Bingham and Casson model fluids are covered as viscoplastic fluid. The Papanastasiou (modified Bingham) model and the modified Casson model are employed in our LB simulations. The Bingham number is an essential physical parameter when considering viscoplastic fluid flows and the modified Bingham number is proposed for modified viscoplastic models. When the value of the modified Bingham number agrees with that of the “normal” Bingham number, viscoplastic fluid flows formulated by modified viscoplastic models strictly reproduce the flow behavior of the ideal viscoplastic fluids. LB simulations are extensively performed for viscoplastic fluid flows through complex flow channels with rectangular and circular obstacles. It is shown that the LB method (LBM) allows us to successfully compute the flow behavior of viscoplastic fluids in various complicated-flow channels with rectangular and circular obstacles. For even low Re and high Bn numbers corresponding to plastic-property dominant condition, it is clearly manifested that the viscosity for both the viscoplastic fluids is largely decreased around solid obstacles. Also, it is shown that the viscosity profile is quite different between both the viscoplastic fluids due to the inherent nature of the models. The viscosity of the Bingham fluid sharply drops down close to the plastic viscosity, whereas the viscosity of the Casson fluid does not rapidly fall. From this study, it is demonstrated that the LBM can be also an effective methodology for computing viscoplastic fluid flows through complex channels including circular obstacles.     

Parametric excitation of waves on a free boundary of a horizontal fluid layer

We consider the dynamical stability of horizontal fluid layer, performing harmonic oscillations in vertical direction. The continued fractions approach allowed us to avoid the conventional restriction to the case of small viscosity and almost-resonant frequencies. Our numerical results cover a wide range of the parameters (viscosity, amplitude and frequency of the oscillation, and depth of the layer). To cite this article: V.I. Yudovich et al., C. R. Mecanique 332 (2004).