Computational fluid dynamics using Graphics Processing Units: Challenges and opportunities

A new paradigm for computing fluid flows is the use of Graphics Processing Units (GPU), which have recently become very powerful and convenient to use. In the past three years, we have implemented five different fluid flow algorithms on GPUs and have obtained significant speed-ups over a single CPU. Typically, it is possible to achieve a factor of 50-100 over a single CPU. In this review paper, we describe our experiences on the various algorithms developed and the speeds achieved

Recent progress and challenges in exploiting graphics processors in computational fluid dynamics

The progress made in accelerating simulations of fluid flow using GPUs, and the challenges that remain, are surveyed. The review first provides an introduction to GPU computing and programming, and discusses various considerations for improved performance. Case studies comparing the performance of CPU- and GPUbased solvers for the Laplace and incompressible Navier–Stokes equations are performed in order to demonstrate the potential improvement even with simple codes. Recent efforts to accelerate CFD simulations using GPUs are reviewed for laminar, turbulent, and reactive flow solvers. Also, GPU implementations of the lattice Boltzmann method are reviewed. Finally, recommendations for implementing CFD codes on GPUs are given and remaining challenges are discussed, such as the need to develop new strategies and redesign algorithms to enable GPU acceleration.Keywords: Graphics processing unit (GPU), Reactive flow, Computational fluid dynamics (CFD), Laminar flows, Turbulent flow, CUD

Doctor of Philosophy

dissertationAccording to a UN report, more than 50% of the total world's population resides in urban areas and this fraction is increasing. Urbanization has a wide range of potential environmental impacts, including those related to the dispersion of potentially dangerous substances emitted from activities such as combustion, industrial processing or from deliberate harmful releases. This research is primarily focused on the investigation of various factors which contribute to the dispersion of certain classes of materials in a complex urban environment and improving both of the fundamental components of a fast response dispersion modeling system - wind modeling and dispersion modeling. Specifically, new empirical parameterizations have been suggested for an existing fast response wind model for street canyon flow fields. These new parameterizations are shown to produce more favorable results when compared with the experimental data. It is also demonstrated that the use of Graphics Processing Unit (GPU) technology can enhance the efficiency of an urban Lagrangian dispersion model and can achieve near real-time particle advection. The GPU also enables real-time visualizations which can be used for creating virtual urban environments to aid emergency responders. The dispersion model based on the GPU architecture relies on the so-called "simplified Langevin equations (SLEs)" for particle advection. The full or generalized form of the Langevin equations (GLEs) is known for its stiffness which tends to generate unstable modes in particle trajectory, where a particle may travel significant distances in a small time step

Development of the marker and cell method for use with unstructured meshes.

The marker and cell method is an efficient co-volume technique suitable for the solution of incompressible flows using a Cartesian mesh. For flows around complex geometries the use of an unstructured mesh is desirable. For geometric flexibility an unstructured mesh implementation is desirable. A co-volume technique requires a dual orthogonal mesh, in the triangular case the Delaunay-Voronoi dual provides the means for determining this dual orthogonal mesh in an unstructured mesh framework. Certain mesh criteria must be placed on the Delaunay-Voronoi to ensure it meets the dual orthogonal requirements. The two dimensional extension of the marker and cell method to an unstructured framework is presented. The requirements of the mesh are defined and methods in their production are discussed. Initially an explicit time stepping scheme is implemented which allows efficient simulation of incompressible fluid flow problems. Limitations of the explicit time stepping scheme that were discovered, mean that high Reynolds number flows that require the use of stretched meshes cannot produce solutions in a reasonable time period. A semi-implicit time stepping routine removes this limitation allowing these types of flows to be successfully modelled. To validate the solvers accuracy and demonstrate its performance, a number of test cases are presented. These include the lid driven cavity, flow over a backward facing step, inviscid flow around a circular cylinder, unsteady flow around a circular cylinder, flow around an SD7003 aerofoil, flow around a NACA0012 aerofoil and flow around a multi element aerofoil. The investigation although revealing a high dependence on the quality of the mesh still demonstrates that accurate results can be obtained efficiently. The efficiency is demonstrated by comparison to the in-house 2D incompressible finite volume solver for flow around a circular cylinder. For this case the unstructured MAC method produced a solution four times faster than the finite volume code

Technology 2000, volume 1

The purpose of the conference was to increase awareness of existing NASA developed technologies that are available for immediate use in the development of new products and processes, and to lay the groundwork for the effective utilization of emerging technologies. There were sessions on the following: Computer technology and software engineering; Human factors engineering and life sciences; Information and data management; Material sciences; Manufacturing and fabrication technology; Power, energy, and control systems; Robotics; Sensors and measurement technology; Artificial intelligence; Environmental technology; Optics and communications; and Superconductivity

NASA Tech Briefs, December 1997

Topics: Design and Analysis Software; Electronic Components and Circuits; Electronic Systems; Physical Sciences; Materials; Computer Software; Mechanics; Manufacturing/Fabrication; Mathematics and Information Sciences; Books and Reports

Bibliography of Lewis Research Center Technical Publications announced in 1991

This compilation of abstracts describes and indexes the technical reporting that resulted from the scientific engineering work performed and managed by the Lewis Research Center in 1991. All the publications were announced in the 1991 issues of STAR (Scientific and Technical Aerospace Reports) and/or IAA (International Aerospace Abstracts). Included are research reports, journal articles, conference presentations, patents and patent applications, and theses

Interactive molecular docking with haptics and advanced graphics

Biomolecular interactions underpin many of the processes that make up life. Molecular docking is the study of these interactions in silico. Interactive docking applications put the user in control of the docking process, allowing them to use their knowledge and intuition to determine how molecules bind together. Interactive molecular docking applications often use haptic devices as a method of controlling the docking process. These devices allow the user to easily manipulate the structures in 3D space, whilst feeling the forces that occur in response to their manipulations. As a result of the force refresh rate requirements of haptic devices, haptic assisted docking applications are often limited, in that they model the interacting proteins as rigid, use low fidelity visualisations or require expensive propriety equipment to use. The research in this thesis aims to address some of these limitations. Firstly, the development of a visualisation algorithm capable of rendering a depiction of a deforming protein at an interactive refresh rate, with per-pixel shadows and ambient occlusion, is discussed. Then, a novel approach to modelling molecular flexibility whilst maintaining a stable haptic refresh rate is developed. Together these algorithms are presented within Haptimol FlexiDock, the first haptic-assisted molecular docking application to support receptor flexibility with high fidelity graphics, whilst also maintaining interactive refresh rates on both the haptic device and visual display. Using Haptimol FlexiDock, docking experiments were performed between two protein-ligand pairs: Maltodextrin Binding Protein and Maltose, and glutamine Binding Protein and Glucose. When the ligand was placed in its approximate binding site, the direction of over 80% of the intra-molecular movement aligned with that seen in the experimental structures. Furthermore, over 50% of the expected backbone motion was present in the structures generated with FlexiDock. Calculating the deformation of a biomolecule in real time, whilst maintaining an interactive refresh rate on the haptic device (> 500Hz) is a breakthrough in the field of interactive molecular docking, as, previous approaches either model protein flexibility, but fail to achieve the required haptic refresh rate, or do not consider biomolecular flexibility at all