Modelling growth and formation of thrombi: a multiphasic approach based on the theory of porous media

Aortic dissection (AD) has a high mortality rate. About 40% of the people with type B AD do not live for more than a month. The prognosis of AD is quite challenging. Hence, we present a triphasic model for the formation and growth of thrombi using the theory of porous media (TPM). The whole aggregate is divided into solid, liquid and nutrient constituents. The constituents are assumed to be materially incompressible and isothermal, and the whole aggregate is assumed to be fully saturated. Darcy’s law describes the flow of fluid in the porous media. The regions with thrombi formation are determined using the solid volume fraction. The velocity- and nutrient concentration-induced mass exchange is defined between the nutrient and solid phases. We introduce the set of equations and a numerical example for thrombosis in type B AD. Here we study the effects of different material parameters and boundary conditions. We choose the values that give meaningful results and present the model’s features in agreement with the Virchow triad. The simulations show that the thrombus grows in the low-velocity regions of the blood. We use a realistic 2-d geometry of the false lumen and present the model’s usefulness in actual cases. The proposed model provides a reasonable approach for the numerical simulation of thrombosis.


Introduction
The aorta is one of the essential arteries in the body.The heart pumps the blood from the left ventricle into the aorta via the aortic valve, which opens and closes with each heartbeat to allow a one-way blood flow.Aortic dissection (AD) begins when a tear occurs in the inner layer (intima) of the aortic wall.This tear allows the blood to flow between the inner and middle layers causing them to separate (dissect).This second blood-filled channel is called a false lumen, where thrombosis (blood clotting) occurs, cf.Fig. 1.Blood clotting or coagulation is the process which prevents excessive bleeding by forming a spatial structure called a thrombus.A thrombus consists of small blood cells (platelets) and fibrous protein (fibrin) that stop the bleeding at the injury site.There are two types of AD depending on the location of the dissection.In type A AD, the dissection happens in the ascending part of the aorta, where the expansion of the false lumen can push other aorta branches and reduce blood flow.In contrast, the dissection occurs in the descending part of the aorta in type B AD, which may extend into the abdomen.The formation of a thrombus involves a complex sequence of biochemical reactions Fig. 1 Illustrations of the true and false lumen in type B AD with entry and exit tears (left), and the initial formation of false lumen (without a thrombus) and formation of a thrombus in false lumen (right) [32] [ 1,2].AD can occur due to high blood pressure leading to increased stress on the aortic wall, weakening of the wall, pre-existing aneurysm or defects in the aortic valve, to name a few factors.Approximately 75% of type B AD patients have hypertension [3,4].Furthermore, Virchow's triad describes three physiological factors that can result in thrombosis.These factors are endothelial injury, hypercoagulability of blood and stasis of blood flow [5,6].
AD is a highly fatal disease.The estimated occurrence is 5-30 cases per million people annually.Among these, the acute cases are 2-−3.5 cases per 100,000 people per year, accounting for 6000 to 10,000 cases per year in the USA alone.To understand the gravity of AD, approximately 75% of patients with a ruptured aortic aneurysm make it to the emergency department alive.However, 40% of AD patients die immediately.The mortality rate for aortic dissection is high, especially in acute cases.The mortality rate is excessive in the first seven days after type B AD due to severe complications, such as malperfusion or rupture in the aorta.In the case of such difficulties, a typical open surgery involves a 14-67% risk of irreversible damage to the spinal cord or, worse, mortality.Looking at the long-term prognosis, the survival rate of patients is 50-80% for five years and 30-60% for ten years [7][8][9][10].The short-term and long-term diagnosis for AD remains unclear, leading to an interest in computational methods to help with the decision-making process for the treatment.We focus here on the formation and growth of a thrombus.
For modelling a thrombus, one has to consider its multiphasic structure and the associated characteristics of the constituents.A highly complicated microscopic model would prevent us from establishing a usable computational model.Therefore, we use a macroscopic continuum mechanical approach of the theory of porous media (TPM).TPM provides an excellent framework to describe the multiphasic microstructure of a thrombus [11,12].It allows to describe macroscopically the complicated microstructure of biological tissues which is almost impossible to determine quantitatively.TPM was developed by combining the theory of mixtures, which was developed using the framework of general thermodynamical considerations, with the concept of volume fractions [13][14][15][16][17].It was continuously improved and developed to the current understanding of the TPM by de Boer & Ehlers [18], and Ehlers [11,19].De Boer has presented an excellent insight into the historical development of TPM in his book [12].Numerous models based on TPM have been presented to describe the behaviour, growth and remodelling of the soft tissues, brain tissue, liver perfusion, intervertebral disc, bone remodelling and tumour growth [20][21][22][23][24][25][26][27][28].The recent works using TPM for modelling biological growth motivate the use of this approach in the presented work.
Chemical, mechanical and metabolic factors drive the growth process of the thrombus.Because of the multiphasic nature of the thrombus, we present a triphasic model consisting of solid, liquid and nutrient phases.Due to the need for more detailed knowledge and parameters to quantify the influence of different factors, the model description is quite challenging.However, the effects of the blood velocity and the nutrients on the growth of thrombi are well researched [29][30][31].Therefore, we present a velocity-and nutrient concentrationinduced growth model based on the theory of porous media.We treat the highly coupled set of differential equations within the framework of the standard Galerkin procedure and implement the weak forms in the nonlinear finite element solver PANDAS.

