Numerical Study of Atrial Fibrillation Effects on Flow Distribution in Aortic Circulation

Role of the vessel morphology on the lenticulostriate arteries hemodynamics during atrial fibrillation: A CFD-based multivariate regression analysis

BACKGROUND AND OBJECTIVE

Atrial fibrillation (AF) is the most common cardiac arrhythmia, inducing accelerated and irregular beating. Beside well-known disabling symptoms - such as palpitations, reduced exercise tolerance, and chest discomfort - there is growing evidence that an alteration of deep cerebral hemodynamics due to AF increases the risk of vascular dementia and cognitive impairment, even in the absence of clinical strokes. The alteration of deep cerebral circulation in AF represents one of the least investigated among the possible mechanisms. Lenticulostriate arteries (LSAs) are small perforating arteries mainly departing from the middle cerebral artery (MCA) and susceptible to small vessel disease, which is one of the mechanisms of subcortical vascular dementia development. The purpose of this study is to investigate the impact of different LSAs morphologies on the cerebral hemodynamics during AF.

METHODS

By combining a computational fluid dynamics (CFD) analysis of LSAs with 7T high-resolution magnetic resonance imaging (MRI), we performed different CFD-based multivariate regression analyses to detect which geometrical and morphological vessel features mostly affect AF hemodynamics in terms of wall shear stress. We exploited 17 cerebral 7T-MRI derived LSA vascular geometries extracted from 10 subjects and internal carotid artery data from validated 0D cardiovascular-cerebral modeling as inflow conditions.

RESULTS

Our results revealed that few geometrical variables - namely the size of the MCA and the bifurcation angles between MCA and LSA - are able to satisfactorily predict the AF impact. In particular, the present study indicates that LSA morphologies exhibiting markedly obtuse LSA-MCA inlet angles and small MCA size downstream of the LSA-MCA bifurcation may be more prone to vascular damage induced by AF.

CONCLUSIONS

The present MRI-based computational study has been able for the first time to: (i) investigate the net impact of LSAs vascular morphologies on cerebral hemodynamics during AF events; (ii) detect which combination of morphological features worsens the hemodynamic response in the presence of AF. Awaiting necessary clinical confirmation, our analysis suggests that the local hemodynamics of LSAs is affected by their geometrical features and some LSA morphologies undergo greater hemodynamic alterations in the presence of AF.

Enhancing the implantation of mechanical circulatory support devices using computational simulations

Introduction: Patients with end-stage heart failure (HF) may need mechanical circulatory support such as a left ventricular assist device (LVAD). However, there are a range of complications associated with LVAD including aortic regurgitation (AR) and thrombus formation. This study assesses whether the risk of developing aortic conditions can be minimised by optimising LVAD implantation technique. Methods: In this work, we evaluate the aortic flow patterns produced under different geometrical parameters for the anastomosis of the outflow graft (OG) to the aorta using computational fluid dynamics (CFD). A three-dimensional aortic model is created and the HeartMate III OG positioning is simulated by modifying (i) the distance from the anatomic ventriculo-arterial junction (AVJ) to the OG, (ii) the cardinal position around the aorta, and (iii) the angle between the aorta and the OG. The continuous LVAD flow and the remnant native cardiac cycle are used as inlet boundaries and the three-element Windkessel model is applied at the pressure outlets. Results: The analysis quantifies the impact of OG positioning on different haemodynamic parameters, including velocity, wall shear stress (WSS), pressure, vorticity and turbulent kinetic energy (TKE). We find that WSS on the aortic root (AoR) is around two times lower when the OG is attached to the coronal side of the aorta using an angle of 45° ± 10° at a distance of 55 mm. Discussion: The results show that the OG placement may significantly influence the haemodynamic patterns, demonstrating the potential application of CFD for optimising OG positioning to minimise the risk of cardiovascular complications after LVAD implantation.

Smoothed particle hydrodynamics based FSI simulation of the native and mechanical heart valves in a patient-specific aortic model

