A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.

Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontalwells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid wellcollapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group, London.

A new approach to the analytical and numerical form-finding of tensegrity structures

We develop a new formulation for the form-finding of tensegrity structures in which the primary variables are the Cartesian components of element lengths. Both an analytical and a numerical implementation of the formulation are described; each require a description of the connectivity of the tensegrity, with the iterative numerical method also requiring a random starting vector of member force densities. The analytical and numerical form-finding of tensegrity structures is demonstrated through six examples, and the results obtained are compared and contrasted with those available in the literature to verify the accuracy and viability of the suggested methods. © 2013 Elsevier Ltd. All rights reserved.

Experimental behaviour of barrel vaulted structure subjected to rocking of the abutments

In this paper, the transverse rocking mechanism of a barrel vaulted structure subjected to horizontal cyclic loads is analysed by means of experimental tests on full scale model and by means of non-linear FE analyses. The study is part of an ongoing experimental and theoretical research program, developed by the University of Brescia, concerning the seismic behaviour of ancient masonry buildings. The scope of the paper is to provide some evidence of the rocking mechanism experienced by barrel vaulted structures under horizontal loading. The understanding of the behaviour of these structural systems is necessary for their seismic vulnerability assessment, as well as for the correct design of possible strengthening techniques. A numeric FE model was validated through comparison with the experimental results and it was used to verify the efficiency of two common strengthening solutions: the technique of the overlaying reinforced concrete slab and the technique of the thin spandrel walls. Experimental and numeric results will be discussed in the paper.

The interaction of elasticity and rocking in flexible structures allowed to uplift

Numerous structures uplift under the influence of strong ground motion. Although many researchers have investigated the effects of base uplift on very stiff (ideally rigid) structures, the rocking response of flexible structures has received less attention. Related practical analysis methods treat these structures with simplified 'equivalent' oscillators without directly addressing the interaction between elasticity and rocking. This paper addresses the fundamental dynamics of flexible rocking structures. The nonlinear equations of motion, derived using a Lagrangian formulation for large rotations, are presented for an idealized structural model. Particular attention is devoted to the transition between successive phases; a physically consistent classical impact framework is utilized alongside an energy approach. The fundamental dynamic properties of the flexible rocking system are compared with those of similar linear elastic oscillators and rigid rocking structures, revealing the distinct characteristics of flexible rocking structures. In particular, parametric analysis is performed to quantify the effect of elasticity on uplift, overturning instability, and harmonic response, from which an uplifted resonance emerges. The contribution of stability and strength to the collapse of flexible rocking structures is discussed. © 2012 John Wiley & Sons, Ltd.

A theoretical approach to the passive control of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this.We study how this spiral mode of vortex breakdown might be suppressed or promoted. We use a Lagrangian approach to identify regions of the flow which are sensitive to small open-loop steady and unsteady (harmonic) forces. We find these regions to be upstream of the vortex breakdown bubble. We investigate passive control using a small axisymmetric control ring. In this case, the steady and unsteady control forces are caused by the drag force on the control ring. We find a narrow region upstream of the bubble where the control ring will stabilise the flow and we verify this using numerical simulations. © 2012 IEEE.

Transverse cracking of a simple lap joint under axial load

The failure mode of axially loaded simple, single lap joints formed between thin adherends which are flexible in bending is conventionally described as one of axial peeling. We have observed - using high-speed photography - that it is also possible for failure to be preceded by the separation front, or crack, moving in a transverse direction, i.e. perpendicular to the direction of the axial load. A simple energy balance analysis suggests that the critical load for transverse failure is the same as that for axial separation for both flexible lap joints, where the bulk of the stored elastic energy lies in the adhesive, and structural lap joints in which the energy stored in the adherends dominates. The initiation of the failure is dependent on a local increases in either stress or strain energy to some critical values. In the case of a flexible joint, this will occur within the adhesive layer and the critical site will be close to one of the corners of the joint overlap from which the separation front can proceed either axially or transversely. These conclusions are supported by a finite element analysis of a joint formed between adherends of finite width by a low modulus adhesive. © 2012 Taylor & Francis.