Theory of porous media
The theory of porous media provides an excellent framework to macroscopically describe the complicated microstructure of the thrombus without knowing its detailed geometry.Therefore, a representative elementary volume (REV) is locally defined, where the individual constituents are considered to be in a state of ideal disarrangement.Using the real or virtual averaging processes over the REV, the microscale information of the overall aggregate and its constituents is homogenised to macro-scale quantities.For the investigated porous body, the immiscible parts lead to a triphasic aggregate ϕ consisting of solid ϕ S , which is saturated by fluid ϕ F .The fluid itself consists of liquid ϕ L and nutrients ϕ N , cf. Figure 2.
The volume fractions n α of the constituents ϕ α , where α ∈ {S, L , N }, are defined as the local ratios of the respective partial volume elements d v α with respect to the bulk volume element d v of the overall aggregate ϕ as [12] where x is the position vector in the current configuration at time t.The volume fractions n α need to permanently fulfil the saturation constraint (1) 2 .In the case of injury (aortic dissection), the nutrients are available in large quantities close to the injury site (false lumen).Therefore, in the presented monograph, the nutrients are described by their volume fractions n α .For small concentrations of nutrients, the approach of molar concentrations should be used [11,28].Moreover, the partial density ρ α = d m α / d v of a constituent ϕ α can be related to its real density ρ α R = d m α / d v α via its volume fraction n α (1) 2 .Due to the volume fraction concept, all geometric and physical quantities, such as motion, deformation and stress, are defined in the total control space.Hence, they can be interpreted as the statistical average values of the real quantities.The overall aggregate body B is defined as the connected manifold of material points P α .At any time t, material points P α of all the constituents ϕ α simultaneously occupy each spatial point x of the current configuration.These particles proceed from different reference positions X α at time t = t o , which leads to individual motion, velocity and acceleration fields for each constituent Moreover, a unique inverse motion function χ χ χ −1 α needs to exist for the motion function χ χ χ α to be unique.The necessary and sufficient condition for this is the existence of non-singular Jacobian J α where det(•) denotes the determinant operator.Following the equations (2) 1 and (3) 1 , any physical quantity can be represented as either Lagrangean (material) or Eulerian (spatial) description.Moreover, the material deformation gradient F α and its inverse F −1 α are defined as During deformation, the Jacobian J α is restricted to J α = det F α > 0. For scalar field functions , the material time derivative is defined as Furthermore, the balance equations for porous media are taken from the balance equations of the constituents ϕ α in mixture theory.The local balance equations of mass for the constituents ϕ α read as which gives volume balance using (1) as The local balance of momentum is given by and the local balance of moment of momentum excluding additional supply term is In Eqs. ( 6) and ( 8), div(•) denotes the spatial divergence operator, T α is the partial Cauchy stress tensor and b is the external volume force per unit mass.ρα represents the total mass production accounting for mass exchange or phase transitions between the constituents ϕ α .Total momentum production ŝα = pα + ρα x α contains the direct momentum exchange pα resulting from the interaction force between the constituents ϕ α as well as indirect parts resulting from the mass exchange ρα .The total production terms are restricted by

