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.

Primary stability of two uncemented acetabular components of different geometry : hemispherical or peripherally enhanced?

Objective This study compared the primary stability of two commercially available acetabular components from the same manufacturer, which differ only in geometry; a hemispherical and a peripherally enhanced design (peripheral self-locking (PSL)). The objective was to determine whether altered geometry resulted in better primary stability. Methods Acetabular components were seated with 0.8 mm to 2 mm interference fits in reamed polyethylene bone substrate of two different densities (0.22 g/cm3 and 0.45 g/cm3). The primary stability of each component design was investigated by measuring the peak failure load during uniaxial pull-out and tangential lever-out tests. Results There was no statistically significant difference in seating force (p = 0.104) or primary stability (pull-out p = 0.171, lever-out p = 0.087) of the two components in the low-density substrate. Similarly, in the high-density substrate, there was no statistically significant difference in the peak pull-out force (p = 0.154) or lever-out moment (p = 0.574) between the designs. However, the PSL component required a significantly higher seating force thanthe hemispherical cup in the high-density bone analogue (p = 0.006). Conclusions Higher seating forces associated with the PSL design may result in inadequate seating and increased risk of component malpositioning or acetabular fracture in the intra-operative setting in high-density bone stock. Our results, if translated clinically, suggest that a purely hemispherical geometry may have an advantage over a peripherally enhanced geometry in high density bone stock.

Evaluation of two-view geometry methods with automatic ground-truth generation

A large number of methods have been published that aim to evaluate various components of multi-view geometry systems. Most of these have focused on the feature extraction, description and matching stages (the visual front end), since geometry computation can be evaluated through simulation. Many data sets are constrained to small scale scenes or planar scenes that are not challenging to new algorithms, or require special equipment. This paper presents a method for automatically generating geometry ground truth and challenging test cases from high spatio-temporal resolution video. The objective of the system is to enable data collection at any physical scale, in any location and in various parts of the electromagnetic spectrum. The data generation process consists of collecting high resolution video, computing accurate sparse 3D reconstruction, video frame culling and down sampling, and test case selection. The evaluation process consists of applying a test 2-view geometry method to every test case and comparing the results to the ground truth. This system facilitates the evaluation of the whole geometry computation process or any part thereof against data compatible with a realistic application. A collection of example data sets and evaluations is included to demonstrate the range of applications of the proposed system.

Analysis of terrain geometry representations for traversability of a Mars rover

For a planetary rover to successfully traverse across unstructured terrain autonomously, one of the major challenges is to assess its local traversability such that it can plan a trajectory to achieve its mission goals efficiently while minimising risk to the vehicle itself. This paper aims to provide a comparative study on different approaches for representing the geometry of Martian terrain for the purpose of evaluating terrain traversability. An accurate representation of the geometric properties of the terrain is essential as it can directly affect the determination of traversability for a ground vehicle. We explore current state-of-the-art techniques for terrain estimation, in particular Gaussian Processes (GP) in various forms, and discuss the suitability of each technique in the context of an unstructured Martian terrain. Furthermore, we present the limitations of regression techniques in terms of spatial correlation and continuity assumptions, and the impact on traversability analysis of a planetary rover across unstructured terrain. The analysis was performed on datasets of the Mars Yard at the Powerhouse Museum in Sydney, obtained using the onboard RGB-D camera.

Are in vitro estimates of cell diffusivity and cell proliferation rate sensitive to assay geometry?

Cells respond to various biochemical and physical cues during wound–healing and tumour progression. In vitro assays used to study these processes are typically conducted in one particular geometry and it is unclear how the assay geometry affects the capacity of cell populations to spread, or whether the relevant mechanisms, such as cell motility and cell proliferation, are somehow sensitive to the geometry of the assay. In this work we use a circular barrier assay to characterise the spreading of cell populations in two different geometries. Assay 1 describes a tumour–like geometry where a cell population spreads outwards into an open space. Assay 2 describes a wound–like geometry where a cell population spreads inwards to close a void. We use a combination of discrete and continuum mathematical models and automated image processing methods to obtain independent estimates of the effective cell diffusivity, D, and the effective cell proliferation rate, λ. Using our parameterised mathematical model we confirm that our estimates of D and λ accurately predict the time–evolution of the location of the leading edge and the cell density profiles for both assay 1 and assay 2. Our work suggests that the effective cell diffusivity is up to 50% lower for assay 2 compared to assay 1, whereas the effective cell proliferation rate is up to 30% lower for assay 2 compared to assay 1.

Monitoring phase behavior of sub- and supercritical CO2 confined in porous fractal silica with 85% porosity

