An efficient heat-based model for solid-liquid-gas phase transition and dynamic interaction

For melting simulation, solid-liquid coupling, liquid-gas interaction, bubble/foam generation, etc., many new methods have been emerging in recent years in computer graphics. To further push advance the technical frontier of the aforementioned phenomena, our novel solution is to focus on an efficient heat-based method towards faithful simulation of physical procedures pertinent to phase transitions and their dynamic interactions. On the methodology aspect, this paper details a simplified temperature-based model to animate the phase transitions and their dynamic interactions, including melting, freezing, and vaporization, by integrating the latent heat model with relevant governing physical laws. On the numerical aspect, our framework supports a new algorithm aiming at tight coupling of heat transfer and multiphase FLIP-based fluids. Specifically for liquid-gas phase transition, we take into account the dissolved gas involved in liquid which further enhances the bubble generation effects. Besides the unique feature of heat transfer, we also devise a SPH-FLIP coupled model to simulate sub-grid bubbles, which enables three-phase dynamic interactions among solid, liquid, and gas. The extensive experiments show that our hybrid approach can simultaneously handle multi-phase transition driven by physics-based heat conditions, as well as the multi-phase dynamic interactions with high fidelity and visual appeal.     

Circuit model for microfluidic bubble generation under controlled pressure

We explore the microfluidic generation of bubbles in a flow-focusing junction using a pressure-controlled device rather than the more common flow rate-controlled devices. This device is a prototype for extending microfluidic drop generation methods to molten polymers. We show that the bubble generation process is highly sensitive to pressure: small changes in pressure induce large changes in bubble size and bubble formation frequency. A simple resistance circuit model can explain this pressure dependence. Briefly, we show that bubble generation is possible only within a finite pressure range. Near the ends of this pressure range, the ratio of the flow rates of the dispersed to continuous phase is highly sensitive to pressure, and therefore so also is the bubble generation process. The circuit model offers a way to use existing models of drop generation (which are based on flow rate-controlled operation) to predict pressure-controlled operation. We also examine drop formation using a highly viscous polymer as the dispersed phase. Drops are formed far downstream of the flow-focusing junction, and they are far smaller than the microfluidic channel dimensions. These results suggest that existing microfluidic drop generation methods may be exploited to make complex particles from thermoplastic polymers.     

Lagrangian particle model for multiphase flows

A Lagrangian particle model for multiphase multicomponent fluid flow, based on smoothed particle hydrodynamics (SPH), was developed and used to simulate the flow of an emulsion consisting of bubbles of a non-wetting liquid surrounded by a wetting liquid. In SPH simulations, fluids are represented by sets of particles that are used as discretization points to solve the Navier-Stokes fluid dynamics equations. In the multiphase multicomponent SPH model, a modified van der Waals equation of state is used to close the system of flow equations. The combination of the momentum conservation equation with the van der Waals equation of state results in a particle equation of motion in which the total force acting on each particle consists of many-body repulsive and viscous forces, two-body (particle-particle) attractive forces, and body forces such as gravitational forces. Similar to molecular dynamics, for a given fluid component the combination of repulsive and attractive forces causes phase separation. The surface tension at liquid-liquid interfaces is imposed through component dependent attractive forces. The wetting behavior of the fluids is controlled by phase dependent attractive interactions between the fluid particles and stationary particles that represent the solid phase. The dynamics of fluids away from the interface is governed by purely hydrodynamic forces. Comparison with analytical solutions for static conditions and relatively simple flows demonstrates the accuracy of the SPH model.     

Particle‐based simulation of bubbles in water–solid interaction

In this paper, a particle‐based multiphase method for creating realistic animations of bubbles in water–solid interaction is presented. To generate bubbles from gas dissolved in the water on the fly, we propose an approximate model for the creation of bubbles, which takes into account the influence of gas concentration in the water, the solid material, and water–solid velocity difference. As the air particle on the bubble surface is treated as a virtual nucleation site, the bubble absorbs air from surrounding water and grows. The density and pressure forces of air bubbles are computed separately using smoothed particle hydrodynamics; then, the two‐way coupling of bubbles with water and solid is solved by a new drag force, so the generated bubbles’ flow on the surface of solid and the deformation in the rising process can be simulated. Additionally, touching bubbles merge together under the cohesion forces weighted by the smoothing kernel and velocity difference. The experimental results show that this method is capable of simulating bubbles in water–solid interaction under different physical conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Boundary Detection in Particle‐based Fluids