Assumptions
The system is investigated under the condition that all the constitutes ϕ α are materially incompressible, i.e.
This leads to the conclusion that volumetric deformations are only a result of a change in volume fractions n α .Moreover, the nutrient and the liquid phases are assumed to be in the fluid phase.For simplification, both phases are assigned the same velocity We assume that the liquid phase is not involved in the mass exchange.Using this assumption and Eq. ( 10) 1 , we get Furthermore, only isothermal processes are considered, energy transfer due to chemical reactions is neglected, accelerations are excluded, and the internal structure of the thrombus is considered to be isotropic.
Along with the assumptions, the Clausius-Planck inequality is written as where d α is the symmetric part of spatial velocity gradient.This equation is essential for developing a thermodynamically consistent model.

Constitutive modelling
The dependencies of the Helmholtz free energy for the solid, liquid and nutrient phases are considered as where C S = F T S F S is the right Cauchy-Green tensor related to solid.Furthermore, the universal dissipation principle must be satisfied by the constitutive relations.Therefore, we evaluate the Clausius-Planck inequality by following the procedure of Coleman and Noll [33].For more information regarding the evaluation of the entropy principle, the reader is referred to Bowen [34] and Ehlers [11].
Initially, an additional saturation constraint (1) 2 is added to the entropy inequality to ensure the fully saturated condition in the overall aggregate at any given time.This is done by introducing a Lagrangean multiplier p as a weight to the saturation condition as [12] p Moreover, we multiply the volume balance of the individual constituent ϕ α with the respective Lagrangean multipliers p α and make use of the relation l l l α : I = d α : Now using the relations ( 16)-( 18), (12) and the summation assumption (10) 1,2 , we can evaluate the inequality (14), following the methodology from Ricken and Bluhm [35].
which must hold for fixed values of the process variables and arbitrary values of freely available quantities d α , (n α ) α [35].Therefore, we obtain the following structure of the entropy inequality In this context, we obtain the necessary and sufficient thermodynamic restrictions.The parts concerning (n α ) α give Furthermore, we obtain the relations for the partial Cauchy stresses using ( 19) and ( 21) We can now introduce the chemical potentials α using ( 19) and ( 21) as For more information on chemical potentials, the reader is referred to Bowen [34].Furthermore, we make use of the assumption ρ L = 0 which gives us the dissipative part as where pF E is called the extra momentum production [36].Moreover, this gives the restrictions for the solid mass production ρS and momentum production pF as postulated by Ricken and Bluhm [35] ρS = δ N S ( N − S ), where δ N S ≥ 0, where S F is the permeability tensor between the solid and fluid phases.These restrictions give us the possibility to further formulate stresses, mass production and interaction forces.

Stress
The stress relations can be introduced following the restrictions from the entropy inequality (22) and neglecting the effective or frictional fluid stress, i.e.T F E ≈ 0. It is assumed that the fluid extra stress is much smaller in comparison with pF E also known as the effective drag force.Therefore, assuming ∂ψ F /∂n F = 0 where ϕ F = ϕ L + ϕ N yields [11,35] where I is the second-order identity tensor.This also leads to p N = p L = p using (21).Furthermore, p is identified as the unspecified pore pressure.The total stress is defined as the sum of partial stresses and using Eq.(1) 2 yields: According to the principle of material objectivity, the constitutive equations should not depend on the observer's position.Thus, the mathematical interpretation of such an objectivity condition states that constitutive equations must be invariant under rigid body rotations of the actual configuration [37,38].Therefore, the free energy functions will be formulated with a dependence on the principal invariants I 1 , I 2 and I 3 of C S using [39] where tr(•) denotes the trace operator.Thereon, the Helmholtz free energy function can be written as follows The structure of the equation satisfies the invariance and polyconvexity conditions, which also implies quasiconvexity.That would ensure the existence of minimisers of the related variational principles in finite elasticity.
For a detailed discussion on convexity conditions, the reader is referred to Dacorogna [40].Now the Helmholtz free energy function can be constructed in the following way where the term in the front with nutrient volume fraction n S represents the change in solid rigidity with respect to the initial volume fraction n S O S .Here (•) O S represents the initial value of (•) with respect to the referential configuration of the solid.This term in the front with solid volume fraction n S accounts for the change in solid rigidity with respect to the initial volume fraction n S O S .The material parameter n helps to define the extent of dependence on change in porosity.Carter and Hayes [41] identified material parameter n = 3, and it has been used in multiple porous media growth applications [24][25][26]35].The rest of the part ψ S neo is the Neo-Hookean material law.μ S and λ S are the macroscopic Lam é constants.
From ( 30) and ( 26), the effective solid Cauchy stress can be obtained μ S (B S − I) + λ S (ln J S )I , (31) where B S is the left Cauchy-Green tensor B S = F S F T S .With equations ( 26), ( 30) and ( 31), total solid Cauchy stress can be written as The effective solid Kirchhoff stress and total solid Kirchhoff stress then read as

