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.

Numerical simulation of railway track supporting system using finite-infinite and thin layer elements under impulsive loads

The present study deals with two dimensional, numerical simulation of railway track supporting system subjected to dynamic excitation force. Under plane strain condition, the coupled finite-infinite elements to represent the near and far field stress distribution and thin layer interface element was employed to model the interfacial behavior between sleepers and ballast. To account for the relative debonding, slipping and crushing that could take place in the contact area between the sleepers and ballast, modified Mohr-Coulomb criterion was adopted. Furthermore an attempt has been made to consider the elasto-plastic material non-linearity of the railway track supporting media by employing different constitutive models to represent steel, concrete and supporting materials. Based on the proposed physical and constitutive modeling a code has been developed for dynamic loads. The applicability of the developed F.E code has been demonstrated by analyzing a real railway supporting structure.

Specificity and context in post-exercise recovery: It is not a one-size-fits-all approach

The concept of specificity of exercise prescription and training is a longstanding and widely accepted foundation of the exercise sciences. Simply, the principle holds that training adaptations are achieved relative to the stimulus applied. That is, the manipulation of training variables (e.g. intensity or loading, mode, volume and frequency) directly influences the acute training stimulus, and so the long-term adaptive response (Young et al., 2001; Bird et al., 2005). Translating this concept to practice then recommends that exercise be prescribed specific to the desired outcomes, and the more closely this is achieved, the greater the performance gain is likely to be. However, the cardiovascular and metabolic adaptations traditionally associated with long, slow distance training types, similarly achieved using high-intensity training methods (for a review see Gibala et al., 2012), highlights understanding of underlying physiology as paramount for effective training program design. Various other factors including illness, sleep and psychology also impact on the training stimulus (Halson, 2014) and must be managed collectively with appropriate post-exercise recovery to continue performance improvements and reduce overtraining and injury risks (Kenttä and Hassmén, 1998).

Wake angle for surface gravity waves on a finite depth fluid

Linear water wave theory suggests that wave patterns caused by a steadily moving disturbance are contained within a wedge whose half-angle depends on the depth-based Froude number $F_H$. For the problem of flow past an axisymmetric pressure distribution in a finite-depth channel, we report on the apparent angle of the wake, which is the angle of maximum peaks. For moderately deep channels, the dependence of the apparent wake angle on the Froude number is very different to the wedge angle, and varies smoothly as $F_H$ passes through the critical value $F_H=1$. For shallow water, the two angles tend to follow each other more closely, which leads to very large apparent wake angles for certain regimes.

A laboratory assessment of the choice of vessel size under individual transferable quota regimes

This paper examines the effect of individual transferable quota regimes on technology choice, such as choice of vessel size, by using the laboratory experiment method. We find that even if vessel sizes change over time, the quota price can converge to the fundamental value conditioned on the vessels chosen. We also find that subjects choose their vessel type to maximise their profits based on the quota trading prices in the previous period. This result implies that the efficiency of quota markets in the beginning period is important because any inefficiency in quota markets may affect vessel sizes in ensuing periods. Moreover, we find that the initial allocations may significantly influence vessel sizes through two channels: first, a higher initial allocation to a subject increases the likelihood that the subject invests in a large-sized vessel; second, the quota price may be higher and more unstable under unequal allocation than under equal allocation; thus, whether the allocation is equal influences subjects' choice of vessel type. © 2014 Australian Agricultural and Resource Economics Society Inc.

Size-resolved particle distribution and gaseous concentrations by real-world road tunnel measurement