This paper presents a novel method to detect free‐surfaces on particle‐based volume representation. In contrast to most particle‐based free‐surface detection methods, which perform the surface identification based on physical and geometrical properties derived from the underlying fluid flow simulation, the proposed approach only demands the spatial location of the particles to properly recognize surface particles, avoiding even the use of kernels. Boundary particles are identified through a Hidden Point Removal (HPR) operator used for visibility test. Our method is very simple, fast, easy to implement and robust to changes in the distribution of particles, even when facing large deformation of the free‐surface. A set of comparisons against state‐of‐the‐art boundary detection methods show the effectiveness of our approach. The good performance of our method is also attested in the context of fluid flow simulation involving free‐surface, mainly when using level‐sets for rendering purposes.     

Lattice-based flow field modeling

We present an approach for simulating the natural dynamics that emerge from the interaction between a flow field and immersed objects. We model the flow field using the lattice Boltzmann model (LBM) with boundary conditions appropriate for moving objects and accelerate the computation on commodity graphics hardware (GPU) to achieve real-time performance. The boundary conditions mediate the exchange of momentum between the flow field and the moving objects resulting in forces exerted by the flow on the objects as well as the back-coupling on the flow. We demonstrate our approach using soap bubbles and a feather. The soap bubbles illustrate Fresnel reflection, reveal the dynamics of the unseen flow field in which they travel, and display spherical harmonics in their undulations. Our simulation allows the user to directly interact with the flow field to influence the dynamics in real time. The free feather flutters and gyrates in response to lift and drag forces created by its motion relative to the flow. Vortices are created as the free feather falls in an otherwise quiescent flow.     

Biochip with integrated pumps for plug-based sequential exchange of solutions

A novel method for exchanging solutions used in biochemical analyses and a device to carry out the exchange are proposed. An array of plugs formed using six injectors was transported in a microflow channel using a main pump located at one end of the main flow channel. The injectors and main pump were operated on the basis of the change in volume caused by the electrolysis of water. Bubbles were produced from working electrodes; these bubbles caused a diaphragm placed below the injectors to inflate and occlude the inlet of the solution reservoir. Increase in the number of bubbles caused the reservoir to inject the solution into the main flow channel in the form of a plug. Each plug was individually transported downstream to the sensing area by the main pump, which was operated in a similar manner to the injector. The device was used for the detection of a tumor marker, α-fetoprotein (AFP). Plugs of necessary solutions were individually transported to the sensing area with immobilized primary antibodies to allow antigen–antibody binding, cleaning, and detection. The fluorescence intensity from the antibodies showed clear dependence on the concentration of AFP. The immobilization of antibodies could also be carried out on-chip.     

A Divergence-free Mixture Model for Multiphase Fluids

We present a novel divergence free mixture model for multiphase flows and the related fluid-solid coupling. The new mixture model is built upon a volume-weighted mixture velocity so that the divergence free condition is satisfied for miscible and immiscible multiphase fluids. The proposed mixture velocity can be solved efficiently by adapted single phase incompressible solvers, allowing for larger time steps and smaller volume deviations. Besides, the drift velocity formulation is corrected to ensure mass conservation during the simulation. The new approach increases the accuracy of multiphase fluid simulation by several orders. The capability of the new divergence-free mixture model is demonstrated by simulating different multiphase flow phenomena including mixing and unmixing of multiple fluids, fluid-solid coupling involving deformable solids and granular materials.     

Distortion-Free Steganography for Polygonal Meshes

We present a technique for steganography in polygonal meshes. Our method hides a message in the indexed rep‐resentation of a mesh by permuting the order in which faces and vertices are stored. The permutation is relative to a reference ordering that encoder and decoder derive from the mesh connectivity in a consistent manner. Our method is distortion‐free because it does not modify the geometry of the mesh. Compared to previous steganographic methods for polygonal meshes our capacity is up to an order of magnitude better. Our steganography algorithm is universal and can be used instead of the standard permutation steganography algorithm on arbitrary datasets. The standard algorithm runs in Ω (n² log² n log log n) time and achieves optimal O(nlog n) bit capacity on datasets with n elements. In contrast, our algorithm runs in O(n) time, achieves a capacity that is only one bit per element less than optimal, and is extremely simple to implement.     

Shadowing Dynamic Scenes with Arbitrary BRDFs

We present a real-time relighting and shadowing method for dynamic scenes with varying lighting, view and BRDFs. Our approach is based on a compact representation of reflectance data that allows for changing the BRDF at run-time and a data-driven method for accurately synthesizing self-shadows on articulated and deformable geometries. Unlike previous self-shadowing approaches, we do not rely on local blocking heuristics. We do not fit a model to the BRDF-weighted visibility, but rather only to the visibility that changes during animation. In this manner, our model is more compact than previous techniques and requires less computation both during fitting and at run-time. Our reflectance product operators can re-integrate arbitrary low-frequency view-dependent BRDF effects on-the-fly and are compatible with all previous dynamic visibility generation techniques as well as our own data-driven visibility model. We apply our reflectance product operators to three different visibility generation models, and our data-driven model can achieve framerates well over 300Hz.     

