Next Article in Journal
Pushing the Limits of Spatial Assay Resolution for Paper-Based Microfluidics Using Low-Cost and High-Throughput Pen Plotter Approach
Previous Article in Journal
Tissue-Engineered Models for Glaucoma Research
Previous Article in Special Issue
Cell Lysis Based on an Oscillating Microbubble Array
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design Applicable 3D Microfluidic Functional Units Using 2D Topology Optimization with Length Scale Constraints

1
Changchun Institute of Optics, Fine Mechanics and Physics (CIOMP), Chinese Academy of Science, Changchun 130033, China
2
University of Chinese Academy of Science, Beijing 100049, China
3
Department of Computing Science, Umeå University, SE-901 87 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Micromachines 2020, 11(6), 613; https://doi.org/10.3390/mi11060613
Submission received: 15 November 2019 / Revised: 1 June 2020 / Accepted: 13 June 2020 / Published: 24 June 2020
(This article belongs to the Special Issue Optimization of Microfluidic Devices)

Abstract

:
Due to the limits of computational time and computer memory, topology optimization problems involving fluidic flow frequently use simplified 2D models. Extruded versions of the 2D optimized results typically comprise the 3D designs to be fabricated. In practice, the depth of the fabricated flow channels is finite; the limited flow depth together with the no-slip condition potentially make the fluidic performance of the 3D model very different from that of the simplified 2D model. This discrepancy significantly limits the usefulness of performing topology optimization involving fluidic flow in 2D—at least if special care is not taken. Inspired by the electric circuit analogy method, we limit the widths of the microchannels in the 2D optimization process. To reduce the difference of fluidic performance between the 2D model and its 3D counterpart, we propose an applicable 2D optimization model, and ensure the manufacturability of the obtained layout, combinations of several morphology-mimicking filters impose maximum or minimum length scales on the solid phase or the fluidic phase. Two typical Lab-on-chip functional units, Tesla valve and fluidic channel splitter, are used to illustrate the validity of the proposed application of length scale control.

1. Introduction

With the trend of miniaturization in recent years, Lab-on-a-chip devices have been widely used in the area of the analysis, synthesis, and separations of biofluidics due to the advantages of high efficiency, portability, and low reagent consumption [1]. Many functions of the conventional analytical laboratory can be achieved on a centimeter-level chip, such as injection, mixing, reaction, cleaning, separation, and detection. In recent years, various microfluidic devices have been designed and applied, such as micropumps, microvalves, micromixers, and microchannels [2,3,4,5]. To pursue more functions and higher efficiency, the design and optimization of microfluidic devices is a continuously studied topic.
Borrvall et al. [6] performed topology optimization for fluidic flow governed by the Stokes equation. Their approach was later generalized to the Navier–Stokes equations by Gersborg–Hansen et al. [7] and Olesen et al. [8]. Topology optimization methods have been used to successfully design lab-on-a-chip microfluidic devices, such as microchannel splitters, micropumps, no-moving-part microvalves, and micromixers [9,10,11,12,13,14,15]. The topology optimization model of fluidic flow usually uses a simplified 2D problem. Extruded versions of the 2D optimization results typically comprise the 3D designs to be fabricated. We remark that if the 2D design is stretched so that the 3D flow path has infinite depth, then the full 3D problem describing the fluidic flow reduces to the 2D problem used in the optimization. However, in practice, the depth of the flow channel is not infinite. The limited flow depth and the no-slip condition, which states that the flow speed is zero along the walls, together make the fluidic performance of the 3D model very different from that of the 2D model. This discrepancy significantly limits the application of the obtained 2D results. Some simplification methods have been proposed by build pseudo-3D models in recent years [16,17,18,19], and these research works are mostly focus on heat sinks problem. Zhao et al. [20] impose length scale control on the cooling channel. Based on our previous work [9,14], the length scale control methods are used to design the microfluidic devices with physical constraints (Diodicity for the example of Tesla valve, or the ratio of flowrate at the outlet for the example of splitter) in this paper.
Electric circuit analogy is a method of calculating and designing microchannel networks by analogizing the foundation parameters of the fluid to those of an electrical circuit [21]. When the height of the channel is constant, the hydraulic resistance is linear with the length and non-linear with the width of the channel. Inspired by this, we limit the widths of the microchannels in the optimization process to reduce the difference between the fluidic performance of the 2D and 3D model. Combinations of morphology mimicking filters [22,23,24] impose maximum or minimum length scales on the solid phase or the fluidic phase. The length scale control also ensures the manufacturability of the obtained layout, which avoids unreasonable small-sized islands inside the fluidic channel or thin-sized width of the fluidic channel. This paper is organized as follows—the topology optimization method and electric circuit analogy method are introduced in Section 2; the length scale control method using the morphology-mimicking filters and several numerical results are presented in Section 3; the discussion and conclusion are stated in Section 4.

