A RBF meshless approach for modeling a fractal mobile/immobile transport model

A fractional differential equation is used to describe a fractal model of mobile/immobile transport with a power law memory function. This equation is the limiting equation that governs continuous time random walks with heavy tailed random waiting times. In this paper, we firstly propose a finite difference method to discretize the time variable and obtain a semi-discrete scheme. Then we discuss its stability and convergence. Secondly we consider a meshless method based on radial basis functions (RBFs) to discretize the space variable. In contrast to conventional FDM and FEM, the meshless method is demonstrated to have distinct advantages: calculations can be performed independent of a mesh, it is more accurate and it can be used to solve complex problems. Finally the convergence order is verified from a numerical example which is presented to describe a fractal model of mobile/immobile transport process with different problem domains. The numerical results indicate that the present meshless approach is very effective for modeling and simulating fractional differential equations, and it has good potential in the development of a robust simulation tool for problems in engineering and science that are governed by various types of fractional differential equations.

Predicting crack initiation due to ratchetting in rail heads using critical element analysis

This paper presents a strategy to predict the lifetime of rails subjected to large rolling contact loads that induce ratchetting strains in the rail head. A critical element concept is used to calculate the number of loading cycles needed for crack initiation to occur in the rail head surface. In this technique the finite element method (FEM) is used to determine the maximum equivalent ratchetting strain per load cycle, which is calculated by combining longitudinal and shear stains in the critical element. This technique builds on a previously developed critical plane concept that has been used to calculate the number of cycles to crack initiation in rolling contact fatigue under ratchetting failure conditions. The critical element concept simplifies the analytical difficulties of critical plane analysis. Finite element analysis (FEA) is used to identify the critical element in the mesh, and then the strain values of the critical element are used to calculate the ratchetting rate analytically. Finally, a ratchetting criterion is used to calculate the number of cycles to crack initiation from the ratchetting rate calculated.

Temperature redistribution modelling during intermittent microwave convective heating

Microwave power is used for heating and drying processes because of its faster and volumetric heating capability. Non-uniform temperature distribution during microwave application is a major drawback of these processes. Intermittent application of microwave potentially reduces the impact of non-uniformity and improves energy efficiency by redistributing the temperature. However, temperature re-distribution during intermittent microwave heating has not been investigated adequately. Consequently, in this study, a coupled electromagnetic with heat and mass transfer model was developed using the finite element method embedded in COMSOL-Multyphysics software. Particularly, the temperature redistribution due to intermittent heating was investigated. A series of experiments were performed to validate the simulation. The test specimen was an apple and the temperature distribution was closely monitored by a TIC (Thermal Imaging Camera). The simulated temperature profile matched closely with thermal images obtained from experiments.

Direct second-order elastic analysis for steel frame design

The traditional structural design procedure, especially for the large-scale and complex structures, is time consuming and inefficient. This is due primarily to the fact that the traditional design takes the second-order effects indirectly by virtue of design specifications for every member instead of system analysis for a whole structure. Consequently, the complicated and tedious design procedures are inevitably necessary to consider the second-order effects for the member level in design specification. They are twofold in general: 1) Flexural buckling due to P-d effect, i.e. effective length. 2) Sway effect due to P-D effect, i.e. magnification factor. In this study, a new system design concept based on the second-order elastic analysis is presented, in which the second-order effects are taken into account directly in the system analysis, and also to avoid the tedious member-by-member stability check. The plastic design on the basis of this integrated method of direct approach is ignored in this paper for simplicity and clarity, as the only emphasis is placed on the difference between the second-order elastic limit-state design and present system design approach. A practical design example, a 57m-span dome steel skylight structure, is used to demonstrate the efficiency and effectiveness of the proposed approach. This skylight structure is also designed by the traditional design approach BS5950-2000 for comparison on which the emphasis of aforementioned P-d and P-D effects is placed.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Performance-based structural fire design and testing of structures

Fire incident in buildings is common in Hong Kong and this could lead to heavy casualties due to its high population density, so the fire safety design of the framed structure is an important research topic. This paper describes a computer tool for determination of capacity of structural safety against various fire scenarios and the well-accepted second-order direct plastic analysis is adopted for simulation of material yielding and buckling. A computer method is developed to predict structural behaviour of bare steel framed structures at elevated temperatures but the work can be applied to structures made of other materials. These effects of thermal expansion and material degradation due to heating are required to be considered in order to capture the actual behavior of the structure under fire. Degradation of material strength with increasing temperature is included by a set of temperature-stress-strain curves according to BS5950 Part 8 mainly, which implicitly allows for creep deformation. Several numerical and experimental verifications of framed structures are presented and compared against solutions by other researchers. The proposed method allows us to adopt the truly performance-based structural fire analysis and design with significant saving in cost and time.

