An official website of the United States government
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
- PAR ID:
- 10423315
- Publisher / Repository:
- Johns Hopkins Turbulence Databases
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract We present the HelioCubed, a high-order magnetohydrodynamic (MHD) code designed for modeling the inner heliosphere. The code is designed to achieve 4th order accuracy both in space and in time. In addition, HelioCubed can perform simulations on mapped grids, such as those based on cubed spheres, which makes it possible to overcome stability limitations caused by the geometrical singularity at the polar axis of a spherical grid, thus enabling substantially larger time steps. HelioCubed has been developed using the high-level Proto library, ensures performance portability across CPU and GPU architectures, and supports back-end implementations, e.g., CUDA, HIP, OpenMP, and MPI. The code is compatible with the HDF5 library, which facilitates seamless data handling for simulations and boundary conditions derived from semi-empirical and MHD models of the solar corona. While presenting the results of preliminary simulations, we demonstrate that our simulations are indeed performed with 4th order of accuracy. Our approach ensures that HelioCubed solves the MHD equations preserving the radial flow to machine round-off error even on cubed-sphere grids. Solar wind simulations are performed using the boundary conditions provided by the Wang–Sheeley–Arge coronal model of the ambient solar wind. It also allows us to to simulate coronal mass ejections using observation-driven flux rope models. These capabilities make HelioCubed a versatile and powerful tool to advance heliophysics research and space weather forecasting.more » « less
-
Direct numerical simulation (DNS) of a transitional boundary layer over a plate with an elliptical leading edge. Navier-Stokes was discretized on a curvilinear grid and solved using a finite volume DNS code. A fractional-step algorithm was adopted, and the spatial discretization was a staggered volume-flux formulation. The viscous terms were integrated in time implicitly using the Crank-Nicolson and the advections terms were treated explicitly using the Adams-Bashforth. Pressure was treated using implicit Euler in the δp-form. The pressure equation was Fourier transformed in the span, and the resulting Helmholtz equation was solved for every spanwise wavenumber using two-dimensional multi-grid. After the simulation has reached a statistical stationary state, 4701 frames 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). Since the grid is staggered, data at the wall are not stored in the database. However, JHTDB provides values in the region between the wall and the first grid point, y∈[0, 0.0036], using 4th-order Lagrange polynomial inter- and extrapolation. The y-locations of the grid points in the vertical direction can be downloaded from this text file.more » « less
-
ABSTRACT A long-standing problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models based on the traditional Chandrasekhar DF formula have been proposed, yet they suffer from fundamental ambiguities in the definition of some terms in Chandrasekhar’s formula when applied to real galaxies, as well as difficulty in evaluating continuous quantities from (spatially) discrete simulation data. In this work, we present a new sub-grid DF estimator based on the discrete nature of N-body simulations, which also avoids the ambiguously defined quantities in Chandrasekhar’s formula. We test our estimator in the gizmo code and find that it agrees well with high-resolution simulations where DF is fully captured, with negligible additional computational cost. We also compare it with a Chandrasekhar estimator and discuss its applications in real galactic simulations.more » « less
