Comparison of three-dimensional profiles over time

In this paper, we describe an analysis for data collected on a three-dimensional spatial lattice with treatments applied at the horizontal lattice points. Spatial correlation is accounted for using a conditional autoregressive model. Observations are defined as neighbours only if they are at the same depth. This allows the corresponding variance components to vary by depth. We use the Markov chain Monte Carlo method with block updating, together with Krylov subspace methods, for efficient estimation of the model. The method is applicable to both regular and irregular horizontal lattices and hence to data collected at any set of horizontal sites for a set of depths or heights, for example, water column or soil profile data. The model for the three-dimensional data is applied to agricultural trial data for five separate days taken roughly six months apart in order to determine possible relationships over time. The purpose of the trial is to determine a form of cropping that leads to less moist soils in the root zone and beyond.We estimate moisture for each date, depth and treatment accounting for spatial correlation and determine relationships of these and other parameters over time.

Investigating the change in three dimensional deformity for idiopathic scoliosis using axially loaded MRI

Background: Adolescent idiopathic scoliosis is a complex three-dimensional deformity, involving a lateral deformity in the coronal plane and axial rotation of the vertebrae in the transverse plane. Gravitational loading plays an important biomechanical role in governing the coronal deformity, however, less is known about how they influence the axial deformity. This study investigates the change in three-dimensional deformity of a series of scoliosis patients due to compressive axial loading. Methods: Magnetic resonance imaging scans were obtained and coronal deformity (measured using the coronal Cobb angle) and axial rotations measured for a group of 18 scoliosis patients (Mean major Cobb angle was 43.4 o). Each patient was scanned in an unloaded and loaded condition while compressive loads equivalent to 50% body mass were applied using a custom developed compressive device. Findings: The mean increase in major Cobb angle due to compressive loading was 7.4 o (SD 3.5 o). The most axially rotated vertebra was observed at the apex of the structural curve and the largest average intravertebral rotations were observed toward the limits of the coronal deformity. A level-wise comparison showed no significant difference between the average loaded and unloaded vertebral axial rotations (intra-observer error = 2.56 o) or intravertebral rotations at each spinal level. Interpretation: This study suggests that the biomechanical effects of axial loading primarily influence the coronal deformity, with no significant change in vertebral axial rotation or intravertebral rotation observed between the unloaded and loaded condition. However, the magnitude of changes in vertebral rotation with compressive loading may have been too small to detect given the resolution of the current technique.

The Galerkin finite element approximation of the fractional cable equation

The cable equation is one of the most fundamental equations for modeling neuronal dynamics. Cable equations with a fractional order temporal derivative have been introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, the fractional cable equation involving two integro-differential operators is considered. The Galerkin finite element approximations of the fractional cable equation are proposed. The main contribution of this work is outlined as follow: • A semi-discrete finite difference approximation in time is proposed. We prove that the scheme is unconditionally stable, and the numerical solution converges to the exact solution with order O(Δt). • A semi-discrete difference scheme for improving the order of convergence for solving the fractional cable equation is proposed, and the numerical solution converges to the exact solution with order O((Δt)2). • Based on the above semi-discrete difference approximations, Galerkin finite element approximations in space for a full discretization are also investigated. • Finally, some numerical results are given to demonstrate the theoretical analysis.

Mixed-dimensional consistent coupling by multi-point constraint equations for efficient multi-scale modeling

Finding an appropriate linking method to connect different dimensional element types in a single finite element model is a key issue in the multi-scale modeling. This paper presents a mixed dimensional coupling method using multi-point constraint equations derived by equating the work done on either side of interface connecting beam elements and shell elements for constructing a finite element multiscale model. A typical steel truss frame structure is selected as case example and the reduced scale specimen of this truss section is then studied in the laboratory to measure its dynamic and static behavior in global truss and local welded details while the different analytical models are developed for numerical simulation. Comparison of dynamic and static response of the calculated results among different numerical models as well as the good agreement with those from experimental results indicates that the proposed multi-scale model is efficient and accurate.

An analytical tool that quantifies cellular morphology changes from three-dimensional fluorescence images

