Optimal anisotropic three-phase conducting composites: Plane problem

The paper establishes tight lower bound for effective conductivity tensor K₁ of two-dimensional three-phase conducting anisotropic composites and defines optimal microstructures. It is assumed that three materials are mixed with fixed volume fractions and that the conductivity of one of the materials is infinite. The bound expands the Hashin–Shtrikman and translation bounds to multiphase structures, it is derived using a combination of translation method and additional inequalities on the fields in the materials; similar technique was used by Nesi, 1995, Cherkaev, 2009 for isotropic multiphase composites. This paper expands the bounds to the anisotropic composites with effective conductivity tensor K₁. The lower bound of conductivity (G-closure) is a piece-wise analytic function of eigenvalues of K₁, that depends only on conductivities of components and their volume fractions. Also, we find optimal microstructures that realize the bounds, developing the technique suggested earlier by Albin et al., 2007a, Cherkaev, 2009. The optimal microstructures are laminates of some rank for all regions. The found structures match the bounds in all but one region of parameters; we discuss the reason for the gap and numerically estimate it.     

Pointwise Bounds for the Solutions of the Smoluchowski Equation with Diffusion

We prove various decay bounds on solutions (f _n : n > 0) of the discrete and continuous Smoluchowski equations with diffusion. More precisely, we establish pointwise upper bounds on n ^? f _n in terms of a suitable average of the moments of the initial data for every positive ?. As a consequence, we can formulate sufficient conditions on the initial data to guarantee the finiteness of ${L^p(\mathbb{R}^d \times [0, T])}$ norms of the moments ${X_a(x, t) := \sum_{m\in\mathbb{N}}m^a f_m(x, t)}$ , ( ${\int_0^{\infty} m^a f_m(x, t)dm}$ in the case of continuous Smoluchowski’s equation) for every ${p \in [1, \infty]}$ . In previous papers [11] and [5] we proved similar results for all weak solutions to the Smoluchowski’s equation provided that the diffusion coefficient d(n) is non-increasing as a function of the mass. In this paper we apply a new method to treat general diffusion coefficients and our bounds are expressed in terms of an auxiliary function ${\phi(n)}$ that is closely related to the total increase of the diffusion coefficient in the interval (0, n].     

Small scale intermittency in a rough wall turbulent boundary layer

Relatively good quality isotropy is observed in the central region of a turbulent boundary layer developing over a mesh screen rough wall. Spectra of velocity and, more especially, vorticity fluctuations satisfy isotropy over a significant wavenumber range. Inertial range scaling exponents ζ _u2 (p) and ζ _u3 (p) of moments of order p(?8) of increments of the transverse velocity fluctuations u ₂ and u ₃ are significantly smaller than the exponents ζ _u1 (p) of increments of the longitudinal velocity fluctuation u ₁. Exponents inferred from the locally averaged values of squared transverse vorticity fluctuations are only slightly smaller than ζ _u1 (p). The difference between ζ _u1 (p) and ζ _u2 (p) [or ζ _u3 (p)] more likely reflects the departure from isotropy of inertial range scales. There is evidence to suggest that the difference decreases with an increase in the Reynolds number and/or a decrease in the magnitude of the mean shear.     

Nonlinear stability of Poiseuille flow of a Bingham fluid: theoretical results and comparison with phenomenological criteria

We present new results on the nonlinear stability of Bingham fluid Poiseuille flows in pipes and plane channels. These results show that the critical Reynolds number for transition, Re_c, increases with Bingham number, B, at least as fast as Re_cB^1/2 as B→∞. Estimates for the rate of increase are also provided. We compare these bounds and existing linear stability bounds with predictions from a series of phenomenological criteria for transition, as B→∞, concluding that only Hanks [AIChE J. 9 (1963) 306; 15 (1) (1963) 25] criteria can possibly be compatible with the theoretical criteria as B→∞. In the more practical range of application, 0≤B≤50, we show that there exists a large disparity between the different phenomenological criteria that have been proposed.     