2. Problem Statement

2.1. Fluid Flow Model

Based on the continuity assumption, Navier–Stokes equations describe the fluid flow
ρ ( u · ) u η Δ u + p = f , · u = 0 ,
where f is the body forces acting on the fluid and u , p, ρ , and η are the velocity, pressure, density, and viscosity of the fluid, respectively. To model the fluid–solid interface, we let f be an artificial friction force, similarly as proposed for the Stokes flow by Borrvall et al. [6]. More precisely, the artificial friction force f = α ( γ ) u , where α is the impermeability of the artificial porous material and γ is a material indicator function. The function γ varies in the interval [0, 1], where 0 and 1 denote the solid and fluid phase, respectively. The impermeability of the porous material is the interpolation function of the function γ :
α ( γ ) = α max q ( 1 γ ) q + γ ,
where α max is the impermeability of the solid phase and q is a positive value used to adjust the convexity of the interpolation function. To obtain perfect impermeability of the solid no-slip boundary, α max should be infinite; but a finite number has to be chosen to ensure numerical stability. For the topology optimization of fluid, the aim of the optimization is usually to minimize
Φ = Ω 2 η ϵ ( u ) : ϵ ( u ) + α ( γ ) | u | 2 d Ω ,
where ϵ ( u ) = u + u T / 2 , | | is L 2 norm. That is, the typical aim is to minimize the viscous dissipation inside the computational domain Ω .

2.2. Electric Circuit Analogy Method

The electric circuit analogy method provides a fast calculation of laminar flow in microchannels based on the analogous behavior of hydraulic and electric circuits with correlations of pressure to voltage, volumetric flowrate to current, and hydraulic resistance to electric resistance [21]. The hydraulic resistance is an important parameter in the design of microfluidic networks by electric circuit analogy. The hydraulic resistance R H of a rectangular microchannel section can be calculated as [25]
R H = 12 η L w h 3 1 h w 192 π 5 n = 1 , 3 , 5 1 n 5 tanh n π w 2 h ,
where η is the viscosity of the fluid and w, h, and L are the width, height, and length of the microchannel section, respectively. Figure 1 illustrates the hydraulic resistance changes with the length (red), height (blue), and width (green) of the channel, where the viscosity of the fluid η is normalized as 1. When the height of the channel is constant, the hydraulic resistance changes linearly with the length and nonlinearly with the width of the channel. By limiting the width of the flow channel, the hydraulic resistance changes linearly with the length. The hydraulic resistance ratio between different parts of the fluid network can remain unchanged in different stretching heights, so the difference of fluidic performance between the 2D model and its 3D counterpart is small. In order to preserve the flexibility of topology optimization, the changes of the channel width is constrained within a certain range.

3. Numerical Experiments

