Thermal compact modeling approach of droplet microreactor based Lab-on-a-Chip devices
Introduction
The rising complexity of Lab-on-a-Chip (LoC) devices for medical diagnostic purposes requires a multi-domain design. Sample handling is usually enabled by microfluidic techniques as the fluids are flowing in micro-sized channels where separation, mixing and bio-reactions may occur. The effects of such reactions can be detected by biosensors which provide electrical signals as output, more often required to be processed on-chip by integrated circuits. The method of co-design of MEMS and integrated circuits shown in Fig. 1 is already utilized in commercial products [1]. A MEMS device is usually modeled by numerical methods (e.g. FEM or CFD simulation) in order to generate a reduced order model (ROM). The ROM is considered as a behavioral description of the MEMS device therefore the electrical and the mechanical functionality can be handled together at a high level of abstraction. Numerical modeling usually provides high accuracy, and requires high computation times, though. In contrast, reduced order models focus on the key behavior of the device and require only moderate computation time. Therefore bypassing the numerical modeling step the design iteration time can be significantly reduced (Fig. 1). This approach can be generalized [2] especially for electro-mechanical microdevices, however the multi-physical nature of LoC devices requires special care as fluid dynamics, electrostatic/magnetic forces, chemical reactions and signal processing need to be handled at the same time. This paper introduces a reduced order thermal model for droplet microreactors which would enhance the design process of LoC devices.
The motivation for designing and manufacturing microreactors is that further reduction in sample volume can be achieved by using separate compartments e.g. liquid droplets encapsulating individual reactions.
The flow of droplets is a moving series of slugs of two phases which could be gas and liquid or two immiscible liquids. The slugs form monodisperse compartments in which species and reactions can be isolated and provide the ability to perform large numbers of experiments using extremely small sample volumes [3]. This technique enables the construction of a high throughput Lab-on-a-Chip platform. In addition to several other techniques, reaction heat measurement based biodetection seems to be one of the promising concepts in such micro-scale (chip size) laboratory devices [4].
Heat measurements constitute a direct means of determining enthalpy changes, which appear in all cases of protein–protein or protein–ligand interactions. The method referred to as calorimetry is capable of detecting antibody–antigen, protein–ligand or enzyme–ligand interactions at the bedside for any biomedically relevant case, such as cancer, neurological disorders, diabetes, and metabolic diseases. Although calorimetry is traditionally a low throughput method, the micro-scale realization enables higher reaction speeds and multiple measurements, therefore higher throughput can be achieved [5]. Continuous flow [6] and droplet calorimetry platforms are already reported [7] and demonstrated in various biological applications such as biomolecular characterization [8]. Thermal design is considered to be the key factor of sensitivity and resolution of such calorimetry based biosensors [9], [10]. The amount of thermal energy to be measured typically falls in the nano-Joule range. Lee et. al. successfully demonstrated the enthalpy measurement of urease enzyme activity with a resolution of 10 nJ in a microfluidic calorimeter with a reactor volume of 3.5 nL [10]. Their measurement results fitted very well with the simple theory of Michaelis–Menten kinetics. Droplet size, reaction kinetics, droplet velocity, material of the channel etc. affect the heat transfer from the droplet towards the ambient. These properties should be considered together with the signal processing solutions (e.g. on-chip integrated analog or digital circuits) at behavioral level. Thermal compact models are widely used in design practice, where such reduced order models (in contrast to detailed numerical models) yield results quickly in the early design phase.
Section snippets
Modeling
Heat transfer models and investigations are traditionally based on the determination of the Nusselt number, which is obviously practical when the heat transfer properties of the individual droplets are neglected and an overall description of the flow is possible. This is widely used e.g. for modeling fluids in electronic cooling applications such as cooling integrated circuits. In contrast, for bio-analytical investigations the thermal interactions should be analyzed in detail as they happen in
Simulation settings
Four simulations were performed by the ANSYS FLUENT (CFD) tool and by our novel compact thermal model (CTM) with constant heat flux (CHF) and with heat generation by enzyme reaction (ENZ) boundary conditions. Further five simulations were done for the model validity range analysis in CHF case. Material properties and simulation parameters are summarized in Table 1 and Table 2, respectively. The CFD simulation is built according to the experimentally validated recommendations of Gupta et. al.
Performance analysis
For the performance test the linear solver of ANSYS 14 was run on an Intel Core i5 4 core CPU with 4 GB RAM under 64 bit Debian GNU/Linux (kernel 3.2.0-2-amd64). The solver has a parallel computing capability but parallelism did not improve the solution speed as the steps must be processed sequentially. In order to achieve reasonable results an axial resolution of was applied for a channel diameter of which yields a sum of 2736 nodes (radial mesh resolution ranges from to .
Conclusion
Based on the general properties of the segmented flow, a simplified model was introduced to describe heat transfer in microchannels. Thermal analysis was performed by a novel compact model and the results were validated by ANSYS FLUENT. The presented compact model provided a 450 fold reduction in execution with little compromise for accuracy. Relative errors of the temperature profiles were calculated in order to compare the results of the two methods. The relative errors of the temperature
References (18)
- et al.
Higher throughput calorimetryopportunities, approaches and challenges
Curr. Opin. Struct. Biol.
(2010) - et al.
Chip calorimetry and its use for biochemical and cell biological investigations
Chem. Eng. Process.: Process Intensif.
(2008) - et al.
Enthalpy array analysis of enzymatic and binding reactions
Anal. Biochem.
(2008) - et al.
Temperature-modulated differential scanning calorimetry in a MEMS device
Sens. Actuators B: Chem.
(2013) - et al.
Chip calorimeters for the investigation of liquid phase reactionsdesign rules
Thermochim. Acta
(2006) - et al.
CFD modelling of flow and heat transfer in the Taylor flow regime
Chem. Eng. Sci.
(2010) - et al.
CFD approaches for the simulation of hydrodynamics and heat transfer in Taylor flow
Chem. Eng. Sci.
(2011) - et al.
Validation of a CFD model of Taylor flow hydrodynamics and heat transfer
Chem. Eng. Sci.
(2012) - et al.
Catalytic behavior of carbonic anhydrase enzyme immobilized onto nonporous silica nanoparticles for enhancing CO2 absorption into a carbonate solution
Int. J. Greenh. Gas Control
(2013)