《Scientific Reports》：Comparative numerical and experimental evaluation of 3D-printed micromixers for biotechnological applications
编辑推荐：
摘要
微混合器在基于微流控的生物技术和诊断中至关重要，在这些应用中，可重复的反应和快速的质量传递是关键。尽管存在许多几何构型，但在尺寸、制造方法和评估技术之间进行一致的比较仍然具有挑战性。在本研究中，研究人员对六种微混合器设计进行了比较分析，重点关注混合效率、
摘要
微混合器在基于微流控的生物技术和诊断中至关重要，在这些应用中，可重复的反应和快速的质量传递是关键。尽管存在许多几何构型，但在尺寸、制造方法和评估技术之间进行一致的比较仍然具有挑战性。在本研究中，研究人员对六种微混合器设计进行了比较分析，重点关注混合效率、压降、剪切应力和能量耗散速率。这些器件采用高清晰度3D打印制造，并通过实验和高分辨率计算流体力学（Computational Fluid Dynamics，CFD）模拟进行研究。测试覆盖了较宽的雷诺数范围（Re = 0.1–200）。在该范围的两端，大多数设计都实现了充分的混合。最大的差异出现在中间区域（Re = 1–50），其中Kenics几何构型的性能优于所有其他设计。网格优化的CFD结果与实验趋势相符，并量化了局部剪切和能量耗散速率。使用工业相关的中国仓鼠卵巢（Chinese hamster ovary，CHO）细胞系进行的活/死测定显示，混合后立即及后续培养期间的细胞存活率未发生改变。通过将经过实验验证的CFD预测、系统性的网格优化评估和生物学验证整合在单一框架内，本研究为比较微混合器设计以及选择符合混合性能、压降、体积、剪切应力和能量耗散速率等预期组合的设计建立了一致的基础。
### 3D打印微混合器比较评价研究解读
#### 一、研究背景与问题
微流控技术在临床诊断、环境监测、细胞培养和药物开发等生物技术领域正获得日益广泛的应用。其紧凑的外形、低试剂消耗以及与自动化工作流程的兼容性，使其在转化研究和临床实施中都具有显著吸引力。在连续流微流控系统中，高效的流体混合是一项关键的功能需求，因为试剂输送、颗粒形成等操作都依赖于层流条件下的快速且可重复的均质化。然而，在微通道中，流动通常由低雷诺数表征，粘性力占主导，湍流不会发展，质量传递主要受分子扩散控制，导致特征混合长度和时间过长。
为克服这些限制，研究者开发了多种特定设计的通道几何构型，通过引入额外的平流输运机制来促进二次流和混沌平流。与依赖外部能量输入的主动微混合器相比，被动微混合器仅通过几何设计和流动诱导效应实现混合，其制造更简单、更易集成、在连续流操作中更稳健。过去二十年中，已发展出多种基于重复单元的被动微混合器几何构型，包括Tesla型结构、Kenics螺旋元件、分割-重组（Split-and-Recombine，SAR）单元、交错人字形混合器、Dean涡旋及螺旋形几何构型、基于障碍物的设计以及扩张-收缩元件等。
尽管许多微混合器几何构型展现出良好的混合性能，但不同研究的结果难以直接比较，因为各项研究在通道尺寸、制造方法、雷诺数范围和评估技术等方面存在显著差异。此外，压降数据未得到一致报告，剪切应力或局部能量耗散速率的定量评估也较少纳入研究，进一步限制了跨研究的可比性。在方法学层面，许多计算流体力学（Computational Fluid Dynamics，CFD）研究依赖粗糙网格或不够严格的网格收敛标准，可能引入数值扩散导致混合效率被高估；不同研究中混合指数的数学定义、空间评估窗口或强度归一化策略也各不相同，即使名义研究条件相似也会产生不一致的结果。这些方法学差异使得比较研究、验证模拟或建立可靠的微混合器性能基准变得十分困难。
高分辨率3D打印技术的日益普及使得制造具有可重复通道尺寸的复杂微混合器几何构型成为可能，这为在相同制造约束下并排比较根本不同的设计提供了条件。然而，在单一实验和计算框架内对不同几何构型进行系统比较的研究仍然有限。此外，在生物技术应用中，微混合器越来越多地用于梯度生成、生物传感、连续低温保护剂-细胞混合或纳米颗粒形成等场景，混合效率必须与机械应力暴露相平衡，因为哺乳动物细胞——特别是悬浮适应的中国仓鼠卵巢（Chinese hamster ovary，CHO）细胞——对流体动力应力敏感。尽管如此，局部剪切应力、能量耗散速率以及后续细胞反应的定量评估尚未常规纳入微混合器研究。
为弥补上述空白，研究人员对六种3D打印被动微混合器几何构型在宽雷诺数范围（0.1–200）内进行了系统比较，涵盖微流控应用中常见的流动区域，通过实验和CFD模拟评估混合效率、压降、剪切应力和局部能量耗散速率，并将CHO细胞活力作为生物验证指标，建立了统一的微混合器性能基准框架。
#### 二、研究结论与意义
本研究的核心结论表明，不同微混合器几何构型之间的性能差异显著。在扩散主导的极低雷诺数条件下，所有设计的混合行为几乎相同；最大差异出现在中间区域（Re = 1–50），即扩散已不足而几何诱导平流尚未完全建立时，仅Kenics和HC混合器在该范围内保持高混合指数，Tesla混合器则表现出明显的最低值且仅在较高雷诺数下恢复。DVBM和T型混合器在Re > 1的任何测试条件下均未达到完全混合。流体动力载荷与混合性能不呈比例关系，本研究中观察到的最大能量耗散率约为8 × 10? W m?3，比已报道的悬浮哺乳动物细胞损伤水平低一至两个数量级，单次通过任何测试几何构型均未在72小时后续培养中显著降低CHO细胞活力或活细胞密度。网格收敛研究表明收敛行为不均匀，中间区域即使在最精细分辨率下仍具网格依赖性。该研究不支持单一的微混合器优选指标，但为根据应用需求权衡混合效率、压降、流体动力载荷、占用面积和可制造性等因素提供了定量基础。该论文发表在《Scientific Reports》。
#### 三、主要关键技术方法
本研究采用高清晰度多喷头3D打印（Projet MJP 2500 Plus，标称层高32 μm，面内分辨率800 × 900 dpi）制造六种微混合器。实验混合性能通过溴酚蓝pH指示剂显色法结合自定义混合指数（Mixing Index，MI）量化。CFD模拟使用COMSOL Multiphysics? v. 6.3在稳态层流条件下进行，系统评估了五种网格密度下的网格收敛性，并采用网格收敛指数量化数值不确定性。剪切应力和能量耗散率（Energy Dissipation Rate，EDR）由速度梯度张量的对称部分计算。细胞实验使用来自德国莱布尼茨汉诺威大学技术化学研究所的CHO-K1悬浮细胞系，以台盼蓝排斥法评估细胞活力和活细胞密度。
#### 四、研究结果
**实验混合性能**
研究人员通过溴酚蓝显色实验在Reynolds数为0.1至100范围内量化了六种设计的混合指数。在Re = 0.1时，混合主要由扩散驱动，所有几何构型在2–4 μL内即实现均质化。在Re = 1–10范围内，大多数设计的混合性能显著下降，仅Kenics和HC混合器实现完全混合（MI > 0.9），Tesla混合器的降幅最明显，DVBM在Re = 1时即表现出极低混合指数。在Re ≈ 10–50的平流主导区域，Kenics混合器持续实现快速均质化。在Re ≥ 50时，Kenics、Tesla、HC和CSSR达到高混合水平，而DVBM和T型混合器未能实现，DVBM的Dean涡旋在该通道尺寸和流速下未能充分发挥。
**CFD模拟与网格收敛研究**
CFD模拟采用稳态层流假设，在CSSR混合器上于四个Reynolds数（0.1、1、10、100）下进行五种网格密度的收敛研究。结果表明在Re = 0.1和100时各网格水平给出几乎一致的混合指数曲线；在Re = 1时呈明显收敛趋势；而在Re = 10的中间区域，预测结果对网格分辨率表现出显著依赖性。网格收敛指数量化显示不确定性在Re = 10处达到峰值（中位数±0.18，最大±0.26），而其他Reynolds数下可忽略不计。研究还验证了交叉扩散稳定的影响可忽略，最终选用最大单元尺寸30 μm的统一网格并关闭交叉扩散。
**模拟与实验结果对比**
所有几何构型的CFD预测均与实验混合进展密切吻合，Lin一致性相关系数范围为0.920–0.982。Kenics、T型混合器和HC的偏差不超过0.07，而CSSR、DVBM和Tesla最大偏差达0.16，对应的分别是各几何构型的过渡区域，该区域扩散和平流输运贡献相当、数值敏感性最高。实验测量在弱混合状态下也表现出较大变异性，模拟假设光滑壁面而未考虑打印侧壁粗糙度，这些因素共同导致偏差。
**流体动力学特征**
在Re = 100条件下，CSSR、Kenics和Tesla混合器表现出最高的整体耗散水平，平均EDR超过数千W m?3，90百分位延伸至10? W m?3范围，源于狭窄收缩结构导致的强烈二次流。HC和T型混合器处于中间水平，T型混合器虽混合性能差但耗散相对较高，因其为达到相同内部体积而采用长而窄的通道产生陡峭壁面法向速度梯度。DVBM的平均和百分位EDR值在所有设计中最低。本研究观察到的最大局部EDR约8 × 10? W m?3，远低于CHO细胞致死性流体动力响应的报道水平（10?–10? W m?3量级）。
**细胞对混合的响应**
CHO悬浮细胞以Re = 1、10、100、200通过各混合器单次后，即刻及后续72小时培养中，各条件下的活细胞密度和活力均与未处理对照相当。暴露后立即检测到1–2%的轻微活力下降，但细胞在72小时内完全恢复，没有任何混合器或Reynolds数条件导致持续活力丧失或生长受损。这表明尽管存在几何收缩、局部剪切应力升高和测量的侧壁粗糙度，单次通过测试混合器在所考察操作条件下不会可测量地影响CHO细胞活力或后续生长。
#### 五、总结讨论与结论
讨论部分指出，微混合器的性能强烈依赖于几何构型与流动区域之间的相互作用。Kenics混合器凭借重复的流体界面分割、旋转和重新定向机制，在过渡条件下表现最稳健。Tesla混合器在中间Reynolds数区域的性能下降源于其在该范围内缺乏充分发展的回流涡旋结构。DVBM在纳米颗粒生产中虽广泛应用，但通常在显著更高的流速下运行，其简单几何构型还便于注塑成型。细胞实验结果与EDR定量分析一致，确认所研究的操作条件远低于机械损伤阈值。数值方面，在单一有利工作点进行的收敛性评估可能高估微混合器模拟的可靠性，中间区域需要专门的数值研究。
研究结论部分翻译如下：实验和计算相结合的分析突出了所研究微混合器几何构型之间的明显性能差异。在扩散主导的极低Reynolds数下，所有设计表现相似。最大差异出现在中间区域（Re = 1–50），此时扩散不再充分而几何诱导平流尚未完全建立，仅Kenics和HC混合器在该范围内保持高混合指数，Tesla混合器表现出明显最低值且仅在更高Reynolds数下恢复，DVBM和T型混合器在Re > 1的所有测试条件下均未达到完全混合。流体动力载荷与混合性能不呈比例关系，本研究观察到的最大EDR约为8 × 10? W m?3，比已报道的悬浮哺乳动物细胞损伤水平低一至两个数量级。与此一致，单次通过任何测试几何构型在72小时后续培养中未显著降低CHO细胞活力或活细胞密度。由于该评估仅限于单一CHO悬浮细胞系，其他细胞类型对流体动力应力的敏感性需另行评估。网格收敛研究显示非均匀收敛行为：极低和极高Reynolds数下的预测对网格密度基本不敏感，而中间区域即使在最精细分辨率下仍具网格依赖性，因此该范围内的定量结论基于实验数据而非模拟。这些结果不支持单一的微混合器选择指标，因为混合效率、压降、流体动力载荷、占用面积和可制造性在不同应用中的权重不同且可能相互制约。在需要宽流速范围高效混合的场景中，Kenics几何构型是最稳健的设计；在背压受限的集成化芯片实验室设备中，HC几何构型以较低水力需求提供相当的混合性能；在需要注塑量产时，平面Tesla和DVBM几何构型更为适合，两者均需高流速才能有效混合且DVBM所需流速显著更高。本研究提出的框架为根据具体应用需求权衡这些取舍提供了定量基础。