The most common software analysis tools available for measuring fluorescence images are for two-dimensional (2D) data that rely on manual settings for inclusion and exclusion of data points, and computer-aided pattern recognition to support the interpretation and findings of the analysis. It has become increasingly important to be able to measure fluorescence images constructed from three-dimensional (3D) datasets in order to be able to capture the complexity of cellular dynamics and understand the basis of cellular plasticity within biological systems. Sophisticated microscopy instruments have permitted the visualization of 3D fluorescence images through the acquisition of multispectral fluorescence images and powerful analytical software that reconstructs the images from confocal stacks that then provide a 3D representation of the collected 2D images. Advanced design-based stereology methods have progressed from the approximation and assumptions of the original model-based stereology(1) even in complex tissue sections(2). Despite these scientific advances in microscopy, a need remains for an automated analytic method that fully exploits the intrinsic 3D data to allow for the analysis and quantification of the complex changes in cell morphology, protein localization and receptor trafficking. Current techniques available to quantify fluorescence images include Meta-Morph (Molecular Devices, Sunnyvale, CA) and Image J (NIH) which provide manual analysis. Imaris (Andor Technology, Belfast, Northern Ireland) software provides the feature MeasurementPro, which allows the manual creation of measurement points that can be placed in a volume image or drawn on a series of 2D slices to create a 3D object. This method is useful for single-click point measurements to measure a line distance between two objects or to create a polygon that encloses a region of interest, but it is difficult to apply to complex cellular network structures. Filament Tracer (Andor) allows automatic detection of the 3D neuronal filament-like however, this module has been developed to measure defined structures such as neurons, which are comprised of dendrites, axons and spines (tree-like structure). This module has been ingeniously utilized to make morphological measurements to non-neuronal cells(3), however, the output data provide information of an extended cellular network by using a software that depends on a defined cell shape rather than being an amorphous-shaped cellular model. To overcome the issue of analyzing amorphous-shaped cells and making the software more suitable to a biological application, Imaris developed Imaris Cell. This was a scientific project with the Eidgenössische Technische Hochschule, which has been developed to calculate the relationship between cells and organelles. While the software enables the detection of biological constraints, by forcing one nucleus per cell and using cell membranes to segment cells, it cannot be utilized to analyze fluorescence data that are not continuous because ideally it builds cell surface without void spaces. To our knowledge, at present no user-modifiable automated approach that provides morphometric information from 3D fluorescence images has been developed that achieves cellular spatial information of an undefined shape (Figure 1). We have developed an analytical platform using the Imaris core software module and Imaris XT interfaced to MATLAB (Mat Works, Inc.). These tools allow the 3D measurement of cells without a pre-defined shape and with inconsistent fluorescence network components. Furthermore, this method will allow researchers who have extended expertise in biological systems, but not familiarity to computer applications, to perform quantification of morphological changes in cell dynamics.

Pulse dynamics in a three-component system : existence analysis

In this article, we analyze the three-component reaction-diffusion system originally developed by Schenk et al. (PRL 78:3781–3784, 1997). The system consists of bistable activator-inhibitor equations with an additional inhibitor that diffuses more rapidly than the standard inhibitor (or recovery variable). It has been used by several authors as a prototype three-component system that generates rich pulse dynamics and interactions, and this richness is the main motivation for the analysis we present. We demonstrate the existence of stationary one-pulse and two-pulse solutions, and travelling one-pulse solutions, on the real line, and we determine the parameter regimes in which they exist. Also, for one-pulse solutions, we analyze various bifurcations, including the saddle-node bifurcation in which they are created, as well as the bifurcation from a stationary to a travelling pulse, which we show can be either subcritical or supercritical. For two-pulse solutions, we show that the third component is essential, since the reduced bistable two-component system does not support them. We also analyze the saddle-node bifurcation in which two-pulse solutions are created. The analytical method used to construct all of these pulse solutions is geometric singular perturbation theory, which allows us to show that these solutions lie in the transverse intersections of invariant manifolds in the phase space of the associated six-dimensional travelling wave system. Finally, as we illustrate with numerical simulations, these solutions form the backbone of the rich pulse dynamics this system exhibits, including pulse replication, pulse annihilation, breathing pulses, and pulse scattering, among others.

Three dimensional simulation of a parabolic trough concentrator thermal collector

