Permeability and fluid flow-induced wall shear stress of bone tissue scaffolds: Computational fluid dynamic analysis using Newtonian and non-Newtonian blood flow models

Enhancing precision in effective regurgitant orifice area estimation by transthoracic echocardiography for functional mitral regurgitation using computational fluid dynamics

Computational fluid dynamics (CFD) was used to identify factors influencing the accuracy of the hemispherical proximal isovelocity surface area (PISA) method in calculating the effective regurgitant orifice area (EROA) for patients with functional mitral regurgitation (FMR). Ninety-nine CFD models were constructed to investigate the impact of regurgitant orifice shape and leaflet tethering on the EROA calculation using the PISA method. The correction factors for regurgitation orifice shape (CFs) and for leaflet tethering (CFt) were derived by comparing the 2D PISA method and the actual orifice area. The correction formula was then tested in vivo via 2D transthoracic echocardiography with 3D transesophageal echocardiography of the vena contracta area (VCA) as a reference method in 62 patients with FMR. Based on the CFD simulation results, the two major factors for correcting the EROA calculation were vena contracta length (VCL) and coaptation depth (CD). The correction formula for the EROA was corrected effective regurgitant orifice area (CEROA) = EROA*CFs*CFt, where CFs = 0.59 × VCL(cm) + 0.6 × MR Vmax(cm/s)-0.63 × PISA R(cm)-1.51 and CFt = 0.4 × CD (cm) + 0.96. The correction formula was applied to FMR patients, and the bias and LOA between the CEROA and VCA (0.01 ± 0.13 cm2) were much smaller than those between the EROA and VCA (0.26 ± 0.32 cm2). The CFD-based correction formula improves the accuracy of the EROA calculation based on the hemispheric PISA method, possibly leading to more accurate and reliable data for treatment decision-making in FMR patients.

An explainable machine learning-based probabilistic framework for the design of scaffolds in bone tissue engineering

Recently, 3D-printed biodegradable scaffolds have shown great potential for bone repair in critical-size fractures. The differentiation of the cells on a scaffold is impacted among other factors by the surface deformation of the scaffold due to mechanical loading and the wall shear stresses imposed by the interstitial fluid flow. These factors are in turn significantly affected by the material properties, the geometry of the scaffold, as well as the loading and flow conditions. In this work, a numerical framework is proposed to study the influence of these factors on the expected osteochondral cell differentiation. The considered scaffold is rectangular with a 0/90 lay-down pattern and a four-layered strut made of polylactic acid with a 5% steel particle content. The distribution of the different types of cells on the scaffold surface is estimated through a scalar stimulus, calculated by using a mechanobioregulatory model. To reduce the simulation time for the computation of the stimulus, a probabilistic machine learning (ML)-based reduced-order model (ROM) is proposed. Then, a sensitivity analysis is performed using the Shapley additive explanations to examine the contribution of the various parameters to the framework stimulus predictions. In a final step, a multiobjective optimization procedure is implemented using genetic algorithms and the ROM, aiming to identify the material parameters and loading conditions that maximize the percentage of surface area populated by bone cells while minimizing the area corresponding to the other types of cells and the resorption condition. The results of the performed analysis highlight the potential of using ROMs for the scaffold design, by dramatically reducing the simulation time while enabling the efficient implementation of sensitivity analysis and optimization procedures.

Evaluating the effect of pore size for 3d-printed bone scaffolds

The present study investigated the influence of pore size of strut-based Diamond and surface-based Gyroid structures for their suitability as medical implants. Samples were made additively from laser powder bed fusion process with a relative density of 0.3 and pore sizes ranging from 300 to 1300 μm. They were subsequently examined for their manufacturability and mechanical properties. In addition, non-Newtonian computational fluid dynamics and discrete phase models were conducted to assess pressure drop and cell seeding efficiency. The results showed that both Diamond and Gyroid had higher as-built densities with smaller pore sizes. However, Gyroid demonstrated better manufacturability as its relative density was closer to the as-designed one. In addition, based on mechanical testing, the elastic modulus was largely unaffected by pore size, but post-yielding behaviors differed, especially in Diamond. High mechanical sensitivity in Diamond could be explained partly by Finite Element simulations, which revealed stress localization in Diamond and more uniform stress distribution in Gyroid. Furthermore, we defined the product of the normalized specific surface, normalized pressure drop, and cell seeding efficiency as the indicator of an optimal pore size, in which this factor identified an optimal pore size of approximately 500 μm for both Diamond and Gyroid. Besides, based on such criterion, Gyroid exhibited greater applicability as bone scaffolds. In summary, this study provides comprehensive assessment of the effect of pore size and demonstrates the efficient estimation of an in-silico framework for evaluating lattice structures as medical implants, which could be applied to other lattice architectures.