To demonstrate the validity of the proposed method, which can reduce the difference of fluidic performance between the 2D model and its 3D counterpart by using the size control of the widths of the fluidic channels, two examples are used in this paper. The viscosity and density of the fluid are 10 3 Pa · s and 998 kg · m 3 , respectively. The inlet velocity is chosen to have a parabolic profile. The intermediate fluid velocity in the height direction at the entrance of the 3D model is selected to be consistent with the 2D model’s inlet velocity. For all numerical experiments, Svanberg’s method of moving asymptotes [26] solves the resulting optimization problem.
We discretize our computational domain Ω by N quadrilateral elements. Comsol multiphysics solves the governing Navier–Stokes equation using bi-linear tensor product elements together with GLS (Galerkin least squares) streamline diffusion with the crosswind diffusion parameter set to 0.1.
We let M denote the number of nodal points in the discretization of the design domain Ω D . We define our material indicator function γ to be an nodal-wise bi-linear function whose restriction to Ω \ Ω D satisfies γ Ω \ Ω D 1 . The M nodal values of the material indicator function in Ω D are obtained by using a harmonic mean based open–close f W -mean filter. For a comprehensive review of mathematical morphology, we refer the reader to the Heijmans article [27]. Figure 2 illustrates the four-basic morphological operators: dilation, erosion, opening, and closing. In colloquial terms, given a brush with shape B, then the opening of S by B holds all points the can be filled by using this brush without touching any point outside S; and given a fully filled space and an eraser with shape B, the closing of S is the region that remains after having erased as much as possible without removing any point in S. Or, in other words, applying the opening operator to a set S removes all small parts from S; similarly, applying the closing operator to a set S removes all small internals holes. The region obtained by the open operation can be written as a union of translates of B, and the complement of the region obtained by the close operation can be written as a union of translates of B. It can be shown that the minimum size of O B ( S ) as well as R d \ C B ( S ) is at least the size of B.
The use of a cascade of filters to approximate the open–close operator was suggested by Sigmund [28], and harmonic mean based filters were introduced by Svanberg et al. [29]. Using the notation for f W -filters [30], we define the equally weighted discrete harmonic erode and dilate operators with radius r > 0 and parameter β > 0 by
E r , β H ( x ) = f E β H 1 D r 1 G r f E β H ( x ) ,
D r , β H ( x ) = f D β H 1 D r 1 G r f D β H ( x ) ,
where x is an M × 1 vector; f E β H ( x ) and f D β H ( x ) represent entry-wise application of the functions f E β H and f D β H , defined as f E β H ( x ) = ( x + β ) 1 and f D β H ( x ) = f E β H ( 1 x ) , respectively; the neighborhood matrix G r has entries g i j = 1 if the distance between the centroids of elements i and j is smaller than r, else g i j = 0 ; and D r = diag { G r 1 } , where 1 is the M × 1 vector with all entries set to one. The vector γ that holds the nodal values of the material indicator function γ in Ω D is
γ = E r c , β H D r c , β H D r o , β H E r o , β H ( ξ ) ,
where ξ is our design vector. In the numerical experiments, we use parameter β = 10 4 . The values for r o and r c dictate the desired minimum length scale for the fluidic and solid phase, respectively. These values are selected separately for each experiment. We remark that the existence of a minimum size scale for both phases cannot be guaranteed when using open–close filters [31]. However, in many cases the optimized designs, obtained by using such filters, possess a minimum size scale on both phases [24,28].
To impose a constraint on the total volume of the fluidic phase as well as to constrain the maximum length scale of the fluidic phase, we add additional two constraints [22,23]
v T γ V ,
v T E r m , β H ( γ ) ε ,
where V is the maximum fraction of Ω to be filled with the fluid, ε = 10 4 , parameter r m specifies the desired maximum length scale of the fluidic channels, and v is a vector with entries
v i = 1 | Ω D | Ω D φ i d Ω ,
where φ i is basis function corresponding to the ith node in Ω D and | Ω D | is the area of Ω D .

3.1. Tesla Valve

