From video to computation of biological fluid–structure interaction problems

This work deals with the techniques necessary to obtain a purely Eulerian procedure to conduct CFD simulations of biological systems with moving boundary flow phenomena. Eulerian approaches obviate difficulties associated with mesh generation to describe or fit flow meshes to body surfaces. The challenges associated with constructing embedded boundary information, body motions and applying boundary conditions on the moving bodies for flow computation are addressed in the work. The overall approach is applied to the study of a fluid–structure interaction problem, i.e., the hydrodynamics of swimming of an American eel, where the motion of the eel is derived from video imaging. It is shown that some first-blush approaches do not work, and therefore, careful consideration of appropriate techniques to connect moving images to flow simulations is necessary and forms the main contribution of the paper. A combination of level set-based active contour segmentation with optical flow and image morphing is shown to enable the image-to-computation process.

[1]  Kevin J. Parker,et al.  Multiple Resolution Bayesian Segmentation of Ultrasound Images , 1995 .

[2]  Andrew F. Laine,et al.  Speckle reduction and contrast enhancement of echocardiograms via multiscale nonlinear processing , 1998, IEEE Transactions on Medical Imaging.

[3]  D. Adkins,et al.  CFD simulation of fish-like body moving in viscous liquid , 2006 .

[4]  Srinivasan Rajagopalan,et al.  Optical flow based volumetric spatio-temporal interpolation , 2008, 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[5]  Berthold K. P. Horn,et al.  Determining Optical Flow , 1981, Other Conferences.

[6]  H. S. Udaykumar,et al.  Adaptively refined, parallelised sharp interface Cartesian grid method for three-dimensional moving boundary problems , 2009 .

[7]  H. S. Udaykumar,et al.  A massively parallel adaptive sharp interface solver with application to mechanical heart valve simulations , 2012 .

[8]  C. Alvino,et al.  Active contours and optical flow for automatic tracking of flying vehicles , 2004, Proceedings of the 2004 American Control Conference.

[9]  H. Liu Simulation-Based Biological Fluid Dynamics in Animal Locomotion , 2005 .

[10]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[11]  Ross T. Whitaker A level-set approach to image blending , 2000, IEEE Trans. Image Process..

[12]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[13]  Christopher V. Alvino,et al.  Multigrid computation of rotationally invariant non-linear optical flow , 2005, IEEE International Conference on Image Processing 2005.

[14]  Seth Dillard,et al.  Techniques to derive geometries for image-based Eulerian computations. , 2014, Engineering computations.

[15]  Thomas Brox,et al.  High Accuracy Optical Flow Estimation Based on a Theory for Warping , 2004, ECCV.

[16]  H. Udaykumar,et al.  A sharp interface method for high-speed multi-material flows: strong shocks and arbitrary materialpairs , 2011 .

[17]  Ulrich Rüde,et al.  3D optical flow computation using a parallel variational multigrid scheme with application to cardiac C-arm CT motion , 2007, Image Vis. Comput..

[18]  L. Costaridou,et al.  A wavelet-based spatially adaptive method for mammographic contrast enhancement. , 2003, Physics in medicine and biology.

[19]  Jun Ohya,et al.  Contrast enhancement of color images based on Wavelet Transform and human visual system , 2007 .

[20]  Liza Shrestha,et al.  From medical images to flow computations without user‐generated meshes , 2014, International journal for numerical methods in biomedical engineering.

[21]  Chunming Li,et al.  Distance Regularized Level Set Evolution and Its Application to Image Segmentation , 2010, IEEE Transactions on Image Processing.

[22]  Nahum Kiryati,et al.  Piecewise-Smooth Dense Optical Flow via Level Sets , 2006, International Journal of Computer Vision.

[23]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[24]  Scott T. Acton,et al.  Speckle reducing anisotropic diffusion , 2002, IEEE Trans. Image Process..

[25]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

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

[27]  Jia Lu,et al.  Fluid–structure interaction methods in biological flows with special emphasis on heart valve dynamics , 2010 .

[28]  Olivier Basset,et al.  Segmentation of ultrasound images--multiresolution 2D and 3D algorithm based on global and local statistics , 2003, Pattern Recognit. Lett..

[29]  John C. Platt,et al.  Elastically deformable models , 1987, SIGGRAPH.

[30]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[31]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[32]  P. Olver,et al.  Conformal curvature flows: From phase transitions to active vision , 1996, ICCV 1995.

[33]  J. Koseff,et al.  A non-staggered grid, fractional step method for time-dependent incompressible Navier-Stokes equations in curvilinear coordinates , 1994 .

[34]  Tony F. Chan,et al.  A Multiphase Level Set Framework for Image Segmentation Using the Mumford and Shah Model , 2002, International Journal of Computer Vision.

[35]  H. S. Udaykumar,et al.  Two-Dimensional Dynamic Simulation of Platelet Activation During Mechanical Heart Valve Closure , 2006, Annals of Biomedical Engineering.

[36]  A. Murat Tekalp,et al.  Motion segmentation by multistage affine classification , 1997, IEEE Trans. Image Process..

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

[38]  C. Lanczos The variational principles of mechanics , 1949 .

[39]  D. Mumford,et al.  Optimal approximations by piecewise smooth functions and associated variational problems , 1989 .

[40]  Michael J. Black,et al.  The Robust Estimation of Multiple Motions: Parametric and Piecewise-Smooth Flow Fields , 1996, Comput. Vis. Image Underst..

[41]  G. Lauder,et al.  The hydrodynamics of eel swimming , 2004, Journal of Experimental Biology.

[42]  J. Sethian Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[43]  Martin D. Fox,et al.  Erratum: Distance regularized level set evolution and its application to image segmentation (IEEE Transactions on Image Processingvol (2010) 19:12 (3243-3254)) , 2011 .

[44]  K. Schulze-Delrieu Visual parameters define the phase and the load of contractions in isolated guinea pig ileum. , 1999, American journal of physiology. Gastrointestinal and liver physiology.

[45]  David E. Breen,et al.  A Level-Set Approach for the Metamorphosis of Solid Models , 2001, IEEE Trans. Vis. Comput. Graph..

[46]  Neil A. Thacker,et al.  Tutorial: Computing 2D and 3D Optical Flow. , 2004 .

[47]  S. Osher,et al.  An improved level set method for incompressible two-phase flows , 1998 .