E. L. Bouzarth's Research Page




My research combines fluid dynamics and numerical analysis to study physical problems often motivated by biology. 

Immersed Semi-Flexible Elastic Filaments

Many physical and biological systems involve inextensible fibers immersed in a fluid; examples include cilia, polymer suspensions, and actin filament transport. In such systems, the dynamics of the immersed fibers may play a significant role in the observed macroscale fluid dynamics. In this study, we simulate the dynamics of an approximately inextensible semi-flexible fiber immersed in a two-dimensional cellular background flow. The system is modeled as an immersed boundary problem with the fluid dynamics described using the Stokes equations. The motion of the immersed fiber is computed by means of the method of regularized Stokeslets, which allows one to calculate fluid velocity, pressure, and stress in the Stokes fluid flow regime due to a collection of regularized point-forces without computing fluid velocities on an underlying grid. Simulations results show that, for some parameter values, the fiber may buckle when approaching a stagnation point. These results provide insight into the stretch-coil transition and macroscale random walk behavior that have been reported in mathematical and experimental literature.


An immersed semi-flexible elastic filament immersed in cellular Stokes flow near a hyperbolic stagnation point.  The white arrows depict fluid velocity while the color contour plot relates to stress.

Hydrodynamics of Cilia

My doctoral research was inspired by the Virtual Lung Project at the University of North Carolina at Chapel Hill. In this endeavor, a diverse interdisciplinary group of researchers is studying a variety of components of the human lung. The goal is to develop a computational model that will allow future scientific breakthroughs pertaining to the prevention, treatment, and cure of pulmonary malfunction, specifically, cystic fibrosis. Of the many research aspects of the Virtual Lung Project, investigating the hydrodynamic effects of cilia on surrounding liquid is the area most closely related to my work. Cilia in the lung are small whip-like protrusions responsible for moving mucus and contaminants out of the lung. Besides pulmonary cilia, primary nodal cilia precess conically and are responsible for the development of left-right asymmetry in embryos.



Sketches of pulmonary (left) and primary nodal (right) cilia motion.
(left image courtesy of http://www.p-c-d.org/en/brochure)

Implementing regularized singularities to model Stokes flow generated by a precessing rod is intended to numerically simulate a situation for which colleagues have exact mathematical solutions and experimentalists have corresponding laboratory studies on the micro- and macro-scales. By strategically distributing regularized singularities in a fluid domain to mimic an immersed boundary, one can compute the velocity and trajectory of the fluid at any point of interest. In addition, a no-slip plane can be created with additional regularized singularities to better match the experimental setup. The exact solution allows for careful error analysis and the experimental studies provide new applications for the numerical model. Spectral deferred correction methods are used to alleviate time stepping restrictions in trajectory calculations. Quantitative and qualitative comparisons to theory and experiment have shown that a numerical simulation of this nature can generate insight into fluid systems that are too complicated to fully understand via experiment or exact numerical solution independently.

Mechanics of Running

Running is a topic that interests researchers in many fields, including sociology, anthropology, sports medicine, anatomy, physics, and mathematics. Running shoe technology and its effect on the evolution of the human running stride is a topic of current conversation in scientific and anthropological studies. In collaboration with an undergraduate student as part of the Duke University PRUV Program, we developed a simple model of the lower leg, foot, and ankle to study various aspects of running. The structure of the model developed from fundamental principles of physics: force, torque, momentum, and energy. This research investigates the force distribution before, during, and after ground impact, which may provide insight into injury prevention. By looking at the transfer and loss of energy during a running stride, one may be able to make observations relating to factors affecting efficiency and performance. Incorporating different foot-strike patterns, physiological parameters, shoes, and running surfaces into our model provides an opportunity to address questions of efficiency and injury prevention.



Research Publications:

  • "Modeling a Semi-Flexible Filament in Cellular Stokes Flow Using Regularized Stokeslets." E.L. Bouzarth, A.T. Layton, and Y.-N. Young. Submitted for publication, 2010.
  • "Modeling Non-Slender Bodies with the Method of Regularized Stokeslets."  E.L. Bouzarth and M.L. Minion.  Submitted for publication, 2010.
  • "Modeling Slender Bodies with the Method of Regularized Stokeslets."  E.L. Bouzarth and M.L. Minion. Submitted for publication, 2010.
  • "A Multi-rate Time Integrator for Regularized Stokeslets."  E.L. Bouzarth and M.L. Minion.  Journal of Computational Physics. Vol. 229, No. 11, pp. 4208-4224, 2010. doi: 10.1016/j.jcp.2010.02.006
  • "Regularized Singularities and Spectral Deferred Correction Methods: A Mathematical Study of Numerically Modeling Stokes Fluid Flow." E.L. Bouzarth, Ph.D. dissertation, University of North Carolina at Chapel Hill, 2008. http://search.lib.unc.edu/search?R=UNCb5803305
  • "Epicyclic orbits in a viscous fluid about a precessing rod: Theory and experiments at the micro and macro scales." E.L. Bouzarth, A. Brooks, R. Camassa, H. Jing, T.J. Leiterman, R.M. McLaughlin, R. Superfine, J. Toledo, and L. Vicci. Physical Review E, Col. 76, 2007, 016313
  • "Helicity conservation under Reidemeister Moves.” E.L. Bouzarth and H. Pfister. American Journal of Physics, Vol. 74, No. 2, 2006. doi: 10.1119/1.2142691
  • "Topological Helicity for Framed Links.” E.L. Bouzarth and D. Richeson. Journal of Knot Theory and its Ramifications, Vol. 13, No. 8, 2004, 1007-1019.  doi: 10.1142/S0218216504003664

Research Presentations:

  • Workshop on Fluid Motion Driven by Immersed Structures, Fields Institute, Toronto, Canada, August 2010
  • Society for Industrial and Applied Mathematics Annual Meeting, Pittsburgh, PA, July 2010
  • Frontiers in Applied and Computational Mathematics, New Jersey Institute of Technology, May 2010
  • Society for Industrial and Applied Mathematics Southeastern Atlantic Sectional Meeting, Raleigh, NC, March 2010
  • American Mathematical Society Southeastern Sectional Meeting, Raleigh, NC, April 2009
  • Duke University, Applied Mathematics and Analysis Seminar, October 2008
  • Virginia Commonwealth University, Mathematics Department Colloquium, February 2008
  • University of Utah, Mathematical Biology Seminar, February 2008
  • Northwestern University, Applied Mathematics Colloquium, January 2008
  • Dickinson College, Sigma Pi Sigma Induction Ceremony Keynote Address, November 2007
  • South Eastern Atlantic Mathematical Sciences Workshop, National Institute of Aerospace, October 2007
  • Graduate Mathematics Association Seminar, University of North Carolina at Chapel Hill, October 2005
  • Math and Computer Science Colloquium Series, Dickinson College, April 2003 

Research Posters:

  • Fluid Dynamics, Analysis, and Numerics, Duke University, June 2010
  • Association for Women in Mathematics Workshop for Women Graduate Students and Recent PhDs, Joint Mathematics Meetings, San Diego, CA, January 2008
  • Workshop for Young Researchers in Mathematical Biology, Mathematical Biosciences Institute at The Ohio State University, September 2007
  • University Research Day, University of North Carolina at Chapel Hill, March 2007
  • South Eastern Atlantic Mathematical Sciences Workshop, College of Charleston, September 2006
  • South Eastern Atlantic Mathematical Sciences Workshop, University of North Carolina at Chapel Hill, September 2005


Contact ELB: bouzarth @ math.duke.edu