A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

External pressure in the hardening of phosphate in tribofilm on iron surfaces

Purpose: In the present work we consider our (in progress) spectroscopy study of zinc and iron phosphates under the influence external high pressure to determine zinc ion change coordination from tetrahedral to octahedral (or hexahedral) structure.----- Design/methodology/approach: The standard equipment is the optical high pressure cell with diamond (DAC). The DAC is assembled and then vibrational or electronic spectra are collected by mounting the cell in an infrared, Raman, EXAFS or UV-visible spectrometer.----- Findings: Mechanism by which zinc and iron methaphosphate material is transformed to glassy meta-phosphate is enhancing mechanical properties of tribofilm. The two decades of intensive study demonstrates that Zn (II) and Fe (III) ions participate to cross-link network under friction, hardening the phosphate.----- Research limitations/implications: Transition metal atoms with d orbital have flexible coordination numbers, for example zinc acts as a cross-linking agent increasing hardness, by changing coordination from tetrahedral to octahedral. Perhaps the external pressure effect on the [Zn–(O-P-)4 ] complex causes a transformation to an [Zn –(O-P-)6] grouping.----- Originality/value: This paper analyses high-pressure spectroscopy which has been applied for the investigation of 3D transition metal ions in solids. When studying pressure effects on coordination compounds structure, we can expect changes in ground electronic state (spin-crossovers), electronic spectra due to structural distortions (piezochromism), and changes in the ligand field causing shifts in the electronic transitions.

Automated road pavement marking detection from high resolution aerial images based on Multi-resolution image analysis and anisotropic Gaussian filtering

Road features extraction from remote sensed imagery has been a long-term topic of great interest within the photogrammetry and remote sensing communities for over three decades. The majority of the early work only focused on linear feature detection approaches, with restrictive assumption on image resolution and road appearance. The widely available of high resolution digital aerial images makes it possible to extract sub-road features, e.g. road pavement markings. In this paper, we will focus on the automatic extraction of road lane markings, which are required by various lane-based vehicle applications, such as, autonomous vehicle navigation, and lane departure warning. The proposed approach consists of three phases: i) road centerline extraction from low resolution image, ii) road surface detection in the original image, and iii) pavement marking extraction on the generated road surface. The proposed method was tested on the aerial imagery dataset of the Bruce Highway, Queensland, and the results demonstrate the efficiency of our approach.

Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.

Comparison between analytical and 3D finite element solutions for borehole stresses in anisotropic elastic rock

We present a rigorous validation of the analytical Amadei solution for the stress concentration around an arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients b11 and b55 are not equal. It is shown from theoretical considerations and published experimental data that the b11 and b55 are not equal for realistic rocks. Second, we develop a 3D finite element elastic model within a hybrid analytical–numerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic, transverse isotropic and orthorhombic symmetries. It is concluded that the analytical Amadei solution is valid with no restriction on the borehole orientation or the symmetry of the elastic anisotropy.

Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

A linear-elastic model of anisotropic tumour growth

This paper examines the effect of anisotropic growth on the evolution of mechanical stresses in a linear-elastic model of a growing, avascular tumour. This represents an important improvement on previous linear-elastic models of tissue growth since it has been shown recently that spatially-varying isotropic growth of linear-elastic tissues does not afford the necessary stress-relaxation for a steady-state stress distribution upon reaching a nutrient-regulated equilibrium size. Time-dependent numerical solutions are developed using a Lax-Wendroff scheme, which show the evolution of the tissue stress distributions over a period of growth until a steady-state is reached. These results are compared with the steady-state solutions predicted by the model equations, and key parameters influencing these steady-state distributions are identified. Recommendations for further extensions and applications of this model are proposed.

Revisiting borehole stresses in anisotropic elastic media : comparison of analytical versus numerical solutions

We present a rigorous validation of the analyticalAmadei solution for the stress concentration around arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients β11 and β55 are not equal. It is shown from theoretical considerations and published experimental data that the β11 and β55 are not equal for realistic rocks. Second, we develop a 3D finite-element elastic model within a hybrid analyticalnumerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic and transverse isotropic symmetries. It is concluded that the analytical Amadei solution is valid with no restrictions on the borehole orientation or elastic anisotropy symmetry.

Numerical investigation of case hardening of plant tissue during dying and its influence on the cellular level shrinkage

Dried plant food materials are one of the major contributors to the global food industry. Widening the fundamental understanding on different mechanisms of food material alterations during drying assists the development of novel dried food products and processing techniques. In this regard, case hardening is an important phenomenon, commonly observed during the drying processes of plant food materials, which significantly influences the product quality and process performance. In this work, a recent meshfree-based numerical model of the authors is further improved and used to simulate the influence of case hardening on shrinkage characteristics of plant tissues during drying. In order to model fluid and wall mechanisms in each cell, Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) are used. The model is fundamentally more capable of simulating large deformation of multiphase materials, when compared with conventional grid-based modelling techniques such as Finite Element Methods (FEM) or Finite Difference Methods (FDM). Case hardening is implemented by maintaining distinct moisture levels in the different cell layers of a given tissue. In order to compare and investigate different factors influencing tissue deformations under case hardening, four different plant tissue varieties (apple, potato, carrot and grape) are studied. The simulation results indicate that the inner cells of any given tissue undergo limited shrinkage and cell wall wrinkling compared to the case hardened outer cell layers of the tissues. When comparing unique deformation characteristics of the different tissues, irrespective of the normalised moisture content, the cell size, cell fluid turgor pressure and cell wall characteristics influence the tissue response to case hardening.

An experimental study on anisotropic mechanical properties of porcine thoracic aorta

Objective: To study the anisotropic mechanical properties of the thoracic aorta in porcine. Methods: Twenty-one porcine thoracic aortas were collected and categorized into three groups. The aortas were then cut through in their axial directions and expanded into two-dimensional planes. Then, by setting the length direction of the planar aortas (i.e., axial directions of the aortas) as 0°, each planar aorta was counterclockwisely cut into 8 samples with orientation of 30°, 45°, 60°, 90°, 120°, 135°, 150° and 180°, respectively. Finally, the uniaxial tensile tests were applied on three groups of samples at the loading rates of 1, 5 and 10 mm/min, respectively, to obtain the elastic modulus and ultimate stress of the aorta in different directions and at different loading rates. Results: The stress-strain curves exhibited different viscoelastic behaviors. With the increase of sample orientations, the elastic modulus gradually increased from 30°, reached the maximum value at 90°, and then gradually decreased till 180°. The variation trend of ultimate stress was similar to that of elastic modulus. Moreover, different loading rates showed a significant influence on the results of elastic modulus and ultimate stress, but a weak influence on the anisotropic degree. Conclusions: The porcine thoracic aorta is highly anisotropic. This research finding provides parameter references for assignment of material properties in finite element modeling, and is significant for understanding biomechanical properties of the arteries.