Posts made by bestucan

@zhe DPMFoam.C 这个的108行在算粒子，下面就开始用算过的信息更新alphac了。

108行的定义在这里 KinematicClud.C 673行。这里用到的solve函数就在这个文件的91行。

然后这solve里一大堆追着函数名找，大概就是particle distribution function的过程了。我猜的

@zhe 我说错了，本来想说alpha1+alpha2不等于1。就是括号里的对立面

（如果alpha1 + alpha2 = 1，那么rho = alpha1*rho1 + alpha2*rho2，未找到alpha1和alpha2的关系。）

这个确实不能说 alpha1 和 alpha2 是未归一化的。
但是这个不影响后面的思路，虽然思路很模糊。
其实绝对正确的平均密度公式是：(alpha1 * rho1 + alpha2 * rho2) / (alpha1 + alpha2)。无论alpha1、alpha2有没有归一化都成立。

还有，刚才又回看了一下。你的第一个帖子说

密度rho應該已經經過alpha1*rho1 + alpha2*rho2得到的整體流體密度

如果这个表达式表达整体流体的密度 ，那是不是就意味着aplha1+alpha2=1?
如果 alpha1+alpha2\=1，那么剩下的体积分数代表的密度不是没算进来。但是代码中并未有限制alpha1和alpha2关系的地方。alphac倒是限制为1-xxx。

不考虑求解器。多相流的平均密度应该是连续相和离散相的平均alphac * rhoc + alphad * rhod，其中alphac + alphad = 1。

MPPIC方法在OpenFOAM里计算好像用的是particle parcel？所有都没有见 alphad。

我说的太外行就不用管了，我没研究过，就是插一句嘴。
附上一个某个帖子的总结：

忘了哪个帖子了，是这个论坛的。

引用自己的两段微信聊天记录：

结构网格和非结构网格的区别在于是否可以根据网格的索引编号确定相邻网格。 结构网格的索引号可以判断网格的相邻网格和面归属，非结构网格不可以，要另存储。是不是完全垂直只会影响算面通量要不要乘以面向量

最后一句换为是不是完全垂直只会影响算面通量要不要乘以带夹角的面向量。

是的，fluent用的也是非结构网格。如果用icem画过图，一般在导出到fluent时都会转换成非结构网格才能导出。因为非结构避不开，结构的可以转换成非结构的，所以基本都是基于非结构网格的。结构网格比较配有限差分，但是现在有限元更常见。约翰。安德森的《计算流体力学基础及其应用》主要讲结构网格配有限差分，要建贴体坐标系，要坐标系变形，主要例子是算飞机翼型，应该说应用范围没有有限元广泛。openfoam用blockmesh生成的网格，或者其他格式转过来能用的网格，在constant/polymesh里，就有定义面归属（owner）和相邻网格（neighbour）的文件，所以是非结构的。总之，fluent和openfoam都是基于非结构的，来了结构网格，先转成非结构才能算。

可能有理解错的，但应该能解除部分疑惑。

说两个猜想。
如果是alphc(fake)*rho = alpha1*rho1 + alpha2*rho2。
那得是alpha1、alpha2都是未归一化的值。
（如果alpha1 + alpha2 = 1，那么rho = alpha1*rho1 + alpha2*rho2，未找到alpha1和alpha2的关系。）
这种情况下rho = (alpha1*rho1 + alpha2*rho2)/alphac(fake)即密度=质量/体积

实际上代码里的rho(code) = alpha1*rho1 + alpha2*rho2，即上面的质量。
即rho(code)是体积为alphac(fake)的混合物的质量
即rho(code)是单位alphac(fake)的质量

下面开始猜：
alphac(code)是归一化过的，在代码里有显示，是1-xxx。
那么猜测：与alpha1和alpha2类似的数值alphac(real)隐含在alphac(code)中。
即alphac(code)=alphac(real)/alphac(fake)

那么rho(code) * alphac(code)代入以上定义就等于：(alpha1*rho1 + alpha2*rho2)/alphac(fake)*alphac(real)
即rho*alphac(real)
(rho(code) * alphac(code) -> rho*alphac(real))

上面是猜rho可能的定义。下面是猜有的地方乘alphac，有的地方不乘。

上面说了alphac是归一化过的，则rho(code)是单位体积内的，这种单位内的量容易搭积分的便车补上计算。对流项没有考虑alphac，因为对流计算过程涉及体积(速度即单位时间内的体积)。

