-
Velocity-Based Monte Carlo Fluids
Authors:
Ryusuke Sugimoto,
Christopher Batty,
Toshiya Hachisuka
Abstract:
We present a velocity-based Monte Carlo fluid solver that overcomes the limitations of its existing vorticity-based counterpart. Because the velocity-based formulation is more commonly used in graphics, our Monte Carlo solver can be readily extended with various techniques from the fluid simulation literature. We derive our method by solving the Navier-Stokes equations via operator splitting and d…
▽ More
We present a velocity-based Monte Carlo fluid solver that overcomes the limitations of its existing vorticity-based counterpart. Because the velocity-based formulation is more commonly used in graphics, our Monte Carlo solver can be readily extended with various techniques from the fluid simulation literature. We derive our method by solving the Navier-Stokes equations via operator splitting and designing a pointwise Monte Carlo estimator for each substep. We reformulate the projection and diffusion steps as integration problems based on the recently introduced walk-on-boundary technique [Sugimoto et al. 2023]. We transform the volume integral arising from the source term of the pressure Poisson equation into a form more amenable to practical numerical evaluation. Our resulting velocity-based formulation allows for the proper simulation of scenes that the prior vorticity-based Monte Carlo method [Rioux-Lavoie and Sugimoto et al. 2022] either simulates incorrectly or cannot support. We demonstrate that our method can easily incorporate advancements drawn from conventional non-Monte Carlo methods by showing how one can straightforwardly add buoyancy effects, divergence control capabilities, and numerical dissipation reduction methods, such as advection-reflection and PIC/FLIP methods.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Reach For the Spheres: Tangency-Aware Surface Reconstruction of SDFs
Authors:
Silvia Sellán,
Christopher Batty,
Oded Stein
Abstract:
Signed distance fields (SDFs) are a widely used implicit surface representation, with broad applications in computer graphics, computer vision, and applied mathematics. To reconstruct an explicit triangle mesh surface corresponding to an SDF, traditional isosurfacing methods, such as Marching Cubes and and its variants, are typically used. However, these methods overlook fundamental properties of…
▽ More
Signed distance fields (SDFs) are a widely used implicit surface representation, with broad applications in computer graphics, computer vision, and applied mathematics. To reconstruct an explicit triangle mesh surface corresponding to an SDF, traditional isosurfacing methods, such as Marching Cubes and and its variants, are typically used. However, these methods overlook fundamental properties of SDFs, resulting in reconstructions that exhibit severe oversmoothing and feature loss. To address this shortcoming, we propose a novel method based on a key insight: each SDF sample corresponds to a spherical region that must lie fully inside or outside the surface, depending on its sign, and that must be tangent to the surface at some point. Leveraging this understanding, we formulate an energy that gauges the degree of violation of tangency constraints by a proposed surface. We then employ a gradient flow that minimizes our energy, starting from an initial triangle mesh that encapsulates the surface. This algorithm yields superior reconstructions to previous methods, even with sparsely sampled SDFs. Our approach provides a more nuanced understanding of SDFs and offers significant improvements in surface reconstruction.
△ Less
Submitted 18 August, 2023;
originally announced August 2023.
-
A Closest Point Method for Surface PDEs with Interior Boundary Conditions for Geometry Processing
Authors:
Nathan King,
Haozhe Su,
Mridul Aanjaneya,
Steven Ruuth,
Christopher Batty
Abstract:
Many geometry processing techniques require the solution of partial differential equations (PDEs) on surfaces. Such surface PDEs often involve boundary conditions prescribed on the surface, at points or curves on its interior or along the geometric (exterior) boundary of an open surface. However, input surfaces can take many forms (e.g., meshes, parametric surfaces, point clouds, level sets, neura…
▽ More
Many geometry processing techniques require the solution of partial differential equations (PDEs) on surfaces. Such surface PDEs often involve boundary conditions prescribed on the surface, at points or curves on its interior or along the geometric (exterior) boundary of an open surface. However, input surfaces can take many forms (e.g., meshes, parametric surfaces, point clouds, level sets, neural implicits). One must therefore generate a mesh to apply finite element-type techniques or derive specialized discretization procedures for each surface representation.
We propose instead to address such problems through a novel extension of the closest point method (CPM) to handle interior boundary conditions specified at surface points or curves. CPM solves the surface PDE by solving a volumetric PDE defined over the Cartesian embedding space containing the surface; only a closest point function is required to represent the surface. As such, CPM supports surfaces that are open or closed, orientable or not, and of any codimension or even mixed-codimension. To enable support for interior boundary conditions, we develop a method to implicitly partition the embedding space across interior boundaries. CPM's finite difference and interpolation stencils are adapted to respect this partition while preserving second-order accuracy. Furthermore, an efficient sparse-grid implementation and numerical solver is developed that can scale to tens of millions of degrees of freedom, allowing PDEs to be solved on more complex surfaces. We demonstrate our method's convergence behaviour on selected model PDEs. Several geometry processing problems are explored: diffusion curves on surfaces, geodesic distance, tangent vector field design, and harmonic map construction. Our proposed approach thus offers a powerful and flexible new tool for a range of geometry processing tasks on general surface representations.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
A Practical Walk-on-Boundary Method for Boundary Value Problems
Authors:
Ryusuke Sugimoto,
Terry Chen,
Yiti Jiang,
Christopher Batty,
Toshiya Hachisuka
Abstract:
We introduce the walk-on-boundary (WoB) method for solving boundary value problems to computer graphics. WoB is a grid-free Monte Carlo solver for certain classes of second order partial differential equations. A similar Monte Carlo solver, the walk-on-spheres (WoS) method, has been recently popularized in computer graphics due to its advantages over traditional spatial discretization-based altern…
▽ More
We introduce the walk-on-boundary (WoB) method for solving boundary value problems to computer graphics. WoB is a grid-free Monte Carlo solver for certain classes of second order partial differential equations. A similar Monte Carlo solver, the walk-on-spheres (WoS) method, has been recently popularized in computer graphics due to its advantages over traditional spatial discretization-based alternatives. We show that WoB's intrinsic properties yield further advantages beyond those of WoS. Unlike WoS, WoB naturally supports various boundary conditions (Dirichlet, Neumann, Robin, and mixed) for both interior and exterior domains. WoB builds upon boundary integral formulations, and it is mathematically more similar to light transport simulation in rendering than the random walk formulation of WoS. This similarity between WoB and rendering allows us to implement WoB on top of Monte Carlo ray tracing, and to incorporate advanced rendering techniques (e.g., bidirectional estimators with multiple importance sampling, the virtual point lights method, and Markov chain Monte Carlo) into WoB. WoB does not suffer from the intrinsic bias of WoS near the boundary and can estimate solutions precisely on the boundary. Our numerical results highlight the advantages of WoB over WoS as an attractive alternative to solve boundary value problems based on Monte Carlo.
△ Less
Submitted 19 May, 2023; v1 submitted 7 May, 2023;
originally announced May 2023.
-
Curl-Flow: Boundary-Respecting Pointwise Incompressible Velocity Interpolation for Grid-Based Fluids
Authors:
Jumyung Chang,
Ruben Partono,
Vinicius C. Azevedo,
Christopher Batty
Abstract:
We propose to augment standard grid-based fluid solvers with pointwise divergence-free velocity interpolation, thereby ensuring exact incompressibility down to the sub-cell level. Our method takes as input a discretely divergence-free velocity field generated by a staggered grid pressure projection, and first recovers a corresponding discrete vector potential. Instead of solving a costly vector Po…
▽ More
We propose to augment standard grid-based fluid solvers with pointwise divergence-free velocity interpolation, thereby ensuring exact incompressibility down to the sub-cell level. Our method takes as input a discretely divergence-free velocity field generated by a staggered grid pressure projection, and first recovers a corresponding discrete vector potential. Instead of solving a costly vector Poisson problem for the potential, we develop a fast parallel sweeping strategy to find a candidate potential and apply a gauge transformation to enforce the Coulomb gauge condition and thereby make it numerically smooth. Interpolating this discrete potential generates a pointwise vector potential whose analytical curl is a pointwise incompressible velocity field. Our method further supports irregular solid geometry through the use of level set-based cut-cells and a novel Curl-Noise-inspired potential ramping procedure that simultaneously offers strictly non-penetrating velocities and incompressibility. Experimental comparisons demonstrate that the vector potential reconstruction procedure at the heart of our approach is consistently faster than prior such reconstruction schemes, especially those that solve vector Poisson problems. Moreover, in exchange for its modest extra cost, our overall Curl-Flow framework produces significantly improved particle trajectories that closely respect irregular obstacles, do not suffer from spurious sources or sinks, and yield superior particle distributions over time.
△ Less
Submitted 26 November, 2023; v1 submitted 1 April, 2021;
originally announced April 2021.
-
Asynchronous Liquids: Regional Time Stepping for Faster SPH and PCISPH
Authors:
Prashant Goswami,
Christopher Batty
Abstract:
This paper presents novel and efficient strategies to spatially adapt the amount of computational effort applied based on the local dynamics of a free surface flow, for both classic weakly compressible SPH (WCSPH) and predictive-corrective incompressible SPH (PCISPH). Using a convenient and readily parallelizable block-based approach, different regions of the fluid are assigned differing time step…
▽ More
This paper presents novel and efficient strategies to spatially adapt the amount of computational effort applied based on the local dynamics of a free surface flow, for both classic weakly compressible SPH (WCSPH) and predictive-corrective incompressible SPH (PCISPH). Using a convenient and readily parallelizable block-based approach, different regions of the fluid are assigned differing time steps and solved at different rates to minimize computational cost. Our approach for WCSPH scheme extends an asynchronous SPH technique from compressible flow of astrophysical phenomena to the incompressible free surface setting, and further accelerates it by entirely decoupling the time steps of widely spaced particles. Similarly, our approach to PCISPH adjusts the the number of iterations of density correction applied to different regions, and asynchronously updates the neighborhood regions used to perform these corrections; this sharply reduces the computational cost of slowly deforming regions while preserving the standard density invariant. We demonstrate our approaches on a number of highly dynamic scenarios, demonstrating that they can typically double the speed of a simulation compared to standard methods while achieving visually consistent results.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
Addressing Troubles with Double Bubbles: Convergence and Stability at Multi-Bubble Junctions
Authors:
Yun Fei,
Christopher Batty,
Eitan Grinspun
Abstract:
In this report we discuss and propose a correction to a convergence and stability issue occurring in the work of Da et al.[2015], in which they proposed a numerical model to simulate soap bubbles.
In this report we discuss and propose a correction to a convergence and stability issue occurring in the work of Da et al.[2015], in which they proposed a numerical model to simulate soap bubbles.
△ Less
Submitted 11 June, 2020; v1 submitted 14 October, 2019;
originally announced October 2019.
-
Constraint Bubbles: Adding Efficient Zero-Density Bubbles to Incompressible Free Surface Flow
Authors:
Ryan Goldade,
Christopher Batty
Abstract:
Liquid simulations for computer animation often avoid simulating the air phase to reduce computational costs and ensure good conditioning of the linear systems required to enforce incompressibility. However, this free surface assumption leads to an inability to realistically treat bubbles: submerged gaps in the liquid are interpreted as empty voids that immediately collapse. To address this shortc…
▽ More
Liquid simulations for computer animation often avoid simulating the air phase to reduce computational costs and ensure good conditioning of the linear systems required to enforce incompressibility. However, this free surface assumption leads to an inability to realistically treat bubbles: submerged gaps in the liquid are interpreted as empty voids that immediately collapse. To address this shortcoming, we present an efficient, practical, and conceptually simple approach to augment free surface flows with negligible density bubbles. Our method adds a new constraint to each disconnected air region that guarantees zero net flux across its entire surface, and requires neither simulating both phases nor reformulating into stream function variables. Implementation of the method requires only minor modifications to the pressure solve of a standard grid-based fluid solver, and yields linear systems that remain sparse and symmetric positive definite. In our evaluations, solving the modified pressure projection system took no more than 10% longer than the corresponding free surface solve. We demonstrate the method's effectiveness and flexibility by incorporating it into commercial fluid animation software and using it to generate a variety of dynamic bubble scenarios showcasing glugging effects, viscous and inviscid bubbles, interactions with irregularly-shaped and moving solid boundaries, and surface tension effects.
△ Less
Submitted 30 November, 2017;
originally announced November 2017.
-
Multimaterial Front Tracking
Authors:
Fang Da,
Christopher Batty,
Eitan Grinspun
Abstract:
We present the first triangle mesh-based technique for tracking the evolution of general three-dimensional multimaterial interfaces undergoing complex topology changes induced by deformations and collisions. Our core representation is a non-manifold triangle surface mesh with material labels assigned to each half-face to distinguish volumetric regions. We advect the vertices of the mesh in a Lagra…
▽ More
We present the first triangle mesh-based technique for tracking the evolution of general three-dimensional multimaterial interfaces undergoing complex topology changes induced by deformations and collisions. Our core representation is a non-manifold triangle surface mesh with material labels assigned to each half-face to distinguish volumetric regions. We advect the vertices of the mesh in a Lagrangian manner, and employ a complete set of collision-safe mesh improvement and topological operations that track and update material labels. In particular, we develop a unified, collision-safe strategy for handling complex topological operations acting on evolving triple- and higher-valence junctions, and a flexible method to merge colliding multimaterial meshes. We demonstrate our approach with a number of challenging geometric flows, including passive advection, normal flow, and mean curvature flow.
△ Less
Submitted 20 January, 2014; v1 submitted 13 June, 2013;
originally announced June 2013.