32 research outputs found
Recommended from our members
Multi-Scale Models to Simulate Interactions between Liquid and Thin Structures
In this dissertation, we introduce a framework for simulating the dynamics between liquid and thin structures, including the effects of buoyancy, drag, capillary cohesion, dripping, and diffusion. After introducing related works, Part I begins with a discussion on the interactions between Newtonian fluid and fabrics. In this discussion, we treat both the fluid and the fabrics as continuum media; thus, the physical model is built from mixture theory. In Part II, we discuss the interactions between Newtonian fluid and hairs. To have more detailed dynamics, we no longer treat the hairs as continuum media. Instead, we treat them as discrete Kirchhoff rods. To deal with the thin layer of liquid that clings to the hairs, we augment each hair strand with a height field representation, through which we introduce a new reduced-dimensional flow model to solve the motion of liquid along the longitudinal direction of each hair. In addition, we develop a faithful model for the hairs' cohesion induced by surface tension, where a penalty force is applied to simulate the collision and cohesion between hairs. To enable the discrete strands interact with continuum-based, shear-dependent liquid, in Part III, we develop models that account for the volume change of the liquid as it passes through strands and the momentum exchange between the strands and the liquid. Accordingly, we extend the reduced-dimensional flow model to simulate liquid with elastoviscoplastic behavior. Furthermore, we use a constraint-based model to replace the penalty-force model to handle contact, which enables an accurate simulation of the frictional and adhesive effects between wet strands. We also present a principled method to preserve the total momentum of a strand and its surface flow, as well as an analytic plastic flow approach for Herschel-Bulkley fluid that enables stable semi-implicit integration at larger time steps.
We demonstrate a wide range of effects, including the challenging animation scenarios involving splashing, wringing, and colliding of wet clothes, as well as flipping of hair, animals shaking, spinning roller brushes from car washes being dunked in water, and intricate hair coalescence effects. For complex liquids, we explore a series of challenging scenarios, including strands interacting with oil paint, mud, cream, melted chocolate, and pasta sauce
Thickness and Thermal Conductivities of the Walls and Fluid Layer Effects on the Onset of Thermal Convection in a Horizontal Fluid Layer Heated from Below
The thermal boundary conditions have important effects on the hydrodynamics of a thermoâconvective fluid layer. These effects are introduced through the Biot number under the Robin type boundary conditions. The thermal conductivity and thicknesses of the walls are key properties to bridge two known ideal situations widely studied: the fluid layer bounded by two insulating walls and the fluid layer bounded by two perfect thermal conducting walls. This chapter is devoted to the physical mechanisms involved in the thermal boundary conditions, its influence on the linear stability of the fluid layer and its implications with the pattern formation. A review of very important investigations on the subject is also given. The role of the thermal conductivities and thicknesses of the walls is explained with help of curves of criticality for the thermoconvection in a horizontal Newtonian fluid layer
Conformation constraints for efficient viscoelastic fluid simulation
The simulation of high viscoelasticity poses important computational challenges. One is the difficulty to robustly measure strain and its derivatives in a medium without permanent structure. Another is the high stiffness of the governing differential equations. Solutions that tackle these challenges exist, but they are computationally slow. We propose a constraint-based model of viscoelasticity that enables efficient simulation of highly viscous and viscoelastic phenomena. Our model reformulates, in a constraint-based fashion, a constitutive model of viscoelasticity for polymeric fluids, which defines simple governing equations for a conformation tensor. The model can represent a diverse palette of materials, spanning elastoplastic, highly viscous, and inviscid liquid behaviors. In addition, we have designed a constrained dynamics solver that extends the position-based dynamics method to handle efficiently both position-based and velocity-based constraints. We show results that range from interactive simulation of viscoelastic effects to large-scale simulation of high viscosity with competitive performance
A Physically Consistent Implicit Viscosity Solver for SPH Fluids
In this paper, we present a novel physically consistent implicit solver for the simulation of highly viscous fluids using the Smoothed Particle Hydrodynamics (SPH) formalism. Our method is the result of a theoretical and practical inâdepth analysis of the most recent implicit SPH solvers for viscous materials. Based on our findings, we developed a list of requirements that are vital to produce a realistic motion of a viscous fluid. These essential requirements include momentum conservation, a physically meaningful behavior under temporal and spatial refinement, the absence of ghost forces induced by spurious viscosities and the ability to reproduce complex physical effects that can be observed in nature. On the basis of several theoretical analyses, quantitative academic comparisons and complex visual experiments we show that none of the recent approaches is able to satisfy all requirements. In contrast, our proposed method meets all demands and therefore produces realistic animations in highly complex scenarios. We demonstrate that our solver outperforms former approaches in terms of physical accuracy and memory consumption while it is comparable in terms of computational performance. In addition to the implicit viscosity solver, we present a method to simulate melting objects. Therefore, we generalize the viscosity model to a spatially varying viscosity field and provide an SPH discretization of the heat equation
Turning point principle for relativistic stars
Upon specifying an equation of state, spherically symmetric steady states of
the Einstein-Euler system are embedded in 1-parameter families of solutions,
characterized by the value of their central redshift. In the 1960's Zel'dovich
[50] and Wheeler [22] formulated a turning point principle which states that
the spectral stability can be exchanged to instability and vice versa only at
the extrema of mass along the mass-radius curve. Moreover the bending
orientation at the extrema determines whether a growing mode is gained or lost.
We prove the turning point principle and provide a detailed description of the
linearized dynamics. One of the corollaries of our result is that the number of
growing modes grows to infinity as the central redshift increases to infinity.Comment: 33 pages, 1 figur
Variational Stokes: A Unified Pressure-viscosity Solver for Accurate Viscous Liquids
© ACM, 2017. This is the author's version of the work. It is posted here by permission of ACM for your personal use. Not for redistribution. The definitive version was published in Larionov, E., Batty, C., & Bridson, R. (2017). Variational Stokes: A Unified Pressure-viscosity Solver for Accurate Viscous Liquids. ACM Trans. Graph., 36(4), 101:1â101:11. https://doi.org/10.1145/3072959.3073628We propose a novel unsteady Stokes solver for coupled viscous and pressure forces in grid-based liquid animation which yields greater accuracy and visual realism than previously achieved. Modern fluid simulators treat viscosity and pressure in separate solver stages, which reduces accuracy and yields incorrect free surface behavior. Our proposed implicit variational formulation of the Stokes problem leads to a symmetric positive definite linear system that gives properly coupled forces, provides unconditional stability, and treats difficult boundary conditions naturally through simple volume weights. Surface tension and moving solid boundaries are also easily incorporated. Qualitatively, we show that our method recovers the characteristic rope coiling instability of viscous liquids and preserves fine surface details, while previous grid-based schemes do not. Quantitatively, we demonstrate that our method is convergent through grid refinement studies on analytical problems in two dimensions. We conclude by offering practical guidelines for choosing an appropriate viscous solver, based on the scenario to be animated and the computational costs of different methods.Natural Sciences and Engineering Research Council of Canad
Turning Point Principle for Relativistic Stars
Upon specifying an equation of state, spherically symmetric steady states of the Einstein-Euler system are embedded in 1-parameter families of solutions, characterized by the value of their central redshift. In the 1960âs Zelâdovich (Voprosy Kosmogonii 9:157â170, 1963) and Harrison et al. (Gravitation Theory and Gravitational Collapse. The University of Chicago press, Chicago, 1965) formulated a turning point principle which states that the spectral stability can be exchanged to instability and vice versa only at the extrema of mass along the mass-radius curve. Moreover the bending orientation at the extrema determines whether a growing mode is gained or lost. We prove the turning point principle and provide a detailed description of the linearized dynamics. One of the corollaries of our result is that the number of growing modes grows to infinity as the central redshift increases to infinity