Residual stress in polyethylene pipes

This paper deals with residual stress in polyethylene (PE) pipes as an important factor that influences their lifetime. Residual hoop stress distribution in a PE pipe was determined using a methodology previously carried out by the authors on polypropylene pipes. Axial residual stress magnitude was determined by comparison of experimental data and finite element modelling. Based on the obtained results, a new simplified methodology for determining the residual hoop stress is proposed. The method needs only one circular specimen made of pipe, but, unlike older methods, provides a more precise estimation of residual hoop stress distribution taking into account its exponential shape. Some older results from literature are recalculated using this method and residual hoop stress of various PE pipe grades and dimensions are then compared. To provide an idea of residual hoop stress influence on lifetime of a pipe, a lifetime estimation is carried out for the examined pipe.     

The intersection marker method for 3D interface tracking of deformable surfaces in finite volumes

Currently, the majority of computational fluid dynamics (CFD) codes use the finite volume method to spatially discretise the computational domain, sometimes as an array of cubic control volumes. The Finite volume method works well with single‐phase flow simulations, but two‐phase flow simulations are more challenging because of the need to track the surface interface traversing and deforming within the 3D grid. Surface area and volume fraction details of each interface cell must be accurately accounted for, in order to calculate for the momentum exchange and rates of heat and mass transfer across the interface. To attain a higher accuracy in two‐phase flow CFD calculations, the intersection marker (ISM) method is developed. The ISM method is a hybrid Lagrangian–Eulerian front‐tracking algorithm that can model an arbitrary 3D surface within an array of cubic control volumes. The ISM method has a cell‐by‐cell remeshing capability that is volume conservative and is suitable for the tracking of complex interface deformation in transient two‐phase CFD simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Palladium‐Catalyzed Direct Cross‐Coupling of Carboranyllithium with (Hetero)Aryl Halides

A palladium‐catalyzed direct C‐arylation reaction of readily available cage carboranyllithium reagents with aryl halides has been developed for the first time. This method is applicable to a wide range of aryl halide substrates including aryl iodides, aryl bromides, and heteroaromatic halides.     

Permeability of microscale fibrous porous media using the lattice Boltzmann method

The permeabilities of microscale fibrous porous media were calculated using the multiple-relaxation-time (MRT) lattice Boltzmann method (LBM). Two models of the microscale fibrous porous media were constructed based on overlapping fibers (simple cubic, body-centered cubic). Arranging the fibers in skew positions yielded two additional models comprising non-overlapping fibers (skewed simple cubic, skewed body-centered cubic). As the fiber diameter increased, the fibers acted as granular inclusions. The effects of the overlapping fibers on the media permeability were investigated. The overlapping fibers yielded permeability values that were a factor of 2.5 larger than those obtained from non-overlapping fibers, but the effects of the fiber arrangement were negligible. Two correlations were obtained for the overlapping and non-overlapping fiber models, respectively. The effects of the rarefaction and slip flow are also discussed. As the Knudsen number increased, the dimensionless permeability increased; however, the increase differed depending on the fiber arrangement. In the slip flow regime, the fiber arrangement inside the porous media became an important factor.     

Dust modelling using a combined CFD and discrete element formulation

The production and dispersal of airborne dust is an important issue in both environmental and industrial contexts. Dust pollution is a major environmental concern, and long exposure in occupational settings has been linked with numerous respiratory health issues. Industrial dust pollution can also present a significant explosion hazard, even in facilities with dust extraction systems. Computational models for dust generation and dispersal have, however, generally been formulated for specific geophysical applications and restricted to static, two‐dimensional, approaches. Here, we present a method for simulating dust production from a dynamic granular bed by using a three‐dimensional coupled discrete element method and Navier–Stokes computational model. Dust production is based on an energy formulation in which micro‐scale dust particles are assumed to overcome cohesion to macro‐scale grains. This model is used over the entire range of energies present within the system, from macro‐scale collisions to aerodynamic entrainment and bombardment of micro‐scale particles. The dust concentration is modelled as a scalar density field, which is advected and diffused through turbulence in the gas flow field. The model is tested against empirical predictions for two test cases, a slug of granular material dropped from a set height and air flow over a granular stockpile. Both give good agreement to the empirical relations, showing that the model can accurately predict the production and subsequent dispersal of dust from a dynamic granular bed. Copyright © 2012 John Wiley & Sons, Ltd.     

FREE VIBRATION ANALYSIS OF LATTICE SANDWICH BEAMS UNDER SEVERAL TYPICAL BOUNDARY CONDITIONS

