Simulation of plant cell shrinkage during drying : A SPH–DEM approach

Plant based dried food products are popular commodities in global market where much research is focused to improve the products and processing techniques. In this regard, numerical modelling is highly applicable and in this work, a coupled meshfree particle-based two-dimensional (2-D) model was developed to simulate micro-scale deformations of plant cells during drying. Smoothed Particle Hydrodynamics (SPH) was used to model the viscous cell protoplasm (cell fluid) by approximating it to an incompressible Newtonian fluid. The visco-elastic characteristic of the cell wall was approximated to a Neo-Hookean solid material augmented with a viscous term and modelled with a Discrete Element Method (DEM). Compared to a previous work [H. C. P. Karunasena, W. Senadeera, Y. T. Gu and R. J. Brown, Appl. Math. Model., 2014], this study proposes three model improvements: linearly decreasing positive cell turgor pressure during drying, cell wall contraction forces and cell wall drying. The improvements made the model more comparable with experimental findings on dried cell morphology and geometric properties such as cell area, diameter, perimeter, roundness, elongation and compactness. This single cell model could be used as a building block for advanced tissue models which are highly applicable for product and process optimizations in Food Engineering.

A coupled SPH-DEM model for micro-scale structural deformations of plant cells during drying

A single plant cell was modeled with smoothed particle hydrodynamics (SPH) and a discrete element method (DEM) to study the basic micromechanics that govern the cellular structural deformations during drying. This two-dimensional particle-based model consists of two components: a cell fluid model and a cell wall model. The cell fluid was approximated to a highly viscous Newtonian fluid and modeled with SPH. The cell wall was treated as a stiff semi-permeable solid membrane with visco-elastic properties and modeled as a neo-Hookean solid material using a DEM. Compared to existing meshfree particle-based plant cell models, we have specifically introduced cell wall–fluid attraction forces and cell wall bending stiffness effects to address the critical shrinkage characteristics of the plant cells during drying. Also, a moisture domain-based novel approach was used to simulate drying mechanisms within the particle scheme. The model performance was found to be mainly influenced by the particle resolution, initial gap between the outermost fluid particles and wall particles and number of particles in the SPH influence domain. A higher order smoothing kernel was used with adaptive smoothing length to improve the stability and accuracy of the model. Cell deformations at different states of cell dryness were qualitatively and quantitatively compared with microscopic experimental findings on apple cells and a fairly good agreement was observed with some exceptions. The wall–fluid attraction forces and cell wall bending stiffness were found to be significantly improving the model predictions. A detailed sensitivity analysis was also done to further investigate the influence of wall–fluid attraction forces, cell wall bending stiffness, cell wall stiffness and the particle resolution. This novel meshfree based modeling approach is highly applicable for cellular level deformation studies of plant food materials during drying, which characterize large deformations.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Secondary flow induced by the rotation of two concentric spheres about their common diameter in micropolar fluids

In this paper, we have studied the secondary flow induced in a micropolar fluid by the rotation of two concentric spheres about a fixed diameter. The secondary flow exhibits behaviour commonly observed in visco-elastic fluids. In particular we have obtained the expressions for microrotation vector. Numerical results have been obtained for a number of values of relative rotations of the two spheres for a chosen set of values of fluid parameters. The results are presented graphically and compared with the previous investigations.

Magnetically actuated artificial cilia: the effect of fluid inertia.

Natural cilia are hairlike microtubule-based structures that are able to move fluid on the micrometer scale using asymmetric motion. In this article, we follow a biomimetic approach to design artificial cilia lining the inner surfaces of microfluidic channels with the goal of propelling fluid. The artificial cilia consist of polymer films filled with superparamagnetic nanoparticles, which can mimic the motion of natural cilia when subjected to a rotating magnetic field. To obtain the magnetic field and associated magnetization local to the cilia, we solve the Maxwell equations, from which the magnetic body moments and forces can be deduced. To obtain the ciliary motion, we solve the dynamic equations of motion, which are then fully coupled to the Navier-Stokes equations that describe the fluid flow around the cilia, thus taking full account of fluid inertial forces. The dimensionless parameters that govern the deformation behavior of the cilia and the associated fluid flow are arrived at using the principle of virtual work. The physical response of the cilia and the fluid flow for different combinations of elastic, fluid viscous, and inertia forces are identified.

Seismic mitigation of building structural systems using passive dampers

This paper treats the seismic mitigation of medium rise frame-shear wall structures and building facade systems using passive damping devices. The frame shear wall structures have embedded viscoelastic and friction dampers in different configurations and placed in various locations in the structure. Influence of damper type, configuration and location are investigated. Results for tip deflections which provide an overall evaluation of the seismic response of the structure, are determined. Seismic mitigation of building facade systems in which visco-elastic dampers are fitted at the horizontal connections between the facades and the frame, instead of the traditional rigid connections, are also treated. Finite element techniques are used to model and analyse the two structural systems under different earthquake loadings, scaled to the same peak ground acceleration for meaningful comparison of responses. Results demonstrate the feasibility of these techniques for seismic mitigation.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Biomimetic molecular organization of naphthalene diimide in the solid state: tunable (chiro-) optical, viscoelastic and nanoscale properties

