Next Article in Journal
A Time Division Multiplexing Inspired Lightweight Soft Exoskeleton for Hip and Ankle Joint Assistance
Next Article in Special Issue
Development of an Automated, Non-Enzymatic Nucleic Acid Amplification Test
Previous Article in Journal
Microfluidic Single-Cell Proteomics Assay Chip: Lung Cancer Cell Line Case Study
Previous Article in Special Issue
Toward Development of a Label-Free Detection Technique for Microfluidic Fluorometric Peptide-Based Biosensor Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Computational Simulations in Advanced Microfluidic Devices: A Review

by
Violeta Carvalho
1,2,*,
Raquel O. Rodrigues
3,
Rui A. Lima
1,4 and
Senhorinha Teixeira
2
1
MEtRICs, Campus de Azurém, University of Minho, 4800-058 Guimarães, Portugal
2
ALGORITMI, Campus de Azurém, University of Minho, 4800-058 Guimarães, Portugal
3
Center for MicroElectromechanical Systems (CMEMS-UMinho), Campus de Azurém, University of Minho, 4800-058 Guimarães, Portugal
4
CEFT, R. Dr. Roberto Frias, Faculty of Engineering of the University of Porto (FEUP), 4200-465 Porto, Portugal
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(10), 1149; https://doi.org/10.3390/mi12101149
Submission received: 9 August 2021 / Revised: 14 September 2021 / Accepted: 21 September 2021 / Published: 23 September 2021

Abstract

:
Numerical simulations have revolutionized research in several engineering areas by contributing to the understanding and improvement of several processes, being biomedical engineering one of them. Due to their potential, computational tools have gained visibility and have been increasingly used by several research groups as a supporting tool for the development of preclinical platforms as they allow studying, in a more detailed and faster way, phenomena that are difficult to study experimentally due to the complexity of biological processes present in these models—namely, heat transfer, shear stresses, diffusion processes, velocity fields, etc. There are several contributions already in the literature, and significant advances have been made in this field of research. This review provides the most recent progress in numerical studies on advanced microfluidic devices, such as organ-on-a-chip (OoC) devices, and how these studies can be helpful in enhancing our insight into the physical processes involved and in developing more effective OoC platforms. In general, it has been noticed that in some cases, the numerical studies performed have limitations that need to be improved, and in the majority of the studies, it is extremely difficult to replicate the data due to the lack of detail around the simulations carried out.

Graphical Abstract

1. Introduction

Over the years, the demand for more reliable and effective preclinical models has grown as we try to reduce and, where possible, replace the utilization of animal models for assessing drug efficacy and safety [1]. Although animal models provide a physiological structure similar to that of human cells and organs, they fail to adequately forecast the behavior of and represent human metabolism. These experiments also require laborious protocols, present ethical issues, and are expensive [2,3,4,5]. Due to the low throughput of in vivo testing, researchers have been using advanced microfluidic devices that integrate biomodels to replicate the physiology and function of cells and/or organs. A microfluidic device is the building block of the microfluidic chip family, defined as a set of technologies embedded in a chip that contains microchannels with a dimension of tens to hundreds of micrometers to process a small volume of fluids. These microdevices first emerged mainly for analytical purposes, and due to the small amount of reagents and sample requirements, they attracted increasing interest thanks to their several advantages over other standard technologies, such as the short time needed for analysis, reduction in fabrication and reagent costs, miniaturization, selectivity, sensitivity, portability, and biocompatibility. In time, the ability to integrate multiple processes in a single device led to the concept of lab-on-a-chip or micro total analysis systems (µTASs), which were first introduced by Andreas Manz in 1990 [6,7]. Their versatility enables their application in a variety of fields, from environmental monitoring and the food industry and microelectronics to biomedical or clinical applications.
Their biocompatibility and ability to precisely study, evaluate, or diagnose individual or multiple cells in a dynamic and accessible way brought forth the concept of cells-on-a-chip. In these biodevices, cells are inserted (mixed with an organic fluid) into the microfluidic chip, which allows assessing their unique function and behavior outside of the body, maintaining their integrity and biological functionality.
Lastly, and more recently, another concept has been introduced: organ-on-a-chip (OoC), also known as microphysiological systems. In these advanced microfluidic devices, 3D cell structures containing mono or multiple cell lines are embedded in an extracellular matrix to emulate the complexity and architecture of the mimicked living organ [8,9]. These devices allow more complex studies and analyses which are physiologically closer to what we can find in a human body system. In these organ-on-a-chip devices, studies can be devoted to single organs and functions, but also to the interaction and response of multiple organs, where the ultimate goal is a body-on-a-chip which is able to replicate the entire human body in a single platform [10,11]. Among the several attributes of these models, their ability to accelerate research, accurately reproduce the micro-environment of human cells, and be mass-produced in a cost-effective way should be stressed [12,13]. Nevertheless, there are also certain drawbacks associated with them, which have to be overcome over time. Frequently, polydimethylsiloxane (PDMS) is used in the fabrication of biomedical and advanced microfluidic devices owing to its outstanding properties, including optical transparency, permeability to gases, and biocompatibility [14,15,16,17,18,19,20,21,22]. However, small hydrophobic molecules, such as biomolecules, proteins, and drugs, can be adsorbed to the surface and consequently compromise the quality of the experiments’ analysis in OoC platforms [23,24,25]. Nonetheless, as the advantages of these devices exceed their limitations, these models continue to be a valuable tool in the study of new health treatments. Many researchers have developed several OoCs to simulate the function of the heart [26], lungs [27], liver [28], kidney [29], brain [30], bone [31], etc., or a combination of several organs, known as body-on-a-chip [32]. The last one allows a better understanding of not only the efficacy of a drug, but also its potential toxicity in other organs, which in turn guarantees a better validity of the results and extrapolation to reality [33]. In addition to the previous models, recently, Si et al. [34] used a human-airway-on-a-chip to investigate the antimalarial drug amodiaquine as a powerful inhibitor of infection with SARS-CoV-2. The research team modeled a human bronchial airway epithelium and pulmonary endothelium infected with pseudotyped severe acute respiratory syndrome coronavirus 2 and verified that the antimalarial drug amodiaquine has the potential to inhibit infection. Hence, the use of OoC helped in the rapid advances witnessed in the testing of new drugs to fight the new coronavirus, which has become an urgent need for the fight against the COVID-19 pandemic.
Along with the progress made in the development of new microphysiological systems, computational tools have played an important role. These allow complementing experimental studies, from an engineering point of view, allowing us to study physical phenomena, assess the viability of a device, and optimize the design at a lower cost [35,36,37,38,39]. Furthermore, their capability to determine critical parameters which are difficult to measure via experimental methods, such as pressure, velocity, shear rate, and temperature, constitute other benefits of using numerical methods. For these reasons, researchers have complemented their studies with numerical methods [40,41,42]. Although a substantial number of review papers exist in the field of microfluidic OoC-based systems, only a small number have addressed the different approaches using numerical simulations as an auxiliary tool in the development of these devices [11,35,37,43]. Hence, the present work aims to provide an overview of the most recent progress in numerical studies in advanced microfluidic devices and how researchers have applied computational tools in their investigations.

2. Applications of Numerical Simulations in Advanced Microfluidic Devices

In this section, the usefulness of conducting numerical simulations in parallel with experimental studies is presented, namely, in the study of fluid flow and mass transfer, nanoparticle development, and also in the optimization of advanced microfluidic devices, i.e., OoC and cell-on-a-chip.

2.1. Fluid Flow and Mass Transfer

