scholarly journals Second-Order Symmetric Smoothed Particle Hydrodynamics Method for Transient Heat Conduction Problems with Initial Discontinuity

Processes ◽  
2018 ◽  
Vol 6 (11) ◽  
pp. 215 ◽  
Author(s):  
Zhanjie Song ◽  
Yaxuan Xing ◽  
Qingzhi Hou ◽  
Wenhuan Lu

To eliminate the numerical oscillations appearing in the first-order symmetric smoothed particle hydrodynamics (FO-SSPH) method for simulating transient heat conduction problems with discontinuous initial distribution, this paper presents a second-order symmetric smoothed particle hydrodynamics (SO-SSPH) method. Numerical properties of both SO-SSPH and FO-SSPH are analyzed, including truncation error, numerical accuracy, convergence rate, and stability. Experimental results show that for transient heat conduction with initial smooth distribution, both FO-SSPH and SO-SSPH can achieve second order convergence rate, which is consistent with the theoretical analysis. However, for one- and two-dimensional conduction with initial discontinuity, the FO-SSPH method suffers from serious unphysical oscillations, which do not disappear over time, and hence it only achieves a first-order convergence rate; while the present SO-SSPH method can avoid unphysical oscillations and has second-order convergence rate. Therefore, the SO-SSPH method is a feasible tool for solving transient heat conduction problems with both smooth and discontinuous distributions, and it is easy to be extended to high dimensional cases.

Author(s):  
Samir Hassan Sadek ◽  
Mehmet Yildiz

This work presents the development of both weakly compressible and incompressible Smoothed Particle Hydrodynamics (SPH) models for simulating two-dimensional transient viscoelastic free surface flow which has extensive applications in polymer processing industries. As an illustration with industrial significance, we have chosen to model the extrudate swell of a second-order polymeric fluid. The extrudate or die swell is a phenomenon that takes place during the extrusion of polymeric fluids. When a polymeric fluid is forced through a die to give a polymer its desired shape, due to its viscoelastic non-Newtonian nature, it shows a tendency to swell or contract at the die exit depending on its rheological parameters. The die swell phenomenon is a typical example of a free surface problem where the free surface is formed at the die exit after the polymeric fluid has been extruded. The swelling process leads to an undesired increase in the dimensions of the extrudate. To be able to obtain a near-net shape product, the flow in the extrusion process should be well-understood to shed some light on the important process parameters behind the swelling phenomenon. To this end, a systematic study has been carried out to compare constitutive models proposed in literature for second-order fluids in terms of their ability to capture the physics behind the swelling phenomenon. The effect of various process and rheological parameters on the die swell such as the extrusion velocity, normal stress coefficients, and Reynolds and Deborah numbers have also been investigated. The models developed here can predict both swelling and contraction of the extrudate successfully. The die swell problem was solved for a wide range of Deborah numbers and for two different Re numbers. The numerical model was validated through the solution of fully developed Newtonian and Non-Newtonian viscoelastic flows in a two-dimensional channel, and the results of these two benchmark problems were compared with analytic solutions, and good agreements were obtained.


2016 ◽  
Vol 14 (05) ◽  
pp. 1750050 ◽  
Author(s):  
Lu Wang ◽  
Fei Xu ◽  
Yang Yang

Compared with the traditional Smoothed Particle Hydrodynamics (SPH), Finite Particle Method (FPM) has higher accuracy for boundary region. However, there are still two inherent defects which are the time consuming and the numerical instability in FPM. In this paper, a high-order algorithm based on the Taylor series expansion and the matrix theory is proposed and the corresponding particles selected mode is discussed. It is validated that the algorithm has higher-order accuracy than the previous low-order improvement algorithm for FPM. Further, transient heat conduction examples have been discussed to verify the feasibility and effectiveness of the new algorithm.


Sign in / Sign up

Export Citation Format

Share Document