Two new methods for the approximate determination of the initial coefficient of the first normal stress difference

Two methods for determining the initial coefficient of the first normal stress difference are presented. They are based on the evaluation of the steady viscosity function η($\text{γ}\text{.}$) and the viscosity function η⁺($\text{γ}\text{.}$, t) at the start-up of a flow with a very small rate of deformation $\text{γ}\text{.}$ < $\text{γ}\text{.}$₀. For the functions η($\text{γ}\text{.}$) and η⁺($\text{γ}\text{.}$), equations are given which can be used for a simple evaluation of the integral relationships obtaiend for ψ₁₀. The values for ψ₁₀ calculated by the two methods are compared with values obtained by the well-known methods via measurement of the ψ₁($\text{γ}\text{.}$) or η″(ω)/ω functions and extrapolation to zero). Both methods give values which are in satisfactory agreement with the experimental values.     

Plane-strain crack tip field and a dual-parameter fracture-criterion for non-linear fracture mechanics

Hancock and Cowling measured the critical crack tip opening displacements, δ_f, at fracture initiation in HY-80 steel specimens of six different configurations. δ_f varied from 90 μm in a deeply double-edge-cracked tensile panel to 900 μm in a single-edge-cracked tensile panel.McMeeking and Parks, and Shih and German have shown by their finite element calculations that the characteristics of the plane strain crack tip fields in both large scale yielding and general yielding are strongly dependent on specimen geometry and load level.In this study, the plane strain crack tip fields in the specimens tested by Hancock and Cowling were calculated using the finite element method. The crack tip triaxial tensile stress field is strongly affected by specimen geometric constraint, and the state of the triaxial tensile stress in a crack tip region is monitored by the ratio between the local tensile stress and the effective stress, i.e., ( ), at a distance x=2δ from the crack tip. The values of ( ) vary from 3.1 for the double-edge-cracked tensile panel to 1.7 for the single-edge-cracked tensile panel. The δ_f measured by Hancock and Cowling correlates very well with the ratio ( ). δ_f is a measure of the fracture ductility of the material ahead of the crack tip, and the ductility decreases with an increase in the triaxial tensile stress, i.e., the ratio ( ).     

Flow structure of wake behind a rotationally oscillating circular cylinder

Flow around a circular cylinder oscillating rotationally with a relatively high forcing frequency has been investigated experimentally. The dominant parameters affecting this experiment are the Reynolds number (Re), oscillation amplitude (θ_A), and frequency ratio F_R=f_f/f_n, where f_f is the forcing frequency and f_n is the natural frequency of vortex shedding. Experiments were carried out under conditions of Re=4.14×10³, 0°θ_A60° and 0.0F_R2.0. Rotational oscillation of the cylinder significantly modified the flow structure in the near-wake. Depending on the frequency ratio F_R, the cylinder wake showed five different flow regimes, each with a distinct wake structure. The vortex formation length and the vortex shedding frequency were greatly changed before and after the lock-on regime where vortices shed at the same frequency as the forcing frequency. The lock-on phenomenon always occurred at F_R=1.0 and the frequency range of the lock-on regime expanded with increasing oscillation amplitude θ_A. In addition, the drag coefficient was reduced when the frequency ratio F_R was less than 1.0 (F_R<1.0) while fixing the oscillation amplitude at θ_A=30°. When the oscillation amplitude θ_A was used as a control parameter at a fixed frequency ratio F_R=1.0 (lock-on regime), the drag reduction effect was observed at all oscillation amplitudes except for the case of θ_A=30°. This type of active flow control method can be used effectively in aerodynamic applications while optimizing the forcing parameters.     

Shear layers of a circular cylinder with rotary oscillation