A key factor that affects the functionality of OoC platforms is fluid flow, which includes the evaluation of several parameters, such as velocity fields, shear stress, and the concentration of vital parameters for cellular function, such as oxygen, glucose, and temperature.
The flow of the culture medium is of great importance to maintain an adequate amount of nutrients for the cell culture as well as the shear stresses exposed to cells, and therefore, understanding flow behavior is vital in the design of OoC microfluidic devices [44]. Kocal et al. [45] simulated flow behavior in a cancer microenvironment to understand the uniformity of fluid shear stress along with microfluidic devices. The results showed that adherent cancer cells experienced a homogenous wall shear stress of up to 25 μPa in the microchannels (Figure 1a). Another interesting study was conducted by Lo et al. [46]. They investigated the influence of shear stress and antioxidant concentrations on the production of reactive oxygen species in lung cancer cells by combining experimental and numerical simulations, and analogous results were obtained for both methods. They verified that an increase in shear stress leads to an increase in the production of reactive oxygen species (Figure 1b). Additionally, two antioxidants, α-tocopherol and ferulic acid, were tested to verify their ability to reduce reactive oxygen species. It was found that a high dose of α-tocopherol was not able to eliminate reactive oxygen species, but a lower dose could. The same outcome was not observed for ferulic acid. This antioxidant could eliminate the production of reactive oxygen species in a concentration-independent way. Kou et al. [47], on the other hand, developed a microfluidic system with the ability to create four different shear flows in only one device to study the cytosolic calcium concentration dynamics of osteoblasts. They observed that cytosolic calcium concentration in the osteoblasts increased proportionally to the intensity of shear stress from 0.03 to 0.30 Pa (Figure 1c). Komen et al. [48] developed a microfluidic device with the ability to expose cancer cells to an in-vivo-like concentration profile of a drug and measure the effectiveness on-chip. This device comprised a cell culture chamber and a drug-dosing channel separated by a transparent membrane so that the drug is not exposed to shear stresses and also to allow a label-free growth quantification. They simulated cell exposure and verified that it followed the blood concentrations measured in vivo.
The diffusion of species is also of utmost importance for cell maintenance, and porous zones are commonly modeled to study these phenomena. Wong et al. [36] created a computational model of a microfluidic device that incorporates flow-induced shear stress and monitors cell-secreted biomolecules. For that, they modeled a bilayer device in which cells are cultured on a porous membrane support. This way, researchers can adjust the device to operate with physiological shear stresses, while ensuring the rapid transport of secreted biomolecules to the biosensor (Figure 2a). In another study [49], Wong et al. used the same model but to simulate the concentration of electroactive tracers, which they validated with experiments. A similar model was proposed by Chen et al. [50], consisting of a porous membrane-separated and a microfluidic chip channel with two layers. In their study, they evaluated the dependences of flow features on the chip geometry and membrane permeability using immersed boundary methods—IBM. The results presented elucidated the dependences of flow flux, wall shear stress, and transmembrane pressure difference on the chip geometry (channel height and length) and membrane permeability, which were consistent with the experimental results. For a membrane with lower permeability, a lengthier channel was required to achieve a permeability-independent stage. Moreover, it was also observed that wall shear stress was dependent on channel length and membrane permeability. Mosavati et al. [51] developed a numerical model of a three-dimensional placenta-on-a-chip model with a porous medium to model the membrane separating the two channels (Figure 2b). The researchers aimed to analyze the unsteady flow into microchannels and glucose concentration profiles in different locations of the microfluidic chip. After validation of the model, Mosavati et al. studied the effects of flow rate and membrane porosity on glucose diffusion through the placental barrier. Similarly, in a recent study [4], Banaeiyan et al. explored the diffusion of glucose, since this is the main element present in the culture medium, through diffusion channels. By contrast, Hu and Li [52] numerically investigated oxygen concentration in the interior of a three-dimensional microchannel containing tumor spheroids and tested two different velocities, 50 µm/s and 100 µm/s. They observed that with lower velocity, oxygen concentration downstream of the spheroid surface is not enough for spheroid growth. By raising the perfusion velocity, oxygen concentration increases and may consequently enhance the growth of the spheroids. Another model comprising a permeable layer was presented by Zahorodny-Burke and Elias [53]. Briefly, they developed a rectangular channel with a polymer layer permeable to oxygen linked to a glass substrate seeded with one layer of oxygen-consuming cells. In contrast, Bhise and coworkers [54] did not simulate a porous membrane where oxygen passes through it. Instead, they modeled hydrogels as a porous medium with a homogenous volumetric oxygen consumption rate related to the total number of encapsulated cells. Other researchers have studied not only the diffusion of species but also the convection of carbon dioxide—CO2 [55]. A computational model was developed to investigate CO2 transport in a cylindrical microfluidic device. The level of CO2 at the base of the water chamber for different CO2 supplying concentrations was investigated, and it was found that by increasing the level of CO2 in the feeding gas to above 10%, CO2 concentration at the bottom of the chamber reaches the desired level of 5%.
The study of the effects of radiation treatment on microcirculation was recently addressed by Cicchetti et al. [56]. The authors developed a numerical model to understand the effects of radiation therapy on normal tissues surrounding a tumor. They modeled the vasodilation and variation of membrane permeability as inputs with consequent fluid extravasation and the deposition of fibrin or the formation of edema. In a second simulation, the researchers modeled the impact of variations in vessel wall elasticity on the distribution of blood flow velocity.
Another important parameter in these types of studies is temperature because it affects not only the environment but also the behavior and function of cells [57]. Hence, its monitoring is of great importance for cell culture. Peng et al. [58] conducted a multiphysics simulation (fluid flow and heat transfer) in a microfluidic device to obtain an accurate and rapid control of on-chip local temperature control. According to the results, the device was able to operate with temperatures ranging between 2 °C and 37 °C using cooling and heating components (Figure 3a). A similar study was conducted by Das and coworkers [59], in which a microfluidic device to study the viability and activity of the cells under a temperature gradient was developed. Through numerical simulations, the temperature gradient that can be reached in the proposed chip was assessed, and it was found that the proposed chip is capable of creating temperature conditions that realistically mimic physiological/biological conditions (Figure 3b).
Although single-phase simulations are mostly used, the use of multiphase models is a valuable option to improve the precision of simulations and obtain more detailed information. To study the behavior of circulating tumor cells in a 3D bioprinted vasculature chip, Hynes et al. [60] performed numerical simulations using both the lattice Boltzmann method (LBM) and finite element method (Figure 4a). The particles were treated as either rigid or deformable spheres, and the results matched closely to the experimental findings, proving that the numerical model used was appropriate. By contrast, Ye et al. [61] used a smoothed dissipative particle dynamics (SDPD) method to model the fluid flow and an IBM for the fluid–cell interaction to study the RBC mechanics in a microfluidic chip (Figure 4b). Likewise, Zhang et al. [62], in order to simulate mechanical interactions between flow and particles (cells) for a cells-on-a-chip device, used a two-way Euler/Lagrange multiphase model. Another example of the use of multiphase models in this type of study was presented by Tanaka et al. [63]. The authors developed a hybrid pump driven by cardiomyocytes that self-organize into bridges between elastic micro-pier structures (Figure 4c). In addition, they performed numerical simulations to understand the flow pattern and distribution of flow velocity, and the results were in good agreement with the experimental results. They used the volume of fluid (VOF) model to better understand and verify the flow generation, flow pattern, and velocity distribution of the ultra-small fluid oscillation unit.
Another key parameter that affects the function of a device is the material. As previously mentioned, PDMS is preferred due to its suitable properties and low price [5,23,64]. Moreover, this material can be obtained with different air permeability, which is interesting for improving the diffusion of species in advanced microfluidic devices. Lamberti et al. [65] demonstrated that by increasing the PDMS mixing ratio, air permeability can be improved up to 300% when compared to the standard 10:1 ratio. However, this material also has certain limitations that encourage researchers to seek alternative materials, one of them being the possibility of absorption of small hydrophobic species [23,64]. Two alternative materials are poly(-methylmethacrylate) (PMMA) [66,67,68,69] and cyclic olefin copolymer (COC) [70]. Zahorodny-Burke et al. [53] studied the nature of oxygen transport within microfluidic cell culture devices considering PDMS, COC, and PMMA through a two-dimensional convection–diffusion mass transfer model. To this end, flow rate, diffusive layer thickness, and material selection were investigated in order to determine their relation to ensure adequate supplies of oxygen for cell growth. However, special attention must be paid when using COC or PMMA. Due to the low oxygen diffusion of these materials, flow rate must be optimized to ensure a sufficient supply of oxygen to the cells by applying high flow rates, at an order of 10−2 m/s (Figure 5). By contrast, Ochs and coworkers [71] investigated not only PDMS and COC, but also a different type of material, polymethylpentene (PMP). They experimentally and numerically studied the influence of PMP on oxygen concentration in cell culture and verified that this is a possible alternative for PDMS. PMP was shown to be able to deliver sufficient oxygen for various types of cells, namely, hepatocytes with a high oxygen consumption rate.

