Operator Splitting for Navier-Stokes and Chorin-Marsden Product Formula

Author(s):  
J. Thomas Beale ◽  
Claude Greengard ◽  
Enrique Thomann
Author(s):  
Guillermo E. Ovando ◽  
Juan C. Prince ◽  
Sandy L. Ovando

Fluid dynamics for a Newtonian fluid in the absence of body forces in a two-dimensional cavity with top and bottom curved walls was studied numerically. The vertical walls are fixed and the curved walls are in motion. The Navier-Stokes equations were solved using the finite element method combined with the operator splitting scheme. We analyzed the behaviour of the velocity fields, the vorticity fields and the velocity profiles of the fluid inside the cavity. The analysis was carried out for two different Reynolds numbers of 50 and 500 with two ratios (R = 1, −1) of the top to the bottom curved lid speed. For these values of parameters the flow is characterized by vortex formation inside the cavity. The spatial symmetry on the flow patterns are also investigated. We found that when the velocities of the top and bottom walls have opposite direction only one cell is formed in the central part of the cavity; however when the velocities of the top and bottom walls have the same direction the vortex formation inside the cavity is more complex.


2013 ◽  
Vol 2013 ◽  
pp. 1-11 ◽  
Author(s):  
Jürgen Geiser

We present some operator splitting methods improved by the use of the Zassenhaus product and designed for applications to multiphysics problems. We treat iterative splitting methods that can be improved by means of the Zassenhaus product formula, which is a sequential splitting scheme. The main idea for reducing the computation time needed by the iterative scheme is to embed fast and cheap Zassenhaus product schemes, since the computation of the commutators involved is very cheap, since we are dealing with nilpotent matrices. We discuss the coupling ideas of iterative and sequential splitting techniques and their convergence. While the iterative splitting schemes converge slowly in their first iterative steps, we improve the initial convergence rates by embedding the Zassenhaus product formula. The applications are to multiphysics problems in fluid dynamics. We consider phase models in computational fluid dynamics and analyse how to obtain higher order operator splitting methods based on the Zassenhaus product. The computational benefits derive from the use of sparse matrices, which arise from the spatial discretisation of the underlying partial differential equations. Since the Zassenhaus formula requires nearly constant CPU time due to its sparse commutators, we have accelerated the iterative splitting schemes.


2011 ◽  
Vol 217 (9) ◽  
pp. 4557-4575 ◽  
Author(s):  
Jürgen Geiser ◽  
Gamze Tanoglu

2012 ◽  
Vol 2012 ◽  
pp. 1-17 ◽  
Author(s):  
K. Poochinapan

The aim of this paper is to study the properties of approximations to nonlinear terms of the 2D incompressible Navier-Stokes equations in the stream function formulation (time-dependent biharmonic equation). The nonlinear convective terms are numerically solved by using the method with internal iterations, compared to the ones which are solved by using explicit and implicit schemes (operator splitting scheme Christov and Marinova; (2001)). Using schemes and algorithms, the steady 2D incompressible flow in a lid-driven cavity is solved up to Reynolds number Re =5000 with second-order spatial accuracy. The schemes are thoroughly validated on grids with different resolutions. The result of numerical experiments shows that the finite difference scheme with internal iterations on nonlinearity is more efficient for the high Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document