https://www.cs.ubc.ca/~mitchell/Papers/myJCP02.pdf

A hybrid particle level set method for improved interface capturing

For more information on the USGS—the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment, visit http://www.usgs.gov or call 1–888–ASK–USGS. For an overview of USGS information products, including maps, imagery, and publications, Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Although this report is in the public domain, permission must be secured from the individual copyright owners to reproduce any copyrighted materials contained within this report. 10. Diagram showing the GMTED2010 layer extents (minimum and maximum latitude and longitude) are a result of the coordinate system inherited from the 1-arc-second SRTM

Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010)

The growth of a vapor bubble in a superheated liquid is controlled by three factors: the inertia of the liquid, the surface tension, and the vapor pressure. As the bubble grows, evaporation takes place at the bubble boundary, and the temperature and vapor pressure in the bubble are thereby decreased. The heat inflow requirement of evaporation, however, depends on the rate of bubble growth, so that the dynamic problem is linked with a heat diffusion problem. Since the heat diffusion problem has been solved, a quantitative formulation of the dynamic problem can be given. A solution for the radius of the vapor bubble as a function of time is obtained which is valid for sufficiently large radius. This asymptotic solution covers the range of physical interest since the radius at which it becomes valid is near the lower limit of experimental observation. It shows the strong effect of heat diffusion on the rate of bubble growth. Comparison of the predicted radius‐time behavior is made with experimental observations in superheated water, and very good agreement is found.

The growth of vapor bubbles in superheated liquids. report no. 26-6

In the present work experimental, numerical, and theoretical investigations of a normal drop impact onto a liquid film of finite thickness are presented. The dynamics of drop impact on liquid surfaces, the shape of the cavity, the formation and propagation of a capillary wave in the crater, and the residual film thickness on the rigid wall are determined and analyzed. The shape of the crater within the film and the uprising liquid sheet formed upon the impact are observed using a high-speed video system. The effects of various influencing parameters such as drop impact velocity, liquid film thickness and physical properties of the liquids, including viscosity and surface tension, on the time evolution of the crater formation are investigated. Complementary to experiments the direct numerical simulations of the phenomena are performed using an advanced free-surface capturing model based on a two-fluid formulation of the classical volume-of-fluid (VOF) model in the framework of the finite volume numerical method. In this model an additional convective term is introduced into the transport equation for phase fraction, contributing decisively to a sharper interface resolution. Furthermore, an analytical model for the penetration depth of the crater is developed accounting for the liquid inertia, viscosity, gravity, and surface tension. The model agrees well with the experiments at the early times of penetration far from the wall if the impact velocity is high. Finally, a scaling analysis of the residual film thickness on the wall is conducted demonstrating a good agreement with the numerical predictions.

Drop impact onto a liquid layer of finite thickness: dynamics of the cavity evolution.

We devise a numerical method for passive advection of a surface, such as the interface between two incompressible fluids, across a computational mesh. The method is called isoAdvector, and is developed for general meshes consisting of arbitrary polyhedral cells. The algorithm is based on the volume of fluid (VOF) idea of calculating the volume of one of the fluids transported across the mesh faces during a time step. The novelty of the isoAdvector concept consists of two parts. First, we exploit an isosurface concept for modelling the interface inside cells in a geometric surface reconstruction step. Second, from the reconstructed surface, we model the motion of the face-interface intersection line for a general polygonal face to obtain the time evolution within a time step of the submerged face area. Integrating this submerged area over the time step leads to an accurate estimate for the total volume of fluid transported across the face. The method was tested on simple two-dimensional and three-dimensional interface advection problems on both structured and unstructured meshes. The results are very satisfactory in terms of volume conservation, boundedness, surface sharpness and efficiency. The isoAdvector method was implemented as an OpenFOAM® extension and is published as open source.

/pdf/a-computational-method-for-sharp-interface-advection-2fzz651ohr.pdf

A computational method for sharp interface advection

Coupling of the interface tracking and the two-fluid models for the simulation of incompressible two-phase flow

Modelling of Two-Phase Flow with Second-Order Accurate Scheme

Microsimulation of electric vehicle energy consumption

One of the important limitations of the interface tracking algorithms is that they can be used only as long as the local computational grid density allows surface tracking. In a dispersed flow, where the dimensions of the particular fluid parts are comparable or smaller than the grid spacing, several numerical and reconstruction errors become considerable. In this paper the analysis of the interface tracking errors is performed for the volume-of-fluid method with the least squares volume of fluid interface reconstruction algorithm. A few simple two-fluid benchmarks are proposed for the investigation of the interface tracking grid dependence. The expression based on the gradient of the volume fraction variable is introduced for the estimation of the reconstruction correctness and can be used for the activation of an adaptive mesh refinement algorithm. Copyright © 2002 John Wiley & Sons, Ltd.

Stojan Petelin

Papers

Coupling of the interface tracking and the two-fluid models for the simulation of incompressible two-phase flow

Modelling of Two-Phase Flow with Second-Order Accurate Scheme

Microsimulation of electric vehicle energy consumption

Numerical errors of the volume‐of‐fluid interface tracking algorithm

Modeling the impact of road network configuration on vehicle energy consumption