Stress intensity factors for two offset parallel circular cracks: Part II — Semi-infinite solid

With the aid of the formulation in [1] (R. Muki, Progress in Solid Mechanics (North-Holland, 1961)) for general three-dimensional asymmetric problems and the superposition principle, Part II of this work makes use of the method in Part I (G.A.C. Graham and Q. Lan, J. Theor. Appl. Fract. Mech. 20, 207–225 (1994) [2]) to examine the interaction of arbitrarily located penny-shaped cracks in an infinite elastic solid to the case of a semi-infinite solid. As in Part I for the infinite body, the problem of a semi-infinite solid containing two penny-shaped cracks is reduced to a system of Fredholm integral equations of the second kind. These integral equations are then solved for some special cases when cracks are far apart and far away from the boundary. Some asymptotic solutions are presented and comparisons are made with the results for the special case where there is only one crack under axisymmetric loading.     

Stress intensity factors for two offset parallel circular cracks: Part III — Elastic layer

Part III of this work is concerned with the interaction of two penny-shaped cracks in the mid-plane of an elastic layer. Two cases, namely the stress free boundary case and the fixed boundary case, are considered. It is shown that these two cases are mathematically equivalent. As in Part I (G.A.C. Graham and Q. Lan, J. Theor. Appl. Fract. Mech. 20, 207–225 (1994) [1])_for the problems of an infinite solid and Part II (G.A.C. Graham and Q. Lan, J. Theor. Appl. Fract. Mech. 20, 227–237 (1994) [2]) of a semi-infinite solid, the problem is reduced to a system of Fredholm integral equations of the second kind. These integral equations are then solved when the crack size is small compared to the distance between them and the cracks are far away from the boundaries. It is also shown that the problem decouples when the cracks are subjected to normal loading and shear loading. Asymptotic solutions are presented for these two loadings.     

On a moderate rotation theory of laminated anisotropic shells—Part 1. Theory

A moderate rotation theory of laminated anisotropic shells, proposed by Schmidt and Reddy [J. appl. Mech. 55, 611–617.1988], is developed and its application is presented. All aspects of the derivations are explicitly developed and specific forms of the equations are derived in this part. The finite-element formulation and its applications are presented in Part 2 of the paper.     

Stress intensity factors for two offset parallel circular cracks: Part I — Infinite elastic solid

The method (W.D. Collins, Proc. R. Soc. London A274, 507–528 (1963); W.S. Fu and L.M. Keer, Int. J. Eng. Sci. 7, 361–372 (1969) [1,2]) used to solve co-planar penny-shaped cracks is generalized to investigate interaction of arbitrarily located penny-shaped cracks. The solution (M.K. Kassir and G.C. Sih, Three-dimensional Crack Problems (Noordhoff International, 1975) [3]) for the problem of an isolated crack in an infinite solid is applied together with the superposition principle to reduce the problem to a system of Fredholm integral equations of the second kind. These integral equations are then solved iteratively when the cracks are far apart. Some asymptotic solutions for the stress intensity factors are presented and comparisons are made whenever possible. Numerical solutions reveal some interesting phenomena.     

Solutions of some classes of second order non-linear ordinary differential equations

A second order non-linear ordinary differential equation satisfied by a homogeneous function of u and v where u is a solution of the linear equation ÿ + p(t)ÿ + r(t)y = 0 and v = ωu, ω being an arbitrary function of t, is obtained. Defining ω suitably in two specific cases, solutions are obtained for a non-linear equation of the form ÿ + p(t)ÿ + q(t)y = μÿ²y^{−1 +} f(t)yⁿ where μ ≠ 1, n≠ 1. Applying our results, some classes of equations of the above type possessing solutions involving two or one or no arbitrary constants are derived. Some illustrative examples are also discussed.     

Numerical solutions of pulsating flow and heat transfer characteristics in a pipe

