Two-dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first-order and second-order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth-order RungeKutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two-dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side-by-side. Results of these simulations were extensively compared with the previous numerical data. Copyright (C) 2011 John Wiley & Sons, Ltd.

VAM applied to dimensional reduction of non-linear hyperelastic plates

This work aims at dimensional reduction of non-linear isotropic hyperelastic plates in an asymptotically accurate manner. The problem is both geometrically and materially non-linear. The geometric non-linearity is handled by allowing for finite deformations and generalized warping while the material non-linearity is incorporated through hyperelastic material model. The development, based on the Variational Asymptotic Method (VAM) with moderate strains and very small thickness to shortest wavelength of the deformation along the plate reference surface as small parameters, begins with three-dimensional (3-D) non-linear elasticity and mathematically splits the analysis into a one-dimensional (1-D) through-the-thickness analysis and a two-dimensional (2-D) plate analysis. Major contributions of this paper are derivation of closed-form analytical expressions for warping functions and stiffness coefficients and a set of recovery relations to express approximately the 3-D displacement, strain and stress fields. Consistent with the 2-D non-linear constitutive laws, 2-D plate theory and corresponding finite element program have been developed. Validation of present theory is carried out with a standard test case and the results match well. Distributions of 3-D results are provided for another test case. (c) 2012 Elsevier Ltd. All rights reserved.

Seismic stability of a long unsupported circular tunnel

The stability of a long unsupported circular tunnel (opening) in a cohesive frictional soil has been assessed with the inclusion of pseudo-static horizontal earthquake body forces. The analysis has been performed under plane strain conditions by using upper bound finite element limit analysis in combination with a linear optimization procedure. The results have been presented in the form of a non-dimensional stability number (gamma H-max/c); where H = tunnel cover, c refers to soil cohesion and gamma(max) is the maximum unit weight of soil mass which the tunnel can support without collapse. The results have been obtained for various values of H/D (D = diameter of the tunnel), internal friction angle (phi) of soil, and the horizontal earthquake acceleration coefficient (alpha(h)). The computations reveal that the values of the stability numbers (i) decrease quite significantly with an increase in alpha(h), and (ii) become continuously higher for greater values of H/D and phi. As expected, the failure zones around the periphery of the tunnel becomes always asymmetrical with an inclusion of horizontal seismic body forces. (c) 2012 Elsevier Ltd. All rights reserved.

A note on the inverse mode shape problem for bars, beams, and plates

Motivated by the idea of designing a structure for a desired mode shape, intended towards applications such as resonant sensors, actuators and vibration confinement, we present the inverse mode shape problem for bars, beams and plates in this work. The objective is to determine the cross-sectional profile of these structures, given a mode shape, boundary condition and the mass. The contribution of this article is twofold: (i) A numerical method to solve this problem when a valid mode shape is provided in the finite element framework for both linear and nonlinear versions of the problem. (ii) An analytical result to prove the uniqueness and existence of the solution in the case of bars. This article also highlights a very important question of the validity of a mode shape for any structure of given boundary conditions.

Minimum power consumption for drag reduction on a circular cylinder by tangential surface motion

