A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

We present a finite volume method to solve the time-space two-sided fractional advection-dispersion equation on a one-dimensional domain. The spatial discretisation employs fractionally-shifted Grünwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes. We demonstrate how the finite volume formulation provides a natural, convenient and accurate means of discretising this equation in conservative form, compared to using a conventional finite difference approach. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

A comparison of finite difference and finite volume methods for solving the space-fractional advection-dispersion equation with variable coefficients

Transport processes within heterogeneous media may exhibit non- classical diffusion or dispersion which is not adequately described by the classical theory of Brownian motion and Fick’s law. We consider a space-fractional advection-dispersion equation based on a fractional Fick’s law. Zhang et al. [Water Resources Research, 43(5)(2007)] considered such an equation with variable coefficients, which they dis- cretised using the finite difference method proposed by Meerschaert and Tadjeran [Journal of Computational and Applied Mathematics, 172(1):65-77 (2004)]. For this method the presence of variable coef- ficients necessitates applying the product rule before discretising the Riemann–Liouville fractional derivatives using standard and shifted Gru ̈nwald formulas, depending on the fractional order. As an alternative, we propose using a finite volume method that deals directly with the equation in conservative form. Fractionally-shifted Gru ̈nwald formulas are used to discretise the Riemann–Liouville fractional derivatives at control volume faces, eliminating the need for product rule expansions. We compare the two methods for several case studies, highlighting the convenience of the finite volume approach.

Multiscale coupling and multiphysics approaches in earth sciences : applications

Geoscientists are confronted with the challenge of assessing nonlinear phenomena that result from multiphysics coupling across multiple scales from the quantum level to the scale of the earth and from femtoseconds to the 4.5 Ga of history of our planet. We neglect in this review electromagnetic modelling of the processes in the Earth’s core, and focus on four types of couplings that underpin fundamental instabilities in the Earth. These are thermal (T), hydraulic (H), mechanical (M) and chemical (C) processes which are driven and controlled by the transfer of heat to the Earth’s surface. Instabilities appear as faults, folds, compaction bands, shear/fault zones, plate boundaries and convective patterns. Convective patterns emerge from buoyancy overcoming viscous drag at a critical Rayleigh number. All other processes emerge from non-conservative thermodynamic forces with a critical critical dissipative source term, which can be characterised by the modified Gruntfest number Gr. These dissipative processes reach a quasi-steady state when, at maximum dissipation, THMC diffusion (Fourier, Darcy, Biot, Fick) balance the source term. The emerging steady state dissipative patterns are defined by the respective diffusion length scales. These length scales provide a fundamental thermodynamic yardstick for measuring instabilities in the Earth. The implementation of a fully coupled THMC multiscale theoretical framework into an applied workflow is still in its early stages. This is largely owing to the four fundamentally different lengths of the THMC diffusion yardsticks spanning micro-metre to tens of kilometres compounded by the additional necessity to consider microstructure information in the formulation of enriched continua for THMC feedback simulations (i.e., micro-structure enriched continuum formulation). Another challenge is to consider the important factor time which implies that the geomaterial often is very far away from initial yield and flowing on a time scale that cannot be accessed in the laboratory. This leads to the requirement of adopting a thermodynamic framework in conjunction with flow theories of plasticity. This framework allows, unlike consistency plasticity, the description of both solid mechanical and fluid dynamic instabilities. In the applications we show the similarity of THMC feedback patterns across scales such as brittle and ductile folds and faults. A particular interesting case is discussed in detail, where out of the fluid dynamic solution, ductile compaction bands appear which are akin and can be confused with their brittle siblings. The main difference is that they require the factor time and also a much lower driving forces to emerge. These low stress solutions cannot be obtained on short laboratory time scales and they are therefore much more likely to appear in nature than in the laboratory. We finish with a multiscale description of a seminal structure in the Swiss Alps, the Glarus thrust, which puzzled geologists for more than 100 years. Along the Glarus thrust, a km-scale package of rocks (nappe) has been pushed 40 km over its footwall as a solid rock body. The thrust itself is a m-wide ductile shear zone, while in turn the centre of the thrust shows a mm-cm wide central slip zone experiencing periodic extreme deformation akin to a stick-slip event. The m-wide creeping zone is consistent with the THM feedback length scale of solid mechanics, while the ultralocalised central slip zones is most likely a fluid dynamic instability.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Web crippling tests of hollow flange channel beams – end one flange and interior one flange load cases

