Researcher profile

Terrance Strother

Terrance Strother contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2015arXiv

Knudsen number dependence of 2D single-mode Rayleigh-Taylor fluid instabilities

We present a study of single-mode Rayleigh-Taylor instabilities (smRTI) with a modified Direct Simulation Monte Carlo (mDSMC) code in two dimensions. The mDSMC code is aimed to capture the dynamics of matter for a large range of Knudsen numbers within one approach. Our method combines the traditional Monte Carlo technique to efficiently propagate particles and the Point-of-Closest-Approach method for high spatial resolution. Simulations are performed using different particle mean-free-paths and we compare the results to linear theory predictions for the growth rate including diffusion and viscosity. We find good agreement between theoretical predictions and simulations and, at late times, observe the development of secondary instabilities, similar to hydrodynamic simulations and experiments. Large mean-free-paths favor particle diffusion, reduce the occurrence of secondary instabilities and approach the non-interacting gas limit.

preprint2014arXiv

Hydrodynamic Shock Wave Studies within a Kinetic Monte Carlo Approach

Kinetic approaches are routinely employed to simulate the dynamics of systems that are too rarified to be described by the Navier-Stokes equations. However, generally they are far too computationally expensive to be applied for systems that are governed by continuum hydrodynamics. In this paper, we introduce a massively parallelized test-particle based kinetic Monte Carlo code that is capable of modeling the phase space evolution of an arbitrarily sized system that is free to move in and out of the continuum limit. Using particle mean free paths which are small with respect to the characteristic length scale of the simulated system, we retrieve continuum behavior, while non-equilibrium effects are observed when the mean free path is increased. To demonstrate the ability of our code to reproduce hydrodynamic solutions, we apply a test-suite of classic hydrodynamic shock problems. Simulations using tens of millions of test-particles are found to reproduce the analytical solutions well.

preprint2014arXiv

Kinetic Simulations of Rayleigh-Taylor Instabilities

We report on an ongoing project to develop a large scale Direct Simulation Monte Carlo code. The code is primarily aimed towards applications in astrophysics such as simulations of core-collapse supernovae. It has been tested on shock wave phenomena in the continuum limit and for matter out of equilibrium. In the current work we focus on the study of fluid instabilities. Like shock waves these are routinely used as test-cases for hydrodynamic codes and are discussed to play an important role in the explosion mechanism of core-collapse supernovae. As a first test we study the evolution of a single-mode Rayleigh-Taylor instability at the interface of a light and a heavy fluid in the presence of a gravitational acceleration. To suppress small-wavelength instabilities caused by the irregularity in the separation layer we use a large particle mean free path. The latter leads to the development of a diffusion layer as particles propagate from one fluid into the other. For small amplitudes, when the instability is in the linear regime, we compare its position and shape to the analytic prediction. Despite the broadening of the fluid interface we see a good agreement with the analytic solution. At later times we observe the development of a mushroom like shape caused by secondary Kelvin-Helmholtz instabilities as seen in hydrodynamic simulations and consistent with experimental observations.

preprint2013arXiv

Building a Hydrodynamics Code with Kinetic Theory

We report on the development of a test-particle based kinetic Monte Carlo code for large systems and its application to simulate matter in the continuum regime. Our code combines advantages of the Direct Simulation Monte Carlo and the Point-of-Closest-Approach methods to solve the collision integral of the Boltzmann equation. With that, we achieve a high spatial accuracy in simulations while maintaining computational feasibility when applying a large number of test-particles. The hybrid setup of our approach allows us to study systems which move in and out of the hydrodynamic regime, with low and high particle densities. To demonstrate our code's ability to reproduce hydrodynamic behavior we perform shock wave simulations and focus here on the Sedov blast wave test. The blast wave problem describes the evolution of a spherical expanding shock front and is an important verification problem for codes which are applied in astrophysical simulation, especially for approaches which aim to study core-collapse supernovae.

preprint2013arXiv

Parallelization of Kinetic Theory Simulations

Numerical studies of shock waves in large scale systems via kinetic simulations with millions of particles are too computationally demanding to be processed in serial. In this work we focus on optimizing the parallel performance of a kinetic Monte Carlo code for astrophysical simulations such as core-collapse supernovae. Our goal is to attain a flexible program that scales well with the architecture of modern supercomputers. This approach requires a hybrid model of programming that combines a message passing interface (MPI) with a multithreading model (OpenMP) in C++. We report on our approach to implement the hybrid design into the kinetic code and show first results which demonstrate a significant gain in performance when many processors are applied.