至少把你的现有进展说一下，进行到哪一步了。particle distribution function 是代码里的还是文献里的，在代码哪里。你看了什么文献、哪个代码却没有找到。好歹让想帮你的人把你踩过的坑给避开了。

用 sudo fdisk -l 看看 swap 分区多大。如果很小或者没有，可能是因为内存。如果很大可以排除内存问题。swap是虚拟内存，物理内存不够用会把硬盘当内存用，只会很慢不会崩。

也可以在编译过程中开着 system monitor 看着内存和 swap 使用状况。

当初装linux的时候分区选择了lvm方式，还能调整swap分区大小。如果不是lvm而是直接分区。只能重装系统修改分区大小。

然后你的源代码放在OpenFOAM1，1??? 不知道会不会影响。

估计你的源代码是用apt下载的，或者是哪位大佬的祖传代码，所以默认安装位置改成了/opt。

各种教程的编译安装位置基本都在$HOME\OpenFOAM，不会有权限问题。
调整 OpenFOAM 的安装位置还是挺麻烦的。

@RebelYoung 这种错误可能是数据本身有问题，不一定是paraview的问题。

如果想确定，可以去paraview的官网下载最新的paraview。如果连它都崩，那就是数据的问题。

https://cfd-china.com/topic/3453/运行pisofoam的时候-卡住不动了

里面有个可以尝试的方法。

在求解器的各处随机插入
cout << "run in line xx"
重新编译这个求解器，看看具体在哪一步之后就不输出了。

根据描述，断点续算可以，但是直接突破不过去，而且只是不输出而不是不计算。应当是某个变量积累超过输出语句的判断范围。所以还在算但不输出。所以上面提到的具体哪一步就是修改了这个变量超出了范围。

更精准的定位错误就要看输出残差的函数怎么写的了。

OF 的 debug 模式编译有更强大和方便的调试功能。

@林之流风 在 能否在集群的个人文件夹内安装OF 中说：

1、用集群的OpenFOAM编译。之前不敢用集群的OpenFOAM编译自己个人文件夹的求解器，因为我不知道 $FOAM_USER开头的这些环境变量具体的是什么目录，

如果是个人帐户，怎么折腾都行，不用怕搞坏系统，因为没权限。

在终端键入env后就可以看到这些环境变量了，其实用集群的个人账户登录，这些目录都在自己的个人文件夹里，所以接下来和在虚拟机操作类似了。

如果只是看某个环境变量 echo $FOAM_INST_DIR 就行。

2、在集群个人文件夹内安装OpenFOAM（如果想用其他版本的OF的话）。开始报错说我缺少flex，我找集群管理员要了root账号，因为我们集群的系统是Centos，所以我键入yum -y install flex 来安装flex。

yum 的 -y 最好不用，这个是默认 yes 。有时候安装软件会更新一些软件再卸载一些软件。还是看看没啥大影响再 yes。

flex安装成功后，登陆个人账户，在终端键入env后可以看到 $FOAM_INST_DIR是在我个人文件夹内的目录。编译OF，经过了很长时间编译成功，

按说编译时候缺包很难只缺一个包，你说编译成功了，怎么确认的呢？编译log有没有保留？有没有搜 log 里的 warning 和 error？ 当然不考虑这些兴许也能往下走，但是出了问题，就比较难定位了，因为可能的坑比较多。编译安装最好按官网的要求把需要的包都装上。

这种东西要看文档，你的ansys安装包如果全的话，应该有三个镜像文件，其中一个是文件名含有“DOCS”。
里面有很多PDF。PBM的手册在新版中没有单独的PDF。而在这里：ANSYS Fluent Advanced Add-On Modules
这里有你要找的东西。

@林之流风 能编译，要动那些影响生成文件的环境变量，和存放库文件的环境变量。他们不只是生成的位置，还是读取的位置。FOAM_USER 开头的应该都要变。没具体实施过，目前想象是这么些。

这些能操作成功的前提是 openfoam 的版本号得对的上。否则比较麻烦。还有一个思路可以考虑：交叉编译。即用A平台编译B平台可以运行的二进制文件。

无超级用户权限还是很麻烦的，要联系有权限的。不然，目前看依赖三个包。一装，每个包又各自依赖三个包，会让人崩溃的。