2.2. Nanoparticle Simulations

Nanoparticles (NPs) have been attracting interest for their potential use as drug delivery systems, but the optimization of their physical properties, such as their size and shape, is a point of interest to accomplish appropriate cell responses [72,73,74]. Due to the multidisciplinarity of computational software, several researchers have taken advantage of its capabilities to broaden the understanding of the transport of NPs and drugs in OoC [75,76,77,78,79]. Soheili et al. [75] performed a numerical and experimental study on the microfluidic preparation and production of chitosan NPs. To develop NPs with the desired size, they implemented a model that relates NPs’ diameter to the flow rate, and sodium tripolyphosphate (TPP) concentration as a critical factor that affects cell flow behavior. Further, the researchers found the optimal mixing efficiency of tripolyphosphate conditions by associating TPP concentration with the length of the channels and flow rate. Kwak and coworkers [79] developed a tumor-microenvironment-on-chip to recapitulate the key features of complex transport of drugs and NPs within a tumor microenvironment. By combining computational simulations and experimental tests, they obtained more detailed information about the dynamic transport behavior of NPs and concluded that NPs should be conceived considering their interactions with the tumor microenvironment.
Other authors have been motivated by the “catch and release” of biomolecules observed in physiological processes and used numerical simulations to design synthetic systems that can catch, transport, and release targeted species inside the surrounding solution, which can be useful for effective separation processes within microfluidic devices [77]. Their setup allows the determination of the conditions where the oscillating fins of the device can selectively bind and ‘‘catch’’ target NPs in the upper channel, and then release them in the lower channel. Their findings have provided insights into fabrication devices to detect and separate target molecules from complex fluids. This topic has been addressed by various researchers [80,81,82,83]. Maleki and coworkers [78] developed a hybrid drug carrier of doxorubicin (DOX) due to the high consumption rate of riboflavin in cancerous cells. In parallel, molecular simulations to study the effects of the microfluidic environment and shift metal dichalcogenide nanolayers on the stability, size, and self-assembly interaction energies of the nanocarriers were conducted (Figure 6a). The results showed that among the studied configurations, the most suitable ones for the adsorption of DOX molecules are MoSe2, MoSSe, MoO3, and MoS2, respectively. Furthermore, the appearance of the new coronavirus has triggered several scientific investigations, and the study conducted by Arefi and colleagues [76] was one of them. They studied the effect of airborne NPs on human health through both numerical simulations and experimental work. Transport and adsorption in a microfluidic lung-on-a-chip device were investigated, and two different multiphase models were used. For the fine NPs, particle transport was estimated using an Eulerian advection-diffusion approach, while for coarse NPs, the Lagrangian particle tracking approach was used (Figure 6b). One of their findings showed that during physical exercise, particle deposition increases, and thus, outdoor exercise on days with poor air quality should be reduced.

3. Numerical Optimization

The use of numerical simulations to optimize devices is one of the main steps that must be carried out during the design process. This helps researchers to improve and accelerate the development of more efficient and representative models for a given problem by subjecting the chip to design changes and improvements. Experimental conditions can also be optimized in accordance with numerical simulations. For instance, Karakas and coworkers [84] developed a microfluidic chip for screening individual cancer cells and performed numerical simulations in order to optimize the exposure time for autophagy, a cellular mechanism where proteins are consumed and recycled to give another source of energy to cells (Figure 7a). This numerical study allowed investigators to determine the minimum exposure time required to ensure the viability of the experiments, saving time and resources during experimental tests. On the other hand, Zhang et al. [81] studied the effect of varying the inlet flow rate of blood and buffer on the separation performance of circulating tumor cells in a microfluidic chip and optimized the operating conditions to increase separation efficiency. Another way to improve a device’s accuracy was presented by Jun-Shan et al. [85], who developed a microfluidic chip with micropillar arrays for 3D cell culture, and using numerical simulations, the space between micropillars was optimized, which allowed the nutrients in the medium to quickly diffuse into the chamber, while cell metabolites could diffuse out of the chamber in a timely manner (Figure 7b). A similar study with micropillars was performed by Chen and coworkers [86]. However, in this case, the authors studied pillars with circular, elliptical, and square cross-sections assembled in both aligned and staggered patterns. Through numerical simulations, the researchers found that in the latter case, fluid flows through both the center of the array and around the pillars. These observations led the research team to address this approach since in staggered arrangements, fluid is better distributed throughout the device. Zhang et al. [62] used a two-way Euler/Lagrange multiphase model to simulate mechanical interactions between flow and particles (cells) for cells-on-a-chip devices. The authors developed three different designs and studied the effect of using different cell densities, inlet flow velocities, and inlet cell numbers. The results showed that regions of lower strain rates had high cell concentrations, and when lower inlet velocities (10 and 20 µm/s) were used, high cell concentrations were observed. However, they observed that cells were not able to reach the outlet, while with a velocity of 40 µm/s, some cells were.
Optimization methods are also of great importance for discovering the relationship between several variables. For instance, Huang and Nguyen [87] simulated a cell-stretching device and optimized important dimensional parameters. They performed a parametrical optimization to determine the relationship between the lateral displacement of the membrane and wall height, wall thickness, membrane thickness, and vacuum strength in order to maximize the output strain of the stretched membrane.
In Table 1, a summary of the numerical considerations used in the previously reviewed papers is presented, including the software used, the flow behavior, flow type, fluid rheology, mesh, flow analysis, and also the validation.
A brief analysis of the results reveals several blank spaces which represent the absence of information regarding the assumptions made when performing the numerical works. Additionally, the preferred software for these types of studies was found to be COMSOL Multiphysics® software, followed by Ansys software®, although researchers have begun to develop their own simulation software to overcome some of the limitations of the existing ones [56,60]. Moreover, a fundamental step in numerical simulations is mesh development. Without guaranteeing that the mesh is appropriate, one cannot obtain truthful and accurate results. Hence, both mesh quality and independence have to be ensured, and this is a step that is commonly missing in the studies presented in this review article. Some authors conducted mesh independence studies, but the majority only indicate the number of elements or do not present any information. One way to carry out such studies with more accuracy is through the calculation of the grid convergence index, as represented in the study conducted by Soheili et al. [75], increasing, in this way, the reliability of the presented results. Moreover, it can be seen that the validation of these types of studies is usually disregarded, but this is an extremely important step that should be addressed whenever possible.

4. Final Remarks and Future Directions