Parabolic Trough Concentrators (PTC) are the most proven solar collectors for solar thermal power plants, and are suitable for concentrating photovoltaic (CPV) applications. PV cells are sensitive to spatial uniformity of incident light and the cell operating temperature. This requires the design of CPV-PTCs to be optimised both optically and thermally. Optical modelling can be performed using Monte Carlo Ray Tracing (MCRT), with conjugate heat transfer (CHT) modelling using the computational fluid dynamics (CFD) to analyse the overall designs. This paper develops and evaluates a CHT simulation for a concentrating solar thermal PTC collector. It uses the ray tracing work by Cheng et al. (2010) and thermal performance data for LS-2 parabolic trough used in the SEGS III-VII plants from Dudley et al. (1994). This is a preliminary step to developing models to compare heat transfer performances of faceted absorbers for concentrating photovoltaic (CPV) applications. Reasonable agreement between the simulation results and the experimental data confirms the reliability of the numerical model. The model explores different physical issues as well as computational issues for this particular kind of system modeling. The physical issues include the resultant non-uniformity of the boundary heat flux profile and the temperature profile around the tube, and uneven heating of the HTF. The numerical issues include, most importantly, the design of the computational domain/s, and the solution techniques of the turbulence quantities and the near-wall physics. This simulation confirmed that optical simulation and the computational CHT simulation of the collector can be accomplished independently.

Improved design rules for fixed ended cold-formed steel columns subject to flexural-torsional buckling

This paper has presented the details of an investigation into the flexural and flexuraltorsional buckling behaviour of cold-formed structural steel columns with pinned and fixed ends. Current design rules for the member capacities of cold-formed steel columns are based on the same non-dimensional strength curve for both fixed and pinned-ended columns. This research has reviewed the accuracy of the current design rules in AS/NZS 4600 and the North American Specification in determining the member capacities of cold-formed steel columns using the results from detailed finite element analyses and an experimental study of lipped channel columns. It was found that the current Australian and American design rules accurately predicted the member capacities of pin ended lipped channel columns undergoing flexural and flexural torsional buckling. However, for fixed ended columns with warping fixity undergoing flexural-torsional buckling, it was found that the current design rules significantly underestimated the column capacities as they disregard the beneficial effect of warping fixity. This paper has therefore proposed improved design rules and verified their accuracy using finite element analysis and test results of cold-formed lipped channel columns made of three cross-sections and five different steel grades and thicknesses.

Ultimate flexural limit states analysis of prestressed concrete sleeper

Railway is one of the most important, reliable and widely used means of transportation, carrying freight, passengers, minerals, grains, etc. Thus, research on railway tracks is extremely important for the development of railway engineering and technologies. The safe operation of a railway track is based on the railway track structure that includes rails, fasteners, pads, sleepers, ballast, subballast and formation. Sleepers are very important components of the entire structure and may be made of timber, concrete, steel or synthetic materials. Concrete sleepers were first installed around the middle of last century and currently are installed in great numbers around the world. Consequently, the design of concrete sleepers has a direct impact on the safe operation of railways. The "permissible stress" method is currently most commonly used to design sleepers. However, the permissible stress principle does not consider the ultimate strength of materials, probabilities of actual loads, and the risks associated with failure, all of which could lead to the conclusion of cost-ineffectiveness and over design of current prestressed concrete sleepers. Recently the limit states design method, which appeared in the last century and has been already applied in the design of buildings, bridges, etc, is proposed as a better method for the design of prestressed concrete sleepers. The limit states design has significant advantages compared to the permissible stress design, such as the utilisation of the full strength of the member, and a rational analysis of the probabilities related to sleeper strength and applied loads. This research aims to apply the ultimate limit states design to the prestressed concrete sleeper, namely to obtain the load factors of both static and dynamic loads for the ultimate limit states design equations. However, the sleepers in rail tracks require different safety levels for different types of tracks, which mean the different types of tracks have different load factors of limit states design equations. Therefore, the core tasks of this research are to find the load factors of the static component and dynamic component of loads on track and the strength reduction factor of the sleeper bending strength for the ultimate limit states design equations for four main types of tracks, i.e., heavy haul, freight, medium speed passenger and high speed passenger tracks. To find those factors, the multiple samples of static loads, dynamic loads and their distributions are needed. In the four types of tracks, the heavy haul track has the measured data from Braeside Line (A heavy haul line in Central Queensland), and the distributions of both static and dynamic loads can be found from these data. The other three types of tracks have no measured data from sites and the experimental data are hardly available. In order to generate the data samples and obtain their distributions, the computer based simulations were employed and assumed the wheel-track impacts as induced by different sizes of wheel flats. A valid simulation package named DTrack was firstly employed to generate the dynamic loads for the freight and medium speed passenger tracks. However, DTrack is only valid for the tracks which carry low or medium speed vehicles. Therefore, a 3-D finite element (FE) model was then established for the wheel-track impact analysis of the high speed track. This FE model has been validated by comparing its simulation results with the DTrack simulation results, and with the results from traditional theoretical calculations based on the case of heavy haul track. Furthermore, the dynamic load data of the high speed track were obtained from the FE model and the distributions of both static and dynamic loads were extracted accordingly. All derived distributions of loads were fitted by appropriate functions. Through extrapolating those distributions, the important parameters of distributions for the static load induced sleeper bending moment and the extreme wheel-rail impact force induced sleeper dynamic bending moments and finally, the load factors, were obtained. Eventually, the load factors were obtained by the limit states design calibration based on reliability analyses with the derived distributions. After that, a sensitivity analysis was performed and the reliability of the achieved limit states design equations was confirmed. It has been found that the limit states design can be effectively applied to railway concrete sleepers. This research significantly contributes to railway engineering and the track safety area. It helps to decrease the failure and risks of track structure and accidents; better determines the load range for existing sleepers in track; better rates the strength of concrete sleepers to support bigger impact and loads on railway track; increases the reliability of the concrete sleepers and hugely saves investments on railway industries. Based on this research, many other bodies of research can be promoted in the future. Firstly, it has been found that the 3-D FE model is suitable for the study of track loadings and track structure vibrations. Secondly, the equations for serviceability and damageability limit states can be developed based on the concepts of limit states design equations of concrete sleepers obtained in this research, which are for the ultimate limit states.

