Predicting occupant risk indicators for frontal impacts with redirective crash cushions

Crash cushions are devices deployed on the road network in order to shield fixed roadside hazards and the non-crashworthy ends of road safety barriers. However crash cushions vary in terms of configuration and operation, meaning that different devices may also vary in terms of ability to mitigate occupant risk. In this study, data derived from crash testing of eleven redirective crash cushions is used as the base input to a numerical procedure for calculation of occupant risk indicators Occupant Impact Velocity (OIV), Occupant Ridedown Acceleration (ORA) and longitudinal Acceleration Severity Index (ASI) for a range of simulated impacting vehicles (mass 800 kg to 2,500 kg) impacting each crash cushion at a range of impact speeds (18 m/s to 32 m/s). The results may be interpreted as demonstrating firstly that enhanced knowledge of the performance of a device over a range of impact conditions, i.e., beyond the crash testing, may assist in determining the crash cushion most suited to a particular application; secondly that a more appropriate conformance test for occupant risk would be a frontal impact by a small (light) vehicle travelling parallel to and aligned with the centreline of the crash cushion; and thirdly that current documented numerical procedures for calculating occupant risk indicators may require review.

Numerical modelling of particle deposition in idealized porous channels

This project provides a steppingstone to comprehend the mechanisms that govern particulate fouling in metal foam heat exchangers. The method is based on development of an advanced Computational Fluid Dynamics model in addition to performing analytical validation. This novel method allows an engineer to better optimize heat exchanger designs, thereby mitigating fouling, reducing energy consumption caused by fouling, economize capital expenditure on heat exchanger maintenance, and reduce operation downtime. The robust model leads to the establishment of an alternative heat exchanger configuration that has lower pressure drop and particulate deposition propensity.

Implementation of parallel tridiagonal solvers for a heterogeneous computing environment

Tridiagonal diagonally dominant linear systems arise in many scientific and engineering applications. The standard Thomas algorithm for solving such systems is inherently serial forming a bottleneck in computation. Algorithms such as cyclic reduction and SPIKE reduce a single large tridiagonal system into multiple small independent systems which can be solved in parallel. We have developed portable cyclic reduction and SPIKE algorithm OpenCL implementations with the intent to target a range of co-processors in a heterogeneous computing environment including Field Programmable Gate Arrays (FPGAs), Graphics Processing Units (GPUs) and other multi-core processors. In this paper, we evaluate these designs in the context of solver performance, resource efficiency and numerical accuracy.

Analytical model of reactive transport processes with spatially variable coefficients

Analytical solutions of partial differential equation (PDE) models describing reactive transport phenomena in saturated porous media are often used as screening tools to provide insight into contaminant fate and transport processes. While many practical modelling scenarios involve spatially variable coefficients, such as spatially variable flow velocity, v(x), or spatially variable decay rate, k(x), most analytical models deal with constant coefficients. Here we present a framework for constructing exact solutions of PDE models of reactive transport. Our approach is relevant for advection-dominant problems, and is based on a regular perturbation technique. We present a description of the solution technique for a range of one-dimensional scenarios involving constant and variable coefficients, and we show that the solutions compare well with numerical approximations. Our general approach applies to a range of initial conditions and various forms of v(x) and k(x). Instead of simply documenting specific solutions for particular cases, we present a symbolic worksheet, as supplementary material, which enables the solution to be evaluated for different choices of the initial condition, v(x) and k(x). We also discuss how the technique generalizes to apply to models of coupled multispecies reactive transport as well as higher dimensional problems.

Numerical studies of gypsum plasterboard and MGO board lined LSF walls exposed to fire

