abarret / clotting_model Goto Github PK

We need an implementation of a CartGridFunction that evaluates the transition rates from unactivated to activated.

We need to calculate the force balance on the wall points.

This involves two things:

1. Moving the wall points with velocity given from the force balance.
2. Applying a force onto the fluid from the force balance.

Note that in some cases, our force balance equations do not make sense. In those instances, the wall points and wall points velocity should be the same as the mesh.

Below is a list of things needed to implement the full bound platelet model.

• Bond formations (#47)
• Bound / Activated platelet sources (#44)
• Cohesion stress formation / relaxation (#45)
• Wall site consumption (#49)
• Implementation in main routines (#49)
• Bound velocity (#42)
• Drag force (not needed because it balances with div(sigma))
• Conversion function for cohesion stress tensor (#48)

It would be useful to have a function that computes the free wall sites (w_max - w). What makes this challenging is the w is computed from W which is a Lagrangian state variable. We need to track the free wall sites on the Lagrangian structures and then spread that to the Eulerian grid.

We need to create an extensions of CFRelaxationOperator that computes the relaxation and stress build up. An initial implementation is given in CohesionStressRHS. We need to verify that this is up-to-date.

We need an implementation of CartGridFunction that computes the bond formation/breakage rates. Note that this should look similar to #10, although because the bond equation is a scalar equation, we don't need to evolve with CFINSForcing.

An implementation of CartGridFunction that evaluates the transition from unbound-activated to bound and bound to unbound-activated. This in part involves a convolution function, which is given in #4.

The bound platelet's move at velocity u_b which is different from the fluid velocity inside the thrombus. We need an implementation of CartGridFunction that evaluates u_b inside the thrombus.

Note that this velocity field is not divergence free so we will need to use an appropriate convective operator for evolving phi_b.