Three-dimensional modeling of steel portal frame buildings

The realistic strength and deflection behavior of industrial and commercial steel portal frame buildings are understood only if the effects of rigidity of end frames and profiled steel claddings are included. The conventional designs ignore these effects and are very much based on idealized two-dimensional (2D) frame behavior. Full-scale tests of a 1212 m steel portal frame building under a range of design load cases indicated that the observed deflections and bending moments in the portal frame were considerably different from those obtained from a 2D analysis of frames ignoring these effects. Three-dimensional (3D) analyses of the same building, including the effects of end frames and cladding, were carried out, and the results agreed well with full-scale test results. Results clearly indicated the need for such an analysis and for testing to study the true behavior of steel portal frame buildings. It is expected that such a 3D analysis will lead to lighter steel frames as the maximum moments and deflections are reduced.

Sub-modelling for the ratchetting failure of insulated rail joints

Insulated rail joints are critical for train safety as they control electrical signalling systems; unfortunately they exhibit excessive ratchetting of the railhead near the endpost insulators. This paper reports a three-dimensional global model of these joints under wheel–rail contact pressure loading and a sub-model examining the ratchetting failures of the railhead. The sub-model employs a non-linear isotropic–kinematic elastic–plastic material model and predicts stress/strain levels in the localised railhead zone adjacent to the endpost which is placed in the air gap between the two rail ends at the insulated rail joint. The equivalent plastic strain plot is utilised to capture the progressive railhead damage adequately. Associated field and laboratory testing results of damage to the railhead material suggest that the simulation results are reasonable.

Deformation of a three-dimensional red blood cell in a stenosed micro-capillary

Red blood cells (RBCs) exhibit different types of motions and deformations when the blood flows through capillaries. Interestingly, due to the complex three-dimensional structure of the RBC membrane, RBCs show three-dimensional motions and deformations in the blood flow. These motions and deformations of the RBCs highly depend on the stiffness of the RBC membrane and on the geometrical parameters of the capillary through which blood flows. However, capillaries always do not have uniform cross sections and some capillaries have stenosed segments, where cross sectional area suddenly reduces. Further, some diseases can alter the stiffness of the RBC membrane drastically. In this study, the deformation behaviour of a single three-dimensional RBC is examined, when it moves through a stenosed capillary. A three-dimensional spring network is used to model the RBC membrane. The RBC’s inside and outside fluids are discretized into a finite number of mass points and treated by smoothed particle hydrodynamics (SPH) method. The capillary is considered as a rigid tube with a stenosed section. The deformation index, mean velocity and total energy of the RBC are analysed when it flows through the stenosed capillary. Further, motion and deformation of the RBCs with different membrane stiffness (KB) are compared when they flow through the stenosed segment of the capillary. The simulation results demonstrate the RBCs are subjected to a larger deformation when they move through the stenosed part of the capillary and the RBCs with lower KBvalues easily pass through the stenosed segment of the capillary. Further, RBCs having higher KBvalues have a lower mean velocity and it leads to slow down the overall blood flow rate

