scholarly journals Development of a Parallel 3D Navier–Stokes Solver for Sediment Transport Calculations in Channels

Computation ◽  
2020 ◽  
Vol 8 (4) ◽  
pp. 84
Author(s):  
Gokhan Kirkil

We propose a method to parallelize a 3D incompressible Navier–Stokes solver that uses a fully implicit fractional-step method to simulate sediment transport in prismatic channels. The governing equations are transformed into generalized curvilinear coordinates on a non-staggered grid. To develop a parallel version of the code that can run on various platforms, in particular on PC clusters, it was decided to parallelize the code using Message Passing Interface (MPI) which is one of the most flexible parallel programming libraries. Code parallelization is accomplished by “message passing” whereby the computer explicitly uses library calls to accomplish communication between the individual processors of the machine (e.g., PC cluster). As a part of the parallelization effort, besides the Navier–Stokes solver, the deformable bed module used in simulations with loose beds are also parallelized. The flow, sediment transport, and bathymetry at equilibrium conditions were computed with the parallel and serial versions of the code for the case of a 140-degree curved channel bend of rectangular section. The parallel simulation conducted on eight processors gives exactly the same results as the serial solver. The parallel version of the solver showed good scalability.

2014 ◽  
Vol 493 ◽  
pp. 215-220
Author(s):  
Vivien Djanali ◽  
Steven W. Armfield ◽  
Michael P. Kirkpatrick ◽  
Stuart Norris

Parallel performance of a fractional step Navier-Stokes solver is investigated. Parallelisation is performed using Message Passing Interface, with domain partitioning. Block preconditioning is applied to the solution of the pressure Poisson equation, which is often the bottleneck in the computation of the fractional step method. Preconditioners tested are classes of incomplete matrix decompositions and sparse approximate inverses. The computational domain is decomposed into eight parts of about equal size in terms of the number of cells, and solved on eight parallel processors. Several aspects of the parallelisation, such as domain splitting directions, speed-up and scalability of the preconditioners, are discussed.


Volume 1 ◽  
2004 ◽  
Author(s):  
Sungsu Lee ◽  
Kyung-Soo Yang ◽  
Jong-Yeon Hwang

Development of geometry-independent computational method and educational codes for simulation of 2D flows around objects of complex geometry is presented. Referred as immersed boundary method, it introduces virtual forcing to governing equations to represent the effect of physical boundaries. The present method is based on a finite-volume approach on a staggered grid with a fractional-step method to solve Navier-Stokes equation and continuity equation. Both momentum and mass forcings are introduced on and inside the object to satisfy no-slip condition and mass conservation. Since Cartesian grid lines in general do not coincide with the immersed boundaries, several interpolation schemes are employed. Several examples are simulated using the method presented in this study and the results agree well with other results. Both user-friendly preprocessor with GUI and FORTRAN-based solver are open to the public for educational purposes.


Author(s):  
Francisco Jose´ de Souza ◽  
Aristeu Silveira Neto

Subgrid-scale modeling, which characterizes Large Eddy Simulation (LES), has been used to predict the behavior of a water-fed hydrocyclone operating without an air core. The governing equations were solved by a fractional step method on a staggered grid. The Smagorinsky subgrid-scale model was employed to account for turbulent effects. Numerical results actually capture the main features of the flow pattern and agree reasonably well with experiments, suggesting that LES represents an interesting alternative to classical turbulence models when applied to the numerical solution of fluid flows within hydrocyclones.


2000 ◽  
Author(s):  
J. Rafael Pacheco ◽  
Arturo Pacheco-Vega ◽  
Sigfrido Pacheco-Vega

Abstract A new approach for the solution of time-dependent calculations of buoyancy driven currents is presented. This method employs the idea that density variation can be pursued by using markers distributed in the flow field. The analysis based on the finite difference technique with the non-staggered grid fractional step method is used to solve the flow equations written in terms of primitive variables. The physical domain is transformed to a rectangle by means of a numerical mapping technique. The problems analyzed include two-fluid flow in a tank with sloping bottom and colliding density currents. The numerical experiments performed show that this approach is efficient and robust.