We investigate the effect of a prescribed tangential velocity on the drag force on a circular cylinder in a spanwise uniform cross flow. Using a combination of theoretical and numerical techniques we make an attempt at determining the optimal tangential velocity profiles which will reduce the drag force acting on the cylindrical body while minimizing the net power consumption characterized through a non-dimensional power loss coefficient (C-PL). A striking conclusion of our analysis is that the tangential velocity associated with the potential flow, which completely suppresses the drag force, is not optimal for both small and large, but finite Reynolds number. When inertial effects are negligible (R e << 1), theoretical analysis based on two-dimensional Oseen equations gives us the optimal tangential velocity profile which leads to energetically efficient drag reduction. Furthermore, in the limit of zero Reynolds number (Re -> 0), minimum power loss is achieved for a tangential velocity profile corresponding to a shear-free perfect slip boundary. At finite Re, results from numerical simulations indicate that perfect slip is not optimum and a further reduction in drag can be achieved for reduced power consumption. A gradual increase in the strength of a tangential velocity which involves only the first reflectionally symmetric mode leads to a monotonic reduction in drag and eventual thrust production. Simulations reveal the existence of an optimal strength for which the power consumption attains a minima. At a Reynolds number of 100, minimum value of the power loss coefficient (C-PL = 0.37) is obtained when the maximum in tangential surface velocity is about one and a half times the free stream uniform velocity corresponding to a percentage drag reduction of approximately 77 %; C-PL = 0.42 and 0.50 for perfect slip and potential flow cases, respectively. Our results suggest that potential flow tangential velocity enables energetically efficient propulsion at all Reynolds numbers but optimal drag reduction only for Re -> infinity. The two-dimensional strategy of reducing drag while minimizing net power consumption is shown to be effective in three dimensions via numerical simulation of flow past an infinite circular cylinder at a Reynolds number of 300. Finally a strategy of reducing drag, suitable for practical implementation and amenable to experimental testing, through piecewise constant tangential velocities distributed along the cylinder periphery is proposed and analysed.

A constant factor approximation algorithm for boxicity of circular arc graphs

Boxicity of a graph G(V, E) is the minimum integer k such that G can be represented as the intersection graph of k-dimensional axis parallel boxes in Rk. Equivalently, it is the minimum number of interval graphs on the vertex set V such that the intersection of their edge sets is E. It is known that boxicity cannot be approximated even for graph classes like bipartite, co-bipartite and split graphs below O(n0.5-ε)-factor, for any ε > 0 in polynomial time unless NP = ZPP. Till date, there is no well known graph class of unbounded boxicity for which even an nε-factor approximation algorithm for computing boxicity is known, for any ε < 1. In this paper, we study the boxicity problem on Circular Arc graphs - intersection graphs of arcs of a circle. We give a (2+ 1/k)-factor polynomial time approximation algorithm for computing the boxicity of any circular arc graph along with a corresponding box representation, where k ≥ 1 is its boxicity. For Normal Circular Arc(NCA) graphs, with an NCA model given, this can be improved to an additive 2-factor approximation algorithm. The time complexity of the algorithms to approximately compute the boxicity is O(mn+n2) in both these cases and in O(mn+kn2) which is at most O(n3) time we also get their corresponding box representations, where n is the number of vertices of the graph and m is its number of edges. The additive 2-factor algorithm directly works for any Proper Circular Arc graph, since computing an NCA model for it can be done in polynomial time.

Asymptotic wavenumber expansions of a strongly orthotropic fluid-filled circular cylindrical shell

In this paper, a suitable nondimensional `orthotropy parameter' is defined and asymptotic expansions are found for the wavenumbers in in vacuo and fluid-filled orthotropic circular cylindrical shells modeled by the Donnell-Mushtari theory. Here, the elastic moduli in the two directions are greatly different; the particular case of E-x >> E-theta is studied in detail, i.e., the elastic modulus in the longitudinal direction is much larger than the elastic modulus in the circumferential direction. These results are compared with the corresponding results for a `slightly orthotropic' shell (E-x approximate to E-theta) and an isotropic shell. The novelty of this presentation lies in obtaining closed-form expansions for the in vacuo and coupled wavenumbers in an orthotropic shell using perturbation methods aiding in a better physical understanding of the problem.

Recent advances in muffler acoustics