The failure of the aortic heart valve is common, resulting in deterioration of the pumping function of the heart. For the end stage valve failure, bi-leaflet mechanical valve (most popular artificial valve) is implanted. However, due to its non-physiological behaviour, a significant alteration is observed in the normal haemodynamics of the aorta. While in-vivo experimentation of a human heart valve (native and artificial) is a formidable task, in-silico study using computational fluid dynamics (CFD) with fluid structure interaction (FSI) is an effective and economic tool for investigating the haemodynamics of natural and artificial heart valves. In the present work, a haemodynamic model of a natural and mechanical heart valve has been developed using meshless particle-based smoothed particle hydrodynamics (SPH). In order to further enhance its clinical relevance, this study employs a patient-specific vascular geometry and presents a successful validation against traditional finite volume method and 4D magnetic resonance imaging (MRI) data. The results have demonstrated that SPH is ideally suited to simulate the heart valve function due to its Lagrangian description of motion, which is a favourable feature for FSI. In addition, a novel methodology for the estimation of the wall shear stress (WSS) and other related haemodynamic parameters have been proposed from the SPH perspective. Finally, a detailed comparison of the haemodynamic parameters has been carried out for both native and mechanical aortic valve, with a particular emphasis on the clinical risks associated with the mechanical valve.

Exergy destruction in atrial fibrillation and a new 'Exergy Age Index'

The concept of exergy in living organisms has been widely used to explore correlations between exergy and different physiological conditions. Atrial fibrillation (AF) is an abnormal physiological condition that takes place inside the heart and is recognised as a common supraventricular arrhythmia. AF can significantly undermine heart function and subsequently circulatory system. Thus, exergy analysis of cardiac flow during AF is a procedure to quantify the long-term impact of persistent AF. The present study adopts the lumped modelling approach for considering cardiovascular circulation and thermoregulation of the body to evaluate the exergy consumption and destruction of the heart in AF. In order to assess the impact of AF, four common AF-associated characteristics including lack of atrial kick, left atrial remodelling, left ventricular systolic dysfunction, and high-frequency fibrillation are examined. The results show that among AF deficiencies, high-frequency fibrillation is the main cause of exergy destruction of the heart during AF. Moreover, a novel 'exergy age index' is proposed which has shown that high fibrillatory conditions in AF can significantly accelerate the heart ageing process, which in turn substantiates AF's adverse impact on the heart.

The Differences of Quantitative Flow Ratio in Coronary Artery Stenosis with or without Atrial Fibrillation

Quantitative flow ratio (QFR) is a new method for the assessment of the extent of coronary artery stenosis. But it may be obscured by the cardiac remodeling and abnormal blood flow of the coronary artery when encountering atrial fibrillation (AF). The present study aimed to examine the impact of these changed structures and blood flow of coronary arteries on QFR results in AF patients. Methods and Results. We evaluated QFR in 223 patients (112 patients with AF; 111 non-AF patients served as controls) who had undergone percutaneous coronary intervention (PCI) due to severe stenoses in coronary arteries. QFR of the target coronary was determined according to the flow rate of the contrast agent. Results showed that AF patients had significantly higher QFR values than control (0.792 ± 0.118 vs. 0.685 ± 0.167, p < 0.001). We further analyzed local QFR around the stenoses (0.858 ± 0.304 vs. 0.756 ± 0.146, p=0.002), residual QFR (0.958 ± 0.055 vs. 0.929 ± 0.093, p=0.005), and index QFR (0.807 ± 0.108 vs. 0.713 ± 0.152, p < 0.001) in these two groups of patients with and without AF. Further analysis revealed that QFR in AF patients was negatively correlated with coronary flow velocity (R = -0.22, p=0.02) and area of stenosis (R = -0.70, p < 0.001) but positively correlated with the minimum lumen area (MLA) (R = 0.47, p < 0.001). Conclusion. AF patients with coronary artery stenosis have higher QFR values, which are associated with decreased blood flow velocity, smaller stenosis, and larger MLA in AF patients upon cardiac remodeling.

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia

Introduction

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH.

Methods

Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TAR_NO).

Results

Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm³/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa^-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa^-1 vs. 1.044, p = 0.001) and deceased TAR_NO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs.

Conclusion

Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TAR_NO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts.

Imaging and biophysical modelling of thrombogenic mechanisms in atrial fibrillation and stroke

Atrial fibrillation (AF) underlies almost one third of all ischaemic strokes, with the left atrial appendage (LAA) identified as the primary thromboembolic source. Current stroke risk stratification approaches, such as the CHA₂DS₂-VASc score, rely mostly on clinical comorbidities, rather than thrombogenic mechanisms such as blood stasis, hypercoagulability and endothelial dysfunction-known as Virchow's triad. While detection of AF-related thrombi is possible using established cardiac imaging techniques, such as transoesophageal echocardiography, there is a growing need to reliably assess AF-patient thrombogenicity prior to thrombus formation. Over the past decade, cardiac imaging and image-based biophysical modelling have emerged as powerful tools for reproducing the mechanisms of thrombogenesis. Clinical imaging modalities such as cardiac computed tomography, magnetic resonance and echocardiographic techniques can measure blood flow velocities and identify LA fibrosis (an indicator of endothelial dysfunction), but imaging remains limited in its ability to assess blood coagulation dynamics. In-silico cardiac modelling tools-such as computational fluid dynamics for blood flow, reaction-diffusion-convection equations to mimic the coagulation cascade, and surrogate flow metrics associated with endothelial damage-have grown in prevalence and advanced mechanistic understanding of thrombogenesis. However, neither technique alone can fully elucidate thrombogenicity in AF. In future, combining cardiac imaging with in-silico modelling and integrating machine learning approaches for rapid results directly from imaging data will require development under a rigorous framework of verification and clinical validation, but may pave the way towards enhanced personalised stroke risk stratification in the growing population of AF patients. This Review will focus on the significant progress in these fields.

Cerebral hemodynamics during atrial fibrillation: Computational fluid dynamics analysis of lenticulostriate arteries using 7 T high-resolution magnetic resonance imaging

Atrial fibrillation (AF) is the most common cardiac arrhythmia, inducing irregular and faster heart beating. Aside from disabling symptoms-such as palpitations, chest discomfort, and reduced exercise capacity-there is growing evidence that AF increases the risk of dementia and cognitive decline, even in the absence of clinical strokes. Among the possible mechanisms, the alteration of deep cerebral hemodynamics during AF is one of the most fascinating and least investigated hypotheses. Lenticulostriate arteries (LSAs)-small perforating arteries perpendicularly departing from the anterior and middle cerebral arteries and supplying blood flow to basal ganglia-are especially involved in silent strokes and cerebral small vessel diseases, which are considered among the main vascular drivers of dementia. We propose for the first time a computational fluid dynamics analysis to investigate the AF effects on the LSAs hemodynamics by using 7 T high-resolution magnetic resonance imaging (MRI). We explored different heart rates (HRs)-from 50 to 130 bpm-in sinus rhythm and AF, exploiting MRI data from a healthy young male and internal carotid artery data from validated 0D cardiovascular-cerebral modeling as inflow condition. Our results reveal that AF induces a marked reduction of wall shear stress and flow velocity fields. This study suggests that AF at higher HR leads to a more hazardous hemodynamic scenario by increasing the atheromatosis and thrombogenesis risks in the LSAs region.

Insights from computational modeling on the potential hemodynamic effects of sinus rhythm versus atrial fibrillation