The Tesla valve is a basic device for controlling fluid motion [14,32]. Figure 3 shows the computational domain Ω for the Tesla valve. We denote the channels connected to the inlet and outlet by Ω C , and the Tesla valve occupies the domain Ω T = Ω C T Ω D , where Ω D is the design domain. The width of inlet L is 200 μ m.
The domain for the Tesla valve is a two-port network, where boundaries Γ in and Γ out represent the two ports. We consider the flow in the two directions from Γ in to Γ out and vice-versa. Henceforth, we denote the two flow directions by the forward and backward direction, respectively. In both the forward and backward direction, the flow profile is specified at the port from which the fluid flows. We let p f and u f denote the pressure and velocity, respectively, for the forward-directed fluid field. Similarly, p r and u r denote the pressure and velocity, respectively, for the backward-directed fluid field. The forward flow can be modeled by Navier–Stokes Equation (1) in Ω together with the boundary condition
u f = u 0 on Γ in , p f I + 2 η ϵ ( u f ) · n = 0 on Γ out ,
where u 0 is the given flow velocity profile at Γ in . The backward flow can be modeled by Navier–Stokes Equation (1) in Ω together with the boundary condition
u r = u 0 on Γ out , p r I + 2 η ϵ ( u r ) · n = 0 on Γ in ,
where u 0 is the given flow velocity profile at Γ out . The performance of Tesla valve is measured by diodicity, which is defined as the ratio of pressure drop of the forward flow to that of the reverse flow. From the kinetic energy viewpoint, the diodicity is also expressed as the ratio of the consumption of the forward pressure work to the reverse pressure work, under the assumption that the flux in the two opposite directions is equal. The diodicity of the Tesla valve can be alternatively calculated as [13,14,32,33]
Di = Γ in p f ( u f · n ) d Γ Γ out p f ( u f · n ) d Γ Γ in p r ( u r · n ) d Γ Γ out p r ( u r · n ) d Γ = Ω T · ( p f u f ) d Ω Ω T · ( p r u r ) d Ω .
The topology optimization problem is
min ξ R M Φ , subject to Di C max , v T γ V , 0 ξ 1 .
For the numerical simulations, each square space with an edge size equal to L is discretized by 20 × 20 quadrilateral elements. Hence 120 × 120 elements discretize the design domain Ω D . The values of α max and q in the topology method are chosen to be 5 × 10 8 and 1, respectively. Here, we set length scale parameters r c and r o small enough to ensure that the neighborhood of each node in Ω D only consists of the node itself. This yields that γ = ξ . The initial value of the design variable ξ is 1 , the Reynolds number is 100, and V is 0.55. Figure 4 shows two optimized Tesla valves. Due to physical problems and geometric constraints, there is a lower limit for the diodicity of the Tesla valve in a particular design domain. And the optimization results depend on the optimization process. Initially, we set a low C max and then increased during the optimization. The left and right valves are obtained for the value of C max is 0.351 and 0.551, respectively. We create body-fitted 2D models that only consider the fluidic phase. The diodicities computed by using these body-fitted 2D models are 0.408 and 0.600 for the left and right design in Figure 4, respectively.
Next, we add a constraint on the maximum length scale in the optimization model and consider the problem
min ξ R M Φ , subject to Di C max , v T γ V , v T E r m , β H ( γ ) ε , 0 ξ 1 .
The maximum and minimum length of the fluidic phase are r m = 200 μ m and r o = 100 μ m, respectively, and the minimum length of the solid phase is r c = 60 μ m. For this test, we set V = 0.25 and C max = 0.5 . Figure 5 shows the optimized Tesla valve with the value of C max = 0.829 . We can see that the optimized design is very similar to the two commonly used Tesla valves in series. We create a body-fitted 2D model that only considers the fluidic phase. The diodicity computed by using the body fitted 2D model is 0.882. It is clear that the design constraint of diodicity is violated because of using size control during the optimization procedure.
Finally, we create 3D models for the results in Figure 4 and Figure 5 by stretching the body-fitted 2D models. Figure 6, Figure 7 and Figure 8 shows the velocity fields of the 2D model (left image) as well as the 3D model (right image) for the optimized Tesla valve in Figure 4 and Figure 5. Having access to the 3D models, we compute the difference
δ 1 = 1 Di 3 D Di 2 D × 100 % ,
where the subscripts 2D and 3D represent the 2D model and its 3D counterpart, respectively. In the ideal case, δ 1 is equal to 0. Figure 9 shows the difference of three optimization results between 2D and 3D models with the ratio of depth to width of the inlet n. The difference in fluidic performance between the optimized 2D model and its 3D counterpart gradually decreases as the height of the 3D model increases. This confirms that the 2D model assumption that the height is infinite, the higher the stretching height of the 3D model, the smaller the difference between the 2D model and the 3D model. The difference of the optimization results using dimension control is significantly smaller than that of unused optimization results, and it is acceptable when n is greater than 1.

3.2. Microfluidic Splitters with Equivalent Outlet Flowrate

