Trying to simulate a case of a Stirred Tank
Hi everyone.
I have been working in a simulation of a stirred tank using DualSPHysics, the simulations are based on experiment realized by Hartmann (Hartmann, H., et al. "Assessment of large eddy and RANS stirred tank simulations by means of LDA."). I want to replicate this simulations using SPH, but, I have a problem with the convergence of the solution.
I have put some "probes" inside the tank and I have been checking the velocity and the pressure over the time. The velocity at the start and middle of the simulation is behave good, but the pressure, has a really curious behavior. The pressure goes up at the start of the simulation (obviously) and then goes down until 0 and at this moment the my simulation start to fail (the density increase and reduce dramatically) and the velocity goes to 0 as if the impeller was not moving.
I have attached some graphs that show what is happening:
Press:
https://www.dropbox.com/s/j21pzym5o85dkv6/Press.png?dl=0
Velocity:
https://www.dropbox.com/s/ho04drd7sdor0yp/Velocity.png?dl=0
Someone has an idea that what could be the problem?
I'm using Wenland Kernel, Sympletic algorithm, CFL=0.2, Sounspeed=280 m/s. If you need more information I can send you the xml file.
Thank so much
Regards.
I have been working in a simulation of a stirred tank using DualSPHysics, the simulations are based on experiment realized by Hartmann (Hartmann, H., et al. "Assessment of large eddy and RANS stirred tank simulations by means of LDA."). I want to replicate this simulations using SPH, but, I have a problem with the convergence of the solution.
I have put some "probes" inside the tank and I have been checking the velocity and the pressure over the time. The velocity at the start and middle of the simulation is behave good, but the pressure, has a really curious behavior. The pressure goes up at the start of the simulation (obviously) and then goes down until 0 and at this moment the my simulation start to fail (the density increase and reduce dramatically) and the velocity goes to 0 as if the impeller was not moving.
I have attached some graphs that show what is happening:
Press:
https://www.dropbox.com/s/j21pzym5o85dkv6/Press.png?dl=0
Velocity:
https://www.dropbox.com/s/ho04drd7sdor0yp/Velocity.png?dl=0
Someone has an idea that what could be the problem?
I'm using Wenland Kernel, Sympletic algorithm, CFL=0.2, Sounspeed=280 m/s. If you need more information I can send you the xml file.
Thank so much
Regards.