Atrial fibrillation (AF) is the most common clinical tachyarrhythmia, posing a significant burden to patients, physicians, and healthcare systems worldwide. With the advent of more effective rhythm control strategies, such as AF catheter ablation, an early rhythm control strategy is progressively demonstrating its superiority not only in symptoms control but also in prognostic terms, over a standard strategy (rate control, with rhythm control reserved only to patients with refractory symptoms). This review summarizes the different impacts exerted by AF on heart mechanics and systemic circulation, as well as on cerebral and coronary vascular beds, providing computational modeling-based hemodynamic insights in favor of pursuing sinus rhythm maintenance in AF patients.

A Finite Element Based Analysis of a Hemodynamics Efficient Flow Stent Suitable for Different Abdominal Aneurysm Shapes

This research aimed to examine the impact of a proposed flow stent (PFS) on different abdominal artery shapes. For that purpose, a finite element-based model using the computational fluid dynamics (CFD) method is developed. The effect of PFS intervention on the hemodynamic efficiency is estimated by all of the significant criteria used for the evaluation of aneurysm occlusion and possible rupture; the flow velocity, pressure, wall shear stress (WSS), and WSS-related indices. Results showed that PFS intervention preserves the effects of high flow rate and decreases irregular flow recirculation in the sac of the aneurysm. The flow velocity decreases inside the aneurysm sac in the range of 55% to 80%. The time-averaged wall shear stress (TAWSS) was reduced from 42% to 53% by FPS deployment. The simulation results implies that PFS could heal an aneurysm efficiently with a mechanism that causes the development of thrombus and ultimately leads to aneurysm resorption.

Development of a computational fluid dynamic model to investigate the hemodynamic impact of REBOA

Background: Resuscitative endovascular balloon occlusion of the aorta (REBOA) is a lifesaving intervention for major truncal hemorrhage. Balloon-tipped arterial catheters are inserted via the femoral artery to create a temporary occlusion of the aorta, which minimizes the rate of internal bleeding until definitive surgery can be conducted. There is growing concern over the resultant hypoperfusion and potential damage to tissues and organs downstream of REBOA. To better understand the acute hemodynamic changes imposed by REBOA, we developed a three-dimensional computational fluid dynamic (CFD) model under normal, hemorrhage, and aortic occlusion conditions. The goal was to characterize the acute hemodynamic changes and identify regions within the aortic vascular tree susceptible to abnormal flow and shear stress. Methods: Hemodynamic data from established porcine hemorrhage models were used to build a CFD model. Swine underwent 20% controlled hemorrhage and were randomized to receive a full or partial aortic occlusion. Using CT scans, we generated a pig-specific aortic geometry and imposed physiologically relevant inlet flow and outlet pressure boundary conditions to match in vivo data. By assuming non-Newtonian fluid properties, pressure, velocity, and shear stresses were quantified over a cardiac cycle. Results: We observed a significant rise in blood pressure (∼147 mmHg) proximal to REBOA, which resulted in increased flow and shear stress within the ascending aorta. Specifically, we observed high levels of shear stress within the subclavian arteries (22.75 Pa). Alternatively, at the site of full REBOA, wall shear stress was low (0.04 ± 9.07E-4 Pa), but flow oscillations were high (oscillatory shear index of 0.31). Comparatively, partial REBOA elevated shear levels to 84.14 ± 19.50 Pa and reduced flow oscillations. Our numerical simulations were congruent within 5% of averaged porcine experimental data over a cardiac cycle. Conclusion: This CFD model is the first to our knowledge to quantify the acute hemodynamic changes imposed by REBOA. We identified areas of low shear stress near the site of occlusion and high shear stress in the subclavian arteries. Future studies are needed to determine the optimal design parameters of endovascular hemorrhage control devices that can minimize flow perturbations and areas of high shear.

In Silico Hemodynamics and Filtering Evaluation of a Commercial Embolic Protection Device

