National Academies Press: OpenBook
« Previous: 4 Modeling with Compartments
Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×

5

From the Compartment to the Fluid

The mathematical modeling presented in Chapters 3 and 4 was based mainly on compartment models, models that lead to differential equations with unknown rate parameters. Another approach, taken by two other workshop presenters, is to investigate the fluid dynamics of the underlying biological system and use those principles to enhance our biological understanding.

George Oster, of the University of California at Berkeley, started his talk by showing a computer-generated movie representative of myxobacteria movement and then presenting a mathematical model that describes the collective behavior during the “ripple phase” of development. The ripple phase (see Figure 5-1) is characterized by a unique pattern of waves that appear to pass through one another; moreover, they can occur without any net movement of bacteria (see Igoshin et al., 2001).

Knowing that myxobacteria move with a combination of two motility mechanisms, labeled A and S, which are controlled by different physiology, Oster sought to combine mathematical models of these mechanisms with a model of the bacterial communication system into a dynamic model that could produce the ripple phase. The resulting model succeeded in emulating and explaining important characteristics of the motion of myxobacteria, such as how the bacteria in crests move with the wave while the bacteria in troughs move against the wave. When waves collide, some bacteria continue moving forward and others reverse, in accordance with experimental observations. The model also captures the ability of ripple waves to move through one another and to propagate without any net transfer of mass.

In another application of fluid dynamics, Charles Peskin, of New York University, described a mathematical model for the heart that considers the muscle tissue to be a time-dependent elastic material, which can be modeled using fluid dynamics. The geometry of his model builds on work in the 1950s by Carolyn Thomas, which described the fiber architecture of the heart as a system of spiraling muscle fibers (see Figure 5-2).

By considering the heart as a composite material of fiber and fluid, Peskin developed equations of a viscous, incompressible fluid to describe the force applied by the fibers to the fluid. Figure 5-3 illustrates the relationship between pressure and volume and how it changes as the valves open and close, while Figure 5-4 shows the force vectors, illustrating the forces on the muscle fiber of the heart.

Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×

FIGURE 5-1 Micrograph and model output showing patterns in the ripple phase. Figure courtesy of George Oster.

FIGURE 5-2 Muscular architecture of hog ventricle. Left, anterior view; right, posterior view. Reprinted with permission from Carolyn Eyster Thomas, Muscular architecture of the ventricles of hog and dog hearts, American Journal of Anatomy, Vol. 101 (1957), 17-57.

Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×

FIGURE 5-3 Pressure/volume relationship in the heart. Figure courtesy of Charles Peskin.

FIGURE 5-4 Illustration of forces on heart muscle fibers. Figure courtesy of Charles Peskin.

Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×

Using the following definitions,

q,r,s = material parameters (s varies along the fiber)

t = time

x = x(q,r,s,t), position

τ = τ(q,r,s,t), force

T = T(q,r,s,t), stress

Peskin derived the following differential equations:

  • The fluid equation, in Cartesian coordinates, which describes the force of the fibers on the fluid:

  • The fiber equations, which describe the stress on the fibers:

  • The interaction equation, which ties them together:

Note that the equations are expressed using the Dirac delta function, an alternative to computing a Jacobean. Use of the delta functions results in an algorithm that is numerically more stable than one that relies on a Jacobean.

Peskin also described a numerical method for solving the equations, using a second-order immersed boundary method derived by M.-C. Lai, which is based on an extension of the Runge-Kutta method. This method can also be used in traditional fluid mechanics problems.

The model output (example frames shown in Plates 1 and 2) compares admirably with MRI scans. Perhaps most impressive, the model captures the swirling movement of blood within the ventricles, a phenomenon of physiological importance that had not emerged in simpler models. This swirling explains why less force is required to exit the chamber than would otherwise be predicted, and it also eliminates leakage in valves from back pressure. However, some improvements are still needed in the model, including refining the representation of the up-down valve motion, getting a larger valve opening, and decreasing the movement of the base of the heart.

It is interesting to note that the Peskin work reflects two of the main themes common to other talks at the workshop. First, the evolution of his heart model over the years has been guided by a clear

Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×

understanding of what physiological features must be captured and what mathematical methods might be suitable to do so, rather than by, say, simply increasing model resolution. For instance, his previous model, although doing an adequate job of reproducing some heart actions, displayed an irregular pattern of blood flow and leakage in the valves that was significantly different from the pattern in a functioning heart. His current version has become more complex in a way that can capture this important phenomenology. Second, the Peskin model illustrates how a mathematical representation can suggest insights (subject to experimental validation) that would not be apparent from current experimental data. The Peskin heart model is an approximate surrogate for a beating heart, a surrogate that can be manipulated and inspected in ways that a living heart cannot.

Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×
Page 18
Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×
Page 19
Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×
Page 20
Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×
Page 21
Suggested Citation:"5 From the Compartment to the Fluid." National Research Council. 2002. Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems. Washington, DC: The National Academies Press. doi: 10.17226/10356.
×
Page 22
Next: 6 Gene Transfer as a Biomedical Tool »
Making Sense of Complexity: Summary of the Workshop on Dynamical Modeling of Complex Biomedical Systems Get This Book
×
Buy Paperback | $21.00 Buy Ebook | $16.99
MyNAP members save 10% online.
Login or Register to save!
Download Free PDF

On April 26-28, 2001, the Board on Mathematical Sciences and Their Applications (BMSA) and the Board on Life Sciences of the National Research Council cosponsored a workshop on the dynamical modeling of complex biomedical systems. The workshop's goal was to identify some open research questions in the mathematical sciences whose solution would contribute to important unsolved problems in three general areas of the biomedical sciences: disease states, cellular processes, and neuroscience. The workshop drew a diverse group of over 80 researchers, who engaged in lively discussions.

To convey the workshop's excitement more broadly, and to help more mathematical scientists become familiar with these very fertile interface areas, the BMSA appointed one of its members, George Casella, of the University of Florida, as rapporteur. He developed this summary with the help of two colleagues from his university, Rongling Wu and Sam S. Wu, assisted by Scott Weidman, BMSA director.

This summary represents the viewpoint of its authors only and should not be taken as a consensus report of the BMSA or of the National Research Council.

  1. ×

    Welcome to OpenBook!

    You're looking at OpenBook, NAP.edu's online reading room since 1999. Based on feedback from you, our users, we've made some improvements that make it easier than ever to read thousands of publications on our website.

    Do you want to take a quick tour of the OpenBook's features?

    No Thanks Take a Tour »
  2. ×

    Show this book's table of contents, where you can jump to any chapter by name.

    « Back Next »
  3. ×

    ...or use these buttons to go back to the previous chapter or skip to the next one.

    « Back Next »
  4. ×

    Jump up to the previous page or down to the next one. Also, you can type in a page number and press Enter to go directly to that page in the book.

    « Back Next »
  5. ×

    Switch between the Original Pages, where you can read the report as it appeared in print, and Text Pages for the web version, where you can highlight and search the text.

    « Back Next »
  6. ×

    To search the entire text of this book, type in your search term here and press Enter.

    « Back Next »
  7. ×

    Share a link to this book page on your preferred social network or via email.

    « Back Next »
  8. ×

    View our suggested citation for this chapter.

    « Back Next »
  9. ×

    Ready to take your reading offline? Click here to buy this book in print or download it as a free PDF, if available.

    « Back Next »
Stay Connected!