A three dimensional virtual medical imaging computed tomography suite: Innovation in education

Aims: The Medical Imaging Training Immersive Environment(MITIE) Computed Tomography(CT) system is an innovative virtual reality (VR) platform that allows students to practice a range of CT techniques. The aim of this pilot study was to harvest user feedback about the educational value of teh application and inform future pedagogical development. This presentation explores the use of this technology for skills training. Background: MITIE CT is a 3D VR environment that allows students to position a patient,and set CT technical parameters including IV contrast dose and dose rate. As with VR initiatives in other health disciplines the software mimics clinical practice as much as possible and uses 3D technology to enhance immersion and realism. The software is new and was developed by the Medical Imaging Course Team at a provider University with funding from a Health Workforce Australia 'Simulated Learning Environments' grant Methods: Current third year medical imaging students were provided with additional 1 hour MITIE laboratory tutorials and studnet feedback was collated with regard to educational value and performance. Ethical approval for the project was provided by the university ethics panel Results: This presentation provides qualitative analysis of student perceptions relating to satisfaction, usability and educational value. Students reported high levels of satisfaction and both feedback and assessment results confirmed the application's significance as a pre-clinical tool. There was a clear emerging theme that MITIE could be a useful learning tool that students could access to consolidate their clinical learning, either on campus or during their clinical placement. Conclusion: Student feedback indicates that MITIE CT has a valuable role to play in the clinial skills training for medical imaging students both in the academic and clinical environment. Future work will establish a framework for an appropriate supprting pedagogy that can cross the boundary between the two environments

Multifunctional three-dimensional T-junction graphene micro-wells: Energy-efficient, plasma-enabled growth and instant water-based transfer for flexible device applications

The “third-generation” 3D graphene structures, T-junction graphene micro-wells (T-GMWs) are produced on cheap polycrystalline Cu foils in a single-step, low-temperature (270 °C), energy-efficient, and environment-friendly dry plasma-enabled process. T-GMWs comprise vertical graphene (VG) petal-like sheets that seemlessly integrate with each other and the underlying horizontal graphene sheets by forming T-junctions. The microwells have the pico-to-femto-liter storage capacity and precipitate compartmentalized PBS crystals. The T-GMW films are transferred from the Cu substrates, without damage to the both, in de-ionized or tap water, at room temperature, and without commonly used sacrificial materials or hazardous chemicals. The Cu substrates are then re-used to produce similar-quality T-GMWs after a simple plasma conditioning. The isolated T-GMW films are transferred to diverse substrates and devices and show remarkable recovery of their electrical, optical, and hazardous NO2 gas sensing properties upon repeated bending (down to 1 mm radius) and release of flexible trasparent display plastic substrates. The plasma-enabled mechanism of T-GMW isolation in water is proposed and supported by the Cu plasma surface modification analysis. Our GMWs are suitable for various optoelectronic, sesning, energy, and biomedical applications while the growth approach is potentially scalable for future pilot-scale industrial production.

Universal crack closure integral for SIF estimation

Numerical analysis of cracked structures often involves numerical estimation of stress intensity factors (SIFs) at a crack tip/front. A newly developed formulation called universal crack closure integral (UCCI) for the evaluation of potential energy release rates (PERRs) and the corresponding SIFs is presented in this paper. Unlike the existing element dedicated forms of crack closure integrals (MCCI, VCCI) with application limited to finite element analysis, this new numerical SIF/PERR estimation technique is independent of the basic stress analysis procedure, making it universally applicable. The second merit of this procedure is that it avoids the generally error-producing zones close to the crack tip/front singularity. The UCCI procedure, based on Irwin's original CCI, is formulated and explored using a simple 2D problem of a straight crack in an infinite sheet. It is then applied to some three-dimensional crack geometries with the stresses and displacements obtained from a boundary element program.