January 2019, 24(1): 127-147. doi: 10.3934/dcdsb.2018111

A dimension splitting and characteristic projection method for three-dimensional incompressible flow

1. 

School of Mechanics and Engineering, Southwest Jiaotong University, Chengdu 610031, China

2. 

Department of Civil and Mechanical Engineering, University of Missouri-Kansas City, Kansas City, MO, 64110, USA

3. 

School of Mathematics and Statistics, Xi'an Jiaotong University, Xi'an 710049, China

4. 

School of Mechanical Engineering, Dongguan University of Technology, Dongguan 523000, China

5. 

Department of Mathematics and Statistics, Missouri University of Science and Technology, Rolla, MO, 64509, USA

* Corresponding author: Yuchuan Chu(ychuan.chu@hit.edu.cn)

Received  May 2017 Revised  September 2017 Published  March 2018

Fund Project: The first author is supported by the Fundamental Research Funds for the Central Universities of China, grant 2682015CX044

A dimension splitting and characteristic projection method is proposed for three-dimensional incompressible flow. First, the characteristics method is adopted to obtain temporal semi-discretization scheme. For the remaining Stokes equations we present a projection method to deal with the incompressibility constraint. In conclusion only independent linear elliptic equations need to be calculated at each step. Furthermore on account of splitting property of dimension splitting method, all the computations are carried out on two-dimensional manifolds, which greatly reduces the difficulty and the computational cost in the mesh generation. And a coarse-grained parallel algorithm can be also constructed, in which the two-dimensional manifold is considered as the computation unit.

Citation: Hao Chen, Kaitai Li, Yuchuan Chu, Zhiqiang Chen, Yiren Yang. A dimension splitting and characteristic projection method for three-dimensional incompressible flow. Discrete & Continuous Dynamical Systems - B, 2019, 24 (1) : 127-147. doi: 10.3934/dcdsb.2018111
References:
[1]

A. Allievi and R. Bermejo, Finite element modified of characteristics for the Navier-Stokes equations, Int. J. Numer. Meth. Fluids, 32 (2000), 439-463. doi: 10.1002/(SICI)1097-0363(20000229)32:4<439::AID-FLD946>3.0.CO;2-Y.

[2]

V. Babu and S. Korpela, Numerical solution of the incompressible three-dimensional Navier-Stokes equations, Comput. Fluids, 22 (1994), 675-691. doi: 10.1016/0045-7930(94)90009-4.

[3]

O. Botella and R. Peyret, Benchmark spectral results on the lid-driven cavity flow, Comput. Fluids, 27 (1998), 421-433. doi: 10.1016/S0045-7930(98)00002-4.

[4]

R. Bouffanais, M. O. Deville and E. Leriche, Large-eddy simulation of the flow in a lid-driven cubical cavity, Phys. Fluids, 19 (2007), 055108. doi: 10.1063/1.2723153.

[5]

J. ChanJ. A. Evans and W. Qiu, A dual Petrov-Galerkin finite element method for the convection-diffusion equation, Comput. Math. Appl., 68 (2014), 1513-1529. doi: 10.1016/j.camwa.2014.07.008.

[6]

H. ChenK. Li and S. Wang, A dimension split method for the incompressible Navier-Stokes equations in three dimensions, Int. J. Numer. Meth. Fluids, 73 (2013), 409-435. doi: 10.1002/fld.3803.

[7]

H. ChenJ. SuK. Li and S. Wang, A characteristic projection method for incompressible thermal flow, Numer. Heat Tr. B-Fund., 65 (2014), 554-590. doi: 10.1080/10407790.2013.836052.

[8]

Z. Chen, Characteristic mixed discontinuous finite element methods for advection-dominated diffusion problems, Comput. Meth. Appl. Mech. Eng., 191 (2002), 2509-2538. doi: 10.1016/S0045-7825(01)00411-X.

[9]

A. J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), 745-762. doi: 10.1090/S0025-5718-1968-0242392-2.

[10]

A. J. Chorin, On the convergence of discrete approximations to the Navier-Stokes equations, Math. Comput., 23 (1969), 341-353. doi: 10.1090/S0025-5718-1969-0242393-5.