The interfacing of aromatic molecules with biomolecules to design functional molecular materials is a promising area of research. Intermolecular interactions determine the performance of these materials and therefore, precise control over the molecular organization is necessary to improve functional properties. Herein we describe the tunable biomimetic molecular engineering of a promising n-type organic semiconductor, naphthalene diimide (NDI), in the solid state by introducing minute structural mutations in the form of amino acids with variable Ca-functionality. For the first time we could achieve all four possible crystal packing modes, namely cofacial, brickwork, herringbone and slipped stacks of the NDI system. Furthermore, amino acid conjugated NDIs exhibit ultrasonication induced organogels with tunable visco-elastic and temperature responsive emission properties. The amino acid-NDI conjugates self-assemble into 0D nanospheres and 1D nanofibers in their gel state while the ethylamine-NDI conjugate forms 2D sheets from its solution. Photophysical studies indicated the remarkable influence of molecular ordering on the absorption and fluorescence properties of NDIs. Interestingly, the circular dichroism (CD) and X-ray diffraction (XRD) studies revealed the existence of helical ordering of NDIs in both solution and solid state. The chiral amino acids and their conformations with respect to the central NDI core are found to influence the nature of the helical organization of NDIs. Consequently, the origin of the preferential handedness in the helical organization is attributed to transcription of chiral information from the amino acid to the NDI core. On account of these unique properties, the materials derived from NDI-conjugates might find a wide range of future interdisciplinary applications from materials to biomedicine.

Dynamics of intermittent force fluctuations in vesicular nanotubulation

Irregular force fluctuations are seen in most nanotubulation experiments. The dynamics behind their presence has, however, been neither commented upon nor modeled. A simple estimate of the mean energy dissipated in force drops turns out to be several times the thermal energy. This coupled with the rate dependent nature of the deformation reported in several experiments point to a dynamical origin of the serrations. We simplify the whole process of tether formation through a three-stage model of successive deformations of sphere to ellipsoid, neck-formation, and tubule birth and extension. Based on this, we envisage a rate-softening frictional force at the neck that must be overcome before a nanotube can be pulled out. Our minimal model includes elastic and visco-elastic deformation of the vesicle, and has built-in dependence on pull velocity, vesicle radius, and other material parameters, enabling us to capture various kinds of serrated force-extension curves for different parameter choices. Serrations are predicted in the nanotubulation region. Other features of force-extension plots reported in the literature such as a plateauing serrated region beyond a force drop, serrated flow region with a small positive slope, an increase in the elastic threshold with pull velocity, force-extension curves for vesicles with larger radius lying lower than those for smaller radius, are all also predicted by the model. A toy model is introduced to demonstrate that the role of the friction law is limited to inducing stick-slip oscillations in the force, and all other qualitative and quantitative features emerging from the model can only be attributed to other physical mechanisms included in the deformation dynamics of the vesicle. (C) 2014 AIP Publishing LLC.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Insight into differences in nanoindentation properties of bone.

Nanoindentation provides the ideal framework to determine mechanical properties of bone at the tissue scale without being affected by the size, shape, and porosity of the bone. However, the values of tissue level mechanical properties vary significantly between studies. Since the differences in the bone sample, hydration state, and test parameters complicate direct comparisons across the various studies, these discrepancies in values cannot be compared directly. The objective of the current study is to evaluate and compare mechanical properties of the same bones using a broad range of testing parameters. Wild type C56BL6 mice tibiae were embedded following different processes and tested in dry and rehydrated conditions. Spherical and Berkovich indenter probes were used, and data analysis was considered within the elasto-plastic (Oliver-Pharr), viscoelastic and visco-elastic-plastic frameworks. The mean values of plane strain modulus varied significantly depending on the hydration state, probe geometry and analysis method. Indentations in dry bone analyzed using a visco-elastic-plastic approach gave values of 34 GPa. After rehydrating the same bones and indenting them with a spherical tip and utilizing a viscoelastic analysis, the mean modulus value was 4 GPa, nearly an order of magnitude smaller. Results suggest that the hydration state, probe geometry and the limitations and assumptions of each analysis method influence significantly the measured mechanical properties. This is the first time that such a systematic study has been carried out and it has been concluded that the discrepancies in the mechanical properties of bone measured by nanoindentation found in the literature should not be attributed only to the differences between the bones themselves, but also to the testing and analysis protocols.

Rheological properties in blends of poly(aryl ether ether ketone) and liquid crystalline poly(aryl ether ketone)