During the last years, several kinds of Embolic Protection Devices (EPD) have been developed, with the aim of minimizing complication caused by thrombi generated during Carotid Artery Stenting (CAS). These devices are capable of capturing small particles generated during the intervention, avoiding cerebral stroke and improving the outcomes of the surgery. However, they have associated complications, like the increase on flow resistance associated by their use or the lack of knowledge on their actual filtration efficiency for thrombi of low size. Current work proposes a validated computational methodology in order to predict the hemodynamic features and filtering efficiency of a commercial EPD. It will be observed how Computational Fluid Dynamics predicts pressure drop with fair agreement with the experimental measurements. Finally, this work analyzes the filtration efficiency and the influence of the distribution of injected particles on this parameter. The capabilities of the filter for retaining particles of diameter below the pore size is, additionally, discussed.

Feasibility and validation of trans-valvular flow derived by four-dimensional flow cardiovascular magnetic resonance imaging in patients with atrial fibrillation

Background: Four-dimensional (4D) flow cardiovascular magnetic resonance imaging (MRI) is an emerging technique used for intra-cardiac blood flow assessment. The role of 4D flow cardiovascular MRI in the assessment of trans-valvular flow in patients with atrial fibrillation (AF) has not previously been assessed. The purpose of this study was to assess the feasibility, image quality, and internal validity of 4D flow cardiovascular MRI in the quantification of trans-valvular flow in patients with AF. Methods: Patients with AF and healthy controls in sinus rhythm underwent cardiovascular MRI, including 4D flow studies. Quality assurance checks were done on the raw data and streamlines. Consistency was investigated by trans-valvular flow assessment between the mitral valve (MV) and the aortic valve (AV). Results: Eight patients with AF (88% male, mean age 62±13 years, mean heart rate (HR) 83±16 beats per minute (bpm)) were included and compared with ten healthy controls (70% male, mean age 41±20 years, mean HR 68.5±9 bpm). All scans were of either good quality with minimal blurring artefacts, or excellent quality with no artefacts. No significant bias was observed between the AV and MV stroke volumes in either healthy controls (-4.8, 95% CI -15.64 to 6.04; P=0.34) or in patients with AF (1.64, 95% CI -4.7 to 7.94; P=0.56). A significant correlation was demonstrated between MV and AV stroke volumes in both healthy controls (r=0.87, 95% CI 0.52 to 0.97; P=0.001) and in AF patients (r=0.82, 95% CI 0.26 to 0.97; P=0.01). Conclusions: In patients with AF, 4D flow cardiovascular MRI is feasible with good image quality, allowing for quantification of trans-valvular flow.

Feasibility and validation of trans-valvular flow derived by four-dimensional flow cardiovascular magnetic resonance imaging in patients with atrial fibrillation

Background: Four-dimensional (4D) flow cardiovascular magnetic resonance imaging (MRI) is an emerging technique used for intra-cardiac blood flow assessment. The role of 4D flow cardiovascular MRI in the assessment of trans-valvular flow in patients with atrial fibrillation (AF) has not previously been assessed. The purpose of this study was to assess the feasibility, image quality, and internal validity of 4D flow cardiovascular MRI in the quantification of trans-valvular flow in patients with AF. Methods: Patients with AF and healthy controls in sinus rhythm underwent cardiovascular MRI, including 4D flow studies. Quality assurance checks were done on the raw data and streamlines. Consistency was investigated by trans-valvular flow assessment between the mitral valve (MV) and the aortic valve (AV). Results: Eight patients with AF (88% male, mean age 62±13 years, mean heart rate (HR) 83±16 beats per minute (bpm)) were included and compared with ten healthy controls (70% male, mean age 41±20 years, mean HR 68.5±9 bpm). All scans were of either good quality with minimal blurring artefacts, or excellent quality with no artefacts. No significant bias was observed between the AV and MV stroke volumes in either healthy controls (-4.8, 95% CI -15.64 to 6.04; P=0.34) or in patients with AF (1.64, 95% CI -4.7 to 7.94; P=0.56). A significant correlation was demonstrated between MV and AV stroke volumes in both healthy controls (r=0.87, 95% CI 0.52 to 0.97; P=0.001) and in AF patients (r=0.82, 95% CI 0.26 to 0.97; P=0.01). Conclusions: In patients with AF, 4D flow cardiovascular MRI is feasible with good image quality, allowing for quantification of trans-valvular flow.