Advanced microfluidic devices have revolutionized drug development processes and in vitro tests by allowing to save time and money but also, more importantly, to present a preclinical platform that can accurately mimic human tissues/organs and thus diminish the use of animals in clinical testing. These devices have become research accelerators, therefore rapidly achieving the status of a promising technology, which has now been proven to be useful for advances in emerging diseases, such as research related to the new coronavirus, Sars-Cov-2, and cancer treatment. Nevertheless, although this technology has been shown to be extremely valuable, it is still far from being able to precisely mimic a proper human organ, which has hindered the passage of these devices to the market. Therefore, research and development of more complex and reliable advanced microfluidic devices will continue to progress for several years, and numerical simulations are expected to boost this process by not only allowing mass transport insights, but also reducing the costs and saving time in research. Hence, computational methods are envisioned to become a fundamental tool for the advance of research in advanced microfluidic platforms.
The present review intended to encourage current and future researchers to include this tool in their studies, and it has shown the importance of using numerical simulations in parallel with experimental works to improve the performance and quality of advanced microfluidic devices. Nevertheless, it was found that most numerical studies did not include details about the simulations used, which makes it unfeasible for other researchers to recreate such simulations and understand what methods were actually used to obtain certain outputs.
In addition to the advantages of using numerical tools, there are also certain challenges and limitations associated with numerical modeling platforms. As biomedical devices become more complicated, mathematical modeling platforms need further improvement and validation to become more precise, and this is a lengthy and difficult process. Moreover, in these biological studies, there are a lot of physical and chemical processes occurring both at a micro- and macro-scale, and it becomes difficult to assess and represent these in one numerical model, making simulations extremely complex and time-consuming. However, given the advancements in this type of technology, the field is expected to grow by leaps and bounds, and more biological features will be included in the existing software. For instance, molecular dynamics simulations may be one good alternative to perform more complex simulations.

Author Contributions

Conceptualization, V.C., R.O.R., R.A.L., and S.T.; methodology, V.C.; validation, R.O.R., R.A.L., and S.T.; formal analysis, R.O.R., R.A.L., and S.T.; writing—original draft preparation, V.C.; writing—review and editing, V.C., R.O.R., R.A.L., and S.T.; supervision, R.O.R., R.A.L., and S.T.; funding acquisition, R.O.R., R.A.L., and S.T. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful for the funding of FCT through projects NORTE-01-0145-FEDER-029394, NORTE-01-0145-FEDER-030171 funded by COMPETE2020, NORTE2020, PORTUGAL2020, and FEDER. This work was also supported by Fundação para a Ciência e a Tecnologia (FCT) under the strategic grants UIDB/04077/2020, UIDB/00319/2020, UIDB/04436/2020, and UIDB/00532/2020.

Acknowledgments

