Fine Grid Numerical Solutions of Triangular Cavity Flow

Numerical solutions of 2-D steady incompressible flow inside a triangular cavity are presented. For the purpose of comparing our results with several different triangular cavity studies with different triangle geometries, a general triangle mapped onto a computational domain is considered. The Navier-Stokes equations in general curvilinear coordinates in streamfunction and vorticity formulation are numerically solved. Using a very fine grid mesh, the triangular cavity flow is solved for high Reynolds numbers. The results are compared with the numerical solutions found in the literature and also with analytical solutions as well. Detailed results are presented

Finite volume simulation of 2-D steady square lid driven cavity flow at high reynolds numbers

In this work, computer simulation results of steady incompressible flow in a 2-D square lid-driven cavity up to Reynolds number (Re) 65000 are presented and compared with those of earlier studies. The governing flow equations are solved by using the finite volume approach. Quadratic upstream interpolation for convective kinematics (QUICK) is used for the approximation of the convective terms in the flow equations. In the implementation of QUICK, the deferred correction technique is adopted. A non-uniform staggered grid arrangement of 768x768 is employed to discretize the flow geometry. Algebraic forms of the coupled flow equations are then solved through the iterative SIMPLE (Semi-Implicit Method for Pressure-Linked Equation) algorithm. The outlined computational methodology allows one to meet the main objective of this work, which is to address the computational convergence and wiggled flow problems encountered at high Reynolds and Peclet (Pe) numbers. Furthermore, after Re > 25000 additional vortexes appear at the bottom left and right corners that have not been observed in earlier studies

Performance and Scalability of Finite Element Analysis for Distributed Parallel Computation

A two dimensional (h; p) finite element scheme for distributed parallel computation is developed. The approach is based on an element - by - element domain decomposition and is implemented on the NCUBE2 system. Example problems are used to demonstrate performance of the algorithm for a range of (h; p) and to validate a run-time model. The run-time model is then used to examine the scaling properties of conjugate gradients and the (h; p) FEM over a range of p. 1 Introduction Recently, considerable attention has been focused on the performance of iterative methods for solving standard finite difference and element discretizations on distributed memory multiprocessors. This includes both gradient type methods, such as conjugate gradients (CG), and multigrid methods (MG). Parallel schemes based on Domain Decomposition and block iteration have also been extensively investigated. Recent studies of such methods can be found in [4, 6, 9]. However, comparatively little has been done to date wi..

A study of application sensitivity to variation in message passing latency and bandwidth

This study measures the effects of changes in message latency and bandwidth for production-level codes on a current generation tightly coupled MPP, the Intel Paragon. Messages are sent multiple times to study the application sensitivity to variations in band - width and latency. This method preserves the effects of contention on the interconnection network. Two applications are studied, PCTH, a shock physics code developed at Sandia National Laboratories, and PSTSWM, a spectral shallow water code developed at Oak Ridge National Laboratory and Argonne National Laboratory. These codes are significant in that PCTH is a {open_quote}full physics{close_quotes} application code in production use, while PSTSWM serves as a parallel algorithm test bed and benchmark for production codes used in atmospheric modeling. They are also significant in that the message-passing behavior differs significantly between the two codes, each representing an important class of scientific message-passing applications

Finite difference and cubic interpolated profile lattice boltzmann method for prediction of two-dimensional lid-driven shallow cavity flow

In this paper, two-dimensional lid-driven cavity flow phenomena at steady state were simulated using two different scales of numerical method: the finite difference solution to the Navier–Stokes equation and the cubic interpolated pseudo-particle lattice Boltzmann method. The aspect ratio of cavity was set at 1, 2/3, 1/2 and 1/3 and the Reynolds number of 100, 400 and 1,000 for every simulation condition. The results were presented in terms of the location of the center of main vortex, the streamline plots and the velocity profiles at vertical and horizontal midsections. In this study, it is found that at the simulation of Reynolds numbers 100 and 400, both methods demonstrate a good agreement with each other; however, small discrepancies appeared for the simulation at the Reynolds number of 1,000. We also found that the number, size and formation of vortices strongly depend on the Reynolds number. The effect of the aspect ratio on the fluid flow behavior is also presented