[11]

J. Douglas Jr and T. F. Russell, Numerical method for convection-dominated diffusion problem based on combining the method of characteristics with finite element of finite difference procedures, SIAM J. Numer. Anal., 19 (1982), 871-885. doi: 10.1137/0719063.

[12]

C. J. FreitasR. L. StreetA. N. Findikakis and J. R. Koseff, Numerical simulation of three-dimensional flow in a cavity, Int. J. Numer. Meth. Fluids, 5 (1985), 561-575. doi: 10.1002/fld.1650050606.

[13]

C. J. Freitas and R. L. Street, Non-linear transient phenomena in a complex recirculating flow: A numerical investigation, Int. J. Numer. Meth. Fluids, 8 (1988), 769-802. doi: 10.1002/fld.1650080704.

[14]

U. GhiaK. N. Ghia and C. T. Shin, High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method, J. Comput. Phys., 48 (1982), 387-411. doi: 10.1016/0021-9991(82)90058-4.

[15]

J. L. GuermondP. Minev and J. Shen, An overview of projection methods for incompressible flows, Comput. Meth. Appl. Mech. Eng., 195 (2006), 6011-6045. doi: 10.1016/j.cma.2005.10.010.

[16]

M. Hermanns, Parallel programming in Fortran 95 using OpenMP, 2002. Available from: http://www.openmp.org/wp-content/uploads/F95_OpenMPv1_v2.pdf.

[17]

C. Johnson, Numerical Solution of Partial Differential Equations by the Finite Element Method, Cambridge University Press, Cambridge, 1987.

[18]

J. R. Koseff and R. L. Street, Visualization studies of a shear driven three-dimensional recirculating flow, J. Fluids Eng., 106 (1984), 21-27. doi: 10.1115/1.3242393.

[19]

J. R. Koseff and R. L. Street, The lid-driven cavity flow: A synthesis of qualitative and quantitative observations, J. Fluids Eng., 106 (1984), 390-398. doi: 10.1115/1.3243136.

[20]

O. A. Ladyzhenskaya, the Mathematical Theory of Viscous Incompressible Flow, Gordon and Breach Science Publishers, New York, 1969.

[21]

K. LiA. Huang and W. Zhang, A dimension split method for the 3-d compressible Navier-Stokes equations in turbomachine, Commun. Numer. Meth. Eng., 18 (2002), 1-14.

[22]

K. W. Morton, A. Priestley and E. Süli, Convergence Analysis of the Lagrange-Galerkin Method with Non-Exact Integration, Technical report, Oxford University Computing Laboratory. Rept. N86/14, Oxford, 1986.

[23]

A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer-Verlag, Berlin, 1994.

[24]

C. ShuX. D. Niu and Y. T. Chew, Taylor series expansion and least squares-based lattice boltzmann method: three-dimensional formulation and its applications, Int. J. Mod. Phys. C, 14 (2003), 925-944. doi: 10.1142/S0129183103005133.

[25]

E. Süli, Convergence and nonlinear stability of the Lagrange-Galerkin method for the Navier-Stokes equations, Numer. Math., 53 (1988), 459-483. doi: 10.1007/BF01396329.

[26]

R. Temam, Sur l'approximation de la solution des equations de Navier-Stokes par la ḿethode des fractionnarires Ⅱ, Arch. Rational Mech. Anal., 33 (1969), 377-385.

[27]

C. Wu, A general theory of three-dimensional flow in subsonic and supersonic turbomachines of axial-, radial-, and mixed-flow types, Tech. Notes. Nat. Adv. Comm. Aeronaut., 1952 (1952), ⅱ+93 pp.

[28]

P. X. Yu and Z. F. Tian, A compact streamfunction-velocity scheme on nonuniform grids for the 2D steady incompressible Navier-Stokes equations, Comput. Math. Appl., 66 (2013), 1192-1212. doi: 10.1016/j.camwa.2013.07.013.

[29]

O. C. Zienkiewicz, P. Nithiarasu and R. L. Taylor, the Finite Element Method for Fluid Dynamics, seventh ed., Elsevier/Butterworth Heinemann, Amsterdam, 2014.

