Aerodynamic flow investigations in an isothermal model of an afterburner

Experimental and numerical investigations of the three-dimensional flowfields in an isothermal model of an afterburner are presented. Five-hole pilot probe measurements along the entire length of the model in three different azimuthal planes, allow the determination of three mean velocity components which provide comprehensive information to aid understanding of such complex flows. The numerical calculations are performed using a SIMPLE based algorithm with staggered grid arrangement. The standardk-ɛ model is used for physical modeling. The numerical results agree quite satisfactorily with the time-mean velocity measurements. The predicted turbulence kinetic energy contours have also been presented.     

Void collapse in an elastic solid

Radial deformations of an infinite medium surrounding a traction-free spherical cavity are considered. The body is composed of an isotropic, incompressible elastic material and is subjected to a uniform pressure at infinity. The possibility of void collapse (i.e. the void radius becoming zero at a finite value of the applied stress) is examined. This does not occur in all materials. The class of materials that do exhibit this phenomenon is determined, and for such materials, an explicit expression for the value of the applied pressure at which collapse occurs is derived. The stability of the deformation and the influence of a finite outer radius are also considered. The results are illustrated for a particular class of power-law materials. In certain respects, the present results for void collapse are complementary to Ball (1982)'s results for cavitation in an incompressible elastic material.Some brief observations on void collapse in compressible materials are made. The collapse of a void under non-symmetric conditions is also discussed by utilizing a solution obtained by Varley and Cumberbatch (1977, 1980).The results reported here were obtained in the course of an investigation supported in part by the U.S. Army Research Office.     

MHD boundary-layer flow of an upper-convected Maxwell fluid in a porous channel

Two-dimensional magnetohydrodynamic (MHD) boundary layer flow of an upper-convected Maxwell fluid is investigated in a channel. The walls of the channel are taken as porous. Using the similarity transformations and boundary layer approximations, the nonlinear partial differential equations are reduced to an ordinary differential equation. The developed nonlinear equation is solved analytically using the homotopy analysis method. An expression for the analytic solution is derived in the form of a series. The convergence of the obtained series is shown. The effects of the Reynolds number Re, Deborah number De and Hartman number M are shown through graphs and discussed for both the suction and injection cases.     

Stability of a two-phase process in an elastic solid

This work investigates the linear stability of an antiplane shear motion which involves a steadily propagating normal planar phase boundary in an arbitrary element of a family of non-elliptic generalized neo-Hookean materials. It is shown that such a process is linearly unstable with respect to a large class of disturbances if and only if the kinetic response function—a constitutively supplied entity which relates the normal velocity of a phase boundary to the driving traction which acts on it—is locally decreasing as a function of the appropriate argument. This result holds whether or not inertial effects are taken into consideration, demonstrating that the linear stability of the relevant process depends entirely upon the transformation kinetics intrinsic to the kinetic response function. The morphological evolution of the interface is then, in an inertia-free setting, tracked for a short time subsequent to the perturbation. It is found that, when the kinetic response function is non-monotonic, the phase boundary can evolve so as to qualitatively resemble the plate-like structures which are found in displacive solid-solid phase transformations.     

Chaotic phenomena in the dynamic buckling of an elastic-plastic column under an impact

The present study is concerned with the dynamic anomalous response of an elastic-plastic column struck axially by a massm with an initial velocityv ₀. This simple example is considered in order to clarify the influence of the impact characteristics and the material plastic properties on the dynamic buckling phenomenon and particularly on the final vibration amplitudes of the column when it shakes down to a wholly elastic behaviour. The material is assumed to have a linear strain hardening with a plastic with a plastic reloading allowed. These material properties are the reason a number of elastic-plastic cycles to be realized prior to any wholly elastic stable behaviour, which causes different amounts of energy to be absorbed due to the plastic deformations.The column exhibits two types of behaviour over the range of the impact masses — a quasi-periodic and a chaotic response. The chaotic behaviour is caused by the multiple equilibrium states of the column when any small changes in the loading parameters cause small changes in the plastic strains which result in large changes in the response behaviour. The two types of behaviour are represented by displacement-time and phase-plane diagrams. The sensitivity to the load parameters is illustrated by the calculation of a Lyapunov-like exponent. Poincaré maps are shown for three particular cases.Notation c elastic wave propagation speed - m impact mass - m _c column mass - s step of the spatial discretization - t time - u(x,t) axial displacement - v ₀ initial velocity - w ₀(x) initial imperfections - w(x,t)+w ₀(x) total lateral displacements - x axial axis - z axis along the column thickness - A cross-section areahb - E Young's modulus - E _t hardening modulus (Figure 2) - M(x,t) bending moment - N(x,t) axial force - impact mass ratiom/m _c - (x,z) strain - Lyapunov-like exponent - material density - (x,z) stress     

