PRF | 宾州州立、南科大杨翔、张雯等:粗糙壁湍流的低维建模新范式

粗糙壁湍流的低维建模新范式

Rough-wall modeling on a low-dimensional manifold

Shyam S. Nair, Robert F. Kunz, Wen Zhang (张雯), Xiang I. A. Yang (杨翔)

引用格式:Nair Shyam S, Kunz Robert F, Zhang Wen, et al. Rough-wall modeling on a low-dimensional manifold[J]. PHYSICAL REVIEW FLUIDS, 2025, 10(6):064606.

编者按

粗糙壁面流动广泛存在于自然与工程场景中,但由于粗糙面结构复杂多变,传统依赖于经验参数的阻力预测方法在泛化能力与精度上存在明显不足,且不具备生成能力。本文提出一种基于低维潜空间学习的粗糙壁面建模框架,不依赖人为选取的粗糙度统计量,而是通过深度卷积自编码器自动提取粗糙面本征特征,并建立与等效砂粒粗糙度高度(ks)的关联模型。该方法不仅预测准确,还具备粗糙面重建与优化设计潜力,为湍流边界层建模提供了全新思路。

一、 研究背景

湍流与壁面粗糙度之间复杂的相互作用,一直是壁面流动研究中的核心挑战。粗糙表面扰乱近壁流动结构(如条带与涡),显著改变边界层动量输运特性,从而增加阻力。如何有效表征粗糙度并准确预测由此产生的阻力,对飞行器设计、舰船水动力、天气预报与地表水力学等领域均具有重要意义。

传统的粗糙度建模依赖于对表面统计量(如均方根粗糙度、高斯偏度、有效坡度等)的经验选取与拟合。然而,这类参数的选择具有较强主观性,难以泛化,尤其是在面对不同类型、多尺度甚至非规则粗糙面的情形下。

二、研究方法

该研究提出了一种新颖的粗糙壁面建模思路——将复杂粗糙面的形貌压缩映射到一个仅包含三个变量的低维流形上。核心方法包括:

  • 自动编码器(Autoencoder)压缩粗糙面形貌:通过深度卷积神经网络将93种实验与DNS获得的粗糙面高度分布,映射到三个无量纲的潜变量,实现低维数据驱动特征提取。

  • 前馈神经网络预测ks值:将潜变量作为输入,通过训练前馈神经网络,预测等效砂粒粗糙度高度ks,这一参数是用于量化粗糙度引起阻力的重要水动力尺度。

  • 构造新粗糙面并验证模型:通过在潜变量空间内生成未见粗糙面样本,使用DNS计算其真实ks,与模型预测值对比,验证模型泛化能力。

图1:卷积自编码器与前馈神经网络粗糙壁面模型的示意图

三、结果与分析

  • 重构效果优秀:Autoencoder能够准确还原原始粗糙面结构,并在保留统计量(如krms、Sk、ka)方面表现良好。

  • 建模准确性优于现有方法:相较于Flack等(2020)与Forooghi等(2017)提出的粗糙度经验公式,本模型在多个测试样本上的预测误差更小,R2拟合优度达0.89。

  • 对未见粗糙面具备预测能力:模型对潜变量空间中新生成的粗糙面进行ks预测,结果与DNS计算吻合良好,表明该低维流形具备一定的生成性与泛化性。

  • 参数物理意义仍在探索中:潜变量与传统粗糙度统计量之间相关性较弱,说明其提取的是非线性高阶特征,这为未来的解释性研究留出空间

图2:粗糙表面样本的输入(左)与重建结果(右)对比
图3:不同粗糙壁面模型对 ks/ka 的预测结果: (a) Flack 等人 [23],(b) Forooghi 等人 [24],(c) Yang 等人 [32],以及 (d) 本研究
图4:对训练集之外的粗糙表面的由 DNS 得到的流向平均速度剖面与前馈神经网络粗糙壁面模型预测的剖面对比

四、研究意义

本研究展示了将粗糙壁面问题映射到“潜空间”的可行性,突破了传统粗糙度统计建模的局限。该方法不仅提升了建模精度与可扩展性,更为未来粗糙面优化设计提供了数据驱动的设计变量空间。潜变量的可控性也为逆向设计(如阻力最小化粗糙面形貌优化)奠定基础。

相关论文:

BLM| 南科大、宾州州立张雯、杨翔等:湍流边界层中大尺度二次涡结构演化的大涡模拟与建模 https://link.springer.com/article/10.1007/s10546-025-00904-1