The behavior of the separated shear layers and the near wake of a circular cylinder with small-amplitude rotary oscillations (Ω₁ = 0.05−0.15 for f _f/f _o ≤ 1.25) were investigated experimentally at Re = 3,700. Measurements of an unforced cylinder were also made for 2,000 ≤ Re ≤ 10,000 to better understand the effects of rotary oscillations. The results show that the shear-layer vortices formed closer to the cylinder and the distance separating them was found to decrease with cylinder oscillations. The shear-layer frequency, however, increased with increasing forcing frequency f _f. The formation-region length l _f decreased significantly with increasing f _f while decreased to a lesser extent with increasing normalized oscillation amplitude Ω₁. The shear layer also diffused to a length L _d larger than that of an unforced cylinder, while the l _f-L _d-Strouhal frequency offsetting mechanism was generally maintained. The near wake was of lower momentum compared to an unforced cylinder, and the transverse velocity fluctuations associated with the unforced vortex-shedding frequency f _o always presented a local peak at f _f/f _o = 0.5, regardless of Ω₁ tested.     

Polar plasticity of thin-walled composites made of ideal fibre-reinforced material

The present study investigates the influence that polar material response has on the plastic behaviour of thin-walled structures made of ideal fibre-reinforced materials (Spencer, 1972); or, equivalently, on the response of thin-walled fibrous composites within the first branch of the matrix dominated form (MDM) of the bimodal theory of plasticity (Soldatos, 2011, Dvorak and Bahei-El-Din, 1987). The plasticity studies mentioned above assume that fibres are infinitely thin and, therefore, perfectly flexible. They possess no bending stiffness and, hence, their negligible bending resistance cannot influence the developed stress state, which is accordingly described by a symmetric stress tensor. In contrast, the present study considers that if fibres resistant in bending are embedded in a material at high volume concentrations, their flexure produces couple-stress and, as a result of this kind of polar material behaviour, the stress tensor becomes non-symmetric. Under plane stress conditions that dominate behaviour of thin-walled structures, the stress-space and, therefore, conditions of plastic yield and relevant yield surfaces are thus four-dimensional. However, shapes and properties of initial yield surfaces relevant to the f₁-branch of MDM are studied comprehensively by considering their projection on particular planes of such a four-dimensional stress-space. It then becomes easier understood that, in the regime of polar material response, a thin-walled structure made of ideal fibre-reinforced material deforms plastically when suitable combinations of shear stress values are reached simultaneously, rather than when only one of two unequal shear stress components reaches some maximum absolute value. Thus, polar material plasticity dismisses the conventional concept of material yield stress in shear and replaces it with a pair of two independent yield moduli. Existence of the latter is perceived as a theoretical justification of the expectation that, due to the presence of fibres, two rather than one shear yield parameters of the composite should be present and accountable for. The non-zero values of those parameters are shown to exert paramount influence on the form of the yield surface of the ideal fibre-reinforced material of interest.     

Sensitivity of probabilistic fracture mechanics analysis to variations in statistical parameters

Probability of failure (p_f) of a structure is usually calculated for a specified set of statistical parameters (mean, standard deviation, and probability distribution) that characterize random variables. This approach may not be efficient in cases where one would like to know the effect of variations in statistical parameters on the probability of failure. A method based on generating and analyzing randomly selected statistical parameters is proposed. The method consists of generating databases of mean and coefficient of variation (COV = mean/standard deviation) values of relevant fracture mechanics variables through a random process. The method was applied to surface cracks in a flat wide plate loaded under elastic conditions. Probability of failure was calculated for each database record using the first-order reliability method (FORM). Multiple linear regression analyses of the database records were performed with p_f as dependent variable and statistical parameters as independent variables. The predicted p_f values were in very good agreement with the directly calculated p_f values for the specified variations of statistical parameters (±10%, ±15%, and ±20%), except those for fracture toughness and tensile stress, where variations should be limited to ±10% and ±15% ranges.     

Estimating the convergence rate for eigenfrequencies of anisotropic plates with variable thicknessEstimation du taux de convergence des valeurs propres des plaques anisotropes,avec épaisseur variable