Measurements of aerosol particle number size distributions (15-700 nm), CO and NOx were performed in a bus tunnel, Australia. Daily mean particle size distributions of mixed diesel/CNG (Compressed Natural Gas) buses traffic flow were determined in 4 consecutive measurement days. EFs (Emission Factors) of Particle size distribution of diesel buses and CNG buses were obtained by MLR (Multiple Linear Regression) methods, particle distributions of diesel buses and CNG buses were observed as single accumulation mode and nuclei-mode separately. Particle size distributions of mixed traffic flow were decomposed by two log-normal fitting curves for each 30 minutes interval mean scans, all the mix fleet PSD emission can be well fitted by the summation of two log-normal distribution curves, and these were composed of nuclei mode curve and accumulation curve, which were affirmed as the CNG buses and diesel buses PN emission curves respectively. Finally, particle size distributions of diesel buses and CNG buses were quantified by statistical whisker-box charts. For log-normal particle size distribution of diesel buses, accumulation mode diameters were 74.5～87.5nm, geometric standard deviations were 1.89～1.98. As to log-normal particle size distribution of CNG buses, nuclei-mode diameters were 21～24 nm, geometric standard deviations were 1.27～1.31.

Field size consistency of nominally matched linacs

Given that there is increasing recognition of the effect that submillimetre changes in collimator position can have on radiotherapy beam dosimetry, this study aimed to evaluate the potential variability in small field collimation that may exist between otherwise matched linacs. Field sizes and field output factors were measured using radiochromic film and an electron diode, for jaw- and MLC-collimated fields produced by eight dosimetrically matched Varian iX linacs (Varian Medical Systems, Palo Alto, USA). This study used nominal sizes from 0.6×0.6 to 10×10 cm215 , for jaw-collimated fields,and from 1×1 to 10×10 cm216 , for MLC-collimated fields, delivered from a zero (head up, beam directed vertically downward) gantry angle. Differences between the field sizes measured for the eight linacs exceeded the uncertainty of the film measurements and the repositioning uncertainty of the jaws and MLCs on one linac. The dimensions of fields defined by MLC leaves were more consistent between linacs, while also differing more from their nominal values than fields defined by orthogonal jaws. The field output factors measured for the different linacs generally increased with increasing measured field size for the nominal 0.6×0.6 and 1×1 cm2 fields, and became consistent between linacs for nominal field sizes of 2×2 cm2 25 and larger. The inclusion in radiotherapy treatment planning system beam data of small field output factors acquired in fields collimated by jaws (rather than the more-reproducible MLCs), associated with either the nominal or the measured field sizes, should be viewed with caution. The size and reproducibility of the fields (especially the small fields) used to acquire treatment planning data should be investigated thoroughly as part of the linac or planning system commissioning process. Further investigation of these issues, using different linac models, collimation systems and beam orientations, is recommended.

Effect of carrier size on the dispersion of salmeterol xinafoate from interactive mixtures

The objective of this study was to determine the influence of lactose carrier size on drug dispersion of salmeterol xinafoate (SX) from interactive mixtures. SX dispersion was measured by using the fine particle fractions determined by a twin stage impinger attached to a Rotahaler1. The particle size of the lactose carrier in the SX interactive mixtures was varied using a range of commercial inhalation-grade lactoses. In addition, differing size fractions of individual lactose samples were achieved by dry sieving. The dispersion ofSXappeared to increase as the particle size of the lactose carrier decreased for the mixtures prepared from different particle size commercial samples of lactose and from different sieve fractions of the same lactose. Fine particles of lactose (<5 mm) associated with the lactose carrier were removed from the carrier surface by a wet decantation process to produce lactose samples with low but similar concentrations of fine lactose particles. The fine particle fractions of SX in mixtures prepared with the decanted lactose decreased significantly (analysis of variance, p<0.001) and the degree of dispersion became independent of the volume mean diameter of the carriers (analysis of variance, p<0.05). The dispersion behavior is therefore associated with the presence of fine adhered particles associated with the carriers and the inherent size of the carrier itself has little influence on dispersion.

Two-dimensional bow and stern flows in a finite-depth fluid with constant vorticity

