Problems in computational mechanics involve higher order nonlinear differential equations with complex boundary conditions, which are difficult to solve analytically and need numerical methods to predict the approximate solution. A large number of mesh points are utilized for better accuracy of the numerical technique which results in storage and operations of a large amount of data. It is of utmost importance that the time taken to perform these calculations is reduced to a realizable scale. General purpose graphical processing units (GPGPUs) provide high number of floating point operations per second (FLOPS) and potentially offer the most efficient architecture to carry out large-scale calculations in computational mechanics. In the present work, an attempt has been made to reduce the computational time for obtaining numerical solution of heat transfer by conduction, laminar flow in a rectangular channel, lid-driven cavity, and flow past square cylinder by programming GPGPUs using compute unified device architecture (CUDA), while maintaining overall second-order accuracy.
[1]
Jonathan Cohen,et al.
Title: A Fast Double Precision CFD Code using CUDA
,
2009
.
[2]
F. Harlow,et al.
Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface
,
1965
.
[3]
Michael J. Flynn,et al.
Some Computer Organizations and Their Effectiveness
,
1972,
IEEE Transactions on Computers.
[4]
I. Orlanski.
A Simple Boundary Condition for Unbounded Hyperbolic Flows
,
1976
.
[5]
Suhas V. Patankar,et al.
Numerical prediction of vortex shedding behind a square cylinder
,
1992
.
[6]
Krishnamurthy Muralidhar,et al.
Computational Fluid Flow and Heat Transfer
,
2003
.
[7]
U. Ghia,et al.
High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method
,
1982
.
[8]
F. Durst,et al.
Accurate computations of the laminar flow past a square cylinder based on two different methods : lattice-Boltzmann and finite-volume
,
2000
.