Reactive-diffusive equation with variable pre-exponential factor

We consider the steady-state solutions for the exothermic chemical reaction, taking the diffusion of the reactant in a slab into account and assuming an Arrhenius temperature dependence with variable pre-exponential factor. We solve the resulting nonlinear boundary value problem using both numerical and analytical methods. This new analytical solution for the Frank-Kamenetskii parameter δ, as it is in terms of Bernoulli’s numbers, is in accordance with the numerical integration provided the activation energy parameter ε (1) is very small and for ε→0 it goes to the well-known Frank-Kamenetskii case. We determine numerically the transitional values of δ, ε and the dimensionless central temperature θ_m.     

The behavior of cracked cross-ply composite laminates under shear loading

An interlaminar-shear-stress analysis developed earlier by Tsai et al. (1990, Micro-cracking-Induced Damage in Composites) for a [φ_m/θ_n], bi-directional composite laminate is used to solve the case of a cross-ply [0_m/90_n]_x laminate with the 90° layer only or both layers cracked under pure shear loading. Strains, forces and laminate shear modulus reduction due to matrix cracking were obtained. Experimental results for shear modulus as a function of crack densities were obtained by a simple shear test and they agree very well with the theoretical prediction.     

MOTIONS, FORCES AND MODE TRANSITIONS IN VORTEX-INDUCED VIBRATIONS AT LOW MASS-DAMPING

These experiments, involving the transverse oscillations of an elastically mounted rigid cylinder at very low mass and damping, have shown that there exist two distinct types of response in such systems, depending on whether one has a low combined mass-damping parameter (low m^*ζ), or a high mass-damping (highm^*ζ ). For our low m^*ζ, we find three modes of response, which are denoted as an initial amplitude branch, an upper branch and a lower branch. For the classical Feng-type response, at highm^*ζ , there exist only two response branches, namely the initial and lower branches. The peak amplitude of these vibrating systems is principally dependent on the mass-damping (m^*ζ), whereas the regime of synchronization (measured by the range of velocity U^*) is dependent primarily on the mass ratio, m^*ζ. At low (m^*ζ), the transition between initial and upper response branches involves a hysteresis, which contrasts with the intermittent switching of modes found, using the Hilbert transform, for the transition between upper–lower branches. A 180° jump in phase angle φ is found only when the flow jumps between the upper–lower branches of response. The good collapse of peak-amplitude data, over a wide range of mass ratios (m^*=1–20), when plotted against (m^*+C_A) ζ in the “Griffin” plot, demonstrates that the use of a combined parameter is valid down to at least (m^*+C_A)ζ 0·006. This is two orders of magnitude below the “limit” that had previously been stipulated in the literature, (m^*+C_A) ζ>0·4. Using the actual oscillating frequency (f) rather than the still-water natural frequency (f_N), to form a normalized velocity (U^*/f^*), also called “true” reduced velocity in recent studies, we find an excellent collapse of data for a set of response amplitude plots, over a wide range of mass ratiosm^* . Such a collapse of response plots cannot be predicted a priori, and appears to be the first time such a collapse of data sets has been made in free vibration. The response branches match very well the Williamson–Roshko (Williamson & Roshko 1988) map of vortex wake patterns from forced vibration studies. Visualization of the modes indicates that the initial branch is associated with the 2S mode of vortex formation, while the Lower branch corresponds with the 2P mode. Simultaneous measurements of lift and drag have been made with the displacement, and show a large amplification of maximum, mean and fluctuating forces on the body, which is not unexpected. It is possible to simply estimate the lift force and phase using the displacement amplitude and frequency. This approach is reasonable only for very low m^*.     

Collinear Central Configurations and Singular Surfaces in the Mass Space