show all references

References:
[1]

A. Allievi and R. Bermejo, Finite element modified of characteristics for the Navier-Stokes equations, Int. J. Numer. Meth. Fluids, 32 (2000), 439-463. doi: 10.1002/(SICI)1097-0363(20000229)32:4<439::AID-FLD946>3.0.CO;2-Y.

[2]

V. Babu and S. Korpela, Numerical solution of the incompressible three-dimensional Navier-Stokes equations, Comput. Fluids, 22 (1994), 675-691. doi: 10.1016/0045-7930(94)90009-4.

[3]

O. Botella and R. Peyret, Benchmark spectral results on the lid-driven cavity flow, Comput. Fluids, 27 (1998), 421-433. doi: 10.1016/S0045-7930(98)00002-4.

[4]

R. Bouffanais, M. O. Deville and E. Leriche, Large-eddy simulation of the flow in a lid-driven cubical cavity, Phys. Fluids, 19 (2007), 055108. doi: 10.1063/1.2723153.

[5]

J. ChanJ. A. Evans and W. Qiu, A dual Petrov-Galerkin finite element method for the convection-diffusion equation, Comput. Math. Appl., 68 (2014), 1513-1529. doi: 10.1016/j.camwa.2014.07.008.

[6]

H. ChenK. Li and S. Wang, A dimension split method for the incompressible Navier-Stokes equations in three dimensions, Int. J. Numer. Meth. Fluids, 73 (2013), 409-435. doi: 10.1002/fld.3803.

[7]

H. ChenJ. SuK. Li and S. Wang, A characteristic projection method for incompressible thermal flow, Numer. Heat Tr. B-Fund., 65 (2014), 554-590. doi: 10.1080/10407790.2013.836052.

[8]

Z. Chen, Characteristic mixed discontinuous finite element methods for advection-dominated diffusion problems, Comput. Meth. Appl. Mech. Eng., 191 (2002), 2509-2538. doi: 10.1016/S0045-7825(01)00411-X.

[9]

A. J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), 745-762. doi: 10.1090/S0025-5718-1968-0242392-2.

[10]

A. J. Chorin, On the convergence of discrete approximations to the Navier-Stokes equations, Math. Comput., 23 (1969), 341-353. doi: 10.1090/S0025-5718-1969-0242393-5.

[11]

J. Douglas Jr and T. F. Russell, Numerical method for convection-dominated diffusion problem based on combining the method of characteristics with finite element of finite difference procedures, SIAM J. Numer. Anal., 19 (1982), 871-885. doi: 10.1137/0719063.

[12]

C. J. FreitasR. L. StreetA. N. Findikakis and J. R. Koseff, Numerical simulation of three-dimensional flow in a cavity, Int. J. Numer. Meth. Fluids, 5 (1985), 561-575. doi: 10.1002/fld.1650050606.

[13]

C. J. Freitas and R. L. Street, Non-linear transient phenomena in a complex recirculating flow: A numerical investigation, Int. J. Numer. Meth. Fluids, 8 (1988), 769-802. doi: 10.1002/fld.1650080704.

[14]

U. GhiaK. N. Ghia and C. T. Shin, High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method, J. Comput. Phys., 48 (1982), 387-411. doi: 10.1016/0021-9991(82)90058-4.

[15]

J. L. GuermondP. Minev and J. Shen, An overview of projection methods for incompressible flows, Comput. Meth. Appl. Mech. Eng., 195 (2006), 6011-6045. doi: 10.1016/j.cma.2005.10.010.

[16]

M. Hermanns, Parallel programming in Fortran 95 using OpenMP, 2002. Available from: http://www.openmp.org/wp-content/uploads/F95_OpenMPv1_v2.pdf.

[17]

C. Johnson, Numerical Solution of Partial Differential Equations by the Finite Element Method, Cambridge University Press, Cambridge, 1987.

[18]

J. R. Koseff and R. L. Street, Visualization studies of a shear driven three-dimensional recirculating flow, J. Fluids Eng., 106 (1984), 21-27. doi: 10.1115/1.3242393.