Numerical studies are made of flow and heat transfer characteristics of a pulsating flow in a pipe. Complete time-dependent laminar boundary-layer equations are solved numerically over broad ranges of the parameter spaces, i.e., the frequency parameter β and the amplitude of oscillation A. Recently developed numerical solution procedures for unsteady boundary-layer equations are utilized. The capabilities of the present numerical model are satisfactorily tested by comparing the instantaenous axial velocities with the existing data in various parameters. The time-mean axial velocity profiles are substantially unaffected by the changes in β and A. For high frequencies, the prominent effect of pulsations is felt principally in a thin layer near the solid wall. Skin friction is generally greateer than that of a steady flow. The influence of oscillation on skin friction is appreciable both in terms of magnitude and phase relation. Numerical results for temperature are analyzed to reveal significant heat transfer characteristics. In the downstream fully established region, the Nusselt number either increases or decreases over the steady-flow value, depending on the frequency parameter, although the deviations from the steady values are rather small in magnitude for the parameter ranges computed. The Nusselt number trend is amplified as A increases and when the Prandtl number is low below unity. These heat transfer characteristics are qualitatively consistent with previous theoretical predictions.     

Rarefactive solitary waves in two-phase fluid flow of compacting media

Rarefactive solitary wave solutions of a third order nonlinear partial differential equation derived by Scott and Stevenson (Geophys. Res. Lett. 11, 1161–1164 (1984)) to describe the one-dimensional migration of melt under the action of gravity through the Earth's mantle are investigated. The partial differential equation contains two parameters, n and m, which are the exponents in power laws relating, respectively, the permeability of the medium and the bulk and shear viscosities of the solid matrix to the voidage. It is proved that, for any value of m, rarefactive solitary wave solutions satisfying certain physically reasonable boundary conditions always exist ifn>1 but do not exist if 0n1. It is also proved that the speed of the solitary wave is an increasing function of the amplitude of the wave. Six new exact rarefactive solitary wave solutions, four of which are expressed in terms of elementary functions and two in terms of elliptic integrals, are derived for six sets of values of n and m. The large amplitude approximation is considered and the results of Scott and Stevenson for n>2, m=0 and n>1, m=1 are extended to n>1 and all m0. It is shown that, for sufficiently large amplitude, larger amplitude solitary waves are broader in width if 0m1 and are narrower in width if m>1.     

Shear-free turbulent diffusion—classical and new scaling laws

We consider the problem of turbulence generation at a vibrating grid in the x₂–x₃ plane. Turbulence diffuses in the x₁-direction. Analyzing the multi-point correlation equation using Lie-group analysis, we find three different invariant solutions (scaling laws): classical diffusion-like solution (heat equation like), decelerating diffusion-wave solution and finite domain diffusion due to rotation. All solutions have been obtained using Lie-group (symmetry) methods. It is shown that if only one spatial dimension is considered, models based on Reynolds averaging are only capable to model either the diffusion-like solution or the decelerating diffusion-wave solution. The latter solution is only admitted under certain algebraic constraints on the model constants; e.g. in case of the K– model the model constants need to obey the relation c₂σ/σ_K=2. Turbulent diffusion on a finite domain induced by rotation is not admitted by any of the classical models. Finally, in the appendix it is shown that Lele's transformation (Phys. Fluids 28(1) (1985) 64) leads to a complete analytic solution of the steady diffusion problem modelled by the K– equation.     

A Navier-Stokes code for S-shaped diffusers—a review

Three-dimensional, compressible, internal flow solutions obtained using a thin-layer Navier-Stokes code are presented. The code, formulated by P.D. Thomas, is based on the Beam-Warming implicit factorization scheme; the boundary conditions also are formulated implicitly. Turbulent flow is treated through the use of the Baldwin-Lomax two-layer, algebraic eddy viscosity model. Steady-state solutions are obtained by solving numerically the time-dependent equations from given initial conditions until the time-dependent terms become negligible. The configuration considered is a rectangular cross-section, S-shaped centreline diffuser duct with an exit/inlet area ratio of 2.25. The Mach number at the duct entrance is 0.9, with a Reynolds number of 5.82 × 10⁵. Convergence to the final results required about 2700 time steps or 11 hours of CPU time on our CRAY-1M computer. The averaged residuals were reduced by about two orders of magnitude during the computations. Several regions of separated flow exist within the diffuser. The separated flow region on the upper wall, downstream of the second bend, is by far the largest and extends to the exit plane.     

On a Korteweg–de Vries-like equation with higher degree non-linearity