JFM | 南科大、宾州州立张雯、杨翔等:摩擦分解积分方法及其在粗糙壁流动中的扩展https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/integral-methods-for-friction-decomposition-and-their-extensions-to-roughwall-flows/2E30FE1D68980F06AA49B9FB464A2914

JFE | 南科大、宾州州立张雯、杨翔等:探索通用粗糙壁模型 https://asmedigitalcollection.asme.org/fluidsengineering/article/145/10/101302/1164127

JFM | 南科大、宾州州立张雯、杨翔等:几何对称粗糙面上方的非对称二次流研究https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/asymmetric-secondary-flows-above-geometrically-symmetric-surface-roughness/4EE0100A56EBC3C3DEAD5F14D87F3D6E

JFM | 南科大、宾州州立张雯、杨翔等:深层均匀城市冠层流动中Raupach 混合层类比的证据https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/evidence-for-raupach-et-als-mixinglayer-analogy-in-deep-homogeneous-urbancanopy-flows/2B00522F7284DB3E69D64EA4E188C4FA

公众号原文链接(文末附论文资源):

PRF | 宾州州立、南科大杨翔、张雯等:粗糙壁湍流的低维建模新范式

注:文章由作者原创供稿,并获得作者授权发布。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/1220279.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Infoseek 媒介投放系统技术实现:基于与辉同行风波的风险防控架构设计

Infoseek 媒介投放系统技术实现:基于与辉同行风波的风险防控架构设计 摘要:本文结合 2026 年初与辉同行 “皖西麻黄鸡” 打假风波,围绕品牌媒介投放的风险防控需求,拆解 Infoseek 字节探索媒介投放系统的技术架构,重点…

疆鸿智能MODBUS TCP转PROFIBUS:网关智构精密组装新脉络

疆鸿智能MODBUS TCP转PROFIBUS:网关智构精密组装新脉络在工业自动化系统的演进中,数据流的无缝贯通是实现智能制造升级的关键。其中,工业通信网关扮演着不可或缺的“协议翻译官”与“数据调度者”角色。本次电子元件精密组装流水线升级项目&a…

2026年评价高的防静电地板公司推荐:水泥纤维网络架空地板、活动架空地板、玻璃防静电地板、硫酸钙防静电地板、通风防静电地板选择指南

2026年医疗机构适用防静电地板厂家推荐行业背景与需求洞察据《2026-2030年中国机房防静电地板行业发展白皮书》显示,2026年国内机房防静电地板市场规模突破120亿元,年复合增长率达8.7%。随着医疗机构数字化转型加速,…

记重要需严格

1、必须7-8小时睡眠 后面租司近,或去深; 2、不炒cai,最多zheng3 3、关于托伐要少喝了 4、qian bianyong bianzan 5 wa ----

2026年评价高的耐火砖公司推荐:耐火材料推荐/耐火材料电话/耐火砖供应厂家/耐火砖厂商/耐火砖厂家/耐火砖厂家电话/选择指南

四川电力行业定制耐火砖优质品牌推荐行业背景与需求洞察据《2026中国耐火材料行业发展白皮书》数据显示,2024年国内耐火材料市场规模突破1200亿元,其中电力行业需求占比达18%,年复合增长率保持在6.2%以上。 电力工业…

计算机毕业设计springboot协同过滤的就业系统的设计与实现 基于Springboot框架的就业推荐系统设计与实现 Springboot协同过滤技术在就业平台中的应用与开发

计算机毕业设计springboot协同过滤的就业系统的设计与实现qd11f(配套有源码 程序 mysql数据库 论文) 本套源码可以在文本联xi,先看具体系统功能演示视频领取,可分享源码参考。随着互联网技术的飞速发展,就业市场也逐渐向数字化、智…

2026美团代运营公司实力榜:这5家真正懂外卖流量!三十六行网络科技(阜阳分公司)领跑内容与策略赛道

在餐饮行业数字化竞争日趋激烈的今天,美团点评平台已成为商家获取客流、塑造品牌的核心战场。然而,复杂的平台规则、多变的营销玩法与精细化的运营需求,让许多餐饮商家倍感压力。选择一家专业、靠谱的代运营服务商,…

JWT(JSON Web Tokens )简洁说明

JWT(JSON Web Tokens )简洁说明实施 Token 验证的方法挺多的,还有一些标准方法,比如 JWT,读作:jot ,表示:JSON Web Tokens 。JWT 标准的 Token 有三个部分:header(头部) payload(数据) signature(签名)…

