非均匀磁场下Maxwell磁纳米流体的拉伸流动与磁扩散分析.pdf

应用数学和力学编委会，：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析吴学珂，刘春燕，白 羽，张 艳，（北京建筑大学 理学院，北京；建筑结构与环境修复功能材料北京市重点实验室，北京）摘要：磁性纳米颗粒可以提升聚合物的导电性和导热性等性能，被广泛应用于机械、生物医学、能源存储等领域 当外界施加非均匀磁场时，感应磁场在高 数的情况下不可忽略 为探究磁性纳米颗粒对层流边界层内黏弹性流体非稳态拉伸流动与磁扩散的影响，将时间分布阶 本构方程与动量方程耦合，建立了二维不可压缩 磁纳米流体的速度与磁扩散偏微分方程组 采用有限差分法进行数值分析，通过控制磁性纳米颗粒种类、体积分数和磁参数大小，分析了流体的速度和感应磁场在边界层内的分布 研究发现：在熔融聚合物中添加 纳米颗粒后，流体的速度、感应磁场最大，速度和磁边界层的厚度最厚；纳米流体的松弛时间参数增大，速度与磁扩散均减小；另外，随着磁参数增大，流体的速度边界层厚度减小，磁边界层厚度增大；纳米颗粒的体积分数越大，流体流动越快，感应磁场越小 因此，非均匀磁场下在聚合物中添加磁性纳米颗粒的研究，为改善材料的性能给予了可参考的数据 关 键 词：流体；磁性纳米颗粒；感应磁场；数值差分格式中图分类号：文献标志码：，（，；，）：，应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金青年科学基金（）；北京市教育委员会科技计划一般项目（）作者简介：吴学珂（），女，硕士生（：）；刘春燕（），女，博士，硕士生导师（通讯作者：）；白羽（），女，教授，博士，硕士生导师（：）；张艳（），女，教授，博士，硕士生导师（：）引用格式：吴学珂，刘春燕，白羽，张艳 非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析 应用数学和力学，（）：，：；引 言聚合物是经聚合反应生成的高分子化合物，而熔融的聚合物在力的作用下可以流动，是一类同时具有黏性和弹性的黏弹性流体 在聚合物中加入纳米颗粒可有效改善聚合物的性能，提高其韧性和传热速率，在工业和生物医学等领域被广泛应用，增加了能源的利用率 磁性纳米颗粒具有磁性特征且具有纳米颗粒的独特效应，常见的磁性纳米颗粒有金属铁、钴、镍、金属氧化物 和 等，在熔融聚合物中加入微量的磁性纳米颗粒，可以使材料的各种性能如导电性、导热性、阻隔性等提高 等应用边界层理论研究了包含电磁流体动力学在拉伸板上复合纳米流体的非稳态常驻点流动 等研究了非均匀磁场对 水基纳米流体强制对流换热的影响 对磁性纳米流体施加外磁场，可以在非接触的条件下对其流动产生影响 这种非接触操控方式在实际操作中更便于实现，且操控范围广，不受 值、离子强度、表面电荷和温度的影响 在多数实际情况下，外加磁场是可变的，近年来，很多研究探讨了非均匀磁场对纳米流体流动的影响 等研究了变磁场下磁场力对纳米流体传热的影响 等研究了非均匀磁场对磁性纳米流体在开腔通道内混合对流换热的影响，发现腔加热壁的上角附近产生了涡流 磁流体力学（）是结合经典流体力学和电动力学的方法研究导电流体和磁场相互作用的学科，在天体物理、地球物理、宇航工程、电磁学以及工程技术中都有广泛的应用 等研究了导电金属流体在感应磁场作用下的流动，发现增加磁参数会使速度提高、感应磁场降低 等采用有限差分格式数值计算了随时间变化的黏性纳米流体在感应磁场中的流动，并讨论了 运动和热泳运动 但在现有的关于非均匀磁场的数值模拟研究中，尚未清楚阐明感应磁场的物理机制 考虑到黏弹性磁纳米流体在流动过程中的复杂特性，将分数阶导数引入本构关系能更灵活地描述黏弹性流体的性质 杨旭等基于分数阶微积分理论，采用空间分数阶导数建立了圆管内分数阶非 流体本构模型，为非 流体的记忆特征提供了一种建模方法 等通过在本构关系中引入分数阶 剪切应力和 热流模型，研究了 流体在平面上的非稳态 对流换热 然而上述模型中分数阶导数参数固定，导致了有限的记忆特性和非局部特性，且无法准确描述一些复杂的动力学过程，例如复合材料的流变特性 分布阶本构模型是分数阶导数在参数值范围内的积分，具有不同的时间和空间特征，作为更有效的工具引起了很多学者的关注 等建立了空间分布阶本构关系来研究边界层中的流动和传热 等基于 流体流动和 传热的传统本构关系，建立了分布阶导数的第 期 吴学珂，等：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析非稳态 对流边界层流动和传热模型 等将分布阶导数引入 流体的本构模型，并分析了相关参数的影响 基于上述研究发现，在感应磁场作用下，非稳态时间分布阶的 磁纳米流体流动的问题很少有人研究 本文将时间分布阶 本构关系代入动量方程，并与磁扩散方程建立流动和感应磁场的控制方程，然后结合有限差分方法与 算法获得控制方程的数值解，最后分析相关参数对流动和感应磁场的影响 数 学 模 型时间分布阶 流体的本构方程为（），（）其中，是随时间变化的剪切应力，和 分别指的是松弛时间参数和黏度，是分数阶导数参数，（）是控制分数阶导数的权重系数且满足条件（），（）是 定义的分数阶导数，如下所示：（）（）（）（），其中，（）为 函数 考虑线性拉伸板上的二维不可压缩非稳态 磁纳米流体的边界层流动问题，如图 所示 建立二维直角坐标系，其中 轴与平板平行，轴垂直于平板 施加非均匀磁场 （），假设（，）和（，）分别是沿着板和垂直于板的速度和感应磁场的分量，并假设感应磁场的法向分量 在壁处消失、平行分量 在边界层边缘接近给定值 则时间分布阶 磁纳米流体的边界层流动和磁扩散控制方程为 ，（）（）（）（）（）（），（）（）满足下列初始条件和边界条件：，：，（）式（）（）中，和分别是纳米流体的黏度和密度，是磁扩散系数，和分别是磁导率和纳米流体的电阻率，为拉伸参数，为特征长度（板长）应 用 数 学 和 力 学 年 第 卷图 物理模型示意图 纳米流体的物理性质参数为（），（），（）其中 是磁性纳米颗粒的体积分数，是流体的黏度，和分别是流体和磁性纳米颗粒的密度，和分别是流体和磁性纳米颗粒的电阻率 磁性纳米颗粒铁、钴、和 主要的物理性质如表 所示 表 磁性纳米颗粒的物理性质 （）（）对方程（）（）进行无量纲化：，（）其中 是磁参数，是磁 数，是 数，得到无量纲控制方程如下（为了方便，后面将省略标记“”）：，（）（）（）（）（）（）第 期 吴学珂，等：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析 （）（），（）（）无量纲化后的初值和边界条件为 ：，：，（）因此，建立了二维非稳态 磁纳米流体的控制方程（）（），方程（）为相应的初始条件和边界条件 数值差分格式采用有限差分法与 算法相结合来求解耦合的二维分布阶控制方程（）（）对时间和空间进行网格划分，设，分别为沿 轴和 轴的空间步长，为时间步长 定义 ，其中，是网格划分的数量 定理 设 ，（），则 分数阶导数为（）（）（）（）（）（）（）（）（）（）（）（），称之为 算法，其基本思想是将被积函数中出现的 的分数阶导数直接用数值微分公式逼近，还可改写为（）（）（）（）（）（）（）（）（）（）（），（）其中 （），根据中点 求积规则，使用多个分数项的加权和变换时间分布阶分数导数，可以获得分布项的数值离散格式：（）（），（）（），（）（），（）应 用 数 学 和 力 学 年 第 卷其中 表示区间，中分数阶 的步长，且 （）（，），其余项的离散格式，和方程（）类似 设在网格点（，）处速度的精确解和数值解分别表示为（，）和，则对于控制方程中整数阶导数的离散格式分别采用向后差分和中心差分 对于离散区域，中 阶的时间分数阶导数，利用定理 中的 数值离散格式（），可得到，（），（）（，）（，）（），（），（）（，）（）（，），（，）（），（）其中，（），将离散格式（）、（）和（）代入控制方程（）（）和初边值条件（），得到求解问题的数值差分格式 结 果 分 析利用第 节中的有限差分方法得到动量方程和磁扩散方程的数值解，并说明相关参数对速度和感应磁场的影响 构造解析解验证了差分格式的收敛性，如图 所示 图 数值解和解析解的比较 图 和图 描绘了磁性纳米颗粒、和 对速度和感应磁场的影响 从图 可以看出，密度相近的 和 纳米颗粒对速度的影响几乎相同，而密度最大的 使流速受到最大限制，速度边界层厚度最薄 从物理角度看，当磁性纳米颗粒的体积分数相同时，密度越大，颗粒间的间隙越小，从而阻碍了第 期 吴学珂，等：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析流体的运动 图 表明，添加的 纳米颗粒电阻率最小，则感应磁场最大，磁边界层厚度最厚；对于电阻率相近的 和，感应磁场的分布几乎相同 由于磁性纳米颗粒的电阻率与电导率呈反比，所以当电阻率增大时，感应磁场减小 图 和图 为变化的磁参数 对 纳米流体速度和感应磁场的影响 图 中的曲线表明流动的模式一致，流速值随着磁参数的增加而降低 在非均匀磁场作用下，边界层内产生与流体运动方向相反的磁场力，所以当 时，磁场力对流体的速度起抑制作用，使得边界层厚度变薄 当 时，无磁场力，纳米流体速度由于内摩擦阻力作用逐渐减小 如图 所示，在边界层内感应磁场随 的增大而增大 因为 增大时磁导率增加，磁扩散效应显著，从而感应磁场增加 图 不同磁性纳米颗粒对速度的影响图 不同磁性纳米颗粒对感应磁场的影响 图 不同 对速度的影响图 不同 对感应磁场的影响 图 和图 显示了速度和感应磁场随磁性纳米颗粒 体积分数 的变化 当 增加时，速度增大，边界层厚度变厚，感应磁场减小 由于拉伸板的运动导致固体颗粒运动增加，从而速度随着磁纳米颗粒的加入而增加 而在外加磁场的作用下，随着 的增加，速度的增大导致磁对流项增大，使得磁扩散效应减小，从而在边界层内感应磁场减小 图 和图 显示了松弛时间参数 对 纳米流体中速度和感应磁场的影响 当选择权重系数（）（），无量纲时间 时，可以观察到 越大，速度和感应磁场越小，这是由于流体的黏性增加，阻碍了流体的运动进而导致磁场减小，这和应力松弛现象是一致的 应 用 数 学 和 力 学 年 第 卷图 不同 对速度的影响图 不同 对感应磁场的影响 图 不同 对速度的影响图 不同 对感应磁场的影响 结 论磁性纳米颗粒可以改善聚合物的性能，线性拉伸板上非稳态流动与磁扩散的研究为其提供了良好的理论基础 本文利用分布阶本构模型，研究了黏弹性基磁纳米流体在非均匀磁场下的二维不可压缩边界层流动，并考虑了高 数下的感应磁场分布 将时间分布阶 方程引入动量方程，建立了速度和磁扩散控制方程组 使用有限差分法和 算法耦合求解偏微分方程组的数值解，并验证了收敛性，结合图像分析了相关物理参数对速度和感应磁场的影响 所得主要结论如下：流体的速度和感应磁场在添加 纳米颗粒时达到最大，此时速度边界层和磁边界层的厚度最厚；由于磁场力的影响，增大的磁参数阻碍了流体的运动，增强了磁扩散；随着 纳米流体体积分数增加，速度增大，感应磁场减小；此外，较高的松弛时间参数会使 流体的流动和磁场的扩散减弱 因此，本文从数值解的角度研究了添加磁性纳米颗粒的熔融聚合物在非均匀磁场下的边界层拉伸流动与磁扩散，对制备出高性能的聚合物材料具有重要意义 致谢 本文作者衷心感谢北京建筑大学金字塔人才培养工程项目（）、北京建筑大学促进高校内涵发展定额项目（）和北京建筑大学研究生创新项目（）对本文的资助 第 期 吴学珂，等：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析参考文献（）：庄昕，刘付军，孙艳萍，等 非等温黏弹性聚合物流体圆柱绕流的高精度数值模拟 应用数学和力学，（）：（，（）：（），：，（）：，：，：，：，：赵芳彪 非均匀磁场下磁液液滴生成与输运的实验研究 昆明：昆明理工大学，（：，（），：，：翟梦情，李琦，郑素佩 求解一维理想磁流体方程的移动网格熵稳定格式 计算力学学报，（）：（，（）：（），（）：，：，（）：杨旭，梁英杰，孙洪广，等 空间分数阶非 流体本构及圆管流动规律研究 应用数学和力学，（）：（，（）：（），：，（）：，：应 用 数 学 和 力 学 年 第 卷 ，：，：，：，：，（）：，（）：，：郭柏灵，蒲学科，黄凤辉 分数阶偏微分方程及其数值解 北京：科学出版社，（，：，（）第 期 吴学珂，等：非均匀磁场下 磁纳米流体的拉伸流动与磁扩散分析