As our second test problem, we study the problem of designing microfluidic splitters originally proposed by Zhou et al. [9]. Figure 10 shows the computational domain Ω = Ω C Ω D for this test case. The width of inlet L is 100 μ m. Here, we have one inlet Γ in and three outlets Γ out ( 1 ) , Γ out ( 2 ) , and Γ out ( 3 ) ; we denote the union of the outlets by Γ out . Given a specified velocity profile, given by u 0 on the inlet, the flow is modeled by Navier–Stokes Equation (1) in Ω together with boundary condition (11).
In this example, we consider two fluidic cases in a four-port splitter. The flow direction is from Γ in to Γ out in both cases. In contrast to the previous example, we will require a specified flowrate,
Q = Γ u · n d Γ ,
at selected ports of the domain. We let p real and u real denote the pressure and velocity, respectively, for the so-called (in this work) real case. The real case is modeled by Navier–Stokes Equation (1) in Ω together with the boundary condition
Γ in u real · n d Γ = Q in on Γ in , p real I + 2 η ϵ ( u real ) · n = 0 on Γ out .
Here Q in is the given flowrate at Γ in . Similarly, we let p ref and u ref denote the pressure and velocity, respectively, for the so-called reference case. The reference case is modeled by Navier–Stokes Equation (1) in Ω together with the boundary condition
p ref I + 2 η ϵ ( u ref ) · n = 0 on Γ in , Γ out ( i ) u ref · n d Γ = Q in / 3 on Γ out ( i ) for i = 1 , 2 , 3 .
The topology optimization problem is
min ξ R M Φ real + Φ ref b + 1 v T γ + θ ( Φ real Φ ref ) subject to 0 ξ 1 .
The first term of the objective function is a fraction, in which the numerator is the sum of the viscous dissipations of the two cases and the denominator is the volume of the solid phase, to obtain an appropriate volume automatically. The second term of the objective function is the difference between the viscous dissipations of the two cases, and θ is the scalar number used to penalize this difference. By this method, the flowrate constraint of outlet in the real model is realized. The values of α max and q in the topology method are 5 × 10 9 and 1, respectively. The initial value of the design variable ξ is 0.5 . During the optimization, θ increases from 0 to 100. The left column of Figure 11 shows designs optimized without length scale control. That is, just as for the results in Figure 4, r c and r o are small enough to ensure γ = ξ . The Reynolds numbers for the simulations were 0.1 (top row), and 10 (bottom row), and the value of b in the objective function was 0.1.
The right column of Figure 11 shows designs optimized under a maximum length scale constraint for the fluid domain. That is, we solve optimization problem (20) appended with constraint (9).
The maximum length of the fluidic phase is 100 μ m, and the value of b used in the objective function is 0.1 and 0.2 for Reynolds numbers 0.1 and 10, respectively.
Having access to the 3D models, we compute the difference
δ 2 = 1 3 i = 1 3 1 3 Q ( i ) Q ( 1 ) + Q ( 2 ) + Q ( 3 ) × 100 % ,
where Q ( i ) is the flowrate at port Γ out ( i ) . In the ideal case, δ 2 is equal to 0. Figure 12 shows the difference for 3D models corresponding to optimized designs in Figure 11 as functions of the ratio of depth to width of the inlet n. The difference in fluidic performance between the optimized 2D model and its 3D counterpart gradually decreases as the height of the 3D model increases. The difference of the optimization results using dimension control is significantly smaller than that of unused optimization results when Re = 10. and it is not significantly smaller than that of unused optimization results when Re = 0.1. The value of Reynolds numbers maybe influences the effect of the proposed method.

4. Discussion and Conclusions

We can see that the difference in fluidic performance between the optimized 2D model and its 3D counterpart gradually decreases as the height of 3D model increases in Figure 9 and Figure 12. This conforms to the 2D model assumption is that the height is infinite.
There is a significant difference in fluidic performance in the case where no length scale control is applied. The difference in the optimization result of length scale control applied is much smaller. Imposing size control of the widths of the fluidic channels can effectively reduce the difference of fluidic performance between the 2D model and the 3D model obtained by stretching the 2D model, so that the 2D optimization results of flow have more practical significance. Moreover, the length scale control also ensures the manufacturability of the obtained. Numerical results demonstrate the validity of the proposed method. We hope that this paper will spur further interest and development on this topic.