Blast response and vulnerability assessment of piled foundations

This research treated the response of pile foundations to blast loads. The influence of important parameters was investigated. The research techniques and the results will enable safer design of pile foundations that are vulnerable to blast loads.

Effect of joint mechanism on vehicle redirectional capability of water-filled road safety barrier systems

Portable water-filled barriers (PWFBs) are roadside appurtenances that prevent vehicles from penetrating into temporary construction zones on roadways. PWFBs are required to satisfy the strict regulations for vehicle re-direction in tests. However, many of the current PWFBs fail to re-direct the vehicle at high speeds due to the inability of the joints to provide appropriate stiffness. The joint mechanism hence plays a crucial role in the performance of a PWFB system at high speed impacts. This paper investigates the desired features of the joint mechanism in a PWFB system that can re-direct vehicles at high speeds, while limiting the lateral displacement to acceptable limits. A rectangular “wall” representative of a 30 m long barrier system was modeled and a novel method of joining adjacent road barriers was introduced through appropriate pin-joint connections. The impact response of the barrier “wall” and the vehicle was obtained and the results show that a rotational stiffness of 3000 kNm/rad at the joints seems to provide the desired features of the PWFB system to re-direct impacting vehicles and restrict the lateral deflection. These research findings will be useful to safety engineers and road barrier designers in developing a new generation of PWFBs for increased road safety.

Unsteady natural convection within a porous enclosure of sinusoidal corrugated side walls

Numerically investigation of free convection within a porous cavity with differential heating has been performed using modified corrugated side walls. Sinusoidal hot left and cold right walls are assumed to receive sudden differentially heating where top and bottom walls are insulated. Air is considered as working fluid and is quiescent, initially. Numerical experiments reveal 3 distinct stages of developing pattern including initial stage, oscillatory intermediate and finally steady state condition. Implicit Finite Volume Method with TDMA solver is used to solve the governing equations. This study has been performed for the Rayleigh numbers ranging from 100 to 10,000. Outcomes have been reported in terms of isotherms, streamline, velocity and temperature plots and average Nusselt number for various Ra, corrugation frequency and corrugation amplitude. The effects of sudden differential heating and its resultant transient behavior on fluid flow and heat transfer characteristics have been shown for the range of governing parameters. The present results show that the transient phenomena are enormously influenced by the variation of the Rayleigh Number with corrugation amplitude and frequency.

Optical properties of silver nanowires grown electrochemically in nanoporous alumina

We report fabrication and optical properties of electrochemically deposited silver nanowires into nanoporous alumina template. A finite element method is used to study plasmonic coupling of dipole emitters with the silver nanowires.

Blast response and sensitivity analysis of segmental tunnel

This research treated the response of underground transportation tunnels to surface blast loads using advanced computer simulation techniques. The influences of important parameters, such as tunnel material, geometrical configuration of segments and surrounding soil were investigated. The findings of this research offer significant new information on the blast performance of underground tunnels and will contribute towards future civil engineering applications.

Numerical investigation of deformation and failure mechanisms of osteopontin-hydroxyapatite interfaces and mineralized collagen fibril arrays

This thesis is a comprehensive study of deformation and failure mechanisms in bone at nano- and micro-scale levels. It explores the mechanical behaviour of osteopontin-hydroxyapatite interfaces and mineralized collagen fibril arrays, through atomistic molecular dynamics and finite element simulations. This thesis shows some main factors contributing to the excellent material properties of bone and provides some guidelines for development of new artificial biological materials and medical implants.

Natural convection heat transfer in a baffled triangular enclosure

To reduce the natural convection heat loss from enclosures many researchers used convection suppression devices in the past. In this study a single baffle is used under the top tip to investigate numerically the natural convection heat loss in an attic shaped enclosure which is a cost effective approach. The case considered here is one inclined wall of the enclosure is uniformly heated while the other inclined wall is uniformly cooled with adiabatic bottom wall. The finite volume method has been used to discretize the governing equations, with the QUICK scheme approximating the advection term. The diffusion terms are discretized using central-differencing with second order accuracy. A wide range of governing parameters are studied (Rayleigh number, aspect ratio, baffle length etc.). It is observed that the heat transfer due to natural convection in the enclosure reduces when the baffle length is increased. Effects of other parameters on heat transfer and flow field are described in this study.