Flow speed and RBC generation

[jounihuo]

I have modified the pipeflow example to have a simple time dependent flow speed amplitude.

setExternalVector(*hemocell.lattice, hemocell.lattice->getBoundingBox(), DESCRIPTOR<T>::ExternalField::forceBeginsAt, plb::Array<T, DESCRIPTOR<T>::d>(poiseuilleForce*amplitude, 0.0, 0.0));

But do I need to take steps to ensure RBC generation has similar amplitude (if yes then how?) or is it automatic?

[gzavo]

Hi @jounihuo,

If you have a periodic simulation (as in the pipeflow case), the number of cells are constant during the run, so you only need to reset the driving force in every iteration, the cell generation happens using the cellPack tool before running the simulation.
If you have a more complex domain, that is not periodic (e.g., a capillary network), then you can use the periodic pre-inlet to drive the domain. See the documentation here, and related publications describing it here and applying it here.