计算机毕业设计springboot校园学生健康管理与服务系统 基于Spring Boot的校园学生健康管理系统设计与开发 Spring Boot框架下的校园学生健康管理服务平台构建

计算机毕业设计springboot校园学生健康管理与服务系统29o2z9 (配套有源码 程序 mysql数据库 论文) 本套源码可以在文本联xi,先看具体系统功能演示视频领取,可分享源码参考。 随着社会的快速发展,人们对健康的关注度不断提升&…

国歌下载地址

https://www.gov.cn/guoqing/guoge/index.htm

Mysql8 允许远程登录

root@vultr:~/mysql8# mysql -uroot -p Enter password: Welcome to the MySQL monitor. Commands end with ; or \g. Your MySQL connection id is 11 Server version: 8.0.29 MySQL Community Server - GPLCopyrig…

2026年10款主流项目管理软件精选推荐|高效落地团队协作

选对项目管理软件,能大幅降低沟通成本、精准把控项目进度,让团队从繁琐的手动跟进中解放出来。精选10款主流易上手的项目管理工具。 进度猫(轻量进度可视化工具) • 核心优势:甘特图思维导图双向联动,关键路…

告别 DIFY 构建卡顿 / 失败!巧用 Makefile 一键构建本地 WEB 镜像,附核心避坑指南

在基于 DIFY 进行二次开发的过程中,我们常常需要修改其前端 WEB 源代码(比如自定义页面样式、新增前端功能),修改完成后就需要重新构建 WEB 镜像来验证效果。 官方 DIFY 在根目录(包含 Makefile 文件的目录)封装了便捷的构建命令,无需我们手动编写冗长的docker build -…

36W电源方案LP3798ESM+LP10R060SD(12V3A)隔离型适配器/充电器电源方案分析(成本低、体积小)

LP3798ESM(ASOP6,准谐振 QR 反激 PWM 控制器) LP10R060SD(SOP7L,60V 同步整流 SR),构成 90–264Vac 全压输入、12V/3A(36W)隔离适配器 / 充电器方案,峰值效率…

基于springboot的机票预定系统的设计与实现 计算机毕业设计选题 计算机毕设项目 前后端分离【源码-文档报告-代码讲解】

🍊作者:计算机毕设匠心工作室 🍊简介:毕业后就一直专业从事计算机软件程序开发,至今也有8年工作经验。擅长Java、Python、微信小程序、安卓、大数据、PHP、.NET|C#、Golang等。 擅长:按照需求定制化开发项目…

java基础-InputStream 类

InputStream 是 Java 中所有字节输入流的抽象基类,位于 java.io 包中。它定义了读取字节数据的基本方法。一、核心特性抽象类 - 不能直接实例化,需要通过子类实现字节流 - 以字节(byte)为单位读取数据单字节读取 - 最基本的读取单…

2026年四川丧葬一条龙哪家强?四川公墓咨询、绵阳公墓推荐——孝爱领衔五大品牌,温情守护生命终章

引言 生命的落幕,需要庄重与温情的守护,殡葬服务的品质直接关乎逝者尊严与生者慰藉。在四川殡葬行业规范化发展的当下,如何挑选正规专业的丧葬一条龙服务,成为许多家庭的难题。为提供客观参考,我们依托行业协会标…

DSU

DSUstruct DSU {int n;vector<int> p, sz;DSU() : n(0) {}DSU(int n_) { init(n_); }void init(int n_) {n = n_;p.resize(n + 1);sz.assign(n + 1, 1);iota(p.begin(), p.end(), 0);}int find(int x) {if (x !=…

2026年架空地板厂家推荐:机房架空地板、水泥纤维网络架空地板、活动架空地板、玻璃防静电地板、硫酸钙防静电地板选择指南

2026西南地区架空地板优质厂家推荐榜根据《2026年中国防静电与架空地板行业发展白皮书》显示,2026年国内架空地板市场规模达180亿元,年复合增长率12.3%,其中西南地区以数据中心、智能工业园为代表的新兴场景需求增速…

BIT

BITstruct BIT {int n;vector<int> tr, vis;int tag;BIT(int _n = 0) { tag = 1; if (_n) init(_n); }void init(int _n) {n = _n;tr.assign(n + 1, 0);vis.assign(n + 1, 0);tag = 1;}int lowbit(int x) { retu…