Mechanical and Computational Fluid Dynamic Models for Magnesium-Based Implants

Today, mechanical properties and fluid flow dynamic analysis are considered to be two of the most important steps in implant design for bone tissue engineering. The mechanical behavior is characterized by Young's modulus, which must have a value close to that of the human bone, while from the fluid dynamics point of view, the implant permeability and wall shear stress are two parameters directly linked to cell growth, adhesion, and proliferation. In this study, we proposed two simple geometries with a three-dimensional pore network dedicated to a manufacturing route based on a titanium wire waving procedure used as an intermediary step for Mg-based implant fabrication. Implant deformation under different static loads, von Mises stresses, and safety factors were investigated using finite element analysis. The implant permeability was computed based on Darcy's law following computational fluid dynamic simulations and, based on the pressure drop, was numerically estimated. It was concluded that both models exhibited a permeability close to the human trabecular bone and reduced wall shear stresses within the biological range. As a general finding, the proposed geometries could be useful in orthopedics for bone defect treatment based on numerical analyses because they mimic the trabecular bone properties.

Understanding of trabecular-cortical transition zone: Numerical and experimental assessment of multi-morphology scaffolds

Applications of additive manufacturing (AM) in tissue engineering develop rapidly. AM offers layer-by-layer creation of complex objects, developed to restore functionality of, or replace, damaged tissues. Porous 3D-printed functional gradient structures are of particular interest: their special architecture makes it possible to simulate the heterogeneity of the replaced tissue and, by continuously changing the mechanical properties, to avoid the concentration of stresses that can be caused by abrupt geometric changes. Such structures also allow combinations of different types of unit cells and a smooth transition between them, making design of personalised scaffolds with optimal parameters for the replacement of damaged host tissue at the interface between tissues possible. This paper presents the results of development of scaffold structures with gradients of porosity and multi-morphology using unit cells based on triply periodic minimal surfaces (TPMS). The mechanical behaviour of additively manufactured scaffold prototypes made of polylactide acid (PLA) was studied under compressive loading. Strain fields on their surface were captured using the Vic-3d Micro-DIC digital image correlation system and compared with those obtained with detailed numerical simulations, employing elastic-plastic properties of PLA, obtained in experiments. The effect of gradient parameters and unit-cell morphology on the stress distribution in scaffolds was analysed. A smooth gradient transition between cells with different morphologies was found to reduce the probability of structural failure under intense compressive loading. A good agreement between numerical results and experimental data was achieved, which justifies application of the developed approach to design of personalised bone scaffolds.

Orthopedic Scaffolds: Evaluation of Structural Strength and Permeability of Fluid Flow via an Open Cell Neovius Structure for Bone Tissue Engineering

The ability of bone to regenerate itself through mechanobiological responses is its dynamic property. Mechanical cues from a neighboring environment produce the structural strain to promote blood flow and bone marrow mobility that in turn aids the bone regeneration process. Occurrences of these phenomena are crucial for the success of metallic scaffolds implanted in the host bone tissue. Thus, permeability and fluid flow-induced wall shear stress (WSS) are two parameters that directly influence cell bioactivities inside a scaffold and are crucial for effective bone tissue regeneration. Given that the scaffolds shall be implanted in the body, permeability assessment was carried out using non-Newtonian fluid. In this work, the triply periodic minimal surface scaffolds with Neovius architectures were fabricated by using selective laser melting technology. The estimation of fluid flow was carried out using computational fluid dynamics (CFD) analysis with a non-Newtonian blood fluid model. Further, the structural strength of various open cell Neovius lattices was evaluated using a static compression test, and in vitro cell culture using Alamar blue assay was evaluated. Results revealed that the values of intrinsic blood flow permeability of the three-dimensional (3D)-printed open cell porous scaffold with Neovius architecture were of the same order of magnitude as those of human bone, ranging from 0.0025 × 10-9 to 0.0152 × 10-9 m2. The structural elastic modulus and compressive strength of NOCL40, NOCL50, and NOCL60 lattices range from 3.27 to 3.71 GPa and 194 to 205 MPa, respectively. All of the values are comparable to the human bone, thus making these lattices a suitable alternative for orthopedic applications.