flex这个命令是缺了的。但是如果你没超级用户权限是无法按装很多依赖软件的，可以编译安装在指定文件夹，但会很麻烦。

修改OpenFOAM的安装位置这个就更麻烦了。可以参考这个：二楼的帖子

一楼的从三楼描述来看其实并没有修改安装位置，只是修改了源代码的位置，安装位置还是默认的。

很多资源只是在国内稀缺，国外很多免费的。
human body

把allrun里的 run application $applications删除。然后执行allrun后然后手动执行fsifoam试试。

本来想装上软件跑个看看，我装都装不上。你都能装上了，肯定能解决，或者找装的人解决

问题描述清楚的好处，除了“使人理解”（描述不清无法使人理解），还有“使更多的人理解”。有价值的问题值得转化成文本与世界讨论。所以下一步才是体现“描述清楚问题”价值的一步：问题翻译成英文，去cfd online找个合适的板块发帖

@Zhoujq 自助者天助，你那么详细认真的描述我才给以同样回答。
侠骨柔情，无须回报

不过给联系方式没啥

我的微博主页
没错，虽然这里面都是舞蹈，但这是我的微博，你看头像都一样

@bestucan 在 关于paraview中VOF的奇怪现象 中说：

alpah.water里的object -> alpha_object_flag

写错一个
alpah.water里的object -> alpha_object_flag.water

@Zhoujq 签名里的链接都是维基百科的，国内登不太稳定，我换了，现在应该没问题了。

换名字当然会影响计算，文件名和 object 名应该是一样的。这种修改方式只是为了放大不同找出来问题。如果你怀疑是文件名影响或者最想orig影响……咦！你的怎么是org

你把文件名加上标识符接着读啊，把怀疑的也读上。

alpha.water                -> alpha_file_flag.water
alpah.water里的object      -> alpha_object_flag
alpha.water.orig           -> alpha_file_flag_o.water.orig
alpah.water.orig里的object -> alpha_object_flag_o.water.orig。

这样无论显始什么都被你打上标签了，你都知道是从哪来的，paraview 能看到的东西也知道是什么东西了。如果 alpha.water 还阴魂不散就一个文件一个文件的把 0 文件夹移空，总会有那么一个，你移走它，就读不出 alpha.water 。它就是那个鬼

你这描述的太详细了
先来谈谈这个描述问题。为什么要描述这么详细，因为环节太多，每个都可能出错，所以要尽量描述可能出错的细节（比如alpha.water文件内容）而忽略不怎么可能会影响的细节（比如cpu型号、硬盘大小）。

当然这些的前提是无法缩小错误范围，而你已经缩小错误范围了。就是在“纯ubuntu”（ubunut单系统，对应的叫双系统或者虚拟机环境）

即使是在纯Ubuntu系统的电脑中进行.msh的网格转换，在不进行计算保存时间步，还是在paraview中显示单相。

所以这个问题实际上可以先限定在这一句话。也就是说前面的问题可能是这这个问题引起的，先解决这个小问题，如果大问题不存在了更好，如果还存在，才需要去考虑那些操作。

就跟调程序一样，先十分确保小部件小程序正常工作。才去找它们联合工作的错误，这样可以忽略小部件错误的可能，把排错精力放在联合上。

然后是错误分析。如果勾选skip zero是无法看见0时间步的alpha.water的，你点了一下，不勾选它，所以能看见alpha.water的。说明软件能理解你的操作意图（显示alpha.water）。

而结果不如你所想，再往下捋就是：

1. 显始的那个文件就是你想要显始的文件么？
2. 你想要显始的文件是你预想的内容么？

1: 修改你认为显示的那个 alpah.water 的显始名称，确定 paraview 里看到的 alpha.water 名称跟着变了，确定就是这个文件的。

操作：为了修改在paraview里的名称不能直接修改文件的名称，而要打开 alpha.water 。把object的名称加上自定义字符比如flag 。（以后这种过程可能会出错，因为有些名字有修改规则，改了就不能正常运行，但这里可以这样改）

    object      alpha_flag.water;

这个过程看起来很蠢很多余，但这确实是一个环节。错误往往出现在忽视的地方。文件多了很容易出现这个问题。

