Vellum Solvers are a type of POP Solver that include support for internal triangle and edge collisions and explicit constraint geometry.
The Geometry data stores the dynamic attributes of the simulation as point attributes. It also stores the collision geometry, so the triangle or poly-line mesh is used to detect and resolve collisions. However, the way the points are held together is determined by the separate Constraint Geometry data. The Constraint Geometry data should have a one-to-one point correspondence with the Geometry. Polygons in the Constraint Geometry are interpreted as rules to maintain certain properties of the source points, such as distance or bend angle.
POP Forces can be wired into the solver-inputs to adjust the behavior of the object.
Controls if any collision detection is done at all.
Controls if self collisions are detected.
Scales the effective time of the Vellum solve. This can be used to create bullet-time like effects where the physics of the Vellum solver run at a different rate than the Houdini playbar. A value of 2 will cause cloth to fall twice as fast, and 0.1 will slow it to a tenth the speed.
Each frame will be broken into this number of substeps. Additional substeps are required for fast moving collisions or sudden forces.
The default substeps can be very aggressive, usually if the Vellum solver is too stretchy, raising substeps to 2 or 5 is a good first start.
Within each substep, this number of passes will be taken by the constraint enforcement operations. Stiff constraints can require more iterations to converge. A good starting point is the diameter of the geometry - the number of edges between the farthest points.
The default constraint iterations use a Gauss-Seidel approach that is fast to converge. However, if it doesn’t fully converge due to too high stiffness, or impossible configurations, it will leave the error as bad looking triangles. The smoothing iterations use a Jacobi approach which is slower to converge but leaves error spread out in a more attractive fashion. The default of ten passes helps smooth out error, but might need to be increased if the overall Constraint Iterations is very high.
Number of collision detection passes to perform. These are interleaved between the constraint iterations. Since collision is expensive, it is best to minimize this. But frequent interleaving helps avoid tent-poling effects where a small collider is fighting with the no-stretch constraints. In practice we find 10 to be correct for most situations, and substeps often being a better solution to increase quality.
Post Collision Passes
After all constraints are performed, a final round of collision detection is done. Collisions are often the most noticeable failure mode, and it is ideal if the next frame can start with non-intersecting geometry. Thus a final cleanup pass can achieve these requirements. We have found that "number of stacked layers + 2" is a good estimate for this number. This allows the effect of the underlying collider to ripple through the stacked layers fully.
In any collision pass, any colliding pair may not be fully resolved. This number of additional collider-pair passes will be run until they are resolved. Since these are only performed on active colliders (and no new collision search is done) this is very cheap.
layer integer point attribute is used to flag point as belonging
to different layers of cloth. Higher numbers refer to higher layers.
Layer Shock will make lower layers this many times heavier during collision
evaluation, ensuring the higher layers will move out of their way. The
rest of the dynamics are unaffected by this, and the difference is fixed
regardless of the number of layers between the two. This can be thought
of as a way to dial between one-way layering of sims and fully coupled
A threshold at which to apply full friction. When the ratio of the tangential velocity and the normal impulse is less than this, the tangential velocity will be fully eliminated through friction. This is roughly tan() of the slope angle that will allow sliding under gravity.
If the static threshold fails, this controls what percentage the tangential velocity will be reduced in the dynamic friction case.
A scale factor on the amount of friction effect to apply for collisions with external geometry.
A scale factor on the amount of friction effect to apply for collisions with self geometry.
Static SDF Enable
A scale factor for the volume collider’s static friction. Useful to create frictionless grounds.
Dynamic SDF Enable
A scale factor for the volume collider’s dynamic friction. Useful to create frictionless grounds.
When integrating points forward by their velocities, curved motions like swinging pendulums will be mispredicted and be damped by the constraint projections. A second order predictor recovers more of this curved motion and maintains more energy in the system.
In any case, if collisions are detected the system will fall back to first order to avoid excessive bouncing.
The Vellum points will receive forces from DOP forces. Turning this off will remove this effect, and the points will only be affected by POP forces.
mass attribute is on the incoming geometry, this will be used as the
default value for the points.
Usually mass should already be set by the Vellum Constraints SOP.
If the incoming geometry has no
pscale attribute, it will be set to this value.
Usually thickness should already be set by the Vellum Constraints SOP.
Disable Failed Collisions
If points fail to resolve their collisions after the post-collision passes, they will be flagged as disabled, allowing them and any primitives attached to them to move collision free.
Reset Collided Points
If a disabled point detects it is no longer in a tangled situation, and is connected to non-disabled point, it will re-activate itself. The hope is that it has moved back to the correct side of the geometry.
Initialize Overlap Distances
overlap_external attributes will be created
and initialized to respect the initial setup of the solver.
Note if the attributes are already present, they will not be initialized. If the configuration has changed due to modeling operations you may need to delete those attributes.
Update Overlap Distances
overlap_external will be updated to reflect
the current configuration. The only decrease with this operation,
so to initialize set to a large initial value.
Tetrahedral Open Faces
The unshared faces of a tetrahedral mesh will generate collision geometry if this is set. Sometimes one wishes to use both a triangle shell and a tetrahedral mesh, making this operation redundant.
Various acceleration limiting options can be used to prevent the simulation from being over-eager to obey non-realistic forces. These can otherwise result in large energy spikes.
Fallback to First Order Integration on Collision
If Max Acceleration is exceeded for a point, assume it signals a sharp, discontinuous collision where the second order prediction will be wrong and add erroneous motion, often in the form of bouncing. In this case, fallback to first order integration for the affected point.
Limit Acceleration in Velocity Update
Cap the amount that velocity of a particle is allowed to change as a result of any of the dynamics. This is useful to prevent some instant motion being mis-identified as a massive force, and thus avoid fly-away particles.
Limit Displacement on Collision
During the collision resolution, if the collision correction moves a particle more than the acceleration amount, cap the effect. The hope is to fail more gracefully when a part of a model snags on the cloth and is pulling it in a surprising fashion.
How frequently welds and constraints are tested to see if they have reached their breaking point. Higher rates will give more accuracy to the solve, but changes to the topology from a broken constraint can slow down the solve considerably.
Do not test for any breaking.
Check only once every frame, regardless of number of substeps.
Check every substep.
Enable Auto Sleep
Particles that remain near-stationary for sufficient time will be
put to sleep by setting their
stopped attribute to
You can use the Pinned visualization to visualize which points are sleeping.
How slowly particles need to be moving before they are eligible
for sleeping. The amount of time they remain below this speed
is accumulated in the
deactivation_time attribute. If they
exceed this speed, the
deactivation_time is reset to zero.
How many seconds particles need to remain at or below the velocity threshold before they are put to sleep. See the Wake Up POP for methods to wake the particles back up again.
Vellum grains use a separate collision pass to deal with collisions between the grains. These parameters control that behavior. They are similar to those on POP Grains
Potential intersection particles are any within this scaled distance
of the average of the two particles
pscale attribute. This
is an overestimate because usually collisions are not updated
during the constraint iterations, so it needs to record not
just the currently colliding particles, but those that may
start to collide due to the earlier iterations.
This also effects the range of the attraction force in clumping.
The maximum number of particles that will be considered when searching for potential collisions over the substep. Capping this is useful to avoid excessive computations, if too many particles are created at one spot.
A universal weight scale applied repulsion weight, attraction weight, and the friction weight. Because constraint averaging isn’t used, this may need to be reduced for stability, especially when high attraction stiffness values are used. A value of 1 / Max Neighbours should provide stability, but at the cost of slower convergence.
Assume Uniform Radius
pscale attribute is used to determine the radius of each
particle. If all particles have the same radius, faster acceleration
structures can be used to find neighbors.
A weighting for how much the particle collision forces are weighted. A value of zero will disable particle collision.
How strongly particles are kept apart. Higher values result in less bouncy repulsion.
A weighting for how much the particles will naturally stick together when close. A value of zero will disable particle clumping.
Scale with the
attractionweight point attribute.
How strongly nearby particles stick to each other. Higher values result in a less bouncy adhesion.
Enable Mass Shock
Artificially scales the mass of particles according to their position with respect to gravity. By making particles higher up lighter, stacks of particles will converge faster and be more stable.
Shock Scaling Power
The amount of scaling to perform. Higher numbers increase the contrast between successive particles. A value of 0 will cause no ratio between particles, a value of 1 a 15% ratio between two particles stacked vertically.
Too high a number makes higher particles extremely light and destabilizes the system.
The up-vector used to define a gradient of particle masses.
Should be in the direction of stacking, so usually is opposite to that of gravity.
Specify any desired compile flags for the kernel.
The Apple OSX OpenCL compiler requires only a single space between kernel options!
When Finish Kernels is disabled, no attempt is to wait for the OpenCL kernels to complete before continuing the next solver. This lets them run in the background until their results are actually needed. To simplify debugging or timing, it is useful to ensure kernels are finished to make sure errors are detected in the right spot.
When loading kernels from disk the kernel is cached to avoid regenerating it every solve. Turning this on forces the re-loading and recompiling of the kernel. This is useful if #include files refer to code that has changed, or the kernel file is changed in an external text editor.
It should always be disabled when prototyping is complete.