Simulating the mechanical stimulation of cells on a porous hydrogel scaffold using an FSI model to predict cell differentiation

3D-structured hydrogel scaffolds are frequently used in tissue engineering applications as they can provide a supportive and biocompatible environment for the growth and regeneration of new tissue. Hydrogel scaffolds seeded with human mesenchymal stem cells (MSCs) can be mechanically stimulated in bioreactors to promote the formation of cartilage or bone tissue. Although in vitro and in vivo experiments are necessary to understand the biological response of cells and tissues to mechanical stimulation, in silico methods are cost-effective and powerful approaches that can support these experimental investigations. In this study, we simulated the fluid-structure interaction (FSI) to predict cell differentiation on the entire surface of a 3D-structured hydrogel scaffold seeded with cells due to dynamic compressive load stimulation. The computational FSI model made it possible to simultaneously investigate the influence of both mechanical deformation and flow of the culture medium on the cells on the scaffold surface during stimulation. The transient one-way FSI model thus opens up significantly more possibilities for predicting cell differentiation in mechanically stimulated scaffolds than previous static microscale computational approaches used in mechanobiology. In a first parameter study, the impact of the amplitude of a sinusoidal compression ranging from 1% to 10% on the phenotype of cells seeded on a porous hydrogel scaffold was analyzed. The simulation results show that the number of cells differentiating into bone tissue gradually decreases with increasing compression amplitude, while differentiation into cartilage cells initially multiplied with increasing compression amplitude in the range of 2% up to 7% and then decreased. Fibrous cell differentiation was predicted from a compression of 5% and increased moderately up to a compression of 10%. At high compression amplitudes of 9% and 10%, negligible areas on the scaffold surface experienced high stimuli where no cell differentiation could occur. In summary, this study shows that simulation of the FSI system is a versatile approach in computational mechanobiology that can be used to study the effects of, for example, different scaffold designs and stimulation parameters on cell differentiation in mechanically stimulated 3D-structured scaffolds.

Computational Fluid Dynamic Analysis of customised 3D-printed bone scaffolds with different architectures

Through the recent years, tissue engineering has been proven as a solid substitute of autografts in the stimulation of bone tissue regeneration, through the development of three dimensional (3D) porous matrices, commonly known as scaffolds. In this work, we analysed two scaffold structures with 500μm pore size, by performing computational fluid dynamics simulations, to compare permeability, Wall Shear Stress (WSS), velocity and pressure distributions. Taking into account those parameters the geometry named as "PCL-50" was the best to anticipate showing a superior performance in supporting cell growth due to the improved flow characteristics in the scaffold.Clinical Relevance- Bone defects that require invasive surgical treatment with high risks in terms of success and effectiveness. Bone tissue engineering (BTE) in combination with the use of computational fluid dynamics (CFD) analysis tools aim to assist in designing optimal scaffolds that better promote bone growth and repair. The fluid dynamic characteristics of a porous scaffold plays a vital role in cell viability and cell growth, affecting the osteogenic performance of the scaffold.

Assessing the effect of manufacturing defects and non-Newtonian blood model on flow behaviors of additively manufactured Gyroid TPMS structures

In the field of medical engineering, Triply Periodic Minimal Surfaces (TPMS) structures have been studied widely owing to their physical attributes similar to those of human bones. Computational Fluid Dynamics (CFD) is often used to reveal the interaction between structural architectures and flow fields. Nevertheless, a comprehensive study on the effect of manufacturing defects and non-Newtonian behavior on the fluid responses in TPMS scaffolds is still lacking. Therefore, the present study fabricated Gyroid TPMS with four relative densities from 0.1 to 0.4. Non-destructive techniques were used to examine surface roughness and geometric deviation. We found that the manufacturing defects had a minor effect on fluid responses. The pressure drop comparison between defect-containing and defect-free models could be differed up to 7%. The same comparison for the average shear stress showed a difference up to 23%, in which greater deviation between both models was observed at higher relative density. On the contrary, the viscosity model played a significant role in flow prediction. By comparing the Newtonian model with Carreau-Yasuda non-Newtonian model, the resulting pressure drop and average wall shear stress from non-Newtonian viscosity could be higher than those of the Newtonian model by more than a factor of two. In addition, we matched the fluid-induced shear stress from both viscosity models with desirable ranges of shear stresses for tissue growth obtained from the literature. Up to 70% from the Newtonian model fell within the desirable range while the matching stress reduced to lower than 8% for the non-Newtonian results. Furthermore, by correlating geometric features with physical outputs, the geometric deviation was seen associated with surface curvature while the local shear stress revealed a strong correlation with inclination angle. Overall, the present work emphasized the importance of the viscosity model for CFD analysis of the scaffolds, especially when resulting fluid-induced wall shear stress is of interest. In addition, the geometric correlation has introduced the alternative consideration of structural architectures from local perspectives, which could assist the further comparison and optimization among different porous scaffolds in the future.

