skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Survey of Recent Applications of the PISALE Code and PDE Framework
We review the basic equations, numerical solution techniques, and new application areas of a novel multi-purpose computer code framework, PISALE, for the solution of complex Partial Differential Equation (PDE) systems on modern computing platforms. We describe how the code solves equations in the fluid approximation using a novel combination of ALE and AMR methods. Sample problems from areas of ground water flow, high-speed impacts, and X-ray Free Electron Laser (XFEL) experiments are given.  more » « less
Award ID(s):
2005259
PAR ID:
10568319
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Think Mind
Date Published:
ISSN:
2308-4499
ISBN:
978-1-68558-107-7
Subject(s) / Keyword(s):
Adaptive Mesh Refinement Computational Fluid Dynamics Arbitrary Lagrangian Eulerian Methods Volume of Fluid High Performance Computing Surface Tension Ground Water Flow.
Format(s):
Medium: X
Location:
https://www.thinkmind.org/library/ADVCOMP/ADVCOMP_2023/advcomp_2023_1_50_20029.html
Sponsoring Org:
National Science Foundation
More Like this
  1. Every numerical general relativistic investigation starts from the solution of the initial value equations at a given time. Astrophysically relevant initial values for different systems lead to distinct sets of equations that obey specific assumptions tied to the particular problem. Therefore, a robust and efficient solver for a variety of strongly gravitating sources is needed. In this work, we present the OpenMP version of the Compact Object CALculator (COCAL) on shared memory processors. We performed extensive profiling of the core COCAL modules in order to identify bottlenecks in efficiency, which we addressed. Using modest resources, the new parallel code achieves speedups of approximately one order of magnitude relative to the original serial COCAL code, which is crucial for parameter studies of computationally expensive systems such as magnetized neutron stars, as well as its further development towards more realistic scenarios. As a novel example of our new code, we compute a binary quark system where each companion has a dimensionless spin of 0.43 aligned with the orbital angular momentum. 
    more » « less
  2. Memcomputing is a novel computing paradigm beyond the von–Neumann one. Its digital version is designed for the efficient solution of combinatorial optimization problems, which emerge in various fields of science and technology. Previously, the performance of digital memcomputing machines (DMMs) was demonstrated using software simulations of their ordinary differential equations. Here, we present the first hardware realization of a DMM algorithm on a low-cost FPGA board. In this demonstration, we have implemented a Boolean satisfiability problem solver. To optimize the use of hardware resources, the algorithm was partially parallelized. The scalability of the present implementation is explored and our FPGA-based results are compared to those obtained using a python code running on a traditional (von–Neumann) computer, showing one to two orders of magnitude speed-up in time to solution. This initial small-scale implementation is projected to state-of-the-art FPGA boards anticipating further advantages of the hardware realization of DMMs over their software emulation. 
    more » « less
  3. Abstract Hyperbolic Navier–Stokes equations replace the heat operator within the Navier–Stokes equations with a damped wave operator. Due to this second-order temporal derivative term, there exist no known bounded quantities for its solution; consequently, various standard results for the Navier–Stokes equations such as the global existence of a weak solution, that is typically constructed via Galerkin approximation, are absent in the literature. In this manuscript, we employ the technique of convex integration on the two-dimensional hyperbolic Navier–Stokes equations to construct a weak solution with prescribed energy and thereby prove its non-uniqueness. The main difficulty is the second-order temporal derivative term, which is too singular to be estimated as a linear error. One of our novel ideas is to use the time integral of the temporal corrector perturbation of the Navier–Stokes equations as the temporal corrector perturbation for the hyperbolic Navier–Stokes equations. 
    more » « less
  4. We present a fully parallel digital memcomputing solver implemented on a field-programmable gate array (FPGA) board. For this purpose, we have designed an FPGA code that solves the ordinary differential equations associated with digital memcomputing in parallel. A feature of the code is the use of only integer-type variables and integer constants to enhance optimization. Consequently, each integration step in our solver is executed in 96 ns. This method was utilized for difficult instances of the Boolean satisfiability (SAT) problem close to a phase transition, involving up to about 150 variables. Our results demonstrate that the parallel implementation reduces the scaling exponent by about 1 compared to a sequential C++ code on a standard computer. Additionally, compared to C++ code, we observed a time-to-solution advantage of about three orders of magnitude. Given the limitations of FPGA resources, the current implementation of digital memcomputing will be especially useful for solving compact but challenging problems. 
    more » « less
  5. Heterogeneous computing environments combining CPU and GPU resources provide a great boost to large-scale scientific computing applications. Code generation utilities that partition the work into CPU and GPU tasks while considering data movement costs allow researchers to develop high-performance solutions more quickly and easily, and make these resources accessible to a larger user base.We present developments for a domain-specific language (DSL) and code generation framework for solving partial differential equations (PDEs). These enhancements facilitate GPU-accelerated solution of the Boltzmann transport equation (BTE) for phonons, which is the governing equation for simulating thermal transport in semiconductor materials at sub-micron scales. The solution of the BTE involves thousands of coupled PDEs as well as complicated boundary conditions and solving a nonlinear equation that couples all of the degrees of freedom at each time step. These developments enable the DSL to generate configurable hybrid GPU/CPU code that couples accelerated kernels with user-defined code. We observed performance improvements of around 18X compared to a CPU-only version produced by this same DSL with minimal additional programming effort. 
    more » « less