High-order consistent SPH with the pressure projection method in 2-D and 3-D

被引:18
|
作者
Nasar, A. M. A. [1 ]
Fourtakas, G. [1 ]
Lind, S. J. [1 ]
King, J. R. C. [1 ]
Rogers, B. D. [1 ]
Stansby, P. K. [1 ]
机构
[1] Univ Manchester, Sch Engn, Manchester, Lancs, England
基金
英国工程与自然科学研究理事会;
关键词
Smoothed particle hydrodynamics; High-order accurate consistency corrections; High-order SPH smoothing error; Projection methods for incompressible SPH; Stabilisation techniques for co-located; projection methods; SMOOTHED PARTICLE HYDRODYNAMICS; FINITE-ELEMENT-METHOD; BOUNDARY-CONDITIONS; INCOMPRESSIBLE SPH; FLOWS; DIFFUSION; ERROR; HYPERVISCOSITY; IMPLEMENTATION; STABILITY;
D O I
10.1016/j.jcp.2021.110563
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Mesh-free methods such a smoothed particle hydrodynamics (SPH) have advantages over mesh-based methods for flow in complex domains but attaining consistent high-order accurate solutions with the conventional form of SPH has yet to be resolved. The high order smoothing error dominates the SPH error until increasingly fine resolutions cause the low-order discretisation error to dominate; this is related to discretising a volume integral into a summation. In this paper, a high-order consistency correction based on modified SPH (MSPH) and the modified finite particle method (FPM) is proposed for improving the order of the limiting discretisation error. The new technique is an arbitrarily high-order extension of these schemes where the complexity of the consistency correction and the required computations are reduced by using simplified versions of the smoothing kernel derivatives. Tested in Eulerian form, the proposed high-order consistent SPH technique (HOCSPH) is combined with new high-order SPH kernel functions, designed to improve the order of the SPH smoothing error, and the resulting hybrid technique is shown to converge according to the smoothing error initially before converging according to the HOCPSH error once the latter becomes dominant. The initial high-order convergence lowers the computational effort required to achieve higher accuracy with hybrid HOCSPH in comparison to HOCSPH with second-order smoothing accuracy kernel functions. However, for highly irregular distributions it is found that the use of kernels with second-order smoothing accuracy provides more consistent convergence properties. A number of flows are simulated in 2-D and 3-D using the new HOCSPH technique in combination with the pressure projection method, and the results show that the method is accurate and able to model highly complex flow patterns. Some issues with stability of the projection method related to pressure-velocity collocation are identified, and several remedies are proposed. While effective for the test cases herein, these remedies are new in this context and require further attention for generalisation in future studies. (C) 2021 Published by Elsevier Inc.
引用
收藏
页数:36
相关论文
共 50 条
  • [1] Efficient optimization-oriented design methodology of high-order 3-D filters using 2-D and 3-D electromagnetic simulators
    Hinojosa, J.
    Quesada-Pereira, F. D.
    Bozzi, M.
    Alvarez-Melcon, A.
    [J]. INTERNATIONAL JOURNAL OF CIRCUIT THEORY AND APPLICATIONS, 2015, 43 (10) : 1431 - 1445
  • [2] A HIGH-ORDER FLUX RECONSTRUCTION METHOD FOR 2-D VORTICITY TRANSPORT
    Gharakhani, Adrin
    [J]. PROCEEDINGS OF ASME 2021 FLUIDS ENGINEERING DIVISION SUMMER MEETING (FEDSM2021), VOL 1, 2021,
  • [3] A Galerkin method for the simulation of the transient 2-D/2-D and 3-D/3-D linear Boltzmann equation
    Gobbert, Matthias K.
    Webster, Samuel G.
    Cale, Timothy S.
    [J]. JOURNAL OF SCIENTIFIC COMPUTING, 2007, 30 (02) : 237 - 273
  • [4] A Galerkin Method for the Simulation of the Transient 2-D/2-D and 3-D/3-D Linear Boltzmann Equation
    Matthias K. Gobbert
    Samuel G. Webster
    Timothy S. Cale
    [J]. Journal of Scientific Computing, 2007, 30 : 237 - 273
  • [5] 3-D/2-D registration by integrating 2-D information in 3-D
    Tomazevic, D
    Likar, B
    Pernus, F
    [J]. IEEE TRANSACTIONS ON MEDICAL IMAGING, 2006, 25 (01) : 17 - 27
  • [6] A High-Order 2-D CPITD Method for Electrically Large Waveguide Analysis
    Kang, Zhen
    Ma, Xikui
    Liu, Qi
    [J]. IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS, 2016, 26 (02) : 83 - 85
  • [7] A unified framework for consistent 2-D/3-D foreground object detection
    Landabaso, Jose-Luis
    Pardas, Montse
    [J]. IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, 2008, 18 (08) : 1040 - 1051
  • [8] Optimal projection of 2-D displacements for 3-D translational motion estimation
    Garcia, C
    Tziritas, G
    [J]. IMAGE AND VISION COMPUTING, 2002, 20 (11) : 793 - 804
  • [9] 3-D simulations of freak waves based on high-order spectral method
    赵西增
    孙昭晨
    梁书秀
    [J]. Journal of Harbin Institute of Technology., 2010, 17 (02) - 291
  • [10] 3-D simulations of freak waves based on high-order spectral method
    赵西增
    孙昭晨
    梁书秀
    [J]. Journal of Harbin Institute of Technology(New series), 2010, (02) : 286 - 291