scholarly journals Spectral Direction Splitting Schemes for the Incompressible Navier-Stokes Equations

2011 ◽  
Vol 1 (3) ◽  
pp. 215-234 ◽  
Author(s):  
Lizhen Chen ◽  
Jie Shen ◽  
Chuanju Xu

AbstractWe propose and analyze spectral direction splitting schemes for the incompressible Navier-Stokes equations. The schemes combine a Legendre-spectral method for the spatial discretization and a pressure-stabilization/direction splitting scheme for the temporal discretization, leading to a sequence of one-dimensional elliptic equations at each time step while preserving the same order of accuracy as the usual pressure-stabilization schemes. We prove that these schemes are unconditionally stable, and present numerical results which demonstrate the stability, accuracy, and efficiency of the proposed methods.

2014 ◽  
Vol 16 (5) ◽  
pp. 1239-1262 ◽  
Author(s):  
Feng Shi ◽  
Guoping Liang ◽  
Yubo Zhao ◽  
Jun Zou

AbstractWe present a new splitting method for time-dependent convention-dominated diffusion problems. The original convention diffusion system is split into two sub-systems: a pure convection system and a diffusion system. At each time step, a convection problem and a diffusion problem are solved successively. A few important features of the scheme lie in the facts that the convection subproblem is solved explicitly and multistep techniques can be used to essentially enlarge the stability region so that the resulting scheme behaves like an unconditionally stable scheme; while the diffusion subproblem is always self-adjoint and coercive so that they can be solved efficiently using many existing optimal preconditioned iterative solvers. The scheme can be extended for solving the Navier-Stokes equations, where the nonlinearity is resolved by a linear explicit multistep scheme at the convection step, while only a generalized Stokes problem is needed to solve at the diffusion step and the major stiffness matrix stays invariant in the time marching process. Numerical simulations are presented to demonstrate the stability, convergence and performance of the single-step and multistep variants of the new scheme.


2015 ◽  
Vol 15 (3) ◽  
pp. 307-330 ◽  
Author(s):  
Nan Jiang ◽  
Hoang Tran

AbstractIn this work, we study Crank–Nicolson leap-frog (CNLF) methods with fast-slow wave splittings for Navier–Stokes equations (NSE) with a rotation/Coriolis force term, which is a simplification of geophysical flows. We propose a new stabilized CNLF method where the added stabilization completely removes the method's CFL time step condition. A comprehensive stability and error analysis is given. We also prove that for Oseen equations with the rotation term, the unstable mode (for which $u^{n+1}+u^{n-1}\equiv 0 $) of CNLF is asymptotically stable. Numerical results are provided to verify the stability and the convergence of the methods.


Author(s):  
Alexander Danilov ◽  
Alexander Lozovskiy ◽  
Maxim Olshanskii ◽  
Yuri Vassilevski

AbstractThe paper introduces a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method is based on a quasi-Lagrangian formulation of the problem and handling the geometry in a time-explicit way. We prove that numerical solution satisfies a discrete analogue of the fundamental energy estimate. This stability estimate does not require a CFL time-step restriction. The method is further applied to simulation of a flow in a model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.


2013 ◽  
Vol 136 (2) ◽  
Author(s):  
Graham Ashcroft ◽  
Christian Frey ◽  
Kathrin Heitkamp ◽  
Christian Weckmüller

This is the first part of a series of two papers on unsteady computational fluid dynamics (CFD) methods for the numerical simulation of aerodynamic noise generation and propagation. In this part, the stability, accuracy, and efficiency of implicit Runge–Kutta schemes for the temporal integration of the compressible Navier–Stokes equations are investigated in the context of a CFD code for turbomachinery applications. Using two model academic problems, the properties of two explicit first stage, singly diagonally implicit Runge–Kutta (ESDIRK) schemes of second- and third-order accuracy are quantified and compared with more conventional second-order multistep methods. Finally, to assess the ESDIRK schemes in the context of an industrially relevant configuration, the schemes are applied to predict the tonal noise generation and transmission in a modern high bypass ratio fan stage and comparisons with the corresponding experimental data are provided.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu

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.


2013 ◽  
Vol 2013 ◽  
pp. 1-6
Author(s):  
Wen-Juan Wang ◽  
Yan Jia

We study the stability issue of the generalized 3D Navier-Stokes equations. It is shown that if the weak solutionuof the Navier-Stokes equations lies in the regular class∇u∈Lp(0,∞;Bq,∞0(ℝ3)),(2α/p)+(3/q)=2α,2<q<∞,0<α<1, then every weak solutionv(x,t)of the perturbed system converges asymptotically tou(x,t)asvt-utL2→0,t→∞.


Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.


Author(s):  
Sakir Amiroudine

The case of a supercritical fluid heated from below (Rayleigh-Bénard) in a rectangular cavity is first presented. The stability of the two boundary layers (hot and cold) is analyzed by numerically solving the Navier-Stokes equations with a van der Waals gas and stability diagrams are derived. The very large compressibility and the very low heat diffusivity of near critical pure fluids induce very large density gradients which lead to a Rayleigh–Taylor-like gravitational instability of the heat diffusion layer and results in terms of growth rates and wave numbers are presented. Depending on the relative direction of the interface or the boundary layer with respect to vibration, vibrational forces can destabilize a thermal boundary layer, resulting in parametric/Rayleigh vibrational instabilities. This has recently been achieved by using a numerical model which does not require any equation of state and directly calculates properties from NIST data base, for instance.


Author(s):  
Jean-Yves Chemin ◽  
Benoit Desjardins ◽  
Isabelle Gallagher ◽  
Emmanuel Grenier

In this chapter we intend to investigate the stability of the Leray solutions constructed in the previous chapter. It is useful to start by analyzing the linearized version of the Navier–Stokes equations, so the first section of the chapter is devoted to the proof of the well-posedness of the time-dependent Stokes system. The study will be applied in Section 3.2 to the two-dimensional Navier–Stokes equations, and the more delicate case of three space dimensions will be dealt with in Sections 3.3–3.5.


Sign in / Sign up

Export Citation Format

Share Document