基于大涡模拟下的三种足球空气动力学分析及运动轨迹探讨
□ 张旋（中国社会科学院大学体育教研部 北京 100000）
摘要： 本文利用计算机数值模拟方式分别对普通足球、团队之星足球、C4D NIKE顶级比赛用球进行空气动力学分析，通过SolidWorks构建模型，PointWise划分网格及计算域，利用WALE亚格子湍流模型封闭方程组以加强计算精度对以往数据进行修正，以对不同表面沟壑类型的足球升力、阻力和基于Q准则的涡量图进行分析。计算出在不考虑足球自转和角度影响下的足球运动轨迹偏离情况，通过研究发现不同足球表面沟壑类型在不同速度飞行下其影响关系及涡脱落造成的飘球产生机理，进一步利用傅里叶变换和对速度对时间进行积分获得的位移量可较为清晰的判断足球运动轨迹和规律，由此得出普通足球香蕉球弧度最大，团队之星足球电梯球弧度最大，而C4D NIKE顶级足球速度最快，直线威胁最大并制作飞行轨迹示意图对守门员预判来球方向做出预警，同时对后续研究者提供更多数据参考。
关键词： 大涡模拟 足球 空气动力学
1、前言
1.1、概述
对足球进行运动空气动力学分析是国外帮助球队获得比赛胜利，帮助球员改进训练水平的一项重要工作，但在我国该方面的研究人才甚少，绝大多数实验室在采用风洞试验对部分运动项目进行相关研究，其成本较高并不能有效普及于绝大多数运动领域。因此本文通过计算仿真模拟方式对足球相关物理量和飞行轨迹进行量化分析，旨在提供一种新的分析方法并提供足球相关数据予后续相关工作的研究者。
1.2、研究目的与任务
（1）研究目的。
通过LES计算模拟仿真获得不同足球空气动力学特性，并借此对不同足球沟壑类型及其变化造成的物理量变化进行分析，通过对比试验量化足球飞行轨迹。
为进行无碱二元复合驱中表面活性剂组分对驱油效果的影响研究，建立9注16采的3层纵向非均质概念模型见图1。注采井距为200m，有效厚度为18m。网格数为 76×76×5，网格尺寸为 20m×20m×6m，地层孔隙度为0.2，三层的渗透率分别为700×10-3μm2、500×10-3μm2和 300×10-3μm2。初始水饱和度为定值0.32，地层压力为定值9.42MPa。油井含水率达到95%，开始注二元，注入0.6PV，后续注水至98%。
（2）研究任务。
①研究模型及部分方程研究；
下面证明ρ是好同余。设aL*E,则存在CS°(a)和CS°(e)使得a°L*e°和Λe=Λa。记e°λe, e°λeixx°和e°λeiyy°分别为λa,(ex)°和(ey)°。则
1.推行法律服务工作机制，实现社区与职工群众和谐共处。如今企业职工的思想观念呈现多元化趋势。为解决职工思想观念和日常工作生活中的主要问题，企业必须要依法履行管理服务职责，要站在规范理顺职工维权行为的角度，完善职工普法工作规划及企业法律服务工作相关规定、法律服务工作考核细则，抓好法制宣传教育阵地建设，做到法律服务和日常管理服务工作同部署、同运行、同考核。同时，企业可通过聘请资深法律专家召开专题法律讲座、开设法律服务专栏、定期组织开展“阳光咨询”活动和区域法律座谈咨询会等形式，面对面解疑释惑，多渠道宣传法律知识，增强职工依法办事、依法维权的意识和法律服务的针对性和实效性。
②足球沟壑类型对升阻力的影响；
2）测量放线。通常，会根据施工图纸所标注的各个坐标点在施工现场进行标注，在一般情况下，需要每隔30cm打1根边桩，但是如果打桩地区的地形变化较大，或是处于转弯路段，需要根据实际情况进行调整。
就整体而言，火法回收耗能较大，部分方法污染严重，不符合现代工艺要求，方法较为陈旧，鉴于大规模生产要求，火法回收是首选方法。
④不同足球偏移轨迹的研究。
2、研究方法
2.1、文献资料法
查阅1980-2018年相关国内外文献研究300余篇，对湍流模型及空气动力学主要问题提取宝贵经验与数据。
如图3所示为普通足球及2015年生产的团队之星足球、2016年生产的C4D NIKE顶级比赛用球在雷诺数为Re≈3.9×105时的阻力系数图，平均阻力系数分别为0.1276、0.1406和0.1150。在以往实验中研究人员已经发现对于光滑圆球在低速时受到的阻力更大，而在高速时受到的阻力更小，足球以22m/s前进时并未达到高速要求，同时可判断其表面粗糙程度。因此实验中的三种足球以C4D NIKE顶级足球表面粗糙程度最小，随后依次为普通足球和团队之星足球。
2.2、专家访谈法
通过与国内外专家的学习论证，确认操作方案的可行性和理论依据的标准。
2.3、数理统计法
通过数理统计方法对后处理中出现的物理量进行定量分析，对部分数据进行散点图和动画处理。
SPOC混合教学模式按照教学目标重新设计课程，课程设计围绕知识点展开。SPOC平台中的在线学习系统设计符合学生个性化学习的特点与要求，教学评价系统设计对学生的在线学习进行督促与激励，同时SPOC课堂能发挥交流与互动多的优势，对学生加强实践指导，帮助提升学生的实践应用能力。通过SPOC教学改革，促进教师从知识的简单传授者转变为知识学习的引导者与自由讨论的组织者，学生也从被动学习的接受者，转变为自我学习的管理者与知识创新的实践者。
2.4、模拟法
建立实验模型，设置计算机工况，确认相关理论，根据三者相互关系，进行物理与数学模拟。
3、结果与分析
3.1、模型的建立及部分方程方案
为保证计算精度采用增长率1.15，15层边界层，同时FLUENT采用双精度模拟。
本文采用三种不同足球模型如图1所示，其表面沟壑位置不同造成不同物理量及飞行轨迹变化，足球沟壑采用1mm圆角以真实还原足球表征，采用直径0.22m国际标准足球，表皮粗糙程度不在本文模拟范围之内，因此与实际情况会产生部分误差，经验证其误差值极小可以忽略不计，足球远场分别为设为5倍足球直径和10倍足球直径。
图1 三种足球模型
（2）网格划分方案及方程。
有人研究有限游戏，热闹于有限游戏；有人研究无限游戏，潇洒于无限游戏；人类才能够延续到今天，才有了更好的可能性。
此外，将足球表面缝隙加深也会产生不同变化，如普通足球表面沟壑加深至2mm圆角，则其阻力系数在相同雷诺数下有所升高，即表面沟壑深度直接影响粗糙程度，造成粗糙度的增加，其在较低出脚速度时足球飞行速度与其成正比同样速度较慢。
（1）Cf=τq/ρu2（2）1/Cf 0.5=1.7+4.15×log10(Re×Cf)（3）
（1）模型的建立。
图2 标准足球混合网格剖面图
3.2、不同足球类型的阻力变化及其应用
（1）足球阻力系数及其对电梯球飞行速度的影响。
图3 不同足球阻力系数变化折线图
总而言之，在小学数学教学的过程中，通过利用多媒体科技，可以更为直观地化解教学过程中出现的重点和难点；并能够将生活中的数学和课堂中的数学结合起来，在学校和家庭之中，为小学数学的教学营造出一种良好的环境和氛围，从而保障学生的学习质量，建立科学的教学方法，达到促进教学效果提升的目的。
本文采用混合网格以减少网格数量，增加质量并提高运算速度与精度，如图2所示为标准足球网格划分方案图，由于利用LES大涡模拟方法是通过低通滤波将流动分为大小两种尺度的涡，大尺度涡使用Navier-Stokes方程进行求解，小尺度涡使用WALE亚格子尺度模型进行方程组的封闭且由于来流速度设定为22m/s在计算中未能引起流体密度的足够变化，因此流场属于不可压缩。
③三维涡量云图对比分析；
而足球电梯球的防守难度主要由两种情况构成，其一为足球飞行速度、其二为足球变化弧度，而阻力系数的大小将决定足球整体速度变化，相对于实验的三种足球而言，粗糙足球表面边界层提前分离后由于卡门涡街效应会降低足球速度，因此可以断定C4D足球飞行速度最快较难防守，而团队之星足球则飞行速度最慢易于防守。
一是学习不足。医院管理人员多从临床医务工作者中选拔，临床业务技能强，但在管理专业知识方面却相对薄弱，加之医院管理的特殊性，不少干部均为“双肩挑”，学习中埋头拉车而疏于抬头看路，重业务、轻管理的现象客观存在，需要平衡临床业务与管理投入。
同时由于采用LES大涡模拟，必须要求网格能够满足求解粘性影响区域的流动，故本文采用Y+=1进行计算第一层网格高度为1e-05m，利用牛顿迭代求解公式如下：
（2）足球涡脱落现象对阻力系数变化的验证。
为进一步验证以上三种足球阻力系数变化具有可信性，我们利用Q准则对三种足球分别作出三维涡量图，根据涡脱落速度判断阻力系数的大小，其定义公式如下：
Q=-1/2×(Velocity u.Gradient X2+Velocity v.Gradient Y2+Velocity w.Gradient Z2)-Velocity u.Gradient Y×Velocity v.Gradient X×Velocity w.Gradient Z
如图4所示为为Q=0.001时的三种足球三维涡量图及其流线图，可以清晰的发现团队之星足球较普通足球、C4D nike涡脱落速度更快，而涡脱落将造成边界层提前分离，进而使不同位置压差随之增大，造成阻力系数的增大由此可证明以上平均阻力系数的变化。
图4 足球三维涡量图及流线图
3.3、足球电梯球弧度变化规律及下坠时间节点分析
图5 不同足球升力系数折线图
在足球飞行过程中，存在一种特殊飞行轨迹，即球不按照平滑曲线进行运转，而是飘忽不定，其肉眼观察似乎并不符合物理规律。事实上，从圆柱绕流的类似问题早已发现，对于平滑物体在一定雷诺数下通过斯特劳哈尔数公式St=fD/v可以得到涡脱落频率其取值为0.14-0.27之间，同时有研究者发现在Re＞2×105以上的情况，St数和Re的相关度很低。对于足球飞行而言，其Re远高于2×105，也因此足球产生特殊飞行轨迹可以猜测为与球体表面粗糙程度相关，而当雷诺数为3×105-3×106时，有规则的涡街便不再存在，因此根据计算足球速度一般位于此雷诺数区间，涡脱落频率不再为常数，造成足球运动轨迹的不规则变化。
如图5所示为三种不同足球升力系数折线图，在低速下，C4D平均升力系数为0.0164，普通足球平均升力系数0.0061，团队之星平均升力系数0.078，由于足球表面粗糙程度造成的升力不规则变化进而出现特殊飞行轨迹且观察曲线变化规律可知，普通足球高度飞行轨迹最为飘忽不定，其次为C4D NIKE顶级足球，而团队之星足球在高度上飞行轨迹更为规律，这也是飘球的产生机理。
同时通过对升力进行傅里叶变换，可以观察其频率分布，而频率分布越接近某一具体值说明卡门涡街效应越明显，规则的卡门涡街将造成足球在某一时间节点产生急速下坠，以便对电梯球进行预判。
图3和图4分别为在无气流辅助和有气流辅助条件下喷雾场内雾滴运动轨迹，图3和图4中雾滴的颜色代表雾滴粒径的大小(单位mm)。比较两图，可以看出在有气流辅助时，气流流场改变了雾滴的运动方向，引导雾滴运动轨迹明显偏向地面，从而有效地减少了雾滴的飘移潜能。
图6 升力系数傅里叶变换
如图6所示分别为普通足球、C4D NIKE顶级足球和团队之星足球，其中频率分布代表了非对称涡脱落的分布，而越接近周期性涡脱落则由涡脱落造成的球的震荡导致足球产生快速下坠愈加明显，因此普通足球电梯球即出现最大拐点时间最快，其次为C4D NIKE顶级足球和团队之星足球。
再者在不考虑飞行过程中足球自转速度带来的影响，即假设足球始终沿Y轴方向飞行时，以22m/s给定足球初速度，则分别计算足球飞行高度和侧向位移幅度如图7所示：
图7 足球飞行轨迹偏移位置图
图中分别为标准足球、C4D NIKE足球和团队之星足球的不同方向偏移轨，由图可知足球沿Y方向飞行受到侧向力影响所造成的偏移分别为 1.06m至 -0.79m、0.1m至 -0.8m、0.13m至 -0.24m；飞行高度偏离值分别为 0.85m至 -0.74m、0.1m至-0.4m、0.04m至-0.89m、由侧向偏移值即香蕉球最大弧度以普通足球居首，C4D NIKE顶级足球次之，团队之星足球弧度最小；且由升力造成的电梯球以团队之星弧度最大，普通足球次之，C4D NIKE顶级足球最小。且其飞行44m时侧向弧线最大拐点分别出现在距球31.9m、13.64m和13.42m；而高度方向最大弧线拐点分别出现在18.04m、24.2m和39.82m与3.3图6所得结果相同。
在使用密度基求解器对空气湿度和压力变化时普通足球升阻力的变化趋势进行研究发现，空气湿度和压力的增加将使球速降低，其平均升阻力系数由于相对粗糙程度的增加而有所降低，也就造成了足球的电梯球和香蕉球产生弧线逐步变小，但是由于仿真模拟的误差及足球本身速度较慢以及大气压的局限性，其并不会产生足够的肉眼可见的弯曲变化，浮动值极小，在此仅做参考。
4、结论
（1）C4D NIKE顶级足球表面粗糙程度最小，随后依次为普通足球和团队之星足球，而粗糙程度变化导致足球表面边界层分离提前，导致压差变化，进而影响到球飞行速度。因此从飞行速度上可以断定C4D足球飞行速度最快较难防守，而团队之星足球则飞行速度最慢易于防守。同时足球表面沟壑深度会增大足球粗糙程度进而使足球在低速下产生较大阻力，减慢球体飞行速度；
（2）由于足球在一定雷诺数下飞行造成涡脱落频率不再为常数，致使足球运动轨迹的不规则变化是飘球产生的主要机理。根据升力变化可知普通足球高度飞行轨迹最为飘忽不定，其次为C4D NIKE顶级足球，而团队之星足球在高度上飞行轨迹更为规律；
（3）通过傅里叶变换和对速度对时间进行积分获得的位移量可较为清晰的判断足球运动轨迹和规律，由此判断出普通足球香蕉球弧度最大，团队之星足球电梯球弧度最大，而C4D NIKE顶级足球速度最快，直线威胁最大；
（4）空气湿度和压力的增加将使球速降低，其平均升阻力系数由于相对粗糙程度的增加而有所降低，致使足球电梯球和香蕉球均不易产生更大弧度，但是由于仿真误差的存在及其变化量较小，在此仅做参考。
参考文献：
[1]杨春卉,袁志华,梁振刚.乒乓球反弹动态特性的仿真研究[J].计算机仿真,2014(10).
[2]王奇志,杨晓晓.乒乓球轨迹预测的研究与仿真[J].计算机工程与科学,2013(02).
[3]杨华,关志明.基于ODE的乒乓球运动轨迹仿真研究[J].计算机仿真,2011(09).
中图分类号： G843
文献标识码： A
文章编号： 1006-8902-（2019）-08-214-3-ZQ
作者简介： 张旋（1989-），男，籍贯：河北，讲师，研究方向：流体力学研究工作，包括物体建模、气动力学分析、数值计算等实际应用问题。
标签：大涡模拟论文; 足球论文; 空气动力学论文; 中国社会科学院大学体育教研部论文;