|Home | About | Journals | Submit | Contact Us | Français|
Biomechanics is broadly defined as the scientific discipline that investigates the effects of forces acting on and within biological structures. The realm of biomechanics includes the circulatory and respiratory systems, tissue mechanics and mechanotransduction, and the musculoskeletal system and motor control. As in many other biological phenomena, many spatial scales are crossed by biomechanics research: intracellular, multicellular, and extracellular matrices; and tissue, organ, and multiorgan systems. It is well established that the effect of forces at higher scales influence behavior at lower scales and that lower-scale properties influence higher-scale response. However, computational methods that incorporate these interactions in biomechanics are relatively rare. In general, computational models that include representation of multiple spatial or temporal scales are loosely defined as multiscale. The fact that multiscale modeling is not well defined lends the term to a variety of scenarios within the computational physiology community. In biomechanics, multiscale modeling may mean establishing a hierarchical link between the spatial and temporal scales, while the output of a larger-scale system is passed through a finely detailed representation at a lower scale (e.g., body-level movement simulations that provide net joint loading for tissue-level stress analysis). In reality, multiscale modeling may require more intricate representation of interactions among scales. A concurrent simulation strategy is inevitable to adequately represent nonlinear associations that have been known for decades .
Multiscale modeling has existed for many years in basic science and engineering areas such as mathematics, material science, chemistry, and fluid dynamics. Computational and organizational issues common to all these disciplines have been explored, including but not limited to, the standardization of methodology, the necessity of reliable data collection procedures, the need for efficient and accurate algorithms, the lack of coupling tools that address multiphysics phenomena, model and data sharing, and public dissemination. Multiscale biomechanics not only shares all of these common problems, but it is also hindered by the restricted amount of data collection possibilities for model development and validation, the large variability in anatomical and functional properties, and the readily nonlinear nature of the underlying physics even at single scales.
In this article, we describe some current multiscale modeling issues in computational biomechanics from the perspective of the musculoskeletal and respiratory systems and mechanotransduction. First, we justify the requirement of multiscale simulations in these individual systems. Then, we summarize challenges inherent to multiscale biomechanics, followed by system-specific computational challenges. We discuss some of the current tools that have been used to aid research in multiscale mechanics simulations and the priorities to further the field of multiscale biomechanics computation. Overall, our goal is to portray our understanding of the highly complicated and time-sensitive discipline: the so-called multiscale biomechanics modeling.
Musculoskeletal modeling can provide the outlining principles of locomotion including movement control and loading on the hard and soft tissues and muscles. Commonly represented at the body level, these models typically use simplified representations of joints (e.g., hip joint as a spherical joint), passive structures (e.g., modeling of ligaments as nonlinear springs), muscles (e.g., hill-type descriptions), and motor control strategies (e.g., calculation of muscle forces using optimization). If the goal is an overall explanation of muscle function and movement at the body level, the added computational and development costs of increasing the level of detail (therefore introducing multiscale modeling) may not be warranted.
However, there are cases that warrant multiscale modeling in the analysis of the musculoskeletal system. For example, one may be interested in individual muscle fiber function  or the stress–strain profile at the joints  during locomotion. There are also scenarios in which models of muscle coordination coupled with the detailed representation of joints and tissues are needed. In these cases, the interdependency of muscle force and tissue response justifies a concurrent multiscale-modeling approach. As an example, patellofemoral pain (PFP) is a common disorder of the knee whose multifactor etiology is not well understood. It is believed that one mechanism of PFP is excessive stress in the patellar cartilage. Both muscle activation  and muscle reflex response times  have been associated with PFP. In addition, the location of pain receptors in the patellar subchondral bone , the influence of joint tissues (such as the medial patellofemoral ligament ), and the complexity of calculating cartilage stress indicate that a multiscale approach would be beneficial. Temporomandibular joint disorders provide a similar example where complex interdependencies exist between the temporomandibular joint disc and activations of the powerful masticatory muscles . It is generally believed that neuromuscular control is a significant factor in noncontact anterior cruciate ligament (ACL) injury . Understanding noncontact injury mechanisms could be enhanced with multiscale models that include detailed representation of muscle (wrapping, activation, and fiber orientation) coupled with accurate representation of the ACL (interaction with surrounding tissue, insertion areas, fiber orientations, viscoelasticity, and damage accumulation) within a body-level framework. Diabetic foot ulceration provides another example where the interaction of muscle coordination and tissue deformation is important. It is well known that diabetic foot ulceration has a mechanical etiology . Patients with diabetes have to perform similar activities of daily living as healthy individuals. The simple task of walking may be harmful because diabetes can affect various levels of biological function from a mechanical perspective. Dysfunctions at these levels manifest themselves in terms of loss of sensation , changes in control of movement , and alteration of tissue  and cell properties . It is not clear how system-level mechanical loads (e.g., contact forces at the foot) reflect to cellular deformations that may cause cell damage and, therefore, ulceration. Higher–organ level forces (e.g., increased foot pressures), redistribution of stress due to changes in tissue composition (e.g., muscular atrophy ) as well as cell distribution within a tissue, and increased mechanical loading of cells or their decreased damage resistance may all have a role in ulceration. A multiscale modeling approach is likely to identify the pathways to cell damage from organ-level mechanical loading to cell-level deformations.
Holistic simulation of all the aforementioned conditions requires models that optimize neuromuscular response concurrently with detailed models of dynamic tissue behavior, thus emphasizing a multiscale approach in musculoskeletal biomechanics. A further requirement is that these multiscale models have sufficient computational efficiency for optimization-type simulations.
Gas exchange at the respiratory surfaces of the lung is dependent on adequate matching of ventilation and perfusion through complex branching structures that are physically tethered to the surrounding parenchymal tissue. Ventilation, perfusion, and gas exchange are, therefore, intimately dependent on the relationship between stress and strain in the lung parenchymal, airway, and vascular tissues; how this varies regionally; and how it changes with disease. Fredberg and Kamm  recently provided a comprehensive review of stress transmission in the lung, from cell to organ. The review highlights the current state of knowledge of the lung as a mechanical organ with organ-specific interdependencies that arguably make it the most complex system of the human body in which to compute solid mechanics.
The lung parenchymal tissue is extremely delicate; yet, it is required to undergo relatively large strain during the repeated action of ventilation. The tissue accommodates change in the tone of the airways or vasculature through a far lower resistance to shear (and therefore to shape change) than to volume change . The lung tissue deforms readily (because of gravity) with a change of posture, and the resulting regional differences in volume expansion of the lung partly determine the distribution of inspired air. The bronchi and blood vessels are elastic structures that are subjected to internal air and blood pressures, respectively, and through parenchymal tethering they are also subjected to dynamic expanding forces transmitted from the pleural surface. Transmission of force to the airways or vessels depends on the integrity of this tethering: respiratory diseases such as asthma or emphysema disrupt mechanical tethering, but through different mechanisms. Tethering of the arteries and veins extends to the level of the pulmonary capillaries that essentially form the walls of the alveoli. Lung tissue expansion causes expansion of the arteries and veins, but the opposite is true for the capillaries, where tissue expansion at the microstructural level equates to expansion of the alveolar septae (and also potentially some rearrangement of alveolar wall conformation) and results in a reduction in the height of the capillaries. The cyclic motion of the breathing lung also imposes dynamic forces at the cellular level, affecting the regulation of structure, function, and metabolism in a variety of cell types. For example, shear stress through airflow is believed to be transmitted to the epithelial cells, with a role in regulating airway surface liquid, and stretch-sensitive ion channels in a number of cell types regulate an increased influx of calcium under mechanical stretch . Research in cellular mechanotransduction, as described in detail in the following section, provides a wealth of information on the cellular response to stretch or shear, but it is difficult to relate this knowledge to the function of the whole organ without a multiscale computational framework in which to interpret cell-level measurements. The previous examples have considered the influence of the organ and tissue on lower-level structures and cell. An example in the opposite direction is the effect of bronchoconstriction on parenchymal mechanics via tethering of the airway wall to the tissue.
The requirement for multiscale representation in lung mechanics is therefore apparent: from lung interaction with heart, chest wall, and diaphragm to the organ and its internal structures; from the mechanical behavior of a complex functional tissue and the major role that surface forces play in determining this behavior, down to the level of the variety of lung cells that respond to dynamic mechanical forces.
The adaptation of tissue properties due to cell function and mediated by the mechanical and biochemical environment has long been recognized. The phenomenon of mechanotransduction is associated with many normal and pathologic processes including bone remodeling, cardiovascular development, and wound healing. However, though the roots of the theory of mechanotransduction are more than 100 years old , mathematical formulations for the theory that can address the complexities of soft- and hard-tissue biomechanics have entered into the modeling domain fairly recently.
Continuum manifestation of mechanotransduction is driven by cellular and subcellular events, which immediately suggests two unique modeling strategies: capture the continuum-level behavioral characteristics using mathematical laws of growth and modeling or model fundamental cell–cell interactions on a local level and allow the continuum properties to evolve accordingly. The latter approach has been used in several applications, including mesenchymal morphogenesis  and trabecular bone adaptation , . These analyses are driven by local, relatively simple differential equations that govern the evolution of, for example, bone or cell density, from which continuum-level patterns may emerge. The former approach relies on solving continuum-level equilibrium equations and is reviewed further below. The multiscale nature of this problem, then, is the efficient interfacing between these two approaches.
The effects of mechanotransduction include growth and remodeling, which are typically considered as unique processes , representing mass–volume changes due to bulk material deposition or resorption versus structural changes including trabecular or fiber realignment, respectively. Development of models of these processes involves both constitutive formulation (constructing mathematical models that govern the evolution of the state of the material) and computational implementation of the constitutive model, typically within a finite element framework with a few exceptions .
Basic constitutive formulations have been posed in terms of evolution laws for the geometry of a basic construct, for example, a cylindrical vessel . In addition to being a continuum formulation for what is inherently a multiscale, multiphysics process, these formulations assume homogeneity of material response. More rigorous constitutive formulations look at the pointwise response of the material, although still in the continuum sense. The continuum mathematical concept of the deformation gradient is decomposed into an elastic component and a growth component . The growth component effectively alters the reference state of the material and is, thus, able to capture phenomena such as residual stress and opening angle in vasculature  and fiber recruitment and alignment in engineered tendon constructs . The constitutive problem is then to pose the evolution law for the growth component of the deformation gradient in terms of some metric of the local mechanical state including stress , strain , stiffness , or strain energy density . This approach is similar mathematically to nonlinear inelastic constitutive models , in which the deformation gradient is decomposed into an elastic component and a viscous–plastic component.
Much of the computational work of geometric and structural adaptation of bone in response to mechanical loading has involved alteration of the mechanical properties (e.g., density and elastic modulus) of a constituent element or voxel in the computational model, and net addition or subtraction of that element or voxel . Computational models of growth or remodeling of soft tissue are also typically finite element based and can be realized by incorporating user-defined material properties reflective of the inelastic processes into general purpose finite element analysis software . These models are again typically continuum based, employing laws for local element adaptation based on the continuum state.
From a different perspective, mathematical and computational models for cell–matrix interactions have been used to predict cell migration within anisotropic tissue, tissue deformation due to the contractile properties of adherent cells, and fiber reorientation . This approach has been used for hard tissue (bone) to predict remodeling  and callus formation  and can also be used in conjunction with fiber-based constitutive models for soft-tissue behavior to predict evolution of material properties due to cell–matrix interactions. However, such an approach is again continuum based.
The case for multiscale modeling of mechanotransduction is based on the physiological underpinnings of the process itself. Although not yet fully understood, mechanical stimuli are transduced into the cell through structural (integrins that mechanically link the extracellular matrix to the cytoskeleton and in turn the nucleus) and biochemical (stress-based activation of transmembrane ion channels or surface growth factor receptors) pathways . Once internalized, a cascade of intracellular processes ensues, which drives cell function including motility, contraction, proliferation, differentiation, and fibrillogenesis. Continuum-based formulations for mechanotransduction treat these cellular and subcellular processes as black boxes. Although such treatment has been shown to successfully capture tissue-level aspects of growth and remodeling, it has been argued that governing the response of cells using mathematical concepts such as stress or strain is inherently unsatisfactory . Additionally, mathematical coupling of fundamental cellular processes with their effects on the geometric, structural, and constitutive environment, and in turn continuum-level boundary conditions on cell function, can be used to drive functional tissue engineering, understand and better treat pathologic growth or remodeling (hypertrophy), and develop combined mechanical and biological solutions to trauma and joint reconstruction.
Using computer models to simulate biological phenomena, regardless of spatial and temporal scales, involve the associated costs in labor and computation. First, any modeling and simulation platform relies on well-developed tools that facilitate defining models and allow simulations to be conducted in reasonable time frames. Developing robust tools is labor intensive and requires expertise in computational science and mathematics. Given model development and numerical solution tools, another labor-intensive step relates to the development of a model, which is directly related to the research discipline. Compiling adequate model input parameters and representing anatomical structures (e.g., generating meshes of complex geometry) are all part of this process. This is usually a hidden cost, not necessarily reported with the results of the studies. After model development, the next step is to use the model to simulate conditions to answer clinical or research problems. Simulations, particularly in multiscale modeling where coupling between physical domains and scales is necessary, are computationally intensive. Solutions may be obtained, but their interpretation may also be challenging, relying on an expert, possibly spending hours to confirm the validity of results and then extracting useful information applicable to the research area. All these general challenges can be addressed in a research setting, if the development cost of tools and models against decreased research output can be afforded. However, the urgent nature of clinical problems increases the burden on model developers, highly demanding easy-to-use and robust models combined with timely solutions.
Modeling in tissue mechanics requires reliable software to prepare models, to solve them, and later to process the results. Preparation of models usually requires data to reconstruct anatomical geometry, to represent characteristics specific to a physical domain of interest, to describe simulation conditions, and for validation. Hardware for model preparation and visualization of results usually requires graphical processors capable of handling large datasets. In general, the simulation process needs hardware with high-speed numerical computation capabilities and a large memory. It is not uncommon to use high-performance computers with multiple processors and parallel processing using shared or distributed memory architectures.
For multiscale modeling and simulations, software, hardware, and data needs increase. In general, specific and specialized computational tools are used at each dimensional scale of the musculoskeletal system. Linking and concurrently passing information across scales, each using the diverse computational methods, presents a challenge. The primary computational tools of musculoskeletal modeling include multibody dynamics at the body level, continuum finite element methods at the organ level, and specialized algorithms and solvers at the tissue–cell level. The multibody method is computationally efficient but lacks the complexity to accurately capture tissue behavior. Finite element models of organs can estimate tissue deformation but are typically intensive in both model preparation and computational time. Combinations of measured electromyographic activity (EMG), ground reaction forces, kinematics, musculotendon dynamics models, muscle activation optimization schemes, and forward or inverse dynamics can predict net joint loading and the forces of individual muscles. Joint loading and muscle forces can then provide input to finite element models that calculate tissue deformation. For example, predicted quasistatic muscle forces and relative bone displacements have been used to provide boundary conditions for finite element prediction of cartilage stress . In this scheme, parameters at the organ and tissue levels are not part of the muscle control strategy. The finite element method and optimization-based force prediction can be coupled , but the computational cost of repeatedly solving the finite element model is prohibitive. In addition, unless time history is also passed between the separate computational domains, viscoelastic behavior and contact friction cannot be represented.
At the body level, predictive simulation of musculoskeletal movements is possible by using forward dynamics and optimal control of muscle activations . Such simulations are able to predict muscle control patterns for performance-related activities such as maximum height jumping  and for efficient movements such as walking with minimum energy expenditure . These simulations are already costly because of repeated integration of the equations of motion to solve for an optimal muscular control pattern. Nevertheless, adding another level of complexity by introducing models at tissue or even cell scales have practical implications. For example, one can design rehabilitation strategies that fine tune system-level loading to promote healing through cell-level remodeling. In addition, safe movement strategies can be predicted to prevent tissue-level failures. Concurrent simulations, coupling forward dynamics of a musculoskeletal multibody system with tissue-level finite element models, are possible (Figure 1). However, this increases the computational cost of movement prediction using single forward dynamics solution from seconds to hours . The optimal control procedure, which may be seeking minimum tissue stress, can require hundreds if not thousands of these simulations.
Validation of multiscale musculoskeletal models depicts another significant challenge. The current state of in vivo data collection for the musculoskeletal system comprises muscle activations, body segment motion, and ground reaction forces. In addition, musculoskeletal organ geometries can be obtained through magnetic resonance imaging (MRI), although high-resolution images require in vitro MRI or computed tomography (CT) scans. With the exception of a few limited instances, muscle forces cannot be measured in vivo. Models that involve estimation of muscle force typically verify model predictions through the measurement of muscle activation (EMG) . Organ, tissue, and subtissue data collection is typically in vitro. This may include mechanical properties of cortical bone, orientation of collagen fibers, and protein content. The lack of in vivo subject-specific data and the complexity and cost associated with experimental measurements emphasize the need for data repositories for musculoskeletal biomechanics. Parameter sensitivity studies coupled with statistical populations of in vivo and primarily in vitro data may provide feasible validation routes.
In the study of the mechanics of the lung tissue, the main challenge remains as it stood nearly 25 years ago: “first, stress–strain relations, based on independent material testing or microstructural modeling with some conformational testing, are needed.” . Early computational studies of lung–tissue mechanics used linear elastic theory and linear material constants, treating the lung as a uniformly inflated structure subjected to an incremental deformation . The approach was extended by solving for successive small displacement increments and using elastic moduli that were dependent on the associated incremental change in transmural pressure . In reality, the lung undergoes relatively large strain during normal breathing, requiring many of these small increments and potentially accumulating numerical error with each increment; however, the appeal of the linear elastic approach is clear: the governing equations are simpler than equations valid for large deformations, enabling analytic solutions for simple shapes and loading , and the elastic moduli are obtained relatively easily .
However, an accurate analysis requires the use of governing equations that are valid for large deformations, i.e., finite deformation elasticity. Relatively few studies have used large deformation theory, and all of these have had to compromise to some degree in their constitutive law . The elasticity of the lung parenchyma comprises a tissue contribution and surface forces due to the alveolar air–liquid interface. The magnitude of the surface forces changes with surface area, and surface tension distorts the alveolar geometry as the tissue expands or recoils . A mathematical model of the tissue microstructure that includes interaction between alveolar structure and surface tension has been developed , but to date, this has not been extended to link with the tissue-continuum level. Such a link could be made using the approach presented in , where moduli were derived based on an analysis of the elastic response of a model of the lung microstructure. Experimental testing of lung tissue to determine the relationship between stress and strain is technically difficult and is performed on tissue without its normal surface forces. When there is no surface tension, the tissue component of the parenchymal elasticity is different to that in the intact, in vivo lung, because surface tension results in distortion of the tissue causing more elastic energy to be stored. Gathering sufficient accurate data on which to base the derivation of a constitutive law is therefore exceedingly difficult. Ideally, an analysis of lung tissue mechanics would also account for viscoelasticity and tissue hysteresis, and the different contributions of (and interaction between) tissue elasticity and surface forces.
Software to develop musculoskeletal models and to simulate or analyze movements is maturing and becoming freely available . These tools promise an open architecture, potentially allowing linking with physiologically realistic simulations of tissue–cell deformations and multiscale muscle models. Finite element analysis packages that focus on biological problems are also provided for free , with developers open to implement customization specific to research fields. Nevertheless, current multiscale analysis of the musculoskeletal system is still based on individualization of such tools. The aforementioned examples from other disciplines also illustrate the challenges associated with model preparation and solution. While both these aspects of multiscale modeling can be computationally intensive, model preparation can, in addition, be quite labor intensive, often to the point where the analyst invests more time in model creation than in model analysis. Overall, this is to the detriment of science (consider the challenges presented by the linking of scales between the level of the organ and the cell). Organs themselves typically span multiple geometric scales (compare the size of cardiac valves to the coronary arteries or the size of the trachea to the terminal bronchioles). Efficient grid generation across these geometric scales needs to be somewhat structured so that the physics of both scales are correctly resolved, for example, the transition from convection-dominated to diffusion-dominated mechanics in the lung. With complex biological domains, this can be daunting. Furthermore, because organs are by definition spatially heterogeneous arrangements of cells, mechanical properties with efficient continuum representations such as tissue elasticity or mass transfer emerge in a spatially heterogeneous fashion from the cellular and extracellular constituents of tissue. Some examples are the spatial arrangement of myofibers in the heart, of collagen and elastin in heart valves and arteries, or of the smooth muscle in the lung. Structural constitutive models, as opposed to phenomenological models, of these continuum fields are highly desirable as they can substantially contract the parameter space. However, this requires the efficient identification and mapping of these cell-level components to a computational grid. In cases where cell activity is modeled together with organ-level fields, this mapping must be a two-way communication. This latter scenario demands efficient geometric databases that can organize time varying organ- and cell-level inputs and outputs. Later, we outline some current and evolving efforts to meet these challenges.
Several excellent isotropic unstructured tetrahedral grid generation algorithms exist for creating finite-element or finite-volume grids from imaging-derived closed triangulated manifolds . Isotropic algorithms attempt to construct tetrahedral elements with nearly equal internal angles and approximately equal edge lengths. However, with medical imaging data, the lengths of these edges are typically more related to the resolution of the image than they are to the physics to be solved and do not take the geometric scale into account. Moreover, the tessellation of the volume is typically disordered. Structured grids, typically structured hexahedral grids, are computationally efficient and can be made to fit the physics of the problem, but structured grids can be laborious to construct if they can be constructed at all. Forecasts for petascale computing envision that a fourfold increase in speedup will come from adaptive gridding, whereas only a one and a half fold speedup will come from increased parallelism . These forecasts forcefully apply to large-scale biomedical problems. A compromise between fully disordered tetrahedral grids and structured hexahedral grids are a new class of algorithms that automatically generate structured, scale-invariant tetrahedral grids from medical imaging data . These grids are well suited to large-scale computations as they can be designed to concentrate computational grid where it is needed while keeping the overall computational cost of the problem tractable. They are also structured in the sense that elements can be automatically arranged in a user-specified number of nearly orthogonal layers.
The approach (Figure 2) consists of defining a feature-size field on single or multiple material input triangulated surface mesh derived from single or multiple material Marching Cubes algorithms. This is accomplished without a background grid and without referencing the medial axis. Thus, determination of the feature-size field is not only computationally efficient but also robust in the sense that it is continuous and does not change unreasonably under perturbation of the surface mesh. Before volume mesh generation, the input surface mesh is modified (refined and derefined) so that edge lengths are proportional to the feature-size field, with the constraint that refinement or derefinement preserves topology and curvature. Surface modification is iterated with a volume-conserving smoothing , with the result that surface triangles are well shaped, well organized, and graded. Volume-conserving smoothing also maintains the enclosed volume to within machine precision with respect to the original voxelated volume. From these, possibly multiple, modified surfaces, a user-defined number of structured layers are constructed across arbitrarily oriented cross-sections in the domain, independent of scale. User intervention consists solely in specifying the desired number of layers and the desired element anisotropy. Additionally, a scale-dependent function may be specified. A constant function, for example, would result in a constant number of layers independent of scale. A linear function would proportionally increase or decrease the number of layers at either the top or bottom scales. Functions may be completely arbitrary.
Once an efficient computable grid has been defined, it is necessary to establish a bridge between measurements of cell-level data and the organ-level model. These data can become the inputs to upscaling approaches in which effective reaction–diffusion equations or elasticity tensors are locally defined by homogenizing cells and extracellular proteins over spatial windows of several cells. Alternatively, they can be stored in the computational cells of the organ-level model to drive complimentary cell-level models. In either case, some averaging of data inevitably occurs. In the case where these data derive from imaging modalities such as MRI, positron emission tomography (PET), or CT, volume-to-volume mapping can be done efficiently  and can be made to tolerate modest geometric mismatching. However, often in multiscale models, upscaling is used to define boundary conditions. This requires a conservative volume-to-surface mapping. Along these lines, more work needs to be done.
If voxel to unstructured grid mappings are presently feasible, a much greater potential source of cell-level data is histology . However, histology-based data (fluorescence, autoradiography, in situ hybridization, proteomics, immunohistology, or simple hematoxylin and eosin staining) pose two fundamental challenges: reconstruction and quantitation. Reconstruction is challenging because sections that are thin enough to reveal cell-level detail under analysis are also prone to distortion during the process of sectioning. Some recent algorithm development has been dedicated to accurate reconstructions without a reference geometry, especially in neuroscience . These and related algorithms attempt to affect local nonlinear warp transformations on a slice-by-slice basis in an effort to rectify these distortions. These efforts are promising, but more work is needed to extend and adapt these approaches to other organs such as the heart and lung. Quantitation can be challenging depending on the modality. For example, with in situ hybridization, messenger ribonucleic acid (mRNA) for a specific marker is hybridized to the tissue, and the signal strength is proportional to the local concentration. The ability to quantitate that signal in at least a semiquantitative fashion is essentially an image-processing problem . Here again, work in neuroscience indicates a possible path for multiscale tissue modeling, where markers can be specialized for collagen, elastin, and muscle fibers, and histology can provide both spatial information and local density .
Once these data have been referenced to a single geometry, whether voxelated or not, they may be efficiently mapped to unstructured grids with the same methods as are applied to magnetic resonance (MR) or CT data . However, the sheer volume of data represented by histological data with near cellular resolution is prohibitive for efficient communication between the computational grid and the three-dimensional (3-D) database. A preferable solution is to adopt a multiresolution, grid-based approach , allowing the resolution of the 3-D database to be adapted to the averaging window that works best for the communication between scales.
In summary, multiscale modeling requires, at its foundation, measurement and communication of data between scales. To be useful to the biomedical engineer, the construction of quality computable grids and the one-way or two-way communication of data to and from these data over multiple scales must be fully automatic and accessible. One persistent challenge will be the registration of multiple datasets from different modalities to a common geometric database. Although some of these concepts overlap with atlasing projects, such as the Allen Brain Atlas , there is a need to adapt these approaches to the needs of the modeling community by developing procedures and algorithms that make the communication of data routine.
Computational biomechanical modeling typically requires a level of customization beyond what is possible with commercial tools. This is even more forcefully true with multiscale modeling where codes that are specialized for different scales often must communicate, more and more often in parallel. At the same time, a plethora of laboratory-specific codes entails unnecessary duplication, and commonly less effective codes due to the substantial investment required. In contrast, the open-source paradigm is a proven model for complex software development that has the capacity to effectively create computational tools that are geared to the needs of multiscale modeling, while enabling unlimited customization through access to the source code. A few examples of highly successful open-source projects of relevance to multiscale modeling are the image processing and registration toolkit, Insight Segmentation and Registration Toolkit (ITK) , another similar tool, 3-D Slicer , and the post-processor Paraview . The National Institutes of Health (NIH) sponsored SimTK project  is another example that aims to bring together multiple components geared toward molecular, neuromuscular, and cardiovascular dynamics. A few common aspects of successful open-source projects are that they are funded, supported by the biomedical community, and well organized. The organization of open-source efforts around task-oriented components that are designed to work together [image segmentation, grid generation and management, inverse parameter estimation, computational continuum mechanics, network models, system-level models, ordinary differential equation (ODE) solvers, one-dimensional (1-D) and two-dimensional (2-D) partial differential equation (PDE) solvers, etc.] would greatly enhance the ability of multiscale modelers to focus on biomedical problem solving and discovery, the ultimate goals of biomedical multiscale modeling.
Model sharing also has the potential to empower the multiscale modeler by providing a library of models that can be combined across scales for a given application. However, one roadblock to effective model sharing is the lack of a standard format. This problem has been recognized by computational practitioners outside of the biomedical field and has led to such projects as computational fluid dynamics general notation system (CGNS) . The goals of CGNS, as an example, provide an interesting model for the biomedical modeling community in that they provide, in collaboration with commercial entities, a standard format that is self-descriptive, machine independent, and well documented. This paradigm has been so successful that virtually all major commercial computational fluid dynamics codes support it. However, CGNS is more than a format. It is a comprehensive, open-source data archiving system with standard naming conventions that provide definitions of boundary conditions, solver specifications, and convergence history, such that a solution obtained in one commercial or research code can be restarted directly in another code. To assure consistent implementation, the project provides easy-to-implement midlevel libraries in all major computer languages. Adoption of a common database and format is complimentary to a common interfaces approach and is foundational to successful open-source projects in multiscale modeling.
Dissemination of a solution database with model distribution may have practical value. Multiscale simulations involve many models at the bottom of the solution hierarchy, which usually have an input–output (I/O) relationship with the models at a higher level. Solution for these models is requested frequently while solving for the higher-level model, and in many cases, individual solutions are costly. This can be seen in multilevel finite element models of tissue–cell interactions  and in musculoskeletal modeling where rigid body movement simulations requesting simulations of a tissue-level finite element model at each time step . If a solution database exists, one can build a fast surrogate model to represent the (I/O) relationship. These surrogate representations can be based on a global fit (e.g., response surface , which were commonly used in optimization problems), local regression (e.g., moving least squares ), or neural networks . For example, a surrogate-type modeling approach uses discrete mass–spring damper models (Figure 3) to represent soft tissue. Parameters for the mass–spring damper network are optimized to fit tissue deformations from either experimental measurements or finite element simulations . Alternatively, adaptive strategies can be adopted, which estimate a fit or interpolation error and simulate the complicated model only when needed. By using adaptive surrogate modeling techniques based on local weighted learning , it is now possible to solve optimal control problems that require concurrent simulation of musculoskeletal movements and tissue deformations. Muscle activations for a jumping simulation (Figure 1) can now be predicted by using a surrogate model of foot deformations. This approach requires costly finite element analysis to be conducted only for approximately 30% of the time.
Funding to support contributions of individual authors were as follows: Ahmet Erdemir, NIH (1R01EB006735-01, principal investigator: Antonie J. van den Bogert); Daniel Einstein, NIH (1R01HL073598-01A1, principal investigator: Richard A. Corley); Trent Guess, National Science Foundation, under the Interagency Modeling and Analysis Group (IMAG) program for Multiscale Modeling (CMS-0506297); Merryn Tawhai, NIH, under the IMAG program for Multiscale Modeling (R01-EB-005823, principal investigator: Ching-Long Lin) and a Bioengineering Research Partnership Grant (R01-HL-064368, principal investigator: Eric A. Hoffman).
Merryn Tawhai received her Ph.D. degree in engineering science from the University of Auckland in 2001. Since then, she has been working in the Auckland Bioengineering Institute at the University of Auckland, establishing a research team to develop integrative computational models of the lung for studying pulmonary physiology and pathophysiology. She also holds an adjunct appointment in the Department of Biomedical Engineering at the University of Iowa. Her main research interests have been in developing anatomically based finite element models of physical structures in the lung and using these individualized models to study lung tissue mechanics, ventilation distribution, pulmonary perfusion, airway thermofluid dynamics, and inert gas mixing.
Jeff Bischoff received his Ph.D. degree in mechanical engineering from the University of Michigan. He is currently a senior research engineer at Zimmer, Inc. His current work involves the development and implementation of novel inelastic constitutive models for soft materials including polymers and soft tissue and the use of computational techniques within joint reconstruction applications. His research interests include constitutive and computational modeling of soft tissue, orthopedic biomechanics, multiscale modeling of growth and remodeling, experimental biomechanics, and material optimization.
Daniel Einstein received his Ph.D. degree in bioengineering from the University of Washington. His work includes the development of computational methods for fluid–solid interactions, large deformation mechanics and inverse analysis. Since joining Pacific Northwest National Laboratory (PNNL) in 2005, he has established a small team of researchers to develop computational approaches for image-based multiscale modeling in the rat and mouse respiratory tract and through the application of these tools to understand the health effects of environmental exposure. Some of his active interests are developing novel approaches for biomedical grid generation, image segmentation and registration, nanoparticle mechanics, biological atlases, finite volume computational physics, and computational morphology.
Ahmet Erdemir received his Ph.D. degree in kinesiology from Pennsylvania State University in 2002. He has established and is currently directing the Computational Biomodeling Core in the Department of Biomedical Engineering at the Cleveland Clinic to promote simulation-based medicine. He has training in biomechanics and mechanical engineering, and he is particularly interested in musculoskeletal biomechanics and soft tissue mechanics at multiple scales. He has well-established research experience in computational biomechanics using optimization techniques and finite element analysis. He also has expertise in cadaver experimentation and human subjects testing. He is an active reviewer for biomechanics-related journals and is a member of many professional societies.
Trent Guess received his Ph.D. degree in mechanical engineering from the University of Kansas in 2003. He is currently the lead of the IMAG multiscale modeling consortium tissue mechanics working group and is the advisor of the University of Missouri-Kansas City (UMKC)’s human-powered vehicle team. He is an assistant professor of mechanical engineering at UMKC and holds an adjunct position in the Department of Orthopaedic Surgery in UMKC’s School of Medicine. He is the director of the Musculoskeletal Biomechanics Research Laboratory at UMKC with research interests in the relationship between neuromuscular response and tissue loading. His research interests have led to the development of novel multiscale modeling techniques in musculoskeletal biomechanics.
Jeff Reinbolt received his Ph.D. degree in mechanical engineering from the University of Florida in 2006. He is an engineering research associate at Stanford University where he develops and applies software to study the dynamics and function of human health and disease. He is with the Center for Physics-Based Simulation of Biological Structures, which is one of the seven National Centers for Biomedical Computation supported by the NIH Roadmap for Bio-informatics and Computational Biology. His research interests include biomedical computation, simulation-based surgical and rehabilitation treatment planning, innovative patient-specific modeling, and optimization methods.
This article is written through collaborative editing using the IMAG wiki . The authors are grateful to IMAG for providing this platform and also appreciate the efforts of the developers of Mediawiki , a free wiki engine upon which this platform is based. An online version of this article in pre-print form can be found in the IMAG wiki .
This article does not intend to be complete nor apply to all areas of multiscale modeling. The text is based on the discussions between the authors while sharing their experience to combine multiple physics, domains and scales in their computational biomechanics research.