This thesis is concerned with two-dimensional free surface flows past semi-infinite surface-piercing bodies in a fluid of finite-depth. Throughout the study, it is assumed that the fluid in question is incompressible, and that the effects of viscosity and surface tension are negligible. The problems considered are physically important, since they can be used to model the flow of water near the bow or stern of a wide, blunt ship. Alternatively, the solutions can be interpreted as describing the flow into, or out of, a horizontal slot. In the past, all research conducted on this topic has been dedicated to the situation where the flow is irrotational. The results from such studies are extended here, by allowing the fluid to have constant vorticity throughout the flow domain. In addition, new results for irrotational flow are also presented. When studying the flow of a fluid past a surface-piercing body, it is important to stipulate in advance the nature of the free surface as it intersects the body. Three different possibilities are considered in this thesis. In the first of these possibilities, it is assumed that the free surface rises up and meets the body at a stagnation point. For this configuration, the nonlinear problem is solved numerically with the use of a boundary integral method in the physical plane. Here the semi-infinite body is assumed to be rectangular in shape, with a rounded corner. Supercritical solutions which satisfy the radiation condition are found for various values of the Froude number and the dimensionless vorticity. Subcritical solutions are also found; however these solutions violate the radiation condition and are characterised by a train of waves upstream. In the limit that the height of the body above the horizontal bottom vanishes, the flow approaches that due to a submerged line sink in a $90^\circ$ corner. This limiting problem is also examined as a special case. The second configuration considered in this thesis involves the free surface attaching smoothly to the front face of the rectangular shaped body. For this configuration, nonlinear solutions are computed using a similar numerical scheme to that used in the stagnant attachment case. It is found that these solution exist for all supercritical Froude numbers. The related problem of the cusp-like flow due to a submerged sink in a corner is also considered. Finally, the flow of a fluid emerging from beneath a semi-infinite flat plate is examined. Here the free surface is assumed to detach from the trailing edge of the plate horizontally. A linear problem is formulated under the assumption that the elevation of the plate is close to the undisturbed free surface level. This problem is solved exactly using the Wiener-Hopf technique, and subcritical solutions are found which are characterised by a train of sinusoidal waves in the far field. The nonlinear problem is also considered. Exact relations between certain parameters for supercritical flow are derived using conservation of mass and momentum arguments, and these are confirmed numerically. Nonlinear subcritical solutions are computed, and the results are compared to those predicted by the linear theory.

Influence of kaolinite particle size on cross-link density, microstructure and mechanical properties of latex blending styrene butadiene rubber composites

Cross-link density, microstructure and mechanical properties of styrene butadiene rubber (SBR) composites filled with different particle sized kaolinites are investigated. With the increase of kaolinite particle size, the cross-link density of the filled SBR composites, the dispersibility and orientation degree of kaolinite particles gradually decrease. Some big cracks in filled rubber composites are distributed along the fringe of kaolinite aggregates, and the absorbance of all the absorption bands of kaolinites gradually increase with the increase of kaolinite particle size. All mechanical property indexes of kaolinite filled SBR composites decrease due to the decrease of cross-linking and reduction of interface interaction between filler and rubber matrix.

A finite element modelling method for thin layer mortared masonry systems

This paper deals with a finite element modelling method for thin layer mortared masonry systems. In this method, the mortar layers including the interfaces are represented using a zero thickness interface element and the masonry units are modelled using an elasto-plastic, damaging solid element. The interface element is formulated using two regimes; i) shear-tension and ii) shearcompression. In the shear-tension regime, the failure of joint is consiedered through an eliptical failure criteria and in shear-compression it is considered through Mohr Coulomb type failure criterion. An explicit integration scheme is used in an implicit finite element framework for the formulation of the interface element. The model is calibrated with an experimental dataset from thin layer mortared masonry prism subjected to uniaxial compression, a triplet subjected to shear loads a beam subjected to flexural loads and used to predict the response of thin layer mortared masonry wallettes under orthotropic loading. The model is found to simulate the behaviour of a thin layer mortated masonry shear wall tested under pre-compression and inplane shear quite adequately. The model is shown to reproduce the failure of masonry panels under uniform biaxial state of stresses.

Understanding the relationship between interfragmentary movement and callus growth in fracture healing: A novel computational approach