LiteSteel beam (LSB) is a new cold-formed steel hollow flange channel section produced using a patented manufacturing process involving simultaneous cold-forming and dual electric resistance welding. The LSBs were commonly used as floor joists and bearers with web openings in residential, industrial and commercial buildings. Due to the unique geometry of LSBs, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, much of the existing research for common cold-formed steel sections is not directly applicable to LSBs. Many research studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions, predominant shear and combined actions. However, to date, no investigation has been conducted into the web crippling behaviour and strength of LSB sections. Hence detailed experimental studies were conducted to investigate the web crippling behaviour and strengths of LSBs under EOF (End One Flange) and IOF (Interior One Flange) load cases. A total of 26 web crippling tests was conducted and the results were compared with current AS/NZS 4600 design rules. This comparison showed that AS/NZS 4600 (SA, 2005) design rules are very conservative for LSB sections under EOF and IOF load cases. Suitable design equations have been proposed to determine the web crippling capacity of LSBs based on experimental results. This paper presents the details of this experimental study on the web crippling behaviour and strengths of LiteSteel beams under EOF and IOF load cases.

Numerical modelling and design of lipped channel beams subject to shear

Cold-formed steel members are increasingly used as primary structural elements in buildings due to the availability of thin and high strength steels and advanced cold-forming technologies. Cold-formed lipped channel beams (LCB) are commonly used as flexural members such as floor joists and bearers. Many research studies have been carried out to evaluate the behaviour and design of LCBs subject to pure bending actions. However, limited research has been undertaken on the shear behaviour and strength of LCBs. Hence a numerical study was undertaken to investigate the shear behaviour and strength of LCBs. Finite element models of simply supported LCBs with aspect ratios of 1.0 and 1.5 were considered under a mid-span load. They were then validated by comparing their results with test results and used in a detailed parametric study based on the validated finite element models. Numerical studies were conducted to investigate the shear buckling and post-buckling behaviour of LCBs. Experimental and numerical results showed that the current design rules in cold-formed steel structures design codes are very conservative for the shear design of LCBs. Improved design equations were therefore proposed for the shear strength of LCBs. This paper presents the details of this numerical study of LCBs and the results.

Numerical studies and design of hollow flange channel beams subject to combined bending and shear actions

LiteSteel beam (LSB) is a cold-formed steel hollow flange channel section produced using a patented manufacturing process involving simultaneous cold-forming and dual electric resistance welding. It is commonly used as floor joists and bearers in residential, industrial and commercial buildings. Design of the LSB is governed by the Australian cold-formed steel structures code, AS/NZS 4600. Due to the geometry of the LSB, as well as its unique residual stress characteristics and initial geometric imperfections resultant of manufacturing processes, currently available design equations for common cold-formed sections are not directly applicable to the LSB. Many research studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions and predominant shear actions. To date, however, no investigation has been conducted into the strength of LSB sections under combined bending and shear actions. Hence experimental and numerical studies were conducted to assess the combined bending and shear behaviour of LSBs. Finite element models of LSBs were developed to simulate their combined bending and shear behaviour and strength of LSBs. They were then validated by comparing the results with available experimental test results and used in a detailed parametric study. The results from experimental and finite element analyses were compared with current AS/NZS 4600 and AS 4100 design rules. Both experimental and numerical studies show that the AS/NZS 4600 design rule based on circular interaction equation is conservative in predicting the combined bending and shear capacities of LSBs. This paper presents the details of the numerical studies of LSBs and the results. In response to the inadequacies of current approaches to designing LSBs for combined bending and shear, two lower bound design equations are proposed in this paper.