Filter Velocity
The seepage velocity w F S = x F − x S determines the motion of the fluid in relation to the solid.The relation pF = pL + pN is considered.From the evaluation of entropy inequality (25) 2 , we obtain the following relation as postulated by Ricken and Bluhm [35] Using ( 34), (8), S F = α F S I for an isotropic material and rearranging the equation, we get [35] n The material parameter α F S can be described either by using initial Darcy's permeability of fluid k F O S [m/s] and effective fluid weight γ F R [N /m 3 ] or by using initial intrinsic permeability of solid K S O S [m 2 ] and dynamic fluid viscosity where m is a dimensionless parameter which accounts for the change of permeability [12,26].Here, (n F /n F 0 S ) m K S 0 S = K S E is referred to as effective permeability, which takes into consideration the change in volume fractions.Eipper [42] proposed this porosity-dependent effective permeability where m ≥ 0. With Eqs. ( 13), ( 36) and ( 35), we finally get the following relation for seepage velocity

Mass exchange
According to (13), the mass exchange occurs between the solid and nutrient phases ρS = − ρ N .Following the evaluation of the entropy inequality (25) 1 , we have ρS ≥ 0. Furthermore, making use of the postulations proposed by Ricken et al. [25,26], and due to no expert knowledge available for the formulation of free energy functions of liquid and nutrient phases, the mass production for the solid phase is formulated.The effects of the blood velocity and the nutrients on the thrombus growth are well researched [29][30][31].Therefore, ρS [kg/m 3 s] is postulated as a function of w F S and n where C represents the maximum mass exchange, and β 1 and β 2 are the material parameters reflecting the dependence of mass exchange on the seepage velocity and nutrient volume fraction, respectively, cf.Fig. 3.

Numerical treatment
Considering the assumptions, balance equations and constitutive relations from the preceding sections, we have a set of six independent variables where u S is the displacement of the solid phase.Using Darcy's formulation for the seepage velocity w F S (35), the set of unknowns could be decreased to five.Furthermore, the saturation condition (1) 2 (n L = 1 − n S − n N ) reduces the set of unknowns to four Once this is concluded, the weak formulation for the governing equations is formulated in the framework of the standard Galerkin procedure (Bubnov-Galerkin).This is achieved by multiplying the momentum balance of the mixture, volume balance of the mixture, solid and nutrients with the test functions δu S , δp, δn S and δn N , respectively.As a result, the weak formulation of the triphasic model reads • Momentum balance of mixture: • Volume balance of mixture: • Volume balance of solid: • Volume balance of nutrients: In the weak formulation, from ( 41) to (44), t is the external load vector acting on the Neumann boundary t , n F w F S • n is the fluid mass efflux on the Neumann boundary q and n N w F S • n is the nutrient mass efflux on the Neumann boundary υ , where n is the outward oriented unit surface normal.Also, the transport equation for nutrients consists of the general structure of an advection equation.It is known that using the given approach for formulating weak forms, equation ( 44) generates large oscillations if not properly stabilised or if the mesh size is not excessively small.Therefore, an artificial diffusion term (r ) is added only to the volume balance of nutrients (44) to stabilise the transport equation.For more information on the mass transport equation and the stabilisation schemes, the reader is referred to Santos et al. [43] and the references therein.

