scholarly journals Second-order accuracy of volume-of-fluid interface reconstruction algorithms, II: An improved constraint on the cell size

2013 ◽  
Vol 8 (1) ◽  
pp. 123-158 ◽  
Author(s):  
Elbridge Puckett
Keyword(s):  
Cell Size ◽  
Volume Of Fluid ◽  
Second Order ◽  
Reconstruction Algorithms ◽  
Fluid Interface ◽  
Order Accuracy ◽  
Interface Reconstruction

Towards 3D Volume-of-Fluid Methods Featuring Subgrid-Scale Capturing of Interface Curvature

2014 ◽  
Author(s):  
Petar Liovic
Keyword(s):  
Volume Of Fluid ◽  
Planar Interface ◽  
Interface Tracking ◽  
Reconstruction Method ◽  
Interfacial Flow ◽  
Subgrid Scale ◽  
Order Accuracy ◽  
Interface Reconstruction ◽  
Interface Curvature ◽  
3D Volume

A new interface reconstruction method for Volume of Fluid (VOF) interface tracking is presented here, based on subgrid-scale planar interface segment reconstruction (SGS-PISR). In the SGS-PISR method implemented here, the centroid of the initial single-surface interface reconstruction is shifted along that normal to enclose the correct volume. An additional step then moves the SGS plane segments laterally outwards, to ameliorate the SGS curvature by blunting the protrusion of the centroid. The SGS-PISR method results in promising tendency towards second-order accuracy and more importantly reduced interface reconstruction errors across a range of mesh resolutions, and is targeted at improving VOF performance in resolving small grid-scale details of the interface topologies in interfacial flow CFD computations.

Volume-of-Fluid Interface Reconstruction Algorithms on Next-Generation Computer Architectures

2014 ◽  
Author(s):  
Marianne M. Francois ◽  
Li-Ta Lo ◽  
Christopher Sewell
Keyword(s):  
Data Structure ◽  
High Performance ◽  
Reconstruction Algorithm ◽  
Volume Of Fluid ◽  
Interfacial Flows ◽  
Reconstruction Algorithms ◽  
Interface Reconstruction ◽  
Computer Scientists ◽  
Performance Results ◽  
Generation Computer

With the increasing heterogeneity and on-node parallelism of high-performance computing hardware, a major challenge to computational physicists is to work in close collaboration with computer scientists to develop portable and efficient algorithms and software. The objective of our work is to implement a portable code to perform interface reconstruction using NVIDIA’s Thrust library. Interface reconstruction is a technique commonly used in volume tracking methods for simulations of interfacial flows. For that, we have designed a two-dimensional mesh data structure that is easily mapped to the 1D vectors used by Thrust and at the same time is simple to work with using familiar data structures terminology (such as cell, vertices and edges). With this new data structure in place, we have implemented a recursive volume-of-fluid initialization algorithm and a standard piecewise interface reconstruction algorithm. Our interface reconstruction algorithm makes use of a table look-up to easily identify all intersection cases, as this design is efficient on parallel architectures such as GPUs. Finally, we report performance results which show that a single implementation of these algorithms can be compiled to multiple backends (specifically, multi-core CPUs, NVIDIA GPUs, and Intel Xeon Phi coprocessors), making efficient use of the available parallelism on each.

Volume of fluid interface reconstruction methods for multi-material problems

2002 ◽  
Vol 55 (2) ◽  
pp. 151-165 ◽  
Author(s):  
David J Benson
Keyword(s):  
Volume Of Fluid ◽  
Review Article ◽  
Arbitrary Lagrangian Eulerian ◽  
Fluid Interface ◽  
Reconstruction Methods ◽  
Interface Reconstruction

Volume of fluid interface reconstruction methods are used to resolve the interfaces between different materials in Eulerian and arbitrary Lagrangian Eulerian calculations. Their accuracy is critical to the overall accuracy of the calculation since the interfaces define the interactions between adjacent materials. The methods have evolved since the early 1960s, and the early criticisms of them no longer hold. In this review article, the differences between the methods and their relative strengths are reviewed, and 38 reference sources are used.

