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.

A coiled-coil-repeat protein 'Ccrp' in Bdellovibrio bacteriovorus prevents cellular indentation, but is not essential for vibroid cell morphology

Bdellovibrio bacteriovorus are small, vibroid, predatory bacteria that grow within the periplasmic space of a host Gram-negative bacterium. The intermediate-filament (IF)-like protein crescentin is a member of a broad class of IF-like, coiled-coil-repeat-proteins (CCRPs), discovered in Caulobacter crescentus, where it contributes to the vibroid cell shape. The B. bacteriovorus genome has a single ccrp gene encoding a protein with an unusually long, stutter-free, coiled-coil prediction; the inactivation of this did not alter the vibriod cell shape, but caused cell deformations, visualized as chiselled insets or dents, near the cell poles and a general 'creased' appearance, under the negative staining preparation used for electron microscopy, but not in unstained, frozen, hydrated cells. Bdellovibrio bacteriovorus expressing 'teal' fluorescent protein (mTFP), as a C-terminal tag on the wild-type Ccrp protein, did not deform under negative staining, suggesting that the function was not impaired. Localization of fluorescent Ccrp-mTFP showed some bias to the cell poles, independent of the cytoskeleton, as demonstrated by the addition of the MreB-specific inhibitor A22. We suggest that the Ccrp protein in B. bacteriovorus contributes as an underlying scaffold, similar to that described for the CCRP protein FilP in Streptomyces coelicolor, preventing cellular indentation, but not contributing to the vibroid shape of the B. bacteriovorus cells.

A coupled SPH-DEM model for fluid and solid mechanics of apple parenchyma cells during drying

A coupled SPH-DEM based two-dimensional (2-D) micro-scale single cell model is developed to predict basic cell-level shrinkage effects of apple parenchyma cells during air drying. In this newly developed drying model, Smoothed Particle Hydrodynamics (SPH) is used to model the low Reynolds Number fluid motions of the cell protoplasm, and a Discrete Element Method (DEM) is employed to simulate the polymer-like cell wall. Simulations results reasonably agree with published experimental drying results on cellular shrinkage properties such as cellular area, diameter and perimeter. These preliminary results indicate that the model is effective for the modelling and simulation of apple parenchyma cells during air drying.

Scanning electron microscopic study of microstructure of Gala apples during hot air drying

Microscopic changes occur in plant food materials during drying significantly influence the macroscopic properties and quality factors of the dried food materials. It is very critical to study microstructure to understand the underlying cellular mechanisms to improve performance of the food drying techniques. However, there is very limited research conducted on such microstructural changes of plant food material during drying. In this work, Gala apple parenchyma tissue samples were studied using a scanning electron microscope for gradual microstructural changes as affected by temperature, time and moisture content during hot air drying at two drying temperatures: 57 ℃ and 70 ℃. For fresh samples, the average cellular parameter values were; cell area: 20000 μm2, ferret diameter: 160 μm, perimeter: 600 μm, roundness: 0.76, elongation: 1.45 and compactness: 0.84. During drying, a higher degree of cell shrinkage was observed with cell wall warping and increase in intercellular space. However, no significant cell wall breakage was observed. The overall reduction of cell area, ferret diameter and perimeter were about 60%, 40% and 30%. The cell roundness and elongation showed overall increments of about 5% and the compactness remained unchanged. Throughout the drying cycle, cellular deformations were mainly influenced by the moisture content. During the initial and intermediate stages of drying, cellular deformations were also positively influenced by the drying temperature and the effect was reversed at the final stages of drying which provides clues for case hardening of the material.

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 meshfree model for plant tissue deformations during drying

Plant tissue has a complex cellular structure which is an aggregate of individual cells bonded by middle lamella. During drying processes, plant tissue undergoes extreme deformations which are mainly driven by moisture removal and turgor loss. Numerical modelling of this problem becomes challenging when conventional grid-based modelling techniques such as Finite Element Methods (FEM) and Finite Difference Methods (FDM) have grid-based limitations. This work presents a meshfree approach to model and simulate the deformations of plant tissues during drying. This method demonstrates the fundamental capabilities of meshfree methods in handling extreme deformations of multiphase systems. A simplified 2D tissue model is developed by aggregating individual cells while accounting for the stiffness of the middle lamella. Each individual cell is simply treated as consisting of two main components: cell fluid and cell wall. The cell fluid is modelled using Smoothed Particle Hydrodynamics (SPH) and the cell wall is modelled using a Discrete Element Method (DEM). During drying, moisture removal is accounted for by reduction of cell fluid and wall mass, which causes local shrinkage of cells eventually leading to tissue scale shrinkage. The cellular deformations are quantified using several cellular geometrical parameters and a favourably good agreement is observed when compared to experiments on apple tissue. The model is also capable of visually replicating dry tissue structures. The proposed model can be used as a step in developing complex tissue models to simulate extreme deformations during drying.

