Engineering topics
Lim, Hyun
Publications and source records attributed to Lim, Hyun.
Nonlinear evolution of quadratic gravity in 3 + 1 dimensions
We present a numerically stable system of ( 3 + 1 ) evolution equations for the nonlinear gravitational dynamics of quadratic-curvature corrections to General Relativity (quadratic gravity). We also report on the numerical implementation of these evolution equations. We recover a well-known linear instability and gather evidence that—aside from said instability—quadratic gravity exhibits a physically stable Ricci-flat subsector. In particular, we demonstrate that Teukolsky-wave perturbations of a Schwarzschild black hole as well as a full binary inspiral (evolved up to merger) remain Ricci flat throughout evolution. Furthermore, this suggests that, at least in vacuum, classical quadratic gravity can mimic General Relativity, even in the fully nonlinear strong-gravity regime.
ARES v1.x - Performance Portable Tool to Simulate Supernovae based on Parthenon Framework
Historically, codes for simulating supernovae (such as Arepo, FLASH or LEAFS) have been at the forefront of scientific high-performance computing to the immense computational resources required for full 3D simulations. However, given the shift towards heterogenous HPC architectures, many current-generation codes are at the risk of losing their competitiveness as they are only designed to run on homogeneous CPU-only systems. There exist several efforts to enable these codes for GPU’s, however, these efforts only consider specific architectures or vendors (e.g., implement only CUDA or HIP), limiting themselves to a small range of exascale computing systems. Frameworks such as Kokkos aim to provide a framework which is agnostic of the targeted architecture, enabling the development of performant and portable code. In the Ares code, we develop a performance portable tool to simulate supernovae based on the Parthenon Framework, which in turn uses Kokkos in the background. Here, the Parthenon Framework provides an interface to the underlying mesh-refinement routines, which form the backbone of our code. In addition, we incorporate the already existing Singularity-EOS toolkit to provide us with various equations of state, primarily the Helmholtz equation of state. We also include the JINA Reaclib as a basis for our nuclear network solver. Finally, we implement a gravity solver to complete the required physics. This setup will provide us with a minimal code base to simulate supernova in a similar style to the tried-and-tested Arepo code, but in a futureproof performance portable framework.
Simulating Neutron Stars with Solid Quark Cores: Rotations, Oscillations, and Binary Mergers [Slides]
Neutron stars feature complexity from GR, QM, particle physics and nuclear physics effects which cannot be investigated in the lab up to the densities in their core. GWs from binary neutron star mergers are on the rise to provide more insight into the physics inside NSs. There is a possibility that involves crystalline, color-superconducting quark matter, with implications towards altering the GW signal if NSs harbor this exotic matter. There are no 3D dynamical simulations of solid quark matter inside Neutron stars. We perform the first 3D dynamical simulations of solid quark matter inside NSs. These are performed in the Newtonian Approximation.
TEE-ACM2
TEE-ACM2 is a library that computes matrix chain multiplication efficiently on GPUs, using a blocking strategy to load data. The library will provide an energy efficient algorithm for chain Matrix Multiplication, by minimizing both computations and off chip data transfers on the GPUs.
Enhancing FleCSPH Features and Validation Tests [Slides]
Abstract not provided.
Nonlinear dynamics of quadratic gravity in spherical symmetry
We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (quadratic gravity) in the spherically-symmetric sector. The formulation relies on (i) the harmonic gauge to cast the evolution system into quasilinear form (ii) the Cartoon method to reduce to spherical symmetry in keeping with the harmonic gauge, and (iii) order reduction to first order (in time) by means of introducing auxiliary variables. The well posedness of the respective initial-value problem is numerically confirmed by evolving randomly perturbed flat-space and black-hole initial data. In conclusion, our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives.
Modeling Neutron Star Oscillations in a Fixed General Relativistic Background Including Solid Crust Dynamics
Explore the source record for details and available documents.
Overview of FleCSPH Solid Material Modeling Capabilities [Slides]
Abstract not provided.
The Ristra Project: FY20/21 Milestone Report
The ASC Advanced Technology Development and Mitigation (ATDM) sub-program was established in 2014 to develop new simulation tools operating on exascale-class computers to serve NNSA (see Appendix B). Over the course of ATDM, LANL management have set a strategy for exascale-class application codes that follows two supportive and mutually risk-mitigating paths: evolution for established production integrated design codes (IDCs) – with a strong pedigree within the user community – based upon existing programming paradigms(MPI+X); and a new start ATDM project, Ristra, a high-risk/high-reward push for a next-generation multi-physics, multi-scale simulation toolkit based on emerging advanced programming systems(with an initial focus on data-flow task-based models exemplified by Legion). The role of Ristra as the high-risk/high-reward path for LANL’s codes was fully consistent with the goals of ATDM as described in Appendix B, in particular its emphasis on evolving ASC capabilities through novel computing programming models and computing technologies.
Nonlinear Dynamics of Modified Gravity Theories [Slides]
Abstract not provided.
What’s Up Gravitational Waves! More Messages from Ripples of Spacetime
Abstract not provided.
Are we alone? Searching Intelligent Life Beyond the Earth [Slides]
The planetarium presentation focused on how we understand and search for ''life'' in the universe.
FleCSPH: The next generation FleCSIble parallel computational infrastructure for smoothed particle hydrodynamics
Not Available
The Other 95% [Slides]
This planetarium presentation shares some basic facts about the universe. It’s very OLD. It's very BIG. The visible Universe contains billions of galaxies. The Universe is EXPANDING. The expansion is speeding up. 95% of the Universe is dark. Dark Energy is discussed in more detail.
Safe Travel Guide to Black Holes [Slides]
This report includes slides for the Pajarito Environmental Education Center (PEEC) public planetarium show.
FleCSPH : A New SPH Code Based on FleCSI Framework for Multi-physics Problems [Slides]
FleCSPH is a parallel and distributed implementation of smoothed particle hydrodynamics (SPH) based on the FleCSI framework.