Areas of Interest
My academic interests fall under the following relatively broad categories:
- Classical theoretical and applied mechanics, Multi-physics phenomena
- Analysis of numerical methods, Advanced finite element methods
- Large scale high performance computing
Biological growth: Reaction, transport and mechanics
H. Narayanan, E. M. Arruda, K. Grosh, K. Garikipati
Received: date / Accepted: date
Abstract
In this paper, we address some modelling issues related to biological growth. Our treatment is based on a recently-proposed, general formulation for growth (Journal of the Mechanics and Physics of Solids, 52, 2004, 1595-1625) within the context of open system continuum thermodynamics. We aim to enhance this treatment by making it more appropriate for the biophysics of growth in soft tissue, specifically tendon. This involves several modifications to the mathematical formulation to represent the reactions, transport and mechanics, and their interactions. We also reformulate the governing differential equations for reaction-transport to represent the incompressibility constraint on the fluid phase of the tissue. This revision enables a straightforward implementation of numerical stabilisation for the hyperbolic, or advection-dominated, limit. A finite element implementation employing a staggered scheme is utilised to solve the coupled nonlinear partial differential equations that arise from the theory. Motivated by our experimental model, an in vitro scaffold-free engineered tendon formed by self-assembly of tendon fibroblasts [Calve et al., 2004], we solve several numerical examples demonstrating biophysical aspects of growth, and the improved numerical performance of the models.1 Introduction
Growth involves the addition or depletion of mass in biological tissue. In biological systems, growth occurs in combination with remodelling, which is a change in microstructure, and possibly with morphogenesis, which is a change in form in the embryonic state. The physics of these processes are quite distinct, and for modelling purposes can, and must, be separated. Our previous work [Garikipati et al., 2004], upon which we now seek to build, drew in some measure from Cowin and Hegedus [1976],Epstein and Maugin [2000],Taber and Humphrey [2001] and Kuhl and Steinmann [2003], and was focused upon a comprehensive account of the coupling between transport and mechanics. The origins of this coupling were traced to the balance equations, kinematics and constitutive relations. A major contribution of that work was the identification and discussion of several driving forces for transport that are thermodynamically-consistent in the sense that specification of these relations does not violate the Clausius-Duhem dissipation inequality. Now, we seek to restrict the range of physically-admissible possibilities in order to gain greater physiological relevance for modelling growth in soft biological tissue. The advection-diffusion equations for mass transport require numerical stabilisation in the advection-dominated regime (the hyperbolic limit). We draw upon the enforcement of the incompressibility limit for the fluid phase to facilitate this process. Below, we briefly introduce each aspect that we have considered, but postpone details until relevant sections in the paper.- For a tissue undergoing finite strains, the transport equations can be formulated, mathematically, in terms of concentrations with respect to either the reference or current (deformed) configuration. However, the physics of fluid-tissue interactions and the imposition of relevant boundary conditions is best understood and represented in the current configuration.
- The state of saturation is crucial in determining whether the tissue swells or shrinks with infusion/expulsion of fluid.
- The fluid phase, whether slightly compressible or incompressible, can develop compressive stress without bound. However, it can develop at most a small tensile stress [Brennen, 1995], having implications for the stiffness and strength of the tissue in tension as against compression. Though this also has implications for void formation through cavitation, the ambient pressure in the tissue under normal physiological conditions ensures that this manifests itself only as a reduction in compressive pressure.
- When modelling transport, it is common to assume Fickean diffusion [Kuhl and Steinmann, 2003]. This implies the existence of a mixing entropy due to the configurations available to molecules of the diffusing species at fixed values of the macroscopic concentration. The state of saturation directly influences this mixing entropy.
- If fluid saturation is maintained, void formation in the pores is disallowed under an increase in the pores' volume. This has implications for the fluid exchanges between a deforming tissue and a fluid bath with which it is in contact.
- Recognising the incompressibility of the fluid phase, it is common to treat soft biological tissue as either incompressible or nearly-incompressible [Fung, 1993]. At the scale of the pores (the microscopic scale in this case), however, a distinction exists in that the fluid is exactly (or nearly) incompressible, while the porous solid network is not.
- In Garikipati et al. [2004], the acceleration of the solid phase was included as a driving force in the constitutive relation for the flux of other phases. However, acceleration is not frame-invariant and its use in constitutive relations is inappropriate.
- Chemical solutes in the extra-cellular fluid are advected by the fluid velocity and additionally undergo transport under a chemical potential gradient relative to the fluid. In the hyperbolic limit, where advection dominates, spatial instabilities emerge in numerical solutions of these transport equations [Brooks and Hughes, 1982,Hughes et al., 1987]. Numerical stabilisation of the equations is intimately tied to the mathematical representation of fluid incompressibility.
2 Balance equations and kinematics of growth
In this section, the coupled, continuum balance equations governing the behaviour of growing tissue are summarised and specialised as outlined in Section 1. For detailed continuum mechanical arguments underlying the equations, the interested reader is directed to Garikipati et al. [2004]. The tissue of interest is an open subset of with a piecewise smooth boundary. At a reference placement of the tissue, , points in the tissue are identified by their reference positions, . The motion of the tissue is a sufficiently smooth bijective map , where . At a typical time , maps a point to its current position, . In its current configuration, the tissue occupies a region . These details are depicted in Figure 1. The deformation gradient is the tangent map of .2.1 Balance of mass for an open system
As a result of mass transport (via the flux terms) and inter-conversion of species (via the source/sink terms) introduced previously, the concentrations, , change with time. In local form, the balance of mass for an arbitrary species in the reference configuration is2.1.1 The role of mass balance in the current configuration
Though it is not mathematically incorrect to solve the initial-boundary-value problem in terms of reference coordinates, it is important to note that as soft tissues deform, the current configuration, , and its boundary, , change in time. Even as the pore structure at the boundary deforms with the tissue, the fluid concentration with respect to remains constant if the boundary is in contact with a fluid bath. Accordingly, this is the appropriate Dirichlet boundary condition to impose under normal physiological conditions. This is shown in an idealised manner in Figure 2.2.2 The kinematics of growth (changes in concentration)
2.2.1 Saturation and tissue swelling
2.3 Balance of momenta
In soft tissues, the species production rate and flux that appear on the right hand-side in Equation (1), are strongly dependent on the local state of stress. To correctly model this coupling, the balance of linear momentum should be solved to determine the local state of strain and stress. The deformation of the tissue is characterised by the map . Recognising that, in tendons, the solid collagen fibrils and fibroblasts do not undergo mass transport, the material velocity of this species, , is used as the primitive variable for mechanics. The motion of each remaining species is split into a deformation along with the solid species, and mass transport relative to it. To this end, it is useful to define the material velocity of a species relative to the solid skeleton as: . Thus, the total material velocity of a species is . The total first Piola-Kirchhoff stress tensor, , is the sum of the partial stresses (borne by a species ) over all the species present. With the introduction of these quantities, the balance of linear momentum in local form for a species in is where is the body force per unit mass, and is an interaction term denoting the force per unit mass exerted upon by all other species present. The final term with the (reference) gradient denotes the contribution of the flux to the balance of momentum. In practise, the relative magnitude of the fluid mobility (and hence flux) is small, so the final term on the right hand side of Equation (6) is negligible, resulting in a more classical form of the balance of momentum. Furthermore, in the absence of significant acceleration of the tissue during growth, the left hand-side can also be neglected, reducing (6) to the quasi-static balance of linear momentum. The balance of momentum of the entire tissue is obtained by summing Equation (6) over all . Additionally, recognising that the rate of change of momentum of the entire tissue is affected only by external agents and is independent of internal interactions, the following relation arises.3 Constitutive framework and modelling choices
As is customary in field theories of continuum physics, the Clausius-Duhem inequality is obtained by multiplying the entropy inequality (the Second Law of thermodynamics) by the temperature field, , and subtracting it from the balance of energy (the First Law of thermodynamics). We assume the internal energy per unit reference volume of species to be of a sufficiently general form: , where is the entropy per unit system volume. Substituting this in the Clausius-Duhem inequality results in a form of this inequality that the specified constitutive relations must not violate. Only the valid constitutive laws relevant to the examples that follow are listed here. For details, see Garikipati et al., [2004].3.1 An anisotropic network model based on entropic elasticity
Each constituent of the tissue has a mass-specific Helmholtz free energy density, . Utilising the material response of a hyperelastic material, the partial first Piola-Kirchhoff stress of collagen is . Here, is the elastic deformation gradient, and is the growth deformation gradient, of collagen. Along the lines of Equation (4), if we were considering unidirectional growth of collagen along a unit vector , we have , with denoting the initial concentration of collagen at the point. The mechanical response (function) of tendons in tension is determined by their dominant structural component, highly oriented fibrils of collagen. In our preliminary formulation, the strain energy density for collagen has been obtained from hierarchical multi-scale considerations based upon an entropic elasticity-based worm-like chain (WLC) model [Kratky and Porod, 1949]. The WLC model has been widely used for long chain single molecules, most prominently for DNA [Marko and Siggia, 1995,Rief et al., 1997,Bustamante et al., 2003], and recently for the collagen monomer [Sun et al., 2002]. The central parameters of this model are the chain's contour length, , and persistence length, . The latter is a measure of its stiffness and given by , where is the bending rigidity, is Boltzmann's constant and is the temperature [Landau and Lifshitz, 1951]. We have fitted the WLC response function derived by Marko and Siggia [1995] to the collagen fibril data of Graham et al. [2004] with nm and nm. This is to be compared with nm and nm, reported by Sun et al. [2002], for a single collagen molecule. Taken together, these results demonstrate that the WLC analysis correctly predicts a collagen fibril to be longer and slightly more compliant than its constituent molecule due to compliant intermolecular cross-links in a fibril. To model the possibility of a collagen network structure, the WLC model has been embedded as a single constituent chain of an eight-chain model [Bischoff et al., 2002a,Bischoff et al., 2002b], depicted in Figure 5. Homogenisation via averaging then leads to a continuum Helmholtz strain energy function, : Here, is the concentration of collagen, is the density of chains, and and are lengths of the unit cell sides aligned with the principal stretch directions. The material model is isotropic only if .3.2 A nearly incompressible ideal fluid
In this preliminary work, the fluid phase is treated as a nearly incompressible, ideal, i.e., inviscid, fluid. The partial Cauchy stress in the fluid is3.2.1 Response of the fluid in tension; cavitation
The response of the ideal fluid, as defined by Equation (9), does not explicitly distinguish between the cases where the fluid is subjected to tension or compression, i.e., whether . When the fluid phase is subjected to compression, being (nearly) incompressible, it can develop compressive stresses without bound and is modelled accurately. Under tension, the fluid can develop at most a small tensile stress [Brennen, 1995], and the bulk of the tensile stiffness arises from the collagen phase.This is not accurately represented by (9), which predicts a tensile response in the fluid similar to the compressive response. Here, we preclude all tensile load carrying by the fluid by limiting . For consistency, we first introduce an additional component to the mixture, a void species, . Denoting its deformation gradient by , we now only require3.3 Constitutive relations for fluxes
From Garikipati et al. [2004], the constitutive relation for the flux of extra-cellular fluid relative to the collagen takes the following form,3.3.1 Saturation and Fickean diffusion
3.3.2 Transport of solute species
The numerous dissolved solute species (proteins, sugars, nutrients, ...), denoted by s, undergo long range transport primarily by being advected by the perfusing fluid. In addition to this, they undergo transport relative to the fluid. This motivates an additional velocity split of the form , where denotes the velocity of the solute relative to the fluid. The constitutive relation for the corresponding flux, denoted by , has the following form, similar to Equation (13) defined for the fluid flux.3.3.3 Objectivity and the contribution from acceleration
In our earlier treatment [Garikipati et al., 2004], the constitutive relation for the fluid flux had a driving force contribution arising from the acceleration of the solid phase, . This term, being motivated by the reduced dissipation inequality, does not violate the Second Law and supports our intuitive understanding that accelerating the solid skeleton in one direction must result in an inertial driving force on the fluid in the opposite direction. However, as defined, this acceleration is obtained by the time differentiation of kinematic quantities, and does not transform in a frame-indifferent manner. Unlike the superficially similar term arising from the gravity vector, the acceleration term presents an improper dependence on the frame of the observer. Thus, its use in constitutive relations is inappropriate, and the term has been dropped in Equation (13).3.3.4 Incompressible fluid in a porous solid
Upon incorporation of the additional velocity split, , described in Subsection 3.3.2, the resulting mass transport equation (3) for the solute species is3.3.5 Stabilisation of the simplified solute transport equation
SUPG methods are a class of finite element methods, originally developed for the scalar advection-diffusion equation, which have proven efficient in the solution of a variety of flow problems [Hughes, 1987]. The standard form of the scalar advection-diffusion is3.4 Nature of the sources
There exists a large body of literature, Cowin and Hegedus [1976],Epstein and Maugin [2000],Ambrosi and Mollica [2002], that addresses growth in biological tissue mainly based upon a single species undergoing transport and production/annihilation. However, when chemistry is accounted for, it is apparent that growth depends on cascades of complex biochemical reactions involving several species, and additionally involves intimate coupling between mass transfer, biochemistry and mechanics. An example of this chemo-mechanical coupling is described in Provenzano et al., [2003]. The modelling approach followed in this work is to select appropriate functional forms of the source terms for collagen, , and the solutes, , that abstract the complexity of the biochemistry. In our earlier exposition [Garikipati et al., 2004], we utilised simple first order chemical kinetics to define . We now replace it with a source that has greater relevance from the standpoint of biochemistry.3.4.1 Enzyme kinetics
Michaelis-Menten enzyme kinetics (see, for e.g., Sengers et al., [2004]) involves a two-step reaction with the collagen and solute production terms given by3.4.2 Strain energy dependent collagen production
A strain energy dependent source term was originally proposed in the context of bone growth [Harrigan and Hamilton, 1993] and induces growth at a point when the energy density deviates from a basal value, suitably weighted by a relative density ratio. Written for collagen, it has the form4 Numerical examples
The theory presented in the preceding sections results in a system of non-linear coupled partial differential equations. A finite element formulation employing a staggered scheme based upon operator splits Armero, [1999,Garikipati and Rao, [2001] has been implemented in FEAP [Taylor, 1999] to solve the coupled problem. In the biphasic-fluid and collagen-problem for instance, the basic solution scheme involves keeping the displacement field fixed while solving for the concentration field using the mass transport equation. The resulting concentration field is then fixed to solve the mechanics problem. This procedure is repeated until the resulting fields satisfy the differential equations within some suitable magnitude of an error norm. The mechanics problem is solved quasi-statically. Backward Euler is used as the time-stepping algorithm for mass transport. Non-linear projection methods [Simo et al., 1985] are used to treat the near-incompressibility imposed by water. Mixed methods, as described in Garikipati and Rao, [2001], are used for stress (and strain) gradient driven fluxes. The following examples aim to demonstrate the mathematical formulation and aspects of the coupled phenomena as the tissue grows. The model geometry, based on our engineered tendon constructs (see Figure 9), is a cylinder 12 mm in length and 1 mm in cross-sectional area. The finite element mesh used for the calculations is depicted in Figure 10.4.1 The constriction induced growth problem
Only two phases-fluid and collagen-are included for the mass transport and mechanics. The collagen is represented by the anisotropic worm-like chain model outlined previously (see Section 3.1) and the fluid phase is modelled as ideal and nearly incompressible. The parameters used in the analysis are as presented in Table 2.1. The values chosen are representative of the kinds of biological systems we are working with. The classes of initial and boundary conditions imposed are also based on physical experiments. Since we only have two species and we want to demonstrate growth, an "artificial" fluid sink is introduced following simple first order kinetics. The collagen source will be the negative of the fluid sink: , where is the reaction rate, and is the initial concentration of fluid. When , this acts as a source for collagen. The mixing entropy of fluid in the mixture with collagen is written as , where is the molecular weight of the fluid. The boundary conditions simply corresponding to immersing the tendon in a nutrient rich bath. The initial collagen concentration is 500 kg/m everywhere and the fluid concentration is 400 kg/m everywhere. This is exposed to a bath where the fluid concentration is 500 kg/m , so with these concentration boundary conditions set, nutrient rich fluid rushes into the tissue, and growth occurs to form more collagen. The following plots present a few results from the analysis.4.2 A swelling problem
4.3 An enzyme-kinetics based multiphasic problem
5 Conclusion
In this paper, we have discussed a number of enhancements to our original growth formulation presented in Garikipati et al. [2004]. That formulation has served as a platform for posing a very wide range of questions on the biophysics of growth. Some issues, such as saturation, incompressibility of the fluid species and its influence upon the tissue response, and the roles of biochemical and strain energy-dependent source terms are specific to soft biological tissues. We note, however, that other issues are also applicable to a number of systems with a porous solid, transported fluid and reacting solutes. Included in these are issues of current versus reference configurations for mass transport, swelling, Fickean diffusion, fluid response in compression and tension, cavitation, and roles of permeabilities and mobilities. These issues have been resolved using arguments posed easily in the framework derived in Garikipati et al. [2004]. The interactions engendered in the coupled reaction-transport-mechanics system are complex, as borne out by the numerical examples in Section 4. We are currently examining combinations of sources defined in Section 3.4, and aim to calibrate our choices from tendon growth experiments. The treatment of these issues has led to a formulation more suited to the biophysics of growing soft tissue, making progress toward our broader goal of applying it to applications such as the study of wound healing, pathological hypertrophy and atrophy, as well as drug efficacy and interaction.References
- [Ambrosi and Mollica 2002]
- Ambrosi D, Mollica F (2002) On the mechanics of a growing tumor. Int J Engr Sci 40:1297-1316
- [Armero 1999]
- Armero F (1999) Formulation and finite element implementation of a multiplicative model of coupled poro-pplasticity at finite strains under fully-saturated conditions. Comp Methods in Applied Mech Engrg 171:205-241
- [Bischoff et al. 2002a]
- Bischoff JE, Arruda EM, Grosh K (2002a) A microstructurally based orthotropic hyperelastic constitutive law. J Applied Mechanics 69:570-579
- [Bischoff et al. 2002b]
- Bischoff JE, Arruda EM, Grosh K (2002b) Orthotropic elasticity in terms of an arbitrary molecular chain model. J Applied Mechanics 69:198-201
- [Brennen 1995]
- Brennen CE (1995) Cavitation and Bubble Dynamics. Oxford University Press
- [Bromberg and Dill 2002]
- Bromberg S, Dill KA (2002) Molecular Driving Forces: Statistical Thermodynamics in Chemistry and Biology. Garland
- [Brooks and Hughes 1982]
- Brooks A, Hughes T (1982) Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comp Methods in Applied Mech Engrg 32:199-259
- [Bustamante et al. 2003]
- Bustamante C, Bryant Z, Smith SB (2003) Ten years of tension: Single-molecule DNA mechanics. Nature 421:423-427
- [Calve et al. 2004]
- Calve S, Dennis R, Kosnik P, Baar K, Grosh K, Arruda E (2004) Engineering of functional tendon. Tissue Engineering 10:755-761
- [Cowin and Hegedus 1976]
- Cowin SC, Hegedus DH (1976) Bone remodeling I: A theory of adaptive elasticity. Journal of Elasticity 6:313-325
- [Epstein and Maugin 2000]
- Epstein M, Maugin GA (2000) Thermomechanics of volumetric growth in uniform bodies. International Journal of Plasticity 16:951-978
- [Fung 1993]
- Fung YC (1993) Biomechanics: Mechanical properties of living tissues, 2nd edn. Springer-Verlag, New York
- [Garikipati and Rao 2001]
- Garikipati K, Rao VS (2001) Recent advances in models for thermal oxidation of silicon. Journal of Computational Physics 174:138-170
- [Garikipati et al. 2004]
- Garikipati K, Arruda EM, Grosh K, Narayanan H, Calve S (2004) A continuum treatment of growth in biological tissue: Mass transport coupled with mechanics. Journal of Mechanics and Physics of Solids 52:1595-1625
- [Graham et al. 2004]
- Graham JS, Vomund AN, Phillips CL, Grandbois M (2004) Structural changes in human type I collagen fibrils investigated by force spectroscopy. Experimental Cell Research 299:335-342
- [Han et al. 2000]
- Han S, Gemmell SJ, Helmer KG, Grigg P, Wellen JW, Hoffman AH, Sotak CH (2000) Changes in ADC caused by tensile loading of rabbit achilles tendon: Evidence for water transport. Journal of Magnetic Resonance 144:217-227
- [Harrigan and Hamilton 1993]
- Harrigan TP, Hamilton JJ (1993) Finite element simulation of adaptive bone remodelling: A stability criterion and a time stepping method. Int J Numer Methods Engrg 36:837-854
- [Hughes et al. 1987]
- Hughes T, Franca L, Mallet M (1987) A new finite element formulation for computational fluid dynamics: VII. Convergence analysis of the generalized SUPG formulation for linear time-dependent multidimensional advective-diffusive systems. Comp Methods in Applied Mech Engrg 63(1):97-112
- [Hughes 1987]
- Hughes TJR (1987) Recent progress in the development and understanding of SUPG methods with special reference to the compressible Euler and Navier-Stokes equations. International Journal for Numerical Methods in Fluids 7:1261-1275
- [Kratky and Porod 1949]
- Kratky O, Porod G (1949) Röntgenuntersuchungen gelöster Fadenmoleküle. Recueil Trav Chim 68:1106-1122
- [Kuhl and Steinmann 2003]
- Kuhl E, Steinmann P (2003) Theory and numerics of geometrically-nonlinear open system mechanics. Int J Numer Methods Engrg 58:1593-1615
- [Kuhl et al. 2005]
- Kuhl E, Garikipati K, Arruda E, Grosh K (2005) Remodeling of biological tissue: Mechanically induced reorientation of a transversely isotropic chain network. Journal of the Mechanics and Physics of Solids 53(7):1552 - 73
- [Landau and Lifshitz 1951]
- Landau LD, Lifshitz EM (1951) A Course on Theoretical Physics, Volume 5, Statistical Physics, Part I. Butterworth Heinemann (reprint)
- [Lee 1969]
- Lee EH (1969) Elastic-Plastic Deformation at Finite Strains. J Applied Mechanics 36:1-6
- [Marko and Siggia 1995]
- Marko JF, Siggia ED (1995) Stretching DNA. Macromolecules 28:8759-8770
- [Nordin et al. 2001]
- Nordin M, Lorenz T, Campello M (2001) Biomechanics of tendons and ligaments. In: Nordin M, Frankel VH (eds) Basic Biomechanics of the Musculoskeletal System, Lippincott Williams and Wilkins, N.Y., pp 102-125
- [Provenzano et al. 2003]
- Provenzano PP, Martinez DA, Grindeland RE, Dwyver KW, Turner J, Vailas AC, Vanderby R (2003) Hindlimb unloading alters ligament healing. Journal of Applied Physiology 94:314-324
- [Rief et al. 1997]
- Rief M, Oesterhelt F, Heymann B, Gaub HE (1997) Single Molecule Force Spectroscopy o Polysaccharides by Atomic Force Microscopy. Science 275:1295-1297
- [Sengers et al. 2004]
- Sengers BG, Oomens CWJ, Baaijens FPT (2004) An integrated finite-element approach to mechanics, transport and biosynthesis in tissue engineering. J Bio Mech Engrg 126:82-91
- [Simo et al. 1985]
- Simo JC, Taylor RL, Pister KS (1985) Variational and projection methods for the volume constraint in finite deformation elasto-plasticity. Comp Methods in Applied Mech Engrg 51:177-208
- [Sun et al. 2002]
- Sun YL, Luo ZP, Fertala A, An KN (2002) Direct quantification of the flexibility of type I collagen monomer. Biochemical and Biophysical Research Communications 295:382-386
- [Swartz et al. 1999]
- Swartz M, Kaipainen A, Netti PE, Brekken C, Boucher Y, Grodzinsky AJ, Jain RK (1999) Mechanics of interstitial-lymphatic fluid transport: Theoretical foundation and experimental validation. J Bio Mech 32:1297-1307
- [Taber and Humphrey 2001]
- Taber LA, Humphrey JD (2001) Stress-modulated growth, residual stress and vascular heterogeneity. J Bio Mech Engrg 123:528-535
- [Taylor 1999]
- Taylor RL (1999) FEAP - A Finite Element Analysis Program. University of California at Berkeley, Berkeley, CA
- [Tezduyar and Sathe 2003]
- Tezduyar T, Sathe S (2003) Stabilization parameters in SUPG and PSPG formulations. Journal of Computational and Applied Mechanics 4:71-88
- [Truesdell and Noll 1965]
- Truesdell C, Noll W (1965) The Non-linear Field Theories (Handbuch der Physik, band III). Springer, Berlin
Footnotes:
At this point, we do not distinguish the solid species further. This is a good approximation to the physiological setting for tendons, which are relatively acellular and whose dry mass consists of up to 75% collagen [Nordin et al., 2001]. Currently, we do not consider physiological processes which involve migration of the cells or surrounding matrix within the tissue-such as the migration of fibroblasts within the extra-cellular matrix during wound healing. The amino acids, nutrients and regulators are in solution at low concentrations, and do not bear any appreciable stress. Under isothermal conditions, the only contribution to is from the strain energy. Where we are referring to the fluid being subject to pure tension, not a reduction in fluid compressive stress from an initial ambient pressure, manifesting itself as tensile stress. The argument above works for all cases other than the one where the system starts off unsaturated and the applied deformation has determinant , necessitating a collapse of voids. The fix is probably right here. Something like . and not in terms of acceleration relative to fixed stars for e.g., as discussed in Truesdell and Noll, [1965]. where every observer has an implicit knowledge of the directionality of the field relative to a fixed frame, allowing it to transform objectively.Previous Research
(Will be populated when I have the time.)