Application of meshfree methods to numerically simulate microscale deformations of different plant food materials during drying

Plant food materials have a very high demand in the consumer market and therefore, improved food products and efficient processing techniques are concurrently being researched in food engineering. In this context, numerical modelling and simulation techniques have a very high potential to reveal fundamentals of the underlying mechanisms involved. However, numerical modelling of plant food materials during drying becomes quite challenging, mainly due to the complexity of the multiphase microstructure of the material, which undergoes excessive deformations during drying. In this regard, conventional grid-based modelling techniques have limited applicability due to their inflexible grid-based fundamental limitations. As a result, meshfree methods have recently been developed which offer a more adaptable approach to problem domains of this nature, due to their fundamental grid-free advantages. In this work, a recently developed meshfree based two-dimensional plant tissue model is used for a comparative study of microscale morphological changes of several food materials during drying. The model involves Smoothed Particle Hydrodynamics (SPH) and Discrete Element Method (DEM) to represent fluid and solid phases of the cellular structure. Simulation are conducted on apple, potato, carrot and grape tissues and the results are qualitatively and quantitatively compared and related with experimental findings obtained from the literature. The study revealed that cellular deformations are highly sensitive to cell dimensions, cell wall physical and mechanical properties, middle lamella properties and turgor pressure. In particular, the meshfree model is well capable of simulating critically dried tissues at lower moisture content and turgor pressure, which lead to cell wall wrinkling. The findings further highlighted the potential applicability of the meshfree approach to model large deformations of the plant tissue microstructure during drying, providing a distinct advantage over the state of the art grid-based approaches.

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 model for individual and collective cell movement in Dictyostelium discoideum

The cellular slime mold Dictyostelium discoideum is a widely used model system for studying a variety of basic processes in development, including cell–cell signaling, signal transduction, pattern formation, cell motility, and the movement of tissue-like aggregates of cells. Many aspects of cell motion are poorly understood, including how individual cell behavior produces the collective motion of cells observed within the mound and slug. Herein, we describe a biologically realistic model for motile D. discoideum cells that can generate active forces, that interact via surface molecules, and that can detect and respond to chemotactic signals. We model the cells as deformable viscoelastic ellipsoids and incorporate signal transduction and cell–cell signaling by using a previously developed model. The shape constraint restricts the admissible deformations but makes the simulation of a large number of interacting cells feasible. Because the model is based on known processes, the parameters can be estimated or measured experimentally. We show that this model can reproduce the observations on the chemotactic behavior of single cells, streaming during aggregation, and the collective motion of an aggregate of cells driven by a small group of pacemakers. The model predicts that the motion of two-dimensional slugs [Bonner, J. T. (1998) Proc. Natl. Acad. Sci. USA 95, 9355–9359] results from the same behaviors that are exhibited by individual cells; it is not necessary to invoke different mechanisms or behaviors. Our computational experiments also suggest previously uncharacterized phenomena that may be experimentally observable.

Effective Cell-Centred Time-Domain Maxwell's Equations Numerical Solvers

This research work analyses techniques for implementing a cell-centred finite-volume time-domain (ccFV-TD) computational methodology for the purpose of studying microwave heating. Various state-of-the-art spatial and temporal discretisation methods employed to solve Maxwell's equations on multidimensional structured grid networks are investigated, and the dispersive and dissipative errors inherent in those techniques examined. Both staggered and unstaggered grid approaches are considered. Upwind schemes using a Riemann solver and intensity vector splitting are studied and evaluated. Staggered and unstaggered Leapfrog and Runge-Kutta time integration methods are analysed in terms of phase and amplitude error to identify which method is the most accurate and efficient for simulating microwave heating processes. The implementation and migration of typical electromagnetic boundary conditions. from staggered in space to cell-centred approaches also is deliberated. In particular, an existing perfectly matched layer absorbing boundary methodology is adapted to formulate a new cell-centred boundary implementation for the ccFV-TD solvers. Finally for microwave heating purposes, a comparison of analytical and numerical results for standard case studies in rectangular waveguides allows the accuracy of the developed methods to be assessed.