Author(s):  
Sung-Eun Kim ◽  
L. Srinivasa Mohan

Large eddy simulations were carried out for the flow around a hydrodynamically smooth, fixed circular cylinder at two Reynolds numbers, one at a subcritical Reynolds number (Re = 1.4 × 105) and the other at a supercritical Reynolds number (Re = 1.0 × 106). The computations were made using a parallelized finite-volume Navier-Stokes solver based on a multidimensional linear reconstruction scheme that allows use of unstructured meshes. Central differencing was used for discretization of both convection and diffusion terms. Time-advancement scheme, based on an implicit, non-iterative fractional-step method, was adopted in conjunction with a three-level, backward second-order temporal discretization. Subgrid-scale turbulent viscosity was modeled by a dynamic Smagorinsky model adapted to arbitrary unstructured meshes with the aid of a test-filter applicable to arbitrary unstructured meshes. The present LES results closely reproduced the flow features observed in experiments at both Reynolds numbers. The time-averaged mean drag coefficient, root-mean-square force coefficients and the frequency content of fluctuating forces (vortex-shedding frequency) are predicted with a commendable accuracy.


2015 ◽  
Vol 2015 ◽  
pp. 1-10 ◽  
Author(s):  
Daniel S. Abdi ◽  
Girma T. Bitsuamlak

A Navier-Stokes equations solver is parallelized to run on a cluster of computers using the domain decomposition method. Two approaches of communication and computation are investigated, namely, synchronous and asynchronous methods. Asynchronous communication between subdomains is not commonly used in CFD codes; however, it has a potential to alleviate scaling bottlenecks incurred due to processors having to wait for each other at designated synchronization points. A common way to avoid this idle time is to overlap asynchronous communication with computation. For this to work, however, there must be something useful and independent a processor can do while waiting for messages to arrive. We investigate an alternative approach of computation, namely, conducting asynchronous iterations to improve local subdomain solution while communication is in progress. An in-house CFD code is parallelized using message passing interface (MPI), and scalability tests are conducted that suggest asynchronous iterations are a viable way of parallelizing CFD code.


Author(s):  
Lucas I Finn ◽  
Bruce M Boghosian ◽  
Christopher N Kottke

We describe a software package designed for the investigation of topological fluid dynamics with a novel algorithm for locating and tracking vortex cores. The package is equipped with modules for generating desired vortex knots and links and evolving them according to the Navier–Stokes equations, while tracking and visualizing them. The package is parallelized using a message passing interface for a multiprocessor environment and makes use of a computational steering library for dynamic user intervention.


Author(s):  
Sotirios S. Sarakinos ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

In this study an academic Computational Fluid Dynamics (CFD) code, named Galatea-I, is described, which employs the Reynolds Averaged Navier–Stokes (RANS) equations along with the artificial compressibility method and the SST (Shear Stress Transport) turbulence model for the prediction of incompressible viscous flows. For the representation of the computational domain unstructured hybrid grids are utilized, composed of tetrahedral, prismatic and pyramidical elements, while for its discretization a node-centered finite-volume scheme is implemented. Galatea-I is enhanced with a parallelization method, which employs spatial domain decomposition, while the data exchange between processors/processes is performed with the use of the Message Passing Interface (MPI) protocol. In addition, a parallel agglomeration multigrid methodology has been incorporated to improve further its computational performance. The proposed code is validated against steady-state flow benchmark test cases, concerning laminar flow over a cubic cavity and a cylindrical surface, as well as turbulent flow over a rectangular wing with a NACA0012 airfoil. The obtained results, compared with these of corresponding reference solvers, reveal Galatea-I’s potential for simulation of inviscid, viscous laminar and turbulent incompressible flows.


Sign in / Sign up

Export Citation Format

Share Document