A semi-implicit numerical method for the free-surface Navier-Stokes equations

2013 ◽  
Vol 74 (8) ◽  
pp. 605-622 ◽  
Author(s):  
Vincenzo Casulli
Keyword(s):  
Free Surface ◽  
Numerical Method ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Implicit Numerical Method

Level-Set Computations of Free Surface Rotational Flows

2005 ◽  
Vol 127 (6) ◽  
pp. 1111-1121 ◽  
Author(s):  
Giuseppina Colicchio ◽  
Maurizio Landrini ◽  
John R. Chaplin
Keyword(s):  
Free Surface ◽  
Numerical Method ◽  
Level Set ◽  
Stokes Equations ◽  
Vertical Surface ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Two Fluids ◽  
Level Set Function ◽  
Rotational Flows

A numerical method is developed for modeling the violent motion and fragmentation of an interface between two fluids. The flow field is described through the solution of the Navier-Stokes equations for both fluids (in this case water and air), and the interface is captured by a Level-Set function. Particular attention is given to modeling the interface, where most of the numerical approximations are made. Novel features are that the reintialization procedure has been redefined in cells crossed by the interface; the density has been smoothed across the interface using an exponential function to obtain an equally stiff variation of the density and of its inverse; and we have used an essentially non-oscillatory scheme with a limiter whose coefficients depend on the distance function at the interface. The results of the refined scheme have been compared with those of the basic scheme and with other numerical solvers, with favorable results. Besides the case of the vertical surface-piercing plate (for which new laboratory measurements were carried out) the numerical method is applied to problems involving a dam-break and wall-impact, the interaction of a vortex with a free surface, and the deformation of a cylindrical bubble. Promising agreement with other sources of data is found in every case.

scholarly journals A SEMI-IMPLICIT SCHEME FOR SOLVING INCOMPRESSIBLE VISCOUS FREE SURFACE FLOWS

2005 ◽  
Vol 4 (2) ◽  
Author(s):  
C. M. Oishi ◽  
J. A. Cuminato ◽  
V. G. Ferreira ◽  
M. F. Tomé ◽  
A. Castelo ◽  
...  
Keyword(s):  
Free Surface ◽  
Numerical Method ◽  
Stokes Equations ◽  
Free Surface Flows ◽  
Numerical Results ◽  
Navier Stokes ◽  
Variable Order ◽  
Navier Stokes Equations ◽  
Reynolds Number Flow ◽  
The Stability

The present work is concerned with a numerical method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in the primitive variables formulation. The diffusive terms are treated by Implicit Backward and Crank-Nicolson methods, and the non-linear convection terms are, explicitly, approximated by the high order upwind VONOS (Variable-Order Non-Oscillatory Scheme) scheme. The boundary conditions for the pressure field at the free surface are treated implicitly, and for the velocity field explicitly. The numerical method is then applied to the simulation of free surface and confined flows. The numerical results show that the present technique eliminates the stability restriction in the original explicit method. For low Reynolds number flow dynamics, the method is robust and produces numerical results that compare very well with the analytical solutions.

Numerical Modeling of Unsteady Flow Through a Turbomachine Stage

1998 ◽  
Author(s):  
R. I. Issa ◽  
M. A. Sadri
Keyword(s):  
Numerical Method ◽  
Periodic Boundary ◽  
Stokes Equations ◽  
Unsteady Flows ◽  
Periodic Boundary Conditions ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Sliding Mesh ◽  
Flow Through ◽  
Implicit Numerical Method

A numerical method is presented for the simulation of unsteady flows through turbomachine stages with unequal numbers of rotor and stator blades. The method solves the two-dimensional incompressible, unsteady, ensemble averaged, Navier-Stokes equations together with transport equations for the k–ϵ turbulence model employed to simulate the effects of turbulence. The method employs an implicit pressure-based finite volume discretisation procedure. In order to simulate the flow in the rotor and stator passages simultaneously, a sliding mesh methodology is developed which allows the mesh mapping the rotor domain to move in a sliding action relative to the static mesh which maps the stator passage. Phase-lagged periodic boundary conditions are implemented in the context of the implicit numerical method developed to handle unequal rotor and stator pitches efficiently. The effectiveness and accuracy of the method are assessed against data for a rotor/stator configuration with unequal pitches in adjacent rows of a low speed turbine.

scholarly journals A SEMI-IMPLICIT SCHEME FOR SOLVING INCOMPRESSIBLE VISCOUS FREE SURFACE FLOWS

2005 ◽  
Vol 4 (2) ◽  
pp. 106
Author(s):  
C. M. Oishi ◽  
J. A. Cuminato ◽  
V. G. Ferreira ◽  
M. F. Tomé ◽  
A. Castelo ◽  
...  
Keyword(s):  
Free Surface ◽  
Numerical Method ◽  
Stokes Equations ◽  
Free Surface Flows ◽  
Numerical Results ◽  
Navier Stokes ◽  
Variable Order ◽  
Navier Stokes Equations ◽  
Reynolds Number Flow ◽  
The Stability

The present work is concerned with a numerical method for solving the two-dimensional time-dependent incompressible Navier-Stokes equations in the primitive variables formulation. The diffusive terms are treated by Implicit Backward and Crank-Nicolson methods, and the non-linear convection terms are, explicitly, approximated by the high order upwind VONOS (Variable-Order Non-Oscillatory Scheme) scheme. The boundary conditions for the pressure field at the free surface are treated implicitly, and for the velocity field explicitly. The numerical method is then applied to the simulation of free surface and confined flows. The numerical results show that the present technique eliminates the stability restriction in the original explicit method. For low Reynolds number flow dynamics, the method is robust and produces numerical results that compare very well with the analytical solutions.

Simulation of the Gap Resonance Between Two Rectangular Barges in Regular Waves by a Free Surface Viscous Flow Solver

2013 ◽  
Author(s):  
B. Elie ◽  
G. Reliquet ◽  
P.-E. Guillerm ◽  
O. Thilleul ◽  
P. Ferrant ◽  
...  
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Viscous Flow ◽  
Stokes Equations ◽  
Resonance Phenomenon ◽  
Navier Stokes ◽  
Regular Waves ◽  
Navier Stokes Equations ◽  
Flow Solver ◽  
Gap Resonance

This paper compares numerical and experimental results in the study of the resonance phenomenon which appears between two side-by-side fixed barges for different sea-states. Simulations were performed using SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach and results are compared with experimental data on two fixed barges with different headings and bilges. Numerical results, obtained using the SWENSE approach, are able to predict both the frequency and the magnitude of the RAO functions.

NUMERICAL SOLUTIONS OF 2D AND 3D SLAMMING PROBLEMS

2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Type Equation ◽  
The Body ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Highly Nonlinear ◽  
2D And 3D

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.

scholarly journals Linear pulse structure and signalling in a film flow on an inclined plane

1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES
Keyword(s):  
Free Surface ◽  
Saddle Point ◽  
Convective Instability ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Numerical Results ◽  
Navier Stokes ◽  
Surface Boundary ◽  
Navier Stokes Equations ◽  
Surface Boundary Conditions

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.

Sign in / Sign up

Export Citation Format

Share Document