python 字节流分段_气液两相流数值模拟分段线性界面模型

本文主要是介绍python 字节流分段_气液两相流数值模拟分段线性界面模型,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

f95c91eb693c7bbe8f8d31bca0a353cd.png

航空航天丨海洋船舶丨核电能源

石油&天然气 | 化工&过程工艺 | 环境水利

 -  多 相 的 世 界  -

先聊一聊,什么是气液两相流?

0cc97d2f7dc0b30d90ceea427a78464f.png

如果你正在喝一杯加了冰的可乐,仔细观察会发现杯子内壁不断有气泡生成、长大、脱落、上升、划过冰块,在液面处破碎……

eb6f010ad78f341181234b3a12e0a848.gif

168a400ce3851f2b77eb2908e2fc794d.png 168a400ce3851f2b77eb2908e2fc794d.png 168a400ce3851f2b77eb2908e2fc794d.png 这就是气液两相流的一个典型示例,化工、核能等领域中的气液两相流问题与可乐中的上升气泡有着相同的机理。

于是,担心可乐糖分的同时,有人还在为计算、复现这个现象,夜不能寐,甚至头秃。

4670b35f1761e7d07b16e0a981c24b38.gif

关于气液界面

气体和液体的交接处被称为气液界面。对气液界面进行描述的重要性体现在,界面两边的流体的物理性质可能差异极大,例如40摄氏度时液态水和饱和水蒸气的密度比接近20000;它的重要性还体现在,有时需要计算发生在界面处的传热、相变或者化学反应。

对气液界面性质的研究可以被追溯到两个世纪之前:

  • Young,Laplace 和Gauss 认为相间界面是一个带有表面张力的、零厚度的面。这一观点推导出一个有趣的结论,表面张力和界面两侧的压差总是处于力平衡状态。因为界面的厚度为零,所以它的质量也为零,它对力的响应是瞬时的。描述的界面被称为尖锐界面(sharp interface)。我们日常能够观察到的界面大多数是尖锐的,比如说可乐里的气泡。

  • 另一个观点由Rayleigh和van der Waals 提出,他们认为当所研究的问题的尺度和界面的尺度接近时,界面可以被认为是有一定厚度的,即扩散边界(diffuse interface)。例如,在一个单组分气液系统接近临界状态时,界面可以被认为是扩散边界。

4670b35f1761e7d07b16e0a981c24b38.gif

如何计算复现尖锐的气液界面

两相流仿真中的挑战

相方程中的对流项离散一直是两相流仿真中的挑战。从Volume of Fluid的观点出发,相方程如下:

07ba80df2f6351afe3bb1ee4411d9132.png

其中α是液相体积分数。在有限体积方法中,它的离散形式,

8b2adc0022b0739c40d7d1ec1c5c4d65.png

Uf·Sf 代表通过网格面的体积通量,αf是插值到网格面心上的液相体积分数。

αf的插值方法就是两相流仿真中的实质挑战,一众大佬为了这个事,努力了很久。

REFERENCE

感兴趣的同学可以查阅:

Flux Corrected Transport, FCT (Boris, 1973; Zalesak 1979);

TotalVariation Diminishing, TVD (Harten, 1984); 

Normalized Variation Diagram, NVD(Lenard, 1988); 

Compressive Interface Capturing Scheme for Arbitrary Meshes, CICSAM (Ubbink, 1997); 

Tangent of Hyperbola for INterface Capturing, THINC (F.Xiao, 2005);

Constrained Interpolation Profile, CIP (T. Yabe, 2001)。

OpenFOAM中的interFoam的思路和上面的文献不同,直接在相方程中加了一个压缩界面项。 以上方法得到的界面存在一定的厚度,也就是说场需要跨过几个单元才能从0变成1。 尖锐气液界面

我们向尖锐气液界面努力的原因很简单:“It’s physical”;很多情况下,它是物理上正确的。笔者基于1980年代Youngs论文中的想法,开发了一个分段线性界面(Piecewise LinearInterface Calculation, PLIC)模型,植入了开源CFD软件OpenFOAM v6。

开发的PLIC算法示意如下。

096c0179b9a3263740615178ac0cf33d.png

利用界面法向量 35fec22bfd40bbeda4b06f1614c6a2ca.png 和液相体积分数α可以重建单元内的界面位置。 界面重构后,就可以利用界面位置计算αf。 很明显,图中只展示了一种界面位置和对流体积的组合。 笔者提出了一种 “类算法(generic algorithm)” ,囊括了所有组合,大幅度简化了算法。 上升气泡算例

采用Hysing的算例来验证求解器。上升气泡的最终形状:

023d9fdcdfdcc80afc4a109c4047fd3d.png

嗯,第一眼平平无奇。

7e2cf8731561062c8516ab5e2606414e.png

最大曲率处的场和界面位置:

3c969282f2aa63edf3fc1beb21f59027.png

图中网格单元中的黑色线段即为分段线性的界面(Piecewise Linear Interface),界面的左侧是气体,右侧是液体,这样的界面可以被称为“尖锐”。所以使用PLIC方法,可以得到气液界面的精确位置。

