A Fire Fighter's Problem

Suppose that a circular fire spreads in the plane at unit speed. A single fire fighter can build a barrier at speed $v>1$. How large must $v$ be to ensure that the fire can be contained, and how should the fire fighter proceed? We contribute two results. First, we analyze the natural curve \mbox{FF}_v that develops when the fighter keeps building, at speed $v$, a barrier along the boundary of the expanding fire. We prove that the behavior of this spiralling curve is governed by a complex function $(e^{w Z} - s \, Z)^{-1}$, where $w$ and $s$ are real functions of $v$. For $v>v_c=2.6144 \ldots$ all zeroes are complex conjugate pairs. If $\phi$ denotes the complex argument of the conjugate pair nearest to the origin then, by residue calculus, the fire fighter needs $\Theta( 1/\phi)$ rounds before the fire is contained. As $v$ decreases towards $v_c$ these two zeroes merge into a real one, so that argument $\phi$ goes to~0. Thus, curve \mbox{FF}_v does not contain the fire if the fighter moves at speed $v=v_c$. (That speed $v>v_c$ is sufficient for containing the fire has been proposed before by Bressan et al. [7], who constructed a sequence of logarithmic spiral segments that stay strictly away from the fire.) Second, we show that any curve that visits the four coordinate half-axes in cyclic order, and in inreasing distances from the origin, needs speed $v>1.618\ldots$, the golden ratio, in order to contain the fire. Keywords: Motion Planning, Dynamic Environments, Spiralling strategies, Lower and upper boundsComment: A preliminary version of the paper was presented at SoCG 201

Numerical modelling of behaviour of reinforced concrete columns in fire and comparison with Eurocode 2

The paper describes a two-step finite element formulation for the thermo-mechanical non-linear analysis of the behaviour of the reinforced concrete columns in fire. In the first step, the distributions of the temperature over the cross-section during fire are determined. In the next step, the mechanical analysis is made in which these distributions are used as the temperature loads. The analysis employs our new strain-based planar geometrically exact and materially non-linear beam finite elements to model the column. The results are compared with the measurements of the full-scale test on columns in fire and with the results of the European building code EC 2. The resistance times of the present method and the test were close. It is also noted that the building code EC 2 might be non-conservative in the estimation of the resistance time. (c) 2005 Elsevier Ltd. All rights reserved

An embedded FE model for modelling reinforced concrete slabs in fire

This is the post-print version of the final paper published in Engineering Structures. The published article is available from the link below. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. Copyright @ 2008 Elsevier B.V.It is evident from a series of tests on simply supported reinforced concrete slabs that the failure of the slabs at large deflections is due to the formation of individual large cracks. This failure mode was also observed in the Cardington full-scale fire tests. Previous research indicates that the global behaviour of concrete slabs subject to large deflections can be well predicted by the smeared cracking model; however, the model cannot quantitatively predict the openings of individual cracks within the slabs at large deflections. For the discrete approach it is usually assumed that the cracks are formed along element edges, therefore continuous re-meshing is required during the analysis. Consequently, the results are mesh-dependent and the computing cost is high. In recent years, mesh independent finite element procedures, such as embedded (EFEM) and extended (XFEM) approaches, were widely used for modelling of the crack initiation and growth in structural members. However, most of the meshless models developed are either based on in-plane loading conditions or confined to thin shells with assumed full-depth cracks, which form apparent displacement jumps within an element. For a reinforced concrete slab, an out-of-plane load causes coupled stretching and bending of the slab, cracks are usually initiated at discrete positions and then propagated, until at last some individual full-depth cracks are formed. Pure stretching or assumed full-depth cracking is inadequate for modelling this kind of failure. Therefore, in this research, a non-linear layered procedure with embedded weak discontinuity is developed to quantitatively model the progressive tensile failure of reinforced concrete slabs subjected to large deflections. The current model inherits the advantage of the smeared approach, and at the same time, introduces the opening width of crack explicitly by taking the advantage of the better description of the kinematic characteristics of the EFEM approach. A series of validations have been conducted against test data at both ambient and elevated temperatures, and the research shows that the model developed in this paper is not sensitive to the FE mesh size and the aspect ratio of the slab. The results predicted by the model developed agreed well with the test data in terms of deflection and crack open width, also agreeing well with those modelled by the smeared model. Hence, this new approach provides a numerical method to predict the load capacity as well as identifying the occurrence and severity of crack failure in reinforced concrete slabs subjected to extreme loading conditions, such as fire

An extended finite element model for modelling localised fracture of reinforced concrete beams in fire