Magnetohydrodynamic flow in an infinite channel

The magnetohydrodynamic (MHD) flow of an incompressible, viscous, electrically conducting fluid in an infinite channel, under an applied magnetic field has been investigated. The MHD flow between two parallel walls is of considerable practical importance because of the utility of induction flowmeters. The walls of the channel are taken perpendicular to the magnetic field and one of them is insulated, the other is partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 200. All the infinite integrals obtained are transformed to finite integrals which contain modified Bessel functions of the second kind. So, the difficulties associated with the computation of infinite integrals with oscillating integrands which arise for large M have been avoided. It is found that, as M increases, boundary layers are formed near the nonconducting boundaries and in the interface region for both velocity and magnetic fields, and a stagnant region in front of the conducting boundary is developed for the velocity field. Selected graphs are given showing these behaviours.     

Magneto-elastic stress induced by an electric current in an infinite thin plate with an elliptical hole

Two-dimensional magnetic field and stress analyses have been presented for soft ferromagnetic, paramagnetic, and diamagnetic materials of an infinite thin plate with an elliptical hole under steady electric current. The magnetic stress has been analyzed in the Maxwell Stress Model. Except for the approximation of the plane stress state since the plate is the thin plate, any assumption is not made for the stress analysis, though the Maxwell stress components are expressed by nonlinear terms. The boundary condition expressed by Maxwell’s stress is completely satisfied without any linear assumptions on the boundary. Two ways for the boundary condition are stated. The analysis of σ _z in the direction of the plate thickness is also carried out. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived, and the magnitude of the stress intensity factor for the magnetic stress and thermal stress due to the Joule heat caused by the electric current is discussed.     

A non-uniformly travelling shock-wave in an LMS-gas

In this paper an LMS-gas is considered. This is an ideal gas with constant specific heats but with a special entropy distribution, so that generalized Riemann invariants r* and s* exist. A flow with r* and s* constant has the property that the gas moves as a rigid body with a constant acceleration or deceleration. It is shown that a flow composed of two domains, in each of which r* and s* are constant, separated by a plane normal shock-wave is admitted by the Rankine-Hugoniot equations. The path of the shock-wave is calculated exactly and a reasonably realistic situation is suggested to generate this composite flow.     

Phenomenological Modeling of Electrorheological Fluids with an Extended Casson-Model

In this paper we propose a phenomenological theory for electrorheological fluids. In general these are suspensions which undergo dramatic changes in their material properties if they are exposed to an electric field. In the context of continuum mechanics these fluids can be modeled as non-Newtonian fluids. Recalling the governing equations of rational thermodynamics and electrodynamics of moving media (Maxwell-Minkowski-equations), we derive suitable governing equations of electrorheology using essentially two assumptions concerning magnetic quantities. Furthermore we introduce a 3-dimensional nonlinear constitutive equation for the Cauchy stress tensor which is an extension of the model proposed by Ružička (see [14]). Assuming a viscometric flow, we compare the shear stress of our model with other well known models and fit the parameters by using measurements that were obtained in a rotational viscometer. Excellent agreement between model and measurements is achieved. On the basis of these results we propose a 3-dimensional model, the so-called extended Casson -model. This model is investigated further for a channel flow configuration with a homogeneous electric field. We determine analytical solutions for the electric field, the velocity and the volumetric flow rate and illustrate the velocity profiles and the predicted pressure drop. The velocity profiles are flattened compared to parabolic profiles and become more flat if the electric field increases. Received March 21, 2000     