The Effect of Tortuosity on Permeability of Porous Scaffold

In designing porous scaffolds, permeability is essential to consider as a function of cell migration and bone tissue regeneration. Good permeability has been achieved by mimicking the complexity of natural cancellous bone. In this study, a porous scaffold was developed according to the morphological indices of cancellous bone (porosity, specific surface area, thickness, and tortuosity). The computational fluid dynamics method analyzes the fluid flow through the scaffold. The permeability values of natural cancellous bone and three types of scaffolds (cubic, octahedron pillar, and Schoen's gyroid) were compared. The results showed that the permeability of the Negative Schwarz Primitive (NSP) scaffold model was similar to that of natural cancellous bone, which was in the range of 2.0 × 10^-11 m² to 4.0 × 10^-10 m². In addition, it was observed that the tortuosity parameter significantly affected the scaffold's permeability and shear stress values. The tortuosity value of the NSP scaffold was in the range of 1.5-2.8. Therefore, tortuosity can be manipulated by changing the curvature of the surface scaffold radius to obtain a superior bone tissue engineering construction supporting cell migration and tissue regeneration. This parameter should be considered when making new scaffolds, such as our NSP. Such efforts will produce a scaffold architecturally and functionally close to the natural cancellous bone, as demonstrated in this study.

How getting twisted in scaffold design can promote bone regeneration: A fluid-structure interaction evaluation

Bone tissue engineering (BTE) uses engineering principles to repair large bone defects, which requires effective mass transport ability of scaffolds to support cellular activities during bone regeneration. Since the implanted BTE scaffolds keep deforming under physiological loading which influences the fluid flow and mass transport within the scaffold and surrounding tissue, thus, scaffold design needs to consider the mass transport behavior under the physiological loading. This work proposed a novel twist scaffold, and its mass transport efficiency under physiological loading conditions was evaluated by a fluid-structure interaction analysis. The results showed that compared to the non-twist scaffold, the twist scaffold could form a rotating flow under the physiological loading, which enhanced the mass transport and generated more appropriate wall shear stress (WSS) to promote bone regeneration. This highlighted the better mass transport efficiency of the twist scaffold. Therefore, getting twist may be a promising design strategy for future BTE scaffolds, and the fluid-structure interaction approach may be a more reliable method for bone regeneration studies in either in vivo or in vitro systems.

Influence of Rigid-Elastic Artery Wall of Carotid and Coronary Stenosis on Hemodynamics

Cardiovascular system abnormalities can result in serious health complications. By using the fluid-structure interaction (FSI) procedure, a comprehensive realistic approach can be employed to accurately investigate blood flow coupled with arterial wall response. The hemodynamics was investigated in both the coronary and carotid arteries based on the arterial wall response. The hemodynamics was estimated based on the numerical simulation of a comprehensive three-dimensional non-Newtonian blood flow model in elastic and rigid arteries. For stenotic right coronary artery (RCA), it was found that the maximum value of wall shear stress (WSS) for the FSI case is higher than the rigid wall. On the other hand, for the stenotic carotid artery (CA), it was found that the maximum value of WSS for the FSI case is lower than the rigid wall. Moreover, at the peak systole of the cardiac cycle (0.38 s), the maximum percentage of arterial wall deformation was found to be 1.9%. On the other hand, for the stenotic carotid artery, the maximum percentage of arterial wall deformation was found to be 0.46%. A comparison between FSI results and those obtained by rigid wall arteries is carried out. Findings indicate slight differences in results for large-diameter arteries such as the carotid artery. Accordingly, the rigid wall assumption is plausible in flow modeling for relatively large diameters such as the carotid artery. Additionally, the FSI approach is essential in flow modeling in small diameters.

Wall Shear Stress Analysis and Optimization in Tissue Engineering TPMS Scaffolds

When designing scaffolds for bone tissue engineering (BTE), the wall shear stress (WSS), due to the fluid flow inside the scaffold, is an important factor to consider as it influences the cellular process involved in new tissue formation. The present work analyzed the average WSS in Schwartz diamond (SD) and gyroid (SG) scaffolds with different surface topologies and mesh elements using computational fluid dynamics (CFD) analysis. It was found that scaffold meshes with a smooth surface topology with tetrahedral elements had WSS levels 35% higher than the equivalent scaffold with a non-smooth surface topology with hexahedral elements. The present work also investigated the possibility of implementing the optimization algorithm simulated annealing to aid in the design of BTE scaffolds with a specific average WSS, with the outputs showing that the algorithm was able to reach WSS levels in the vicinity of 5 mPa (physiological range) within the established limit of 100 iterations. This proved the efficacy of combining CFD and optimization methods in the design of BTE scaffolds.

Portable hand-held bioprinters promote in situ tissue regeneration

Three-dimensional bioprinting, as a novel technique of fabricating engineered tissues, is positively correlated with the ultimate goal of regenerative medicine, which is the restoration, reconstruction, and repair of lost and/or damaged tissue function. The progressive trend of this technology resulted in developing the portable hand-held bioprinters, which could be used quite easily by surgeons and physicians. With the advent of portable hand-held bioprinters, the obstacles and challenges of utilizing statistical bioprinters could be resolved. This review attempts to discuss the advantages and challenges of portable hand-held bioprinters via in situ tissue regeneration. All the tissues that have been investigated by this approach were reviewed, including skin, cartilage, bone, dental, and skeletal muscle regeneration, while the tissues that could be regenerated via this approach are targeted in the authors' perspective. The design and applications of hand-held bioprinters were discussed widely, and the marketed printers were introduced. It has been prospected that these facilities could ameliorate translating the regenerative medicine science from the bench to the bedside actively.

Geometry-Based Computational Fluid Dynamic Model for Predicting the Biological Behavior of Bone Tissue Engineering Scaffolds

The use of biocompatible and biodegradable porous scaffolds produced via additive manufacturing is one of the most common approaches in tissue engineering. The geometric design of tissue engineering scaffolds (e.g., pore size, pore shape, and pore distribution) has a significant impact on their biological behavior. Fluid flow dynamics are important for understanding blood flow through a porous structure, as they determine the transport of nutrients and oxygen to cells and the flushing of toxic waste. The aim of this study is to investigate the impact of the scaffold architecture, pore size and distribution on its biological performance using Computational Fluid Dynamics (CFD). Different blood flow velocities (BFV) induce wall shear stresses (WSS) on cells. WSS values above 30 mPa are detrimental to their growth. In this study, two scaffold designs were considered: rectangular scaffolds with uniform square pores (300, 350, and 450 µm), and anatomically designed circular scaffolds with a bone-like structure and pore size gradient (476-979 µm). The anatomically designed scaffolds provided the best fluid flow conditions, suggesting a 24.21% improvement in the biological performance compared to the rectangular scaffolds. The numerical observations are aligned with those of previously reported biological studies.

Revealing the compressive and flow properties of novel bone scaffold structure manufactured by selective laser sintering technique

Additive manufacturing is revolutionizing the field of medical sciences through its key application in the development of bone scaffolds. During scaffold fabrication, achieving a good level of porosity for enhanced mechanical strength is very challenging. The bone scaffolds should hold both the porosity and load withstanding capacity. In this research, a novel structure was designed with the aim of the evaluation of flexible porosity. A CAD model was generated for the novel structure using specific input parameters, whereas the porosity was controlled by varying the input parameters. Poly Amide (PA 2200) material was used for the fabrication of bone scaffolds, which is a biocompatible material. To fabricate a novel structure for bone scaffolds, a Selective Laser Sintering machine (SLS) was used. The displacement under compression loads was observed using a Universal Testing Machine (UTM). In addition to this, numerical analysis of the components was also carried out. The compressive stiffness found through the analysis enables the verification of the load withstanding capacity of the specific bone scaffold model. The experimental porosity was compared with the theoretical porosity and showed almost 29% to 30% reductions when compared to the theoretical porosity. Structural analysis was carried out using ANSYS by changing the geometry. Computational Fluid Dynamics (CFD) analysis was carried out using ANSYS FLUENT to estimate the blood pressure and Wall Shear Stress (WSS). From the CFD analysis, maximum pressure of 1.799 Pa was observed. Though the porosity was less than 50%, there was not much variation of WSS. The achievement from this study endorses the great potential of the proposed models which can successfully be adapted for the required bone implant applications.

Challenges in computational fluid dynamics applications for bone tissue engineering

Bone injuries or defects that require invasive surgical treatment are a serious clinical issue, particularly when it comes to treatment success and effectiveness. Accordingly, bone tissue engineering (BTE) has been researching the use of computational fluid dynamics (CFD) analysis tools to assist in designing optimal scaffolds that better promote bone growth and repair. This paper aims to offer a comprehensive review of recent studies that use CFD analysis in BTE. The mechanical and fluidic properties of a given scaffold are coupled to each other via the scaffold architecture, meaning an optimization of one may negatively affect the other. For example, designs that improve scaffold permeability normally result in a decreased average wall shear stress. Linked with these findings, it appears there are very few studies in this area that state a specific application for their scaffolds and those that do are focused on in vitro bioreactor environments. Finally, this review also demonstrates a scarcity of studies that combine CFD with optimization methods to improve scaffold design. This highlights an important direction of research for the development of the next generation of BTE scaffolds.

Effect of guidewire insertion in fractional flow reserve procedure for real geometry using computational fluid dynamics

BACKGROUND

Coronary artery disease is an abnormal contraction of the heart supply blood vessel. It limits the oxygenated blood flow to the heart. Thus, diagnosing its severity helps physicians to select the appropriate treatment plan. Fractional flow reserve (FFR) is the most accurate method to pinpoint the stenosis severity. However, inserting the guidewire across stenosis may cause a false overestimation of severity.

METHODS

To estimate the errors due to guidewire insertion, reconstructed three-dimensional coronary artery geometry from a patient-specific scan is used. A comprehensive three-dimensional blood flow model is developed. Blood is considered non-Newtonian and the flow is pulsatile. The model is numerically simulated using realistic boundary conditions.

RESULTS

The FFR value is calculated and compared with the actual flow ratio. Additionally, the ratio between pressure drop and distal dynamic pressure (CDP) is studied. The obtained results for each case are compared and analyzed with the case without a guidewire. It was found that placing the guidewire leads to overestimating the severity of moderate stenosis. It reduces the FFR value from 0.43 to 0.33 with a 23.26% error compared to 0.44 actual flow ratio and the CDP increases from 5.31 to 7.2 with a 35.6% error. FFR value in mild stenosis does not have a significant change due to placing the guidewire. The FFR value decreases from 0.83 to 0.82 compared to the 0.83 actual flow ratio.

CONCLUSION

Consequently, physicians should consider these errors while deciding the treatment plan.

Porous Geometry Guided Micro-mechanical Environment Within Scaffolds for Cell Mechanobiology Study in Bone Tissue Engineering

Mechanobiology research is for understanding the role of mechanics in cell physiology and pathology. It will have implications for studying bone physiology and pathology and to guide the strategy for regenerating both the structural and functional features of bone. Mechanobiological studies in vitro apply a dynamic micro-mechanical environment to cells via bioreactors. Porous scaffolds are commonly used for housing the cells in a three-dimensional (3D) culturing environment. Such scaffolds usually have different pore geometries (e.g. with different pore shapes, pore dimensions and porosities). These pore geometries can affect the internal micro-mechanical environment that the cells experience when loaded in the bioreactor. Therefore, to adjust the applied micro-mechanical environment on cells, researchers can tune either the applied load and/or the design of the scaffold pore geometries. This review will provide information on how the micro-mechanical environment (e.g. fluid-induced wall shear stress and mechanical strain) is affected by various scaffold pore geometries within different bioreactors. It shall allow researchers to estimate/quantify the micro-mechanical environment according to the already known pore geometry information, or to find a suitable pore geometry according to the desirable micro-mechanical environment to be applied. Finally, as future work, artificial intelligent - assisted techniques, which can achieve an automatic design of solid porous scaffold geometry for tuning/optimising the micro-mechanical environment are suggested.

Numerical modelling of osteocyte growth on different bone tissue scaffolds

The most common solution for the regeneration or replacement of damaged bones is the implantation of prostheses comprising ceramic or metallic materials. However, these implants are known to cause problems such as post-operative infections, collapse of the prosthesis, and lack of osseointegration. Consequently, bone tissue engineering was established because of the limitations of such implants. Osteogenic implants offer promising solutions for bone regeneration; however, three-dimensional scaffolds should be used as supportive structures. It is challenging to correctly design these structures and their compositions or properties to provide a microenvironment that promotes tissue regeneration and expedites bone formation. Computational fluid dynamics can be used to model the main phenomena that occur in bioreactors, such as cell metabolism, nutrient transport, and cell culture growth, or to model the influence of several key mechanisms related to the fluid medium, in particular, the wall shear stress. In this work, a new numerical bone cell growth model was developed, which considered the oxygen and nutrient consumption as well as the wall shear stress effect on cell proliferation. The model was implemented using 35 three-dimensional scaffolds of different porosities, and the effect of the main geometrical parameters involved in each scaffold type was analysed. The porosity plays an important role, however, a similar porosity did not guarantee similar shear stress or cell growth among the scaffolds. Randomised trabecular scaffolds, that more closely resembled trabecular bone, showed the highest cell growth values, so these are the best candidates for cell growth in a bioreactor.

Benchmarking Cellulose Nanocrystals Part II: New Industrially Produced Materials

The demand for industrially produced cellulose nanocrystals (CNCs) has been growing since 2012, when CelluForce Inc. opened its inaugural demonstration plant with a production capacity of 1 tonne per day. Currently, there are 10 industrial CNC producers worldwide, each producing a unique material. Thus, academic researchers and commercial users alike must consider the properties of all available CNCs and carefully select the material which will optimize the performance of their desired application. To support these efforts, this article presents a thorough characterization of four new industrially produced CNCs including sulfated CNCs from NORAM Engineering and Constructors Ltd. (in cooperation with InnoTech Alberta and Alberta-Pacific Forest Industries Inc.) and Melodea Ltd., as well as carboxylated CNCs from Anomera Inc. and Blue Goose Biorefineries Inc. These materials were benchmarked against typical lab-made, sulfated CNCs. While all CNCs were similar in size, shape, crystallinity, and suspension quality, the sulfated CNCs had a higher surface charge density than their carboxylated counterparts, leading to higher colloidal stability. Additionally, significant differences in the rheological profiles of aqueous CNC suspensions, as well as CNC thermal stability and self-assembly behavior, were observed. As such, this article highlights both the subtle and significant differences between five CNC types and acts as a guide for end-users looking to optimize the performance of CNC-based materials.

Numerical-experimental analysis of the permeability-porosity relationship in triply periodic minimal surfaces scaffolds

Bone Tissue Engineering has been focusing on improving the current methods for bone repair, being the use of scaffolds presented as an upgrade to traditional surgery techniques. Scaffolds are artificially porous matrices, meant to promote cell seeding and proliferation, being these properties influenced by the permeability of the structure. This work employed experimental pressure drop tests and Computational Fluid Dynamics models to assess permeability (and fluid streamlines) within different triply periodic minimal surfaces scaffold geometries (Schwarz D, Gyroid and Schwarz P). The pressure outputs from the computational analysis presented a good correlation with the experimental results, with R² equal to 0.903; they have also shown that a lower porosity may not mean a lower permeability if the geometry is altered, such as the difference between 60% porous Gyroid scaffolds (8.1*10^-9 mm²) and 70% porous Schwarz D scaffolds (7.1*10^-9 mm²). Fluid streamlines revealed how the Gyroid geometries are the most appropriate design for most bone tissue engineering applications, due to their consistent fluid permeation, followed by Schwarz D. The Schwarz P geometries have shown flat streamlines and significant variation of the permeability with the porosity (an increase of 10% in their porosity lead to an increase in the permeability from 5.1*10^-9 mm² to 11.7*10^-9 mm²), which would imply a poor environment for cell seeding and proliferation.

Analysis of the effect of guidewire position on stenosis diagnosis using computational fluid dynamics

BACKGROUND

Fractional flow reserve is an accurate method for diagnosing stenosis. The difficulty in using this procedure lies in placing the guidewire precisely at the blood vessel centerline. Owing to the long distance between the insertion point and the stenosis, a guidewire inclination can occur. Therefore, the main objective of this study is to investigate how the measured pressure in a blood vessel varies with the guidewire position.

METHODS

A three-dimensional model of blood flow is developed and numerically simulated. Two positions and two inclination angles from the blood vessel centerline and three throat diameters are investigated. The predicted results are validated using the available experimental data. The predicted results and actual measurements are observed to agree well with each other.

RESULTS

The pressure drop coefficient (CDP) increases because of guidewire insertion. When the guidewire is placed at inclined positions in moderate stenosis, the values of CDP are 66 and 68, depending on the inclination angle; the errors in CDP are 69% and 76%, respectively. At a high flow rate, the errors are reduced to 67% and 70%, respectively. The error in the CDP ranges from 42% to 61% when the guidewire is placed parallel to the centerline. For severe stenosis, the CDP is nearly the same at all positions and varies between 240 and 250; without a guidewire, the CDP is 163.

CONCLUSIONS

The findings confirmed that practitioners should be aware of the guidewire position during the operation. The displacement of the guidewire should be estimated, and the corresponding error must be considered.

Mechanical performance of highly permeable laser melted Ti6Al4V bone scaffolds

Critically engineered stiffness and strength of a scaffold are crucial for managing maladapted stress concentration and reducing stress shielding. At the same time, suitable porosity and permeability are key to facilitate biological activities associated with bone growth and nutrient delivery. A systematic balance of all these parameters are required for the development of an effective bone scaffold. Traditionally, the approach has been to study each of these parameters in isolation without considering their interdependence to achieve specific properties at a certain porosity. The purpose of this study is to undertake a holistic investigation considering the stiffness, strength, permeability, and stress concentration of six scaffold architectures featuring a 68.46-90.98% porosity. With an initial target of a tibial host segment, the permeability was characterised using Computational Fluid Dynamics (CFD) in conjunction with Darcy's law. Following this, Ashby's criterion, experimental tests, and Finite Element Method (FEM) were employed to study the mechanical behaviour and their interdependencies under uniaxial compression. The FE model was validated and further extended to study the influence of stress concentration on both the stiffness and strength of the scaffolds. The results showed that the pore shape can influence permeability, stiffness, strength, and the stress concentration factor of Ti6Al4V bone scaffolds. Furthermore, the numerical results demonstrate the effect to which structural performance of highly porous scaffolds deviate, as a result of the Selective Laser Melting (SLM) process. In addition, the study demonstrates that stiffness and strength of bone scaffold at a targeted porosity is linked to the pore shape and the associated stress concentration allowing to exploit the design freedom associated with SLM.

A multiscale computational fluid dynamics approach to simulate the micro-fluidic environment within a tissue engineering scaffold with highly irregular pore geometry

Mechanical stimulation can regulate cellular behavior, e.g., differentiation, proliferation, matrix production and mineralization. To apply fluid-induced wall shear stress (WSS) on cells, perfusion bioreactors have been commonly used in tissue engineering experiments. The WSS on cells depends on the nature of the micro-fluidic environment within scaffolds under medium perfusion. Simulating the fluidic environment within scaffolds will be important for gaining a better insight into the actual mechanical stimulation on cells in a tissue engineering experiment. However, biomaterial scaffolds used in tissue engineering experiments typically have highly irregular pore geometries. This complexity in scaffold geometry implies high computational costs for simulating the precise fluidic environment within the scaffolds. In this study, we propose a low-computational cost and feasible technique for quantifying the micro-fluidic environment within the scaffolds, which have highly irregular pore geometries. This technique is based on a multiscale computational fluid dynamics approach. It is demonstrated that this approach can capture the WSS distribution in most regions within the scaffold. Importantly, the central process unit time needed to run the model is considerably low.

Effect of scaffold architecture on cell seeding efficiency: A discrete phase model CFD analysis

Within perfusion cell culture systems, scaffold architecture is able to control important biological parameters such as permeability and fluid flow-induced shear stress. As well, one of the main factors affecting the final fate of this process as well as optimal cell differentiation and proliferation in these systems is initial adhesion of cells to scaffolds. In this study, the effect of scaffold architecture on the adhesion of the cells was computationally investigated. For this purpose, four scaffold models including double-diamond, gyroid, FR-D, and Schwarz-primitive were designed using triply periodic minimal surface (TPMS) geometry with a constant porosity of 80%. As well, the inlet velocity of zero to simulate static cell culture and three different inlet velocities for modeling the dynamic cell culture conditions were also selected. The results showed that cell culture efficiency of scaffolds could be changed up to seven times from architecture to architecture under the same conditions. The efficiency of cell culture in scaffolds with tortuous architecture was also reported higher than those with relatively straight microchannels. In terms of culture methods, unlike dynamic cell culture model in which almost a homogeneous cell distribution was observed in static cell culture simulation, more cells adhered, but they had agglomerated in the scaffold entrance regions and had failed to reach all regions. The results of this study shed more light on the selection and design of scaffold architecture for optimal cell culture in tissue engineering.