We deal with a non-linear partial differential equation which has been widely investigated owing to its applications in quantum field theory, as well as plasma and solid-state physics. It is the matter of a third order KdV-like equation with higher degree non-linearity in the coefficient of the transport term; it can be derived from a Lagrangian or an Hamiltonian density. In the current literature specific attention has been devoted to the search for traveling-wave solutions, depending upon a positive parameter v, which assesses the speed of the solitary wave. The velocity v is always assumed to be constant, as its dependence on the wave-amplitude is neglected in the mathematical model. In this context, Coffey [On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592] exploits an algebraic recursive technique to obtain these solutions in closed form for particular values of v. The aim of this paper is to extend these results by showing that closed-form solutions are achievable for every value of v: to this purpose we supply a proper mathematical framework for these issues by taking into account a suitable special function, namely an elliptic function in the sense of Weierstraß. Furthermore we obtain two classes of the so-called kink solutions, see [M.W. Coffey, On series expansions giving closed-form solutions of Korteweg–de Vries-like equations, SIAM J. Appl. Math. 50 (6) (1990) 1580–1592; B. Dey, Domain wall solutions of KdV-like equations with higher order non-linearity, J. Phys. A 19 (1) (1986) L9–L12], as well as an exponential development of the general solution, for which we prove the convergence. Eventually we show how to implement the resulting functions by means of a symbolic manipulation program.     

A comparative study on low‐order,amplitude‐equation and perturbation approaches in thermal convection

A comparison among three weakly nonlinear approaches for thermo‐gravitational instability in a Newtonian fluid layer heated from below is presented. First, the dynamical systems describing the time evolution of the problem from different weakly nonlinear approaches, namely, the Lorenz model, the amplitude equations and the perturbation expansion approaches are obtained. Next, the steady states and their stability, as well as the transient behaviour are obtained from each dynamical system. The similarity and difference among the three models are emphasized. The role of each of the nondimensional groups, the Rayleigh number and the Prandtl number is compared for the three models. The different approaches lead to similar behaviours when the Rayleigh number is just above its critical value and Prandtl number is high. However, only the dynamical system obtained from the amplitude equations is able to reflect the role of the Prandtl number. On the other hand, the amplitude equations and perturbation expansion techniques are not suitable for predicting the uniform oscillatory behaviour observed frequently in Rayleigh–Bénard convection. The novelty of the current work lies in studying the critical differences in the findings of the three popular approaches to investigate weakly nonlinear thermal convection for the first time. Copyright © 2011 John Wiley & Sons, Ltd.     

Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure

A three‐dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds‐averaged Navier–Stokes equations with a non‐hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co‐ordinate system, with a semi‐implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five‐diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non‐hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind‐induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

Experimental investigation of cylindrically diverging solitons in an electric lattice

The results of an experimental investigation of cylindrical solitons in a two-dimensional electric LC-lattice are given. It is shown that in the continuum limit, propagation of cylindrical waves far from the center of symmetry in such a lattice may be described for each ray tube by a known modification of the Korteweg-de Vries equation which takes account of the cylindrical divergence. The dispersion term in this equation depends on the direction of wave propagation relative to the direction of the main axes of the lattice. Formation of solitons from non-soliton-shaped pulses was observed. The variations of soliton amplitude and duration with distance have been determined. They agree well with the numerical calculations by Maxon & Viecelli [2] and Dorfman [9]. Comparison of the obtained experimental data with the known theoretical laws of amplitude attenuation for diverging solitons [2, 12, 14] seems to favor the validity of the law A r^−2/3 rather than A r^−1/2.     

Numerical studies of slow viscous rotating flow past a sphere—1

The Navier-Stokes equations for a steady, viscous rotating fluid, rotating about the z-axis with angular velocity ω are linearized using the Stokes approximation. The linearized Navier-Stokes equations governing the axisymmetric flow can be written as three coupled partial differential equations for the stream function, vorticity and rotational velocity components. One parameter, R_eω = 2ωa²/v, enters the resulting equations. For R_eω « 1, the coupled equations are solved by the Peaceman-Rachford A.D.I. (Alternating Direction Implicit) method and the resulting algebraic equations are solved by the ‘method of sweeps’. Stream lines for ψ = 0·05, 0·2, 0·5 and magnitude of the vorticity vector z = 0·2 are plotted for R_eω = 0·1, 0·3, 0·5. Correction to the Stokes drag due to the rotation of fluid is calculated.     

Applications of the integral equation method to delay differential equations