注:图中的黑色线段使用paraview中的python脚本画出。

4670b35f1761e7d07b16e0a981c24b38.gif

结 语

感谢小伙伴们没有被令人头秃的公式和文献吓跑,一直看到这里。这篇短文主要想介绍一下两相流界面的描述,展示一下分段线性界面(PLIC)模型的能力,也想为两相流中的“尖锐界面模型”正名。在笔者读过的所有有关PLIC的文献中,只有M. van Sint Annaland曾经展示过分段线性的界面,文献附在结尾,有兴趣的朋友可以读一下。很多算法声称“尖锐界面”,但是并没有展示他们的“尖锐界面”结果。

78d04ec0f464287e9dab18a577e7f694.png

REFERENCE

作者博士工作简介:

https://www.tue.nl/en/news/news-overview/01-05-2020-stay-cool-by-phase-transition/

Mvan Sint Annaland, NG Deen, and JAM Kuipers. Numerical simulation of gasbubbles behaviour using a three-dimensional volume of fluid method. Chemical Engineering Science, 60(11):2999–3011, 2005.

863d91e8370acdb99e007a4a02477bae.gif

CFD · 分享

CFD史话

CFD中的预处理方法

OpenFOAM初学资源!前方囤积狂极度舒适预警

OpenFOAM湍流模型的实现

Kelvin–Helmholtz不稳定性

界面追踪:Level Set 与 VOF

涡结构提取

解锁新技能:OpenFOAM三维网格生成器

CFD×AI | OpenFOAM也能深度学习

气液两相流声速的直接数值模拟

关于流固耦合(FSI)模拟

CFD · 计算仿真案例

基于OpenFOAM的水翼空化数值模拟

船舶与浮冰的碰撞仿真模拟

基于重叠网格方法模拟物体高速入水

基于神经网络算法的OpenFOAM燃烧求解器

并联四喷管发动机流场结构形态及CFD计算

多相流模型在火箭发射喷水降噪系统中的应用

OpenFOAM湍流模型的实现

关于稠密气固两相流,你get到多少?

搞懂这些,你就掌握了气固两相流传热模拟的核心

flameletFoam详解及SandiaD湍流射流火焰模拟

54b9bdf516d28bdcba0b91cbc09b00b6.png

 - 仿 真 成 就 智 造 -

更多精彩内容

请关注「多相流在线」

4a16acffe0bb72b664a45fd08dcd0dae.png

工业流体仿真软件

CFD软件开发 | 行业解决方案 | 多相流测试设备

一站式多相流解决方案平台

· Multiphase Flow World ·

38eee516cd27f27fe3cd37e2bdee5ad6.png

多相流在线

网址 : www.cmfd.com.cn

在线客服:QQ 3375113520

微信客服:ID cmfd2008

这篇关于python 字节流分段_气液两相流数值模拟分段线性界面模型的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/283580

相关文章

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.

使用Python实现快速搭建本地HTTP服务器

《使用Python实现快速搭建本地HTTP服务器》:本文主要介绍如何使用Python快速搭建本地HTTP服务器,轻松实现一键HTTP文件共享,同时结合二维码技术,让访问更简单,感兴趣的小伙伴可以了... 目录1. 概述2. 快速搭建 HTTP 文件共享服务2.1 核心思路2.2 代码实现2.3 代码解读3.

Python使用自带的base64库进行base64编码和解码

《Python使用自带的base64库进行base64编码和解码》在Python中,处理数据的编码和解码是数据传输和存储中非常普遍的需求,其中,Base64是一种常用的编码方案,本文我将详细介绍如何使... 目录引言使用python的base64库进行编码和解码编码函数解码函数Base64编码的应用场景注意

Python基于wxPython和FFmpeg开发一个视频标签工具

《Python基于wxPython和FFmpeg开发一个视频标签工具》在当今数字媒体时代,视频内容的管理和标记变得越来越重要,无论是研究人员需要对实验视频进行时间点标记,还是个人用户希望对家庭视频进行... 目录引言1. 应用概述2. 技术栈分析2.1 核心库和模块2.2 wxpython作为GUI选择的优

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python+PyQt5实现多屏幕协同播放功能

《Python+PyQt5实现多屏幕协同播放功能》在现代会议展示、数字广告、展览展示等场景中,多屏幕协同播放已成为刚需,下面我们就来看看如何利用Python和PyQt5开发一套功能强大的跨屏播控系统吧... 目录一、项目概述:突破传统播放限制二、核心技术解析2.1 多屏管理机制2.2 播放引擎设计2.3 专

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

Java的IO模型、Netty原理解析

《Java的IO模型、Netty原理解析》Java的I/O是以流的方式进行数据输入输出的,Java的类库涉及很多领域的IO内容:标准的输入输出,文件的操作、网络上的数据传输流、字符串流、对象流等,这篇... 目录1.什么是IO2.同步与异步、阻塞与非阻塞3.三种IO模型BIO(blocking I/O)NI

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