[19]

J. R. Koseff and R. L. Street, The lid-driven cavity flow: A synthesis of qualitative and quantitative observations, J. Fluids Eng., 106 (1984), 390-398. doi: 10.1115/1.3243136.

[20]

O. A. Ladyzhenskaya, the Mathematical Theory of Viscous Incompressible Flow, Gordon and Breach Science Publishers, New York, 1969.

[21]

K. LiA. Huang and W. Zhang, A dimension split method for the 3-d compressible Navier-Stokes equations in turbomachine, Commun. Numer. Meth. Eng., 18 (2002), 1-14.

[22]

K. W. Morton, A. Priestley and E. Süli, Convergence Analysis of the Lagrange-Galerkin Method with Non-Exact Integration, Technical report, Oxford University Computing Laboratory. Rept. N86/14, Oxford, 1986.

[23]

A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, Springer-Verlag, Berlin, 1994.

[24]

C. ShuX. D. Niu and Y. T. Chew, Taylor series expansion and least squares-based lattice boltzmann method: three-dimensional formulation and its applications, Int. J. Mod. Phys. C, 14 (2003), 925-944. doi: 10.1142/S0129183103005133.

[25]

E. Süli, Convergence and nonlinear stability of the Lagrange-Galerkin method for the Navier-Stokes equations, Numer. Math., 53 (1988), 459-483. doi: 10.1007/BF01396329.

[26]

R. Temam, Sur l'approximation de la solution des equations de Navier-Stokes par la ḿethode des fractionnarires Ⅱ, Arch. Rational Mech. Anal., 33 (1969), 377-385.

[27]

C. Wu, A general theory of three-dimensional flow in subsonic and supersonic turbomachines of axial-, radial-, and mixed-flow types, Tech. Notes. Nat. Adv. Comm. Aeronaut., 1952 (1952), ⅱ+93 pp.

[28]

P. X. Yu and Z. F. Tian, A compact streamfunction-velocity scheme on nonuniform grids for the 2D steady incompressible Navier-Stokes equations, Comput. Math. Appl., 66 (2013), 1192-1212. doi: 10.1016/j.camwa.2013.07.013.

[29]

O. C. Zienkiewicz, P. Nithiarasu and R. L. Taylor, the Finite Element Method for Fluid Dynamics, seventh ed., Elsevier/Butterworth Heinemann, Amsterdam, 2014.

