Numerical Analysis of a Finite Element Approximation to a Level Set Model for Free-Surface Flows

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Tomás Chacón Rebollo ◽  
Macarena Gómez Mármol ◽  
Isabel Sánchez Muñoz

Abstract In this paper, we study a finite element discretization of a Level Set Method formulation of free-surface flow. We consider an Euler semi-implicit discretization in time and a Galerkin discretization of the level set function. We regularize the density and viscosity of the flow across the interface, following the Level Set Method. We prove stability in natural norms when the viscosity and density vary from one to the other layer and optimal error estimates for smooth solutions when the layers have the same density. We present some numerical tests for academic flows.

2020 ◽  
Vol 2020 ◽  
pp. 1-18
Author(s):  
Puyang Gao

In this paper, we develop a new computational framework to investigate the sloshing free surface flow of Newtonian and non-Newtonian fluids in the rectangular tanks. We simulate the flow via a two-phase model and employ the fixed unstructured mesh in the computation to avoid the mesh distortion and reconstruction. As for the solution of Navier–Stokes equation, we utilize the SUPG finite element method based on the splitting scheme. The same order interpolation functions are then used for velocity and pressure. Moreover, the moving interface is captured via the concise level set method. We take advantage of the implicit discontinuous Galerkin method to handle the solution of level set and its reinitialization equations. A mass correction technique is also added to ensure the mass conservation property. The dam break-free surface flow is simulated firstly to demonstrate the validity of our mathematical model. In addition, the sloshing Newtonian fluid in the tank with flat and rough bottoms is considered to illustrate the feasibility and robustness of our computational scheme. Finally, the development of free surface for non-Newtonian fluid is also studied in the two tanks, and the influence of power-law index on the sloshing fluid flow is analyzed.


Author(s):  
Pablo Go´mez ◽  
Julio Herna´ndez ◽  
Joaqui´n Lo´pez ◽  
Fe´lix Faura

A numerical study of the initial stages of wave breaking processes in shallow water is presented. The waves considered are assumed to be generated by moving a piston in a two-dimensional channel, and may appear, for example, in the injection chamber of a high-pressure die casting machine under operating conditions far from the optimal. A numerical model based on a finite-difference discretization of the Navier-Stokes equations in a Cartesian grid and a second-order approximate projection method has been developed and used to carry out the simulations. The evolution of the free surface is described using a level set method, with a reinitialization procedure of the level set function which uses a local grid refinement near the free surface. The ability of different algorithms to improve mass conservation in the reinitialization step of the level set function has been tested in a time-reversed single vortex flow. The results for the breaking wave profiles show the flow characteristics after the impact of the first plunging jet onto the wave’s forward face and during the subsequent splash-up.


2018 ◽  
Vol 52 (5) ◽  
pp. 2065-2082 ◽  
Author(s):  
Erik Burman ◽  
Jonathan Ish-Horowicz ◽  
Lauri Oksanen

We consider a finite element discretization for the reconstruction of the final state of the heat equation, when the initial data is unknown, but additional data is given in a sub domain in the space time. For the discretization in space we consider standard continuous affine finite element approximation, and the time derivative is discretized using a backward differentiation. We regularize the discrete system by adding a penalty on the H2-semi-norm of the initial data, scaled with the mesh-parameter. The analysis of the method uses techniques developed in E. Burman and L. Oksanen [Numer. Math. 139 (2018) 505–528], combining discrete stability of the numerical method with sharp Carleman estimates for the physical problem, to derive optimal error estimates for the approximate solution. For the natural space time energy norm, away from t = 0, the convergence is the same as for the classical problem with known initial data, but contrary to the classical case, we do not obtain faster convergence for the L2-norm at the final time.


2014 ◽  
Vol 6 (2) ◽  
pp. 1523-1554 ◽  
Author(s):  
B. Hillebrand ◽  
C. Thieulot ◽  
T. Geenen ◽  
A. P. van den Berg ◽  
W. Spakman

Abstract. The level set method allows for tracking material surfaces in 2-D and 3-D flow modeling and is well suited for applications of multi-material flow modeling. The level set method utilizes smooth level set functions to define material interfaces, which makes the method stable and free of oscillations that are typically observed in case step-like functions parameterize interfaces. By design the level set function is a signed distance function and gives for each point in the domain the exact distance to the interface and on which side it is located. In this paper we present four benchmarks which show the validity, accuracy and simplicity of using the level set method for multi-material flow modeling. The benchmarks are simplified setups of dynamical geophysical processes such as a Rayleigh–Taylor instability, post glacial rebound, subduction and slab detachment. We also demonstrate the benefit of using the level set method for modeling a free surface with the sticky air approach. Our results show that the level set method allows for accurate material flow modeling and that the combination with the sticky air approach works well in mimicking Earth's free surface. Since the level set method tracks material interfaces instead of materials themselves, it has the advantage that the location of these interfaces is accurately known and that it represents a viable alternative to the more commonly used tracer method.


2013 ◽  
Author(s):  
Gabriel Reliquet ◽  
Aurélien Drouet ◽  
Pierre-Emmanuel Guillerm ◽  
Erwan Jacquin ◽  
Lionel Gentaz ◽  
...  

The purpose of this paper is to present combination of the SWENSE (Spectral Wave Explicit Navier-Stokes Equations – [1]) method — an original method to treat fully nonlinear wave-body interactions — and a free surface RANSE (Reynolds Averaged Navier-Stokes Equations) solver using a single-phase Level Set method to capture the interface. The idea is to be able to simulate wave-body interactions under viscous flow theory with strong deformations of the interface (wave breaking in the vicinity of the body, green water on ship decks…), while keeping the advantages of the SWENSE scheme. The SWENSE approach is based on a physical decomposition by combining incident waves described by a nonlinear spectral scheme based on potential flow theory and an adapted Navier-Stokes solver where only the diffracted part of the flow is solved, incident flow parameters seen as forcing terms. In the single-phase Level Set method [2, 3], the air phase is neglected. Thus, only the liquid phase is solved considering a fluid with uniform properties. The location of the free surface is determined by a Level Set function initialised as the signed distance. The accuracy of simulation depends essentially on the pressure scheme used to impose free surface dynamic boundary condition. Comparisons of numerical results with experimental and numerical data for US navy combatant DTMB 5415 in calm water and in head waves are presented.


2006 ◽  
Vol 18 (6) ◽  
pp. 742-747 ◽  
Author(s):  
Lan-hao Zhao ◽  
Tong-chun Li ◽  
Ling Wang ◽  
M. I. Herreros ◽  
M. Pastor

Sign in / Sign up

Export Citation Format

Share Document