• 中文核心期刊
  • 中国科技核心期刊
  • ISSN 1007-6336
  • CN 21-1168/X

MessageBoard

Dear readers, authors and reviewers, if you have any questions about the contribution, review, editing and publication of this magazine, you can add comments on this page. We'll get back to you as soon as possible. Thank you for your support!

UserName
Email
Phone
Title
Content
Code
Volume 40 Issue 1
Jan.  2021
Article Contents

Citation:

Numerical simulation of convective diffusion of pollutants on the wall surface of ocean pillar structure

  • Received Date: 2019-06-13
    Accepted Date: 2019-12-12
  • For the problem of flow around the cylinder, many domestic and foreign scholars have done a lot of research, but there is limited research on the pollution caused by the surface coating of the column to the marine environment. In this paper, the concentration field caused by the wall surface contaminant flow under the action of two-dimensional cylinder is numerically simulated. Under the oscillating eddy current field generated by the flow around a cylinder with a Reynolds number of 75~175, the cylindrical surface pollutants are used as the wall condition, the concentration field under various oscillating eddy current fields is analyzed, and the correlation and distribution law between the concentration field and the eddy current field are analyzed. The numerical simulation results show that under the action of the oscillating eddy current field, the concentration distribution of the pollutants is synchronized with the distribution of the vortex, the concentration of the center of the vortex is relatively high, and the concentration at the edge of the vortex is relatively low. As the Reynolds number increases, the concentration of the vortex position decreases.
  • 加载中
  • [1] XU J S, LI G X, HORRILLO J J, et al. Calculation of maximum allowable free span length and safety assessment of the DF1-1 submarine pipeline[J]. Journal of Ocean University of China, 2010, 9(1): 1-10. doi: 10.1007/s11802-010-0001-4
    [2] LIANG D F, CHENG L, LI F J. Numerical modeling of flow and scour below a pipeline in currents: part II. Scour simulation[J]. Coastal Engineering, 2005, 52(1): 43-62. doi: 10.1016/j.coastaleng.2004.09.001
    [3] 陈薛浩, 朱志夏. 基于Fluent的海底管线附近流场分析[J]. 上海交通大学学报, 2012, 46(3): 458-462.
    [4] SONG J N, LU L, TENG B, et al. Laboratory tests of vortex-induced vibrations of a long flexible riser pipe subjected to uniform flow[J]. Ocean Engineering, 2011, 38(11/12): 1308-1322.
    [5] 缪国平, 朱仁传, 程建生, 等. 海上风电场建设与海洋工程装备研发中若干水动力学关键技术问题[J]. 上海造船, 2009, (1): 19-22, 25.
    [6] 王 松. 低碳经济环境下的海上风力发电发展探析[J]. 现代经济信息, 2017, (21): 282-283. doi: 10.3969/j.issn.1001-828X.2017.21.225
    [7] 李 军, 王洪彬, 李 燕. 影响海底管道寿命的主要因素及防范建议[J]. 石油工程建设, 2007, 33(2): 35-38. doi: 10.3969/j.issn.1001-2206.2007.02.010
    [8] 黄宏杰. 低雷诺数疏水圆柱绕流减阻机制研究[D]. 合肥: 中国科学技术大学, 2018.
    [9] 王旭东, 缴立立. 海底管线节点防腐涂层施工[J]. 涂料工业, 2015, 45(6): 80-83.
    [10] LIU C, ZHENG X, SUNG C H. Preconditioned multigrid methods for unsteady incompressible flows[J]. Journal of Computational Physics, 1998, 139(1): 35-57. doi: 10.1006/jcph.1997.5859
    [11] RAJANI B N, KANDASAMY A, MAJUMDAR S. Numerical simulation of laminar flow past a circular cylinder[J]. Applied Mathematical Modelling, 2009, 33(3): 1228-1247. doi: 10.1016/j.apm.2008.01.017
    [12] 黄钰期, 邓 见, 任安禄. 黏性非定常圆柱绕流的升阻力研究[J]. 浙江大学学报: 工学版, 2003, 37(5): 596-601.
    [13] 郝 鹏, 李国栋, 杨 兰, 等. 圆柱绕流流场结构的大涡模拟研究[J]. 应用力学学报, 2012, 29(4): 437-443. doi: 10.11776/cjam.29.04.A013
    [14] 桂洪斌, 卢腾超. 海洋立管涡激振动的数值模拟研究[J]. 船舶力学, 2016, 20(1/2): 93-98.
    [15] 胡 琴, 万长成, 肖 祥, 等. 高雷诺数下海底悬跨管道绕流数值模拟[J]. 机械强度, 2018, 40(1): 154-159.
    [16] 顾恩凯. 考虑流固耦合的大位移柔性立管涡激振动研究[D]. 青岛: 中国海洋大学, 2012.
    [17] EL AKOURY R, BRAZA M, PERRIN R, et al. The three-dimensional transition in the flow around a rotating cylinder[J]. Journal of Fluid Mechanics, 2008, 607: 1-11. doi: 10.1017/S0022112008001390
    [18] 魏志理, 孙德军, 尹协远. 圆柱尾迹流场中横向振荡翼型绕流的数值模拟[J]. 水动力学研究与进展, 2006, 21A(3): 299-308.
    [19] 方 正, 李 炜. 粘性流体圆柱绕流的混合有限分析解[J]. 水动力学研究与进展, 1998, 13A(2): 133-139.
    [20] 孟元元. 圆柱绕流的数值模拟研究[D]. 兰州: 甘肃农业大学, 2010.
  • 加载中