Figure 1.  Splitting the flow domain $\Omega$
Figure 2.  Grid structure of two-dimensional manifold $D$
Figure 3.  Sketch of three-dimensional lid-driven cavity flow
Figure 4.  Velocity profiles on middle plane z = 0.5 for Re = 100
Figure 5.  Velocity profiles on middle plane z = 0.5 for Re = 400
Figure 6.  Velocity profiles on middle plane z = 0.5 for Re = 1000
Figure 7.  Streamline profile for various Reynolds numbers: Re = 100(A, B, C); Re = 400(D, E, F); x = 0.5(A, D); z = 0.5(B, E); y = 0.5(C, F)
Figure 8.  Streamline profile for various Reynolds numbers: Re = 1000(A, B, C); Re = 2000(D, E, F); x = 0.5(A, D); z = 0.5(B, E); y = 0.5(C, F)
Figure 9.  Three dimensional streamline for different Reynolds numbers: Re = 100(A, B, C); Re = 400(D, E, F)
Figure 10.  Three dimensional streamline for different Reynolds numbers: Re = 1000(A, B, C); Re = 2000(D, E, F)
Table 1.  Error of numerical solution with different mesh sizes
$\frac{1}{h}$$\|\vec u-\vec u_h\|_{L^2}$$\alpha$$\|p-p_h\|_{L^2}$$\alpha$$\kappa_{div}$
$4$1.149E-002-2.931E-001-4.741E-002
$8$3.513E-0031.7101.394E-0011.0726.953E-003
$16$8.765E-0041.8566.172E-0021.1242.304E-003
$32$1.927E-0041.9661.961E-0021.3013.826E-004
$\frac{1}{h}$$\|\vec u-\vec u_h\|_{L^2}$$\alpha$$\|p-p_h\|_{L^2}$$\alpha$$\kappa_{div}$
$4$1.149E-002-2.931E-001-4.741E-002
$8$3.513E-0031.7101.394E-0011.0726.953E-003
$16$8.765E-0041.8566.172E-0021.1242.304E-003
$32$1.927E-0041.9661.961E-0021.3013.826E-004
Table 2.  Convergence rate with different mesh sizes
$\frac{1}{h}$DSM-CDSM-D
$U_{L^2}$ rate$P_{L^2}$ rateCPU(s)$U_{L^2}$ rate$P_{L^2}$ rateCPU(s)
4--44.5--43.7
81.7101.072138.71.420.876162.4
161.8561.124206.31.491.075383.2
321.9661.301957.61.530.9711996.3
$\frac{1}{h}$DSM-CDSM-D
$U_{L^2}$ rate$P_{L^2}$ rateCPU(s)$U_{L^2}$ rate$P_{L^2}$ rateCPU(s)
4--44.5--43.7
81.7101.072138.71.420.876162.4
161.8561.124206.31.491.075383.2
321.9661.301957.61.530.9711996.3
Table 3.  Parallel performance of DSM-C at $1/h = 8, 16$
$p$ $1/h=8$ $1/h=16$
$T_p$$S_{p}$$E_{p}$ $T_p$$S_{p}$$E_{p}$
152.35--451.69--
237.931.380.69303.141.490.75
421.632.420.61170.442.650.66
616.723.130.52123.073.670.61
815.173.450.43102.424.410.55
1015.313.420.3491.814.920.49
1222.762.300.1994.104.800.40
$p$ $1/h=8$ $1/h=16$
$T_p$$S_{p}$$E_{p}$ $T_p$$S_{p}$$E_{p}$
152.35--451.69--
237.931.380.69303.141.490.75
421.632.420.61170.442.650.66
616.723.130.52123.073.670.61
815.173.450.43102.424.410.55
1015.313.420.3491.814.920.49
1222.762.300.1994.104.800.40
Table 4.  Parallel performance of DSM-C at $1/h = 32, 64$
$p$$1/h=32$$1/h=64$
$T_p$$S_{p}$$E_{p}$$T_P$$S_{p}$$E_{p}$
13847.32--32861.04--
22171.171.770.8917077.321.920.96
41183.793.250.819225.063.560.89
6875.914.390.736744.554.870.81
8717.345.360.675402.696.080.76
10648.785.930.594627.377.100.71
12616.416.240.524147.197.920.66
$p$$1/h=32$$1/h=64$
$T_p$$S_{p}$$E_{p}$$T_P$$S_{p}$$E_{p}$
13847.32--32861.04--
22171.171.770.8917077.321.920.96
41183.793.250.819225.063.560.89
6875.914.390.736744.554.870.81
8717.345.360.675402.696.080.76
10648.785.930.594627.377.100.71
12616.416.240.524147.197.920.66
[1]

Madalina Petcu, Roger Temam, Djoko Wirosoetisno. Averaging method applied to the three-dimensional primitive equations. Discrete & Continuous Dynamical Systems - A, 2016, 36 (10) : 5681-5707. doi: 10.3934/dcds.2016049

[2]

Caterina Calgaro, Meriem Ezzoug, Ezzeddine Zahrouni. Stability and convergence of an hybrid finite volume-finite element method for a multiphasic incompressible fluid model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 429-448. doi: 10.3934/cpaa.2018024

[3]

Cornel M. Murea, H. G. E. Hentschel. A finite element method for growth in biological development. Mathematical Biosciences & Engineering, 2007, 4 (2) : 339-353. doi: 10.3934/mbe.2007.4.339

[4]

Martin Burger, José A. Carrillo, Marie-Therese Wolfram. A mixed finite element method for nonlinear diffusion equations. Kinetic & Related Models, 2010, 3 (1) : 59-83. doi: 10.3934/krm.2010.3.59

[5]

Leonardi Filippo. A projection method for the computation of admissible measure valued solutions of the incompressible Euler equations. Discrete & Continuous Dynamical Systems - S, 2018, 11 (5) : 941-961. doi: 10.3934/dcdss.2018056