For a given m=(m₁,...,m_n)(R⁺)ⁿ, let p and q(R³)ⁿ be two central configurations for m. Then we call p and q equivalent and write pq if they differ by an SO(3) rotation followed by a scalar multiplication as well as by a permutation of bodies. Denote by L(n,m) the set of equivalent classes of n-body collinear central configurations in R³ for any given mass vector m=(m₁,...,m_n)(R⁺)ⁿ. The main discovery in this paper is the existence of a union H₃ of three non-empty algebraic surfaces in the mass half space (m₁,m₂–m₁,m₃–m₂)R⁺×R² besides the planes generated by equal masses, which decreases the number of collinear central configurations. The union H₃ in R⁺×R ² is explicitly constructed by three 6-degree homogeneous polynomials in three variables such that, for any mass vector m=(m₁,m₂,m₃)(R⁺)³, ^# L(3,m)=3, if m₁, m₂, and m₃ are mutually distinct and (m₁,m₂–m₁,m₃–m₂)H₃, ^# L(3,m)=2, if m₁, m₂, and m₃ are mutually distinct and (m₁,m₂–m₁,m₃–m₂)H₃, ^# L(3,m)=2, if two of m₁, m₂, and m₃ are equal but not the third, ^# L(3,m)=1, if m₁=m₂=m₃. We give also a sharp upper bound on ^#L(n,m) for any positive mass vector m(R⁺)ⁿ.     

Effect of material nonhomogeneities on the HRR dominance

This work studies the asymptotic stress and displacement fields near the tip of a stationary crack in an elastic–plastic nonhomogeneous material with the emphasis on the effect of material nonhomogeneities on the dominance of the crack tip field. While the HRR singular field still prevails near the crack tip if the material properties are continuous and piecewise continuously differentiable, a simple asymptotic analysis shows that the size of the HRR dominance zone decreases with increasing magnitude of material property gradients. The HRR field dominates at points that satisfy |α⁻¹ ∂α/∂x_δ|1/r, |α⁻¹ ∂²α/(∂x_δ ∂x_γ)|1/r², |n⁻¹ ∂n/∂x_δ|1/[r|ln(r/A)|] and |n⁻¹ ∂²n/(∂x_δ ∂x_γ)|1/[r²|ln(r/A)|], in addition to other general requirements for asymptotic solutions, where α is a material property in the Ramberg–Osgood model, n is the strain hardening exponent, r is the distance from the crack tip, x_δ are Cartesian coordinates, and A is a length parameter. For linear hardening materials, the crack tip field dominates at points that satisfy |E_tan⁻¹ ∂E_tan/∂x_δ|1/r, |E_tan⁻¹ ∂²E_tan/(∂x_δ ∂x_γ)|1/r², |E⁻¹ ∂E/∂x_δ|1/r, and |E⁻¹ ∂²E/(∂x_δ ∂x_γ)|1/r², where E_tan is the tangent modulus and E is Young’s modulus.     

Understanding the Piche evaporimeter as a simple integrating mass transfer meter

Results are presented of a comparison of measured and calculated evaporation rates of the Piche evaporimeter under indoor and outdoor (within a meteorological screen) conditions. In both cases, application of mass transfer formulae in use for horizontal (turbulent) flow to the evaporating blotting paper of the instrument yield very good results under pure forced convection conditions. For mixed convection regimes, comparisons using either pure free (combined heat and mass transfer) or pure forced convection equations give as expected too low calculated values. Reasons for such differences with measured values are reviewed. Our forced convection results confirm that main stream turbulence is only of influence on mass transfer to a zero incidence flow in combination with pressure gradient (bluff body) effects, which under our conditions appear to be absent around the Piche surfaces. The same results prove absence of any influence of the particular temperature distribution over the blotting paper on the mass transfer. The understanding and importance of these conclusions in relation to the use of the Piche evaporimeter as a simple integrating mass transfer meter under actual farming conditions are discussed. The importance to obtain such mass transfer data is explained in the introduction.Nomenclature A Numerical constant in free convection Sherwood number - Coefficient of thermal expansion (K^–1) - C _{(s, b)} Water vapour concentration average at the evaporating surface (s) and in the bulk air (b) (g m^–3) - D Coefficient of molecular diffusion of water vapour in air (m² s^–1) - d Characteristic dimension of the paper disc in the direction of flow (m) - E _{(c, m)} Evaporation rates of the Piche evaporimeter, calculated (c) and measured (m) (units in text) - e _{(s, b)} Partial water vapour pressure average at the evaporating surface (s) and in the bulk air (b) (mbar) - Gr Grashof number - g Acceleration of gravity (m s^–2) - m Number of measuring periods - n Numerical constant in free convection Sherwood number - Coefficient of kinematic viscosity of air (m² s^–1) - P Atmospheric pressure (mbar) - Re Reynolds number - _{(s, b)} Air density average at the evaporating surface (s) and of the bulk air (b) (g m^–3) - Sh Sherwood number - T _{(s, b)} Temperature average of the evaporating surface (s) and in the bulk air (b) (K) - T _{(vs, vb)} Virtual temperature average at the evaporating surface (s) and in the bulk air (b) (K) - U Wind speed (air movement) average of the bulk air (m s^–1)     