Introduction & Aims Optimising fracture treatments requires a sound understanding of relationships between stability, callus development and healing outcomes. This has been the goal of computational modelling, but discrepancies remain between simulations and experimental results. We compared healing patterns vs fixation stiffness between a novel computational callus growth model and corresponding experimental data. Hypothesis We hypothesised that callus growth is stimulated by diffusible signals, whose production is in turn regulated by mechanical conditions at the fracture site. We proposed that introducing this scheme into computational models would better replicate the observed tissue patterns and the inverse relationship between callus size and fixation stiffness. Method Finite element models of bone healing under stiff and flexible fixation were constructed, based on the parameters of a parallel rat femoral osteotomy study. An iterative procedure was implemented, to simulate the development of callus and its mechanical regulation. Tissue changes were regulated according to published mechano-biological criteria. Predictions of healing patterns were compared between standard models, with a pre-defined domain for callus development, and a novel approach, in which periosteal callus growth is driven by a diffusible signal. Production of this signal was driven by local mechanical conditions. Finally, each model’s predictions were compared to the corresponding histological data. Results Models in which healing progressed within a prescribed callus domain predicted that greater interfragmentary movements would displace early periosteal bone formation further from the fracture. This results from artificially large distortional strains predicted near the fracture edge. While experiments showed increased hard callus size under flexible fixation, this was not reflected in the standard models. Allowing the callus to grow from a thin soft tissue layer, in response to a mechanically stimulated diffusible signal, results in a callus shape and tissue distribution closer to those observed histologically. Importantly, the callus volume increased with increasing interfragmentary movement. Conclusions A novel method to incorporate callus growth into computational models of fracture healing allowed us to successfully capture the relationship between callus size and fixation stability observed in our rat experiments. This approach expands our toolkit for understanding the influence of different fixation strategies on healing outcomes.

Evaluation of a patient-specific finite-element model to simulate conservative treatment in adolescent idiopathic scoliosis

Study design Retrospective validation study. Objectives To propose a method to evaluate, from a clinical standpoint, the ability of a finite-element model (FEM) of the trunk to simulate orthotic correction of spinal deformity and to apply it to validate a previously described FEM. Summary of background data Several FEMs of the scoliotic spine have been described in the literature. These models can prove useful in understanding the mechanisms of scoliosis progression and in optimizing its treatment, but their validation has often been lacking or incomplete. Methods Three-dimensional (3D) geometries of 10 patients before and during conservative treatment were reconstructed from biplanar radiographs. The effect of bracing was simulated by modeling displacements induced by the brace pads. Simulated clinical indices (Cobb angle, T1–T12 and T4–T12 kyphosis, L1–L5 lordosis, apical vertebral rotation, torsion, rib hump) and vertebral orientations and positions were compared to those measured in the patients' 3D geometries. Results Errors in clinical indices were of the same order of magnitude as the uncertainties due to 3D reconstruction; for instance, Cobb angle was simulated with a root mean square error of 5.7°, and rib hump error was 5.6°. Vertebral orientation was simulated with a root mean square error of 4.8° and vertebral position with an error of 2.5 mm. Conclusions The methodology proposed here allowed in-depth evaluation of subject-specific simulations, confirming that FEMs of the trunk have the potential to accurately simulate brace action. These promising results provide a basis for ongoing 3D model development, toward the design of more efficient orthoses.

Genetics of anisotropy asymmetry: Registration and sample size effects

Brain asymmetry has been a topic of interest for neuroscientists for many years. The advent of diffusion tensor imaging (DTI) allows researchers to extend the study of asymmetry to a microscopic scale by examining fiber integrity differences across hemispheres rather than the macroscopic differences in shape or structure volumes. Even so, the power to detect these microarchitectural differences depends on the sample size and how the brain images are registered and how many subjects are studied. We fluidly registered 4 Tesla DTI scans from 180 healthy adult twins (45 identical and fraternal pairs) to a geometrically-centered population mean template. We computed voxelwise maps of significant asymmetries (left/right hemisphere differences) for common fiber anisotropy indices (FA, GA). Quantitative genetic models revealed that 47-62% of the variance in asymmetry was due to genetic differences in the population. We studied how these heritability estimates varied with the type of registration target (T1- or T2-weighted) and with sample size. All methods consistently found that genetic factors strongly determined the lateralization of fiber anisotropy, facilitating the quest for specific genes that might influence brain asymmetry and fiber integrity.