Implementation of a Semi-Implicit Pressure-Based Multigrid Fluid Flow Algorithm on a Graphics Processing Unit

2009 ◽  
Author(s):  
Aaron F. Shinn ◽  
S. P. Vanka
Keyword(s):  
Stokes Equations ◽  
Graphics Processing Unit ◽  
Navier Stokes ◽  
Processing Unit ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Multigrid Algorithm ◽  
Computational Speed ◽  
Speed Up ◽  
Graphics Processing

A semi-implicit pressure based multigrid algorithm for solving the incompressible Navier-Stokes equations was implemented on a Graphics Processing Unit (GPU) using CUDA (Compute Unified Device Architecture). The multigrid method employed was the Full Approximation Scheme (FAS), which is used for solving nonlinear equations. This algorithm is applied to the 2D driven cavity problem and compared to the CPU version of the code (written in Fortran) to assess computational speed-up.

Fast Fluid Thermodynamics Simulation by Solving Heat Diffusion Equation

2021 ◽  
Vol 11 (04) ◽  
pp. 1-11
Author(s):  
Wanwan Li
Keyword(s):  
Diffusion Equation ◽  
Real Time ◽  
Stokes Equations ◽  
Fluid Simulation ◽  
Heat Diffusion ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Heat Diffusion Equation ◽  
Acceleration Techniques ◽  
Speed Up

In mechanical engineering educations, simulating fluid thermodynamics is rather helpful for students to understand the fluid’s natural behaviors. However, rendering both high-quality and realtime simulations for fluid dynamics are rather challenging tasks due to their intensive computations. So, in order to speed up the simulations, we have taken advantage of GPU acceleration techniques to simulate interactive fluid thermodynamics in real-time. In this paper, we present an elegant, basic, but practical OpenGL/SL framework for fluid simulation with a heat map rendering. By solving Navier-Stokes equations coupled with the heat diffusion equation, we validate our framework through some real-case studies of the smoke-like fluid rendering such as their interactions with moving obstacles and their heat diffusion effects. As shown in Fig. 1, a group of experimental results demonstrates that our GPU-accelerated solver of Navier-Stokes equations with heat transfer could give the observers impressive real-time and realistic rendering results.

scholarly journals Finite element method completely implemented for graphic processor units using parallel algorithm libraries

2017 ◽  
Vol 33 (1) ◽  
pp. 53-66 ◽  
Author(s):  
Franz Pichler ◽  
Gundolf Haase
Keyword(s):  
Finite Element ◽  
Graphics Processing Unit ◽  
Computational Cost ◽  
Processing Unit ◽  
Time Step ◽  
Device Architecture ◽  
Transient Problems ◽  
Speed Up ◽  
Automotive Batteries ◽  
Graphics Processing

A finite element code is developed in which all of the computationally expensive steps are performed on a graphics processing unit via the THRUST and the PARALUTION libraries. The code focuses on the simulation of transient problems where the repeated computations per time-step create the computational cost. It is used to solve partial and ordinary differential equations as they arise in thermal-runaway simulations of automotive batteries. The speed-up obtained by utilizing the graphics processing unit for every critical step is compared against the single core and the multi-threading solutions which are also supported by the chosen libraries. This way a high total speed-up on the graphics processing unit is achieved without the need for programming a single classical Compute Unified Device Architecture kernel.

scholarly journals On the Importance of the Influence of both Velocity and Aspect Ratios on the Occurrence of Bifurcation Phenomena within a Two-Sided Lid-Driven Enclosure

2015 ◽  
Vol 55 ◽  
pp. 160-172
Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Numerical Results ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Left Wall ◽  
Driven Cavity ◽  
Aspect Ratios ◽  
The Stability ◽  
The Right

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.

scholarly journals Numerical Simulation of Slip effect on Lid-Driven Cavity Flow Problem for High Reynolds Number: Vorticity – Stream Function Approach

2021 ◽  
Vol 8 (3) ◽  
pp. 418-424
Author(s):  
Syed Fazuruddin ◽  
Seelam Sreekanth ◽  
G. Sankara Sekhar Raju
Keyword(s):  
Reynolds Number ◽  
Stream Function ◽  
Stokes Equations ◽  
Cavity Flow ◽  
Partial Slip ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Slip Conditions ◽  
Driven Cavity Flow

Incompressible 2-D Navier-stokes equations for various values of Reynolds number with and without partial slip conditions are studied numerically. The Lid-Driven cavity (LDC) with uniform driven lid problem is employed with vorticity - Stream function (VSF) approach. The uniform mesh grid is used in finite difference approximation for solving the governing Navier-stokes equations and developed MATLAB code. The numerical method is validated with benchmark results. The present work is focused on the analysis of lid driven cavity flow of incompressible fluid with partial slip conditions (imposed on side walls of the cavity). The fluid flow patterns are studied with wide range of Reynolds number and slip parameters.