Experimental and computational study of gas bubble removal in a microfluidic system using nanofibrous membranes

This paper presents a simple and efficient method for removing gas bubbles from a microfluidic system. This bubble removal system uses a T-junction configuration to generate gas bubbles within a water-filled microchannel. The generated bubbles are then transported to a bubble removal region and vented through a hydrophobic nanofibrous membrane. Four different hydrophobic Polytetrafluorethylene membranes with different pore sizes ranging from 0.45 to 3 μm are tested to study the effect of membrane structure on the system performance. The fluidic channel width is 500 μm and channel height ranges from 100 to 300 μm. Additionally, a 3D computational fluid dynamics model is developed to simulate the bubble generation and its removal from a microfluidic system. Computational results are found to be in a good agreement with the experimental data. The effects of various geometrical and flow parameters on bubble removal capability of the system are studied. Furthermore, gas–liquid two-phase flow behaviors for both the complete and partial bubble removal cases are thoroughly investigated. The results indicate that the gas bubble removal rate increases with increasing the pore size and channel height but decreases with increasing the liquid flow rate.

     

High-Resolution Volumetric Computation of Offset Surfaces with Feature Preservation

We present a new algorithm for the efficient and reliable generation of offset surfaces for polygonal meshes. The algorithm is robust with respect to degenerate configurations and computes (self‐)intersection free offsets that do not miss small and thin components. The results are correct within a prescribed ε‐tolerance. This is achieved by using a volumetric approach where the offset surface is defined as the union of a set of spheres, cylinders, and prisms instead of surface‐based approaches that generally construct an offset surface by shifting the input mesh in normal direction. Since we are using the unsigned distance field, we can handle any type of topological inconsistencies including non‐manifold configurations and degenerate triangles. A simple but effective mesh operation allows us to detect and include sharp features (shocks) into the output mesh and to preserve them during post‐processing (decimation and smoothing). We discretize the distance function by an efficient multi‐level scheme on an adaptive octree data structure. The problem of limited voxel resolutions inherent to every volumetric approach is avoided by breaking the bounding volume into smaller tiles and processing them independently. This allows for almost arbitrarily high voxel resolutions on a commodity PC while keeping the output mesh complexity low. The quality and performance of our algorithm is demonstrated for a number of challenging examples.     

整流型加油枪在柴油加注过程中的数值模拟与试验研究

应用计算流体力学软件Fluent,基于多相流VOF模型,进行加油枪加油时有、无整流管两种状态的模拟,得出的油箱内油气体积分数云图表明：加装整流管加注时油箱内流场更加稳定,产生的气泡更少。通过柴油加注试验,对比分析加油枪有、无整流管两种情况,得到不同加油阶段油液状态的图像,由此表明整流管有较好的整流效果,可有效减少油液面上油沫的产生。最后总结研究结论并对整流管的使用提出了建议。     

Evolving Self-Organizing Behaviors for a Swarm-Bot

In this paper, we introduce a self-assembling and self-organizing artifact, called a swarm-bot, composed of a swarm of s-bots, mobile robots with the ability to connect to and to disconnect from each other. We discuss the challenges involved in controlling a swarm-bot and address the problem of synthesizing controllers for the swarm-bot using artificial evolution. Specifically, we study aggregation and coordinated motion of the swarm-bot using a physics-based simulation of the system. Experiments, using a simplified simulation model of the s-bots, show that evolution can discover simple but effective controllers for both the aggregation and the coordinated motion of the swarm-bot. Analysis of the evolved controllers shows that they have properties of scalability, that is, they continue to be effective for larger group sizes, and of generality, that is, they produce similar behaviors for configurations different from those they were originally evolved for. The portability of the evolved controllers to real s-bots is tested using a detailed simulation model which has been validated against the real s-bots in a companion paper in this same special issue.     

Error handling strategies in multiphase inverse modeling

Parameter estimation by inverse modeling involves the repeated evaluation of a function of residuals. These residuals represent both errors in the model and errors in the data. In practical applications of inverse modeling of multiphase flow and transport, the error structure of the final residuals often significantly deviates from the statistical assumptions that underlie standard maximum likelihood estimation using the least-squares method. Large random or systematic errors are likely to lead to convergence problems, biased parameter estimates, misleading uncertainty measures, or poor predictive capabilities of the calibrated model. The multiphase inverse modeling code iTOUGH2 supports strategies that identify and mitigate the impact of systematic or non-normal error structures. We discuss these approaches and provide an overview of the error handling features implemented in iTOUGH2.     