Free vibration problems of lattice sandwich beams under several typical boundary conditions are investigated in the present paper. The lattice sandwich beam is transformed to an equivalent homogeneous three-layered sandwich beam. Unlike the traditional analytical model in which the rotation angles of the face sheets and the core are assumed the same, different rotation angles are considered in this paper to characterize the real response of sandwich beams. The analytical solutions of the natural frequencies for several typical boundary conditions are obtained. The effects of material properties and geometric parameters on the natural frequencies are also investigated.     

An adaptive enrichment algorithm for advection‐dominated problems

We are interested in developing a numerical framework well suited for advection–diffusion problems when the advection part is dominant. In that case, given Dirichlet type boundary condition, it is well known that a boundary layer develops. To resolve correctly this layer, standard methods consist in increasing the mesh resolution and possibly increasing the formal accuracy of the numerical method. In this paper, we follow another path: we do not seek to increase the formal accuracy of the scheme but, by a careful choice of finite element, to lower the mesh resolution in the layer. Indeed the finite element representation we choose is locally the sum of a standard one plus an enrichment. This paper proposes such a method and with several numerical examples, we show the potential of this approach. In particular, we show that the method is not very sensitive to the choice of the enrichment and develop an adaptive algorithm to automatically choose the enrichment functions.Copyright © 2012 John Wiley & Sons, Ltd.     

Relationship between Euler buckling and unstable equilibria of buckled beams

Classic snap-through of curved beams, plates, and shells has long been an object of attention in structural engineering. Euler buckling under axial loading is perhaps an even more entrenched part of the canon of engineering education and practice. In this paper we introduce a relationship between the two phenomena, that to our knowledge has not been directly addressed before. The relationship shows that Euler buckling configurations are connected by the force–displacement curve under transverse loading. The results are used to develop a very simple metric to estimate the number of unstable static equilibria of a buckled structure based only on its geometry with no need for static or dynamic solvers. The study is focused on beams as this allows for an unambiguous discussion of the idea on the simplest possible structure.     

A hybrid level set method for modelling detonation and combustion problems in complex geometries

We present an accurate and fast wave tracking method that uses parametric representations of tracked fronts, combined with modifications of level set methods that use narrow bands. Our strategy generates accurate computations of the front curvature and other geometric properties of the front. We introduce data structures that can store discrete representations of the location of the moving fronts and boundaries, as well as the corresponding level set fields, that are designed to reduce computational overhead and memory storage. We present an algorithm we call stack sweeping to efficiently sort and store data that is used to represent orientable fronts. Our implementation features two reciprocal procedures, a forward ‘front parameterization’ that constructs a parameterization of a front given a level set field and a backward ‘field construction’ that constructs an approximation of the signed normal distance to the front, given a parameterized representation of the front. These reciprocal procedures are used to achieve and maintain high spatial accuracy. Close to the front, precise computation of the normal distance is carried out by requiring that displacement vectors from grid points to the front be along a normal direction. For front curves in two dimensions, a cubic interpolation scheme is used, and G ¹ surface parameterization based on triangular patches is used for the three-dimensional implementation to compute the distances from grid points near the front. To demonstrate this new, high accuracy method we present validations and show examples of combustion-like applications that include detonation shock dynamics, material interface motions in a compressible multi-material simulation and the Stephan problem associated with dendrite solidification.     

Distorted wave response of ultrasonic annular stator incorporating non-uniform geometry

The traveling wave ultrasonic stator is normally fabricated with teeth. The tooth geometry improves the driving speed, but it creates natural frequency splitting and mode contamination, especially a distorted traveling wave. A dynamic model of a stepped-plate periodic stator is developed to examine the distortion. The stator is treated as an annular supported by a thin mid plate, and the support stiffness is formulated by using equivalent energy principle. The effects of the tooth and mid plate on the natural frequency and vibration mode are examined by using the perturbation method. The rules governing the frequency splitting, frequency perturbation as well as mode contamination are also identified. The traveling wave response and elliptical trace on stator surface are obtained by using the mode superposition method and they are proved to be distorted due to the tooth geometry. The response at the repeated doublets becomes coupled forward and backward traveling waves, but that at the split doublets becomes coupled forward traveling, standing and backward traveling waves. The results indicate that the tooth mass instead of the stiffness decreases the vibration amplitude and driving speed of the dominant wave, but their effects are different at the repeated and split doublets. Inspection of the model implies that the distortion can be suppressed by using a suitable combination of the wavenumber, tooth count, tooth height and occupying fraction. Numerical calculations are carried out to demonstrate the tooth geometry effect on the transient waveform, driving speed and elliptical trace. The optimization of the tooth geometry that can help achieve a purer traveling wave is discussed.