2.已经找到了就是这个文件，paraview 显始的就是它，但不如我所想，那么是paraview读错了？还是读对了但实际上只是文件的内容不如你所想？到这基本就解决了。

如果你眼睁睁看着alpha.water里的场里的值一大片1，而paraview里就是没有显示。你可以找paraview的开发团队报bug了

我不是大神，我只是耗的时间多，熟能生巧，你投入时间长了，你也可以

@ibelief 网盘丫

见过生物体管道仿真的，很少见这种带弹性的。可能生物体内约束关系比较复杂。
displacement
pressure

你这个，new是你自己起得函数名用来代替 phaseChangTwoPhaseMixture 这个函数的？

替代就直接替代，不能改名字啊。如果我没猜错，原始的调用代码是

phaseChangeTwoPhaseMixture(mesh);

这个不是函数调用，是类变量的隐性初始化（隐性初始化也用括号，容易和函数混），这是个声明类变量并赋值的过程，隐性初始化会调用构造函数并传递参数。那个和类名称同名的函数叫构造函数，必须和类名称一模一样。在声明该类变量时会自动调用构造函数。

也就是说，看着你在调用new函数，但是new函数是属于某个类对象的，而你没有声明这样的类对象，所以就找不到new函数。

改就直接改构造函数，不要换名字。要么复刻一个一模一样，只有类名不一样的类。

建议先看看C++，不用多么精，至少都浏览一遍有印象。你调试的这些时间用来浏览早就浏览完了…

描述不详细。无法判断。
这台电脑什么系统、文件从哪来的、来文件的那个什么系统、操作步骤是什么（1234）、文件细节展示（ls -la）。

可能是paraview、可能是文件、可能是系统差别、可能是……

都tutorials了，又没啥机密，使劲往详细了写啊，还是要一直来回猜，不过能看出来你已经尽量往详细了写了

看我签名里有讲怎么写的。

一般服务器的显卡性能都比较弱，一是因为服务器很少用显卡，二（主要的）是因为服务器上的linux系统对独立显卡支持很差，很不稳定，这违反了服务器的初衷（稳定）。服务器上的图形界面对显卡的应用也很受限， xorg 已经不行了，祖传代码成屎山了，wayland作为接班的还没发展起来。所以同样的卡在服务器上也会性能很差。不如在个人计算机上。用服务器加速渲染也一般是作为远程节点，而不是直接显示（比如我最后提到的那种）

想看是什么卡了，可以这样看（system monitor 也可以，不过不能一下看完）：

用conky看，随时都可以看。就像桌面上的背景，无法操作，只显示。左边的CPU，内存，磁盘读写。谁卡就会飙到99%，看谁敢卡，谁卡kill谁

这是我配置过的，默认显示的信息比较少。配置文件在这里conky 配置文件

把这些文件放到家目录下的.conky文件夹，然后从终端运行这些文件里的conkym脚本就行了。（可以把conkym放到自己的路径里，然后做个desktop运行它，然后设定开机运行这个desktop）

注：最近发现一个用服务器加速渲染的比较简易的方法，仅限于 paraview。可以开在服务器端开多线程 pvserver。然后再在服务器上开 paraview，远程连接（其实是本地）pvserver，多线程渲染。直接多线程开paraview好像不行，会开出来一大堆单线程而不是一个多线程。这种可以用于渲染动画，需要很长时间，放那不管它就行。只是看看，远程开了 pvserver，在哪用同版本的 paraview 都可以看。

@小考拉 为什么有三个文件?

我将phaseChangeTwoPhaseMixture.C phaseChangeTwoPhaseMixture.H 和 newphaseChangeTwoPhaseMixture.C文件对应部分都改成了

定义一个类应该两个文件就行了。看样子类型也对上了。

@cccrrryyy 惭愧惭愧，您的理解比我深入多了，

我再猜测一下，可能是求解动量守恒方程的时候（夹在物质守恒和能量守恒中间？）要用密度。我记得求解器里的压力是除过密度的压力。这个时候开头算不算密度应该会有影响？

好了，肚里没货了，再和您讨论要明天充充电了

@cccrrryyy 第二个我知道，求解密度方程是为了求解质量守恒方程，连续性残差是密度残差的反映，具体什么关系没研究过，当做相等也没什么，顶多有个常系数关系。

质量守恒方程所运输的物理量就是密度，所以一直在围着密度转，所以时间步开头就算密度，然后通过迭代减小残差。

