质量流量 — CAE术语解释
质量流量
老师,质量流量是CFD的基本量吧。怎么使用呢?
质量流量的理论基础
质量流量的定义和基本概念
"质量流量"就是"单位时间内流动的流体质量"对吧?和体积流量有什么区别呢?
理解完全正确。定义式是
密度变化意味着气体分析时,体积流量比质量流量更本质吗?CAE软件中通常如何输入?
完全正确。气体的密度因压力、温度大幅变动,因此用质量流量给定边界条件在物理上更稳定。Ansys Fluent和Siemens STAR-CCM+等主流CAE软件,可以直接将"Mass Flow Inlet"指定为进口边界条件。例如,在汽车发动机进气口分析中,可以根据气门升程输入瞬时质量流量曲线。相反,对于送风机出口等压力基本恒定的情况,通常设置"Pressure Outlet",将质量流量作为计算结果输出。
质量流量守恒的连续方程具体怎么写?CAE计算中如何使用?
微分形式的质量守恒方程(连续方程)为
质量流量的数值计算手法
离散化和求解器中的处理
我无法想象求解器内部如何处理"质量流量边界条件"。和指定速度有什么根本区别?
根本区别在于作为未知数处理的变量。速度边界条件将速度矢量固定,求解器调整压力和密度来确定质量流量。而质量流量边界条件则固定
离散化时,单元界面的质量流量怎么计算?简单平均吗?
不是简单算术平均。取决于对流项的离散化格式。例如,一阶迎风格式会用上游单元的值计算
求解器的"Pressure-Based"和"Density-Based"对质量流量的处理方式会改变吗?
改变很大。Pressure-Based求解器(适用于低速~亚音速流动)采用压力修正法(SIMPLE、PISO),连续方程以压力修正方程形式满足。质量流量是这个修正过程的中心变量。而Density-Based求解器(适用于超音速流动)耦合求解密度、动量、能量守恒方程。两者都严格处理质量守恒,但求解方法不同。例如,在Fluent中求解马赫数0.3以上的压缩流时,选择Density-Based求解器,质量流量边界条件会得到更直接的处理,冲击波分辨率可能会提高。
质量流量的实务应用
边界条件设置和结果验证
实际分析中,在进口设置质量流量时,用什么单位比较保险?SI单位kg/s感觉数值太小,容易输入错误。
实务中根据项目和行业惯例来选择。汽车行业常用g/s或kg/h。例如,发动机进气量通常在数十到数百kg/h的数量级。航空航天领域则常见lb/s(磅每秒)。重要的是在软件内保持单位系的一致性。Ansys Workbench中设置"单位制"后,输入栏会显示对应的单位,很方便。另外,数值过小或过大会导致求解器的数值缩放问题,例如将1e-5 kg/s改为10 mg/s输入,调整到数量级1前后,也是一个有用的技巧。
计算完成后,如何检查质量流量是否守恒?具体看哪里?
主流CAE软件都有专门的监视器功能或报告功能。Ansys Fluent中,进入"Reports"→"Fluxes",可以计算选定的边界面(进口、出口、壁面等)的质量流量总和。在定常计算的收敛解中,所有流入面的合计与流出面的合计应该几乎一致(例如差异在1e-6 kg/s以下),在非定常计算中,这个差应等于计算域内的质量变化率。同时,在求解器的收敛历史中添加"Net Mass Flow Rate"或"Mass Imbalance"监视器,实时观察它在反复计算中趋近于零,是标准的检查步骤。
在多个进口或出口的系统中,部分边界条件为质量流量,其他为压力,有什么设置顺序或注意事项吗?
设置顺序本身影响不大,但要避免"过约束"或"欠约束"的物理矛盾。典型的过约束是在所有开口设置质量流量,这样压力无法确定,求解器会发散。至少需要一处采用压力基准(Pressure Inlet/Outlet)。反之,若全部设置压力条件,虽然不会欠约束,但当实际驱动力(例如风机)无法再现时,可能无法获得正确的物理过程。例如,有风机的房间分析中,在风机面设置质量流量(或体积流量),在外部开口设置压力条件,是最符合实际的做法。
质量流量的软件对比
各软件中的实现和特点
Ansys Fluent和Siemens STAR-CCM+在质量流量的输入界面和功能上有差别吗?
基本概念相同,但界面和高级功能有差异。Fluent在边界条件对话框中直接输入"Mass Flow Rate",用"Direction Specification Method"定义流入方向。STAR-CCM+中选择"Mass Flow",同时指定数值和"Flow Direction"("Normal to Boundary"或"Components"等)。STAR-CCM+的优势是可以用"Field Function"灵活地定义随位置或时间变化的质量流量。例如,可以在GUI上直接用数学表达式定义随活塞运动变化的正弦波质量流量曲线。Fluent虽然也能通过UDF(用户定义函数)实现,但需要编程知识。
COMSOL Multiphysics这类多物理场软件中呢?"质量流量"设置会因物理场而改变吗?
在COMSOL中,不同物理场界面的流入条件名称有差异,但本质相同。在"CFD模块"或"微观流"中,可以在"进口"(Inlet)边界指定"质量流量"。在"非等温流"或"反应流"中也是如此。COMSOL的优势在于能直接将这些与"固体力学"、"电磁场"等耦合。例如,在MEMS微泵分析中,可以将压电驱动的膜片变形(固体力学)在流路内自动耦合为流体质量源项。在这种情况下,用户无需直接指定质量流量,而是由耦合物理场间接确定。
开源的OpenFOAM中,质量流量边界条件怎么设置?
在OpenFOAM中,边界条件在"0/"目录下各个场(U、p、T等)的文件中按补丁指定。直接指定质量流量的边界类型是"flowRateInletVelocity"。在这个设置中,用`massFlowRate`或`volumetricFlowRate`关键字写入数值,求解器会内部调整边界速度。例如,`inlet { type flowRateInletVelocity; massFlowRate 0.01; // kg/s value uniform (0 0 0); }` 这样的形式。相比商用软件,需要直接编辑配置文件,初期学习成本较高,但自定义灵活性极强,特殊的质量流量曲线可以通过编程自由实现。
质量流量故障排除
常见错误和对策
设置质量流量进口后计算时,求解器出现"reversed flow"警告。这是问题吗?怎么处理?
"逆流"警告本身在计算初期或流动复杂的区域很常见,不一定是错误。但在定常计算收敛后仍然出现,就有问题了。原因可能是:设置的质量流量相对于出口压力过大,或流路内有大型障碍物导致压力损失超预期等。对策包括:(1) 改进初始条件,例如在进口附近的单元设置合理的初始速度;(2) 修改求解器设置,在Fluent中将"Pressure Outlet"的"Backflow Direction Specification Method"从"Normal to Boundary"改为"From Neighboring Cell"等;(3) 作为非定常计算求解。根本上,需要重新检查边界条件设置是否物理合理。
计算发散,出现"mass flow rate at inlet is too large"这样的错误。除了检查数值,还有其他原因吗?
有的。质量流量数值本身正确,但其他设置可能是原因。主要原因和对策如下:
1. **网格质量**:进口面网格过粗或扭曲,会导致界面流量计算不稳定。改进面网格质量。
2. **密度设置**:材料密度输入错误,或使用理想气体时操作条件(压力、温度)不现实。例如,温度应为300K却输入3K,密度会变成100倍,实际质量流量会大幅增加。
3. **求解器格式**:高阶精度格式配合大质量流量可能发振。先用一阶迎风差分使其稳定,收敛后再提高格式精度,这种"格式松弛"策略有效。
4. **时间步长(非定常)**:时间步太大,单步流入的质量超过单元容量(Courant数过大)。减小时间步长。
在多条流路汇合的"集流管"分析中,在各进口设置细致的质量流量,但出口总质量流量与设置值不符。是网格相关性吗?
网格相关性是一个因素,但更直接的原因是"数值流动分配"。在汇合部的复杂涡旋和逆流下,设定的流量不一定会按预期分配到各支路。对策包括:
1. **细化汇合部网格**:特别是各流路碰撞角的分辨率要提高。
2. **工程化初始条件**:在所有进口设置相同的速度(与质量流量成正比)作为初始条件,稳定计算初期。
3. **调整求解器参数**:缩小压力-速度耦合求解器(SIMPLE、PISO)的松弛系数,使各反复中的流量调整更温和。
4. **改进模型建立**:如可能,将集流管上游的公共管道也包含在模型中,在那里设置单一质量流量,让下游分配由流动阻力(压力损失)自然决定。这种做法往往更符合现实,数值上也更稳定。
相关主题
细节
错误