Time-dependent fractional advection–diffusion equations by an implicit MLS meshless method

Recently, many new applications in engineering and science are governed by a series of fractional partial differential equations (FPDEs). Unlike the normal partial differential equations (PDEs), the differential order in a FPDE is with a fractional order, which will lead to new challenges for numerical simulation, because most existing numerical simulation techniques are developed for the PDE with an integer differential order. The current dominant numerical method for FPDEs is Finite Difference Method (FDM), which is usually difficult to handle a complex problem domain, and also hard to use irregular nodal distribution. This paper aims to develop an implicit meshless approach based on the moving least squares (MLS) approximation for numerical simulation of fractional advection-diffusion equations (FADE), which is a typical FPDE. The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless strong-forms. The stability and convergence related to the time discretization of this approach are then discussed and theoretically proven. Several numerical examples with different problem domains and different nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of the FADE.

Nanostructured mesoporous bioglass coated zirconia scaffolds for bone tissue engineering

For the filling and reconstruction of non-healing bone defects, the application of porous ceramic scaffold as bone substitutes is considered to be a reasonable choice. In bone tissue engineering, an ideal scaffold must satisfy several criterias such as open porosity, having high compressive strength (it depends where in body, and if external fixatures are used) and the practicability for cell migration. Many researchers have focused on enhancing the mechanical properties of hydroxyapatite scaffolds by combining it with other biomaterials, such as bioglass and polymers. Nevertheless, there is still a lack of suitable scaffolds based on porous biomaterials. In this study, zirconia scaffolds from two different templates (polyurethane (PU) and Acrylonitrile Butadiene Styrene (ABS) templates) were successfully fabricated with dissimilar fabrication techniques. The scaffold surfaces were further modified with mesoporous bioglass for the purpose of bone tissue engineering. In the study of PU template scaffold, high porosity (~88%) sol-gel derived yttria-stabilized zirconia (YSZ) scaffold was prepared by a polyurethane (PU) foam replica method using sol-gel derived zirconia for the first time, and double coated with Mesoporous Bioglass (MBGs) coating. For the ABS template scaffold, two types of templates (cube and cylinder) with different strut spacings were used and fabricated by a 3D Rapid Prototyper. Subsequently, zirconia scaffolds with low porosity (63±2.8% to 68±2.5%) were fabricated by embedding the zirconia powder slurry into the ABS templates and burning out the ABS to produce a uniform porous structure. The zirconia scaffolds were double coated with mesoporous bioglass by dip coating for the first time. The porosities of the scaffolds were calculated before and after coating. The microstructures were then examined using scanning electron microscopy and the mechanical properties were evaluated using compressive test. Accordingly, relationships between microstructure, processing and mechanical behaviour of the porous zirconia was discussed. Scaffold biocompatibility and bioactivity was also evaluated using a bone marrow stromal cell (BMSC) proliferation test and a simulated body fluid test.

Preparation and characterisation of tri-calcium phosphate scaffolds with tunnel-like macro-pores for bone tissue engineering

Calcium Phosphate ceramics have been widely used in tissue engineering due to their excellent biocompatibility and biodegradability. In the physiological environment, they are able to gradually degrade, absorbed and promote bone growth. Ultimately, they are capable of replacing damaged bone with new tissue. However, their low mechanical properties limit calcium phosphate ceramics in load-bearing applications. To obtain sufficient mechanical properties as well as high biocompatibility is one of the main focuses in biomaterials research. Therefore, the current project focuses on the preparation and characterization of porous tri-calcium phosphate (TCP) ceramic scaffolds. Hydroxapatite (HA) was used as the raw material, and normal calcium phosphate bioglass was added to adjust the ratio between calcium and phosphate. It was found that when 20% bioglass was added to HA and sintered at 1400oC for 3 hours, the TCP scaffold was obtained and this was confirmed by X-ray diffraction (XRD) analysis. Test results have shown that by applying this method, TCP scaffolds have significantly higher compressive strength (9.98MPa) than those made via TCP powder (<3MPa). Moreover, in order to further increase the compressive strength of TCP scaffolds, the samples were then coated with bioglass. For normal bioglass coated TCP scaffold, compressive strength was 16.69±0.5MPa; the compressive strength for single layer mesoporous bioglass coated scaffolds was 15.03±0.63MPa. In addition, this project has also concentrated on sizes and shapes effects; it was found that the cylinder scaffolds have more mechanical property than the club ones. In addition, this project performed cell culture within scaffold to assess biocompatibility. The cells were well distributed in the scaffold, and the cytotoxicity test was performed by 3-(4,5)-dimethylthiahiazo(-z-y1)-3,5-di- phenytetrazoliumromide (MTT) assay. The Alkaline Phosphatase (Alp) activity of human bone marrow mesenchymal stem cell system (hBMSCs) seeded on scaffold expressed higher in vitro than that in the positive control groups in osteogenic medium, which indicated that the scaffolds were both osteoconductive and osteoinductive, showing potential value in bone tissue engineering.

