Numerical accuracy comparison of two boundary conditions commonly used to approximate shear stress distributions in tissue engineering scaffolds cultured under flow perfusion

Abstract

Introduction Flow‐induced shear stresses have been found to be a stimulatory factor in pre‐osteoblastic cells seeded in 3D porous scaffolds and cultured under continuous flow perfusion. However, due to the complex internal structure of the scaffolds, whole scaffold calculations of the local shear forces are computationally intensive. Instead, representative volume elements (RVEs), which are obtained by extracting smaller portions of the scaffold, are commonly used in literature without a numerical accuracy standard. Objective Hence, the goal of this study is to examine how closely the whole scaffold simulations are approximated by the two types of boundary conditions used to enable the RVEs: “wall boundary condition” (WBC) and “periodic boundary condition” (PBC). Method To that end, lattice Boltzmann method fluid dynamics simulations were used to model the surface shear stresses in 3D scaffold reconstructions, obtained from high‐resolution microcomputed tomography images. Results It was found that despite the RVEs being sufficiently larger than 6 times the scaffold pore size (which is the only accuracy guideline found in literature), the stresses were still significantly under‐predicted by both types of boundary conditions: between 20% and 80% average error, depending on the scaffold's porosity. Moreover, it was found that the error grew with higher porosity. This is likely due to the small pores dominating the flow field, and thereby negating the effects of the unrealistic boundary conditions, when the scaffold porosity is small. Finally, it was found that the PBC was always more accurate and computationally efficient than the WBC. Therefore, it is the recommended type of RVE.

In vivo measurement of blood clot mechanics from computational fluid dynamics based on intravital microscopy images

“…Moreover, it has recently gained popularity for modeling biological flows. [64][65][66][67][68][69][70][71] It was chosen for this work because of the ease with which it handles complex boundary conditions, and due to its efficient parallelizability on supercomputers. [72] Both of these aspects make the method especially attractive for dealing with high-resolution images, such as the ones used in our study.…”

Scaffolds with a High Surface Area-to-Volume Ratio and Cultured Under Fast Flow Perfusion Result in Optimal O2 Delivery to the Cells in Artificial Bone Tissues

“…A previously developed custom-written, in-house code was used in this work [16][17][18]20,26,29,31]. The D3Q15 lattice [32], in conjunction with the single-relaxation time Bhatnagar, Gross, and Krook [33] collision term approximation, was used to perform simulations.…”

Scaffolds with a High Surface Area-to-Volume Ratio and Cultured Under Fast Flow Perfusion Result in Optimal O₂ Delivery to the Cells in Artificial Bone Tissues

“…In addition, the LBM method has gained popularity within the scientific computing community because of the ease with which it can be parallelized on supercomputers [30]. A previously developed custom-written, in-house code was used in this work [16][17][18]20,26,29,31]. The D3Q15 lattice [32] in conjunction with the single-relaxation time Bhatnagar, Gross and Krook [33] collision term approximation was used to perform simulations.…”