Constrained predictive control of ship fin stabilizers to prevent dynamic stall

In moderate to high sea states, the effectiveness of ship fin stabilizers can severely deteriorate due to nonlinear effects arising from unsteady hydrodynamic characteristics of the fins: dynamic stall. These nonlinear effects take the form of a hysteresis, and they become very significant when the effective angle of attack of the fins exceeds a certain threshold angle. Dynamic stall can result in a complete loss of control action depending on how much the fins exceed the threshold angle. When this is detected, it is common to reduce the gain of the controller that commands the fins. This approach is cautious and tends to reduce performance when the conditions leading to dynamic stall disappear. An alternative approach for preventing the effects while keeping high performance, consists of estimating the effective angle of attack and set a conservative constraint on it as part of the control objectives. In this paper, we investigate the latter approach, and propose the use of a model predictive control (MPC) to prevent the development of these nonlinear effects by considering constraints on both the mechanical angle of the fins and the effective angle of attack.

Lumbar epidural steroids are not cost effective in sciatica

Background Lumbar Epidural Steroids Injections (ESI’s) have previously been shown to provide some degree of pain relief in sciatica. Number Needed To Treat (NNT) to achieve 50% pain relief has been estimated at 7 from the results of randomised controlled trials. Pain relief is temporary. They remain one of the most commonly provided procedures in the UK. It is unknown whether this pain relief represents good value for money. Methods 228 patients were randomised into a multi-centre Double Blind Randomised Controlled Trial. Subjects received up to 3 ESI’s or intra-spinous saline depending on response and fall off with the first injection. All other treatments were permitted. All received a review of analgesia, education and physical therapy. Quality of life was assessed using the SF36 at 6 points and compared using independent sample t-tests. Follow up was up to 1 yr. Missing data was imputed using last observation carried forward (LOCF). QALY’s (Quality of Life Years) were derived from preference based heath values (summary health utility score). SF-6D health state classification was derived from SF-36 raw score data. Standard gambles (SG) were calculated using Model 10. SG scores were calculated on trial results. LOCF was not used for this. Instead average SG were derived for a subset of patients with observations for all visits up to week 12. Incremental QALY’s were derived as the difference in the area between the SG curve for the active group and placebo group. Results SF36 domains showed a significant improvement in pain at week 3 but this was not sustained (mean 54 Active vs 61 Placebo P<0.05). Other domains did not show any significant gains compared with placebo. For derivation of SG the number in the sample in each period differed. In week 12, average SG scores for active and placebo converged. In other words, the health gain for the active group as measured by SG was achieved by the placebo group by week 12. The incremental QALY gained for a patient under the trial protocol compared with the standard care package was 0.0059350. This is equivalent to an additional 2.2 days of full health. The cost per QALY gained to the provider from a patient management strategy administering one epidural as suggested by results was £25 745.68. This result was derived assuming that the gain in QALY data calculated for patients under the trial protocol would approximate that under a patient management strategy based on the trial results (one ESI). This is above the threshold suggested by some as a cost effective treatment. Conclusions The transient benefit in pain relief afforded by ESI’s does not appear to be cost-effective. Further work is needed to develop more cost-effective conservative treatments for sciatica.

Studies on the in-plane shear response of confined masonry shear walls