The anisotropy of simple shearing flow

Summary The optical analogue of the second normal stress difference, namely,n ₂₂ –n ₃₃ has been measured for several polymer systems. The experimental method is direct and makes use of capillary and slit dies. For the polymer melts examined — ranging from a polystyrene with a viscosity of 10⁴ Ns/m² to a PIB of viscosity 50 Ns/m² — a non-zero (n ₂₂ –n ₃₃) was found which could not be accounted for by any systematic errors. Except for one case, where (n ₁₁ –n ₂₂) could not be evaluated, the ratio (n ₂₂ –n ₃₃)/(n ₁₁ –n ₂₂) was found to be negative with an absolute value of 0.05 to about 0.1. For solutions the results are not clear-cut: in several casesn ₂₂ –n ₃₃ was immeasurably small but in one case (23% PIB in oil) the above ratio was the same as for the melts. A modification of the (dilute solution) molecular theories is suggested which is qualitatively able to describe most of the experimental results here outlined.

---

Zusammenfassung Die optische Analogie der zweiten Normalspannungs-Differenz, nämlichn ₂₂ –n ₃₃ wurde für verschiedene Polymersysteme gemessen. Dabei wurden Kapillaren und Schlitzdüsen benutzt. Bei allen untersuchten polymeren Schmelzen, die sich von einem Polystyrol mit einer Viskosität von 10⁴ Ns/m² bis zu einem Polyisobuten mit einer Viskosität von 50 Ns/m² erstreckten, wurde n₂₂ –n ₃₃ signifikant ungleich Null gefunden. Mit Ausnahme eines Falles, in dem (n ₁₁ –n ₂₂) nicht bestimmt werden konnte, ergab sich die Größe (n ₂₂ –n ₃₃)/(n ₁₁ –n ₂₂) negativ mit einem Absolutwert zwischen 0,05 und 0,1.Für Lösungen waren die Ergebnisse nicht so eindeutig: In mehreren Fällen warn ₂₂ –n ₃₃ unmeßbar klein, in einem Fall aber (23% PIB in Öl) war diese Größe von der gleichen Ordnung wie für die Schmelzen.Es wurde eine Modifikation der molekularen Theorie verdünnter Lösungen vorgeschlagen, die die meisten der hier dargestellten experimentellen Ergebnisse zu beschreiben gestattet.

---

---

With 12 figures     

Effect of transverse cracks on the mechanical properties of angle-ply composites laminates

A modified shear lag analysis, taking into account the notion of stress perturbation function, is employed to evaluate the effect of transverse cracks on the stiffness reduction in [±θ_n/90_m]_S angle-ply laminated composites. Effects of number of 90° layers and number of ±θ layers on the laminate stiffness have also been studied. The present results represent well the dependence of the degradation of mechanical properties on the fibre orientation angle of the outer layers, the number of cracked cross-ply layers and the number of uncracked outer ±θ layers in the laminate.     

Control of a submerged jet in a thin rectangular cavity

