Modeling the interaction of biological cells with a solidifying interface

In this article, we develop a modified level set method for modeling the interaction of particles with a solidifying interface. The dynamic computation of the van der Waals and drag forces between the particles and the solidification front leads to a problem of multiple length scales, which we resolve using adaptive grid techniques. We present a variety of example problems to demonstrate the accuracy and utility of the method. We also use the model to interpret experimental results obtained using directional solidification in a cryomicroscope.

[1]  K. Diller,et al.  Intracellular ice formation in glycerolized red cells , 1975 .

[2]  A. Karma,et al.  Quantitative phase-field model of alloy solidification. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

[4]  B. Rubinsky,et al.  Mechanical interactions between ice crystals and red blood cells during directional solidification. , 1994, Cryobiology.

[5]  F. Fonseca,et al.  The high viscosity encountered during freezing in glycerol solutions: effects on cryopreservation. , 2006, Cryobiology.

[6]  J. D. de Pablo,et al.  Molecular simulation study on the influence of dimethylsulfoxide on the structure of phospholipid bilayers. , 2003, Biophysical journal.

[7]  C. Körber,et al.  Encapsulation of human erythrocytes by growing ice crystals. , 1994, Cryobiology.

[8]  D. Stefanescu,et al.  An analytical model for the interaction between an insoluble particle and an advancing solid/liquid interface , 1992 .

[9]  K W Cole,et al.  Roles of unfrozen fraction, salt concentration, and changes in cell volume in the survival of frozen human erythrocytes. , 1988, Cryobiology.

[10]  C. Körber,et al.  Survival of directionally solidified B-lymphoblasts under various crystal growth conditions. , 1992, Cryobiology.

[12]  C. Körber,et al.  Phenomena at the advancing ice–liquid interface: solutes, particles and biological cells , 1988, Quarterly Reviews of Biophysics.

[13]  J. Walsh,et al.  Measurement of Cell Volume Loss in the Liquid Region Preceding an Advancing Phase Change Interface a , 1998, Annals of the New York Academy of Sciences.

[14]  Doru M. Stefanescu,et al.  A dynamic model for the interaction between a solid particle and an advancing solid/liquid interface , 1992 .

[15]  Mark S. Shephard,et al.  The versatility of automatic mesh generators based on tree structures and advanced geometric constructs , 1988 .

[16]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[17]  H. S. Udaykumar,et al.  Simulation of micro-scale interaction between ice and biological cells , 2003 .

[18]  Numerical simulation of the interaction of biological cells with an ice front during freezing , 2000 .

[19]  J. Baust.  Advances in Biopreservation , 2006 .

[20]  H. Udaykumar,et al.  Sharp interface Cartesian grid method III: Solidification of pure materials and binary solutions , 2005 .

[21]  K. W. Cole,et al.  Influence of cell concentration on the contribution of unfrozen fraction and salt concentration to the survival of slowly frozen human erythrocytes. , 1985, Cryobiology.

[22]  H. Udaykumar,et al.  Sharp interface Cartesian grid method I: An easily implemented technique for 3D moving boundary computations , 2005 .

[23]  Department of Physics,et al.  EFFICIENT COMPUTATION OF DENDRITIC MICROSTRUCTURES USING ADAPTIVE MESH REFINEMENT , 1998 .

[24]  Stephen H. Davis,et al.  Theory of Solidification , 2001 .

[25]  J. K. Spelt,et al.  Determination of the surface tension of biological cells using the freezing front technique , 1982, Cell Biophysics.

[26]  Britta Nestler,et al.  Phase-field model for solidification of a monotectic alloy with convection , 2000 .

[27]  Laxmikant V. Kale,et al.  Parallelization of a level set method for simulating dendritic growth , 2006, J. Parallel Distributed Comput..

[28]  P. Colella,et al.  An Adaptive Level Set Approach for Incompressible Two-Phase Flows , 1997 .

[29]  R. Trivedi,et al.  On the role of confinement on solidification in pure materials and binary alloys , 2006, cond-mat/0606388.

[30]  Dantzig,et al.  Computation of dendritic microstructures using a level set method , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[31]  Ronald Fedkiw,et al.  A Level Set Approach for the Numerical Simulation of Dendritic Growth , 2003, J. Sci. Comput..

[32]  D. Stefanescu,et al.  The influence of buoyant forces and volume fraction of particles on the particle pushing/entrapment transition during directional solidification of Al/SiC and Al/graphite composites , 1990 .

[33]  J. W. Garvin,et al.  Effect of a premelted film on the dynamics of particle–solidification front interactions , 2006 .

[34]  J J de Pablo,et al.  Stabilization and preservation of Lactobacillus acidophilus in saccharide matrices. , 2000, Cryobiology.

[35]  Grant,et al.  Directional solidification in two and three dimensions. , 1993, Physical review letters.

[36]  Jonathan A. Dantzig,et al.  An adaptive mesh refinement scheme for solidification problems , 1996 .

[37]  J. Crowe,et al.  From anhydrobiosis to freeze-drying of eukaryotic cells. , 2002, Comparative biochemistry and physiology. Part A, Molecular & integrative physiology.

[38]  P. Mazur.  Limits to life at low temperatures and at reduced water contents and water activities , 1980, Origins of life.

[39]  J. Pötschke,et al.  On the behaviour of foreign particles at an advancing solid-liquid interface , 1989 .

[40]  Peter Mazur,et al.  Kinetics of Water Loss from Cells at Subzero Temperatures and the Likelihood of Intracellular Freezing , 1963, The Journal of general physiology.

[41]  A. Hubel,et al.  Freezing Characteristics of Genetically Modified Lymphocytes for the Treatment of MPS II , 1999, Cell transplantation.

[42]  Nikolas Provatas,et al.  Seaweed to dendrite transition in directional solidification. , 2003, Physical review letters.

[43]  P. Mazur,et al.  Relative contributions of the fraction of unfrozen water and of salt concentration to the survival of slowly frozen human erythrocytes. , 1981, Biophysical journal.

[44]  V. Bronstein,et al.  Rejection and capture of cells by ice crystals on freezing aqueous solutions , 1981 .

[45]  Bruce T. Murray,et al.  Adaptive phase-field computations of dendritic crystal growth , 1997 .

[46]  M. Shannon,et al.  The effect of laser light propagation through a self-induced inhomogeneous process gas on temperature dependent laser-assisted chemical etching , 2003 .

[47]  Juan J de Pablo,et al.  Molecular simulation study of phospholipid bilayers and insights of the interactions with disaccharides. , 2003, Biophysical journal.

[48]  A. Karma,et al.  Quantitative phase-field modeling of dendritic growth in two and three dimensions , 1996 .

[49]  J. W. Garvin,et al.  Particle-solidification front dynamics using a fully coupled approach, part II: comparison of drag expressions , 2003 .

[50]  H. Udaykumar,et al.  Sharp interface Cartesian grid method II: A technique for simulating droplet interactions with surfaces of arbitrary shape , 2005 .

[51]  H. C. Hamaker.  The London—van der Waals attraction between spherical particles , 1937 .

[52]  P. Mazur,et al.  Contributions of unfrozen fraction and of salt concentration to the survival of slowly frozen human erythrocytes: influence of warming rate. , 1983, Cryobiology.

[53]  P. Mazur,et al.  Osmotic shrinkage as a factor in freezing injury in plant tissue cultures. , 1976, Plant physiology.

[54]  N. Goldenfeld,et al.  Phase field model for three-dimensional dendritic growth with fluid flow. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.