Estimates of the differences between rescaled eigenvalues of the spectral problem for a thin anisotropic plate and eigenvalues of its two-dimensional models are obtained with bounds expressed in terms of the plate's thickness and attributes of the limit eigenvalue. To cite this article: S.A. Nazarov, C. R. Mecanique 330 (2002) 603–607.     

Bifurcations of Heteroclinic Orbits

The search for traveling wave solutions of a semilinear diffusion partial differential equation can be reduced to the search for heteroclinic solutions of the ordinary differential equation ü − cu̇ + f(u) = 0, where c is a positive constant and f is a nonlinear function. A heteroclinic orbit is a solution u(t) such that u(t) → γ ₁ as t → −∞ and u(t) → γ ₂ as t → ∞ where γ ₁, γ ₂ are zeros of f. We study the existence of heteroclinic orbits under various assumptions on the nonlinear function f and their bifurcations as c is varied. Our arguments are geometric in nature and so we make only minimal smoothness assumptions. We only assume that f is continuous and that the equation has a unique solution to the initial value problem. Under these weaker smoothness conditions we reprove the classical result that for large c there is a unique positive heteroclinic orbit from 0 to 1 when f(0) = f(1) = 0 and f(u) > 0 for 0 < u < 1. When there are more zeros of f, there is the possibility of bifurcations of the heteroclinic orbit as c varies. We give a detailed analysis of the bifurcation of the heteroclinic orbits when f is zero at the five points −1 < −θ < 0 < θ < 1 and f is odd. The heteroclinic orbit that tends to 1 as t → ∞ starts at one of the three zeros, −θ, 0, θ as t → −∞. It hops back and forth among these three zeros an infinite number of times in a predictable sequence as c is varied.     

Quasi-static behavior of Mg-alloys with and without short-fiber reinforcement

Magnesium alloys AE42 and AZ91 reinforced with 23 vol.% carbon short fibers (D_f ≈ 7 μm, L_f ≈ 100 μm) were tested under quasi-static loading. The carbon fibers were quasi-isotropically distributed in the horizontal plane (reinforced plane) of the casting. Compression and tensile tests were carried out on both the matrix alloys and the composites at temperatures between 20 °C and 300 °C. Specimens were machined to be loaded either parallel or normal to the reinforced plane. Due to the reinforcement, the compression yield stress of the composite AE42-C increased to a value approximately three-fold greater than the yield strength of the matrix; for composite AZ91-C this parameter was approximately 2.5-fold greater than that of the AZ91 matrix. The improvement in tensile strength was less than that in compression, which could be related to early tensile fracture through decohesion at the matrix–fiber interface, as detected by SEM investigations conducted on failed tensile specimens. Flow curves for the matrix alloys at different temperatures were described by a modified Kocks–Mecking material law. An idealization of a 2-D mesomodel was used for finite-element simulation of the mechanical behavior of the composites. The fibers were first considered as elastic bodies and the behavior of the matrix material was set according to the material law determined from the flow curves for the matrix alloys. Other calculations were carried out by considering elasto-plastic behavior of the fibers for application of a failure initiation technique to simulate the behavior of the composite materials beyond the ultimate stress.     

A strength reliability model by Markov process of unidirectional composites with fibers placed in hexagonal arrays

This paper proposes a strength reliability model based on a Markov process for unidirectional composites with fibers in a hexagonal array. The model assumes that a group of fiber breaking points, a so-called cluster, evolves with increased stress. The cluster evolution process branches because of various fiber-breakage paths. Load-sharing structure of intact fibers around clusters was estimated from geometric and mechanical local load-sharing rules. Composites fracture if a cluster achieves a critical size, so the model expresses a fracture criterion by setting an absorbing state. Next, the author constituted a state transition diagram concerning cluster evolutions of 1-fiber to 7-fiber breaks and analytically solved simultaneous differential equations obtained from the diagram. Results showed that, as critical cluster size increases, slope of the fracture probability distribution is given in a Weibull probability scale as follows: m_c=i×m_f (i, the number of broken fibers in a cluster; m_c and m_f, Weibull shape parameters for fracture probabilities of a critical cluster and fiber strength, respectively). This relation between m_c and m_f had been shown by Smith et al. [Proc. R. Soc. London, A 388 (1983) 353–391], but the present study demonstrated it analytically without any lower tail of the Weibull distribution used in that paper. In addition, the present model can be approximated by a one-state birth model.     