In this paper, we compare two approaches for determining the amplitude equations; namely, the integral equation method and the method of multiple scales. To describe and compare the methods, we consider three examples: the parametric resonance of a Van der Pol oscillator under state feedback control with a time delay, the primary resonance of a harmonically forced Duffing oscillator under state feedback control with a time delay, and the primary resonance together with 1:1 internal resonance of a two degree-of-freedom model. Using the integral equation method and the method of multiple scales, the amplitude equations are obtained. The stability of the periodic solution is examined by using the Floquet theorem together with the Routh–Hurwitz criterion (without time delay) and the Nyquist criterion (with time delay). By comparison with the solution obtained by the numerical integration, we find that the accuracy of the integral equation method is much better.     

Predicting turbulent flow in a staggered tube bundle

This paper presents the results of calculations performed for the turbulent, incompressible flow around a staggered array of tubes for which carefully obtained experimental results are available as part of an established ERCOFTAC-IAHR test case. The Reynolds-averaged Navier–Stokes equations are solved using a pressure-based finite volume algorithm, using collocated cell vertex store on an unstructured and adaptive mesh of tetrahedra. Turbulence closure is obtained with a truncated form of a low-Reynolds number k– model developed by Yang and Shih. The computational domain covers all seven rows of tubes used in the experimental study and periodic flow is allowed to develop naturally. The results of the computations are surprisingly good and compare favourably with results obtained by others using a wide range of alternative k– models for a single cylinder with periodic inflow and outflow boundaries on structured meshes.     

An iterative method for the computation of the response of linearised Navier–Stokes equations to harmonic forcing and application to forced cylinder wakes

The central aim of this paper is the development and application of an efficient, iterative methodology for the computation of the perturbation fields induced by harmonic forcing of the linearised Navier–Stokes equations. The problem is formulated directly in the frequency domain, and the resulting system of equations is solved iteratively until convergence. The method is easily implemented to any implicit code that can solve iteratively the steady‐state Navier–Stokes equations. In this paper, it is applied to investigate the flow around a static cylinder with pulsating approaching flow and a cylinder undergoing forced stream‐wise oscillations. All terms of the perturbation kinetic energy equation are computed, and it is shown that perturbations grow by extracting energy from two sources: the underlying base flow field and the externally provided energy that maintains the imposed oscillation. The periodic drag force acting on the cylinder is also computed, and it is demonstrated that Morrison's equation is a simple model that can estimate with good accuracy the amplitude and phase of this force with respect to the approaching flow. The perturbation fields induced by periodic inlet flow (static cylinder) and forced stream‐wise cylinder oscillation are closely related: the velocity fields are identical in the appropriate reference frames, and a simple expression is derived, which links the pressures in the two flow cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Renormalization group theory for turbulence: Assessment of the Yakhot-Orszag-Smith theory

The purpose of this paper is to estimate the renormalization group theory for turbulence developed by Yakhot and Orszag [J. Sci. Comput. 1 (1986) 3] and reformulated by Yakhot and Smith [J. Sci. Comput. 7 (1992) 35]. We go into details of their basic theory for the Navier-Stokes equations, the transport equations for the turbulent kinetic energy K, and its dissipation rate . As a result, it becomes evident that their theory bears no relationship to Wilson's renormalization group theory for critical phenomena. Their model is not directly obtained from the renormalization group theory. They evaluated the Kolmogorov constant by setting the expansion parameter ε = 0 and ε = 4 in the same equations. Furthermore, all the constants in their model are invalid because of the same problem.     

三阶模态耦合效应下悬索的1:1主共振与稳定性分析

针对悬索的振动,研究了模态耦合效应对悬索振动特征的影响。首先基于哈密顿原理推导了考虑抗弯刚度影响的悬索的偏微分振动方程,采用Galerkin方法得到了悬索的前三阶模态耦合振动常微分方程组。采用多尺度法分析了悬索的一阶、二阶和三阶主共振,得到了一阶、二阶和三阶主共振的幅-频响应方程,接着基于Lyapunov稳定性理论进行了稳定性分析,最后进行了数值算例分析。算例分析表明,当1:1主共振发生时,一阶主共振产生的幅值远大于二阶和三阶主共振产生的幅值,即当悬索振动时,能量主要以一阶模态幅值的形式散发;在同阶次幅值-σ曲线中,随着F的增加,1:1主共振产生的幅值有所增加;在幅值-V曲线中,随着σ的增加,临界跳跃点有向右偏移的趋势,σ增加会导致幅值增加;档距越大,一阶、二阶和三阶1:1主共振产生的幅值越大,但一阶主共振产生的幅值增加最为明显。