An innovative method is presented for control of an oscillatory turbulent jet in a thin rectangular cavity with a thickness to width ratio of 0.16. Jet flow control is achieved by mass injection of a secondary jet into the region above the submerged primary jet nozzle exit and perpendicular to the primary nozzle axis. An experimental model, a 2-D and a 3-D computational fluid dynamics (CFD) model are used to investigate the flow characteristics under various secondary injection mass flow rates and injection positions. Two-dimensional laser Doppler anemometry (LDA) measurements are compared with results from the CFD models, which incorporate a standard k–ε turbulence model or a 2-D and 3-D realisable k–ε model. Experimental results show deflection angles up to 23.3° for 24.6% of relative secondary mass flow are possible. The key to high jet control sensitivity is found to be lateral jet momentum with the optimum injection position at 12% of cavity width (31.6% of the primary nozzle length) above the primary nozzle exit. CFD results also show that a standard k–ε turbulence closure with nonequilibrium wall functions provides the best predictions of the flow.     

Heteroclinic Orbits and Bernoulli Shift¶for the Elliptic Collision Restricted¶Three-Body Problem

We consider two mass points of masses m ₁=m ₂= moving under Newton's law of gravitational attraction in a collision elliptic orbit while their centre of mass is at rest. A third mass point of mass m ₃≈ 0, moves on the straight line L, perpendicular to the line of motion of the first two mass points and passing through their centre of mass. Since m ₃≈ 0, the motion of the masses m ₁ and m ₂ is not affected by the third mass, and from the symmetry of the motion it is clear that m ₃ will remain on the line L. So the three masses form an isosceles triangle whose size changes with the time. The elliptic collision restricted isosceles three-body problem consists in describing the motion of m ₃. In this paper we show the existence of a Bernoulli shift as a subsystem of the Poincaré map defined near a loop formed by two heteroclinic solutions associated with two periodic orbits at infinity. Symbolic dynamics techniques are used to show the existence of a large class of different motions for the infinitesimal body. Accepted July 6, 2000?Published online February 14, 2001     

Double diffusive natural convection in a partially heated enclosure with Soret and Dufour effects

The effect of double-diffusive natural convection of water in a partially heated enclosure with Soret and Dufour coefficients around the density maximum is studied numerically. The right vertical wall has constant temperature θ_c, while left vertical wall is partially heated θ_h, with θ_h > θ_c. The concentration in right wall is maintained higher than left wall (C_c < C_h) for case I, and concentration is lower in right wall than left wall (C_h > C_c) for case II. The remaining left vertical wall and the two horizontal walls are considered adiabatic. Water is considered as the working fluid. The governing equations are solved by control volume method using SIMPLE algorithm with QUICK scheme. The effect of the various parameters (thermal Rayleigh number, center of the heating location, density inversion parameter, Buoyancy ratio number, Schmidt number, and Soret and Dufour coefficients) on the flow pattern and heat and mass transfer has been depicted. Comprehensive Nusselt and Sherwood numbers data are presented as functions of the governing parameters mentioned above.     

3-D kinematical conservation laws (KCL): Evolution of a surface in – in particular propagation of a nonlinear wavefront

3-D KCL are equations of evolution of a propagating surface (or a wavefront) Ω_t in 3-space dimensions and were first derived by Giles, Prasad and Ravindran in 1995 assuming the motion of the surface to be isotropic. Here we discuss various properties of these 3-D KCL. These are the most general equations in conservation form, governing the evolution of Ω_t with singularities which we call kinks and which are curves across which the normal n to Ω_t and amplitude w on Ω_t are discontinuous. From KCL we derive a system of six differential equations and show that the KCL system is equivalent to the ray equations of Ω_t. The six independent equations and an energy transport equation (for small amplitude waves in a polytropic gas) involving an amplitude w (which is related to the normal velocity m of Ω_t) form a completely determined system of seven equations. We have determined eigenvalues of the system by a very novel method and find that the system has two distinct nonzero eigenvalues and five zero eigenvalues and the dimension of the eigenspace associated with the multiple eigenvalue 0 is only 4. For an appropriately defined m, the two nonzero eigenvalues are real when m>1 and pure imaginary when m<1. Finally we give some examples of evolution of weakly nonlinear wavefronts.     