Effects of Ageing on Aortic Circulation During Atrial Fibrillation; a Numerical Study on Different Aortic Morphologies

Atrial fibrillation (AF) can alter intra-cardiac flow and cardiac output that subsequently affects aortic flow circulation. These changes may become more significant where they occur concomitantly with ageing. Aortic ageing is accompanied with morphological changes such as dilation, lengthening, and arch unfolding. While the recognition of AF mechanism has been the subject of numerous studies, less focus has been devoted to the aortic circulation during the AF and there is a lack of such investigation at different ages. The current work aims to address the present gap. First, we analyse aortic flow distribution in three configurations, which attribute to young, middle and old people, using geometries constructed via clinical data. We then introduce two transient inlet flow conditions representative of key AF-associated defects. Results demonstrate that both AF and ageing negatively affect flow circulation. The main consequence of concomitant occurrence is enhancement of endothelial cell activation potential (ECAP) throughout the vascular domain, mainly at aortic arch and descending thoracic aorta, which is consistent with some clinical observations. The outcome of the current study suggests that AF exacerbates the vascular defects occurred due to the ageing, which increases the possibility of cardiovascular diseases per se.

A computational simulation platform for designing real-time monitoring systems with application to COVID-19

With the aim of contributing to the fight against the coronavirus disease 2019 (COVID-19), numerous strategies have been proposed. While developing an effective vaccine can take months up to years, detection of infected patients seems like one of the best ideas for controlling the situation. The role of biosensors in containing highly pathogenic viruses, saving lives and economy is evident. A new competitive numerical platform specifically for designing microfluidic-integrated biosensors is developed and presented in this work. Properties of the biosensor, sample, buffer fluid and even the microfluidic channel can be modified in this model. This feature provides the scientific community with the ability to design a specific biosensor for requested point-of-care (POC) applications. First, the validation of the presented numerical platform against experimental data and then results and discussion, highlighting the important role of the design parameters on the performance of the biosensor is presented. For the latter, the baseline case has been set on the previous studies on the biosensors suitable for SARS-CoV, which has the highest similarity to the 2019 nCoV. Subsequently, the effects of concentration of the targeted molecules in the sample, installation position and properties of the biosensor on its performance were investigated in 11 case studies. The presented numerical framework provides an insight into understanding of the virus reaction in the design process of the biosensor and enhances our preparation for any future outbreaks. Furthermore, the integration of biosensors with different devices for accelerating the process of defeating the pandemic is proposed.

Numerical framework for simulating bio-species transport in microfluidic channels with application to antibody biosensors

Diagnosis is a fundamental stage in health care and medical treatment. Microfluidic biosensors and lab-on-a-chip devices are amongst the few practical tools for achieving this goal. A new computational code, specifically for designing microfluidic-integrated biosensors is developed, the details of which is presented in this work. This new approach is developed using control-volume based finite-element (CVFEM) method and solves bio-recognition chemical reactions and full Navier–Stokes equations. The results of the proposed platform are validated against the experimental data for a microfluidic based biosensor, where excellent agreement is achieved. The properties of the biosensor, sample, buffer fluid and even the microfluidic channel can easily be modified in this platform. This feature provides the scientific community with the ability to design a specific biosensor for requested point-of-care applications.•

A new approach is developed using control-volume based finite-element (CVFEM) method for investigating flow inside a microfluidic-integrated biosensor. It is also used to study the influence of surface functionalization on binding cycle.

•

The proposed model solves bio-recognition chemical reactions as well as full Navier–Stokes and energy equations. Experimental-based or personalized equations of the chemical reactions and flow behaviour are adoptable to this code.

•

The developed model is Fortran-based and has the potential to be used in both industry and academia for biosensing technology.