The data is from a direct numerical simulation of forced isotropic turbulence on a 4096-cubed periodic grid, using a pseudo-spectral parallel code. The simulations are documented in Ref. 1. Time integration uses second-order Runge-Kutta. The simulation is de-aliased using phase-shifting and truncation. Energy is injected by keeping the energy density in the lowest wavenumber modes prescribed following the approach of Donzis & Yeung. After the simulation has reached a statistical stationary state, a frame of data, which includes the 3 components of the velocity vector and the pressure, are generated and written in files that can be accessed directly by the database (FileDB system).
more »
« less
Rotating Stratified Turbulence Dataset on 40963 Grid
The data is from a direct numerical simulation of rotating stratified turbulence on a 4096-cubed periodic grid using a pseudo-spectral parallel code, GHOST. The simulations are documented in Ref. 1. The relative strength of stratification vs. rotation is characterized by the ratio of the Brunt-Väisälä to inertial wave frequency, N/f = 4.95. The code solves the Boussinesq equations with a solid body rotation force acting as the only external forcing mechanism. Time integration uses fourth-order Runge-Kutta. The simulation is initialized with large-scale isotropic conditions on a coarser grid. As the simulation progresses resolution is increased, peaking with 4096-cubed at maximum dissipation. After the simulation has reached a statistical stationary state, 5 frames of data, which includes the 3 components of the velocity vector and the temperature fluctuations, are generated and written in files that can be accessed directly by the database (FileDB system).
more »
« less
- Award ID(s):
- 2103874
- NSF-PAR ID:
- 10423315
- Publisher / Repository:
- Johns Hopkins Turbulence Databases
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract We present a global kinetic plasma simulation of an axisymmetric pulsar magnetosphere with self-consistent
e ±pair production. We use the particle-in-cell method and log-spherical coordinates with a grid size 4096 × 4096. This allows us to achieve a high voltage induced by the pulsar rotation and investigate pair creation in a young pulsar far from the death line. We find the following: (1) The energy release ande ±creation are strongly concentrated in the thin, Y-shaped current sheet, with a peak localized in a small volume at the Y-point. (2) The Y-point is shifted inward from the light cylinder by ∼15% and “breathes” with a small amplitude. (3) The densee ±cloud at the Y-point is in ultrarelativistic rotation, which we call superrotation, because it exceeds corotation with the star. The cloud receives angular momentum flowing from the star along the poloidal magnetic field lines. (4) Gamma-ray emission peaks at the Y-point and is collimated in the azimuthal direction, tangent to the Y-point circle. (5) The separatrix current sheet between the closed magnetosphere and the open magnetic field lines is sustained by the electron backflow from the Y-point cloud. Its thickness is self-regulated to marginal charge starvation. (6) Only a small fraction of dissipation occurs in the separatrix inward of the Y-point. A much higher power is released in the equatorial plane, including the Y-point where the created densee ±plasma is spun up and intermittently ejected through the nozzle between the two open magnetic fluxes. -
The data is from a direct numerical simulation (DNS) of homogeneous buoyancy driven turbulence on a 1024-cubed periodic grid. (See README-HBDT.pdf linked document for equations and details.) The simulation was performed with the variable-density version of the petascale CFDNS code. The database covers both the buoyancy driven increase in turbulence intensity as well as the buoyancy mediated turbulence decay.more » « less
-
more » « less
Abstract Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observations. Binary population synthesis with full simulation of stellar structure and evolution is computationally expensive, requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code
POSYDON incorporates grids ofMESA binary star simulations that are interpolated to model large-scale populations of massive binaries. The traditional method of computing a high-density rectilinear grid of simulations is not scalable for higher-dimension grids, accounting for a range of metallicities, rotation, and eccentricity. We present a new active learning algorithm,psy-cris , which uses machine learning in the data-gathering process to adaptively and iteratively target simulations to run, resulting in a custom, high-performance training set. We testpsy-cris on a toy problem and find the resulting training sets require fewer simulations for accurate classification and regression than either regular or randomly sampled grids. We further applypsy-cris to the target problem of building a dynamic grid ofMESA simulations, and we demonstrate that, even without fine tuning, a simulation set of only ∼1/4 the size of a rectilinear grid is sufficient to achieve the same classification accuracy. We anticipate further gains when algorithmic parameters are optimized for the targeted application. We find that optimizing for classification only may lead to performance losses in regression, and vice versa. Lowering the computational cost of producing grids will enable new population synthesis codes such asPOSYDON to cover more input parameters while preserving interpolation accuracies. -
more » « less
Abstract Calibrating with detailed 2D core-collapse supernova (CCSN) simulations, we derive a simple CCSN explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behaviour of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ∼90 per cent of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code F ornax adheres to Nature and issues concerning the neutrino–matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.