A Simplified Parallel Two-Level Iterative Method for Simulation of Incompressible Navier-Stokes Equations

2015 ◽  
Vol 7 (6) ◽  
pp. 715-735
Author(s):  
Yueqiang Shang ◽  
Jin Qin
Keyword(s):  
Stokes Equations ◽  
Stokes Problem ◽  
Coarse Grid ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Step Flow ◽  
Grid Solution ◽  
Oseen Problem ◽  
Parallel Iterative

AbstractBased on two-grid discretization, a simplified parallel iterative finite element method for the simulation of incompressible Navier-Stokes equations is developed and analyzed. The method is based on a fixed point iteration for the equations on a coarse grid, where a Stokes problem is solved at each iteration. Then, on overlapped local fine grids, corrections are calculated in parallel by solving an Oseen problem in which the fixed convection is given by the coarse grid solution. Error bounds of the approximate solution are derived. Numerical results on examples of known analytical solutions, lid-driven cavity flow and backward-facing step flow are also given to demonstrate the effectiveness of the method.

Fekete-Gauss Spectral Elements for Incompressible Navier-Stokes Flows: The Two-Dimensional Case

2013 ◽  
Vol 13 (5) ◽  
pp. 1309-1329 ◽  
Author(s):  
Laura Lazar ◽  
Richard Pasquetti ◽  
Francesca Rapetti
Keyword(s):  
Degrees Of Freedom ◽  
Stokes Equations ◽  
Schur Complement ◽  
Spectral Element ◽  
Navier Stokes ◽  
Two Dimensional ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Flow Past A Cylinder ◽  
Accuracy Study

AbstractSpectral element methods on simplicial meshes, say TSEM, show both the advantages of spectral and finite element methods, i.e., spectral accuracy and geometrical flexibility. We present aTSEM solver of the two-dimensional (2D) incompressible Navier-Stokes equations, with possible extension to the 3D case. It uses a projection method in time and piecewise polynomial basis functions of arbitrary degree in space. The so-called Fekete-Gauss TSEM is employed,i.e., Fekete (resp. Gauss) points of the triangle are used as interpolation (resp. quadrature) points. For the sake of consistency, isoparametric elements are used to approximate curved geometries. The resolution algorithm is based on an efficient Schur complement method, so that one only solves for the element boundary nodes. Moreover, the algebraic system is never assembled, therefore the number of degrees of freedom is not limiting. An accuracy study is carried out and results are provided for classical benchmarks: the driven cavity flow, the flow between eccentric cylinders and the flow past a cylinder.

scholarly journals The use of dual reciprocity method for 2D laminar viscous flow

2020 ◽  
Vol 310 ◽  
pp. 00044
Author(s):  
Juraj Mužík
Keyword(s):  
Viscous Flow ◽  
Flow Problem ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Singular Boundary ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Boundary Method ◽  
2D Flow ◽  
Singular Boundary Method

The paper presents the use of the dual reciprocity multidomain singular boundary method (SBMDR) for the solution of the laminar viscous flow problem described by Navier-Stokes equations. A homogeneous part of the solution is solved using a singular boundary method with the 2D Stokes fundamental solution - Stokeslet. The dual reciprocity approach has been chosen because it is ideal for the treatment of the nonhomogeneous and nonlinear terms of Navier-Stokes equations. The presented SBMDR approach to the solution of the 2D flow problem is demonstrated on a standard benchmark problem - lid-driven cavity.

Flow Pattern on the Surface of Airfoil With Lower Reynolds Number

2009 ◽  
Author(s):  
Xu Sun ◽  
Jia-Zhong Zhang
Keyword(s):  
Reynolds Number ◽  
Flow Pattern ◽  
Stokes Equations ◽  
Separation Bubble ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Naca0012 Airfoil ◽  
Lower Reynolds Number ◽  
Characteristic Based Split

In this paper, aerodynamic performance of the NACA0012 airfoil in the incompressible flow with a lower Reynolds number (Re) is investigated numerically from the viewpoints of flow pattern and nonlinear dynamics. First, the characteristic-based split (CBS) finite element method is introduced for the approximation of the incompressible Navier-Stokes equations, and then the lid-driven cavity flow and flow around a circular cylinder are calculated for varification. Then, at Re = 1000, flow fields around the NACA0012 airfoil at a series of angles of attack are simulated. With the increase of the attack angle, great change of the flow pattern appears, and the flow structures such as trailing edge vortex, separation bubble and shedding vortex are observed. Moreover, it is found that the separation bubble plays an important role in the deterioration of the flow stability at higher attack angles, and the vortex shedding can be taken as the result of a Hopf bifurcation while the bifurcation parameter is the angle of attack.

Sign in / Sign up

Export Citation Format

Share Document