Rheological properties of the blends of poly(aryl ether ether ketone) (PEEK) with liquid crystalline poly(aryl ether ketone) containing substituted 3-trifluoro-methylbenzene side group (F-PAEK), prepared by solution precipitation, have been investigated by rheometer. Dynamic rheological behaviors of the blends under the oscillatory shear mode are strongly dependent on blend composition. For PEEK-rich blends, the systems show flow curves similar to those of the pure PEEK, i.e., dynamic storage modulus G' is larger than dynamic loss modulus G", showing the feature of elastic fluid. For F-PAEK-rich systems, the rheological behavior of the blends has a resemblance to pure F-PAEK, i.e., G" is greater than G', showing the characteristic of viscous fluid. When the PEEK content is in the range of 50-70%, the blends exhibit an unusual rheological behavior, which is the result of phase inversion between the two components. Moreover, as a whole, the complex viscosity values of the blends are between those of two pure polymers and decrease with increasing F-PAEK content. However, at 50% weight fraction of PEEK, the viscosity-composition curves exhibit a local maximum, which may be mainly attributed to the phase separation of two components at such a composition.

Rheological characterisation and printing performance of Sn/Ag/Cu solder pastes

Lead-free solder paste printing process accounts for majority of the assembly defects in the electronic manufacturing industry. The study investigates rheological behaviour and stencil printing performance of the lead-free solder pastes (Sn/Ag/Cu). Oscillatory stress sweep test was carried out to study the visco-elastic behaviour of the lead-free solder pastes. The visco-elastic behaviour of the paste encompasses solid and liquid characteristic of the paste, which could be used to study the flow behaviour experienced by the pastes during the stencil printing process. From this study, it was found that the solid characteristics (G0) is higher than the liquid characteristic (G0 0) for the pastes material. In addition, the results from the study showed that the solder paste with a large G0 = G0 0 has a higher cohesiveness resulting in poor withdrawal of the paste during the stencil printing process. The phase angles (d) was used to correlate the quality of the dense suspensions to the formulation of solder paste materials. This study has revealed the value of having a rheological measurement for explaining and characterising solder pastes for stencil printing. As the demand for lead free pastes increases rheological measurements can assist with the formulation or development of new pastes.

Rheological characterisation of solder pastes and isotropic conductive adhesives used for flip-chip assembly

Solder pastes and isotropic conductive adhesives (ICAs) are widely used as a principal bonding medium in the electronic industry. This study investigates the rheological behaviour of the pastes (solder paste and isotropic conductive adhesives) used for flip-chip assembly. Oscillatory stress sweep test are performed to evaluate solid characteristic and cohesiveness of the lead-free solder pastes and isotropic conductive adhesive paste materials. The results show that the G' (storage modulus) is higher than G '' (loss modulus) for the pastes material indicating a solid like behaviour. It result shows that the linear visco-elastic region for the pastes lies in a very small stress range, below 10 Pa. in addition, the stress at which the value of storage modulus is equal to that of loss modulus can be used as an indicator of the paste cohesiveness. The measured cross-over stress at G'=G '' shows that the solder paste has higher stress at G'=G '' compared to conductive adhesives. Creep-recovery test method is used to study the slump behaviour in the paste materials. The conductive adhesive paste shows a good recovery when compared to the solder pastes. (C) 2008 Elsevier B.V. All rights reserved.

An investigation into the rheological properties of different lead-free solder pastes for surface mount applications

Purpose – The purpose of this paper is to investigate the rheological behaviour of three different lead-free solder pastes used for surface mount applications in the electronic industry.Design/methodology/approach – This study concerns the rheological measurements of solder paste samples and is made up of three parts. The first part deals with the measurement of rhelogical properties with three different measuring geometries, the second part looks into the effect of frequencies on oscillatory stress sweep measurements and the final part reports on the characterisation and comparison of three different types of Pb-free solder pastes. Findings – Among the three geometries, the serrated parallel plate was found effective in minimising the wall-slip effect. From the oscillatory stresssweep data with different frequencies; it was observed that the linear visco-elastic region is independent of frequency for all the solder paste samples. To understand the shear thinning behaviour of solder paste, the well known Cross and Carreau models were fitted to the viscosity data. Moreover,creep-recovery and dynamic frequency-sweep tests were also carried out without destroying the sample’s structure and have yielded useful information on the pastes behaviour.Research limitations/implications – More extensive research is needed to fully characterise the wall-slip behaviour during the rheological measurements of solder pastes. Practical implications – The rheological test results presented in this paper will be of important value for research and development, quality control and facilitation of the manufacturing of solder pastes and flux mediums. Originality/value – This paper shows how wall-slip effects can be effectively avoided during rheological measurements of solder pastes. The paper also outlines how different rheological test methods can be used to characterise solder paste behaviours