Author Contributions

Y.G., H.P. and Z.L. wrote the topology optimization code; E.W. and Z.L. wrote the size control code; Y.G. and Z.L. performed the numerical examples; Y.G., E.W. and Z.L. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation of China (grant Nos. 51675506); National Science and Technology Major Project of the Ministry of Science and Technology of China (No. 2017ZX10304403); STINT (Stiftelsen för internationalisering av högre utbildning och forskning), The Swedish Foundation for International Cooperation in Research and Higher Education (No. IB2018-7470); and by the Swedish strategic research programme eSSENCE. The authors thank Professor Krister Svanberg for supplying the Matlab codes of MMA.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Whitesides, G. The lab finally comes to the chip! Lab Chip 2014, 14, 3125–3126. [Google Scholar] [CrossRef] [PubMed]
  2. Laser, D.J.; Santiag, J.G. A review of micropumps. J. Micromechan. Microeng. 2004, 14, R35–R64. [Google Scholar] [CrossRef]
  3. Oh, K.W.; Ahn, C.H. A review of microvalves. J. Micromechan. Microeng. 2006, 16, R13–R39. [Google Scholar] [CrossRef]
  4. Cai, G.; Xue, L.; Zhang, H.; Lin, J. A Review on Micromixers. Micromachines 2017, 8, 274. [Google Scholar] [CrossRef] [PubMed]
  5. Feng, X.; Ren, Y.; Hou, L.; Tao, Y.; Jiang, T.; Li, W.; Jiang, H. Tri-fluid mixing in a microchannel for nanoparticle synthesis. Lab Chip 2019, 19. [Google Scholar] [CrossRef]
  6. Borrvall, T.; Petersson, J. Topology optimization of fluids in Stokes flow. Int. J. Num. Methods Fluids 2003, 41, 77–107. [Google Scholar] [CrossRef]
  7. Gersborg–Hansen, A.; Sigmund, O.; Haber, R.B. Topology optimization of channel flow problems. Struct. Multidiscip. Optim. 2005, 30, 181–192. [Google Scholar] [CrossRef]
  8. Olesen, L.H.; Okkels, F.; Bruus, H. A high-level programming-language implementation of topology optimization applied to steady-state Navier—Stokes flow. Int. J. Num. Methods Eng. 2006, 65, 975–1001. [Google Scholar] [CrossRef] [Green Version]
  9. Zhou, T.; Liu, T.; Deng, Y.; Chen, L.; Qian, S.; Liu, Z. Design of microfluidic channel networks with specified output flow rates using the CFD-based optimization method. Microfluid. Nanofluid. 2017, 21. [Google Scholar] [CrossRef]
  10. Deng, Y.; Liu, Z.; Zhang, P.; Liu, Y.; Wu, Y. Topology optimization of unsteady incompressible Navier—Stokes flows. J. Comput. Phys. 2011, 230, 6688–6708. [Google Scholar] [CrossRef]
  11. Alexandersen, J.; WAage, N.; Andreasen, C.S.; Sigmund, O. Topology optimisation for natural convection problems. Int. J. Num. Methods Fluids 2014, 76, 699–721. [Google Scholar] [CrossRef] [Green Version]
  12. Deng, Y.; Wu, Y.; Xuan, M.; Korvink, J.G.; Liu, Z. Dynamic Optimization of Valveless Micropump. In Solid-State Sensors, Proceedings of the Actuators and Microsystems Conference, Beijing, China, 5–9 June 2011; IEEE: Piscataway, NJ, USA, 2011. [Google Scholar] [CrossRef]
  13. Deng, Y.; Liu, Z.; Zhang, P.; Wu, Y.; Korvink, J.G. Optimization of no-moving part fluidic resistance microvalves with low Reynolds number. In Proceedings of the 2010 IEEE 23rd International Conference on Micro Electro Mechanical Systems (MEMS), Hong Kong, China, 24–28 January 2010. [Google Scholar] [CrossRef]
  14. Lin, S.; Zhao, L.; Guest, J.K.; Weihs, T.P.; Liu, Z. Topology Optimization of Fixed-Geometry Fluid Diodes. J. Mech. Des. 2015, 13, 081402. [Google Scholar] [CrossRef]
  15. Guo, Y.; Xu, Y.; Deng, Y.; Liu, Z. Topology optimization of passive micromixers based on Lagrangian mapping method. Micromachines 2018, 9, 137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zeng, S.; Kanargi, B.; Lee, P.S. Experimental and numerical investigation of a mini channel forced air heat sink designed by topology optimization. Int. J. Heat Mass Transf. 2018, 121, 663–679. [Google Scholar] [CrossRef]
  17. Haertel, J.; Engelbrecht, K.; Lazarov, B.; Sigmund, O. Topology optimization of a pseudo 3D thermofluid heat sink model. Int. J. Heat Mass Transf. 2018, 121, 1073–1088. [Google Scholar] [CrossRef] [Green Version]
  18. Kobayashi, H.; Yaji, K.; Yamasaki, S.; Fujita, K. Freeform winglet design of fin-and-tube heat exchangers guided by topology optimization. Appl. Therm. Eng. 2019. [Google Scholar] [CrossRef]
  19. Yan, S.; Wang, F.; Hong, J.; Sigmund, O. Topology optimization of microchannel heat sink using a two-layer model. Int. J. Heat Mass Transf. 2019, 143, 118462. [Google Scholar] [CrossRef]
  20. Zhao, X.; Zhou, M.; Sigmund, O.; Andreasen, C.S. A “poor man’s approach” to topology optimization of cooling channels based on a Darcy flow model. Int. J. Heat Mass Transf. 2018, 116, 1108–1123. [Google Scholar] [CrossRef] [Green Version]
  21. Oh, K.W.; Lee, K.; Ahn, B.; Furlani, E.P. Design of pressure-driven microfluidic networks using electric circuit analogy. Lab Chip 2012, 12, 515–545. [Google Scholar] [CrossRef]
  22. Lazarov, B.S.; Wang, F. Maximum length scale in density based topology optimization. Comput. Methods Appl. Mech. Eng. 2017, 318, 826–844. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, F.; Lazarov, B.S. On projection methods, convergence and robust formulations in topology optimization. Struct. Multidisc. Optim. 2011, 43, 767–784. [Google Scholar] [CrossRef]
  24. Hägg, L.; Wadbro, E. Nonlinear filters in topology optimization: Existence of solutions and efficient implementation for minimum compliance problems. Struct. Multidiscip. Optim. 2017, 55, 1017–1028. [Google Scholar] [CrossRef] [Green Version]
  25. Cornish, R.J. Flow in a pipe of rectangular cross-section. Proc. R. Soc. Lond. Ser. A 1928, 120, 691–700. [Google Scholar]
  26. Svanberg, K. The method of moving asymptotes—A new method for structural optimization. Int. J. Num. Methods Eng. 1987, 24, 359–373. [Google Scholar] [CrossRef]
  27. Heijmans, H. Mathematical morphology: A modern approach in image processing based on algebra and geometry. SIAM Rev. 1995, 37, 1–36. [Google Scholar] [CrossRef] [Green Version]
  28. Sigmund, O. Morphology-based black and white filters for topology optimization. Struct. Multidiscip. Optim. 2007, 33, 401–424. [Google Scholar] [CrossRef] [Green Version]
  29. Svanberg, K.; Svärd, H. Density filters for topology optimization based on the Pythagorean means. Struct. Multidiscip. Optim. 2013, 48, 859–875. [Google Scholar] [CrossRef]
  30. Wadbro, E.; Hägg, L. On quasi-arithmetic mean based filters and their fast evaluation for large-scale topology optimization. Struct. Multidiscip. Optim. 2015, 52, 879–888. [Google Scholar] [CrossRef]
  31. Schevenels, M.; Sigmund, O. On the implementation and effectiveness of morphological close—Open and open— Close filters for topology optimization. Struct. Multidiscip. Optim. 2016, 54, 15–21. [Google Scholar] [CrossRef] [Green Version]
  32. Sato, Y.; Yaji, K.; Izui, K.; Yamada, T.; Nishiwaki, S. Topology optimization of a no-moving-part valve incorporating Pareto frontier exploration. Struct. Multidiscip. Optim. 2017. [Google Scholar] [CrossRef]
  33. Liu, Z.; Deng, Y.; Lin, S.; Xuan, M. Optimization of micro Venturi diode in steady flow at low Reynolds number. Eng. Optim. 2012, 44, 1389–1404. [Google Scholar] [CrossRef]