The simulation of oscillatory pipe flow on an analogue computer

Summary A formulation of equations governing oscillatory flow of a fluid in a circular pipe is given. The equations are then cast in a form suitable for analogue computation and a corresponding programming diagram suitable for use with an EAI 680 Analogue Computer is provided. Results have been obtained forNewtonian flow for a range ofReynolds numbers and for combinations of the realReynolds number, and the ratio of elastic to viscous forces in the case of linear viscoelastic fluids.The displacement profiles given in this work do not depend upon specific rheological models and are therefore of universal validity for linear fluids.     

Logarithmic singularity of an elastic composite wedge subjected to out-of-the-plane extensional strain

When an elastic composite wedge is not under a plane strain deformation, an out-of-the-plane extensional strain exists. The singularity analysis for the stresses at the apex of the composite wedge reduces to a system of non-homogeneous linear equations. When the composite wedge consists of two anisotropic elastic materials, it is shown that the stresses have the (ln r) term for all combinations of wedge angles with few exceptions. The same is true when the materials are isotropic except that the (ln r) term may appear in the form of r(ln r) in the displacements only. For these isotropic composite wedges therefore the stresses are bounded, though not continuous, at the apex. However, there are isotropic composite wedges for which the stress singularity is logarithmic. Conditions are given for isotropic composite wedges for which the stresses are (a) uniform, (b) non-uniform but bounded and (c) logarithmic. Unlike the r^−λ singularity, the existence of the (ln r) term does not depend on the complete boundary conditions.     

Waveless Gravity Flow over an Inclined Step

Waveless gravity flows over an inclined step of heightH are studied. The casesH > 0 andH < 0 are considered. It is always assumed that the flow is supercritical at infinity on the right, which ensures the existence of a solution with no downstream wave. For the case of subcritical flow, the relation between the Froude number and the step height is determined that ensures a waveless regime. An approximate analytical relation between the step height and the Froude number is obtained from an analysis of numerical data. This dependence is shown to be nearly identical for steps of any slopes. For the case of supercritical flow, it is established that the problem has a twoparameter set of solutions. For the caseH <0, approximate analytical formulas for the freesurface shape are obtained.     

Diffraction by a grating of cylinders with an arbitrary cross-section

The diffraction of a monochromatic plane electromagnetic wave by a grating consisting of parallel, electrically perfectly conducting cylinders is investigated. Two separate two-dimensional scalar problems are dealt with viz. the case ofE-polarization and the case ofH-polarization. A new Green's function formulation of the problem is employed leading to integral equations for the unknown field distributions on an elementary scatterer of the grating. Numerical results are presented for a grating of cylinders with either square or circular cross-section.     

An internal crack problem for an infinite elastic layer

The elastostatic plane problem of an infinite elastic layer with an internal crack is considered. The elastic layer is subjected to two different loadings, (a) the elastic layer is loaded by a symmetric transverse pair of compressive concentrated forces P/2, (b) it is loaded by a symmetric transverse pair of tensile concentrated forces P/2. The crack is opened by an uniform internal pressure p ₀ along its surface and located halfway between and parallel to the surfaces of the elastic layer. It is assumed that the effect of the gravity force is neglected. Using an appropriate integral transform technique, the mixed boundary value problem is reduced to a singular integral equation. The singular integral equation is solved numerically by making use of an appropriate Gauss–Chebyshev integration formula and the stress-intensity factors and the crack opening displacements are determined according to two different loading cases for various dimensionless quantities.     

On unsteady flow of an elastico-viscous fluid past an infinite plate with variable suction