Analysis of rural activity spaces and transport disadvantage using a multi-method approach

Current knowledge about the relationship between transport disadvantage and activity space size is limited to urban areas, and as a result, very little is known about this link in a rural context. In addition, although research has identified transport disadvantaged groups based on their size of activity space, these studies have, however, not empirically explained such differences and the result is often a poor identification of the problems facing disadvantaged groups. Research has shown that transport disadvantage varies over time. The static nature of analysis using the activity space concept in previous research studies has lacked the ability to identify transport disadvantage in time. Activity space is a dynamic concept; and therefore possesses a great potential in capturing temporal variations in behaviour and access opportunities. This research derives measures of the size and fullness of activity spaces for 157 individuals for weekdays, weekends, and for a week using weekly activity-travel diary data from three case study areas located in rural Northern Ireland. Four focus groups were also conducted in order to triangulate quantitative findings and to explain the differences between different socio-spatial groups. The findings of this research show that despite having a smaller sized activity space, individuals were not disadvantaged because they were able to access their required activities locally. Car-ownership was found to be an important life line in rural areas. Temporal disaggregation of the data reveals that this is true only on weekends due to a lack of public transport services. In addition, despite activity spaces being at a similar size, the fullness of activity spaces of low-income individuals was found to be significantly lower compared to their high-income counterparts. Focus group data shows that financial constraint, poor connections both between public transport services and between transport routes and opportunities forced individuals to participate in activities located along the main transport corridors.

Efficient hybrid-game strategies coupled to evolutionary algorithms for robust multidisciplinary design optimization in aerospace engineering

A number of game strategies have been developed in past decades and used in the fields of economics, engineering, computer science, and biology due to their efficiency in solving design optimization problems. In addition, research in multiobjective and multidisciplinary design optimization has focused on developing a robust and efficient optimization method so it can produce a set of high quality solutions with less computational time. In this paper, two optimization techniques are considered; the first optimization method uses multifidelity hierarchical Pareto-optimality. The second optimization method uses the combination of game strategies Nash-equilibrium and Pareto-optimality. This paper shows how game strategies can be coupled to multiobjective evolutionary algorithms and robust design techniques to produce a set of high quality solutions. Numerical results obtained from both optimization methods are compared in terms of computational expense and model quality. The benefits of using Hybrid and non-Hybrid-Game strategies are demonstrated.

Adaptive wing/aerofoil design optimisation using MOEA coupled to uncertainty design method