Bulk modulus of particulate composites: A comparison between a statistical and micro-mechanical approach

Two models are compared. One is based on the theory of elastic continua, and describes the interaction between filler and matrix in terms of an interfacial layer of varying volume fraction and elastic properties. The other derives from an equation of state for the constituents and the composite, based on molecular considerations. The filler-matrix interaction is then expressed in terms of segmental attractions and repulsions. We examine the dependence of the bulk modulusK _c ( _f ) on the volume fraction _f of filler and then show the correspondence between the two theories in terms of the infinite dilution limit of the ratio [K _c ( _f ) –K _m ]/(K _{m f} ) where the indexm refers to the matrix.Dedicated to Prof. Dr. F. R. Schwarzl on the occasion of his 60th birthdayOn leave from Rajdhani College, University of Delhi, India     

Rheology of suspensions of spherical particles in a newtonian and a non-newtonian fluid

Properties of suspensions of spherical glass beads (25–38 μm dia.) in a Newtonian fluid and a non-Newtonian (NBS Fluid 40) fluid were measured at volume fractions, φ, of 0%, 10%, 20% and 30%. Measurements were made using a modified and computerized Weissenberg Rheogoniometer. Properties measured included steady shear viscosity, η($\text{γ}\text{.}$), first normal stress difference, N₁($\text{γ}\text{.}$), linear viscoelastic properties, η′(ω) and G′(ω), shear stress relaxation, σ^? ($\text{γ}\text{.}$, t), and growth, σ⁺($\text{γ}\text{.}$, t) and normal stress relaxation, N₁^?($\text{γ}\text{.}$, t).For a the Newtonian fluid, increasing φ causes both η and η′ to increase, with η′ showing a slight frequency dependence. Both N₁ and G′ are zero and stress relaxation and growth occur essentially instantaneously. For the NBS fluid, both η and η′ increse with φ at all $\text{γ}\text{.}$ and ω, respectively, the increase being greater as $\text{γ}\text{.}$ and ω approach zero. N₁ and G′ are less affected by the presence of the particles than η and η′ with the effect on G′ being more pronounced than on N₁. For fixed $\text{γ}\text{.}$, stress relaxation and growth exhibit greater non-linear effects as φ is increased. A model for predicting a priori the linear viscoelastic properties for suspensions was found to yeild reasonable estimates up to φ = 20%.     

Collision of water drops with a plane water surface

The dependence of the outcome of the collision of uncharged water drops with a plane water surface on the impact angle , the velocity v₁ and the radius r₁ of the drops has been investigated experimentally. The impact parameters were varied over the intervals: v₁=0.40–1.05 m/sec, r₁=75–150m, and =16-85°. The method employed made it possible to avoid having to monitor the individual high-speed impact process. A stream of drops, produced in a vibrating reed type monodisperse droplet generator, was directed at the target. The impact parameters were measured by means of pulsed illumination. The results are expressed in the form of the dependence of the rebound probability and the coalescence coefficient E_S on the impact parameters. The existence of alternating conditional rebound-coalescence-rebound zones for different impact angles is established, together with a decrease in E_S with increase in r₁ and v₁. The data obtained generalize the results of previous experiments.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 165–168, May–June, 1990.     

Lock-on characteristics of a cavity shear layer

