trixi-framework / trixi-framework/TrixiParticles.jl
Sticking artifacts and negative pressure of dummy boundary particles
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 83
- Forks
- 24
- Avg merge
- 8d 18h
- Merged PRs (30d)
- 3
Description
With ContinuityDensity, we allow negative particle pressure and corresponding attractive forces between particles. This is a good thing, as clipping the pressure results in unphysical rarefaction of the fluid in some situations. A realistic incompressible or weakly compressible fluid can't just be pulled apart to a significantly lower density.
However, we usually (always?) don't want attractive forces between fluid and boundary particles on (free-slip?) walls. But with dummy particle boundaries, we do get attractive forces if the pressure becomes negative.
To demonstrate the severity of the problem, here is a dam break with dummy boundary particles using SummationDensity (i.e. their density is calculated through kernel summation and their pressure from the state equation). This is the result when the problem is not treated at all:
With SummationDensity, all boundary particles away from the fluid end up with a large negative pressure. Therefore, the fluid is sucked up the walls of the tank (I should've used a tank without ceiling, that would look even funnier).
Currently, we just limit the pressure of dummy boundary particles to be non-negative. This eliminates the sucking, but still produces some sticking artifacts, as there are still attractive forces when a fluid particle close to the boundary has a negative pressure. The example looks like this, then.
As can be seen very clearly, there are still sticking artifacts in the form of fluid particles gliding along the boundary.
Another strategy that came to my mind is to limit the coefficient in front of the kernel gradient in the pressure acceleration to be non-negative for fluid-boundary interaction, eliminating all attractive forces between fluid and boundary particles. This is the result:
Obviously, without attractive forces, the fluid will bounce off the boundary after the impact.
I created this issue to discuss options and see if we can find anything about this in the literature (I couldn't).
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Review the ContinuityDensity and SummationDensity boundary-pressure behavior described in the issue, along with the three dam-break results. Compare the existing non-negative dummy-boundary pressure limit with the proposed pressure-acceleration coefficient limit and relevant literature; done means selecting and documenting a physically justified approach to eliminate the sticking artifacts.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- julia
- Domain
- hpc
- Issue type
- Bug
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100