A quantitative sub-grid air entrainment model for bubbly flows - plunging jets

The accurate prediction of air entrainment is critical in simulating various important multiphase (air/water) flows. In this paper, we present a sub-grid air entrainment model that quantitatively predicts the rate of air entrainment and subsequent disperse bubbly flow for a plunging jet. The derivation of this model is based on the two-stage (i.e., low and high liquid jet velocity) air entrainment mechanisms suggested by Sene [Sene KJ. Air entrainment by plunging jets. Chem Eng Sci 1988;43(10):2615-23]. This model was validated against extensive experimental data for water jets in air over a wide range of liquid velocities (from around 1 to 10 m/s) for the total rate of air entrainment. It was then implemented into an Eulerian/Eulerian two-fluid computational multiphase fluid dynamics (CMFD) model, wherein the liquid and the bubbles are modeled as two distinct continua. This multiphase model, supplemented by the new sub-grid air entrainment model, was used to predict the void fraction distribution underneath plunging water jets at different depths and water jet velocities. It was found that this approach yields results that match the experimental observations very well.     

GPU-based Fast Ray Casting for a Large Number of Metaballs

Metaballs are implicit surfaces widely used to model curved objects, represented by the isosurface of a density field defined by a set of points. Recently, the results of particle‐based simulations have been often visualized using a large number of metaballs, however, such visualizations have high rendering costs. In this paper we propose a fast technique for rendering metaballs on the GPU. Instead of using polygonization, the isosurface is directly evaluated in a per‐pixel manner. For such evaluation, all metaballs contributing to the isosurface need to be extracted along each viewing ray, on the limited memory of GPUs. We handle this by keeping a list of metaballs contributing to the isosurface and efficiently update it. Our method neither requires expensive precomputation nor acceleration data structures often used in existing ray tracing techniques. With several optimizations, we can display a large number of moving metaballs quickly.     

CFD simulations of wave-current-body interactions including greenwater and wet deck slamming

A numerical method for the solution of unsteady Navier-Stokes equations has been employed in conjunction with an interface-preserving level-set method for the simulation of greenwater effect on offshore structures and ships. In this method, the free surface flows are modeled as immiscible air-water two-phase flows and the free surface itself is represented by the zero level-set function. The Navier-Stokes equations for both the water and air flows are formulated in moving curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered multi-block grid system. Large eddy simulation (LES) approach is used with Smagorinsky model to account for the effects of turbulence induced by violent free surface motions. A chimera domain decomposition approach is implemented using overlapping, embedding, or matching grids to facilitate the simulation of complex flow around practical configurations. The overset grid system also greatly simplified the simulation of arbitrary translational and rotational motions among various computational blocks. Calculations were performed first for dam-breaking flow and free jet problems involving violent free surface motions. The level-set Navier-Stokes method was then employed for the simulation of slamming of a hemisphere, greenwater on offshore structure and ships, and wet deck slamming of an X-Craft in pitch and heave motions. The numerical results clearly demonstrated the capability of the level-set method to deal with violent free surface flows involving breaking waves, water droplets, trapped air bubbles, and wave-current-body interactions.     

Modeling shear-induced diffusion force in particulate flows

In this paper, a novel approach is presented to mechanistically model the dispersion of solid particles due to shear-induced diffusion effect. The model is based on the interfacial force concept, so that it is compatible with the ensemble averaged multifield model of multicomponent flows. Although the new model has been developed with membrane filtration processes in mind, the proposed shear-induced diffusion force formulation can be utilized to model a variety of industrial processes where particulates are present. The model has been verified against the other phenomenological models currently in use.The first part of this paper is concerned with theoretical aspects of model derivation. Then, a numerical analysis is presented to illustrate the application of the new model to dilute liquid/particle two-phase flows in tube membrane systems similar to those that are used for micro- and nano-filtration processes. The numerical simulations have been performed using a state of the art multiphase/multicomponent CFD code, NPHASE [Antal SP, Ettorre SM, Kunz RF, Podowski MZ. Development of a next generation computer code for the prediction of multicomponent multiphase flows. In: International meeting on trends in numerical and physical modeling for industrial multiphase flow, Cargese, France; 2000; Kunz RF, Yu WS, Antal SP, Ettorre SM. An unstructured two-fluid method based on the coupled phasic exchange algorithm. AIAA Paper 2001; 2001:2672]. The numerical consistency of the results of computer calculations have been verified against simplified semi-analytical approximations.