SPH boundary & discontinuity correction terms
I'm working on a CFD Matlab code using smoothing particle hydrodynamics (SPH) approach.
I have a great problem with the consistency at the boundaries & the discontinuities due to Shock waves. They make the code results totally wrong. I've tried to add terms for normalization & correction defined by (Liu) & (Chen) in the governing equations but I fail to get right answers. Could anyone help me how to use the normalization & correction techniques in Navier-Stokes equations especially in discontinuities |
What are you using to model your boundary conditions??
|
I used virtual exterior particles to mirror the fluid particles near the boundary. I didn't use boundary virtual particles.
Now I achieved the right results near the boundary. But in the wave discontinuity, I have an inaccuracy in result (represented in form of varying non-smooth results for a flow). |
So you have used ghost particles...?
And are you also using Artificial Viscosity SPH? |
Yes I used Monaghan Artificial Viscosity
|
As far as I know with AV SPH you will always encounter some problem with your results. If you want more quality and precision you should look for preciser variants of SPH like Riemann-SPH and delta-SPH...
|
OK I'll check & try that....many thanks
|
All times are GMT -4. The time now is 05:40. |