Approximate solutions of the Navier-Stokes equations are derived through the Laplace transform for two dimensional, incompressible, elastico-viscous flow past a flat porous plate. The flow is assumed to be independent of the distance parallel to the plate. General formulae for the velocity distribution, skin friction and displacement thickness as functions of the given free stream velocity and suction velocity are obtained. The response of skin friction to the impulsive perturbations in the stream and suction velocities is studied. It is found that the order of singularity in the skin friction at t=0 increases due to the elastic property of the fluid in the impulsive case. When the stream is accelerated the skin friction still anticipates the velocity but the time of anticipation is reduced from 1/4 to (1/4) (1—k), where k is the elastic parameter of the fluid. It is found that in general the resistance of the elastico-viscous fluids to an impulsive increase in the stream velocity is greater than the viscous fluids, the elasticoviscous fluids also reach the steady state earlier than the viscous fluids.     

On the dynamics of an elastic body containing a moving crack

It is proved that the energy release rate and the rate of entropy production in the dynamics of an elastic body containing a moving crack are proportional. Moreover, a theorem of the domain of influence type and a uniqueness theorem for solutions to the boundary-initial-value problem of brittle fracture mechanics are proved.     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

Transmission Dynamics of an Influenza Model with Age of Infection and Antiviral Treatment

Age of infection (the time lapsed since infection) is an important factor to consider when modeling the transmission dynamics of influenza under the influence of antiviral treatment and drug-resistance. In this paper, we consider an influenza model which includes an age of infection. The model includes partial differential equations (PDEs) in order to describe the variable infectiousness and effect of antivirals during the infectious period. We derived the formulas for various reproduction numbers (RN) including R_SC{\mathcal R_{SC}} (the controlled RN by one sensitive case), R_TC{\mathcal R_{TC}} (the controlled total RN by one sensitive case), and R_R{\mathcal R_R} (the RN by one resistant case). The model analysis shows that R_SC{\mathcal R_{SC}} and R_R{\mathcal R_R} determine both the global stability of the disease free equilibrium and the existence of the non-trivial equilibria. Local stabilities of the non-trivial equilibria are also discussed. Numerical simulations are conducted to not only confirm or extend the analytic results on qualitative behaviors of the system, but also reveal important quantitative properties of the disease dynamics influenced by antiviral treatment. These results are then used to assess the effectiveness of treatment programs in terms of both the RNs and the epidemic size. Our findings illustrate possibility that a higher level of antiviral use may lead to an increase of the epidemic size, despite the fact that there is a fitness cost of the drug-resistant strains. This suggests that programs for antiviral use should be designed carefully to avoid the adverse effect.     

Calculations and experiments on an optimal adjustment of function generators

Summary Formulae are derived for the approximation of a given arbitrary functionr(x) defined forx ₀xx _n by a broken line consisting ofn linear parts. Numerical data, obtained for any arbitrary value ofn, are important for analogue computer applications involving the use of function generators.Experiments have shown that the determination of an optimal adjustment of a function generator forr(x) is simple and accurate if performed with the aid of an analogue computer.Shell Research N.V.     

Arrest characteristics of an impacted crack dissipating plastic energy

When a structural member is accidentally struck, an initial defect may grow and then arrest. An estimate of its size increase is made by considering the geometry of a centrally cracked panel subjected to a step function load in time. Two load amplitudes differing by a factor of five are considered. Under impact, the crack accelerates and then decelerates prior to arrest. The dynamic characteristics depend on the interaction of the elastic-plastic stress waves intervening with the physical boundaries. This effect is assessed quantitatively by computing for the energy stored in a unit volume of material and by incorporating sliding nodes in the finite element method. The energy dissipated by plastic deformation must be accounted for as it is no longer available for creating new macrocrack surface. Obtained are the near tip normal stresses that are found to change from compression to tension. Their magnitude is considerably larger than the corresponding static values. Increase in crack length changes from 0.87% to 20.6% when the magnitude of the impact load is raised five times. The rate of change of the strain energy density factor ΔS with crack growth Δa is found to be governed by the condition ΔS/Δa = const. during loading while dynamic relaxation corresponded to a nonlinear behavior. The physical implication of this remains to be clarified in view of the fact that plasticity theory may not adequately explain the near tip crack behavior.