This investigation focuses on defining the lock-on regions of a cavity shear layer subject to local periodic excitations. A circular cylinder of small diameter (d=4 mm), located very close to the upstream edge of cavity, is used to generate the local periodic excitations in the form of oscillatory rotation about its center with various angular amplitudes (Δθ) and frequencies (f_e). All the experiments were conducted in a recirculating water channel at three different Reynolds numbers that are based on the momentum thickness at the upstream edge of cavity (Re_θ0=152, 216 and 278). The LDV system and the laser-sheet technique are employed to perform the quantitative velocity measurements and the qualitative flow visualization, respectively. For cavity flows at three Reynolds numbers studied, the resonant lock-on is found to be the primary lock-on region within the range of frequency ratio (f_e/f₀=0.28–2.0). Here f₀ denotes the natural instability frequency of an unexcited cavity shear layer. The frequency bandwidth of resonant lock-on region does increase with increasing excitation amplitudes (Δθ). While the excitation amplitudes are smaller than 5° (Δθ5°), the resonant lock-on region, at Reynolds numbers 216 and 278, distributes asymmetrically about f_e/f₀=1.0 and biases to the high frequency (or large f_e/f₀) side. However, the sidewise expansion of resonant lock-on region is enlarged and the degree of asymmetric distribution is alleviated at large excitation amplitudes (Δθ>5°). The amount of sidewise expansion of the resonant lock-on region biased toward the high-frequency side is more significant at the lowest Reynolds number (152) than those at two higher Reynolds numbers (216 and 278). Besides, there exists a sub-harmonic lock-on region only at the lowest Reynolds number 152. The existence of a sub-harmonic lock-on region clearly reveals that the differential equation governing the self-excited oscillation within a cavity contains the quadratic nonlinear term. Further, at the lowest Reynolds number (152), the sidewise expansion of the sub-harmonic lock-on region is much narrower than that of the resonant lock-on region.     

Reynolds' analogy for the thermal convection driven by nonisothermal stretching surfaces

In the present paper the steady boundary-layer flows induced by permeable stretching surfaces with variable temperature distribution are investigated under the aspect of Reynolds' analogy r = St _x /C _f(x). It is shown that for certain stretching velocities and wall temperature distributions, “Reynolds' function”r, i.e. the ratio of the local Stanton number St _x and the skin friction coefficient C _f(x) equals −1/2 for any value of the Prandtl number Pr and of the dimensionless suction/injection velocity f _w. In all of these cases, the dimensionless temperature field ϑ is connected to the dimensionless downstream velocity f ^′ by the simple relationship ϑ=(f ^′)^Pr. It is also shown that in the general case, Reynolds' function r may possess several singularities in f _w. The largest of them represents a critical value, so that for f _w<f _w,crit the solutions of the energy equation (although they still satisfy all the boundary conditions) become nonphysical.     

The interlayer shear effect on graphene multilayer resonators

Graphene nanostrips with single or few layers can be used as bending resonators with extremely high sensitivity to environmental changes. In this paper we report molecular dynamics (MD) simulation results on the fundamental and secondary resonant frequencies f of cantilever graphene nanostrips with different layer number n and different nanostrip length L. The results deviate significantly from the prediction of not only the Euler-Bernoulli beam theory (f∝nL⁻²), but also the Timoshenko's model. Since graphene nanostrips have extremely high intralayer Young's modulus and ultralow interlayer shear modulus, we propose a multibeam shear model (MBSM) that neglects the intralayer stretch but accounts for the interlayer shear. The MBSM prediction of the fundamental and secondary resonant frequencies f can be well expressed in the form f−f_mono∝[(n-1)/n]^bL^−2(1−b), where f_mono denotes the corresponding resonant frequency as the layer number is 1, with b=0.61 and 0.77 for the fundamental and secondary resonant modes. Without any additional parameters fitting, the prediction from MBSM agrees excellently with the MD simulation results. The model is thus of importance for designing multilayer graphene nanostrips based applications, such as resonators, sensors and actuators, where interlayer shear has apparent impacts on the mechanical deformation, vibration and energy dissipation processes therein.