Figure 1. The hydraulic resistance changes with the length (red), height (blue), and width (green) of the channel, η = 1 .
Figure 1. The hydraulic resistance changes with the length (red), height (blue), and width (green) of the channel, η = 1 .
Micromachines 11 00613 g001
Figure 2. The four basic morphological operators defined by the ball B acting on the set S.
Figure 2. The four basic morphological operators defined by the ball B acting on the set S.
Micromachines 11 00613 g002
Figure 3. The computational domain Ω for the Tesla valve.
Figure 3. The computational domain Ω for the Tesla valve.
Micromachines 11 00613 g003
Figure 4. Two optimized Tesla valves with the value of C max is 0.351 (a) and 0.551 (b), respectively.
Figure 4. Two optimized Tesla valves with the value of C max is 0.351 (a) and 0.551 (b), respectively.
Micromachines 11 00613 g004
Figure 5. Tesla valve optimized with the full length scale control optimization Formulation (15).
Figure 5. Tesla valve optimized with the full length scale control optimization Formulation (15).
Micromachines 11 00613 g005
Figure 6. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve with the value of C max = 0.408 in Figure 4.
Figure 6. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve with the value of C max = 0.408 in Figure 4.
Micromachines 11 00613 g006
Figure 7. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve with the value of C max = 0.600 in Figure 4.
Figure 7. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve with the value of C max = 0.600 in Figure 4.
Micromachines 11 00613 g007
Figure 8. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve in Figure 5.
Figure 8. Reverse flow velocity field of the 2D model (a) and its 3D counterpart (b) for the optimized Tesla valve in Figure 5.
Micromachines 11 00613 g008
Figure 9. The difference δ 1 of three optimization results between 2D and 3D models with the ratio of depth to width of the inlet n.
Figure 9. The difference δ 1 of three optimization results between 2D and 3D models with the ratio of depth to width of the inlet n.
Micromachines 11 00613 g009
Figure 10. Computational domain for the second test problem. Here, the computational domain Ω comprises the design domain Ω D as well as four channels Ω C that are connected to the inlet and outlets of the domain.
Figure 10. Computational domain for the second test problem. Here, the computational domain Ω comprises the design domain Ω D as well as four channels Ω C that are connected to the inlet and outlets of the domain.
Micromachines 11 00613 g010
Figure 11. Optimized designs with aim to obtain equivalent outlet flowrate.
Figure 11. Optimized designs with aim to obtain equivalent outlet flowrate.
Micromachines 11 00613 g011
Figure 12. The difference δ 2 of four optimization results between 2D and 3D models with the ratio of depth to width of the inlet n.
Figure 12. The difference δ 2 of four optimization results between 2D and 3D models with the ratio of depth to width of the inlet n.
Micromachines 11 00613 g012

Share and Cite

MDPI and ACS Style

Guo, Y.; Pan, H.; Wadbro, E.; Liu, Z. Design Applicable 3D Microfluidic Functional Units Using 2D Topology Optimization with Length Scale Constraints. Micromachines 2020, 11, 613. https://doi.org/10.3390/mi11060613

AMA Style

Guo Y, Pan H, Wadbro E, Liu Z. Design Applicable 3D Microfluidic Functional Units Using 2D Topology Optimization with Length Scale Constraints. Micromachines. 2020; 11(6):613. https://doi.org/10.3390/mi11060613

Chicago/Turabian Style

Guo, Yuchen, Hui Pan, Eddie Wadbro, and Zhenyu Liu. 2020. "Design Applicable 3D Microfluidic Functional Units Using 2D Topology Optimization with Length Scale Constraints" Micromachines 11, no. 6: 613. https://doi.org/10.3390/mi11060613

APA Style

Guo, Y., Pan, H., Wadbro, E., & Liu, Z. (2020). Design Applicable 3D Microfluidic Functional Units Using 2D Topology Optimization with Length Scale Constraints. Micromachines, 11(6), 613. https://doi.org/10.3390/mi11060613

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