High order resolution and parallel implementation on unstructured grids

The numerical solution of the two-dimensional inviscid Euler flow equations is given. The unstructured mesh is generated by the advancing front technique. A cell-centred upwind finite volume method has been adopted to discretize the Euler equations. Both explicit and point implicit time stepping algorithms are derived. The flux calculation using Roe's and Osher's approximate Riemann solvers are studied. It is shown that both the Roe and Osher's schemes produce an accurate representation of discontinuities (e.g. shock wave). It is also shown that better convergence performance has been achieved by the point implicit scheme than that by the explicit scheme. Validations have been done for subsonic and transonic flow over airfoils, supersonic flow past a compression corner and hypersonic flow past cylinder and blunt body geometries. An adaptive remeshing procedure is also applied to the numerical solution with the objective of getting improved results. The issue of high order reconstruction on unstructured grids has been discussed. The methodology of the Taylor series expansion is adopted. The calculation of the gradient at a reference point is carried out by the use of either the Green-Gauss integral formula or the least-square methods. Some recently developed limiter construction methods have been used and their performance has been demonstrated using the test example of the transonic flow over a RAE 2822 airfoil. It has been shown that similar pressure distributions are obtained by all limiters except for shock wave regions where the limiter is active. The convergence problem is illustrated by the mid-mod type limiter. It seems only the Venkatakrishnan limiter provides improved convergence. Other limiters do not appear to work as well as that shown in their original publications. Also the convergence history given by the least-square method appears better than that by the Green-Gauss method in the test

Proceedings, MSVSCC 2018

Proceedings of the 12th Annual Modeling, Simulation & Visualization Student Capstone Conference held on April 19, 2018 at VMASC in Suffolk, Virginia. 155 pp

Adaptive heterogeneous parallelism for semi-empirical lattice dynamics in computational materials science.

With the variability in performance of the multitude of parallel environments available today, the conceptual overhead created by the need to anticipate runtime information to make design-time decisions has become overwhelming. Performance-critical applications and libraries carry implicit assumptions based on incidental metrics that are not portable to emerging computational platforms or even alternative contemporary architectures. Furthermore, the significance of runtime concerns such as makespan, energy efficiency and fault tolerance depends on the situational context. This thesis presents a case study in the application of both Mattsons prescriptive pattern-oriented approach and the more principled structured parallelism formalism to the computational simulation of inelastic neutron scattering spectra on hybrid CPU/GPU platforms. The original ad hoc implementation as well as new patternbased and structured implementations are evaluated for relative performance and scalability. Two new structural abstractions are introduced to facilitate adaptation by lazy optimisation and runtime feedback. A deferred-choice abstraction represents a unified space of alternative structural program variants, allowing static adaptation through model-specific exhaustive calibration with regards to the extrafunctional concerns of runtime, average instantaneous power and total energy usage. Instrumented queues serve as mechanism for structural composition and provide a representation of extrafunctional state that allows realisation of a market-based decentralised coordination heuristic for competitive resource allocation and the Lyapunov drift algorithm for cooperative scheduling

Research summary, January 1989 - June 1990

The Research Institute for Advanced Computer Science (RIACS) was established at NASA ARC in June of 1983. RIACS is privately operated by the Universities Space Research Association (USRA), a consortium of 62 universities with graduate programs in the aerospace sciences, under a Cooperative Agreement with NASA. RIACS serves as the representative of the USRA universities at ARC. This document reports our activities and accomplishments for the period 1 Jan. 1989 - 30 Jun. 1990. The following topics are covered: learning systems, networked systems, and parallel systems

Distributed solution of Laplacian eigenvalue problems

The purpose of this article is to approximately compute the eigenvalues of the symmetric Dirichlet Laplacian within an interval $(0,\Lambda)$. A novel domain decomposition Ritz method, partition of unity condensed pole interpolation method, is proposed. This method can be used in distributed computing environments where communication is expensive, e.g., in clusters running on cloud computing services or networked workstations. The Ritz space is obtained from local subspaces consistent with a decomposition of the domain into subdomains. These local subspaces are constructed independently of each other, using data only related to the corresponding subdomain. Relative eigenvalue error is analysed. Numerical examples on a cluster of workstations validate the error analysis and the performance of the method.Comment: 28 page

Mesh generation for voxel -based objects

A new physically-based approach to unstructured mesh generation via Monte-Carlo simulation is proposed. Geometrical objects to be meshed are represented by systems of interacting particles with a given interaction potential. A new way of distributing nodes in complex domains is proposed based on a concept of dynamic equilibrium ensemble, which represents a liquid state of matter. The algorithm is simple, numerically stable and produces uniform node distributions in domains of complex geometries and different dimensions. Well-shaped triangles or tetrahedra can be created by connecting a set of uniformly-spaced nodes. The proposed method has many advantages and potential applications.;The new method is applied to the problem of meshing of voxel-based objects. By customizing system potential energy function to reflect surface features, particles can be distributed into desired locations, such as sharp corners and edges. Feature-preserved surface mesh can then be constructed by connecting the node set.;A heuristic algorithm using an advancing front approach is proposed to generate triangulated surface meshes on voxel-based objects. The resultant surface meshes do not inherit the anisotropy of the underlying hexagonal grid. However, the important surface features, such as edges and corners may not be preserved in the mesh.;To overcome this problem, surface features such as edges, corners need to be detected. A new approach of edge capturing is proposed and demonstrated. The approach is based on a Laplace solver with incomplete Jacobi iterations, and as such is very simple and efficient. This edge capturing approach combined with the mesh generation methods above forms a simple and robust technique of unstructured mesh generation on voxel-based objects.;A graphical user interface (GUI) capable of complex geometric design and remote simulation control was implemented. The GUI was used in simulations of large fuel-cell stacks. It enables one to setup, run and monitor simulations remotely through secure shell (SSH2) connections. A voxel-based 3D geometrical modeling module is built along with the GUI. The flexibility of voxel-based geometry representation enables one to use this technique for both geometric design and visualization of volume data

Semiannual report

This report summarizes research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, fluid mechanics, and computer science during the period 1 Oct. 1994 - 31 Mar. 1995