The prime aim of this research project is to evaluate the performance of confined masonry walls under in-plane shear with a view to contributing to the national masonry design standard through a set of design clauses. This aim stems from the criticisms of the current provisions of the in-plane shear capacity equations in the Australian Masonry Standard AS3700 (2011) being highly non-conservative. This PhD thesis is an attempt to address this gap in the knowledge through systematic investigation of the key parameters that affects the in-plane shear strength of the masonry walls through laboratory experiments and extensive finite element analyses.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early onset scoliosis (EOS). Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option for EOS is fusionless scoliosis surgery. Similar to bracing, this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods are one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semiconstrained growing rods (Medtronic, Sofamor, Danek, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard "constrained / rigid" rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into a 7 level thoracolumbar multi-segment unit (MSU), removing all non-ligamentous soft tissues and leaving 3cm of ribs either side. Pure nondestructive axial rotation moments of ±4Nm at a constant rotation rate of 8deg.s-1 were applied to the mounted MSU spines using a biaxial Instron testing machine. Displacement of each vertebral level was captured using a 3D motion tracking system (Optotrak 3020, Northern Digital Inc, Waterloo, ON). Each spine was tested in an un-instrumented state first and then with appropriately sized semi-constrained growing rods and rigid rods in alternating sequence. The rods were secured by multi-axial pedicle screws (Medtronic CD Horizon) at levels 2 and 6 of the construct. The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm.deg-1) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of the order of testing, rigid rods significantly reduced the total ROM compared with semi-constrained rods (p<0.05) with in a significantly stiffer spine for both left and right axial rotation (p<0.05). Analysing the intervertebral motion within the instrumented levels 2-6, rigid rods showed reduced ROM compared with semi-constrained growing rods and compared with un-instrumented motion segments. CONCLUSION Semi-constrained growing rods maintain similar stiffness in axial rotation to un-instrumented spines, while dual rigid rods significantly reduce axial rotation. Clinically the effect of semi-constrained growing rods as observed in this study is that they would be expected to allow growth via the telescopic rod components while maintaining the axial flexibility of the spine, which may reduce occurrence of the crankshaft phenomenon.

A biomechanical analysis of growing rods used in the management of early onset scoliosis

INTRODUCTION Managing spinal deformities in young children is challenging, particularly early-onset scoliosis (EOS). Any progressive spinal deformity particularly in early life presents significant health risks for the child and a challenge for the treating surgeon. Surgical intervention is often required if EOS has been unresponsive to conservative treatment particularly with rapidly progressive curves. An emerging treatment option particularly for EOS is fusionless scoliosis surgery. Similar to bracing this surgical option potentially harnesses growth, motion and function of the spine along with correcting spinal deformity. Dual growing rods is one such fusionless treatment, which aims to modulate growth of the vertebrae. The aim of this study was to ascertain the extent to which semi-constrained growing rods (Medtronic, Memphis, TN) with a telescopic sleeve component, reduce rotational constraint on the spine compared with standard rigid rods and hence potentially provide a more physiological mechanical environment for the growing spine. METHODS Six 40-60kg English Large White porcine spines served as a model for the paediatric human spine. Each spine was dissected into 7 level thoracolumbar multi-segment unit (MSU) spines, removing all non-ligamentous soft tissues. Appropriately sized semi-constrained growing rods and rigid rods were secured by multi-axial screws (Medtronic) prior to testing in alternating sequences for each spine. Pure nondestructive moments of +/4Nm at a constant rotation rate of 8deg/s was applied to the mounted MSU spines. Displacement of each level was captured using an Optotrak (Northern Digital Inc, Waterloo, ON). The range of motion (ROM), neutral zone (NZ) size and stiffness (Nm/deg) were calculated from the Instron load-displacement data and intervertebral ROM was calculated through a MATLAB algorithm from Optotrak data. RESULTS Irrespective of sequence order rigid rods significantly reduced the total ROM (deg) than compared to semi-constrained rods (p<0.05) and resulted in a significantly stiffer (Nm/deg) spine for both left and right axial rotation testing (p<0.05). Analysing the intervertebral motion within the instrumented levels, rigid rods showed reduced ROM (Deg) than compared to semi-constrained growing rods and the un-instrumented (UN-IN) test sequences. CONCLUSION The semi-constrained growing rods maintained rotation similar to UN-IN spines while the rigid rods showed significantly reduced axial rotation across all instrumented levels. Clinically the effect of semi-constrained growing rods evaluated in this study is that they will allow growth via the telescopic rod components while maintaining the axial rotation ability of the spine, which may also reduce the occurrence of the crankshaft phenomenon.

Shear capacities of rivet fastened Rectangular Hollow Flange Channel Beams