scholarly journals A New Nutation Model Of Nonrigid Earth With Ocean And Atmosphere

2000 ◽  
Vol 180 ◽  
pp. 242-247
Author(s):  
Cheng-li Huang ◽  
Wen-jing Jin ◽  
Xing-hao Liao
Keyword(s):  
Transfer Functions ◽  
Second Order ◽  
Earth Model ◽  
Surface Boundary ◽  
Complex Frequency ◽  
Order Accuracy ◽  
Complex Scalar ◽  
Motion Equations ◽  
Surface Boundary Conditions ◽  
Nutation Series

AbstractBy integrating the truncated complex scalar gravitational motion equations for an anelastic, rotating, slightly elliptical Earth, the complex frequency dependent Earth transfer functions are computed directly. Unlike the conventional method, the effects of both oceanic loads and tidal currents are included via outer surface boundary conditions, all of which are expanded to second order in ellipticity. A modified ellipticity profile in second order accuracy for the non-hydrostatic Earth is obtained from Clairaut’s equation and the PREM Earth model by adjusting both the ellipticity of the core-mantle boundary and the global dynamical ellipticity to modern observations. The effects of different Earth models, anelastic models, and ocean models are computed and compared. The atmospheric contributions to prograde annual, retrograde annual and retrograde semiannual nutation are also included as oceanic effects. Finally, a complete new nutation series of more than 340 periods, including in-phase and out-of-phase parts of longitude and obliquity terms, for a more realistic Earth, is obtained and compared with other available nutation series and observations.

scholarly journals 3-D numerical modelling of crustal polydiapirs with volume-of-fluid methods

2020 ◽  
Vol 222 (1) ◽  
pp. 474-506
Author(s):  
Aurélie Louis-Napoléon ◽  
Muriel Gerbault ◽  
Thomas Bonometti ◽  
Cédric Thieulot ◽  
Roland Martin ◽  
...  
Keyword(s):  
Volume Of Fluid ◽  
Vof Method ◽  
Earth’S Crust ◽  
Gravitational Instabilities ◽  
Interface Reconstruction ◽  
Rayleigh Taylor Instability ◽  
Earth's Crust ◽  
Two Stages ◽  
Rayleigh Benard ◽  
Earth Dynamics

SUMMARY Gravitational instabilities exert a crucial role on the Earth dynamics and in particular on its differentiation. The Earth’s crust can be considered as a multilayered fluid with different densities and viscosities, which may become unstable in particular with variations in temperature. With the specific aim to quantify crustal scale polydiapiric instabilities, we test here two codes, JADIM and OpenFOAM, which use a volume-of-fluid (VOF) method without interface reconstruction, and compare them with the geodynamics community code ASPECT, which uses a tracking algorithm based on compositional fields. The VOF method is well-known to preserve strongly deforming interfaces. Both JADIM and OpenFOAM are first tested against documented two and three-layer Rayleigh–Taylor instability configurations in 2-D and 3-D. 2-D and 3-D results show diapiric growth rates that fit the analytical theory and are found to be slightly more accurate than those obtained with ASPECT. We subsequently compare the results from VOF simulations with previously published Rayleigh–Bénard analogue and numerical experiments. We show that the VOF method is a robust method adapted to the study of diapirism and convection in the Earth’s crust, although it is not computationally as fast as ASPECT. OpenFOAM is found to run faster than, and conserve mass as well as JADIM. Finally, we provide a preliminary application to the polydiapiric dynamics of the orogenic crust of Naxos Island (Greece) at about 16 Myr, and propose a two-stages scenario of convection and diapirism. The timing and dimensions of the modelled gravitational instabilities not only corroborate previous estimates of timing and dimensions associated to the dynamics of this hot crustal domain, but also bring preliminary insight on its rheological and tectonic contexts.

Sign in / Sign up

Export Citation Format

Share Document