Violeta Carvalho acknowledges the PhD scholarship UI/BD/151028/2021 attributed by FCT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Filipovic, N.; Nikolic, M.; Sustersic, T. Simulation of organ-on-a-chip systems. In Biomaterials for Organ and Tissue Regeneration; Vrana, N.E., Knopf-Marques, H., Barthes, J., Eds.; Woodhead Publishing Series in Biomaterial; Woodhead Publishing: Sawston, UK, 2020. [Google Scholar]
  2. Soker, S.; Carolina, N. A reductionist metastasis-on-a-chip platform for in vitro tumor progression modeling and drug screening. Biotechnol. Bioeng. 2016, 113, 2020–2032. [Google Scholar] [CrossRef]
  3. Metiner, P.S.; Gulce-Iz, S.; Biray-Avci, C. Bioengineering-inspired three-dimensional culture systems: Organoids to create tumor microenvironment. Gene 2018, 686, 203–212. [Google Scholar] [CrossRef]
  4. Banaeiyan, A.A.; Theobald, J.; Paukstyte, J.; Wölfl, S.; Adiels, C.B.; Goksör, M. Design and fabrication of a scalable liver-lobule-on-a-chip microphysiological platform. Biofabrication 2017, 9, 015014. [Google Scholar] [CrossRef] [PubMed]
  5. Carvalho, V.; Maia, I.; Souza, A.; Ribeiro, J.; Costa, P.; Puga, H.; Teixeira, S.; Lima, R.A. In vitro biomodels in stenotic arteries to perform blood analogues flow visualizations and measurements: A review. Open Biomed. Eng. J. 2020, 14, 87–102. [Google Scholar] [CrossRef]
  6. Rodrigues, R.O.; Lima, R.; Gomes, H.T.; Silva, A.M.T. Polymer microfluidic devices: An overview of fabrication methods. U. Porto. J. Eng. 2017, 1, 67–79. [Google Scholar] [CrossRef] [Green Version]
  7. Manz, A.; Graber, N.; Widmer, H. Miniaturized total chemical analysis systems: A novel concept for chemical sensing. Sens. Actuators B Chem. 1990, 1, 244–248. [Google Scholar] [CrossRef]
  8. Huh, D.; Kim, H.J.; Fraser, J.P.; E Shea, D.; Khan, M.; Bahinski, A.; Hamilton, G.A.; Ingber, D.E. Microfabrication of human organs-on-chips. Nat. Protoc. 2013, 8, 2135–2157. [Google Scholar] [CrossRef]
  9. Fabre, K.M.; Livingston, C.; Tagle, D.A. Organs-on-chips (microphysiological systems): Tools to expedite efficacy and toxicity testing in human tissue. Exp. Biol. Med. 2014, 239, 1073–1077. [Google Scholar] [CrossRef]
  10. Wu, Q.; Liu, J.; Wang, X.; Feng, L.; Wu, J.; Zhu, X.; Wen, W.; Gong, X. Organ-on-a-chip: Recent breakthroughs and future prospects. Biomed. Eng. Online 2020, 19, 9. [Google Scholar] [CrossRef] [Green Version]
  11. Carvalho, V.; Gonçalves, I.; Lage, T.; Rodrigues, R.; Minas, G.; Teixeira, S.; Moita, A.; Hori, T.; Kaji, H.; Lima, R. 3D printing techniques and their applications to organ-on-a-chip platforms: A systematic review. Sensors 2021, 21, 3304. [Google Scholar] [CrossRef]
  12. Huh, D.; Hamilton, G.A.; Ingber, D.E. From 3D cell culture to organs-on-chips. Trends Cell Biol. 2011, 21, 745–754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Polini, A.; Prodanov, L.; Bhise, N.S.; Manoharan, V.; Dokmeci, M.R.; Khademhosseini, A. Organs-on-a-chip: A new tool for drug discovery. Expert Opin. Drug Discov. 2014, 9, 335–352. [Google Scholar] [CrossRef] [PubMed]
  14. Pinho, D.; Carvalho, V.; Gonçalves, I.M.; Teixeira, S.; Lima, R. Visualization and measurements of blood cells flowing in microfluidic systems and blood rheology: A personalized medicine perspective. J. Pers. Med. 2020, 10, 249. [Google Scholar] [CrossRef] [PubMed]
  15. Faustino, V.; Rodrigues, R.O.; Pinho, D.; Costa, E.; Santos-Silva, A.; Miranda, V.; Amaral, J.S.; Lima, R. A microfluidic deformability assessment of pathological red blood cells flowing in a hyperbolic converging microchannel. Micromachines 2019, 10, 645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Catarino, S.O.; Rodrigues, R.O.; Pinho, D.; Miranda, J.M.; Minas, G.; Lima, R. Blood cells separation and sorting techniques of passive microfluidic devices: From fabrication to applications. Micromachines 2019, 10, 593. [Google Scholar] [CrossRef] [Green Version]
  17. Quirós-Solano, W.F.; Gaio, N.; Stassen, O.M.J.A.; Arik, Y.B.; Silvestri, C.; Van Engeland, N.C.A.; van der Meer, A.; Passier, R.; Sahlgren, C.; Bouten, C.; et al. Microfabricated tuneable and transferable porous PDMS membranes for organs-on-chips. Sci. Rep. 2018, 8, 13524. [Google Scholar] [CrossRef]
  18. Boas, L.V.; Faustino, V.; Lima, R.; Miranda, J.M.; Minas, G.; Fernandes, C.S.V.; Catarino, S.O. Assessment of the deformability and velocity of healthy and artificially impaired red blood cells in narrow polydimethylsiloxane (PDMS) microchannels. Micromachines 2018, 9, 384. [Google Scholar] [CrossRef] [Green Version]
  19. Lima, R.; Wada, S.; Tanaka, S.; Takeda, M.; Ishikawa, T.; Tsubota, K.-I.; Imai, Y.; Yamaguchi, T. In vitro blood flow in a rectangular PDMS microchannel: Experimental observations using a confocal micro-PIV system. Biomed. Microdevices 2008, 10, 153–167. [Google Scholar] [CrossRef] [Green Version]
  20. Pinho, D.; Muñoz-Sánchez, B.N.; Anes, C.F.; Vega, E.J.; Lima, R. Flexible PDMS microparticles to mimic RBCs in blood particulate analogue fluids. Mech. Res. Commun. 2019, 100, 103399. [Google Scholar] [CrossRef] [Green Version]
  21. Mata, A.; Fleischman, A.J.; Roy, S. Characterization of polydimethylsiloxane (PDMS) properties for biomedical micro/nanosystems. Biomed. Microdevices 2005, 7, 281–293. [Google Scholar] [CrossRef]
  22. Leclerc, E.; Sakai, Y.; Fujii, T. Cell culture in 3-dimensional microfluidic structure of PDMS (polydimethylsiloxane). Biomed. Microdevices 2003, 5, 109–114. [Google Scholar] [CrossRef]
  23. Halldorsson, S.; Lucumi, E.; Gomez-Sjoberg, R.; Fleming, R.M. Advantages and challenges of microfluidic cell culture in polydimethylsiloxane devices. Biosens. Bioelectron. 2015, 63, 218–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wang, Y.; Chen, S.; Sun, H.; Li, W.; Hu, C.; Ren, K. Recent progresses in microfabricating perfluorinated polymers (Teflons) and the associated new applications in microfluidics. Microphysiol. Syst. 2018, 2, 6. [Google Scholar] [CrossRef]
  25. Sosa-Hernández, J.E.; Villalba-Rodríguez, A.M.; Romero-Castillo, K.D.; Aguilar-Aguila-Isaías, M.A.; García-Reyes, I.E.; Hernández-Antonio, A.; Ahmed, I.; Sharma, A.; Parra-Saldívar, R.; Iqbal, H.M.N. Organs-on-a-chip module: A review from the development and applications perspective. Micromachines 2018, 9, 536. [Google Scholar] [CrossRef] [Green Version]
  26. Sakamiya, M.; Fang, Y.; Mo, X.; Shen, J.; Zhang, T. A heart-on-a-chip platform for online monitoring of contractile behavior via digital image processing and piezoelectric sensing technique. Med. Eng. Phys. 2019, 75, 36–44. [Google Scholar] [CrossRef]
  27. Zamprogno, P.; Wüthrich, S.; Achenbach, S.; Thoma, G.; Stucki, J.D.; Hobi, N.; Schneider-Daum, N.; Lehr, C.-M.; Huwer, H.; Geiser, T.; et al. Second-generation lung-on-a-chip with an array of stretchable alveoli made with a biological membrane. Commun. Biol. 2021, 4, 168. [Google Scholar] [CrossRef]
  28. Deng, J.; Wei, W.; Chen, Z.; Lin, B.; Zhao, W.; Luo, Y.; Zhang, X. Engineered liver-on-a-chip platform to mimic liver functions and its biomedical applications: A review. Micromachines 2019, 10, 676. [Google Scholar] [CrossRef] [Green Version]
  29. Benam, K.H.; Dauth, S.; Hassell, B.; Herland, A.; Jain, A.; Jang, K.-J.; Karalis, K.; Kim, H.J.; MacQueen, L.; Mahmoodian, R.; et al. Engineered in vitro disease models. Annu. Rev. Pathol. Mech. Dis. 2015, 10, 195–262. [Google Scholar] [CrossRef] [Green Version]
  30. Ndyabawe, K.; Cipriano, M.; Zhao, W.; Haidekker, M.; Yao, K.; Mao, L.; Kisaalita, W.S. Brain-on-a-chip device for modeling multiregional networks. ACS Biomater. Sci. Eng. 2021, 7, 350–359. [Google Scholar] [CrossRef]
  31. Hao, S.; Ha, L.; Cheng, G.; Wan, Y.; Xia, Y.; Sosnoski, D.M.; Mastro, A.M.; Zheng, S.-Y. A spontaneous 3D bone-on-a-chip for bone metastasis study of breast cancer cells. Small 2018, 14, e1702787. [Google Scholar] [CrossRef]
  32. Sung, J.H.; Wang, Y.I.; Sriram, N.N.; Jackson, M.; Long, C.; Hickman, J.J.; Shuler, M.L. Recent advances in body-on-a-chip systems. Anal. Chem. 2019, 91, 330–351. [Google Scholar] [CrossRef] [PubMed]
  33. Shuler, M.L. Organ-, body- and disease-on-a-chip systems. Lab Chip 2017, 17, 2345–2346. [Google Scholar] [CrossRef] [PubMed]
  34. Si, L.; Bai, H.; Rodas, M.; Cao, W.; Oh, C.Y.; Jiang, A.; Moller, R.; Hoagland, D.; Oishi, K.; Horiuchi, S.; et al. A human-airway-on-a-chip for the rapid identification of candidate antiviral therapeutics and prophylactics. Nat. Biomed. Eng. 2021, 5, 815–829. [Google Scholar] [CrossRef]
  35. Sung, J.H.; Wang, Y.; Shuler, M.L. Strategies for using mathematical modeling approaches to design and interpret multi-organ microphysiological systems (MPS). APL Bioeng. 2019, 3, 021501. [Google Scholar] [CrossRef]
  36. Wong, J.F.; Young, E.W.K.; Simmons, C.A. Computational analysis of integrated biosensing and shear flow in a microfluidic vascular model. AIP Adv. 2017, 7, 115116. [Google Scholar] [CrossRef] [Green Version]
  37. Sheidaei, Z.; Akbarzadeh, P.; Kashaninejad, N. Advances in numerical approaches for microfluidic cell analysis platforms. J. Sci. Adv. Mater. Devices 2020, 5, 295–307. [Google Scholar] [CrossRef]
  38. Carvalho, V.; Pinho, D.; Lima, R.A.; Teixeira, J.C.; Teixeira, S. Blood flow modeling in coronary arteries: A review. Fluids 2021, 6, 53. [Google Scholar] [CrossRef]
  39. Lopes, D.; Agujetas, R.; Puga, H.; Teixeira, J.; Lima, R.; Alejo, J.; Ferrera, C. Analysis of finite element and finite volume methods for fluid-structure interaction simulation of blood flow in a real stenosed artery. Int. J. Mech. Sci. 2021, 207, 106650. [Google Scholar] [CrossRef]
  40. Carvalho, V.; Rodrigues, N.; Ribeiro, R.; Costa, P.F.; Teixeira, J.C.F.; Lima, R.A.; Teixeira, S.F.C.F. Hemodynamic study in 3D printed stenotic coronary artery models: Experimental validation and transient simulation. Comput. Methods Biomech. Biomed. Eng. 2020, 24, 623–636. [Google Scholar] [CrossRef]
  41. Carvalho, V.; Rodrigues, N.; Ribeiro, R.; Costa, P.F.; Lima, R.A.; Teixeira, S.F. 3D printed biomodels for flow visualization in stenotic vessels: An experimental and numerical study. Micromachines 2020, 11, 549. [Google Scholar] [CrossRef]
  42. Doutel, E.; Viriato, N.; Carneiro, J.; Campos, J.B.; Miranda, J.M. Geometrical effects in the hemodynamics of stenotic and non-stenotic left coronary arteries—Numerical and in vitro approaches. Int. J. Numer. Methods Biomed. Eng. 2019, 35, e3207. [Google Scholar] [CrossRef]
  43. Rodrigues, R.O.; Sousa, P.C.; Gaspar, J.; Bañobre-López, M.; Lima, R.; Minas, G. Organ-on-a-chip: A preclinical microfluidic platform for the progress of nanomedicine. Small 2020, 16, e2003517. [Google Scholar] [CrossRef]
  44. Nguyen, N.T.; Wereley, S.T.; Shaegh, S.A.M. Fundamentals and Applications of Microfluidics, 3rd ed.; Artech House: Norwood, MA, USA, 2019; Volume 110. [Google Scholar]
  45. Kocal, G.C.; Guven, S.; Foygel, K.; Goldman, A.; Chen, P.; Sengupta, S.; Paulmurugan, R.; Baskin, Y.; Demirci, U. Dynamic microenvironment induces phenotypic plasticity of esophageal cancer cells under flow. Sci. Rep. 2016, 6, 38221. [Google Scholar] [CrossRef] [Green Version]
  46. Lo, K.-Y.; Zhu, Y.; Tsai, H.-F.; Sun, Y.-S. Effects of shear stresses and antioxidant concentrations on the production of reactive oxygen species in lung cancer cells. Biomicrofluidics 2013, 7, 64108. [Google Scholar] [CrossRef] [Green Version]
  47. Kou, S.; Pan, L.; van Noort, D.; Meng, G.; Wu, X.; Sun, H.; Xu, J.; Lee, I. A multishear microfluidic device for quantitative analysis of calcium dynamics in osteoblasts. Biochem. Biophys. Res. Commun. 2011, 408, 350–355. [Google Scholar] [CrossRef]
  48. Komen, J.; Westerbeek, E.Y.; Kolkman, R.W.; Roesthuis, J.; Lievens, C.; Berg, A.V.D.; Van Der Meer, A.D. Controlled pharmacokinetic anti-cancer drug concentration profiles lead to growth inhibition of colorectal cancer cells in a microfluidic device. Lab Chip 2020, 20, 3167–3178. [Google Scholar] [CrossRef]
  49. Wong, J.F.; Simmons, C.A. Microfluidic assay for the on-chip electrochemical measurement of cell monolayer permeability. Lab Chip 2019, 19, 1060–1070. [Google Scholar] [CrossRef] [PubMed]
  50. Chen, S.; Xue, J.; Hu, J.; Ding, Q.; Zhou, L.; Feng, S.; Cui, Y.; Lü, S.; Long, M. Flow field analyses of a porous membrane-separated, double-layered microfluidic chip for cell co-culture. Acta Mech. Sin. 2020, 36, 754–767. [Google Scholar] [CrossRef]
  51. Mosavati, B.; Oleinikov, A.V.; Du, E. Development of an organ-on-a-chip-device for study of placental pathologies. Int. J. Mol. Sci. 2020, 21, 8755. [Google Scholar] [CrossRef]
  52. Hu, G.; Li, D. Three-dimensional modeling of transport of nutrients for multicellular tumor spheroid culture in a microchannel. Biomed. Microdevices 2007, 9, 315–323. [Google Scholar] [CrossRef] [PubMed]
  53. Zahorodny-Burke, M.; Nearingburg, B.; Elias, A. Finite element analysis of oxygen transport in microfluidic cell culture devices with varying channel architectures, perfusion rates, and materials. Chem. Eng. Sci. 2011, 66, 6244–6253. [Google Scholar] [CrossRef]
  54. Bhise, N.S.; Manoharan, V.; Massa, S.; Tamayol, A.; Ghaderi, M.; Miscuglio, M.; Lang, Q.; Zhang, Y.S.; Shin, S.R.; Calzone, G.; et al. A liver-on-a-chip platform with bioprinted hepatic spheroids. Biofabrication 2016, 8, 014101. [Google Scholar] [CrossRef] [PubMed]
  55. Mäki, A.-J.; Peltokangas, M.; Kreutzer, J.; Auvinen, S.; Kallio, P. Modeling carbon dioxide transport in PDMS-based microfluidic cell culture devices. Chem. Eng. Sci. 2015, 137, 515–524. [Google Scholar] [CrossRef]
  56. Cicchetti, A.; Laurino, F.; Possenti, L.; Rancati, T.; Zunino, P. In silico model of the early effects of radiation therapy on the microcirculation and the surrounding tissues. Phys. Medica 2020, 73, 125–134. [Google Scholar] [CrossRef]
  57. Kim, M.; Kim, G.; Kim, D.; Yoo, J.; Kim, D.-K.; Kim, H. Numerical study on effective conditions for the induction of apoptotic temperatures for various tumor aspect ratios using a single continuous-wave laser in photothermal therapy using gold nanorods. Cancers 2019, 11, 764. [Google Scholar] [CrossRef] [Green Version]
  58. Peng, J.; Fang, C.; Ren, S.; Pan, J.; Jia, Y.; Shu, Z.; Gao, D. Development of a microfluidic device with precise on-chip temperature control by integrated cooling and heating components for single cell-based analysis. Int. J. Heat Mass Transf. 2019, 130, 660–667. [Google Scholar] [CrossRef]
  59. Das, S.K.; Chung, S.; Zervantonakis, I.; Atnafu, J.; Kamm, R.D. A microfluidic platform for studying the effects of small temperature gradients in an incubator environment. Biomicrofluidics 2008, 2, 034106. [Google Scholar] [CrossRef] [Green Version]
  60. Hynes, W.F.; Pepona, M.; Robertson, C.; Alvarado, J.; Dubbin, K.; Triplett, M.; Adorno, J.J.; Randles, A.; Moya, M.L. Examining metastatic behavior within 3D bioprinted vasculature for the validation of a 3D computational flow model. Sci. Adv. 2020, 6, eabb3308. [Google Scholar] [CrossRef]
  61. Ye, T.; Shi, H.; Phan-Thien, N.; Lim, C.T.; Li, Y. Numerical design of a microfluidic chip for probing mechanical properties of cells. J. Biomech. 2019, 84, 103–112. [Google Scholar] [CrossRef]
  62. Zhang, M.; Zheng, A.; Zheng, Z.C.; Wang, M.Z. Multiphase flow experiment and simulation for cells-on-a-chip devices. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2019, 233, 432–443. [Google Scholar] [CrossRef]
  63. Tanaka, N.; Yamashita, T.; Yalikun, Y.; Amaya, S.; Sato, A.; Vogel, V.; Tanaka, Y. An ultra-small fluid oscillation unit for pumping driven by self-organized three-dimensional bridging of pulsatile cardiomyocytes on elastic micro-piers. Sens. Actuators B Chem. 2019, 293, 256–264. [Google Scholar] [CrossRef]
  64. Subramaniam, A.; Sethuraman, S. Biomedical applications of nondegradable polymers. In Natural and Synthetic Biomedical Polymers; Elsevier Science: Amsterdam, The Netherlands, 2014. [Google Scholar]
  65. Lamberti, A.; Marasso, S.L.; Cocuzza, M. PDMS membranes with tunable gas permeability for microfluidic applications. RSC Adv. 2014, 4, 61415–61419. [Google Scholar] [CrossRef]
  66. Li, J.M.; Liu, C.; Dai, X.D.; Chen, H.H.; Liang, Y.; Sun, H.L.; Tian, H.; Ding, X.P. PMMA microfluidic devices with three-dimensional features for blood cell filtration. J. Micromechanics Microengineering 2008, 18, 95021. [Google Scholar] [CrossRef]
  67. Faustino, V.; Catarino, S.O.; Lima, R.; Minas, G. Biomedical microfluidic devices by using low-cost fabrication techniques: A review. J. Biomech. 2016, 49, 2280–2292. [Google Scholar] [CrossRef] [Green Version]
  68. Singhal, J.; Pinho, D.; Lopes, R.; Sousa, P.C.; Garcia, V.; Schutte, H.; Lima, R.; Gassmann, S. Blood flow visualization and measurements in microfluidic devices fabricated by a micromilling technique. Micro Nanosyst. 2015, 7, 148–153. [Google Scholar] [CrossRef] [Green Version]
  69. Lopes, R.; Rodrigues, R.O.; Pinho, D.; Garcia, V.; Schütte, H.; Lima, R.; Gassmann, S. Low cost microfluidic device for partial cell separation: Micromilling approach. In Proceedings of the 2015 IEEE International Conference on Industrial Technology (ICIT), Seville, Spain, 17–19 March 2015; pp. 3347–3350. [Google Scholar]
  70. Aghvami, S.A.; Opathalage, A.; Zhang, Z.; Ludwig, M.; Heymann, M.; Norton, M.; Wilkins, N.; Fraden, S. Rapid prototyping of cyclic olefin copolymer (COC) microfluidic devices. Sens. Actuators B Chem. 2017, 247, 940–949. [Google Scholar] [CrossRef]
  71. Ochs, C.J.; Kasuya, J.; Pavesi, A.; Kamm, R.D. Oxygen levels in thermoplastic microfluidic devices during cell culture. Lab Chip 2014, 14, 459–462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Santiago, G.T.-D.; Flores-Garza, B.G.; Negrete, J.T.; Lara-Mayorga, I.M.; González-Gamboa, I.; Zhang, Y.S.; Rojas-Martínez, A.; Ortiz-López, R.; Álvarez, M.M. The tumor-on-chip: Recent advances in the development of microfluidic systems to recapitulate the physiology of solid tumors. Materials 2019, 12, 2945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Hua, S.; De Matos, M.B.C.; Metselaar, J.M.; Storm, G. Current trends and challenges in the clinical translation of nanoparticulate nanomedicines: Pathways for translational development and commercialization. Front. Pharmacol. 2018, 9, 790. [Google Scholar] [CrossRef]
  74. Ahn, J.; Ko, J.; Lee, S.; Yu, J.; Kim, Y.; Jeon, N.L. Microfluidics in nanoparticle drug delivery; From synthesis to pre-clinical screening. Adv. Drug Deliv. Rev. 2018, 128, 29–53. [Google Scholar] [CrossRef]
  75. Soheili, S.; Mandegar, E.; Moradikhah, F.; Doosti-Telgerd, M.; Javar, H.A. Experimental and numerical studies on microfluidic preparation and engineering of chitosan nanoparticles. J. Drug Deliv. Sci. Technol. 2020, 61, 102268. [Google Scholar] [CrossRef]
  76. Arefi, S.M.A.; Yang, C.W.T.; Sin, D.D.; Feng, J.J. Simulation of nanoparticle transport and adsorption in a microfluidic lung-on-a-chip device. Biomicrofluidics 2020, 14, 044117. [Google Scholar] [CrossRef] [PubMed]
  77. Liu, Y.; Bhattacharya, A.; Kuksenok, O.; He, X.; Aizenberg, M.; Aizenberg, J.; Balazs, A.C. Computational modeling of oscillating fins that “catch and release” targeted nanoparticles in bilayer flows. Soft Matter 2016, 12, 1374–1384. [Google Scholar] [CrossRef] [PubMed]
  78. Maleki, R.; Khedri, M.; Malekahmadi, D.; Mohaghegh, S.; Jahromi, A.M.; Shahbazi, M.A. Simultaneous doxorubicin encap-sulation and in-situ microfluidic micellization of bio-targeted polymeric nanohybrids using dichalcogenide monolayers: A Molecular in-silico study. Mater. Today Commun. 2021, 26, 101948. [Google Scholar] [CrossRef]
  79. Kwak, B.; Ozcelikkale, A.; Shin, C.; Park, K.; Han, B. Simulation of complex transport of nanoparticles around a tumor using tumor-microenvironment-on-chip. J. Control. Release 2014, 194, 157–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Nasiri, R.; Shamloo, A.; Akbari, J.; Tebon, P.; Dokmeci, M.R.; Ahadian, S. Design and simulation of an integrated centrifugal microfluidic device for CTCs separation and cell lysis. Micromachines 2020, 11, 699. [Google Scholar] [CrossRef]
  81. Zhang, X.; Xu, X.; Ren, Y.; Yan, Y.; Wu, A. Numerical simulation of circulating tumor cell separation in a dielectrophoresis based Y-Y shaped microfluidic device. Sep. Purif. Technol. 2021, 255, 117343. [Google Scholar] [CrossRef]
  82. Ma, J.-T.; Xu, Y.-Q.; Tang, X.-Y. A numerical simulation of cell separation by simplified asymmetric pinched flow fractionation. Comput. Math. Methods Med. 2016, 2016, 2564584. [Google Scholar] [CrossRef] [Green Version]
  83. Shamloo, A.; Boodaghi, M. Design and simulation of a microfluidic device for acoustic cell separation. Ultrasonics 2018, 84, 234–243. [Google Scholar] [CrossRef]
  84. Karakas, H.E.; Kim, J.; Park, J.; Oh, J.M.; Choi, Y.; Gozuacik, D.; Cho, Y.-K. A microfluidic chip for screening individual cancer cells via eavesdropping on autophagy-inducing crosstalk in the stroma niche. Sci. Rep. 2017, 7, 2050. [Google Scholar] [CrossRef] [Green Version]
  85. Liu, J.-S.; Zhang, Y.-Y.; Wang, Z.; Deng, J.-Y.; Ye, X.; Xue, R.-Y.; Ge, D.; Xu, Z. Design and validation of a microfluidic chip with micropillar arrays for three-dimensional cell culture. Chin. J. Anal. Chem. 2017, 45, 1109–1114. [Google Scholar] [CrossRef]
  86. Chen, W.; Li, J.; Wan, X.; Zou, X.; Qi, S.; Zhang, Y.; Weng, Q.; Li, J.; Xiong, W.; Xie, C.; et al. Design of a microfluidic chip consisting of micropillars and its use for the enrichment of nasopharyngeal cancer cells. Oncol. Lett. 2019, 17, 1581–1588. [Google Scholar] [CrossRef] [PubMed]
  87. Huang, Y.; Nguyen, N.-T. A polymeric cell stretching device for real-time imaging with optical microscopy. Biomed. Microdevices 2013, 15, 1043–1054. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) CFD simulation of shear stress on the XY, YZ, and ZX planes. Adapted from [45]; (b) numerical simulation of H 2 O 2 concentration inside the microfluidic chip. Adapted from [46]; (c) contours of wall shear stress distribution in designed chambers. From left to right, central areas of each chamber, with different horizontal distances. The orange and red dashed lines indicate the uniform wall shear stress area boundary. Adapted from [47].