Fire resistance of cold-formed light gauge steel frame (LSF) wall systems is enhanced by lining them with single or multiple layers of wall boards with varying thermal properties. These wall boards are gypsum plasterboards or Magnesium Oxide (MgO) boards produced by different manufacturers. Thermal properties of these boards appear to show considerable variations and this can lead to varying fire resistance levels (FRL) for their wall systems. Currently FRLs of wall systems are determined using full scale fire tests, but they are time consuming and expensive. Recent research studies on the fire performance of LSF wall systems have used finite element studies to overcome this problem, but they were developed based on 1-D and 2-D finite element platform capable of performing either heat transfer or structural analysis separately. Hence in this research a 3-D finite element model was developed first for LSF walls lined with gypsum plasterboard and cavity insulation materials. Accurate thermal properties of these boards are essential for finite element modelling, and thus they were measured at both ambient and elevated temperatures. This experimental study included specific heat, relative density and thermal conductivity of boards. The developed 3-D finite element model was then validated using the available fire tests results of LSF walls lined with gypsum plasterboard, and is being used to investigate the fire performance of different LSF wall configurations. The tested MgO board exhibited significant variations in their thermal properties in comparison to gypsum plasterboards with about 50% loss of its initial mass at about 500 ºC compared to 16% for gypsum plasterboards. Hence the FRL of MgO board lined LSF wall systems is likely to be significantly reduced. This paper presents the details of this research study on the fire performance of LSF wall systems lined with gypsum plasterboard and MgO board including the developed 3-D finite element models, thermal property tests and the results.

Numerical modelling of deformation behaviour of red blood cells in microvessels using the coupled SPH-DEM method

This thesis developed an advanced computational model to investigate the motion and deformation properties of red blood cells in capillaries. The novel model is based on the meshfree particle methods and is capable of modelling the large deformation of red blood cells moving through blood vessels. The developed model was employed to simulate the deformation behaviour of healthy and malaria infected red blood cells as well as the motion of red blood cells in stenosed capillaries.

The clinical effectiveness and cost-effectiveness of exercise referral schemes: A systematic review and economic evaluation

- Background Exercise referral schemes (ERS) aim to identify inactive adults in the primary-care setting. The GP or health-care professional then refers the patient to a third-party service, with this service taking responsibility for prescribing and monitoring an exercise programme tailored to the needs of the individual. - Objective To assess the clinical effectiveness and cost-effectiveness of ERS for people with a diagnosed medical condition known to benefit from physical activity (PA). The scope of this report was broadened to consider individuals without a diagnosed condition who are sedentary. - Data sources MEDLINE; EMBASE; PsycINFO; The Cochrane Library, ISI Web of Science; SPORTDiscus and ongoing trial registries were searched (from 1990 to October 2009) and included study references were checked. - Methods Systematic reviews: the effectiveness of ERS, predictors of ERS uptake and adherence, and the cost-effectiveness of ERS; and the development of a decision-analytic economic model to assess cost-effectiveness of ERS. - Results Seven randomised controlled trials (UK, n = 5; non-UK, n = 2) met the effectiveness inclusion criteria, five comparing ERS with usual care, two compared ERS with an alternative PA intervention, and one to an ERS plus a self-determination theory (SDT) intervention. In intention-to-treat analysis, compared with usual care, there was weak evidence of an increase in the number of ERS participants who achieved a self-reported 90-150 minutes of at least moderate-intensity PA per week at 6-12 months' follow-up [pooled relative risk (RR) 1.11, 95% confidence interval 0.99 to 1.25]. There was no consistent evidence of a difference between ERS and usual care in the duration of moderate/vigorous intensity and total PA or other outcomes, for example physical fitness, serum lipids, health-related quality of life (HRQoL). There was no between-group difference in outcomes between ERS and alternative PA interventions or ERS plus a SDT intervention. None of the included trials separately reported outcomes in individuals with medical diagnoses. Fourteen observational studies and five randomised controlled trials provided a numerical assessment of ERS uptake and adherence (UK, n = 16; non-UK, n = 3). Women and older people were more likely to take up ERS but women, when compared with men, were less likely to adhere. The four previous economic evaluations identified suggest ERS to be a cost-effective intervention. Indicative incremental cost per quality-adjusted life-year (QALY) estimates for ERS for various scenarios were based on a de novo model-based economic evaluation. Compared with usual care, the mean incremental cost for ERS was £169 and the mean incremental QALY was 0.008, with the base-case incremental cost-effectiveness ratio at £20,876 per QALY in sedentary people without a medical condition and a cost per QALY of £14,618 in sedentary obese individuals, £12,834 in sedentary hypertensive patients, and £8414 for sedentary individuals with depression. Estimates of cost-effectiveness were highly sensitive to plausible variations in the RR for change in PA and cost of ERS. - Limitations We found very limited evidence of the effectiveness of ERS. The estimates of the cost-effectiveness of ERS are based on a simple analytical framework. The economic evaluation reports small differences in costs and effects, and findings highlight the wide range of uncertainty associated with the estimates of effectiveness and the impact of effectiveness on HRQoL. No data were identified as part of the effectiveness review to allow for adjustment of the effect of ERS in different populations. - Conclusions There remains considerable uncertainty as to the effectiveness of ERS for increasing activity, fitness or health indicators or whether they are an efficient use of resources in sedentary people without a medical diagnosis. We failed to identify any trial-based evidence of the effectiveness of ERS in those with a medical diagnosis. Future work should include randomised controlled trials assessing the cinical effectiveness and cost-effectivenesss of ERS in disease groups that may benefit from PA. - Funding The National Institute for Health Research Health Technology Assessment programme.