Figures(14) / Tables(4)

Article Metrics

Article views(337) PDF downloads(4) Cited by()

Related
Proportional views
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Numerical simulation of convective diffusion of pollutants on the wall surface of ocean pillar structure

  • College of Environment Science and Engineering, Dalian Maritime University, Dalian 116026, China

Abstract: For the problem of flow around the cylinder, many domestic and foreign scholars have done a lot of research, but there is limited research on the pollution caused by the surface coating of the column to the marine environment. In this paper, the concentration field caused by the wall surface contaminant flow under the action of two-dimensional cylinder is numerically simulated. Under the oscillating eddy current field generated by the flow around a cylinder with a Reynolds number of 75~175, the cylindrical surface pollutants are used as the wall condition, the concentration field under various oscillating eddy current fields is analyzed, and the correlation and distribution law between the concentration field and the eddy current field are analyzed. The numerical simulation results show that under the action of the oscillating eddy current field, the concentration distribution of the pollutants is synchronized with the distribution of the vortex, the concentration of the center of the vortex is relatively high, and the concentration at the edge of the vortex is relatively low. As the Reynolds number increases, the concentration of the vortex position decreases.

    HTML

  • 随着海洋工程的发展,海底管道技术日趋成熟,广泛应用于海洋工程的石油和天然气工业[1],如海底运输管线[2-3]、海洋立管[4]等。绿色能源理念的兴起,使人们的目光开始投向清洁的可再生能源——风能。海上风电场是近年来海洋工程的一个重要开发领域[5]。海上风电场要面对风和波浪的双重负荷考验,对机组支撑结构要求很高[6]。风电机组支撑结构通常选用圆柱体作为塔基,在一定的流速条件下,水流流经柱体时,柱体背流面两侧会产生交替脱离柱体表面的漩涡[7],这种流场被称为振荡涡流场。

    国内外大量学者已对海水中柱体周边的局部冲刷问题和涡激振动对柱体磨损冲击问题进行了详尽的研究[8],在实际情况中,还要考虑腐蚀问题。海洋结构物通常会在表面喷涂防腐材料作为保护,延缓腐蚀[9]。对于海洋柱体防腐有两种最常用的方法,一种是由常用的海洋涂料组成的外涂层,例如,3PE材料和3PP材料等,可以避免海水与柱体的直接接触;另一种是金属涂层,金属微粒表面形成的致密氧化膜有防腐蚀的作用。由于圆柱绕流产生的振荡涡流场效应,柱体的保护涂料在一定条件下会从表面脱落入海,对海洋环境造成污染,导致污染物质在海洋生物体内累积,进而影响人类健康。

    本文在圆柱绕流产生的涡流场条件下,对柱体结构壁面污染物在海水中的对流扩散过程进行数值模拟,得出其浓度场,分析其扩散及分布规律,为下一步的防护与治理提供科学支持。

  • 1.   材料与方法

      1.1.   流体流控制方程

    • 本文采用黏性不可压缩流体Navier-Stokes方程,求解低雷诺数下圆柱绕流问题的流场,方程如下。

      连续性方程:

      动量方程:

      式中:${{u}_{i}}(i=1{\text{、}}2) $xy方向的流体速度;$\rho $为流体密度;p为流体压力;v为流体粘性系数。

    • 1.2.   浓度方程

    • 假设在浓度方程中,扩散系数各向同性,方程如下。

      式中: C为污染物浓度;D为扩散系数。

    • 1.3.   计算参数

      1.3.1.   计算区域
    • 计算域为一个二维矩形域,圆柱体位于矩形域的左侧,计算域尺寸为10.0 d×25.0 dd为圆柱直径,本文设d=0.1 m),网格总数为20245个,为多边形结构化网格。圆柱体结构的圆心距上、下边界均为5 d,距左侧入流边界5 d,计算域如图1所示。

      Figure 1.  Calculation area

    • 1.3.2.   初、边值条件
    • 边界条件:在计算域左侧边界给定来流速度u1=0.15~0.35 m/s、u2=0;上、下边界设定为对称面边界;右侧设为压力出口边界,并设为非反射条件;圆柱表面设置为无滑移壁面。

      初值:除边界外,内部节点速度均为0。

      在方程(1)(2)基础上,模拟浓度场,设置圆柱体边界浓度条件为1,左侧及上下边界为0。

      模型及参数设置:根据本文研究内容,通过控制来流速度设置雷诺数。

      式中:uρμ分别为流体的流速、密度与黏性系数(选择水作为模拟对象);d为特征长度,此处为圆柱体直径。来流速度区间为0.15~0.35 m/s,故而雷诺数范围为75~175。

    • 1.3.3.   数值方法
    • 本文数值求解方法为有限体积法,将求解域细分为有限数量的小控制体积,作为计算网格单元。选用隐式时间积分求解,将时间步长设置为0.02 s,求解时间设置为8 s。

    2.   结果与讨论

      2.1.   流场分析

    • 本文计算了雷诺数为200工况下的涡流场、升力、阻力系数及斯特劳哈尔数,并与其他文献进行对比。图2Re=200,求解时间为8 s时圆柱绕流的涡流图,图3图4为文献Liu[10]和Rajani[11]Re=200时圆柱绕流的涡流图。经对比可得,本文与文献中的模拟结果规律一致。圆柱下游均形成两排交错脱落的漩涡,紧靠圆柱的尾流处,漩涡幅度较小,涡量较大。距离圆柱越远,漩涡幅度越大,涡量越小。流场整体呈梯状振荡分布,具有较好的周期性。

      Figure 2.  Eddy current diagram when the Reynolds number is 200

      Figure 3.  Eddy current diagram when the Reynolds number is 200 from paper [10]

      Figure 4.  Eddy current diagram when the Reynolds number is 200 from paper [11]

      升力系数CL、阻力系数CD和斯托劳哈尔数St是圆柱绕流数值模拟中重要的特征参数[12]。升力、阻力系数可反映圆柱升力、阻力的变化规律,斯托劳哈尔数反映了绕流对柱体作用的周期性[13],对比特征参数可以验证计算的速度场结果是否合理。

      圆柱绕流的升力、阻力系数均呈周期性振荡[14],升力系数的振荡值明显大于阻力系数[15]图5为本文雷诺数为200时的升力、阻力系数变化曲线,由图5可得,在低雷诺数的层流状态下,圆柱表面的升力、阻力系数均呈周期性振荡,并且升力系数的振荡频率远大于阻力系数,符合振荡规律。

      Figure 5.  Rising resistance curve when the Reynolds number is 200

      通过计算求解,得出雷诺数为200时阻力系数和升力系数的均值分别为1.26和0.64,计算结果与1987年至2016年各文献中数据[16-18]的对比如表1所示,本文计算数据与文献中的结果基本一致,这说明该模拟方法具有可行性。

      Re=200阻力系数升力系数
      桂洪斌[14] 1.31 0.65
      Loceinte[16] 1.29 0.6
      Zhang J F[16] 1.25 0.54
      Marinez[16] 1.3 0.7
      Rockwell[16] 1.2
      顾恩凯[16] 1.29 0.74
      Braza [17] 1.35 0.5
      魏志理[18] 1.45 0.8
      本文 1.26 0.69

      Table 1.  The mean value of the drag coefficient and lift coefficient when the Reynolds number is 200

      在低雷诺数阶段,随着雷诺数的增大,涡街脱落的特征频率逐渐增大,相应的斯特劳哈尔数也呈增长趋势[16]表2为本文不同雷诺数下涡街脱落特征频率及斯特劳哈尔数数据,由表2可知,斯特劳哈尔数随雷诺数的增大而增大。方正[19]和孟元元[20]计算了Re=100时旋涡的形成和脱落,斯特劳哈尔数均为0.16,与本文雷诺数为100时的结果基本一致,说明本文计算的圆柱绕流结果是合理的、可靠的。

      雷诺数75100125150175
      脱落频率2.252253.153154.054054.954965.80541
      斯特劳哈尔数0.150150.15765750.1621620.1651653330.165868857

      Table 2.  The frequency of vortex shedding and the Strouhal number when the Reynolds number is 75~175

    • 2.2.   浓度场分析

    • 为观察低雷诺数下,圆柱绕流流场对污染物浓度对流扩散的影响,在方程(1)(2)计算结果的基础上,模拟浓度场,对比分析其涡流图和浓度场图。

      图6图9分别为雷诺数为75,求解时间为2 s和8 s时的涡流图和浓度场图。由图6图9可得,受涡流影响,污染物的对流扩散与涡流呈同步振荡分布。在求解时间为2 s时,涡流还未成型,此时表面污染物质的浓度扩散情况受涡流影响还不明显,圆柱体尾流中污染物的浓度扩散基本呈泄状分布。在求解时间为8 s时,涡流已基本成型,流场逐渐趋于稳定,通过对比可以发现,在圆柱尾流区域中,污染物质的扩散与涡流场的变化基本一致,均呈周期性的振荡分布。同时,受涡流影响,污染物质在流场中不易对流,每个漩涡中心的浓度值相对较高,而在漩涡边缘的浓度值相对较低,形成一定的浓度梯度。

      Figure 6.  Eddy current field diagram at t=2 s

      Figure 7.  Concentration field diagram at t=2 s

      Figure 8.  Eddy current field diagram at t=8 s

      Figure 9.  Concentration field diagram at t=8 s

    • 2.3.   浓度场和雷诺数的关系

    • 为了观察在低雷诺数下,雷诺数的变化对表面污染物质浓度对流扩散的影响,设置求解时间为8 s,雷诺数分别为75、100、125、150、175时,得到无量纲的平均质量值、浓度等值线图以及尾流各漩涡浓度峰值。

      表3为计算区域内不同雷诺数下的无量纲平均质量值,由表3可得,区域内平均质量值随速度和雷诺数的增大而减小。在模拟中,圆柱表面污染物质的量是一定的,当入口边界来流速度增大时,就会加快污染物质的扩散。本文的计算域为矩形,速度增大会加快污染物质从边界流失,故而整个区域内的平均质量会随着雷诺数的增大而减小。

      由于流场尾流涡街脱落,污染物质扩散的浓度场会形成交替出现的漩涡,随着雷诺数的增加,涡街脱落的频率发生变化,因此污染物质的对流扩散会受到一定的影响,浓度场发生变化。

      雷诺数75100125150175
      来流速度0.150.20.250.30.35
      平均质量0.0261560.0221490.0198480.0183150.017215

      Table 3.  Changes in the average quality with speed and Reynolds number

      图10图14为浓度等值线图,由图可得,在污染物随涡流扩散的过程中,随着雷诺数的增大,尾流的振荡幅度变小,涡流逐渐趋于平缓,影响范围缩小。

      Figure 10.  Concentration contour map with the Reynolds number of 75

      Figure 11.  Concentration contour map with the Reynolds number of 100

      Figure 12.  Concentration contour map with the Reynolds number of 125

      Figure 13.  Concentration contour map with the Reynolds number of 150

      Figure 14.  Concentration contour map with the Reynolds number of 175

      表4为雷诺数为75~175时,圆柱尾流区域浓度场从左边界到右边界各漩涡的浓度峰值,由表可得,涡流中心污染物浓度峰值的平均值随雷诺数的增大而减小。由此可得,随着雷诺数的增大,涡流效应减小,尾流漩涡趋于平缓,涡流中心的浓度峰值也逐渐降低。这说明,在Re=75~175的范围内,随雷诺数的增大,污染物质的累积效应呈下降趋势。

      雷诺数峰值1峰值2峰值3峰值4峰值5均值
      750.151920.120510.0890990.0890990.0576890.101663
      1000.118390.0869020.0869020.0554160.0554160.080605
      1250.117990.0864850.0549850.0549850.0549850.073886
      1500.0862090.0546990.0546990.0546990.0546990.061001
      1750.0866000.0551030.0551030.0551030.0236070.055103

      Table 4.  Peak value of each vortex concentration in the concentration field when the Reynolds number is 75~175

    3.   结 论
    • (1)本文模拟了低雷诺数下的圆柱绕流问题,所得流场、升力、阻力系数及脱落频率等参数与其他文献基本一致,证明通过本文模型进行数值模拟可以有效地分析圆柱绕流问题,验证了数值模型的准确性。

      (2)雷诺数为75时,在流场的作用下,污染物质的对流扩散会受到影响,污染物质的对流扩散与涡流场的变化基本一致。受涡流影响,表面附着物不易对流,并且表面附着物在每个漩涡中心的浓度值相对较高,而在漩涡边缘的浓度值相对较低。

      (3)在低雷诺数下,Re为75~175时,雷诺数的变化会对污染物质的浓度场产生影响。随着雷诺数的增大,浓度场涡流效应减小,尾流漩涡趋于平缓,涡流中心的浓度均值也逐渐降低。这就说明,在低雷诺数下,随着雷诺数的增大,涡心位置的浓度值呈下降趋势。

Reference (20)

Catalog

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return