• 可能错误1 new函数没有定义，如果是定义了，那就是定义的形参和调用时的实参不是一个数据类型。

因为c++默认的函数重载，参数不同就是不同的函数。

• 可能错误2 新的class没有包含进编译时的链接库，就是你在尝试解决的方向。
  用这个命令进入

cd $FOAM_USER_LIBBIN

看看新的class的so文件在不在。如果在，那基本就是可能错误1了。

• 另：编译class能过这点没有太大的说明意义。因为是个库文件，里面都是类的定义和成员函数，除了成员函数之间互相调用，不会有其他调用，所以很难出错。如果细看，也能看见myFoam.o是已经编译出来了的，错误是出在链接各.o文件。这时候就要检查调用了，就出错了。

总感觉缺一个点。
得证明大密度瓜皮可以等效替代低密度瓜瓤。
因为体积、质量都一样，然后质心也在同一点上。

@七辆战车 船就是这样的，如果把船看成整体，平均密度比水小，所以才浮着，但是船体是铁，密度比水大。

实现变密度不倒翁不一定非得是实心的。事实上玩具不倒翁也不是实心的。其实真的不倒翁，放大好多倍就是“上一半密度比水小，下一半密度比水大”。

就像一个西瓜被掏空了一半，它就是个不倒翁西瓜，而且上下密度不等。西瓜皮加空气的质量除以半个西瓜的体积就是上半部的密度（可以调整瓜皮厚度来调整上半部密度）。下半部就是西瓜的密度。

@七辆战车 你好，我也是年轻的研究生

照你这么一描述，特别像一个圆的木船嘛

可以搜搜船舶类的仿真，给船底整厚点，然后加个盖…

从这入手，然后在论坛搜 两相流，多相流。然后看看我签名里的链接 https://www.cfd-china.com/topic/1152/多相流与矩方法的数学模型/1?_=1592748369580

https://www.jianguoyun.com/p/DS1BNA8QsLv2BRjPnqUD

或者把漂浮体用可压缩变形的？压缩了，密度就上去了。

@七辆战车 中文搜不到？不用怀疑了！你就是国内第一人！担负起中国在此领域与世界沟通的桥梁作用！

把要搜的东西翻译成英文，增删替换关键词，多搜几次，一边搜一边换更合适的关键词。https://www.hatarilabs.com/ih-en/floating-object-stability-modeling-with-openfoam

@东岳
肺动脉、肺动脉2、脑动脉
里面有分享许可证，好像是可以随便用的。

这个没找到许可证，但是好像是免费分享模型的网站
血管

李老师也不清楚，那更有研究价值了
倾情推荐《纳米颗粒两相流体动力学》，于老师和林老师的大作

于老师越来越忙了，现在见上一面都困难

OpenFOAM的网格是纯文本的，已经很难藏信息了。但是想藏还是行得通的，只是不知到他们会不会这么办，比如对可计算分布的网格，每个节点偏移本应该的位置一点点，这点偏移量可以编码。

正交的还是比较有特点的，较可能看出。多边形的都一样乱糟糟……

限制华为的不是安卓，是google全家桶。不只是国人对QQ、微信的依赖。国外非常依赖个人的管理和规划，以提升个人效率（国内达不到也跟啥优劣没关，纯粹的穷太久，顾不上）。个人的管理和规划又很依赖邮件、日历、待办、协同工具。这些个人效率套装能三端同步的就苹果、微软、谷歌。谷歌是用起来最“便宜”的、功能最全的（苹果短腿office；微软短腿开发环境；苹果、微软都短腿开源、搜索；而且单论邮箱gmail是独一档，对比其他的有领先一代的优势）、最不挑平台的。这还没提google play、google book、google map，手机电脑的密码可以自动同步填写，这些细节功能。所以某种情况下没有google的安卓对大部分外国人就像砖一样。华为还是可以用开源版本的安卓的，不过已经没啥意义了。

所以应该还没有开源被限制的先例。

即使限制，怎么限制啊，学术出版商使劲限制论文盗版下载不也没啥用么。牛叉的人翻手为云覆手为雨，难免会有能控制一切的错觉（比如控制大自然）。脑子里东西、思想的传播是多大的权威都无法限制的。属于可以被毁灭但不会被打败的存在。