Symmetry and secondary bifurcation of elastic structures

We consider non-linear bifurcation problems for elastic structures modeled by the operator equation F[w;α]=0 where F:X×R^k→Y,X,Y are Banach spaces and XY. We focus attention on problems whose bifurcation equations are of the form

f_i(α₁,α₂;λ,μ)=(a_iμ+b_iλ)α_i+p_iα_i³+q_iα_i∑_j=1,j≠i^kα_j+1²+α_ih_i(λ,μ;α₁,α₂,…α_k) i=1,2,…k

which emanates from bifurcation problems for which the linearization of F is Fredholm operators of index 0. Under the assumption of F being odd we prove an important theorem of existence of secondary bifurcation. Under this same assumption we prove a symmetry condition for the reduced equations and consequently we got an existence result for secondary bifurcation. We also include a stability analysis of the bifurcating solutions.     

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.     

Oblique stagnation-point flow of a viscoelastic fluid with heat transfer

The two-dimensional forced convection stagnation-point flow and heat transfer of a viscoelastic second grade fluid obliquely impinging on an infinite plane wall is considered as an exact solution of the full partial differential equations. This oblique flow consists of an orthogonal stagnation-point flow to which a shear flow whose vorticity is fixed at infinity is added. The relative importance of these flows is measured by a parameter γ. The viscoelastic problem is reduced to two ordinary differential equations governed by the Weissenberg number W_e, two parameters α and β, the later being a free parameter β, introduced by Tooke and Blyth [A note on oblique stagnation-point flow, Physics of Fluids 20 (2008) 033101-1–3], and the Prandtl number Pr. The two cases when α=β and α≠β are, respectively, considered. Physically the free parameter may be viewed as altering the structure of the shear flow component by varying the magnitude of the pressure gradient. It is found that the location of the separation point x_s of the boundary layer moves continuously from the left to the right of the origin of the axes (x_s<0).     

Some plane adiabatic ideal gas flows containing shocks

We obtain the solution describing adiabatic flows of an ideal gas characterized by the two parameters a and b such that [a]=L ^m+1 T ^–1, [b]=ML ^–2–2m where m is arbitrary (m > 0).h This solution permits the construction of flows containing shocks.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 10, No. 3, pp. 71–73, May–June, 1969.     

The global flow of the hyperbolic restricted three-body problem

Two mass points of equal masses m ₁ = m ₂ > 0 move under Newton's law of attraction in a non-collision hyperbolic orbit while their center of mass is at rest. We consider a third mass point, of mass m ₃ = 0, moving on the straight line L perpendicular to the plane of motion of the first two mass points and passing through their center of mass. Since m ₃ = 0, the motion of m ₁ and m ₂ is not affected by the third, and from the symmetry of the motion it is clear that m ₃ remains on the line L. The hyperbolic restricted 3-body problem is to describe the moton of m ₃. Our main result is the characterization of the global flow of this problem.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

Crack repair using an elastic filler

The effect of repairing a crack in an elastic body using an elastic filler is examined in terms of the stress intensity levels generated at the crack tip. The effect of the filler is to change the stress field singularity from order 1/r^1/2 to 1/r^(1-λ) where r is the distance from the crack tip, and λ is the solution to a simple transcendental equation. The singularity power (1-λ) varies from (the unfilled crack limit) to 1 (the fully repaired crack), depending primarily on the scaled shear modulus ratio γ_r defined by G₂/G₁=γ_rε, where 2πε is the (small) crack angle, and the indices (1, 2) refer to base and filler material properties, respectively. The fully repaired limit is effectively reached for γ_r≈10, so that fillers with surprisingly small shear modulus ratios can be effectively used to repair cracks. This fits in with observations in the mining industry, where materials with G₂/G₁ of the order of 10^-3 have been found to be effective for stabilizing the walls of tunnels. The results are also relevant for the repair of cracks in thin elastic sheets.     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.