Experimental and numerical investigation of the behaviour of CFRP strengthened CHS beams subjected to bending

This paper presents the results of an experimental and numerical program to investigate the circular hollow section (CHS) beams, strengthened using Carbon Fibre Reinforced Polymer (CFRP) sheets. The circular hollow shaped steel beams bonded with different CFRP layer orientations were tested under four-point bending. The mid-span deflection, service load and failure load were recorded. The LHL (where L, first inner longitudinal layer, H, second hoop layer and L, third outer longitudinal layer) and LLH (where L, first inner longitudinal layer, L, second longitudinal layer and H, third outer hoop layer) layer oriented strengthened beams perform slightly better than HHL (where H, first inner hoop layer, H, second hoop layer and L, third outer longitudinal layer) layer oriented strengthened beams. The LHL and LLH layer oriented treated beams showed very similar structural behaviour. Numerical analyses were then conducted on the CFRP strengthened steel CHS beams. The validity of the models has been assessed by comparing the failure loads and mid-span deflections. The effects of various parameters such as bond length, section types, tensile modulus of CFRP, adhesive layer thickness and adhesive types have been studied.

Experimental and numerical investigation of strain-rate-dependent behavior of kangaroo shoulder cartilage

This thesis introduces a new animal model, kangaroo, to biomechanical investigations of shoulder cartilage research. It examines the effect of cartilage structure and constituents on tissue behavior and its adaptation to mechanical loading. In doing so, the study explains the relationship of tissue's functional behaviors to its structure and constituents which has important implications for tissue engineering strategies catering joint specific cartilage tissue generation.

Mathematical modelling and numerical simulation of phosphine flow during grain fumigation in leaky cylindrical silos

The phosphine distribution in a cylindrical silo containing grain is predicted. A three-dimensional mathematical model, which accounts for multicomponent gas phase transport and the sorption of phosphine into the grain kernel is developed. In addition, a simple model is presented to describe the death of insects within the grain as a function of their exposure to phosphine gas. The proposed model is solved using the commercially available computational fluid dynamics (CFD) software, FLUENT, together with our own C code to customize the solver in order to incorporate the models for sorption and insect extinction. Two types of fumigation delivery are studied, namely, fan- forced from the base of the silo and tablet from the top of the silo. An analysis of the predicted phosphine distribution shows that during fan forced fumigation, the position of the leaky area is very important to the development of the gas flow field and the phosphine distribution in the silo. If the leak is in the lower section of the silo, insects that exist near the top of the silo may not be eradicated. However, the position of a leak does not affect phosphine distribution during tablet fumigation. For such fumigation in a typical silo configuration, phosphine concentrations remain low near the base of the silo. Furthermore, we find that half-life pressure test readings are not an indicator of phosphine distribution during tablet fumigation.

Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

Numerical study on the surface depression of the concrete runway pavement under impact load