Exhaust noise in engines has always been a major source of automotive noise. Challenges for muffler design have been constraints on size, back pressure, and, of course, the cost. Designing for sufficient insertion loss at the engine firing frequency and the first few harmonics has been the biggest challenge. Most advances in the design of efficient mufflers have resulted from linear plane wave theory, making use of the transfer matrix method. This review paper deals with evaluating approximate source characteristics required for prediction of the unmuffled intake and exhaust noise, making use of the electroacoustical analogies. In the last few years, significant advances have been made in the analysis of variable area perforated ducts, transverse plane wave analysis of short elliptical as well as circular chambers, double-tuned expansion chambers and concentric tube resonators, catalytic converters, diesel particulate filters, air cleaners, etc. The development of long strand fibrous materials that can be used in hot exhaust systems without binders has led to the use of combination mufflers in exhaust systems. Breakthroughs have been achieved in the prediction and control of breakout noise from the elliptical and circular muffler shell as well as the end plates of typical mufflers. Diesel particulate filters and inlet air cleaners have also been modeled acoustically. Some of these recent advances are the subject of this review paper.

Asymptotic expansions for the structural wavenumbers of isotropic and orthotropic fluid-filled circular cylindrical shells in the intermediate frequency range

We consider wavenumbers in in vacuo and fluid-filled isotropic and orthotropic shells. Using the Donnell-Mushtari (DM) theory we find compact and elegant asymptotic expansions for the wavenumbers in the intermediate frequency range, i.e., around the ring frequency. This frequency range corresponds to the frequencies where there is a rapid change in the values of bending wavenumbers and is found to exist in isotropic and orthotropic shells (in vacua and fluid-filled) for low circumferential orders n only. The same is first identified using the n=0 mode of an orthotropic shell. Following this, using the expression for the intermediate frequency, asymptotic expansions are found for other cases. Here, in order to get compact expansions we consider slight orthotropy (epsilon << 1) and light fluid loading (mu << 1). Thus, the orthotropy parameter epsilon and the fluid loading parameter mu are used as asymptotic parameters along with the non-dimensional thickness parameter beta. The methodology can be extended to any order of epsilon, only the expansions become unwieldy. The expansions are matched with the numerical solutions of the corresponding dispersion relation. The match is found to be good.

Stability of a long unsupported circular tunnel in soils with seismic forces

By using the lower-bound finite element limit analysis, the stability of a long unsupported circular tunnel has been examined with an inclusion of seismic body forces. The numerical results have been presented in terms of a non-dimensional stability number (gamma H/c) which is plotted as a function of horizontal seismic earth pressure coefficient (k (h)) for different combinations of H/D and I center dot; where (1) H is the depth of the crest of the tunnel from ground surface, (2) D is the diameter of the tunnel, (3) k (h) is the earthquake acceleration coefficient and (4) gamma, c and I center dot define unit weight, cohesion and internal friction angle of soil mass, respectively. The stability numbers have been found to decrease continuously with an increase in k (h). With an inclusion of k (h), the plastic zone around the periphery of the tunnel becomes asymmetric. As compared to the results reported in the literature, the present analysis provides a little lower estimate of the stability numbers. The numerical results obtained would be useful for examining the stability of unsupported tunnel under seismic forces.

Effect of steady rotation on low Reynolds number vortex shedding behind a circular cylinder

In this paper control of oblique vortex shedding in the wake behind a straight circular cylinder is explored experimentally and computationally. Towards this, steady rotation of the cylinder about its axis is used as a control device. Some limited studies are also performed with a stepped circular cylinder, where at the step the flow is inevitably three-dimensional irrespective of the rotation rate. When there is no rotation, the vortex shedding pattern is three dimensional as described in many previous studies. With a non-zero rotation rate, it is demonstrated experimentally as well as numerically that the shedding pattern becomes more and more two-dimensional. At sufficiently high rotation rates, the vortex shedding is completely suppressed.

On Uniform Approximate Solutions in Bending of Symmetric Laminated Plates