Open Access funded by Engineering and Physical Sciences Research Council under a Creative Commons license.A robust finite element procedure for modelling the localised fracture of reinforced concrete beams at elevated temperatures is developed. In this model a reinforced concrete beam is represented as an assembly of 4-node quadrilateral plain concrete, 3-node main reinforcing steel bar, and 2-node bond-link elements. The concrete element is subdivided into layers for considering the temperature distribution over the cross-section of a beam. An extended finite element method (XFEM) has been incorporated into the concrete elements in order to capture the localised cracks within the concrete. The model has been validated against previous fire test results on the concrete beams.The Engineering and Physical Sciences Research Council of Great Britain under Grant No. EP/I031553/1

Prediction of smoke filling in large volumes by means of data assimilation-based numerical simulations

The concept of numerical simulations for real-time Numerical Fire Forecasting is illustrated for the case of natural smoke filling of a large-scale atrium in case of fire. The numerical simulations are performed within the Inverse Zone Modelling framework. The technique consists of assimilating collected data for a certain parameter, in casu the smoke layer height, into the zone model in order to estimate an unknown of the problem ('model invariant'), mainly the fire heat release rate. A forecast in terms of evolution of smoke level and temperature can then be produced. Because zone model calculations are very fast, positive lead times of several minutes are obtained. The developed model produces reliable forecasts for the cases considered. Equally important, the robustness of the technique is illustrated: the sensitivity of the results to the 'initial guess' of the model invariants is small (i.e. the method converges easily); one model invariant is sufficient to obtain reliable predictions for smoke layer height evolution; the data assimilation window length does not affect the results significantly. The method automatically provides a different value for the plume entrainment constant, depending on the position of the fire (in the middle of the atrium or in a corner)

Discontinuity induced bifurcations of non-hyperbolic cycles in nonsmooth systems

We analyse three codimension-two bifurcations occurring in nonsmooth systems, when a non-hyperbolic cycle (fold, flip, and Neimark-Sacker cases, both in continuous- and discrete-time) interacts with one of the discontinuity boundaries characterising the system's dynamics. Rather than aiming at a complete unfolding of the three cases, which would require specific assumptions on both the class of nonsmooth system and the geometry of the involved boundary, we concentrate on the geometric features that are common to all scenarios. We show that, at a generic intersection between the smooth and discontinuity induced bifurcation curves, a third curve generically emanates tangentially to the former. This is the discontinuity induced bifurcation curve of the secondary invariant set (the other cycle, the double-period cycle, or the torus, respectively) involved in the smooth bifurcation. The result can be explained intuitively, but its validity is proven here rigorously under very general conditions. Three examples from different fields of science and engineering are also reported

Non-linear fire-resistance analysis of reinforced concrete beams

The non-linear structural analysis of reinforced concrete beams in fire consists of three separate steps: (i) The estimation of the rise of surrounding air temperature due to fire; (ii) the determination of the distribution of the temperature within the beam during fire; (iii) the evaluation of the mechanical response due to simultaneous time-dependent thermal and mechanical loads. Steps (ii) and (iii) are dealt with in the present paper. We present a two-step computational procedure where a 2D transient thermal analysis over the cross-sections of beams are made first, followed by mechanical analysis of the structure. Fundamental to the accuracy of the mechanical analysis is a new planar beam finite element. The effects of plasticity in concrete, and plasticity and viscous creep in steel are taken into consideration. The properties of concrete and steel along with the values of their thermal and mechanical parameters are taken according to the European standard ENV 1992-1-2 (1995). The comparison of our numerical and full-scale experimental results shows that the proposed mechanical and 2D thermal computational procedure is capable to describe the actual response of reinforced concrete beam structures to fire

Stochastic mean field formulation of the dynamics of diluted neural networks

We consider pulse-coupled Leaky Integrate-and-Fire neural networks with randomly distributed synaptic couplings. This random dilution induces fluctuations in the evolution of the macroscopic variables and deterministic chaos at the microscopic level. Our main aim is to mimic the effect of the dilution as a noise source acting on the dynamics of a globally coupled non-chaotic system. Indeed, the evolution of a diluted neural network can be well approximated as a fully pulse coupled network, where each neuron is driven by a mean synaptic current plus additive noise. These terms represent the average and the fluctuations of the synaptic currents acting on the single neurons in the diluted system. The main microscopic and macroscopic dynamical features can be retrieved with this stochastic approximation. Furthermore, the microscopic stability of the diluted network can be also reproduced, as demonstrated from the almost coincidence of the measured Lyapunov exponents in the deterministic and stochastic cases for an ample range of system sizes. Our results strongly suggest that the fluctuations in the synaptic currents are responsible for the emergence of chaos in this class of pulse coupled networks.Comment: 12 Pages, 4 Figure