CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

How to tackle volume loss in FLIP?

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 9, 2016, 18:09
Smile How to tackle volume loss in FLIP?
  #1
New Member
 
Join Date: Nov 2016
Posts: 4
Rep Power: 10
usirev is on a distinguished road
Hi,

I created a primitive FLIP simulation - a starting point which I'll further enhance in terms of performance.

Right now I have a problem with volume loss - when external forces are
applied to particles, the volume of the liquid gets significantly smaller after a
while (30% of starting fluid volume).

Next follows description my current implementation of FLIP. I'll be glad for
any advice on how to improve my implementation and most importantly,
what plays the major role in the volume loss.

(How can I constrain the fluid volume in FLIP?)


Info about my primitive solver:
The main loop:
  1. Mark cells containing fluid.
  2. Convert velocity from particles -> grid. (Inverse trilinear interpolation - particle transfers its velocity to the closest 4 grid points in each X, Y and Z velocity grids)
  3. Apply external forces to grid (just add gravity force to every vertex on the Y-velocity grid).
  4. Compute divergence.
  5. Solve the (pseudo)pressure Poisson equation using Gauss-Seidel. Subtract pressure gradient.
  6. Convert velocity from grid -> particles. (Trilinear interpolation - similar to step 2.).
  7. Integrate using Forward Euler.

Other info:
  • I do not extrapolate velocity outside fluid cells.
  • I do not use CFL condition - the time step is always 16ms.
  • Boundary conditions are of type free-slip.
  • I do not remove or add particles.
  • I clamp the integrated positions of particles to stay withing fluid domain ([0;1]x[0;1]x[0;1] cube).
    (Can this be the potential source of volume loss?)
usirev is offline   Reply With Quote

Reply

Tags
flip, volume change


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
multiphaseEulerFoam high Courant number Frenk_T OpenFOAM 5 November 24, 2016 04:23
[blockMesh] error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Meshing & Mesh Conversion 2 March 14, 2012 10:56
On the damBreak4phaseFine cases paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14
[Commercial meshers] CuBit t42 OpenFOAM Meshing & Mesh Conversion 6 July 10, 2008 08:51
Volume Loss in Level Set Method Zi-Wei Chiou Main CFD Forum 6 September 19, 2001 10:37


All times are GMT -4. The time now is 14:15.