A layer-wise theory with the analysis of face ply independent of lamination is used in the bending of symmetric laminates with anisotropic plies. More realistic and practical edge conditions as in Kirchhoff's theory are considered. An iterative procedure based on point-wise equilibrium equations is adapted. The necessity of a solution of an auxiliary problem in the interior plies is explained and used in the generation of proper sequence of two dimensional problems. Displacements are expanded in terms of polynomials in thickness coordinate such that continuity of transverse stresses across interfaces is assured. Solution of a fourth order system of a supplementary problem in the face ply is necessary to ensure the continuity of in-plane displacements across interfaces and to rectify inadequacies of these polynomial expansions in the interior distribution of approximate solutions. Vertical deflection does not play any role in obtaining all six stress components and two in-plane displacements. In overcoming lacuna in Kirchhoff's theory, widely used first order shear deformation theory and other sixth and higher order theories based on energy principles at laminate level in smeared laminate theories and at ply level in layer-wise theories are not useful in the generation of a proper sequence of 2-D problems converging to 3-D problems. Relevance of present analysis is demonstrated through solutions in a simple text book problem of simply supported square plate under doubly sinusoidal load.

Stability of long unsupported twin circular tunnels in soils

The stability of two long unsupported circular parallel tunnels aligned horizontally in fully cohesive and cohesive-frictional soils has been determined. An upper bound limit analysis in combination with finite elements and linear programming is employed to perform the analysis. For different clear spacing (S) between the tunnels, the stability of tunnels is expressed in terms of a non-dimensional stability number (gamma H-max/c); where H is tunnel cover, c refers to soil cohesion, and gamma(max) is maximum unit weight of soil mass which the tunnels can bear without any collapse. The variation of the stability number with tunnels' spacing has been established for different combinations of H/D, m and phi; where D refers to diameter of each tunnel, phi is the internal friction angle of soil and m accounts for the rate at which the cohesion increases linearly with depth. The stability number reduces continuously with a decrease in the spacing between the tunnels. The optimum spacing (S-opt) between the two tunnels required to eliminate the interference effect increases with (i) an increase in H/D and (ii) a decrease in the values of both m and phi. The value of S-opt lies approximately in a range of 1.5D-3.5D with H/D = 1 and 7D-12D with H/D = 7. The results from the analysis compare reasonably well with the different solutions reported in literature. (C) 2013 Elsevier Ltd. All rights reserved.

Detonation-driven-shock wave interactions with perforated plates

The study of detonations and their interactions is vital for the understanding of the high-speed flow physics involved and the ultimate goal of controlling their detrimental effects. However, producing safe and repeatable detonations within the laboratory can be quite challenging, leading to the use of computational studies which ultimately require experimental data for their validation. The objective of this study is to examine the induced flow field from the interaction of a shock front and accompanying products of combustion, produced from the detonation taking place within a non-electrical tube lined with explosive material, with porous plates with varying porosities, 0.7-9.7%. State of the art high-speed schlieren photography alongside high-resolution pressure measurements is used to visualise the induced flow field and examine the attenuation effects which occur at different porosities. The detonation tube is placed at different distances from the plates' surface, 0-30 mm, and the pressure at the rear of the plate is recorded and compared. The results indicate that depending on the level of porosity and the Mach number of the precursor shock front secondary reflected and transmitted shock waves are formed through the coalescence of compression waves. With reduced porosity, the plates act almost as a solid surface, therefore the shock propagates faster along its surface.

A parametric study on the buckling of functionally graded material plates with internal discontinuities using the partition of unity method

In this paper, the effect of local defects, viz., cracks and cutouts on the buckling behaviour of functionally graded material plates subjected to mechanical and thermal load is numerically studied. The internal discontinuities, viz., cracks and cutouts are represented independent of the mesh within the framework of the extended finite element method and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction and the effective properties are estimated using the Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation theory. The influence of various parameters, viz., the crack length and its location, the cutout radius and its position, the plate aspect ratio and the plate thickness on the critical buckling load is studied. The effect of various boundary conditions is also studied. The numerical results obtained reveal that the critical buckling load decreases with increase in the crack length, the cutout radius and the material gradient index. This is attributed to the degradation in the stiffness either due to the presence of local defects or due to the change in the material composition. (C) 2013 Elsevier Masson SAS. All rights reserved.