The intermittently rivet fastened Rectangular Hollow Flange Channel Beam (RHFCB) is a new cold-formed hollow section proposed as an alternative to welded hollow flange beams. It is a monosymmetric channel section made by intermittently rivet fastening two torsionally rigid rectangular hollow flanges to a web plate. This method will allow the development of optimum sections by choosing appropriate combinations of web and flange plate widths and thicknesses. RHFCBs can be commonly used as flexural members in buildings. Many experimental and numerical studies have been carried out in the past to investigate the shear behaviour of lipped channel beams. However, no research has been undertaken on the shear behaviour of rivet fastened RHFCBs. Therefore a detailed experimental study involving 19 shear tests was undertaken to investigate the shear behaviour and capacities of rivet fastened RHFCBs. Simply supported test specimens of RHFCB with aspect ratios of 1.0 and 1.5 were loaded at mid-span until failure. Comparison of experimental results with corresponding predictions from the current Australian cold-formed steel design rules showed that the current design rules are very conservative for the shear design of rivet fastened RHFCBs. Significant improvements to web shear buckling occurred due to the presence of rectangular hollow flanges while considerable post-buckling strength was also observed. Appropriate improvements have been proposed for the design rules of shear strength of rivet fastened RHFCBs within the Direct Strength Method format. This paper presents the details of this study and the results.

Web crippling behaviour of C-section beams with web ribs

SupaCee section is one of the cold-formed steel members which is increasingly used in the construction sector. It is characterized by unique ribbed web and curved lip elements, and is claimed to be more economical with extra strength than the traditional channel sections. SupaCee sections are widely used in Australia as floor joists, bearers, purlins and girts. Many experimental and numerical studies have been carried out to evaluate the behaviour and design of conventional channel beams subject to web crippling. To date, however, no investigation has been conducted into the web crippling behaviour and strength of SupaCee sections. Current cold-formed steel design equations do not include any design procedures for SupaCee sections. Hence experimental studies were conducted to assess the web crippling behaviour and strengths of SupaCee sections under ETF and ITF load cases. Thirty six web crippling tests were conducted and the capacity results were compared with the predictions from the AS/NZS 4600 and AISI design rules developed for conventional channel sections. Comparison of ultimate web crippling capacities from tests showed that AS/NZS 4600 and AISI design equations are unconservative for SupaCee sections under ETF load case, but are overly conservative for ITF load case. Hence new equations were proposed to determine the web crippling capacities of SupaCee sections based on the experimental results from this study. Suitable design rules were also developed within the direct strength method format. This paper presents the details of this experimental study of SupaCee sections subject to web crippling and the results.

Experimental investigation of post-fire mechanical properties of cold-formed steels

Cold-formed steel members are widely used in residential, industrial and commercial buildings as primary load-bearing elements. During fire events, they will be exposed to elevated temperatures. If the general appearance of the structure is satisfactory after a fire event then the question that has to be answered is how the load bearing capacity of cold-formed steel members in these buildings has been affected. Hence after such fire events there is a need to evaluate the residual strength of these members. However, the post-fire behaviour of cold-formed steel members has not been investigated in the past. This means conservative decisions are likely to be made in relation to fire exposed cold-formed steel buildings. Therefore an experimental study was undertaken to investigate the post-fire mechanical properties of cold-formed steels. Tensile coupons taken from cold-formed steel sheets of three different steel grades and thicknesses were exposed to different elevated temperatures up to 800 oC, and were then allowed to cool down to ambient temperature before they were tested to failure. Tensile coupon tests were conducted to obtain their post-fire stress-strain curves and associated mechanical properties (yield stress, Young’s modulus, ultimate strength and ductility). It was found that the post-fire mechanical properties of cold-formed steels are reduced below the original ambient temperature mechanical properties if they had been exposed to temperatures exceeding 300 oC. Hence a new set of equations is proposed to predict the post-fire mechanical properties of cold-formed steels. Such post-fire mechanical property assessments allow structural and fire engineers to make an accurate prediction of the safety of fire exposed cold-formed steel buildings. This paper presents the details of this experimental study and the results of post-fire mechanical properties of cold-formed steels. It also includes the results of a post-fire evaluation of cold-formed steel walls.