The use of adaptive wing/aerofoil designs is being considered as promising techniques in aeronautic/aerospace since they can reduce aircraft emissions, improve aerodynamic performance of manned or unmanned aircraft. The paper investigates the robust design and optimisation for one type of adaptive techniques; Active Flow Control (AFC) bump at transonic flow conditions on a Natural Laminar Flow (NLF) aerofoil designed to increase aerodynamic efficiency (especially high lift to drag ratio). The concept of using Shock Control Bump (SCB) is to control supersonic flow on the suction/pressure side of NLF aerofoil: RAE 5243 that leads to delaying shock occurrence or weakening its strength. Such AFC technique reduces total drag at transonic speeds due to reduction of wave drag. The location of Boundary Layer Transition (BLT) can influence the position the supersonic shock occurrence. The BLT position is an uncertainty in aerodynamic design due to the many factors, such as surface contamination or surface erosion. The paper studies the SCB shape design optimisation using robust Evolutionary Algorithms (EAs) with uncertainty in BLT positions. The optimisation method is based on a canonical evolution strategy and incorporates the concepts of hierarchical topology, parallel computing and asynchronous evaluation. Two test cases are conducted; the first test assumes the BLT is at 45% of chord from the leading edge and the second test considers robust design optimisation for SCB at the variability of BLT positions and lift coefficient. Numerical result shows that the optimisation method coupled to uncertainty design techniques produces Pareto optimal SCB shapes which have low sensitivity and high aerodynamic performance while having significant total drag reduction.

An improved genetic algorithm and graph theory based method for optimal sectionalizer switch placement in distribution networks with DG

In this paper a new graph-theory and improved genetic algorithm based practical method is employed to solve the optimal sectionalizer switch placement problem. The proposed method determines the best locations of sectionalizer switching devices in distribution networks considering the effects of presence of distributed generation (DG) in fitness functions and other optimization constraints, providing the maximum number of costumers to be supplied by distributed generation sources in islanded distribution systems after possible faults. The proposed method is simulated and tested on several distribution test systems in both cases of with DG and non DG situations. The results of the simulations validate the proposed method for switch placement of the distribution network in the presence of distributed generation.

An advanced meshless method for time fractional diffusion equation

Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of fractional partial differential equations.

A node-based smoothed conforming point interpolation method (NS-CPIM) for elasticity problems

This paper formulates a node-based smoothed conforming point interpolation method (NS-CPIM) for solid mechanics. In the proposed NS-CPIM, the higher order conforming PIM shape functions (CPIM) have been constructed to produce a continuous and piecewise quadratic displacement field over the whole problem domain, whereby the smoothed strain field was obtained through smoothing operation over each smoothing domain associated with domain nodes. The smoothed Galerkin weak form was then developed to create the discretized system equations. Numerical studies have demonstrated the following good properties: NS-CPIM (1) can pass both standard and quadratic patch test; (2) provides an upper bound of strain energy; (3) avoid the volumetric locking; (4) provides the higher accuracy than those in the node-based smoothed schemes of the original PIMs.

A nanoscratch method for measuring hardness of thin films

Thin solid films were extensively used in the making of solar cells, cutting tools, magnetic recording devices, etc. As a result, the accurate measurement of mechanical properties of the thin films, such as hardness and elastic modulus, was required. The thickness of thin films normally varies from tens of nanometers to several micrometers. It is thus challenging to measure their mechanical properties. In this study, a nanoscratch method was proposed for hardness measurement. A three-dimensional finite element method (3-D FEM) model was developed to validate the nanoscratch method and to understand the substrate effect during nanoscratch. Nanoindentation was also used for comparison. The nanoscratch method was demonstrated to be valuable for measuring hardness of thin solid films.

An enriched radial point interpolation method based on weak-form and strong-form

In this article, an enriched radial point interpolation method (e-RPIM) is developed for computational mechanics. The conventional radial basis function (RBF) interpolation is novelly augmented by the suitable basis functions to reflect the natural properties of deformation. The performance of the enriched meshless RBF shape functions is first investigated using the surface fitting. The surface fitting results have proven that, compared with the conventional RBF, the enriched RBF interpolation has a much better accuracy to fit a complex surface than the conventional RBF interpolation. It has proven that the enriched RBF shape function will not only possess all advantages of the conventional RBF interpolation, but also can accurately reflect the deformation properties of problems. The system of equations for two-dimensional solids is then derived based on the enriched RBF shape function and both of the meshless strong-form and weak-form. A numerical example of a bar is presented to study the effectiveness and efficiency of e-RPIM. As an important application, the newly developed e-RPIM, which is augmented by selected trigonometric basis functions, is applied to crack problems. It has been demonstrated that the present e-RPIM is very accurate and stable for fracture mechanics problems.