Phase behavior of CO2 confined in porous fractal silica with volume fraction of SiO2 φs = 0.15 was investigated using small-angle neutron scattering (SANS) and ultrasmall-angle neutron scattering (USANS) techniques. The range of fluid densities (0<(FCO2)bulk<0.977 g/cm3) and temperatures (T=22 °C, 35 and 60 °C) corresponded to gaseous, liquid, near critical and supercritical conditions of the bulk fluid. The results revealed formation of a dense adsorbed phase in small pores with sizes D<40 A° at all temperatures. At low pressure (P <55 bar, (FCO2)bulk <0.2 g/cm3) the average fluid density in pores may exceed the density of bulk fluid by a factor up to 6.5 at T=22 °C. This “enrichment factor” gradually decreases with temperature, however significant fluid densification in small pores still exists at temperature T=60°C, i.e., far above the liquid-gas critical temperature of bulk CO2 (TC=31.1 °C). Larger pores are only partially filled with liquid-like adsorbed layer which coexists with unadsorbed fluid in the pore core. With increasing pressure, all pores become uniformly filled with the fluid, showing no measurable enrichment or depletion of the porous matrix with CO2.

Object tracking via non-Euclidean geometry : a Grassmann approach

A robust visual tracking system requires an object appearance model that is able to handle occlusion, pose, and illumination variations in the video stream. This can be difficult to accomplish when the model is trained using only a single image. In this paper, we first propose a tracking approach based on affine subspaces (constructed from several images) which are able to accommodate the abovementioned variations. We use affine subspaces not only to represent the object, but also the candidate areas that the object may occupy. We furthermore propose a novel approach to measure affine subspace-to-subspace distance via the use of non-Euclidean geometry of Grassmann manifolds. The tracking problem is then considered as an inference task in a Markov Chain Monte Carlo framework via particle filtering. Quantitative evaluation on challenging video sequences indicates that the proposed approach obtains considerably better performance than several recent state-of-the-art methods such as Tracking-Learning-Detection and MILtrack.

Surface magnetoplasma waves at the interface between a plasma-like medium and a metal in the Voigt geometry

Theoretical and experimental results associated with the studies of different properties of surface-type waves (SW) in plasma-like medium-metal structures are reviewed. The propagation of surface waves in the Voigt geometry (the SW propagate across the external magnetic field, which is parallel to the interface) is considered. Various problems dealing with the linear properties of the SW (dispersion characteristics, electromagnetic fields topography, influence of the inhomogeneity of the medium, etc.); excitation mechanisms of the plasma-metal waveguide structures (parametric, drift, diffraction, etc. mechanisms); nonlinear effects associated with SW propagation (higher harmonics generation, self-interaction, nonlinear damping, nonlinear interactions, etc.) are presented. In many cases the results are valid for both gaseous and solid-state plasmas. © 1999 Elsevier Science B.V. All rights reserved.

Three-dimensional geological modelling of the Galilee and central Eromanga basins, Australia: New insights into aquifer/aquitard geometry and potential influence of faults on inter-connectivity

Study region The Galilee and Eromanga basins are located in central Queensland, Australia. Both basins are components of the Great Artesian Basin which host some of the most significant groundwater resources in Australia. Study focus This study evaluates the influence of regional faults on groundwater flow in an aquifer/aquitard interbedded succession that form one of the largest Artesian Basins in the world. In order to assess the significance of regional faults as potential barriers or conduits to groundwater flow, vertical displacements of the major aquifers and aquitards were studied at each major fault and the general hydraulic relationship of units that are juxtaposed by the faults were considered. A three-dimensional (3D) geological model of the Galilee and Eromanga basins was developed based on integration of well log data, seismic surfaces, surface geology and elevation data. Geological structures were mapped in detail and major faults were characterised. New hydrological insights for the region Major faults that have been described in previous studies have been confirmed within the 3D geological model domain and a preliminary assessment of their hydraulic significance has been conducted. Previously unknown faults such as the Thomson River Fault (herein named) have also been identified in this study.

Formation of the three-dimensional geometry of the red blood cell membrane

Red blood cells (RBCs) are nonnucleated liquid capsules, enclosed in deformable viscoelastic membranes with complex three dimensional geometrical structures. Generally, RBC membranes are highly incompressible and resistant to areal changes. However, RBC membranes show a planar shear deformation and out of plane bending deformation. The behaviour of RBCs in blood vessels is investigated using numerical models. All the characteristics of RBC membranes should be addressed to develop a more accurate and stable model. This article presents an effective methodology to model the three dimensional geometry of the RBC membrane with the aid of commercial software COMSOL Multiphysics 4.2a and Fortran programming. Initially, a mesh is generated for a sphere using the COMSOL Multiphysics software to represent the RBC membrane. The elastic energy of the membrane is considered to determine a stable membrane shape. Then, the actual biconcave shape of the membrane is obtained based on the principle of virtual work, when the total energy is minimised. The geometry of the RBC membrane could be used with meshfree particle methods to simulate motion and deformation of RBCs in micro-capillaries