Airport runway pavement always subjected to huge impact loading due to the hard landing of aircraft on the pavement surface. Therefore runway pavements should have sufficient impact resistance capability to avoid damage causing by hard impact like surface deflection in downward or penetration since the repair works is cumbersome within the operating condition of airport and also increases the service life cost of the pavement structure. Several research works have been carried out on airport runway pavement to measure the present condition of pavement and also to predict future performance of it. However, most of the works are confined by pavement response under moving aircraft loading. Nevertheless, no comprehensive research work is yet conducted to identify the controlling factors which might have significant effect in changing the common pavements damage like surface penetration depth under impact of aircraft. Therefore, a 3D FE study is conducted to determine some effective factors in controlling the top surface penetration depth of runway pavement. Among the exterior factors, mass of the impactor, velocity of the impactor, impact angle and boundary conditions are selected and as interior factors, thickness of the runway pavement, compressive strength and density of materials used in the runway pavement are selected.

Experimental and numerical studies of hollow flange channel beams subject to web crippling under ETF and ITF load cases

This paper presents the details of experimental and numerical studies on the web crippling behaviour of hollow flange channel beams, known as LiteSteel beams (LSB). The LSB has a unique shape of a channel beam with two rectangular hollow flanges, made using a unique manufacturing process. Experimental and numerical studies have been carried out to evaluate the behaviour and design of LSBs subject to pure bending actions, predominant shear actions and combined actions. To date, however, no investigation has been conducted into the web crippling behaviour and strength of LSB sections under ETF and ITF load conditions. Hence experimental studies consisting of 28 tests were first conducted in this research to assess the web crippling behaviour and strengths of LSBs under two flange load cases (ETF and ITF). Experimental web crippling capacity results were then compared with the predictions from AS/NZS 4600 and AISI S100 design rules, which showed that AS/NZS 4600 and AISI S100 design equations are very unconservative for LSBs under ETF and ITF load cases. Hence improved equations were proposed to determine the web crippling capacities of LSBs. Finite element models of the tested LSBs were then developed, and used to determine the elastic buckling loads of LSBs under ETF and ITF load cases. New equations were proposed to determine the corresponding elastic buckling coefficients of LSBs. Finally suitable design rules were also developed under the Direct Strength Method format using the test results and buckling analysis results from finite element analyses.

Numerical modelling of rivet fastened rectangular hollow flange channel beams subject to shear

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. 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 experimental and numerical studies were undertaken to investigate the shear behaviour and strength of rivet fastened RHFCBs. In this research finite element models of rivet fastened RHFCBs were developed to investigate their nonlinear shear behaviour including their buckling characteristics and ultimate shear strength. This paper presents the details of the finite element models of rivet fastened RHFCBs and the results. Both finite element analysis and experimental results showed that the current design rules are very conservative for the shear design of rivet fastened RHFCBs. Appropriate improvements have been proposed for the design rules of shear strength of rivet fastened RHFCBs within the Direct Strength Method format.

Thermal performance of magnesium oxide wall board using numerical modelling

Fire resistance of light-gauge steel frame (LSF) walls can be enhanced by lining them with single or multiple layers of wall boards. This research is focused on the thermal per-formance of Magnesium Oxide (MgO) wall boards in comparison to the conventional gypsum plasterboards exposed to standard fire on one side. Thermal properties of MgO board and gypsum plasterboard were measured first and then used in the finite element heat transfer models of the two types of panels. The measured thermal property results show that MgO board will perform better than the gypsum plasterboards due to its higher specific heat values at elevated temperatures. However, MgO board loses 50% of its ini-tial mass at about 500 °C compared to 16% for gypsum plasterboard. The developed finite element models were validated using the fire test results of gypsum plasterboards and then used to study the thermal performance of MgO board panels. Finite element analysis re-sults show that when MgO board panels are exposed to standard fire on one side the rate of temperature rise on the ambient side is significantly reduced compared to gypsum plas-terboard. This has the potential to improve the overall thermal performance of MgO board lined LSF walls and their fire resistance levels (FRL). However, full scale fire tests are needed to confirm this. This paper presents the details of this investigation and the results.