Numerical example
In this section, we present a numerical example of the formation and growth of a thrombus in type B aortic dissection using a realistic two-dimensional geometry of false lumen.The constitutive equations for the solid stress T S , the mass production term ρS and the seepage velocity w F S provide the thrombus-specific material laws.In addition, the coupled set of governing equations presents the capabilities of the model.We implement the weak forms of balance equations in the FE package PANDAS, where time and space adaptive methods are widely applied.However, modelling the growth of living tissues has its challenges.In the case of thrombosis, it is not straightforward to obtain the data and perform experiments on living tissues.Therefore, the parameters are chosen that give reasonable results for specific cases.
In the numerical example, we want to use a 2-d cross section of a realistic geometry of the false lumen to model thrombosis.To obtain the cross section, we use a 3-d model of an aorta consisting of a true and false lumen [44], cf.Fig. 4 (left).We cut this model in the x − y plane represented as the cutting plane in Fig. 4 (top right).This gives us the 2-d geometry in the x − y plane, cf.Fig. 4 (bottom right).However, because we are modelling thrombosis in the false lumen, the geometry of the false lumen is of interest to us.Therefore, we create and discretise the false lumen's geometry using CUBIT, consisting of 1063 elements, cf.Fig. 5.The geometry consists of a solid matrix saturated with fluid, consisting of liquid and nutrient phases.Furthermore, we arbitrarily choose the entry and exit tears position.The boundary conditions at the entry tear include n S = 0.4 and q = 0.1 m/s.The bottom exit tear (dotted line) is the drained surface with p = 0.0 N/m 2 .The rest of the boundaries are undrained surfaces, cf.Fig. 5. Also, the right side of the false lumen is fixed in both the x and y directions.We use the Taylor Hood elements for spatial discretisation, where quadratic approximation is used for the solid displacements u S and linear approximation for the pressure p, solid volume fraction n S and nutrient volume fraction n N .The simulation is performed using the parameters given in Table 1 and with a time step size of 100 s.Here, we assume zero Poisson's ration (λ S = 0) for the solid skeleton.This is done for simplicity and lack of material data [35,[45][46][47].
We can now discuss the results while drawing an analogy with the process of thrombosis.At time t = 0 hours, the initial solid volume fraction n S is 0.2 because of the presence of subendothelial collagen, wall cells and activated platelets on the formation of the false lumen.As the fluid enters via the entry tear, it creates highand low-velocity regions, cf.Fig. 6.Because of different velocity profiles and the availability of nutrients, the process of thrombosis begins, and the solid volume fractions start increasing, cf.Fig. 7.This can be compared to primary haemostasis, where the platelets accumulate at the injury site and form a platelet plug.The solid volume fraction n S increases further due to the mass exchange rate (38) 3 dependence on seepage velocity w F S and nutrient volume fraction n N , which can be compared to secondary haemostasis.During secondary haemostasis, the clotting factors interact in a complicated series of chemical reactions leading to the formation of fibrin fibre.The platelets and the fibrin fibre form a mesh leading to the development of a stable plug.This process continues further to form a permanent solid plug called a thrombus.Furthermore, Fig. 8 shows that the permeability in the false lumen decreases with an increase in the solid volume fraction.This means that as the thrombus grows and there is more solid, the ability of the fluid to move through porous media decreases.Moreover, we observe a singularity at the exit tear due to the sharp edges, which is a numerical artefact [48,49].
Further, we can explore the features of the model.We can vary the material parameters and the boundary conditions to adapt the model for specific cases.

Influence of material parameters in mass exchange
The material parameters, β 1 and β 2 , in Eq. ( 38) can be varied to change the dependence of mass exchange rate on the seepage velocity and nutrient volume fraction.In the above example, the default values of 0.05 and 5.0 are used for β 1 and β 2 , respectively.To begin with, we vary the material parameter β 1 leading to a change in From this, we can see that as we increase the value of β 1 , forming a thrombus is easier due to growth happening for the wider range of seepage velocities, cf.Fig. 9. Furthermore, we vary the values of the material parameter β 2 and change the influence of nutrient volume fraction on growth.As the value of β 2 increases, we allow the thrombus growth for a more extensive range of nutrient volume fraction n N , leading to higher growth because of the increased availability of nutrients, cf.Fig. 10.

Change in Neumann boundary condition
Moreover, we can see the effects due to variation in the fluid mass influx q on the Neumann boundary q .An increase in the fluid mass influx q at the entry tear results in higher seepage velocity.Because of the mass exchange dependence on the seepage velocity, forming a thrombus in the false lumen, especially in the middle section, is difficult, cf.Fig. 11.This also fits well with the physiological understanding of thrombosis and the Virchow triad, where it is difficult to form blood clots when the blood velocity is high [5].
The above analysis shows that the triphasic model accommodates the well-known Virchow triad, which describes three physiological factors that can result in thrombosis.The first one, endothelial injury, is included in the form of the presence of a false lumen.On the formation of a false lumen, the endothelium is damaged, which lines the inner layer of the blood vessels.The endothelial injury stimulates the platelets and coagulation process.The second factor is hypercoagulability, which is an increased tendency of coagulation in the body due to inherited or acquired disorders.The material parameter β 2 and C can be used to adapt the model for such a scenario.Also, the material parameter β 2 can be used to include this increased tendency of coagulation.The third factor, the stasis of blood, is present in the form of mass exchange dependency on the seepage velocity [5,50].Here, the material parameter β 1 can be used to adapt the model for the specific case.Finally, the fluid mass efflux q boundary condition can be used to incorporate the factor of high blood pressure, which is known to be the major cause of aortic dissection.

Conclusions
A triphasic model has been developed for the thrombus's growth, capable of describing its growth and nonlinear behaviour.The theory of porous media (TPM) provides an excellent framework to consider the multiphasic nature of a thrombus.Therefore, we used TPM to develop a thermodynamically consistent model using a smeared model of solid and nutrient-rich liquid phases.The constitutive relations are proposed based on the (a) q = 0.05 (b) q = 0.1 (c) q = 0.5 Fig. 11 Change in solid volume fractions n S for different values of q at t = 83 hr restrictions obtained by evaluating the entropy inequality.The mass exchange between the solid and nutrient phases is formulated, which depends on the nutrient concentration and the seepage velocity.The simulations show the growth of the thrombus in low-blood-velocity regions.Here, we see that the material parameters play an important role in incorporating the physiological factors that can result in thrombosis, according to the Virchow triad.We have an additional parameter to incorporate the factor of inherited or acquired disorders leading to hypercoagulability.The stasis of blood and endothelial injury, along with hypertension, are also incorporated into the model.Moreover, the triphasic model gives the advantage of including the mass exchange between the nutrient and solid phases without altering the amount of liquid which fits well with the physiological understanding of thrombosis.The model also proves its usefulness in actual cases.
However, biological modelling is challenging.There is a need to quantify the different factors responsible for thrombi growth and determine the material parameters from additional experiments.This would also give us more accurate boundary conditions and possibilities to guide the growth process in a much more accurate way.This would also lead to the potential splitting of the mass production formulation into two approaches: primary and secondary haemostasis.These two processes consist of the major part of the thrombus formation.Moreover, the non-Newtonian behaviour of blood should be considered in the model.Also, there is a lack of availability of enough medical data.With enough medical data (CT/MRI scans) taken at different stages of thrombosis, it would be possible to validate the model with further improvements.
Moreover, the model can be extended to different cases of AD, such as type A AD or false lumen with multiple tears.A larger set of medical and experimental data can aid in developing, validating and training such models and performing patient-focused simulations.Additionally, because the short-term and long-term diagnosis of AD, especially type B AD, is unclear, all the models for different AD cases can be combined to develop a numerical laboratory and help in decision-making.Furthermore, the model of thrombosis could be extended to conditions where the formation of the blood clot is critical, such as deep venous thrombosis, hypercoagulability disorders and disseminated intravascular coagulation (DIC).Understanding the mechanics of growth in such chronic conditions can open new directions in medical device design, personalised medicine, prognosis and controlling disease progression [51].

Fig. 2
Fig. 2 Illustration of the microstructure of the porous false lumen (left), macro-model obtained by volumetric homogenisation process (centre) and superimposed continua (right)

Fig. 3
Fig. 3 Mass productions ρS w F S and ρS n N dependence on the seepage velocity w F S and nutrient volume fraction n N , respectively

Fig. 4 3
Fig.4 3-d model of the aorta with a true and false lumen (left)[44].Cutting plane and the resulting 2-d cross section of the lumens (right)

Fig. 5
Fig. 5 Discretisation and boundary conditions of the false lumen geometry

wFig. 6 Fig. 7
Fig. 6 Norm of the seepage velocity w F S at different stages in time

Fig. 8 1 Fig. 9
Fig. 8 Effective permeability K S E at different stages in time

Fig. 10
Fig. 10 Change in solid volume fractions n S for different values of β 2 at t = 83 hr

Table 1
Parameters for thrombi growth