退一万步，OpenFOAM只是最火的开源CFD，不是唯一的。就算国外的都用不成。国内的也不是开发不出来。这不是树底下长不成树麻。先驱的虹吸效应，后来的没成长的养料。

condenstaion 凝结（用“冷”，觉得限定凝结的原因）
coagulation 聚并（硬把它称为聚并，不用一个与“凝结”重合的字，提醒自己不混，用“碰”，觉得偏固体且偏两个碰撞对象）
nucleation 成核（用结晶感觉更像是固体）

coagulation ?=? aggregation

aggregation<---相对--->breakage

限于目前的视野范围，着这样一一对应他们。感觉中文论文提到的词太乱了，凝结、冷凝、凝聚、聚沉、碰并、结晶。然后方程也揉到一起，盹儿~的一下就变得特别简单……去某术语委员会查，一个词还按领域不同有不同的译名。交叉领域真乱套了……

fox的大作里没用coagulation，这个词只在引用里有。通篇全用的aggregation。但是aggregation又是只针对固体颗粒用，对液滴和气泡用coalescence。

对着aggregation查.查着查着又查出来agglomeration。这回彻底搞不清了。

留个坑，以后搞清了来填。虽然整英文的更高大上，还是希望中文也能更清晰。

paraview 可以直接读 fluent 的 case 。读的时候有两个选项，其中一个可以读对。

问题是你没说明你的代码适用的 OpenFOAM 版本

虽然我们用 emacs 的和用 vim 的向来不和，但是为了合伙掐死用 vs code 的还是可以站在一起的๑乛◡乛๑ 

• emacs插件：疑似OpenFOAM的开发者之一 Henry 的配置文件。
  配置里还有很多 OpenFOAM 自带的命令，这下破案了，那些命令根本不是记住然后手动在bash里用的，是写到脚本里，绑到快捷键上调用的。

• vim插件：vim三种插件。

躲在角落里擦一擦自己的机械硬盘，还能用emacs再续两年……

你给的命名只是几何或网格的，并不是严格的边界条件类型名称。在转成求解器可读的格式的时候（ openfoam 的转换或 fluent 的 read mesh ），它需要一个边界条件类型，这就得靠软件（fluent3DMeshToFoam 和 read mesh）自己猜了。你用inlet和outlet，它还能猜出来，猜不出来的就默认wall了。这在fluent里也是的，导入网格之后除了非常明显表明边界条件类型的part名，其他的要手动设边界条件（就像你上面的命名过程一样）。

那啥，你把"."去了就行了눈_눈

@yuhxFoamer 这哪里是编译不过去，这都没开始编译。wmake就是腾挪源文件的，g++给的才是编译错误。libso 只是 wmake 的运行参数，这个错误只是bash脚本里写的输出。 不过你的命令提示符真是绝了，才看件那是个命令

@yuhxFoamer 那你的炉灶起的太干净了……都是手动的，很容易出问题。
比如wmakeLnInclude的描述

# Description
#     Link all the source files in the <dir> directory into <dir>/lnInclude

应该是 link 到某个 lnInclud 的，wmake是个加壳的make。在foam文件夹下的wmake里有操作过程。wmake调用wmakeLnIncludeAll，wmakeLnIncludeAll调用wmakeLnInclude。硬调也不是不能调。问题可能在option。

植入代码这种事情，如果不是轻车熟路。有个笨方法，但很稳。先把自带的A复制成A1，除了名字什么都不动，基本不会出错。然后再趁计算机不注意，把A1里的东西悄咪咪的、慢慢的换成你的。期间常编译看效果。

@麦迪文 这个里面装X的部分不要看。看看如何描述清楚自己的问题：提问的智慧

看你那么爱用省略号，一看就是囧于表达的。大胆的说，没人吃你，
毕竟你看起来不是那么好吃

https://cfd.direct/openfoam/user-guide/v7-snappyhexmesh/#x27-1970005.4
OpenFOAM7里有，能搜到locationinmesh。可能是你没装snappyHexMesh所需的第三方？esi 版本的默认安装第三方很多，社区的很少。

我去，还能这么调……不是应该找出液相消失的原因么 (:3」∠)
不会是没定义好出口边界吧，进的气体把液体从入口逼出了。

另外你对问题的描述真的挺详细的，再多一句我就能复现你的试验了o(*≧▽≦)ツ┏━┓