[6]

Arun K. Kulshreshth, Andreas Alpers, Gabor T. Herman, Erik Knudsen, Lajos Rodek, Henning F. Poulsen. A greedy method for reconstructing polycrystals from three-dimensional X-ray diffraction data. Inverse Problems & Imaging, 2009, 3 (1) : 69-85. doi: 10.3934/ipi.2009.3.69

[7]

Corinna Burkard, Roland Potthast. A time-domain probe method for three-dimensional rough surface reconstructions. Inverse Problems & Imaging, 2009, 3 (2) : 259-274. doi: 10.3934/ipi.2009.3.259

[8]

Masaru Ikehata, Mishio Kawashita. An inverse problem for a three-dimensional heat equation in thermal imaging and the enclosure method. Inverse Problems & Imaging, 2014, 8 (4) : 1073-1116. doi: 10.3934/ipi.2014.8.1073

[9]

Wangtao Lu, Shingyu Leung, Jianliang Qian. An improved fast local level set method for three-dimensional inverse gravimetry. Inverse Problems & Imaging, 2015, 9 (2) : 479-509. doi: 10.3934/ipi.2015.9.479

[10]

Binjie Li, Xiaoping Xie, Shiquan Zhang. New convergence analysis for assumed stress hybrid quadrilateral finite element method. Discrete & Continuous Dynamical Systems - B, 2017, 22 (7) : 2831-2856. doi: 10.3934/dcdsb.2017153

[11]

Kun Wang, Yinnian He, Yueqiang Shang. Fully discrete finite element method for the viscoelastic fluid motion equations. Discrete & Continuous Dynamical Systems - B, 2010, 13 (3) : 665-684. doi: 10.3934/dcdsb.2010.13.665

[12]

Junjiang Lai, Jianguo Huang. A finite element method for vibration analysis of elastic plate-plate structures. Discrete & Continuous Dynamical Systems - B, 2009, 11 (2) : 387-419. doi: 10.3934/dcdsb.2009.11.387

[13]

So-Hsiang Chou. An immersed linear finite element method with interface flux capturing recovery. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2343-2357. doi: 10.3934/dcdsb.2012.17.2343

[14]

Donald L. Brown, Vasilena Taralova. A multiscale finite element method for Neumann problems in porous microstructures. Discrete & Continuous Dynamical Systems - S, 2016, 9 (5) : 1299-1326. doi: 10.3934/dcdss.2016052

[15]

Qingping Deng. A nonoverlapping domain decomposition method for nonconforming finite element problems. Communications on Pure & Applied Analysis, 2003, 2 (3) : 297-310. doi: 10.3934/cpaa.2003.2.297

[16]

Runchang Lin. A robust finite element method for singularly perturbed convection-diffusion problems. Conference Publications, 2009, 2009 (Special) : 496-505. doi: 10.3934/proc.2009.2009.496

[17]

Christos V. Nikolopoulos, Georgios E. Zouraris. Numerical solution of a non-local elliptic problem modeling a thermistor with a finite element and a finite volume method. Conference Publications, 2007, 2007 (Special) : 768-778. doi: 10.3934/proc.2007.2007.768

[18]

Liejune Shiau, Roland Glowinski. Operator splitting method for friction constrained dynamical systems. Conference Publications, 2005, 2005 (Special) : 806-815. doi: 10.3934/proc.2005.2005.806

[19]

Alexander Zlotnik, Ilya Zlotnik. Finite element method with discrete transparent boundary conditions for the time-dependent 1D Schrödinger equation. Kinetic & Related Models, 2012, 5 (3) : 639-667. doi: 10.3934/krm.2012.5.639

[20]

Zhenlin Guo, Ping Lin, Guangrong Ji, Yangfan Wang. Retinal vessel segmentation using a finite element based binary level set method. Inverse Problems & Imaging, 2014, 8 (2) : 459-473. doi: 10.3934/ipi.2014.8.459

2017 Impact Factor: 0.972

Article outline

Figures and Tables

[Back to Top]