Figure 1. (a) CFD simulation of shear stress on the XY, YZ, and ZX planes. Adapted from [45]; (b) numerical simulation of H 2 O 2 concentration inside the microfluidic chip. Adapted from [46]; (c) contours of wall shear stress distribution in designed chambers. From left to right, central areas of each chamber, with different horizontal distances. The orange and red dashed lines indicate the uniform wall shear stress area boundary. Adapted from [47].
Micromachines 12 01149 g001
Figure 2. (a) Computational modeling of fluid flow and analyte transport in microchannels and porous membrane. White arrows indicate flow direction in the membrane and the upper microfluidic channel. The upper surface of the cell layer has a fixed concentration equal to 1 nM. Analyte flux streamlines are shown in black. Adapted from [36]; (b) chip design and glucose concentration profiles at a flow rate of 50 µL/h in the bare membrane, respectively. Adapted from [51].
Figure 2. (a) Computational modeling of fluid flow and analyte transport in microchannels and porous membrane. White arrows indicate flow direction in the membrane and the upper microfluidic channel. The upper surface of the cell layer has a fixed concentration equal to 1 nM. Analyte flux streamlines are shown in black. Adapted from [36]; (b) chip design and glucose concentration profiles at a flow rate of 50 µL/h in the bare membrane, respectively. Adapted from [51].
Micromachines 12 01149 g002
Figure 3. (a) Numerical simulation of on-chip temperature profiles including spatial temperature distribution (i), planar temperature distribution (ii), spatial temperature gradient distribution (iii), and planar temperature gradient distribution (iv). Adapted from [58]; (b) temperature distribution in the chip in the incubator environment (370 K ambient and bath at 900 K). Adapted from [59].
Figure 3. (a) Numerical simulation of on-chip temperature profiles including spatial temperature distribution (i), planar temperature distribution (ii), spatial temperature gradient distribution (iii), and planar temperature gradient distribution (iv). Adapted from [58]; (b) temperature distribution in the chip in the incubator environment (370 K ambient and bath at 900 K). Adapted from [59].
Micromachines 12 01149 g003
Figure 4. (a) Simulation of relatively stiff beads with a shear elastic modulus, Gs, of 10−3 N/m, compared to deformable cells Gs = 10−5 N/m. Adapted from [60]; (b) deformation of an RBC passing through a narrow microchannel. Adapted from [61]; (c) vector of velocity on the XY plane and position of the selected particles in experimental tests. Adapted from [63].
Figure 4. (a) Simulation of relatively stiff beads with a shear elastic modulus, Gs, of 10−3 N/m, compared to deformable cells Gs = 10−5 N/m. Adapted from [60]; (b) deformation of an RBC passing through a narrow microchannel. Adapted from [61]; (c) vector of velocity on the XY plane and position of the selected particles in experimental tests. Adapted from [63].
Micromachines 12 01149 g004
Figure 5. COC and PMMA at various average flow rates with a polymer capping layer with a thickness of 1.9 mm. At this thickness, for corresponding speeds, both polymers have virtually identical profiles. Adapted from [53].
Figure 5. COC and PMMA at various average flow rates with a polymer capping layer with a thickness of 1.9 mm. At this thickness, for corresponding speeds, both polymers have virtually identical profiles. Adapted from [53].
Micromachines 12 01149 g005
Figure 6. (a) Nanocarrier penetration into the cancer cell membrane without MoSe2 and with MoSe2, respectively. Adapted from [78]; (b) Distribution and instantaneous velocity of non-Brownian particles at t = 2, 4 and 6 s for particle diameter d = 100 nm and 500 nm. Adapted from [76].
Figure 6. (a) Nanocarrier penetration into the cancer cell membrane without MoSe2 and with MoSe2, respectively. Adapted from [78]; (b) Distribution and instantaneous velocity of non-Brownian particles at t = 2, 4 and 6 s for particle diameter d = 100 nm and 500 nm. Adapted from [76].
Micromachines 12 01149 g006
Figure 7. (a) Contour plot of the TGFβ1 exposure ( J TGF β 1 ) in (distance, exposure time)-space, where the numerical simulation was conducted. Adapted from [84]; (b) glucose concentration profile at different instants (1, 18, and 220 s) in the chip. Adapted from [85].
Figure 7. (a) Contour plot of the TGFβ1 exposure ( J TGF β 1 ) in (distance, exposure time)-space, where the numerical simulation was conducted. Adapted from [84]; (b) glucose concentration profile at different instants (1, 18, and 220 s) in the chip. Adapted from [85].
Micromachines 12 01149 g007
Table 1. Summary of the considerations made in the numerical studies.
Table 1. Summary of the considerations made in the numerical studies.
SoftwareFlow RegimePhases NumberFluid RheologyMeshTime-
Dependence
ValidationRef
COMSOLLaminar and turbulentSingle-phaseNewtoniann/dSteady×[45]
CFD-ACE +LaminarSingle-phasen/dn/dSteady[46]
AnsysLaminarSingle-phaseNewtonian78,000 nodesSteady×[47]
COMSOLn/dSingle-phasen/dn/dTransient[48]
COMSOLLaminarSingle-phaseNewtonian182,900 elementsTransient×[36]
n/dLaminarSingle-phasen/dn/dTransient[49]
IBM, MATLAB and AnsysLaminarMultiphaseNewtonianMesh independence test (1 µm)Transient×[50]
COMSOLLaminarSingle-phaseNewtonianMesh independence test (n/d)Transient[51]
COMSOLLaminarSingle-phaseNewtonianFine mesh (n/d)Transient×[4]
COMSOLn/dSingle-phasen/dn/dTransient[52]
COMSOLLaminarSingle-phasen/d7000 elementsSteady[53]
COMSOLn/dSingle-phasen/dMesh independence test (800,000 nodes)Transient×[54]
COMSOLn/dMultiphaseNewtonian24,300 elementsTransient[55]
User-Defined SoftwareLaminarSingle-phasen/dn/dSteady×[56]
COMSOLLaminarSingle-phasen/dn/dSteady×[58]
COMSOLLaminarSingle-phasen/d123,334 elementsSteady[59]
User-Defined Softwaren/dMultiphaseNewtoniann/dTransient[60]
SDPD and IBMn/dMultiphasen/dn/dTransient[61]
AnsysLaminarMultiphaseNewtonianMesh independence test (220,000 elements)Transient[62]
AnsysLaminarMultiphaseNewtoniann/dSteady[63]
COMSOLn/dSingle-phasen/d 1.8 × 10 6 elementsTransient[71]
COMSOLLaminarMultiphasen/dMesh independence test (6000 elements)Transient×[75]
n/dn/dMultiphasen/dn/dTransient[79]
n/dn/dMultiphasen/dn/dTransient×[77]
GROMACSn/dMultiphaseNewtoniann/dTransient[78]
COMSOLLaminarMultiphasen/dMesh independence test (168,000 elements)Steady and Transient[76]
COMSOLn/dSingle-phasen/dn/dTransient×[84]
n/dn/dMultiphasen/dn/dTransient×[81]
COMSOL and Ansysn/dSingle-phasen/dn/dTransient[85]
Ansysn/dSingle-phasen/dn/dSteady×[86]
Ansysn/dSingle-phasen/d2–8 μm element sizeTransient[87]
n/d—non-defined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Carvalho, V.; Rodrigues, R.O.; Lima, R.A.; Teixeira, S. Computational Simulations in Advanced Microfluidic Devices: A Review. Micromachines 2021, 12, 1149. https://doi.org/10.3390/mi12101149

AMA Style

Carvalho V, Rodrigues RO, Lima RA, Teixeira S. Computational Simulations in Advanced Microfluidic Devices: A Review. Micromachines. 2021; 12(10):1149. https://doi.org/10.3390/mi12101149

Chicago/Turabian Style

Carvalho, Violeta, Raquel O. Rodrigues, Rui A. Lima, and Senhorinha Teixeira. 2021. "Computational Simulations in Advanced Microfluidic Devices: A Review" Micromachines 12, no. 10: 1149. https://doi.org/10.3390/mi12101149

APA Style

Carvalho, V., Rodrigues, R. O., Lima, R. A., & Teixeira, S. (2021). Computational Simulations in Advanced Microfluidic Devices: A Review. Micromachines, 12(10), 1149. https://doi.org/10.3390/mi12101149

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop