Please wait a minute...
浙江大学学报（理学版）
浙江大学学报（理学版）  2023, Vol. 50 Issue (1): 49-55    DOI: 10.3785/j.issn.1008-9497.2023.01.008
数学与计算机科学     
三维非牛顿流体充填过程的有限元-间断有限元数值模拟研究
高普阳()
长安大学 理学院，陕西 西安 710064
The numerical investigation of non-Newtonian fluid filling process via finite element and discontinuous Galerkin method
Puyang GAO()
School of Science，Chang'an University，Xi'an 710064，China
 全文: PDF(2056 KB)   HTML( 6 )
摘要：

针对三维非牛顿流体充填问题，建立了有限元-间断有限元耦合算法。对于两相Navier-Stokes方程，基于压力增量修正格式分三步求解，分别采用二次和一次拉格朗日插值多项式求解速度和压力，以确保计算过程稳定。采用守恒型水平集（level set）方法追踪运动界面，并依据间断有限元方法求解水平集和重新初始化方程。以三维圆球剪切流动及非牛顿流体三维平板型腔充填过程为例，并与已有文献的数值和实验结果进行比较，以验证数值算法的稳定性、准确性以及流体的质量守恒性。
关键词： 有限元;  间断有限元;  水平集;  非牛顿充填过程    
Abstract:

In this paper, we develop a coupled finite element and discontinuous Galerkin method in three dimension and study the non-Newtonian fluid filling process. To solve two phase Navier-Stokes equations, we employ the incremental pressure correction scheme to accomplish it in three steps. In order to guarantee the computational stability, we take the second order and first order interpolation polynomials for the velocity and pressure, respectively. In addition, the conservative Level Set method is employed to capture the moving interface. The discontinuous Galerkin method is used to solve the Level Set and its re-initialization equations. We take the three dimensional vortex shearing problem and the three dimensional non-Newtonian fluid filling process to verify the proposed approach, compare the result with the numerical results and existing experimental data to illustrate the stability, accuracy and the mass conservation property of the coupled scheme.
Key words: finite element    discontinuous Galerkin    level set    non-Newtonian filling process
收稿日期: 2021-12-13 出版日期: 2023-01-13
CLC:  O 242.1  
基金资助: 国家自然科学基金资助项目(11901051);陕西省自然科学基础研究计划青年项目(2020JQ-338);长安大学中央高校基本科研业务费项目(300102122107);陕西省科学技术协会青年人才托举计划项目(20220504)
作者简介: 高普阳（1991—），ORCID：https：//orcid.ord/0000-0001-8620-1783，男，博士，讲师，主要从事非牛顿流动问题的数值算法研究，E-mail： gaopuyang@chd.edu.cn.
服务　　
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章　　
高普阳

引用本文:

高普阳. 三维非牛顿流体充填过程的有限元-间断有限元数值模拟研究[J]. 浙江大学学报（理学版）, 2023, 50(1): 49-55.

Puyang GAO. The numerical investigation of non-Newtonian fluid filling process via finite element and discontinuous Galerkin method. Journal of Zhejiang University (Science Edition), 2023, 50(1): 49-55.

链接本文:

https://www.zjujournals.com/sci/CN/10.3785/j.issn.1008-9497.2023.01.008    或    https://www.zjujournals.com/sci/CN/Y2023/V50/I1/49

图1  算法流程
图2  初始状态下的圆球
图3  当t =3.0时不同网格的自由界面形态
图4  不同时刻的自由界面形态
图5  网格3上相对质量误差随时间的变化
图6  三维型腔示意
图7  计算网格3示意
图8  当t =1.9 s时不同网格的熔体前沿界面形态
图9  不同时刻的实验结果（左）、文献［9］的数值结果（中）和本文的数值结果（右）
图10  网格3熔体质量相对误差随时间的变化
图11  不同时刻前沿界面的三维视图
13 HE L P， LU G， CHEN D C， et al. Three-dimensional smoothed particle hydrodynamics simulation for injection molding flow of short fiber-reinforced polymer composites［J］. Modelling and Simulation in Materials Science and Engineering， 2017， 25（5）： 055007. doi:10.1088/1361-651X/aa6dc9
doi: 10.1088/1361-651X/aa6dc9
14 BORZENKO E I， HEGAJ E I. Three-dimensional simulation of a tank filling with a viscous fluid using the VOF method［J］. Journal of Siberian Federal University （Mathematics & Physics）， 2020， 13（6）： 670-677. DOI：10.17516/1997-1397-2020-13-6-670-677
doi: 10.17516/1997-1397-2020-13-6-670-677
15 GU Z H， WEN H L， YU C H， et al. Interface-preserving level set method for simulating dam-break flows［J］. Journal of Computational Physics， 2018， 374： 249-280. DOI：10.1016/j.jcp.2018.07.057
doi: 10.1016/j.jcp.2018.07.057
16 AGUIRRE A， CASTILLO E， CRUCHAGA M， et al. Stationary and time-dependent numerical approximation of the lid-driven cavity problem for power-law fluid flows at high Reynolds numbers using a stabilized finite element formulation of the VMS type［J］. Journal of Non-Newtonian Fluid Mechanics， 2018， 257： 22-43. DOI：10.1016/j.jnnfm.2018.03.014
doi: 10.1016/j.jnnfm.2018.03.014
17 GUERMOND J L， MINEV P， SHEN J. An overview of projection methods for incompressible flows［J］. Computer Methods in Applied Mechanics and Engineering， 2006， 195（44-47）： 6011-6045. DOI：10.1016/j.cma.2005.10.010
doi: 10.1016/j.cma.2005.10.010
18 GUERMOND J L， de LUNA M Q， THOMPSON T. An conservative anti-diffusion technique for the level set method［J］. Journal of Computational and Applied Mathematics， 2017， 321： 448-468. DOI：10.1016/j.cam.2017.02.016
doi: 10.1016/j.cam.2017.02.016
1 XU X Y， OUYANG J， YANG B X， et al. SPH simulations of three-dimensional non-Newtonian free surface flows［J］. Computer Methods in Applied Mechanics and Engineering， 2013， 256： 101-116. DOI：10.1016/j.cma.2012.12.017
doi: 10.1016/j.cma.2012.12.017
2 FIGUEIREDO R A， OISHI C M， CUMINATO J A， et al. Three-dimensional transient complex free surface flows： Numerical simulation of XPP fluid［J］. Journal of Non-Newtonian Fluid Mechanics， 2013， 195： 88-98. DOI：10.1016/j.jnnfm.2013.01.004
doi: 10.1016/j.jnnfm.2013.01.004
3 ZHANG C H， CHEN S G， SUN Q C， et al. Free-surface simulations of Newtonian and non-Newtonian fluids with the lattice Boltzmann method［J］. Acta Geologica Sinica（English Edition）， 2016， 90（3）： 999-1010. DOI：10.1111/1755-6724.12740
doi: 10.1111/1755-6724.12740
4 周文， 欧阳洁， 杨斌鑫， 等. 三维非等温非牛顿流体充模过程的建模与模拟［J］. 化工学报， 2011， 62（3）： 618-627. doi:10.3788/gzxb20114002.0199
ZHOU W， OUYANG J， YANG B X， et al. Modeling and simulation of 3D mold filling process for non-isothermal non-Newtonian fluid［J］. CIESC Journal， 2011， 62（3）： 618-627. doi:10.3788/gzxb20114002.0199
doi: 10.3788/gzxb20114002.0199
5 BAUM M， ANDERS D. A numerical simulation study of mold filling in the injection molding process［J］. Computer Methods in Materials Science， 2021， 21（1）： 25-34. DOI：10.1002/fld.166
doi: 10.1002/fld.166
19 LUO K， SHAO C X， YANG Y， et al. A mass conserving level set method for detailed numerical simulation of liquid atomization［J］. Journal of Computational Physics， 2015， 298： 495-519. DOI：10.1016/j.jcp.2015.06.009
doi: 10.1016/j.jcp.2015.06.009
6 LI Q， QU F C. A level set based immersed boundary method for simulation of non-isothermal viscoelastic melt filling process［J］. Chinese Journal of Chemical Engineering， 2021， 32： 119-133. DOI：10.1016/j.cjche.2020.09.057
doi: 10.1016/j.cjche.2020.09.057
7 XU X Y， TIAN L Y， PENG S， et al. Development of SPH for simulation of non-isothermal viscoelastic free surface flows with application to injection molding［J］. Applied Mathematical Modelling， 2022， 104： 782-805. DOI：10.1016/j.apm.2021.12.015
doi: 10.1016/j.apm.2021.12.015
8 TOMÉ M F， CASTELO A， NÓBREGA J M， et al. Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids［J］. Computers & Fluids， 2014， 90： 172-185. DOI：10.1016/j.compfluid.2013.11.015
doi: 10.1016/j.compfluid.2013.11.015
9 MUKRAS S M S， AL-MUFADI F A. Simulation of HDPE mold filling in the injection molding process with comparison to experiments［J］. Arabian Journal for Science and Engineering， 2016， 41（5）： 1847-1856. DOI：10.1007/s13369-015-1970-9
doi: 10.1007/s13369-015-1970-9
10 LIU Q S， OUYANG J， ZHOU W， et al. Numerical simulation of the sequential coinjection molding process based on level set method［J］. Polymer Engineering & Science， 2015， 55（8）： 1707-1719. DOI：10.1002/pen.24009
doi: 10.1002/pen.24009
11 LIU Q S， OUYANG J， LI W M， et al. Three-dimensional modelling and simulation of sequential co-injection moulding with application to a toothbrush handle［J］. The Canadian Journal of Chemical Engineering， 2016， 94（2）： 382-390. DOI：10.1002/cjce.22394
doi: 10.1002/cjce.22394
[1] 孙美玲,江山. 多尺度有限元法结合分层网格模拟二维奇异摄动的两端边界层问题[J]. 浙江大学学报（理学版）, 2022, 49(5): 564-569.